基于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命令施加沿某个坐标轴方向荷载值变化的荷载。ansys中加变化的面载荷将在下文中介绍。