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

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

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

3、dd鍵，定義第一個(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)變和塑性應(yīng)變累加的總應(yīng)變。為使其實(shí)現(xiàn)相加，還需進(jìn)行以下操作：拾取主菜單：Main Menu>Time Hi

4、st 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)用菜單：Utility Menu>plot ctrls>styles>Graphs>Modify ax

5、is.將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)指正！/prep7 et,1,solid65 et,2,link8 et,3,solid45r,1r,2,236e-6r,3,75.4e-6 r,4,151.2e-6 !混凝土本構(gòu)關(guān)系mp,ex,1,2.522e10mp,pr

6、xy,1,0.2 mp,dens,1,2500tb,concr,1,1 tbdata,0.5,0.95,5e6,-1tb,miso,1,1,7tbtemp,0 tbpt,0.00068,17.15e6tbpt,0.00136,28.03e6tbpt,0.00204,33.62e6tbpt,0.00272,35.96e6tbpt,0.0034,36.53e6tbpt,0.00408,36.53e6!縱筋的本構(gòu)關(guān)系mp,ex,2,2.0e11mp,prxy,2,0.28 tb,bkin,2 tbtemp,0 tbdata,362e6,2e9！箍筋的本構(gòu)關(guān)系mp,EX,3,2.1E11mp,PRXY

7、,3,.27tb,bkin,3tbtemp,0tbdata,272.5e6,2.1e9 !建立模型local,11,1,0,0,0,0,-90,0 ！建立局部坐標(biāo)系，采用柱坐標(biāo)，繞y軸順時(shí)針轉(zhuǎn)90度csys,11*do,j,1,8n,j,0.185,j*360/8,-0.6*enddo ngen,3,100,1,8,1,0,0,0.3 ！復(fù)制節(jié)點(diǎn)，3次，編號(hào)增量，源節(jié)點(diǎn)18，0,0,0.3ngen,4,100,201,208,1,0,0,0.15ngen,7,100,501,508,1,0,0,0.3type,2 ！指定單元屬性，建立基礎(chǔ)箍筋單元real,3mat,3*do,i,1,6*do,

8、j,1,7e,(i-1)*100+j,(i-1)*100+j+1*enddoe,(i-1)*100+1,(i-1)*100+8*enddotype,2 ！建立柱上箍筋單元real,4mat,3*do,i,1,5*do,j,1,7e,(i+5)*100+j,(i+5)*100+j+1*enddoe,(i+5)*100+1,(i+5)*100+8*enddotype,2 ！建立縱向鋼筋real,2mat,2*do,i,1,11*do,j,1,8e,(i-1)*100+j,i*100+j*enddo*enddoallsel,all ！合并節(jié)點(diǎn)，壓縮節(jié)點(diǎn)nummrg,allnumcmp,allwpcs

9、ys,1,11cylind,0,0.185,-0.6,0,0,360 ！建立圓柱，內(nèi)外徑i-or，高度z1-z2，起終角度cylind,0,0.185,0,0.45,0,360cylind,0.185,0.3,-0.6,0,0,360cylind,0,0.185,0.45,2.25,0,360cylind,0,0.185,2.25,2.55,0,360cylind,0.185,0.20,2.25,2.55,0,360cylind,0,0.185,2.55,2.565,0,360wprota,90 ！繞y軸逆時(shí)針轉(zhuǎn)90度vsel,s,1,7,1 ！選擇體，選擇類型，item，comb，最小編號(hào)，

10、最大編號(hào)，增量vsbw,all ！體相減numcmp,all wprota,90 ！繞z軸逆時(shí)針轉(zhuǎn)90度vsel,s,1,14,1 vsbw,all numcmp,all vsel,s,1,28,1 vglue,all numcmp,all wprota,-90 wprota,-90 allsel,allvsel,s,1,16,1vsel,a,21,24,1vatt,1,1,1 vsel,s,17,20,1vsel,a,25,28,1vatt,2,1,3 allsel,all lsel,s,1,8,1 ！選擇線lsel,a,11,22,1 ！補(bǔ)選（a）線lsel,a,25,40,1 lsel,

