ANSYS基坑弹性地基梁全程序即详解.doc

/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 !将工作平面原点定义在3号点RECTNG, 0, L2, 0, -h,AGLUE, all !粘结所有面ET, 1, SHELL43 !ET，ITYPE，Ename，KOPT1，KOPT6，INOPR（定义单元） !KOPT1KOPT6为元素特性编码 !shell43 4 节点塑性大应变单元ET, 2, COMBIN14 !COMBIN14弹簧-阻尼器Spring-DamperMPTEMP, !删除系统中已存在的温度表MPTEMP, 1, 0 !定义一个温度表MPDATA, EX, 1, , 2.4E10 !指定与温度相应的材料性能数据 弹性模量MPDATA, PRXY, 1, , 0.15 !主泊松比 ESIZE, 1, 0 !指定单元边长AMESH, ALL !划分面生成面单元NSEL, S, LOC, Z, 0 !选择一组节点子集创建新集ESLN, S !选择已选节点上的单元 NSEL, S, LOC, Z, -1 !选择z坐标值为-1的- ESLN, U !从已选集中删除此时剩下只支撑板 CM, STRUT, ELEM !将选择集命名STRUT生成元件alls !all sel 全选 CMSEL, U, STRUT !去除STRUT元件CM, WALL, ELEM !将选择集命名wall生成元件NSEL, S, LOC, X, 0.1, L1-0.1 !选择一组节点子集创建新集NPLOT !显示节点NSEL, R, LOC, Y, 0 !从当前集选择一组节点子集ESLN, S !从已选集中选择NSEL, S, LOC, Y, 1 !从当前集选择一组节点子集ESLN, U !从已选集中删除ENSYM, , , , ALL !反转壳单元法线方向NSEL, S, LOC, Y, 0.1, L2-0.1 !选择一组节点子集创建新集 NPLOT !显示节点NSEL, R, LOC, X, 0 !从当前集选择一组节点子集ESLN, S !从已选集中选择NSEL, S, LOC, X, 1 !从当前集选择一组节点子集ESLN, U !从已选集中删除ENSYM , , , , ALL !反转壳单元法线方向ALLSNUMCMP, ALL !所有实体进行重新编号 !直接生成节点*DO, i, 1, L1-1 ! 从1到29进行循环CSYS, 0 !激活默认笛卡尔坐标系N, 100000+2*i-1, i, 0, -1 ! 节点编号 后面为坐标N, 100000+2*i, i, 2.5, -1 ! 节点编号 后面为坐标*enddo*DO, i, 1, L1-1 ! 从1到29进行循环 !ModelingCreatElementsElem AttributesTYPE, 2 !设置单元类型属性指示器 MAT , 1 !MP命令中的MAT即材料性能REAL, 0 !材料实常数ESYS, 0 !材料坐标系统属性指示器EN, 100000+i, 100000+2*i-1, 100000+2*i !根据给定的单元号和节点号生成单元*enddoNUMCMP , 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, 2 MAT , 1 REAL, 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, i, 1, L2-1TYPE, 2 MAT , 1 REAL, 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 !所有实体进行重新编号ESEL, S, TYPE, , 2 !选择单元类型号为2的单元EGEN, 25, 100000, ALL, , , , , , , , , , -1,ESEL, S, TYPE, , 2 !选择单元类型号为2的单元CM, SPRING, ELEM !生成一个叫SPRING由单元组成的元件ALLSALLSEL, ALL !选择所有实体NUMMRG, NODE, 0.01, 0.01, LOW !节点合并 距离小于0.01则同保留编号底的点NUMCMP, ALL !所有实体进行重新编号*DO, i, 1, 25NSEL, S, LOC, Z, -1*i !从当前集选择一组节点子集 深度-1以下ESLN, S !从已选集中选择NSEL, R, TYPE, , 2 !选择单元类型号为2的单元R, i, m1*i*b*h, , , !单位面积内受的力随深度增加而增加 EMODIF, ALL, REAL, i, !对已存在单元进行修改*enddo R, 101, 0.6, 0.6, 0.6, 0.6, , !识别号+实常数R, 102, 0.1, 0.1, 0.1, 0.1, ,CMSEL, S, WALLEMODIF, ALL, REAL, 101, !将实常数101组赋给墙CMSEL, S, STRUTEMODIF, ALL, REAL, 102, !赋值给支撑NSEL, S, LOC, Z-25 !约束墙底竖向位移 D, ALL, , , , , , UZ, , , , , CMSEL, S, SPRING !选择土弹簧单元 NSLE, S !以下命令从已选弹簧集合中选 CMSEL, S, WALL !在已选集中选墙单元NSLE, U !从集合中删除刚选择的单元，即与墙有关的单元 D, ALL, , , , , , ALL, , , , , !约束土弹簧单元端点的所有位移 ALLSSAVE/SOLUALLSANTYPE, STATIC, NEW !分析模式 静力NROPT, FULL !指定计算模式*AFUN, DEG !指定角度单位为 度Q=2.0E4 !Q为超载GAMA1=0.9E4 !浮重度FAI1=20 !内摩擦角C1=10e3 !粘聚力m1=1500e3 !比例系数b=1 !单元宽度h=1 !墙体单元高度 !无支撑开挖1m 坑内水位-1.0 坑外0NSEL, S, LOC, Z, 0, -0.9 !0到0.9的位置即开挖面以上 ESLN, S !以下命令在当前集里选择CMSEL, R, wall !在当前集选wll 单元NSLE, S !以下命令在当前集里选择 *GET, ZMIN, NODE, , MNLOC, Z, , , , !墙最浅节点处的位置*GET, ZMAX, NODE, , MXLOC, Z, , , , !墙最深节点处的位置LOCZ1=abs (ZMAX) !取绝对值LOCZ2=abs (ZMIN)KA= (TAN(45.0 - FAI1/2)*2 !主动土压力系数 !采用水土分算 PA1=(Q + GAMA1*LOCZ1)*KA-2.0*C1*SQRT(KA) !最深主动土压力公式* IF , PA1, LT, 0, THEN !去除小于零值的可能PA1=0*ENDIFPA2=10*1E3*LOCZ1 !水压力计算PA=PA1+PA2 !总应力SA1=(Q+GAMA1*LOCZ2)*KA-2*C1*SQRT(KA) SA2=10*1E3*LOCZ2WA=SA2SA=SA1+SA2pressure=SA SLZER=ZMAXSLOPE=(SA-PA)/(LOCZ2-LOCZ1) !增长率SFGRAD, PRES, 0, z, SLZER, SLOPE !沿z方向从-1m开始的面力减少SFE, all, 1, PRES, , -PA, , , ! 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, SLOPESFE, all, 1, PRES, , -pressure, , , !alls !将集合扩大到全集cmsel, S, strut !选择内支撑 全部杀死EKILL, ALLallsNSEL, S, LOC, Z, 0.1, -1.1ESLN, SESEL, R, TYPE, , 2 !选择单元类型号为2的单元 EKILL, ALL*do, i, 2, 25z=i-1 !开挖了1m 要在被动区减去相应的土压力K=m1*z*b*hR, i, K, , ,*enddoTIME, 1 !载荷步1allssolveNSEL, S, LOC, Z, 0, -4.9ESLN, SCMSEL, R, wallNSLE, S*GET, ZMIN, NODE, , MNLOC, Z, , , ,*GET, ZMAX, NODE, , MXLOC, 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, THEN PA1=0*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, SLOPESFE, all, 1, PRES, , -PA, , ,allsNSEL, S, LOC, Z, -5.1, -24.9ESLN, SGMSEL, R, wallNSLE, S*GET, ZMIN, NODE, , MNLOC, Z, , , ,*GET, ZMAX, NODE, , MXLOC, Z, , , ,S