前一篇文章介绍了表面效应单元家族的153单元，今天来谈谈154单元。Surf154称之为3D结构表面效应单元，主要用于施加各种荷载和表面效应，可覆盖在3D单元面积表面，例如Solid单元和Shell单元，通过esurf命令创建。
单元形状如下所示：
该单元可通过4~8个节点定义，对于三角形的单元形状，可通过合并节点实现，建立单元时，该单元的I、J节点默认与单元坐标系的X轴平行。
该单元也可以考虑单元厚度，通过定义材料的密度来考虑单元的重量，也可以考虑附加质量得作用，附加质量通过实常数定义。
该单元的荷载方向编号如下所示：
采用该单元进行加载得时候，选择合适的坐标系非常重要。单元坐标系的选择通过Keyopt(2)进行设置。
不管是采用单元坐标系还是局部坐标系，单元主要荷载方向为1、2、3，分别沿着当前坐标系的Z、X以及Y方向。
至于4号、5号荷载方向则同Surf153的3号、4号荷载方向类似，一般情况下很少采用4号荷载方向，5号荷载方向则是用于施加斜向荷载，方向通过单元向量决定。
下面以几个简单的例子来说明利用该单元的使用。（部分例子来源于ANSYS结构单元分析与应用）
例题一：直角坐标系下的施加，施加如图所示Q1、Q2、Q3、Q4荷载
命令流如下：
finish
/clear
/prep7
l=800
b=100
h=200
q1=3
q2=2
q3=6
q4=0.5
ref=40*acos(-1)/180 !Q1的角度
et,1,solid95
et,2,surf154
keyopt,2,11,2
mp,ex,1,2.1e5
mp,prxy,1,0.3
blc4,,,l,h,b
vatt,1,,1
esize,25
vmesh,all
nsel,s,loc,x,0
d,all,all
nsel,s,loc,y,h
type,2
esurf
nsel,s,loc,x,l
esurf
allsel,all
!施加顶面的各种荷载
nsel,s,loc,y,h
esln,s,1
sfe,all,2,pres,,-q3
sfe,all,3,pres,,-q2
sfe,all,5,pres,,q1,-cos(ref),-sin(ref),0
!施加端面的切向荷载
nsel,s,loc,x,l
esln,s,1
sfe,all,2,pres,,-q4
allsel,all
例题二：切向荷载施加，对某圆柱体施加如下所示的切向荷载
该例子在使用表面效应单元时，要注意转换坐标系，新建一局部坐标系，且该局部坐标系为极坐标系，将单元坐标系旋转到局部坐标系下， 也可以设置Keyopt(2)=1,根据局部坐标系来施加，施加荷载的方向为Y方向，也即是2号荷载方向。命令流如下：
finish
/clear
/prep7
l=400
r=50
et,1,solid95
et,2,surf154
mp,ex,1,2.1e5
mp,prxy,1,0.0
cyl4,,,r,,,,l
wprota,,90
vsbw,all
wprota,,,90
vsbw,all
vatt,1,,1
esize,15
vmesh,all
nsel,s,loc,z,0
d,all,all
csys,1
nsel,s,loc,x,r
local,11,1 !定义局部柱坐标系
esys,11 !将单元坐标系指定在局部坐标系
type,2
esurf
esel,s,type,,2
sfe,all,2,pres,,10
例题三：某网壳均布荷载的施加。
在某种情况下，需要按照投影面积来施加荷载，一个很典型的例子便是坡屋顶时，均布荷载的施加。例如，如下的网壳结构，我们在计算时需要将屋面均布荷载2KN/m^2等效为节点荷载，由于整个网壳结构采用杆梁系单元计算，因此均布荷载的施加可以通过表面效应单元Surf154来实现，通过约束所有节点自由度，计算结构在均布荷载作用下的支座反力，从而实现等效节点荷载的计算。
该过程部分命令流如下：
/prep7
!定义表面效应单元
local,11,0,0,0,0
ET,3,SURF154,,1,,1
!定义面单元surf154 ，局部坐标系，无中间节点
R,3,,,,,,,0.1,0.1,0.1 !定义实常熟
MP,DENS,3,0 !定义材料密度
TYPE,3 !设置单元类型
MAT,3 !设置材料属性
REAL,3 !设置实常熟
esys,11
…….
!================
!创建单元
!===============
!均布荷载施加
P=2000
allsel,all
SFE,all,3,PRES,, P
!===============
计算完后，获取节点反力，便是等效节点荷载，将节点荷载施加到相应节点上，并同时删掉表面效应单元，进行下一步的计算。
暂无评论内容