如何为随时间变化的折叠物体实现重力模拟？

核心思路 #

因为你需要关节无限刚性且旋转由正弦函数强制控制，而非物理引擎自动计算，直接用Pymunk这类通用物理引擎会受其内置约束限制，自定义轻量物理逻辑反而更可控。核心是：将每个部件视为刚体，用正弦函数固定关节旋转角度，手动处理重力、地面碰撞与摩擦驱动。

具体实现步骤 #

• 部件与关节建模：每个部件用矩形/线段表示，记录位置、角度、质量、关节锚点。关节旋转角度由sine(time)直接控制，关联的两个部件需根据锚点坐标同步位置（保证刚性）。
• 重力与碰撞处理：给每个部件施加重力加速度g（比如(0, 9.8)），实时检测部件是否接触地面（y坐标超过地面高度），碰撞后将垂直速度置0，同时保留水平方向的摩擦减速。
• 摩擦驱动逻辑：当部件绕关节旋转时，接触地面的部件边缘会因摩擦产生水平推力。计算部件接触点的线速度，根据摩擦系数修正整体的水平速度，实现“行走”效果。

代码示例（核心片段） #

import pygame
import math

pygame.init()
screen = pygame.display.set_mode((800, 600))
clock = pygame.time.Clock()
ground_y = 500
g = 0.5  # 重力加速度
friction = 0.95  # 地面摩擦系数

# 定义部件类
class Segment:
    def __init__(self, x, y, length, angle, mass=1):
        self.x = x
        self.y = y
        self.length = length
        self.angle = angle
        self.mass = mass
        self.vx = 0
        self.vy = 0
        # 关节锚点（部件末端）
        self.joint_x = self.x + math.cos(self.angle) * self.length
        self.joint_y = self.y + math.sin(self.angle) * self.length

    def update(self, ground_y):
        # 施加重力
        self.vy += g
        # 地面碰撞检测
        if self.y + self.length * math.sin(self.angle) > ground_y:
            self.vy = 0
            self.y = ground_y - self.length * math.sin(self.angle)
            # 摩擦减速
            self.vx *= friction
        # 更新位置
        self.x += self.vx
        self.y += self.vy
        # 更新关节锚点
        self.joint_x = self.x + math.cos(self.angle) * self.length
        self.joint_y = self.y + math.sin(self.angle) * self.length

    def draw(self, screen):
        end_x = self.x + math.cos(self.angle) * self.length
        end_y = self.y + math.sin(self.angle) * self.length
        pygame.draw.line(screen, (255,0,0), (self.x, self.y), (end_x, end_y), 5)
        pygame.draw.circle(screen, (0,255,0), (int(self.joint_x), int(self.joint_y)), 5)

# 初始化两个关联部件
seg1 = Segment(300, 400, 100, -math.pi/2)
seg2 = Segment(seg1.joint_x, seg1.joint_y, 100, 0)

running = True
while running:
    screen.fill((255,255,255))
    pygame.draw.line(screen, (0,0,0), (0, ground_y), (800, ground_y), 3)
    
    for event in pygame.event.get():
        if event.type == pygame.QUIT:
            running = False

    # 用正弦函数控制关节旋转（seg2的角度随时间变化）
    t = pygame.time.get_ticks() / 1000
    seg2.angle = math.sin(t) * math.pi/3
    # 同步seg2的位置到seg1的关节锚点（保证刚性）
    seg2.x = seg1.joint_x
    seg2.y = seg1.joint_y

    # 计算摩擦驱动：seg2接触地面时，旋转产生的水平推力
    seg2_end_y = seg2.y + math.sin(seg2.angle) * seg2.length
    if seg2_end_y >= ground_y:
        # 计算接触点的水平线速度
        rot_speed = math.cos(t) * math.pi/3  # 正弦函数的导数是余弦
        contact_vx = rot_speed * seg2.length * math.cos(seg2.angle + math.pi/2)
        seg1.vx += contact_vx * 0.1  # 推力系数，可调整

    seg1.update(ground_y)
    seg2.update(ground_y)

    seg1.draw(screen)
    seg2.draw(screen)

    pygame.display.flip()
    clock.tick(60)

pygame.quit()

关键说明 #

• 关节刚性通过强制同步部件位置实现：seg2的起点始终绑定seg1的末端，角度由正弦函数直接设置，不受物理影响。
• 行走驱动依赖旋转部件与地面的摩擦：当seg2旋转时，接触地面的端点产生水平方向的速度分量，通过摩擦转化为整体的移动推力。
• 可调整参数：g（重力）、friction（摩擦系数）、正弦函数的频率/振幅，以及推力系数，来优化行走效果。

内容的提问来源于stack exchange，提问作者Material Scientist