我刚刚浏览了 Fenics 书以获取节点位移并尝试了此代码....

u_nodal_values = U.vector()
u_array = u_nodal_values.array()
coor = mesh.coordinates()
center = (0.5,0.5,0.5)
print U(center)
if mesh.num_vertices()==len(u_array):
   for i in range(mesh.num_vertices()): 
       print ’U(%8g,%8g,%8g) = %g’% (coor[i][0], coor[i][1], u_array[i])

但我没有得到预期的结果..有人可以建议我修改吗？