資源描述:
《ansys鄧肯-張材料模型》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、ANSYS鄧肯-張材料模型樓主給的在ANSYS上實(shí)現(xiàn)鄧肯-張模型的方法很有用,但其中還有幾點(diǎn)需要修正的,這也是樓上的兄弟們有疑問的原因。我把樓主的代碼運(yùn)行了一下,然后對照作了修改,現(xiàn)在上傳一下,有問題的兄弟可以仔細(xì)對照一下,在這里我對其中幾個比較明顯的問題說明一下:1.MP命令不能直接給單元加材料,這是對的。在這里,樓主遺漏了一下命令:MPCHG,具體見下面的修改過的代碼。2.關(guān)于密度的問題。這些要在宏中定義,每修改一種材料(即調(diào)用一次鄧肯-張子程序)就要修改一次材料的密度,其他有關(guān)材料的問題可以類推。3.關(guān)于施加重力的問題。要在調(diào)用宏后,在同一個循環(huán)
2、中重新定義一下重力。以下是我修改過的樓主的代碼,希望對兄弟們有所幫助。!用APDL得到初步成果,貼于此供感興趣的朋友參考,不當(dāng)之處敬請指正,!歡迎加以完善。!基本思路:!鄧肯-張模型的關(guān)鍵點(diǎn)是材料的彈性模量隨大小主應(yīng)力差!及小主應(yīng)力(圍壓)的變化而變化,用APDL實(shí)現(xiàn)之的基本思路是:!給每個單元定義一個材料號,分級施加荷載,在每個荷載步結(jié)束時提取出各!單元的大小主應(yīng)力,據(jù)此計(jì)算出下個荷載步的彈性模量Et,修改各單元之MP,!用于下一步計(jì)算。!以下是一個簡單算例,copy出去可直接運(yùn)行。!!!常規(guī)三軸試驗(yàn)?zāi)M!**********************
3、************************************FINISH/CLEAR/TITLE,NumericalSimulationofthreeaxestestingofsoils/PREP7*dim,SUy,array,50!Settlementrecords*dim,MaxPs,array,120!Maxhistoryp1-p3*dim,MaxDs,array,120!MaxhistoryDs!*dim,EEt,array,50!Etofelememt!!!Duncan-ChangModel!!!Symbols:c-粘滯力,F(xiàn)ai
4、-內(nèi)摩擦角,Sf-破壞強(qiáng)度(p1-p3)f,!Ds-應(yīng)力水平,Pa-大氣壓,P3-圍壓!***********************************************************************CREATE,Duncan-Chang!CreatMacrofile*afun,deg!Unitofangle*set,Pa,1e5*set,P1,-ArrS3(i)!注意:巖土工程中應(yīng)力為拉負(fù)壓正*set,P3,-ArrS1(i)*if,P3,LT,0.1*Pa,thenP3=0.1*Pa !圍壓最小取值*endifSf
5、0=2*(c0*cos(Fai)+P3*sin(Fai))/(1-sin(Fai))!Mohr-Coulomb破壞強(qiáng)度(p1-p3)fDs=(P1-P3)/Sf0!應(yīng)力水平,*if,Ds,GT,0.95,thenDs=0.95 !應(yīng)力水平最大取值*endif!判斷加卸荷,如果(P1-P3)小于歷史最大值視為卸荷-再加荷過程*if,MaxPs(i),LT,P1-P3,thenEi=k0*Pa*(P3/Pa)**n0Et0=Ei*(1-Rf*Ds)**2 !加荷情況的切線模量MaxPs(i)=P1-P3!保存歷史最大應(yīng)力*elseif,MaxPs(i
6、),GE,P1-P3Et0=Kur*Pa*(P3/Pa)**n0 !卸荷模量*endifmp,ex,i,Et0!修改材料i的Etmp,nuxy,i,Mu0mp,dens,I,1800!重新定義第I種材料的密度*END!***********************單元類型***********************************et,1,42 !平面四節(jié)點(diǎn)單元KEYOPT,1,3,2 !平面應(yīng)變!*********************以下定義材料初始模量*************************mp,ex,1,3.728e7
7、!砂土的彈性模量mp,nuxy,1,0.33mp,dens,1,1800!******************************************************************blc4,0,0,0.08,0.15 !建立幾何模型/PNUM,AREA,1/REPLOTaesize,all,0.01!網(wǎng)格劃分mat,1amesh,allnsel,s,loc,y,0 !邊界條件d,all,Uy!底邊界豎向約束nsel,s,loc,x,0d,all,Ux!左側(cè)邊界水平向約束nsel,all/replotfini/SOLUtime,
8、0.01!施加圍壓ACEL,0,9.8,0sfl,all,pres,2e5!200kPasol