地層結(jié)構(gòu)法命令流

finish/clear/prep7! !隧道輪廓k,1001,,,k,1002,0,2.217k,1,0,10.258k,2,-5.832,6.629larc,1,2,1001,6.5 !創(chuàng)建圓弧，第一端點，第二端點，圓心所在一側(cè)任一關(guān)鍵點，半徑k,3,-6.5,3.758larc,2,3,1001,6.5k,4,-6.429,2.629larc,3,4,1001,9k,5,-6.367,2.217larc,4,5,1001,9k, 6,-6.179,1.1345,6k, 7,-5.617,1.1346,7k,8,-5.174,0.874larc,7,8,1001,1.802larc,8,1001,1002,15.75lcomb,2,3 !合并線!LCOMB,NL1,NL2,KEEP -其中：NL1,NL2---擬合并的兩條線號。NL1可為ALL，或元件名。!KEEP---是否保留輸入的線及其公共關(guān)鍵點控制參數(shù)。!KEEP=0則刪除NL1和NL2及其公共關(guān)鍵點，如果已經(jīng)劃分網(wǎng)格則不能刪除，或者依附于其它圖素也不能刪除!KEEP=1則保留線及其公共關(guān)鍵點，但公共關(guān)鍵點不依附于新創(chuàng)建的線。!該命令可以合并獨立線或依附于同面上的線，合并后便于網(wǎng)格劃分。!可合并的線可為直線或曲線，以及直線與曲線，可共線或不共線。!當(dāng)為多條時，應(yīng)為多條首尾相連的線。!無論在何種坐標(biāo)系下執(zhí)行合并，合并后的線不改變合并前的空間位置。lcomb,4,5lcomb,7,8lsymm,x,all!對稱，鏡像!LSYMM,Ncomp,NL1,NL2,NINC,KINC,NOELEM,IMOVE!Ncomp---對稱控制選項，可選X（缺?。?，Y,Z值。!NL1,NL2,NINC---按增量NINC從NL1到NL2定義關(guān)鍵點的范圍（缺省為NL1）,NINC缺省為1。NL1也可為ALL或元件名，此時NP2和NINC將被忽略。!KINC---要創(chuàng)建的關(guān)鍵點編號增量，缺省時由系統(tǒng)自動指定（不會覆蓋）。!NOELEM---是否創(chuàng)建單元和節(jié)點控制參數(shù)。NOELEM=0（缺?。┤绻嬖趩卧凸?jié)點則生成；NOELEM=1不生成單元和節(jié)點。!IMOVE---線是否被移動或重新創(chuàng)建。IMOVE=0（缺?。┰瓉淼木€不動，重新創(chuàng)建新線；當(dāng)IMOVE=1不創(chuàng)建新線，原來的線移動到新位!置，此時編號不變（即ITIME、KINC和NOELEM均無效），且單元和節(jié)點一并移動。!該命令要求當(dāng)前坐標(biāo)系為直角坐標(biāo)系，線可以在任意象限。同KSYMM相同，可通過設(shè)定當(dāng)前坐標(biāo)系為工作平面或局部坐標(biāo)系而改變鏡像位置。nummrg,all!粘接!nummrg是一個使用率極高的命令，使用起來非常的方便。它通常用來實現(xiàn)“粘接”的效果，即，合并坐標(biāo)完全一致的兩個圖素，具體做法就是!刪除編號大的圖素，并用編號小的圖素占據(jù)原來的位置。這樣操作以后，會實現(xiàn)在不同圖素之間形成共用部分，從而在一定程度上取代Vglue命令。!但是，對于復(fù)雜模型，在進(jìn)行圖素合并命令之前，對操作結(jié)果的預(yù)期不明朗之時，還是盡量使用vglue為好，因為這樣能最大限度地保障有限元模型的連續(xù)性。numcmp,all!壓縮編號alls! !建立地層范圍blc4,-40,-40,80,90!命令:BLC4,XCORNER,YCORNER,WIDTH,HEIGHT,DEPTH 通過一角點坐標(biāo)和尺寸創(chuàng)建矩形面!XCORNER,YCORNER---矩形面或塊體第1個角點在工作平面上的X和Y坐標(biāo)。!WIDTH---平行于工作平面X軸方向離XCORNER的距離。!HEIGHT--平行于工作平面Y軸方向離YCORNER的距離!DEPTH---離工作平面的垂直距離，即平行于Z軸。!DEPTH=0（缺?。﹦t生成面。如WIDTH或HEIGHT或DEPTH為負(fù)值，則為反方向距離。asbl,all,all!從面中減去線nummrg,allnumcmp,allalls! !細(xì)分隧道分析模型wprota,,,90!工作平面的旋轉(zhuǎn)!命令：WPROTA,THXYTHYZ,THZX!其中THXY,THYZ,THZX為繞工作平面坐標(biāo)系Z軸、X軸和Y軸的旋轉(zhuǎn)角度。asbw,all !從面中減去與工作面相交的部分wpoffs,,,20!將工作平面沿其自身坐標(biāo)軸移動!命令：WPOFFS,XOFF,YOFF,ZOFF!其中XOFF,YOFF,ZOFF為工作平面坐標(biāo)系內(nèi)沿其X軸、Y軸和Z軸的偏移增量。asbw,all!用工作面切分面wpoffs,,,-40asbw,allwpoffs,,2.139wprota,,90asbw,allwpoffs,,,-4asbw,allwpoffs,,,-20asbw,allwpoffs,,,46asbw,allwpoffs,,,-20.995wprota,,90wpoffs,,,-13.821asel,s,,,24,29,5asbw,allwpoffs,,,-0.562asel,s,,,10,17,7asel,a,,,13asbw,allwpoffs,,,-11.234asel,s,,,23,30,7asel,a,,,14asbw,allwpoffs,,,-0.562asel,s,,,35,37,2asbw,allwprota,,90asel,s,,,18,38,20asel,a,,,27,28,1asbw,alllcomb,1,8lcomb,3,14lcomb,54,37lcomb,57,39wpcsys,-1!使工作平面的位置與當(dāng)前坐標(biāo)系的X-Y平面或R-0平面一致，但不改變視圖方向。!WPCSYS，WN，KCN!WN是窗口號，該窗口中的視圖方向修改為與工作平面垂直。默認(rèn)值是1。如果WN為負(fù)值，則視圖方向不變。!KCN是坐標(biāo)系號，默認(rèn)為當(dāng)前激活的坐標(biāo)系。如果KCN是直角坐標(biāo)系，工作平面與其X-Y平面一致；如果是圓柱或球坐標(biāo)系，則與R-0平面一致。allsnummrg,allnumcmp,all!設(shè)置實體單元類型ET,1,PLANE42KEYOPT,1,1,0 !參見ElementsReference部分的ElementLibrary中的plane42KEYOPT,1,2,0KEYOPT,1,3,2 !設(shè)置為平面應(yīng)變模式KEYOPT,1,5,0KEYOPT,1,6,0mp,ex,1,1.5e9 !!地層!!加固地層 mp--定義線性材料屬性mp,prxy,1,0.4mp,dens,1,2000tb,dp,1 !激活非線性材料的數(shù)據(jù)表tbdata,1,0.15e6,25 !定義tb數(shù)據(jù)表中的數(shù)據(jù)et,2,beam3*set,th1,0.25r,2,th1,th1*th1*th1/12,th1mp,ex,2,2.555e10 !!初期支護(hù)mp,prxy,2,0.2mp,dens,2,2500mp,ex,3,1.5e9 !!上臺階內(nèi)部土體mp,prxy,3,0.4mp,dens,3,2000tb,dp,3tbdata,1,0.15e6,25mp,ex,4,1.5e9 !!下臺階內(nèi)部土體mp,prxy,4,0.4mp,dens,4,2000tb,dp,4tbdata,1,0.15e6,25mp,ex,5,1.5e9 !!仰拱內(nèi)部土體mp,prxy,5,0.4mp,dens,5,2000tb,dp,5tbdata,1,0.15e6,25! !!!初次劃分單元網(wǎng)格lsel,s,,,1,3,2lsel,a,,,11,13,2lsel,a,,,4lesize,all,,,12 !確定劃分指定線的網(wǎng)格數(shù)和網(wǎng)格尺寸alls!LESIZE,NL1,SIZE,ANGSIZ,NDIVSPACE,KFORC,LAYER1,LAYER2,KYNDIV!NL1---Numberofthelinetobemodified.IfALL,modifyallselectedlines[LSEL].IfNL1=P,graphical!pickingisenabledandallremainingcommandfieldsareignored(validonlyintheGUI).AcomponentnamemayalsobesubstitutedforNL1.!SIZE---IfNDIVisblank,SIZEisthedivision(elementedge)length.Thenumberofdivisionsisautomaticallycalculatedfromthelinelength(roundedupwardto! nextinteger).IfSIZEiszero(orblank),useANGSIZorNDIV.!ANGSIZ---Thedivisionarc(indegrees)spannedbytheelementedge(exceptforstraightlines,whichalwaysresultinonedivision).Thenumberofdivisionsis! automaticallycalculatedfromthelinelength(roundedupwardtonextinteger).!NDIV---Ifpositive,NDIVisthenumberofelementdivisionsperline.If-1(andKFORC=1),NDIVisassumedtobezeroelementdivisionsperline.TARGE169with!arigidspecificationignoresNDIVandwillalwaysmeshwithoneelementdivision.!SPACE---Spacingratio.Ifpositive,nominalratiooflastdivisionsizetofirstdivisionsize(if>1.0,sizesincrease,if<1.0,sizesdecrease).Ifnegative,!|SPACE|isnominalratioofcenterdivision(s)sizetoenddivisionssize.Ratiodefaultsto1.0(uniformspacing).Forlayer-meshing,avalueof1.0! normallyisused.IfSPACE=FREE,ratioisdeterminedbyotherconsiderations!KFORC---KFORC0-3areusedonlywithNL1=ALL.Specifieswhichselectedlinesaretobemodified.! 0 — Modifyonlyselectedlineshavingundefined(zero)divisions.! 1 — Modifyallselectedlines.! 2 — Modifyonlyselectedlineshavingfewerdivisions(includingzero)thanspecifiedwiththiscommand.! 3 — Modifyonlyselectedlineshavingmoredivisionsthanspecifiedwiththiscommand.! 4 — ModifyonlynonzerosettingsforSIZE,ANGSIZ,NDIV,SPACE,LAYER1,andLAYER2.IfKFORC=4,blankor0settingsremainunchanged.!LAYER1---Layer-meshingcontrolparameter.Distancewhichdefinesthethicknessoftheinnermeshlayer,LAYER1.Elementsinthislayerareuniformly-sizedwith! edgelengthsequaltothespecifiedelementsizefortheline(eitherthroughSIZEorline-length/NDIV).ApositivevalueforLAYER1isinterpretedas! anabsolutelength,whileanegativevalueininterpretedasamultiplieronthespecifiedelementsizefortheline.Asageneralrule,theresulting! thicknessoftheinnermeshlayershouldbegreaterthanorequaltothespecifiedelementsizefortheline.IfLAYER1=OFF,layer-meshingcontrol! settingsareclearedfortheselectedlines.Thedefaultvalueis0.0!LAYER2---Layer-meshingcontrolparameter.Distancewhichdefinesthethicknessoftheoutermeshlayer,LAYER2.Elementsinthislayertransitioninsizefrom! thoseinLAYER1totheglobalelementsize.ApositivevalueofLAYER2isinterpretedasanabsolutelength,whileanegativevalueisinterpretedasa! meshtransitionfactor.AvalueofLAYER2=-2wouldindicatethatelementsshouldapproximatelydoubleinsizeasthemeshprogressesnormaltoLAYER1.! Thedefaultvalueis0.0.!KYNDIV---! 0,No,andOffmeansthatSmartSizingcannotoverridespecifieddivisionsandspacingratios.Mappedmeshfailsifdivisionsdonotmatch.Thisdefines! thespecificationas"hard".! 1,Yes,andOnmeansthatSmartSizingcanoverridespecifieddivisionsandspacingratiosforcurvatureorproximity.Mappedmeshingcanoverride! divisionstoobtainrequiredmatchingdivisions.Thisdefinesthespecificationas"soft".lsel,s,,,35,37,2lsel,a,,,48,51,1lsel,a,,,52lesize,all,,,8allslsel,s,,,36,70,34lsel,a,,,7,9,2lsel,a,,,71,76,5lsel,a,,,74,80,6lesize,all,,,11allslsel,s,,,33,65,32lsel,a,,,6,8,2lsel,a,,,38,69,31lsel,a,,,72,83,11lesize,all,,,1allslsel,s,,,39,41,1lsel,a,,,25,32,7lsel,a,,,87,91,4lsel,a,,,86,90,4lesize,all,,,2allslsel,s,,,55,56,1lsel,a,,,53,54,1lsel,a,,,58,60,2lsel,a,,,61lsel,a,,,42,43,1lsel,a,,,77,88,11lsel,a,,,27,81,54lsel,a,,,15,64,49lsel,a,,,21,24,3lesize,all,,,6allslsel,s,,,12,19,7lsel,a,,,17,30,13lsel,a,,,2,5,3lsel,a,,,34lsel,a,,,16,31,15lsel,a,,,57,59,2lsel,a,,,44,47,3lsel,a,,,26,29,3lsel,a,,,85,92,7lsel,a,,,89,93,4lsel,a,,,66,68,1lsel,a,,,22,75,53lsel,a,,,78,84,6lsel,a,,,23,62,39lsel,a,,,10,18,8lsel,a,,,45,46,1lsel,a,,,28,73,45lsel,a,,,79,82,3lsel,a,,,14,20,6lsel,a,,,63lesize,all,,,6alls!定義單元的材料屬性及單元范圍，并最終劃分單元網(wǎng)格type,1!選取單元類型mat,1 !單元材料屬性amap,24,31,36,1,26 !簡化面映射網(wǎng)格劃分，amap,area,kp1,kp2,kp3,kp4——中，area---擬劃分的面號，kp1、kp2、kp3、kp4---指定的角點，3或4個都可以，順序可任意。amap,25,29,36,1,27amap,19,31,26,23,21amap,21,29,27,24,19amap,38,21,23,2,54amap,39,19,24,4,52allstype,1mat,1asel,s,,,3,9,1asel,a,,,11,12,1asel,a,,,14,17,1asel,a,,,21,22,1asel,a,,,24,27,1asel,a,,,29,40,1asel,a,,,19amesh,all !AMESH,NA1,NA2,NINC---在幾個面上生成面單元，na1,na2,ninc--面編號范圍和編號增量，na1可取all、p、元件名或組件名。allstype,1mat,3amap,1,26,1,25amap,2,27,1,25allstype,1mat,4amap,18,26,25,22,23amap,20,27,25,22,24allstype,1mat,5amap,23,23,44,3,2amap,28,44,22,3,6amap,13,47,22,6,5amap,10,24,47,5,4allstype,2mat,2real,2lsel,s,,,1,3,2lsel,a,,,35,37,2lsel,a,,,39,41,2lsel,a,,,6,9,1lmesh,all !LMESH,NA1,NA2,NINC---在幾個線上生成線單元，na1,na2,ninc--線編號范圍和編號增量，na1可取all、p、元件名或組件名alls! !添加約束antype,0 !指定分析類型為靜態(tài)分析 ANTYPE,Antype,Status,LDSTEP,SUBSTEP,Action---指定分析類型和重新啟動的狀態(tài)nsel,s,loc,x,-40 !選擇所有x坐標(biāo)為-40的結(jié)點創(chuàng)建子集 NSEL,Type,Item,Comp,VMIN,VMAX,VINC,KABSd,all,ux,0 !對所有所選結(jié)點的x放下添加約束，位移值為0nsel,s,loc,x,40d,all,ux,0nsel,s,loc,y,-40d,all,uy,0allsacel,,10 !線性加速度，! !第一步，初始狀態(tài)分析/solu!!1antype,static !!聲明分析類型：靜力分析esel,s,type,,2 !選定第2類單元創(chuàng)建子集ekill,all !殺死所有所選單元allsTIME,1!時間，“1”為荷載步結(jié)束的時間值!deltim,0.2,0.02,0.4 !設(shè)定時間步長，0.2是當(dāng)前荷載步時間步長，0.02,0.4分別為采用自動時間步時代最小時間步長和最大時間步長，此命令可以和NSUBST命令互相替代NSUBST,100,10000000 !設(shè)定子步數(shù)和時間步長，100為當(dāng)前荷載步的子步數(shù)，100000000為最大子步數(shù)，后面還有兩個選項分別是最小子步數(shù)和時間步長繼承控制參數(shù)NEQIT,100 !定義平衡迭代的最大次數(shù)，100為每個子步中平衡迭代的最大次數(shù)autots,on!使用自動時間步若為off則不適應(yīng)自動時間步pred,on!打開時間步長預(yù)測器 pred,sskey,--,lskey其中sskey為子步預(yù)測控制參數(shù)，off為關(guān)閉，on為開，默認(rèn)為on，lskey為荷載步預(yù)測控制參數(shù)，默認(rèn)為off。lnsrch,on!打開線性搜索cutcontrol,plslimit,0.5!回退控制，plslimit為回退準(zhǔn)則控制，0.5為回退準(zhǔn)則的數(shù)值。nlgeom,on!打開大位移效果nropt,full,,on!設(shè)定牛頓-拉普森選項cnvtol,f,,0.01,2!設(shè)定非線性收斂值(收斂準(zhǔn)則)，此處出現(xiàn)警告。OUTRES,basic,1 !數(shù)據(jù)寫入KBC,0 !線性增加荷載savesolve!rescontrol,file_summany!控制文件編寫多幀重新啟動! ! —!提取應(yīng)力/post1set,1 !定義數(shù)據(jù)集asel,s,,,1,2,1!提取上臺階節(jié)點力esla,s,1!選擇與所選面相關(guān)的元素nsle,s !選擇已選單元的相關(guān)節(jié)點nsel,r,ext !建立一個新的節(jié)點子集，選擇的節(jié)點為體單元表面上的節(jié)點，面單元外輪廓上的節(jié)點nplot !顯示節(jié)點*get,n,node,,count !提取節(jié)點數(shù)量， *status,n命令可獲得n的數(shù)值列表*dim,n_num2,array,1,n!定義大小為n=48的一維數(shù)組*do,i,1,n !從小到大排序*get,n_num2(1,i),node,,num,min!此命令為取當(dāng)前節(jié)點子集中最小的節(jié)點號，在此循環(huán)中連取48次nsel,u,,,n_num2(1,i)!從當(dāng)前子集中去掉一組子集，此命令是為了去掉上一行命令*get命令中第i此獲取的最小節(jié)點號，以保證*get命令在第i+1次獲取節(jié)點號時，所要進(jìn)行選擇的!子集中沒有第i次已經(jīng)被選擇的最小節(jié)點號，從而選取到第i+1此選擇中最小節(jié)點號，但卻比第i次選擇的節(jié)點號大的節(jié)點號，由此使得數(shù)組n_num2中節(jié)點從小到大排列*enddo*status,n_num2*dim,n_2fx,array,1,n !定義一維數(shù)組大小為n=48存儲x方向節(jié)點反力*dim,n_2fy,array,1,n*do,i,1,nnsel,s,,,n_num2(1,i) !從小到大選擇節(jié)點fsum!單元節(jié)點力或矩的綜合*get,n_2fx(1,i),fsum,,item,fx!獲取所有節(jié)點x方向的節(jié)點反力*enddo*status,n_2fx!顯示數(shù)組n_2fx*do,i,1,nnsel,s,,,n_num2(1,i)!從小到大選擇節(jié)點fsum !單元節(jié)點力的綜合*get,n_2fy(1,i),fsum,,item,fy!獲取所有節(jié)點y方向的節(jié)點反力*enddo*status,n_2fy!顯示數(shù)組n_2fyparsav,all,001!將所有參數(shù)寫入001.ext文件中! !第二步，挖去上臺階土體，釋放40%荷載/solu!!2antype,,rest,1!在第一個子步重新啟動進(jìn)行靜態(tài)分析esel,s,mat,,3 !選擇材料屬性3的單元，建立新子集ekill,all!殺死所選單元allsesel,s,type,,2 !選擇單元類型2，建立新子集ekill,allalls!應(yīng)力釋放parres,,001!讀取文件001.ext的數(shù)據(jù)allsasel,s,,,1,2,1!選擇面1和面2esla,s,1!選擇與所選面相關(guān)的元素nsle,s!選擇已選單元的相關(guān)節(jié)點nsel,r,ext!建立一個新的節(jié)點子集*do,i,1,nf,n_num2(1,i),fx,0.6*n_2fx(1,i) !？？？？？荷載釋放率怎么確定？？？f,n_num2(1,i),fy,0.6*n_2fy(1,i)*enddoallsTIME,2!deltim,0.2,0.02,0.4 !設(shè)定時間步長NSUBST,100,10000000 !設(shè)定每一步子步數(shù)NEQIT,100autots,on!使用自動時間步pred,on!打開時間步長預(yù)測器lnsrch,on!打開線性搜索cutcontrol,plslimit,0.5nlgeom,on!打開大位移效果nropt,full,,on!設(shè)定牛頓-拉普森選項cnvtol,f,,0.01,2!設(shè)定非線性收斂值OUTRES,basic,1 !數(shù)據(jù)寫入KBC,0 !線行增加荷載，即斜坡荷載savesolve!rescontrol,file_summany! ! !第三步，繼第二步后對上臺階支護(hù)/solu!!3antype,,rest,2 !在第二個子步重新啟動進(jìn)行靜態(tài)分析esel,s,type,,2 !選擇單元類型2建立子集ekill,allallsesel,s,mat,,3 !選擇材料屬性3的單元，建立子集ekill,allallslsel,s,,,1,3,2esll,s!選擇與所選線相關(guān)的單元ealive,all !激活所選單元allsasel,s,,,1,2,1!應(yīng)力釋放esla,s,1 !選擇與所選面相關(guān)的元素nsle,s !選擇已選單元的相關(guān)節(jié)點nsel,r,ext!!建立一個新的節(jié)點子集，選擇的節(jié)點為體單元表面上的節(jié)點，面單元外輪廓上的節(jié)點fdele,all,all!去掉節(jié)點力alls! -TIME,3!deltim,0.2,0.02,0.4 !設(shè)定時間步長NSUBST,100,10000000 !設(shè)定每一步子步數(shù)NEQIT,100autots,on!使用自動時間步pred,on!打開時間步長預(yù)測器lnsrch,on!打開線性搜索cutcontrol,plslimit,0.5nlgeom,on!打開大位移效果nropt,full,,on!設(shè)定牛頓-拉普森選項cnvtol,f,,0.01,2!設(shè)定非線性收斂值OUTRES,basic,1!數(shù)據(jù)寫入KBC,0 !線行增加荷載savesolve!rescontrol,file_summany!提取應(yīng)力/post1set,1asel,s,,,18,20,2!應(yīng)力釋放esla,s,1 !選擇與所選面相關(guān)的元素nsle,s !選擇已選單元的相關(guān)節(jié)點nsel,r,ext !建立一個新的節(jié)點子集*get,n,node,,count !提取節(jié)點數(shù)量*dim,n_num3,array,1,n !定義一維數(shù)組大小為n*do,i,1,n !從小到大排序*get,n_num3(1,i),node,,num,min!選取最小的節(jié)點nsel,u,,,n_num3(1,i) !從當(dāng)前子集中去掉所選節(jié)點*enddo*dim,n_3fx,array,1,n !定義一維數(shù)組大小為n存儲節(jié)點反力*dim,n_3fy,array,1,n*do,i,1,nnsel,s,,,n_num3(1,i) !從小到大選擇節(jié)點fsum*get,n_3fx(1,i),fsum,,item,fx!獲取x方向的節(jié)點反力*enddo*do,i,1,nnsel,s,,,n_num3(1,i) !從小到大選擇節(jié)點fsum*get,n_3fy(1,i),fsum,,item,fy!獲取y方向的節(jié)點反力*enddoparsav,all,002 !將所有參數(shù)寫入002.ext文件中! !第四步，繼第三步后再開挖下臺階土體，并釋放30%的荷載/solu!!4antype,,rest,3 !在第三個子步重新啟動進(jìn)行靜態(tài)分析esel,s,type,,2!選擇單元類型2ekill,allallsesel,s,mat,,3!選擇材料3的單元ekill,allallsesel,s,mat,,4 !選擇材料4的單元ekill,allallslsel,s,,,1,3,2esll,s !選擇與所選線相關(guān)的單元ealive,all !激活所選單元allsparres,,002 !讀取文件002.ext的數(shù)據(jù)allsasel,s,,,18,20,2!應(yīng)力釋放esla,s,1nsle,snsel,r,ext*do,i,1,nf,n_num3(1,i),fx,0.7*n_3fx(1,i)f,n_num3(1,i),fy,0.7*n_3fy(1,i)*enddo!去掉下臺階上邊緣的節(jié)點力*do,i,1281,1292fdele,i,all !去掉節(jié)點力*enddo*do,i,1395,1405fdele,i,all!去掉節(jié)點力*enddoalls! TIME,4!deltim,0.2,0.02,0.4 !設(shè)定時間步長NSUBST,100,10000000 !設(shè)定每一步子步數(shù)NEQIT,100autots,on!使用自動時間步pred,on!打開時間步長預(yù)測器lnsrch,on!打開線性搜索cutcontrol,plslimit,0.5nlgeom,on!打開大位移效果nropt,full,,on!設(shè)定牛頓-拉普森選項cnvtol,f,,0.01,2!設(shè)定非線性收斂值OUTRES,basic,1 !數(shù)據(jù)寫入KBC,0 !線行增加荷載savesolve!rescontrol,file_summany! !第五步，下臺階支護(hù)/solu!!5antype,,rest,4 !在第四個子步重新啟動進(jìn)行靜態(tài)分析esel,s,type,,2!選擇單元類型2ekill,allallsesel,s,mat,,3 !選擇材料3的單元ekill,allallsesel,s,mat,,4 !選擇材料4的單元ekill,allallslsel,s,,,1,3,2!選擇與所選線相關(guān)的單元!選擇與所選線相關(guān)的單元ealive,allallslsel,s,,,35,37,2esll,s !選擇與所選線相關(guān)的單元ealive,allallsasel,s,,,18,20,2!應(yīng)力釋放esla,s,1 !選擇與所選面相關(guān)的元素nsle,s !選擇已選單元的相關(guān)節(jié)點nsel,r,ext!建立一個新的節(jié)點子集fdele,all,all !去掉所有節(jié)點力alls! ! TIME,5!deltim,0.2,0.02,0.4 !設(shè)定時間步長NSUBST,100,10000000 !設(shè)定每一步子步數(shù)NEQIT,100autots,on!使用自動時間步pred,on!打開時間步長預(yù)測器lnsrch,on!打開線性搜索cutcontrol,plslimit,0.5nlgeom,on!打開大位移效果nropt,full,,on!設(shè)定牛頓-拉普森選項cnvtol,f,,0.01,2!設(shè)定非線性收斂值OUTRES,basic,1!數(shù)據(jù)寫入KBC,0 !線行增加荷載savesolve!rescontrol,file_summany!提取應(yīng)力/post1set,1asel,s,,,23,28,5!應(yīng)力釋放asel,a,,,10,13,3esla,s,1nsle,snsel,r,ext*get,n,node,,count !提取節(jié)點數(shù)量*dim,n_num4,array,1,n !定義一維數(shù)組大小為n*do,i,1,n !從小到大排序*get,n_num4(1,i),node,,num,min!從小到大提取節(jié)點號nsel,u,,,n_num4(1,i)*enddo*dim,n_4fx,array,1,n !定義一維數(shù)組大小為n存儲節(jié)點力*dim,n_4fy,array,1,n*do,i,1,nnsel,s,,,n_num4(1,i) !從小到大選擇節(jié)點fsum*get,n_4fx(1,i),fsum,,item,fx*enddo*do,i,1,nnsel,s,,,n_num4(1,i)fsum*get,n_4fy(1,i),fsum,,item,fy*enddoparsav,all,003! !第六步，開挖仰拱土體，并釋放20%對荷載/solu!!6antype,,rest,5esel,s,type,,2ekill,allallsesel,s,mat,,3ek