11、a,43,46,1 lesize,all,2 ！線尺寸為2lsel,s,47,50,1 lsel,a,53,56,1lsel,a,58,73,1lesize,all,2lsel,s,98,107,1 lsel,a,110,112,2lesize,all,1lsel,s,23,24,1 lsel,a,9,10,1lesize,all,2lsel,s,78,82,1 lesize,all,3lsel,s,83,87,1 lesize,all,6lsel,s,41,42,1 lsel,a,88,92,1lesize,all,2lsel,s,93,97,1 lesize,all,1mshkey,1 m

12、shape,0,3d vsel,s,1,28vmesh,all allsel,all nummrg,allnumcmp,allwprota,90 csys,0allsel,all /solu ！進(jìn)入加載、求解csys,11nsel,s,loc,x,0.3 ！選擇半徑在x=3處的節(jié)點(diǎn)，并施加約束d,all,allallsel,all OUTRES,ALL,2 ！輸出獲得的解csys,0pred,on ！打開(kāi)預(yù)應(yīng)力選項(xiàng)cnvtol,f,0.05,2 !定義收斂條件,使用缺省的VALUE autots,1 !打開(kāi)自動(dòng)時(shí)間步控制 lnsrch,1 !打開(kāi)線性搜索 ncnv,2 !如果不收斂時(shí)結(jié)束而不退

13、出 neqit,50 !每一子步中方程的迭代次數(shù)限值 savetime,1 ！定義第1載荷步nsubst,10asel,s,29,32,1sfa,all,pres,5468699 ！在已選擇的面上施加壓力荷載allsel,allacel,9.8 ！重力加速度y=9.8kbc,1 ！階躍加載方式lswrite,1 time,10 ！定義第2荷載步cp,1,ux,298,300,303 !將節(jié)點(diǎn)298,300,303ux向位移耦合d,298,ux,0.015 ！再借點(diǎn)298處施加ux位移0.015allsel,all KBC,0allsel,all NSUBST,20,1000000,10 lsw

14、rite,2time,20 ！定義第3載荷步d,298,ux,-0.015allsel,all KBC,0allsel,all NSUBST,100,1000000,10 lswrite,3time,30 ！定義第4載荷步d,298,ux,0.015allsel,all KBC,0allsel,all NSUBST,100,1000000,10 lswrite,4time,40 ！定義第5載荷步d,298,ux,-0.015 allsel,all KBC,0allsel,all NSUBST,100,1000000,10 lswrite,5time,50 ！定義第6載荷步d,298,ux,0.

15、015allsel,all KBC,0allsel,all NSUBST,100,1000000,10 lswrite,6time,60 ！定義第7載荷步d,298,ux,-0.015allsel,all KBC,0allsel,all NSUBST,100,1000000,10 lswrite,7time,70 ！定義第8載荷步d,298,ux,0.015allsel,all KBC,0allsel,all NSUBST,100,1000000,10 lswrite,8time,80 d,298,ux,-0.015allsel,all KBC,0allsel,all NSUBST,100,1

16、000000,10 lswrite,9time,90d,298,ux,0.015allsel,all KBC,0allsel,all NSUBST,100,1000000,10 lswrite,10time,100d,298,ux,-0.015allsel,all KBC,0allsel,all NSUBST,100,1000000,10 lswrite,11time,110d,298,ux,0.03allsel,all KBC,0allsel,all NSUBST,200,1000000,10 lswrite,12time,120d,298,ux,-0.03allsel,all KBC,0a

17、llsel,all NSUBST,200,1000000,10 lswrite,13time,130d,298,ux,0.03allsel,all KBC,0allsel,all NSUBST,200,1000000,10 lswrite,14time,140d,298,ux,-0.03allsel,all KBC,0allsel,all NSUBST,200,1000000,10 lswrite,15time,150d,298,ux,0.03allsel,all KBC,0allsel,all NSUBST,200,1000000,10 lswrite,16time,160d,298,ux,

