《ANSYS基坑弹性地基梁全程序即详解.doc》由会员分享，可在线阅读，更多相关《ANSYS基坑弹性地基梁全程序即详解.doc（8页珍藏版）》请在三一文库上搜索。
1、(/ =20kPaAGLUE, all!粘结所有面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
2、, 0, L2, 0, -h,KWPAVE, 3!将工作平面原点定义在3号点RECTNG, 0, L2, 0, -h,AGLUE, all!粘结所有面ET, 1, SHELL43 元）!ET ，ITYPE，Ename，KOPT1，KOPT6 ，INOPR （定义单ET, 2, COMBIN14!K0PT1K0PT6为元素特性编码!shell43 4 节点塑性大应变单元 !C0MBIN14 弹簧 -阻尼器 Spring-DamperMPTEMP,! 删除系统中已存在的温度表MPTEMP, 1, 0! 定义一个温度表MPDATA, EX, 1, , 2.4E10! 指定与温度相应的材料性能数据弹性
3、模量MPDATA, PRXY , 1, , 0.15!主泊松比ESIZE, 1, 0!指定单元边长AMESH, ALL! 划分面生成面单元NSEL, S, L0C, Z, 0ESLN, SNSEL, S, L0C, Z, -1ESLN, UCM, STRUT, ELEM allsCMSEL, U, STRUT CM, WALL, ELEM! 选择一组节点子集创建新集! 选择已选节点上的单元选择 z 坐标值为 -1 的 -! 从已选集中删除此时剩下只支撑板! 将选择集命名 STRUT 生成元件 !all sel 全选! 去除 STRUT 元件! 将选择集命名 wall 生成元件NSEL, S,
4、L0C, X, 0.1, L1-0.1NPL0TNSEL, R, L0C, Y, 0ESLN, SNSEL, S, L0C, Y, 1ESLN, UENSYM, , , , ALLNSEL, S, L0C, Y, 0.1, L2-0.1NPL0TNSEL, R, L0C, X, 0ESLN, SNSEL, S, L0C, X, 1ESLN, UENSYM , , , , ALL! 选择一组节点子集创建新集! 显示节点! 从当前集选择一组节点子集! 从已选集中选择! 从当前集选择一组节点子集! 从已选集中删除! 反转壳单元法线方向! 选择一组节点子集创建新集! 显示节点! 从当前集选择一组节点子
5、集! 从已选集中选择! 从当前集选择一组节点子集! 从已选集中删除! 反转壳单元法线方向! 所有实体进行重新编号ALLS NUMCMP , ALL从1 到 29进行循环! 直接生成节点 *D0, i, 1, L1-1CSYS, 0!激活默认笛卡尔坐标系N, 100000+2*i-1, i, 0, -1N, 100000+2*i, i, 2.5, -1 *enddo! 节点编号 后面为坐标! 节点编号 后面为坐标*DO, i, 1, L1-1! 从1到 29进行循环!Modeling>Creat>Elements>Elem AttributesTYPE, 2MAT , 1REA
6、L, 0ESYS, 0!设置单元类型属性指示器 !MP 命令中的 MAT 即材料性能 ! 材料实常数 !材料坐标系统属性指示器EN, 100000+i, 100000+2*i-1, 100000+2*i! 根据给定的单元号和节点号生成单元*enddoNUMCMP , ALL!所有实体进行重新编号*DO, i, 1, L1-1CSYS, 0N, 100000+2*i-1, i, L2, -1N, 100000+2*i, i, L2-2.5, *enddo! 激活默认笛卡尔坐标系1*DO, i, 1, L1-1TYPE, 2MAT , 1REAL, 0ESYS, 0EN, 100000+i, 10
7、0000+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, 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*
8、i, L1-2.5, i, -1*enddo*DO, i, 1, L2-1TYPE, 2 $ MA T , 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
9、, ALLNUMMRG , NODE, 0.01, 0.01, LOWNUMCMP , ALL!选择所有实体! 节点合并 距离小于 0.01 则同保留编 !所有实体进行重新编号底的点*DO, i, 1, 25NSEL, S, LOC, Z, -1*iESLN, SNSEL, R, TYPE, , 2R, i, m1*i*b*h, , , EMODIF, ALL, REAL, i, *enddo!识别号 +实常数! 从当前集选择一组节点子集 深度-1 以下 !从已选集中选择!选择单元类型号为 2 的单元 !单位面积内受的力随深度增加而增加 !对已存在单元进行修改R, 101, 0.6, 0.6,
10、 0.6, 0.6, ,R, 102, 0.1, 0.1, 0.1, 0.1, ,CMSEL, S, W ALLEMODIF, ALL, REAL, 101,CMSEL, S, STRUTEMODIF, ALL, REAL, 102,!将实常数 101 组赋给墙! 赋值给支撑NSEL, S, LOC, Z-25D, ALL, , , , , , UZ, , , , ,!约束墙底竖向位移CMSEL, S, SPRINGNSLE, SCMSEL, S, W ALLNSLE, UD, ALL, , , , , , ALL, , , , ,!选择土弹簧单元!以下命令从已选弹簧集合中选! 在已选集中选墙
11、单元! 从集合中删除刚选择的单元，即与墙有关的单元!约束土弹簧单元端点的所有位移ALLSSAVE/SOLUALLSANTYPE, STA TIC, NEWNROPT, FULL*AFUN, DEGQ=2.0E4GAMA1=0.9E4FAI1=20C1=10e3m1=1500e3b=1h=1! 分析模式 静力!指定计算模式! 指定角度单位为 度 !Q 为超载!浮重度! 内摩擦角!粘聚力!比例系数!单元宽度!墙体单元高度! 无支撑开挖 1m 坑内水位 -1.0 坑外 0NSEL, S, LOC, Z, 0, -0.9ESLN, SCMSEL, R, wallNSLE, S!0 到 0.9 的位置即
12、开挖面以上 !以下命令在当前集里选择 ! 在当前集选 wll 单元 !以下命令在当前集里选择*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, THENPA1=0!去除小于
13、零值的可能!水压力计算! 总应力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,乙SLZER, SLOPE !沿z方向从-1m开始的面力减少 SFE, all, 1, PRES, , -PA, , ,!alls! 开挖面以下加载NSEL, S, LOC, Z, -1.1, -24.9ESLN, SGMSEL,
14、 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, A
15、LL*do, i, 2, 25z=i-1!开挖了 1m 要在被动区减去相应的土压力K=m1*z*b*h! 载荷步 1R, i, K, , , *enddoTIME, 1 alls solveNSEL, 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
16、.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=SA SLZER=ZMAX SLOPE=(SA-PA)/(LOCZ2-LOCZ1) SFGRAD, PRES, 0, z, SLZER, SLOPE SFE, all, 1, PRES, , -PA, , , allsNSEL, S, LOC, Z, -5.1, -24.9ESLN, SGM
17、SEL, 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, , ,allscmsel, s, strutealive, ALLallsNSEL, S, LOC, Z, 0.1, -5.12 的单元ESLN, SESEL, R, TYPE, , 2!选择单元类型号为EKILL, ALL*do, i, 6, 25z=i-5K=m1*z*b*hR, i, K, , ,*enddoTIME, 2allssolve