1、AGLUE, all!粘結所有面/prep7L1=30!設置變量L2=30h=-25K, 1,0, 0, 0,K, 2, L1, 0, 0,K, 3, L1, L2, 0,K, 4, 0, L2, 0,KWPAVE, 1!將工作平面原點定義在1號點RECTNG, 0, L1,0, L2,wpro, , -90,!將工作平面繞X軸Z到Y方向90度RECTNG , 0, L1,0, -h,KWPAVE, 4!將工作平面原點定義在4號點RECTNG , 0, L1,0, -h,wpro, ,90!將工作平面繞y軸x到z方向90度RECTNG , 0, L2, 0, -h,KWPAVE, 3!將工作平

2、面原點定義在3號點RECTNG , 0, L2, 0, -h,(/ =20kPaET, 1, SHELL43元）!ET，ITYPE，Ename，KOPT1，KOPT6，INOPR（定義單!K0PT1K0PT6為元素特性編碼ET, 2, COMBIN14!C0MBIN14彈簧-阻尼器Spring-DamperMPTEMP,MPTEMP, 1, 0!刪除系統中已存在的溫度表!定義一個溫度表MPDATA, EX, 1, , 2.4E10MPDATA, PRXY, 1, , 0.15!指定與溫度相應的材料性能數據彈性模量!主泊松比ESIZE, 1, 0AMESH, ALL!指定單元邊長!劃分面生成面單

3、元NSEL, S, LOC, Z, 0ESLN, SNSEL, S, LOC, Z, -1ESLN, UCM, STRUT, ELEM allsCMSEL, U, STRUT CM,WALL, ELEM!選擇一組節點子集創建新集!選擇已選節點上的單元!選擇z坐標值為-1的-!從已選集中刪除此時剩下只支撐板!將選擇集命名STRUT生成元件!all sel全選!去除STRUT元件!將選擇集命名wall生成元件NSEL, S, LOC, X, 0.1, L1-0.1NPLOTNSEL, R, LOC, Y, 0ESLN, SNSEL, S, LOC, Y, 1ESLN, UENSYM, , , ,

4、ALL!選擇一組節點子集創建新集!顯示節點!從當前集選擇一組節點子集!從已選集中選擇!從當前集選擇一組節點子集!從已選集中刪除!反轉殼單元法線方向NSEL, S, LOC, Y, 0.1, L2-0.1NPLOTNSEL, R, LOC, X, 0ESLN, SNSEL, S, LOC, X, 1ESLN, UENSYM , , , , ALL!選擇一組節點子集創建新集!顯示節點!從當前集選擇一組節點子集!從已選集中選擇!從當前集選擇一組節點子集!從已選集中刪除!反轉殼單元法線方向ALLSNUMCMP, ALL!所有實體進行重新編號!直接生成節點*D0, i, 1, L1-1!從1到29進行循

5、環CSYS, 0!激活默認笛卡爾坐標系N, 100000+2*i-1, i, 0,N, 100000+2*i, i, 2.5,*enddo-1!節點編號 后面為坐標-1!節點編號 后面為坐標*DO, i, 1, L1-1!從1到29進行循環!ModelingCreatElementsElem AttributesTYPE, 2MAT , 1REAL, 0ESYS, 0!設置單元類型屬性指示器!MP命令中的MAT即材料性能!材料實常數!材料坐標系統屬性指示器EN, 100000+i, 100000+2*i-1, 100000+2*i!根據給定的單元號和節點號生成單元*enddoNUMCMP ,

6、ALL!所有實體進行重新編號*DO, i, 1, L1-1CSYS, 0!激活默認笛卡爾坐標系N, 100000+2*i-1, i, L2, -1N, 100000+2*i, i, L2-2.5, -1*enddo*DO, i, 1, L1-1TYPE, 2MAT , 1REAL, 0ESYS, 0EN, 100000+i, 100000+2*i-1, 100000+2*i*enddoNUMCMP , ALL!所有實體進行重新編號*DO, i, 1, L2-1CSYS, 0N, 100000+2*i-1, 0, i, -1N, 100000+2*i, 2.5, i, -1 *enddo*DO,

7、 i, 1, L2-1TYPE, 2MAT , 1REAL, 0ESYS, 0EN, 100000+i, 100000+2*i-1, 100000+2*i *enddoNUMCMP , ALL!所有實體進行重新編號*DO, i, 1, L2-1CSYS, 0N,100000+2*i-1, L1, i, -1N,100000+2*i, L1-2.5, i, -1*enddo*DO, i, 1, L2-1TYPE, 2 $ MAT , 1 $ REAL , 0ESYS, 0EN, 100000+i, 100000+2*i-1, 100000+2*i*enddoNUMCMP , ALL!所有實體進行

8、重新編號ESEL, S, TYPE, , 2!選擇單元類型號為2的單元EGEN, 25, 100000, ALL, , , , , , , , , , -1,ESEL, S, TYPE, , 2!選擇單元類型號為2的單元CM, SPRING, ELEM!生成一個叫SPRING由單元組成的元件ALLSALLSEL, ALLNUMMRG, NODE, 0.01, 0.01, LOWNUMCMP, ALL!選擇所有實體!節點合并 距離小于0.01則同保留編!所有實體進行重新編號!從當前集選擇一組節點子集 深度-1以下!從已選集中選擇!選擇單元類型號為2的單元!單位面積內受的力隨深度增加而增加!對已存

9、在單元進行修改R, 101, 0.6, 0.6, 0.6, 0.6, ,R, 102, 0.1, 0.1, 0.1, 0.1, ,CMSEL, S, WALL底的點*DO, i, 1, 25NSEL, S, LOC, Z, -1*iESLN, SNSEL, R, TYPE, , 2R, i, m1*i*b*h, , ,EMODIF, ALL, REAL, i,*enddo!識別號+實常數ALLSANTYPE, STATIC, NEWNROPT, FULL*AFUN, DEGQ=2.0E4GAMA1=0.9E4FAI1=20C1=10e3m1=1500e3b=1h=1!分析模式 靜力!指定計算模

10、式!指定角度單位為 度!Q為超載!浮重度!內摩擦角!粘聚力!比例系數!單元寬度!墻體單元高度EMODIF, ALL, REAL, 101,CMSEL, S, STRUTEMODIF, ALL, REAL, 102,!將實常數101組賦給墻NSEL, S, LOC, Z-25D, ALL, , , , , , UZ, , , , ,CMSEL, S, SPRINGNSLE, SCMSEL, S, WALLNSLE, UD, ALL, , , , , , ALL, , , , ,ALLSSAVE/SOLU!賦值給支撐!約束墻底豎向位移!選擇土彈簧單元!以下命令從已選彈簧集合中選!在已選集中選墻單元

11、!從集合中刪除剛選擇的單元，即與墻有關的單元!約束土彈簧單元端點的所有位移!無支撐開挖1m坑內水位-1.0坑外0NSEL, S, LOC, Z, 0, -0.9ESLN, SCMSEL, R, wallNSLE, S!0到0.9的位置即開挖面以上!以下命令在當前集里選擇!在當前集選wll單元!以下命令在當前集里選擇*GET, ZMIN, NODE, , MNLOC, Z, , , ,*GET, ZMAX, NODE, , MXLOC, Z, , , ,!墻最淺節點處的位置!墻最深節點處的位置LOCZ1=abs (ZMAX)LOCZ2=abs (ZMIN)KA= (TAN(45.0 - FAI1

12、/2)*2!取絕對值!主動土壓力系數!采用水土分算PA1=(Q + GAMA1*LOCZ1)*KA-2.0*C1*SQRT(KA) !最深主動土壓力公式ENDIFPA2=10*1E3*LOCZ1PA=PA1+PA2SA1=(Q+GAMA1*LOCZ2)*KA-2*C1*SQRT(KA)SA2=10*1E3*LOCZ2WA=SA2SA=SA1+SA2pressure=SASLZER=ZMAXSLOPE=(SA-PA)/(LOCZ2-LOCZ1)!增長率SFGRAD, PRES, 0, z, SLZER, SLOPE!沿z方向從-1m開始的面力減少SFE, all, 1, PRES, , -PA,

13、 , ,!alls!開挖面以下加載NSEL, S, LOC, Z, -1.1, -24.9ESLN, SGMSEL, R, wallNSLE, S*GET, ZMIN, NODE, , MNLOC, Z, , , ,*GET, ZMAX, NODE, , MXLOC, Z, , , , SLZER=ZMAXSLOPE=0SFGRAD, PRES, O, z, SLZER, SLOPE SFE, all, 1, PRES, , -pressure, , ,!alls!將集合擴大到全集cmsel, S, strut!選擇內支撐 全部殺死EKILL, ALLallsNSEL, S, LOC, Z,

14、0.1, -1.1ESLN, SESEL, R, TYPE, , 2!選擇單元類型號為2的單元EKILL, ALL*do, i, 2, 25z=i-1!開挖了1m要在被動區減去相應的土壓力* IF , PA1, LT, 0, THENPA1=0!去除小于零值的可能!水壓力計算!總應力K=m1*z*b*hR, i, K, , , *enddoTIME, 1 alls solveNSEL, S, LOC, Z, 0, -4.9ESLN, SCMSEL, R, wallNSLE, S*GET, ZMIN, NODE, , MNLOC, Z, , , ,*GET, ZMAX, NODE, , MXLO

15、C, Z, , , ,LOCZ1=abs (ZMAX)LOCZ2=abs (ZMIN)KA= (TAN(45.0 - FAI1/2)*2PA1=(Q + GAMA1*LOCZ1)*KA-2.0*C1*SQRT(KA)* IF , PA1, LT, 0, THENPA1=0*ENDIFPA2=10*1E3*LOCZ1PA=PA1+PA2SA1=(Q+GAMA1*LOCZ2)*KA-2*C1*SQRT(KA)SA2=10*1E3*LOCZ2WA=SA2SA=SA1+SA2 pressure=SASLZER=ZMAXSLOPE=(SA-PA)/(LOCZ2-LOCZ1)SFGRAD, PRES, 0, z, SLZER, SLOPESFE, all, 1, PRES, , -PA, , , allsNSEL, S, LOC, Z, -5.1, -24.9ESLN, SGMSEL, R, wall!載荷步1NSLE, S*GET, ZMIN, NODE, , MNLOC, Z, , , ,*GET, ZMAX, NODE, , MXLOC, Z, , , , SLZER=ZM