!用表面效應(yīng)單元加任意方向的荷載
finish
/PREP7
et,1,45 !定義實(shí)體單元solid45
et,2,154 !定義三維表面效應(yīng)單元

KEYOPT,2,2,0 !指定表面效應(yīng)單元的K2=0，所加荷載與單元坐標(biāo)系方向相同
KEYOPT,2,4,1 !指定表面效應(yīng)單元的K4=0，去掉邊中點(diǎn)，成為四結(jié)點(diǎn)表面單元

block,-5,5,-5,5,0,5 !建實(shí)體模型

mp,dens,1,2000
mp,ex,1,10e9
mp,prxy,1,0.2

asel,s,loc,z,5.0,5.0 !選中實(shí)體上表面
AATT, 1, , 2, 0, !指定實(shí)體上表面用154號單元
MSHAPE,0,2D
MSHKEY,1
esize,,5
amesh,all !對上表面劃分網(wǎng)格

allsel,all
VATT, 1, , 1, 0 !指定實(shí)體用45號單元
MSHAPE,0,3D
MSHKEY,1
vmesh,all

/PSYMB,ESYS,1 !顯示單元坐標(biāo)系
esel,s,type,,2 !選中實(shí)體上表面的表面效應(yīng)單元以方便加荷載

sfe,all,1,pres,,50 !在面內(nèi)加Z向荷載，大小為50，荷載方向可通過值的正負(fù)控制
sfe,all,2,pres,,100 !在面內(nèi)加X向荷載，大小為100
sfe,all,3,pres,,150 !在面內(nèi)加Y向荷載，大小為150

/psf,pres,,2,0,1 !以箭頭方式顯示所加荷載

!如果已經(jīng)知道荷載在整體坐標(biāo)系內(nèi)的方向失量為(0,1,1)，可以用如語句加該方向的荷載
sfe,all,5,pres,,100,0,1,1 !荷載值100后的三個數(shù)為方向失量

allsel,all
eplot

通過以上命令流得到的荷載圖如下

需要注意的時圖中(0,1,1)方向的荷載值為70.71=100*sqrt(2)/2，剛好是命令流中的荷載值乘以方向余弦。可以用sfelist命令查看單元上的荷載值。

另外，可以再結(jié)合sfgrad命令施加沿某個坐標(biāo)軸方向荷載值變化的荷載。可以參考“[url=http://blog.sina.com.cn/s/blog_47569d4601000aap.html]ANSYS中加變化的面荷載的方法”

·

在ANSYS中如果要在一個面上施加沿某個方向變化的面荷載，需要有兩步來完成：

這里以一個在圓筒內(nèi)表面加內(nèi)水壓力的例子進(jìn)行說明。

第一步，設(shè)置面荷載變化規(guī)律。如果面荷載沿 Z 向變化，后面指定面荷載從 Z=100 開始變化，并按斜率為 -9800 進(jìn)行變化，可用如下語句

sfgrad,pres,,z,100,-9800 ! 也就是準(zhǔn)備在高 100 米的圓柱加內(nèi)水壓力吧
第二步，施加面荷載。在指定的面上施加按第一步設(shè)置的面荷載變化規(guī)律的面荷載。
SFA,P51X,1,PRES,0
這個語句相當(dāng)于在指定面上施加法向荷載 ( 選圓筒體內(nèi)表面 ) ，在 Z=100 時荷載值為 0 ，隨 Z 坐標(biāo)變化荷載值以變化率 -9800 進(jìn)行變化，這樣在 Z=0 時荷載值為 -9800*100
每次用 sfgrad 進(jìn)行設(shè)置后僅對隨后的 sfa 命令有效，直倒下次再用 sfgrad 進(jìn)行設(shè)置。
在面上施加荷載后，對模型剖分后可以執(zhí)行以下命令來查看加的面荷載是否正確
/PSF,PRES,NORM,2,0,1 以箭頭方式顯示面荷載
sftran 將面荷載轉(zhuǎn)化到有限元模型上

文章引用自：