ANSYS基坑弹性地基梁全程序即详解word文档良心出品

1、i/-2(3tPaAGLUE, all!粘结所有面1-r 规 y.I I , LJ：滋g枯I . 00 vJi N CreatElementsElem Attributes!设置单元类型属性指示器!MP 命令中的 MAT 即材料性能!材料实常数!材料坐标系统属性指示器!根据给定的单元号和节点号生成单元NUMCMP , ALL!所有实体进行重新编号!激活默认笛卡尔坐标系*DO, i, 1, L1-1CSYS, 0N, 100000+2*i-1, i, L2, -1N, 100000+2*i, i, L2-2.5, -1 *enddo *DO, i, 1, L1-1 TYPE, 2MAT , 1

2、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-1 TYPE, 2MAT , 1REAL, 0ESYS, 0EN, 100000+i, 100000+2*i-1, 100000+2*i *enddoNUMCMP , ALL!所有实体进行重新编号*DO, i, 1, L2-1CSYS, 0N,100

3、000+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 EGEN, 25, 100000, ALL, , , ,ESEL, S, TYPE, , 2CM, SPRING, ELEM! 选择单元类型号为 2 的单元, , , -1,! 选择单元类型号为 2 的单元! 生成一

4、个叫 SPRING 由单元组成的元件ALLSALLSEL, ALLNUMMRG, NODE, 0.01, 0.01, NUMCMP, ALLLOW!选择所有实体!节点合并 距离小于 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,

5、101, 0.6, 0.6, 0.6, 0.6, ,R, 102, 0.1, 0.1, 0.1, 0.1, ,CMSEL, S, WALL! 识别号 +实常数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, , , ,!选择土弹簧单元!以下命

6、令从已选弹簧集合中选!在已选集中选墙单元!从集合中删除刚选择的单元，即与墙有关的单元!约束土弹簧单元端点的所有位移ALLSSAVE /SOLUALLSANTYPE, STATIC, NEWNROPT, FULL*AFUN, DEGQ=2.0E4GAMA1=0.9E4FAI1=20C1=10e3 m1=1500e3b=1h=1! 分析模式 静力 !指定计算模式 ! 指定角度单位为 度 !Q 为超载!浮重度!内摩擦角!粘聚力!比例系数!单元宽度! 墙体单元高度!无支撑开挖1m 坑内水位 -1.0 坑外 0NSEL, S, LOC, Z, 0, -0.9ESLN, SCMSEL, R, wallNS

7、LE, 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/2)*2!取绝对值!主动土压力系数!采用水土分算最深主动土压力公式PA1=(Q + GAMA1*LOCZ1)*KA-2.0*C1*SQRT(KA) !* IF , PA1, LT, 0,

8、 THENPA1=0!去除小于零值的可能*ENDIFPA2=10*1E3*LOCZ1!水压力计算PA=PA1+PA2!总应力SA1=(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, SLOPE SFE, all, 1, PRES, , -PA, , ,!alls! 增长率! 沿 z 方向从 -1m 开始的面力减少!开挖面以下加载NSEL, S, LOC, Z, -

9、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, , ,!allscmsel, S, strutEKILL, ALL!将集合扩大到全集!选择内支撑 全部杀死allsNSEL, S, LOC, Z, 0.1, -1.1!选择单元类型号为 2 的单元ESLN, SESEL,

10、 R, TYPE, , 2EKILL, ALL*do, i, 2, 25z=i-1! 开挖了 1m 要在被动区减去相应的土压力K=m1*z*b*hR, i, K, , , *enddo! 载荷步 1TIME, 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

11、=(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, SLOPE SFE, 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, , , , SLZER=ZMAXSLOPE=0SFGRAD, PRES, O, z, SLZER, SLOPESFE, all, 1, PRES, , -pressure, ,