最新ANSYS基坑弹性地基梁全程序即详解汇总.docx
《最新ANSYS基坑弹性地基梁全程序即详解汇总.docx》由会员分享,可在线阅读,更多相关《最新ANSYS基坑弹性地基梁全程序即详解汇总.docx(11页珍藏版)》请在冰豆网上搜索。
![最新ANSYS基坑弹性地基梁全程序即详解汇总.docx](https://file1.bdocx.com/fileroot1/2022-12/8/8f2c431f-b11b-489b-a9da-b1709cbf003a/8f2c431f-b11b-489b-a9da-b1709cbf003a1.gif)
最新ANSYS基坑弹性地基梁全程序即详解汇总
ANSYS
基坑弹性地基
梁全程序即详解
(厂20kPa
-rMW▼JU▽
土体参数:
连续墙参数:
▼从N-
^t-19kN/m3
F-2.4Xld°Pa
c-10kPa
y-0.15
E0=20
宀w-1500kN/m4
•
•
1楼板歩数:
地卜水位:
•
厚度0.M
坑外水位0m不变
v/5-2.4XlOl0Pa
境内水位位于开挖
;0.6m|
4
J55.
面
•VZ././//<•.*///•;*,*zz•.,•..
开挖参数:
第一步:
无支撐开挖至Jn第二步:
支设楼板支撐,并开挖至
/prep7
Ll=30!
设置变量
L2=30
h=-25
K,1,0,0,0,
K,2,Ll,0,0,
K,3,L1,L2,0,
K,4,0,L2,0,
KWPAVE,1!
将工作平面原点定义在1号点
RECTNG,0,Ll,0,L2,
wpro,,-90,!
将工作平面绕X轴Z到Y方向90度
RECTNG,0,Ll,0,-h,
KWPAVE,4!
将工作平面原点定义在4号点
RECTNG,0,Ll,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
(定义单元)
!
KOPT1~KOPT6为元素特性编码
!
shell434节点塑性大应变单元
ET,2,COMBIN14!
COMBIN14弹簧-阻尼器Spring-Damper
MPTEMP,,,,,,,,!
删除系统中已存在的温度表
MPTEMP.1,0!
定义一个温度表
MPDATA,EX,1,,2.4E10!
指定与温度相应的材料性能数据弹性模量
MPDATA,PRXY,1,,0.15!
主泊松比
ESIZE,1,0!
指定单元边长
AMESH,ALL!
划分面生成面单元
NSEL,S,LOC,乙0!
选择一组节点子集创建新集
ESLN,S!
选择已选节点上的单元
NSEL,S,LOC,乙-1!
选择z坐标值为-1的一
ESLN,U!
从已选集中删除此时剩下只支撑板
CM,STRUT,ELEM!
将选择集命名STRUT生成元件
alls!
allsei全选
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丄L2-0.1!
选择一组节点子集创建新集
NPLOT
!
显示节点
NSEL,R,LOC,X,0
!
从当前集选择一组节点子集
ESLN,S
!
从已选集中选择
NSEL,S,LOC,X,1
!
从当前集选择一组节点子集
ESLN,U
!
从已选集中删除
ENSYM,,,,ALL
!
反转壳单元法线方向
ALLS
NUMCMP,ALL
!
所有实体进行重新编号
!
直接生成节点
水DO,i,1,LM
!
从1到29进行循环
CSYS,0
!
激活默认笛卡尔坐标系
N,100000+2怜1,i,0.-1!
节点编号后面为坐标
N,100000+2%i,2.5,-1!
节点编号后面为坐标
*enddo
水DO,i,1,LM
!
从1到29进行循环
!
Modeling>Creat>EIements>ElemAttributes
TYPE,2
!
设置单元类型属性指示器
MAT,1
!
MP命令中的MAT即材料性能
REAL,0
!
材料实常数
ESYS,0
!
材料坐标系统属性指示器
EN,100000+i,100000+2怜1,100000+2*i!
根据给定的单元号和节点号生成单
元
*enddo
NUMCMP,ALL!
所有实体进行重新编号
*DO,i,1,L1-1
CSYS,0
!
激活默认笛卡尔坐标系
N,100000+2*i-l,i,L2,-l
N,100000+2*i,i,L2-2.5,-l
*enddo
*DO,i,1,L1-1
TYPE,2
MAT,1
REAL,0
ESYS,0
EN,100000+i,100000+2*i-l,100000+2*i*enddo
NUMCMP.ALL!
所有实体进行重新编号
*DO.i,1,L2-1
CSYS,0
N,100000+2*i-l,0,i,-l
N,100000+2*i,2.5,i,-l
*enddo
*DO,i,1,L2-1
TYPE,2
MAT,1
REAL,0
ESYS,0
EN,100000+i,100000+2*i-l,100000+2*i*enddo
NUMCMP,ALL!
所有实体进行重新编号
*DO,i,1,L2-1
CSYS,0
N,100000+2*i-l,LI,i,-l
N,100000+2%LI-2.5,i,-1
*enddo
*DO,i,1,L2-1
TYPE,2$MAT,1$REAL,0
ESYS,0
EN,100000+i,100000+2*i-l,100000+2*i*enddo
NUMCMP.ALL!
所有实体进行重新编号
ESEL,S,TYPE,,2!
选择单元类型号为2的单元
EGEN,25,100000,ALL,
ESEL,S,TYPE,,2!
选择单元类型号为2的单元
CM,SPRING,ELEM!
生成一个叫SPRING由单元组成的元件
ALLS
ALLSEL,ALL!
选择所有实体
NUMMRG,NODE,0.01,0.01,LOW!
节点合并距离小于0.01则同保留编号
底的点
NUMCMP,ALL!
所有实体进行重新编号
*DO,i,1,25
NSEL,S,LOC,乙-l*i
!
从当前集选择一组节点子集深度-1以下
ESLN,S
!
从已选集中选择
NSEL,R,TYPE,,2
!
选择单元类型号为2的单元
R,i,ml*i*b*h,,,
!
单位面积内受的力随深度増加而増加
EMODIEALL,REAL,i,!
对已存在单元进行修改
*enddo
NSEL,S,LOC,Z-25!
约束墙底竖向位移
D,ALL,.U乙…,
CMSEL,S,SPRING
!
选择土弹簧单元
NSLE,S
!
以下命令从已选弹簧集合中选
CMSEL,S,WALL
!
在已选集中选墙单元
NSLE,U
!
从集合中删除刚选择的单兀,即与墙有关的单兀
ALLS
SAVE
/SOLU
ALLS
!
无支撑开挖Im坑内水位・1・0坑外0
NSEL,S,LOC,乙0,-0.9!
0到0.9的位置即开挖面以上
*GET,ZMIN,NODE…MNLOC,乙…!
墙最浅节点处的位置
*GET,ZMAX,NODE,MXLOC,乙…!
墙最深节点处的位置
LOCZl=abs(ZMAX)!
取绝对值
LOCZ2=abs(ZMIN)
KA=(TAN(45.0-FAIl/2))**2!
主动土压力系数
!
采用水土分算
PA1=(Q+GAMA1*LOCZ1)*KA-2.0*C1*SQRT(KA)!
最深主动土压力公式
SA1=(Q+GAMA1*LOCZ2)*KA-2*C1*SQRT(KA)
SA2=10*1E3*LOCZ2
WA=SA2
SA=SA1+SA2
pressure=SA
SLZER=ZMAX
SLOPE=(SA-PA)/(LOCZ2-LOCZ1)!
增长率
SFGRAD,PRES,0,z,SLZER,SLOPE!
沿z方向从-lm开始的面力减少
SFE,all,1,PRES,,-PA,,,!
alls
!
开挖面以下加载
NSEL,S,LOC,Z,-1.1,-24.9
ESLN,S
GMSEL,R,wall
NSLE,S
*GET,ZMIN,NODE,,MNLOC,乙,,,
*GET,ZMAX,NODE,,MXLOC,乙…
SLZER二ZMAX
SLOPE=0
SFGRAD,PRES,O,乙SLZER,SLOPE
SFE,all,1,PRES,,-pressure,,,!
alls!
将集合扩大到全集
cmsel,S,strut!
选择内支撑全部杀死
EKILL,ALL
alls
NSEL,S,LOC,乙0.1,-1.1
ESLN,S
ESEL,R.TYPE,,2!
选择单元类型号为2的单元
EKILL,ALL
*do,i,2,25
z=i-l!
开挖了lm要在被动区减去相应的土压力
K=ml*z*b*h
R,i,K,,,
*enddo
TIME,1!
载荷步1
alls
solve
NSEL,S,LOC,乙0,-4.9
ESLN,S
CMSEL,R,wall
NSLE,S
*GET,ZMIN,NODE,,MNLOC,Z,,,,
*GET,ZMAX,NODE,,MXLOC,Z,,,,
LOCZl=abs(ZMAX)
LOCZ2=abs(ZMIN)
KA=(TAN(45.0-FAIl/2))**2
PA1=(Q+GAMA1*LOCZ1)*KA-2.0*C1*SQRT(KA)
*IF,PA1,LT,0,THEN
PA1=O
*ENDIF
PA2=10*1E3*LOCZ1
PA=PA1+PA2
SA1=(Q+GAMA1*LOCZ2)*KA-2*C1*SQRT(KA)
SA2=10*1E3*LOCZ2
WA=SA2
SA=SA1+SA2
pressure=SA
SLZER=ZMAX
SLOPE=(SA-PA)/(LOCZ2-LOCZ1)
SFGRAD,PRES,0,z,SLZER,SLOPE
SFE,all,1,PRES,,-PA,,,
alls
NSEL,S,LOC,Z,-5.1,-24.9
ESLN,S
GMSEL,R,wall
NSLE,S
*GET,ZMIN,NODE,,MNLOC,乙,,,
*GET,ZMAX,NODE,,MXLOC,乙,,,
SLZER=ZMAX
SLOPE=0
SFGRAD,PRES,O,乙SLZER,SLOPE
SFE,all,1,PRES.,-pressure,,,
alls
cmsel,s,strut
ealive,ALL
alls
NSEL,S,LOC,Z,0.1,-5.1
ESLN,S
ESEL,R,TYPE,,2!
选择单元类型号为2的单元
EKILL,ALL
*do,i,6,25
z=i-5
K=ml*z*b*h
R,i,K…*enddo
TIME,2allssolve