混凝土非線性分析-應(yīng)力應(yīng)變曲線輸出概要

1、各位同學(xué)，這里幾乎已經(jīng)是最完整的命令流了，關(guān)于這方面工作，初步建議如下：（1）初步的模擬最好先不要考慮混凝土下降段，不要考慮混凝土的壓碎，同時(shí)單元尺寸不 要過小，以免影響收斂。（2 ）反復(fù)試算結(jié)果表明，混凝土的MISO 和 KINH 模型、 DP 和 MISES 準(zhǔn)則、剪力傳遞系數(shù)對(duì)模擬結(jié)果影響均不大，但在一定程度上影響了收斂性的好壞。（3）當(dāng)解決了收斂問題之后，最大的難點(diǎn)是下降段，對(duì)于一個(gè)橋墩（柱）的反復(fù)荷載試驗(yàn)， 引起承載力下降的原因大體可歸結(jié)為： 塑性鉸區(qū)混凝土壓碎脫落、 縱筋屈曲， 對(duì)于以剪切破 壞為主的試件， 還包括剪切滑移的影響。 對(duì)于以黏結(jié)破壞為主的試件， 還包括縱筋混凝土 黏結(jié)

2、滑移破壞。因此，能否正確模擬這些將是模型正確與否的關(guān)鍵。眾所周知，以 link8 單 元模擬鋼筋自然無法模擬縱筋屈曲，這樣就把一個(gè)因素忽略了。（4）作者們最后采用了擴(kuò)大的破壞面，坦率的講，這也是為了模擬混凝土的壓碎破壞而采 取的沒有辦法的辦法， 鑒于作者水平有限， 至今也無法對(duì)此做出更為合理的解釋。 破壞面到 底擴(kuò)大多少，仍無規(guī)律可循。且一旦打開壓碎開關(guān)，會(huì)出現(xiàn)諸如難以收斂、 “假壓碎 ”等一系 列問題。所以我們建議大家起初還是先不要考慮壓碎為好。定義應(yīng)力應(yīng)變曲線1， 定義變量：拾取主菜單： Main Menu>Time Hist postproc>Define V

3、ariables> 在隨之彈出的對(duì)話 框中點(diǎn)擊 Add 鍵，定義第一個(gè)變量序號(hào)為 2，選取第一個(gè)變量 stress, 確定與之對(duì)應(yīng)的下一 級(jí)選項(xiàng)（如 Y-direction SY 等）；返回定義變量對(duì)話框，再點(diǎn)擊 add 鍵，定義第二個(gè)變量序 號(hào)為 3，選取第二個(gè)變量 strain-elastic 及以及對(duì)贏得下一級(jí)選項(xiàng)（如 Y-dir'n EPEL Y 等，在 應(yīng)力 - 應(yīng)變圖中，其向量的取向應(yīng)相同）。同理再定義變量4，選取變量 strain-plastic 及與之對(duì)應(yīng)的下一級(jí)選項(xiàng)如 Y-dir'n EPEL Y 等），在應(yīng)力 - 應(yīng)變圖中，應(yīng)變是彈性應(yīng)變和塑性

4、應(yīng) 變累加的總應(yīng)變。為使其實(shí)現(xiàn)相加，還需進(jìn)行以下操作：拾取主菜單： Main Menu>Time Hist postproc>math operation>add,定義計(jì)算變量序號(hào)為 5，同時(shí)在相應(yīng)交互框內(nèi)輸入 3 和 4。點(diǎn)擊確認(rèn)鍵，則由變量 3，4代表的應(yīng)變之和就存在變量 5 中。2，繪制應(yīng)力 -應(yīng)變曲線：拾取主菜單： Main Menu>Time Hist postproc>setting>graph.設(shè)置 x 軸向變量為單變量，并將其變量序號(hào)定義為 5。點(diǎn)擊確定鍵退出退化框。 拾取應(yīng)用菜單： Utili

