!用表面效應單元加任意方向的荷載

　　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

　　通過以上命令流得到的荷載圖如附件所示。