18、-0.03allsel,all KBC,0allsel,all NSUBST,200,1000000,10 lswrite,17time,170d,298,ux,0.03allsel,all KBC,0allsel,all NSUBST,200,1000000,10 lswrite,18time,180d,298,ux,-0.03allsel,all KBC,0allsel,all NSUBST,200,1000000,10 lswrite,19time,190d,298,ux,0.03allsel,all KBC,0allsel,all NSUBST,200,1000000,10 lswri

19、te,20time,200d,298,ux,-0.03allsel,all KBC,0allsel,all NSUBST,200,1000000,10 lswrite,21time,210d,298,ux,0.045allsel,all KBC,0allsel,all NSUBST,300,1000000,10 lswrite,22lssolve,1,22,1命令流：1、單調(diào)水平荷載作用下的命令流：建模：!mono brick wall!solid65-整體式建模!*Finish/clear/units,si/com,Structural/config,nres,400000/prep7/ti

20、tle,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,1700mp,ex,1,0.24e10!砌體屈服準(zhǔn)則TB,MKIN,1TBTEMP,strainTBdata,1,0.2E-3,1.0E-3,2.

21、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.13E6hntrl=0.21E6tb,concr,1,4,tbdata,0.2,0.9,hntrl,-1 !不考慮壓碎!*-彈性支撐-!彈性支撐暫取鋼筋材料屬性mp,ex,2,2.0e11mp,dens,2,7.600 mp,prxy,2,0.30!real constantsr,1!*-

22、modeling-!geometry modelblock,0,3.12,0,0.24,0,3.12block,-0.24,3.12,0,0.24,3.12,3.36Vsel,s,volu,1,2,1$aslv,s$lsla,slesize,all,0.24Vsel,s,volu,1$vatt,1,1,1, !brick wallVsel,s,volu,2$vatt,2,1,2, !brick wallvsel,all !必須再次選擇所有單元VSWEEP,ALL savenumcmp,allallsel!*/device,vector,1/eshape,1!*-耦合-allsel,all!el

23、astic support nodes setesel,s,type,2$nsle,s,all$nsel,r,loc,z,3.12nsel,u,loc,x,-0.24$cm,slavenode,node!other node setesel,s,type,1$nsle,s,all$cm,masternode,node/input,ucouple,macsavefinish耦合文件ucouple,mac：!*-開(kāi)始-allsel !最好保留這句命令!*將從屬節(jié)點(diǎn)編號(hào)依次存入數(shù)組*cmsel,s,slavenode*get,count1_node,node,0,count*del,slave_no