5、ty Menu>plot ctrls>styles>Graphs>Modify axis.將 x,y 坐標(biāo)軸分別命名為 Y-strains,Y-stress, 拾取主菜單： Main Menu>Time Hist postproc>graph variables. 在對(duì)話框上 "the first variable" 對(duì)應(yīng)的交互框中輸入 2。點(diǎn)擊確定鍵，則預(yù)想的應(yīng)力 -應(yīng)變曲線就顯示在屏幕上。 ok! 試試看！歡迎大家繼續(xù)批評(píng)指正！/prep7et,1,solid65et,2,link8et

6、,3,solid45r,1r,2,236e-6r,3,75.4e-6r,4,151.2e-6!混凝土本構(gòu)關(guān)系mp,ex,1,2.522e10mp,prxy,1,0.2mp,dens,1,2500tb,concr ,1,1 tbdata,0.5,0.95,5e6,-1tb,miso,1,1,7 tbtemp,0 tbpt,0.00068,17.15e6 tbpt,0.00136,28.03e6 tbpt,0.00204,33.62e6 tbpt,0.00272,35.96e6 tbpt,0.0034,36.53e6 tbpt,0.00408,36.53e6!縱筋的本構(gòu)關(guān)系mp,ex,2,2.0e

7、11mp,prxy,2,0.28tb,bkin,2 tbtemp,0 tbdata,362e6,2e9！箍筋的本構(gòu)關(guān)系mp,EX,3,2.1E11mp,PRXY,3,.27tb,bkin,3 tbtemp,0 tbdata,272.5e6,2.1e9！建立局部坐標(biāo)系，采用柱坐標(biāo)，繞 y 軸順時(shí)針轉(zhuǎn) 90 度！復(fù)制節(jié)點(diǎn)， 3 次，編號(hào)增量，源節(jié)點(diǎn) 18 ， 0,0,0.3!建立模型 local,11,1,0,0,0,0,-90,0 csys,11*do,j,1,8n,j,0.185,j*360/8,-0.6 *enddo ngen,3,100,1,8,1,0,0,0.3ngen,4,100,20

8、1,208,1,0,0,0.15 ngen,7,100,501,508,1,0,0,0.3type,2real,3mat,3*do,i,1,6*do,j,1,7e,(i-1)*100+j,(i-1)*100+j+1 *enddoe,(i-1)*100+1,(i-1)*100+8*enddotype,2real,4mat,3*do,i,1,5*do,j,1,7 e,(i+5)*100+j,(i+5)*100+j+1 *enddoe,(i+5)*100+1,(i+5)*100+8*enddotype,2real,2mat,2*do,i,1,11*do,j,1,8e,(i-1)*100+j,i*10

9、0+j*enddo*enddoallsel,allnummrg,allnumcmp,all！指定單元屬性，建立基礎(chǔ)箍筋單元！建立柱上箍筋單元！建立縱向鋼筋！合并節(jié)點(diǎn)，壓縮節(jié)點(diǎn)wpcsys,1,11 cylind,0,0.185,-0.6,0,0,360 cylind,0,0.185,0,0.45,0,360 cylind,0.185,0.3,-0.6,0,0,360 cylind,0,0.185,0.45,2.25,0,360 cylind,0,0.185,2.25,2.55,0,360 cylind,0.185,0.20,2.25,2.55,0,360 cylind,0,0.185,2.55

10、,2.565,0,360！建立圓柱，內(nèi)外徑 i-or ，高度 z1-z2 ，起終角度wprota,90 vsel,s,1,7,1 vsbw,all numcmp,all wprota,90 vsel,s,1,14,1 vsbw,all numcmp,all！繞 y 軸逆時(shí)針轉(zhuǎn) 90 度！選擇體，選擇類型， item ， comb ，最小編號(hào)，最大編號(hào)，增量 ！體相減！繞 z 軸逆時(shí)針轉(zhuǎn) 90 度vsel,s,1,28,1vglue,allnumcmp,all wprota,-90 wprota,-90allsel,allvsel,s,1,16,1 vsel,a,21,24,1 vatt,1,1

