今日和各位分享一个后处理小技巧，即杆梁系结构在ANSYS中如何像设计软件那般显示构件应力比的云图。
用过设计软件诸如3D3S、Midas、YJK的同学都知道，在复核杆件强度和稳定性的时候，这些软件可以一键生成构件应力比云图或者稳定性云图，这样可以非常方便的知道在整个结构体系中哪些杆件处于薄弱，需要适当加强，从而指导后续的结构设计工作。
然而对于通用有限元软件来讲，由于没有套我国设计规范，有限元软件只能做整体性结构分析，杆件复核做不到，也即后续的杆件强度、稳定复核需要人工去实现。
本文以某单层网壳结构为例，具体演示如何实现强度应力比云图在ANSYS中的显示。
根据《钢标》，杆件强度复核公式如下：
首先在做完工况分析后，根据上述公式，需要手动提取杆件的轴力、杆件的Mx、My,并按公式8.1.1-2计算杆件应力，两边同时除以f，则可到强度应力比，命令流如下：
!根据《钢标》3.5.1确定宽厚比等级，截面等级为S1
Wn1=3.14*325**3*(1-Ab1**4)/32
An1=3.14/4*325**2*(1-Ab1**2)
Wn2=3.14*325**3*(1-Ab2**4)/32
An2=3.14/4*325**2*(1-Ab2**2)
Wn3=3.14*325**3*(1-Ab3**4)/32
An3=3.14/4*325**2*(1-Ab3**2)
*get,ele_max,elem,0,count
*dim,Stress_ratio,,ele_max
*dim,Ele_xuhao,char,ele_max
*get,Real_num,elem,i,attr,real
*elseif,Real_num,eq,2,then
*elseif,Real_num,eq,3,then
*get,Mz_i(i),elem,i,smisc,3
*get,Mz_j(i),elem,i,smisc,16
*get,My_i(i),elem,i,smisc,2
*get,My_j(i),elem,i,smisc,15
Mzy_i=(Mz_i(i)**2+My_i(i)**2)**0.5
Mzy_j=(Mz_j(i)**2+My_j(i)**2)**0.5
*get,Axial_i(i),elem,i,smisc,1
*get,Axial_j(i),elem,i,smisc,14
!计算应力比，注意Q345对应的抗拉强度设计值为305Mpa
Stress_ratio_i=abs(Axial_i(i))/An/fn+Mzy_i/rm/Wn/fn
Stress_ratio_j=abs(Axial_j(i))/An/fn+Mzy_j/rm/Wn/fn
Stress_ratio(i)=max(Stress_ratio_i,Stress_ratio_j)
上述数组 Stress_ratio 即为根据单元序号存储的单元强度应力比，接下来需要将该应力比反应到云图当中。
esel,s,ename,,188
etable,ZL,smisc,1
pletab,ZL
*do,i,1,ele_max
*vput,Stress_ratio(i),elem,i,etab,ZL
*enddo
pletab,ZL
上述命令流中，首先定义了一个单元表，这个单元表可随意定义，这里定义的是轴力单元表，然后用vput命令覆盖单元表里面的数据，将原来轴力的数据替换为了计算的应力比，最后重新绘制单元表即可显示的是应力比云图。
从应力比云图则可清楚的了解杆件的应力分布，最大应力比，最小应力比等，从而对后续的结构优化提供指导意义。
稳定系系数云图的显示工作原理同理，这里就不在单独演示。