24、de*dim,slave_node,array,count1_node*get,slave_node(1),node,0,num,min*do,i,2,count1_nodeslave_node(i)=ndnext(slave_node(i-1)*enddo !*將主節(jié)點(diǎn)編號(hào)依次存入數(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_nod

25、emaster_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)=0*enddo!*開(kāi)始選擇程序*allselcmsel,s,masternode*do,i,1,count1_node kk=1 k=1 *dowhile,kk k=nnear(slave_node(i)nsel,s,cp,allkk=nsel(k)+0.001 allselcmsel,s,masternodensel,u,no

26、de,kcm,masternode,node *enddo cp_node(i)=k *enddo!*選擇完畢*!*開(kāi)始耦合*allsel,all/prep7*do,i,1,count1_nodecp,next,UX,slave_node(i),cp_node(i)cp,next,UY,slave_node(i),cp_node(i)cp,next,UZ,slave_node(i),cp_node(i)*enddo!*耦合完畢*求解方法NR法：brick wall NR solution!displacement load!用NR法!*-solution-/soluantype,0!-cons

27、traint-nsel,s,loc,z,0d,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,-0.24!d,all,uzlumpm,onsolcontrol,onpred,onlnsrch,onautots,1 !打開(kāi)自動(dòng)時(shí)間步控制ncnv,2 !如果不收斂時(shí)結(jié)束而不退出compression=0.24E7!*-load steps-!*-1st-acel,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.

28、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 step-nsel,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,1.0E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNS

29、UBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,0.05,2,0.05allsellswrite,2!*-3rd step-nsel,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-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,0.09,2,0.25allsellswrite,3!*-4th ste

30、p-nsel,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,4.0E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,4!*allsel lssolve,1,4,1FINISH也可以用弧長(zhǎng)法，不過(guò)沒(méi)有一定經(jīng)驗(yàn)這比較難控制：!mono brick wall-1 arclength solution!displa

31、cement load!*-solution-/soluantype,0!-constraint-nsel,s,loc,z,0d,all,all!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uy!*lumpm,onsolcontrol,onncnv,2 !如果不收斂時(shí)結(jié)束而不退出compression=0.24E7!*-load steps-!*-first step-acel,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.12esln,s,0,allsfe,all,6,pres,0,0.03*com

32、pressionallsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,1E-5, !使用弧長(zhǎng)法NSUBST,200neqit,15KBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsel,alllswrite,1!*-second step-esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6d,all,ux,2.50E-3!d,all,ux,3.51E-3all

33、sel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,1E-5, !使用弧長(zhǎng)法kbc,0NSUBST,250!arctrm,u,0.0028,313,uxneqit,15nlgeom,1!cnvtol,f,0.005,2,0.01!cnvtol,f,0.05,2,0.5cnvtol,f,0.1,2,0.5allsellswrite,2!*-third step-esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,

34、0,60*0.0593E6!d,all,ux,2.79E-3!d,all,ux,3.18E-3d,all,ux,3.00E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,1E-7, !使用弧長(zhǎng)法kbc,0NSUBST,2000!arctrm,u,0.0028,313,uxneqit,15nlgeom,1!cnvtol,f,0.005,2,0.01!cnvtol,f,0.05,2,0.5cnvtol,f,0.5,2,0.8allsellswrite,3allsel!*-fourth step-esel,s,type,2$nsle,s,all$n

35、sel,r,loc,x,-0.24d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6d,all,ux,3.51E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,1E-7, !使用弧長(zhǎng)法kbc,0NSUBST,300!arctrm,u,0.0028,313,uxneqit,15nlgeom,1!cnvtol,f,0.005,2,0.01!cnvtol,f,0.05,2,0.5cnvtol,f,0.6,2,0.9allsellswrite,4allsel!*

36、lssolve,1,4,1FINISH后處理：/post26esel,s,type,1eplotlines,5000nsol,2,326,u,x,ux!nsol,2,313,u,x,ux*do,i,1,14jj=node(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,fxadd,6,3,6rforce,4,kk,f,x,fxadd,6,4,6! rforce,5,ll,f,x,fx!add,6,5,6*enddoprod,7,2,1000prod,8,6,-1/

37、1000/axlab,x,Displacement(mm)/axlab,y,F(kN)/XRANGE,0,4/YRANGE,0,100xvar,7plvar,8!rforce,5,340,f,x,fx!esol,5,87,221,f,x,fx2、循環(huán)水平荷載作用下的命令流：出求解文件不同外，其他都相同；!cyclic NR solution!displacement load!*-solution-/soluantype,0!-constraint-nsel,s,loc,z,0d,all,all!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all

38、,uy!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,3.36d,all,uy!*lumpm,onsolcontrol,onpred,onlnsrch,onautots,1 !打開(kāi)自動(dòng)時(shí)間步控制ncnv,2 !如果不收斂時(shí)結(jié)束而不退出compression=0.24E7!*-load steps-!*-1st step-acel,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,ERASEOUTR

39、ES,ALL,ALLNSUBST,200,2000,10neqit,15KBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsel,alllswrite,1!*-初裂荷載20%做1次循環(huán)-!*-用NR法-!*-2nd step-esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uz!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,3.36d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6d,all,ux,1.48

40、E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,20,2000,10neqit,15NROPT ,FULL,OFFKBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,2!*-3rd step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-1.48E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,40,2000,20neqit,15NROPT ,FULL,OFFKBC,0nlgeom,1cnvtol

41、,f,0.1,2,0.5allsellswrite,3!*-4th step-nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,1.48E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,40,2000,20neqit,15NROPT ,FULL,OFFKBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,4!*-one cycle-!*-初裂荷載50%做1次循環(huán)-!*-用NR法-!*-5th step-nsel,s,loc,x,-0.24ddele,a

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