11、,1vsel,s,17,20,1 vsel,a,25,28,1 vatt,2,1,3allsel,alllsel,s,1,8,1 lsel,a,11,22,1 lsel,a,25,40,1 lsel,a,43,46,1 lesize,all,2！選擇線!補(bǔ)選（a）線！線尺寸為 2lsel,s,47,50,1lsel,a,53,56,1lsel,a,58,73,1lesize,all,2lsel,s,98,107,1lsel,a,110,112,2lesize,all,1lsel,s,23,24,1lsel,a,9,10,1lesize,all,2 lsel,s,78,82,1lesize,al

12、l,3lsel,s,83,87,1lesize,all,6lsel,s,41,42,1 lsel,a,88,92,1 lesize,all,2lsel,s,93,97,1lesize,all,1mshkey,1 mshape,0,3d vsel,s,1,28 vmesh,allallsel,allnummrg,all numcmp,all wprota,90 csys,0 allsel,all/solu！進(jìn)入加載、求解csys,11 nsel,s,loc,x,0.3 d,all,all allsel,all！選擇半徑在 x=3 處的節(jié)點(diǎn)，并施加約束OUTRES,ALL,2！輸出獲得的解csys

13、,0pred,on cnvtol,f,0.05,2 autots,1 lnsrch,1 ncnv,2！打開預(yù)應(yīng)力選項(xiàng)! 定義收斂條件 ,使用缺省的 VALUE !打開自動(dòng)時(shí)間步控制 !打開線性搜索!如果不收斂時(shí)結(jié)束而不退出neqit,50!每一子步中方程的迭代次數(shù)限值savetime,1 nsubst,10 asel,s,29,32,1 sfa,all,pres,5468699 allsel,all acel,9.8kbc,1lswrite,1！定義第 1 載荷步！在已選擇的面上施加壓力荷載！重力加速度 y=9.8！階躍加載方式time,10cp,1,ux,298,300,303d,298,u

14、x,0.015！定義第 2 荷載步!將節(jié)點(diǎn) 298,300,303ux 向位移耦合！再借點(diǎn) 298 處施加 ux 位移 0.015allsel,allKBC,0allsel,allNSUBST,20,1000000,10lswrite,2time,20 ！定義第 3 載荷步 d,298,ux,-0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,3time,30 ！定義第 4 載荷步 d,298,ux,0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,4time

15、,40 ！定義第 5 載荷步 d,298,ux,-0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,5time,50d,298,ux,0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,6time,60d,298,ux,-0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,7！定義第 6 載荷步！定義第 7 載荷步time,70d,298,ux,0.015allsel,allKBC,0allsel

16、,allNSUBST,100,1000000,10lswrite,8time,80d,298,ux,-0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,9time,90d,298,ux,0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,10！定義第 8 載荷步time,100d,298,ux,-0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,11time,110d,298,ux,0.03al

17、lsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,12time,120d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,13time,130d,298,ux,0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,14time,140d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,15time,150d,298

18、,ux,0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,16time,160d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,17time,170d,298,ux,0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,18time,180d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,19tim

19、e,190d,298,ux,0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,20time,200d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10lswrite,21 time,210 d,298,ux,0.045allsel,allKBC,0allsel,allNSUBST,300,1000000,10lswrite,22lssolve,1,22,1命令流：1、單調(diào)水平荷載作用下的命令流：建模：!mono brick wall!solid65- 整體式建模

