杆梁壳单元内力的输出
计算完成后,就可以对杆梁壳单元进行后处理,显示位移、应力(需打开单元开关显示开关)的操作和普通的单元的后处理是一样的。
但其内力的输出和显示需要进行特殊的处理。
首先要保证在定义单元类型时,打开了内力输出的选项;其次,在完成计算后需要通过定义单元表的方法来提取。
一、定义输出选项
Beam4单元说明的input summary中可以看到关于内力输出的开关选项说明:
二、定义单元表
查看对应项和序号:在单元说的后的单元输出表中查到MFOR(X,Y,Z)和MMOM(X,Y,Z)对应的项和序列号如下表,注意I和J节点上对应的序列号不同。
GUI方式如下图:
命令流过程如下:
(I节点处结果)
Etable,forx_i,smisc,1 !x方向剪力(即轴力)
Etable,fory_i,smisc,2 !y方向剪力
Etable,forz_i,smisc,3 !z方向剪力
Etable,momx_i,smisc,4 !x方向弯矩
Etable,momy_i,smisc,5 !y方向弯矩
Etable,momz-i,smisc,6 !z方向弯矩
(J节点处结果)
Etable,forx_j,smisc,7 !x方向剪力(即轴力)
Etable,fory_j,smisc,8 !y方向剪力
Etable,forz_j,smisc,9 !z方向剪力
Etable,momx_j,smisc,10 !x方向弯矩
Etable,momy_j,smisc,11 !y方向弯矩
Etable,momz-j,smisc,12 !z方向弯矩
三、绘制内力图
对线单元和面单元,其绘制单元表的结果的命令是不一样的。
线单元的内力云图的绘制为:PLLS. 面单元的内力云图的绘制为:PLETAB
四、注意事项
1.弯矩图的绘制与显示比例
在绘制弯矩图时(其他内力图也一样),可以给定一个显示比例系数Fact(默认等于1),当为负值时,反向显示,数值为放大系数。
同时还须注意的是,弯矩图的显示还与单元坐标系的方向有关,因此,如果同一段梁中各个单元的单元坐标系不一致(可以通过改变线的方向来改变线单元的单元坐标系),还将导致弯矩图显示不连续。
壳单元的后处理中更加要注意单元坐标系的方向,而且如果面上的壳单元划分不整齐,还将导致提取出来的内力图方向不统一,没有参考价值。
2.划分数的影响
要特别注意梁单元划分数对弯矩图的影响,除了因为有限元法是一种把连续结构体离散化的近似求解方法,还因为弯矩图是以梁单元的I和J节点结果进行连线绘制的。
可通过例子中用10份,2份,100份进行对比。