[转]ANSYS中用表面效应单元加任意方向的荷载

(本文转自:徐瑞有限元

 !用表面效应单元加任意方向的荷载
finish
/PREP7
et,1,45    !定义实体单元solid45
et,2,154   !定义三维表面效应单元

KEYOPT,2,2,0    !指定表面效应单元的K2=0,所加荷载与单元坐标系方向相同
KEYOPT,2,4,1    !指定表面效应单元的K4=0,去掉边中点,成为四结点表面单元

block,-5,5,-5,5,0,5      !建实体模型
mp,dens,1,2000
mp,ex,1,10e9
mp,prxy,1,0.2

asel,s,loc,z,5.0,5.0                !选中实体上表面
AATT,       1, ,   2,       0,      !指定实体上表面用154号单元
MSHAPE,0,2D
MSHKEY,1
esize,,5
amesh,all                   !对上表面划分网格

allsel,all
VATT,       1, ,   1,       0       !指定实体用45号单元
MSHAPE,0,3D
MSHKEY,1
vmesh,all

/PSYMB,ESYS,1  !显示单元坐标系
esel,s,type,,2                !选中实体上表面的表面效应单元以方便加荷载

sfe,all,1,pres,,50   !在面内加Z向荷载,大小为50,荷载方向可通过值的正负控制
sfe,all,2,pres,,100   !在面内加X向荷载,大小为100
sfe,all,3,pres,,150   !在面内加Y向荷载,大小为150

/psf,pres,,2,0,1   !以箭头方式显示所加荷载
!如果已经知道荷载在整体坐标系内的方向失量为(0,1,1),可以用如语句加该方向的荷载
sfe,all,5,pres,,100,0,1,1 !荷载值100后的三个数为方向失量

allsel,all
eplot


通过以上命令流得到的荷载图如下图所示:

表面单元施加荷载

需要注意的时图中(0,1,1)方向的荷载值为70.71=100*sqrt(2)/2,刚好是命令流中的荷载值乘以方向余弦。可以用sfelist命令查看单元上的荷载值。

另外,可以再结合sfgrad命令施加沿某个坐标轴方向荷载值变化的荷载。

P.S 在Workbench中对实体面施加的边界条件,一般都是通过这种方法实现的,只是创建这个面单元的工作由Workbench后台自动执行罢了。