20、!*Finish/clear/units,si/com,Structural/config,nres,400000/prep7/title,mono-brick-wall!*depth_sup=0.24 ! 彈性支墊 240mm!單元屬性!for brick, mortar and leveling coatet,1,solid65!*KEYOPT,1,1,0KEYOPT,1,5,0KEYOPT,1,6,0KEYOPT,1,7,1!for elastic supportet,2,solid45!* 砌體材料 !MU10,M5 砌體材料屬性mp,prxy,1,0.15mp,dens,1,1700

21、mp,ex,1,0.24e10!砌體屈服準(zhǔn)則TB,MKIN,1 TBTEMP,strain TBdata,1,0.2E-3,1.0E-3,2.0E-3,3.0E-3,4.8E-3TBTEMP, !TBDATA,1,0.48E6,0.82E6,1.32E6,1.5E6,1.32E6TBDATA,1,0.48E6,1.31E6,2.11E6,2.4E6,2.11E6 !/XRANGE,0,0.02!TBPLOT,MKIN,1 !砌體破壞準(zhǔn)則 !hntrl=0.13E6 hntrl=0.21E6 tb,concr,1,4, tbdata,0.2,0.9,hntrl,-1 ! 不考慮壓碎!* 彈性支撐

22、 !彈性支撐暫取鋼筋材料屬性mp,ex,2,2.0e11 mp,dens,2,7.600 mp,prxy,2,0.30 !real constantsr,1!*modeling!geometry modelblock,0,3.12,0,0.24,0,3.12 block,-0.24,3.12,0,0.24,3.12,3.36 Vsel,s,volu,1,2,1$aslv,s$lsla,s lesize,all,0.24Vsel,s,volu,1$vatt,1,1,1, !brick wallVsel,s,volu,2$vatt,2,1,2, !brick wall vsel,all ! 必須再

23、次選擇所有單元VSWEEP,ALL save numcmp,all allsel!*/device,vector,1/eshape,1!* 耦合 allsel,all!elastic support nodes set esel,s,type,2$nsle,s,all$nsel,r,loc,z,3.12 nsel,u,loc,x,-0.24$cm,slavenode,node!other node set esel,s,type,1$nsle,s,all$cm,masternode,node /input,ucouple,macsave finish耦合文件 ucouple,mac ：!dow

