本文介紹了Ansys APDL語(yǔ)言方式的非線性方程求解相關(guān)內(nèi)容。

　　在學(xué)習(xí)編寫ABAQUS的UMAT子程序時(shí)，在受拉部分考慮了混凝土損傷的影響，需要確定屈服應(yīng)力和屈服應(yīng)變。這個(gè)方程簡(jiǎn)化后屬于非線性方程，格式如：A*(x-1)^1.7+x-B=0

　　為了驗(yàn)算方程求解的正確性，使用APDL來求這個(gè)方程的解。假定A=B=2，分別給定X初始，并給出迭代次數(shù);

　　x初始值為2時(shí)，迭代5次;初始值為1，迭代6次;初始值為1000，迭代13次;

　　結(jié)果均無限逼近1.461942.

　　下面給出APDL的命令，有興趣的朋友也可以試一下。(注意X的初始不要取負(fù)值，這是方程成立的前提)

　　A=2

　　B=2

　　x0=1

　　*cfopen,aaaaa,txt

　　*do,i,1,1000

　　xn=x0-(A*(x0-1)**1.7+x0-B)/(1.7*A*(x0-1)**0.7+1)

　　i=i+1

　　*if,abs((xn-x0)/x0),gt,1.0e-6,then

　　x0=xn

　　*else

　　*vwrite(*,*) xn,i

　　(G15.5,I3)

　　*cfclose

　　*exit

　　*endif

　　*enddo

　　簡(jiǎn)單求解非線性，APDL確實(shí)好用，比excel更具公式化，便于檢查;相對(duì)一般的編程軟件，又方便，不易出錯(cuò)。

　　1.與求解一相比，直接使用條件循環(huán)，*dowhile實(shí)現(xiàn)。

　　2.該命令流的使用：通過界面/input輸入。

　　ft=1.43

　　E0=3e4

　　E1=ft/(ft**0.54*65e-6)

　　At=0.312*ft**2

　　A=E0/E1-1

　　B=A+1

　　x0=1.0

　　D=x0

　　c=1

　　i=0

　　*cfopen,aaaaa,txt

　　*dowhile,c

　　TERM1=1/(At*(x0-1)**1.7+x0)

　　TERM2=1.7*At*(x0-1)**0.7+1

　　TERM3=x0*B-A*x0*TERM1-D

　　TERM4=B-A*(TERM1-x0*TERM1**2*TERM2)

　　x=x0-TERM3/TERM4

　　*vwrite,x,x,i

　　(2(F15.5,","),F15.0)

　　*if,abs(x-x0)/x0,lt,1e-6,then

　　c=0

　　*else

　　c=1

　　*endif

　　x0=x

　　i=i+1

　　*enddo

　　*cfclose