用Table数组参数实现线性荷载的施加

ANSYS中，命令"SFL"可以实现线上施加面荷载的功能，但其只对2D面单元的边界（线）、轴对称单元本身、壳单元边界（线）有效，对3D梁单元无效。命令"SF"和"SFBEAM"一次也只可以对一个单元进行线性荷载的施加。本实例通过使用具有线性插值功能的Table数组参数，快速实现对多个梁单元线性渐变荷载的施加。修改Table数组参数可实现任意线性变化的荷载。（注：本例系在汤子扬网友的命令流基础上修改而成，在此表示感谢。）

命令流如下：

FINI
/CLEAR
/PREP7
ELEM_L=0.25
ET,1,BEAM3
MP,EX,1,207E9
MP,PRXY,1,0.3
MP,DENS,1,7.8E3
R,1,0.03,0.25E-4,0.1
K,1
K,2,6
K,3,7
K,4,8
L,1,2
L,2,3
L,3,4
LSEL,ALL
LATT,1,1,1
LESIZE,ALL,ELEM_L
LMESH,ALL
D,1,UX,0,,,,UY
D,26,UY,0
!定义Table数组参数表示线性荷载
*DIM,XXHZ,TABLE,3,1
XXHZ(1,0)=0,0.5,1
XXHZ(1)=800,400,600

ALLSEL,ALL
L_LOAD=1    !加载直线编号

*GET,LL,LINE,L_LOAD,LENG
*GET,K1,LINE,L_LOAD,KP,1
N1=NODE(KX(K1),KY(K1),KZ(K1))
LSEL,S,LINE,,L_LOAD
ESLL,S

*GET,MAXN,ELEM,,NUM,MAX
*GET,MINN,ELEM,,NUM,MIN
*SET,I,MINN
*DO,I,MINN,MAXN
   *IF,ESEL(I),EQ,1,THEN
     KKI=DISTND(NELEM(I,1),N1)/LL
     KKJ=DISTND(NELEM(I,2),N1)/LL
     FI=XXHZ(KKI)
     FJ=XXHZ(KKJ)
     SFBEAM,I,1,PRES,FI,FJ
   *ENDIF
*ENDDO

ALLSEL,ALL
/PBC,ALL,,1
EPLOT

/SOLU
!...

运行结果如下图：

---

【新鲜一下】：厌倦了Google和Baidu的搜索框了没？看看这个网页版的命令行搜索goosh吧，搜索只需要简单几个命令就可以了，ft的是俺用它来搜索"鲁班人"关键词，第2条竞然出来“以建筑施工为主并从事房地产等多元经营的国有大型建筑施工企业”，强列声明，与偶无关。纳闷儿中...