24、hile,kk k=nnear(slave_node(i) nsel,s,cp,all kk=nsel(k)+0.001 開始 allsel ! 最好保留這句命令!* 將從屬節(jié)點(diǎn)編號(hào)依次存入數(shù)組 * cmsel,s,slavenode *get,count1_node,node,0,count*del,slave_node*dim,slave_node,array,count1_node*get,slave_node(1),node,0,num,min*do,i,2,count1_node slave_node(i)=ndnext(slave_node(i-1) *enddo將主節(jié)點(diǎn)編號(hào)依次存

25、入數(shù)組 *allselcmsel,s,masternode*get,count2_node,node,0,count*del,master_node*dim,master_node,array,count2_node *get,master_node(1),node,0,num,min*do,i,2,count2_node master_node(i)=ndnext(master_node(i-1) *enddo!* 將與從屬節(jié)點(diǎn)耦合的節(jié)點(diǎn)數(shù)組初始化 * *del,cp_node*dim,cp_node,array,count1_node*do,i,1,count1_nodecp_node(i

26、)=0 *enddoI*開始選擇程序*allsel cmsel,s,masternode*do,i,1,count1_node kk=1 k=1allsel cmsel,s,masternode nsel,u,node,k cm,masternode,node *enddocp_node(i)=k*enddo選擇完畢 *開始耦合 *allsel,all/prep7*do,i,1,count1_node cp,next,UX,slave_node(i),cp_node(i) cp,next,UY,slave_node(i),cp_node(i) cp,next,UZ,slave_node(i),

27、cp_node(i) *enddo耦合完畢 *求解方法 NR 法： brick wall NR solution !displacement load !用 NR 法!*solution/soluantype,0!constraintnsel,s,loc,z,0 d,all,all!* esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uy!*!esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24!d,all,uz lumpm,on solcontrol,on pred,on lnsrch,on autots,

28、1 ! 打開自動(dòng)時(shí)間步控制ncnv,2 ! 如果不收斂時(shí)結(jié)束而不退出compression=0.24E7!*load steps!*1stacel,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.119esln,s,0,allsfe,all,6,pres,0,0.03*compressionallsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,2,200,1neqit,15KBC,0nlgeom,1cnvtol,f,0.01,2,0.05allsellswrite,1!*2nd stepnsel,s,loc,x,-0.24

29、!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3d,all,ux,1.0E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,0.05,2,0.05allsellswrite,2!*3rd stepnsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3d,all,ux,2.5E-4 allsel,all

30、 !*OUTRES,ERASE OUTRES,ALL,ALL NSUBST,200,2000,100 neqit,15KBC,0 nlgeom,1 cnvtol,f,0.09,2,0.25 allsel lswrite,3!*4th stepnsel,s,loc,x,-0.24 !esln,s,0,all !sfe,all,4,pres,0,60*0.0593E6 !d,all,ux,2.79E-3 d,all,ux,4.0E-3 allsel,all!* OUTRES,ERASE OUTRES,ALL,ALL NSUBST,200,2000,100 neqit,15KBC,0 nlgeom,

31、1 cnvtol,f,0.1,2,0.5 allsel lswrite,4 !* allsel lssolve,1,4,1 FINISH也可以用弧長法，不過沒有一定經(jīng)驗(yàn)這比較難控制：!mono brick wall-1 arclength solution!displacement load!*solution/soluantype,0!constraintnsel,s,loc,z,0d,all,all!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uy!* lumpm,on solcontrol,on ncnv,2 ! 如果不收斂時(shí)結(jié)束

32、而不退出 compression=0.24E7!*load steps!*first stepacel,9.8 nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.12 esln,s,0,all sfe,all,6,pres,0,0.03*compression allsel,all!*OUTRES,ERASE OUTRES,ALL,ALL arclen,on,1E-5, ! 使用弧長法 NSUBST,200 neqit,15KBC,0 nlgeom,1 cnvtol,f,0.1,2,0.5 allsel,all lswrite,1!*second stepesel,

33、s,type,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uz!* nsel,s,loc,x,-0.24 !esln,s,0,all !sfe,all,4,pres,0,60*0.0593E6 d,all,ux,2.50E-3 !d,all,ux,3.51E-3 allsel,all!* OUTRES,ERASE OUTRES,ALL,ALL arclen,on,1E-5, ! 使用弧長法 kbc,0NSUBST,250 !arctrm,u,0.0028,313,ux neqit,15nlgeom,1!cnvtol,f,0.005,2,0.01!cnvtol,f

34、,0.05,2,0.5 cnvtol,f,0.1,2,0.5 allsel lswrite,2!*third stepesel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uz!*nsel,s,loc,x,-0.24 !esln,s,0,all !sfe,all,4,pres,0,60*0.0593E6 !d,all,ux,2.79E-3!d,all,ux,3.18E-3 d,all,ux,3.00E-3 allsel,all!*OUTRES,ERASE OUTRES,ALL,ALL arclen,on,1E-7, ! 使用弧長法 kbc,0NSU

35、BST,2000 !arctrm,u,0.0028,313,ux neqit,15 nlgeom,1!cnvtol,f,0.005,2,0.01!cnvtol,f,0.05,2,0.5 cnvtol,f,0.5,2,0.8 allsel lswrite,3 allsel!*fourth stepesel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6 d,all,ux,3.51E-3 allsel,all!*OUTRES,E

36、RASE OUTRES,ALL,ALL arclen,on,1E-7, ! 使用弧長法kbc,0NSUBST,300 !arctrm,u,0.0028,313,ux neqit,15nlgeom,1 !cnvtol,f,0.005,2,0.01 !cnvtol,f,0.05,2,0.5 cnvtol,f,0.6,2,0.9 allsellswrite,4allsel!*lssolve,1,4,1FINISH后處理：/post26 esel,s,type,1 eplot lines,5000 nsol,2,326,u,x,ux !nsol,2,313,u,x,ux *do,i,1,14 jj=n

37、ode(i-1)*0.24,0.24,0) kk=node(i-1)*0.24,0,0)! ll=node(225.5+(i-1)*225.5,3300,0) rforce,3,jj,f,x,fx add,6,3,6rforce,4,kk,f,x,fxadd,6,4,6! rforce,5,ll,f,x,fx!add,6,5,6*enddo prod,7,2,1000 prod,8,6,-1/1000 /axlab,x,Displacement(mm) /axlab,y,F(kN) /XRANGE,0,4 /YRANGE,0,100xvar,7plvar,8!rforce,5,340,f,x,

38、fx!esol,5,87,221,f,x,fx2、循環(huán)水平荷載作用下的命令流： 出求解文件不同外，其他都相同； !cyclic NR solution !displacement load!*solution/soluantype,0!constraintnsel,s,loc,z,0 d,all,all!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uy !*esel,s,type,2$nsle,s,all$nsel,r,loc,x,3.36d,all,uy !* lumpm,on solcontrol,on pred,on lnsrch,

39、on autots,1 ! 打開自動(dòng)時(shí)間步控制 ncnv,2 ! 如果不收斂時(shí)結(jié)束而不退出 compression=0.24E7!*load steps!*1st stepacel,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.119 esln,s,0,allsfe,all,6,pres,0,0.03*compression allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,10neqit,15KBC,0 nlgeom,1 cnvtol,f,0.1,2,0.5 allsel,all lswrite,

40、1!* 初裂荷載 20%做 1 次循環(huán)!* 用 NR 法 !*2nd stepesel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uz!* esel,s,type,2$nsle,s,all$nsel,r,loc,x,3.36 d,all,uz!* nsel,s,loc,x,-0.24 !esln,s,0,all !sfe,all,4,pres,0,60*0.0593E6 d,all,ux,1.48E-4 allsel,all!*OUTRES,ERASE OUTRES,ALL,ALL NSUBST,20,2000,10 neqit,15NROPT

41、,FULL,OFFKBC,0 nlgeom,1 cnvtol,f,0.1,2,0.5 allsel lswrite,2!*3rd stepnsel,s,loc,x,-0.24 ddele,all,uxnsel,s,loc,x,3.36 d,all,ux,-1.48E-4 allsel,all!*OUTRES,ERASE OUTRES,ALL,ALL NSUBST,40,2000,20 neqit,15NROPT ,FULL,OFFKBC,0 nlgeom,1cnvtol,f,0.1,2,0.5 allsel lswrite,3!*4th stepnsel,s,loc,x,3.36 ddele,

42、all,uxnsel,s,loc,x,-0.24 d,all,ux,1.48E-4 allsel,all!*OUTRES,ERASE OUTRES,ALL,ALL NSUBST,40,2000,20 neqit,15NROPT ,FULL,OFFKBC,0 nlgeom,1 cnvtol,f,0.1,2,0.5 allsel lswrite,4!*one cycle!* 初裂荷載 50%做 1 次循環(huán)!* 用 NR 法 !*5th stepnsel,s,loc,x,-0.24 ddele,all,uxnsel,s,loc,x,-0.24 d,all,ux,2.85E-4 allsel,all!

43、*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,100,2000,10 neqit,15NROPT ,FULL,OFFKBC,0 nlgeom,1 cnvtol,f,0.1,2,0.5 allsel lswrite,5!*6th stepnsel,s,loc,x,-0.24 ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-2.85E-4 allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,30 neqit,15NROPT ,FULL,OFFKBC,0 nlgeom,1 cnvtol,f,0

44、.1,2,0.5 allsel lswrite,6!*7th stepnsel,s,loc,x,3.36 ddele,all,uxnsel,s,loc,x,-0.24 d,all,ux,2.85E-4 allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,30 neqit,15NROPT ,FULL,OFFKBC,0 nlgeom,1 cnvtol,f,0.1,2,0.5 allsel lswrite,7!*two cycle!* 初裂荷載 80%做 1 次循環(huán)!* 用 NR 法 !*8th stepnsel,s,loc,x,-0.24 d

45、dele,all,uxnsel,s,loc,x,-0.24d,all,ux,3.5E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, ! 使用弧長法 kbc,0NSUBST,100,2000,20 neqit,15KBC,0 nlgeom,1 cnvtol,f,0.5,2,0.9 allsel lswrite,8!*9th stepnsel,s,loc,x,-0.24 ddele,all,uxnsel,s,loc,x,3.36 d,all,ux,-3.5E-4 allsel,all!*OUTRES,ERASE OUTRES,ALL

46、,ALL !arclen,on,1E-7, ! 使用弧長法 kbc,0NSUBST,200,2000,30 neqit,15KBC,0 nlgeom,1 cnvtol,f,0.5,2,0.9 allsel lswrite,9!*10th step-nsel,s,loc,x,3.36 ddele,all,uxnsel,s,loc,x,-0.24 d,all,ux,3.5E-4 allsel,all!*OUTRES,ERASE OUTRES,ALL,ALL !arclen,on,1E-7, ! 使用弧長法 kbc,0NSUBST,200,2000,30 neqit,15KBC,0 nlgeom,1

47、 cnvtol,f,0.5,2,0.9 allsel lswrite,10!*three cycle!* 初裂荷載 100% 做 1 次循環(huán)!* 用 NR 法 !*11th stepnsel,s,loc,x,-0.24 ddele,all,uxnsel,s,loc,x,-0.24 d,all,ux,4.42E-4 allsel,all!*OUTRES,ERASE OUTRES,ALL,ALL !arclen,on,1E-7, ! 使用弧長法 kbc,0NSUBST,100,2000,20 neqit,15KBC,0 nlgeom,1 cnvtol,f,0.5,2,0.9 allsel lswr

48、ite,11!*12th stepnsel,s,loc,x,-0.24 ddele,all,uxnsel,s,loc,x,3.36 d,all,ux,-4.42E-4 allsel,all!*OUTRES,ERASE OUTRES,ALL,ALL !arclen,on,1E-7, ! 使用弧長法 kbc,0NSUBST,200,2000,30 neqit,15KBC,0 nlgeom,1 cnvtol,f,0.5,2,0.9 allsel lswrite,12!*13th stepnsel,s,loc,x,3.36 ddele,all,uxnsel,s,loc,x,-0.24 d,all,ux

49、,4.42E-4 allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, ! 使用弧長法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,13!*four cycle!*61.4KN 的荷載做 1 次循環(huán) !* 用 NR 法 !*14th stepnsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,6.72E-4allsel,all!*OUTRES,ERASEOUTRES,ALL

50、,ALL!arclen,on,1E-7, ! 使用弧長法kbc,0NSUBST,100,2000,20neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,14!*15th stepnsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-6.72E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, ! 使用弧長法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,0.5,

51、2,0.9allsellswrite,15!*16th stepnsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,6.72E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, ! 使用弧長法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,16!*five cycle!*74.1KN 的荷載做 1 次循環(huán) !* 用 NR 法 !*17th stepnsel,s,loc,x

52、,-0.24 ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,9.01E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, ! 使用弧長法kbc,0NSUBST,100,2000,20neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,17!*18th step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-9.01E-4allsel,all!*OUTRES,ERASEOUTRES,

53、ALL,ALL!arclen,on,1E-7, ! 使用弧長法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,18!*19th step-nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,9.01E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, ! 使用弧長法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9