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

/prep7

L1=30 !设置变量

L2=30

h=-25

K, 1, 0, 0, 0,

K, 2, L1, 0, 0,

K, 3, L1, L2, 0,

K, 4, 0, L2, 0,

KWPA VE, 1 !将工作平面原点定义在1号点RECTNG, 0, L1, 0, L2,

wpro, , -90, !将工作平面绕X轴Z到Y方向90度RECTNG, 0, L1, 0, -h,

KWPA VE, 4 !将工作平面原点定义在4号点RECTNG, 0, L1, 0, -h,

wpro, , ,90 !将工作平面绕y轴x到z方向90度RECTNG, 0, L2, 0, -h,

KWPA VE, 3 !将工作平面原点定义在3号点RECTNG, 0, L2, 0, -h,

AGLUE, all !粘结所有面

ET, 1, SHELL43 !ET，ITYPE，Ename，KOPT1，~，KOPT6，INOPR（定义单元）

!KOPT1～KOPT6为元素特性编码

!shell43 4 节点塑性大应变单元

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, 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, W ALL, 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 !反转壳单元法线方向

ALLS

NUMCMP, 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进行循环

!Modeling>Creat>Elements>Elem Attributes

TYPE, 2 !设置单元类型属性指示器

MAT , 1 !MP命令中的MA T即材料性能

REAL, 0 !材料实常数

ESYS, 0 !材料坐标系统属性指示器

EN, 100000+i, 100000+2*i-1, 100000+2*i !根据给定的单元号和节点号生成单元*enddo

NUMCMP , ALL !所有实体进行重新编号

*DO, i, 1, L1-1

CSYS, 0 !激活默认笛卡尔坐标系

N, 100000+2*i-1, i, L2, -1

N, 100000+2*i, i, L2-2.5, -1

*enddo

*DO, i, 1, L1-1

TYPE, 2

MAT , 1

REAL, 0

ESYS, 0

EN, 100000+i, 100000+2*i-1, 100000+2*i

*enddo

NUMCMP , ALL !所有实体进行重新编号

*DO, i, 1, L2-1

CSYS, 0

N, 100000+2*i-1, 0, i, -1

N, 100000+2*i, 2.5, i, -1

*enddo

*DO, i, 1, L2-1

TYPE, 2

MAT , 1

REAL, 0

ESYS, 0

EN, 100000+i, 100000+2*i-1, 100000+2*i

*enddo

NUMCMP , ALL !所有实体进行重新编号

*DO, i, 1, L2-1

CSYS, 0

N,100000+2*i-1, L1, i, -1

N,100000+2*i, L1-2.5, i, -1

*enddo

*DO, i, 1, L2-1

TYPE, 2 $ MAT , 1 $ REAL , 0

ESYS, 0

EN, 100000+i, 100000+2*i-1, 100000+2*i

*enddo

NUMCMP , ALL !所有实体进行重新编号

ESEL, S, TYPE, , 2 !选择单元类型号为2的单元

EGEN, 25, 100000, ALL, , , , , , , , , , -1,

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, 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, WALL