近期学习到的一个方法，可供参考，也欢迎探讨：

1、在 LS-Prepost 中先通过 F11 回到老版界面，Setting 板块中，Local Coord System >>Whole>>cyl>>Create 创建柱坐标系；

2、在柱坐标创建面板中采用 2Position 方式建立，其中 OriginP 的 x0,y0,z0 代表柱坐标系原点，两个 Position（x1,y1,z1）（x2,y2,z2）分别与（x0,y0,z0) 组合表示柱坐标系的轴以及过轴面（暂且这样称）。例：假如柱坐标系轴线为 z 轴，则（x1=x0,y1=y0,z1=z0+/-1)，对应的 axis 选择 Z_axis；过轴面选择 YZ_Axis 时，（x2=x0,y2=y0+1,z2= z0)，如图定义了一个在（0，71.6249，0）位置的柱坐标系。

3、然后选中定义的坐标系，点击 Apply，再回到 Fcomp 中，选中单元，输出单元应力，x-stress 就变成了径向应力，y-stress 就变成了环向应力。

通过柱坐标查看应力分布，首先需要明确的是数值仿真模型中使用的材料本构模型是否支持柱坐标，部分材料卡（*MAT) 中有一个 AOPT 选项，需要将其设置为：

AOPT = 4

计算模型切换到圆柱坐标。
此外，还需要在下一行（XP,YP,ZP) 中定义原点 P 的坐标，并定义沿 z 轴方向的矢量（V1,V2,V3)。

您可以通过：

Element Tool -> Measure -> Inertia -> Active Elements Only -> All -> Apply

数值模型的重心设置为柱坐标的原点。

如果在 LS-DYNA 中无法执行此操作，则仍可以使用一些方程将应力从笛卡尔坐标系转换为圆柱坐标系。