《ANSYS基坑弹性地基梁全程序即详解》由会员分享，可在线阅读，更多相关《ANSYS基坑弹性地基梁全程序即详解（8页珍藏版）》请在金锄头文库上搜索。
1、l/20klJa/prep7L1=30L2=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,S
2、HELL43元）ET,2,COMBIN14!ET,ITYPE,Ename,K0PT1,，KOPT6,INOPR（定义单!KOPT1KOPT6为元素特性编码!shell434节点塑性大应变单元!COMBIN14弹簧-阻尼器Spring-DamperMPTEMP,!删除系统中已存在的温度表MPTEMP,1,0!定义一个温度表MPDATA,EX,1,2.4E10MPDATA,PRXY,1,0.15!指定与温度相应的材料性能数据弹性模量!主泊松比ESIZE,1,0AMESH,ALLNSEL,S,LOC,Z,0ESLN,SNSEL,S,LOC,Z,-1ESLN,UCM,STRUT,ELEMallsCMS
3、EL,U,STRUTCM,WALL,ELEM!指定单元边长!划分面生成面单元!选择一组节点子集创建新集!选择已选节点上的单元!选择z坐标值为-1的-!从已选集中删除此时剩下只支撑板!将选择集命名STRUT生成元件!allsel全选!去除STRUT元件!将选择集命名wall生成元件!所有实体进行重新编号!直接生成节点!从1到29进行循环NSEL,S,LOC,X,0.1,L1-0.1NPLOTNSEL,R,LOC,Y,0ESLN,SNSEL,S,LOC,Y,1ESLN,UENSYM,ALLNSEL,S,LOC,Y,0.1,L2-0.1NPLOTNSEL,R,LOC,X,0ESLN,SNSEL,S,
4、LOC,X,1ESLN,UENSYM,ALLALLSNUMCMP,ALL*DO,i,1,L1-1!选择一组节点子集创建新集!显示节点!从当前集选择一组节点子集!从已选集中选择!从当前集选择一组节点子集!从已选集中删除!反转壳单元法线方向!选择一组节点子集创建新集!显示节点!从当前集选择一组节点子集!从已选集中选择!从当前集选择一组节点子集!从已选集中删除!反转壳单元法线方向CSYS,0!激活默认笛卡尔坐标系N,100000+2*i-1,i,0,-1N,100000+2*i,i,2.5,-1*enddo!节点编号后面为坐标!节点编号后面为坐标*DO,i,1,L1-1!从1到29进行循环!Mode
5、lingCreatElementsElemAttributesTYPE,2MAT,1REAL,0ESYS,0!设置单元类型属性指示器!MP命令中的MAT即材料性能!材料实常数!材料坐标系统属性指示器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,2MAT,1REAL,0ESYS,0EN,100
6、000+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,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-
7、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,
8、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,0.6,0.6,R,102,0.1,0.1,0.1,0.1,CMSEL,S,WALLEMODIF,ALL,REAL,101,CMSEL,S,STRUTEMODIF,ALL,REAL,102,!将实常数101组赋给墙!赋值给支撑NSEL,S,LOC,Z-25D,ALL,UZ,!约束墙底
9、竖向位移CMSEL,S,SPRINGNSLE,SCMSEL,S,WALLNSLE,UD,ALL,ALL,!选择土弹簧单元!以下命令从已选弹簧集合中选!在已选集中选墙单元!从集合中删除刚选择的单元，即与墙有关的单元!约束土弹簧单元端点的所有位移ALLSSAVE/SOLUALLSANTYPE,STATIC,NEWNROPT,FULL*AFUN,DEGQ=2.0E4GAMA1=0.9E4FAI1=20C1=10e3m1=1500e3b=1h=1!分析模式静力!指定计算模式!指定角度单位为度!Q为超载!浮重度!内摩擦角!粘聚力!比例系数!单元宽度!墙体单元高度!无支撑开挖1m坑内水位-1.0坑外0NS
10、EL,S,LOC,Z,0,-0.9!0到0.9的位置即开挖面以上ESLN,SCMSEL,R,wallNSLE,S!以下命令在当前集里选择!在当前集选wll单元!以下命令在当前集里选择*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,THENPA1=0!去除小于零
11、值的可能!水压力计算!总应力*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方向从-lm开始的面力减少SFE,all,1,PRES,-PA,!alls!开挖面以下加载NSEL,S,LOC,Z,-l.l,-24.9ESLN,SGMSEL,R,wallNSLE,S*GET,ZMIN,NO
12、DE,MNLOC,Z,*GET,ZMAX,NODE,MXLOC,Z,SLZER=ZMAXSLOPE=0SFGRAD,PRES,O,z,SLZER,SLOPESFE,all,l,PRES,-pressure,!alls!将集合扩大到全集cmsel,S,strut!选择内支撑全部杀死EKILL,ALLallsNSEL,S,LOC,Z,0.l,-l.lESLN,SESEL,R,TYPE,2!选择单元类型号为2的单元EKILL,ALL*do,i,2,25z=i-l!开挖了lm要在被动区减去相应的土压力K=m1*z*b*hR,i,K,*enddoTIME,1!载荷步1allssolveNSEL,S,LO
13、C,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,THENPA1=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,SLZER=ZMAXSLOPE=0SFGRAD,PRES,O,z,SLZER,SLOPESFE,all,1,PRES,-pressure,
