计算机图形学之动画和模拟算法：头发和毛发的物理模型.docx

#

模拟循环

for t in range(1000):

#

计算拉伸力

forces = np.zeros((n_nodes, 3))

for i in range(1, n_nodes):

delta_pos = positions[i] - positions[i-1]

stretch_force = elastic_coeff * (np.linalg.norm(delta_pos) - rest_length) * delta_pos / np.lina

lg.norm(delta_pos)

forces[i] -= stretch_force

forces[i-1] += stretch_force