本帖最后由 civilman狼人 于 2012-10-26 22:28 编辑
首先根据表中的初始条件求出S(x)=S0*e(tx)中的两个系数S0和t(看不清该字母,用t来表示),即120/2=S0*e(t*0),500/2=S0*e(t*950),求得S0=60,t=(ln(250/60))/950=1.502227743E-03。
最后应用ANSYS的APDL语言进行编程即可实现。
具体见如下步骤:
finish
/clear
S0=60
t=(log(250/60))/950
/prep7
*do,i,0,95,1
k,i+1,i*10,S0*((exp(1))**(t*i*10))
*enddo
*do,i,1,95,1
l,i,i+1
*enddo
k,1000,0,0,0
k,1001,950,0,0
AROTAT,all, , , , , , 1000, 1001,360,12,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!上面命令流为几何图形,以下为划分为单元的命令流
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
finish
/clear
/prep7
ET,1,SHELL63
R,1,0.1, , , , , ,
MP,EX,1,2.1e11
MP,PRXY,1,0.3
MP,DENS,1,7850
*do,i,0,95,1
k,i+1,i*10,S0*((exp(1))**(t*i*10))
*enddo
*do,i,1,95,1
l,i,i+1
*enddo
k,1000,0,0,0
k,1001,950,0,0
AROTAT,all, , , , , , 1000, 1001,360,36,
LESIZE,ALL, , ,1, ,1, , ,1,
MSHKEY,0
mshape,0,3d
amesh,all
|