FLAC3D整体坐标系的变换

z511367698

附命令：如x与y互换
def rotation
p_gp=gp_head
loop while p_gp # null
w_temp1=gp_xpos(p_gp)
w_temp2=gp_ypos(p_gp)
w_temp3=gp_zpos(p_gp)
gp_xpos(p_gp)=-w_temp2
gp_ypos(p_gp)=w_temp1
gp_zpos(p_gp)=w_temp3
p_gp = gp_next(p_gp)
endloop
end
@rotation

cumt_em17

我也像你这样尝试过，但网格可以转变，计算的时候有时会提示体积为负，想要旋转的话需要遵守坐标的轮换对称性，即x->y->z->x的顺序才行。