ANSYS在结构中 的仿真模拟与分析.docx

1 楼层浇筑过程仿真分析

/PREP7
!定义单元类型、实常数、截面参数
ET,1,BEAM188
ET,2,SHELL63
MP,EX,1,3.25E10 !定义 C40 混泥土材料属性
MP,PRXY,1,0.17
MP,DENS,1,2800
MP,EX,2,3.0E10 !定义 C30 混泥土材料属性
MP,PRXY,2,0.17
MP,DENS,2,2800
MP,EX,3,3.0E10
MP,PRXY,3,0.167
MP,DENS,3,2800
SECTYPE,1,BEAM,RECT !框架柱截面特性，设定为矩形截面梁
SECDATE,1,1 ！梁宽 1，高 1
SECTYPE,2,BEAM,RECT !外框架梁截面特性
SECDATA,0.3,0.5
SECTYPE,3,BEAM,RECT !内框架梁截面特性
SECDATA,0.4,0.7
SECTYPE,4,BEAM,RECT !次梁截面特性
SECDATE,0.2,0.4

R,2,0.2
R,3,0.3

!建模
K,3000,22,8,72
k,4000,1000
K,,

!定义楼板及外墙厚度
!定义筒体单元实常数

K,12,44
KFILL,1,12
KGEN,5,1,12,,,4
KGEN,11,1,60,,,,3
!建立框架柱模型
*DO,I,1,541,60
L,I,I+60
*ENDDO

LGEN,3,1,10,1,8,,,2
LGEN,2,1,30,1,,16,,48
LGEN,2,1,10,1,,8,,24
LGEN,2,21,30,1,,8,,24
LGEN,2,1,80,1,28,,,7
LATT,1,,1,,4000,,1
LESIZE,ALL,2
LMESH,ALL
LSEL,U,,,ALL

!建立外环梁模型
L,61,65
L,65,89
L,89,92
L,92,68
L,68,72
L,72,120
L,120,109
L,109,61
LGEN,10,161,168,1,,,3,60
LATT,1,,1,,4000,,2
LESIZE,ALL,2
LMESH,ALL
LSEL,U,,,ALL

!建立内框架梁模型
L,63,111
L,85,89
L,89,113
L,70,118
L,92,96
L,92,116
LGEN,10,241,246,1,,,3,60
LATT,1,,1,,3000,,3
LESIZE,ALL,2
LMESH,ALL
LSEL,U,,,ALL

！建立次梁模型
L,62,110
L,64,112
L,69,117

L,71,119
L,97,108
L,73,77
L,80,84
L,90,102
L,91,103
LGEN,10,301,309,1,,,3,60
LATT1,,1,,3000,,4
LESIZE,ALL,2
LMESH,ALL
ALLS

!建立楼层板模型
*DO,I,61,64
A,I,I+1,I+13,I+12
*ENDDO
AGEN,4,1,4,1,,4,,12

*DO,I,68,71,1
A,I,I+1,I+13,I+12
*ENDDO
AGEN,4,17,20,1,,4,,12

*DO,I,89,91,1
A,I,I+1,I+13,I+12
*ENDDO
AGEN,10,1,35,1,,,3,60
A,653,656,644,641

AATT,2,2,2
AESIZE,ALL,2

AMESH,ALL
ASEL,U,,,ALL

!建立外墙模型并划分网格
*DO,I,61,64,1
A,I,I+1,I+61,I+60
*ENDDO

*DO,I,89,91,1
A,I,I+1,I+61,I+60
*ENDDO

*DO,I,68,71,1
A,I,I+1,I+61,I+60
*ENDDO

*DO,I,109,112,1
A,I,I+1,I+61,I+60
*ENDDO

*DO,I,116,119,1
A,I,I+1,I+61,I+60
*ENDDO

*DO,I,61,97,12
A,I,I+12,I+72,I+60
*ENDDO

*DO,I,65,77,12
A,I,I+12,I+72,I+60
*ENDDO

*DO,I,68,80,12
A,I,I+12,I+72,I+60
*ENDDO

*DO,I,72,108,12
A,I,I+12,I+72,I+60
*ENDDO

AGEN,9,352,382,,,,3,60
AATT,3,2,2
AESIZE,ALL,2
AMESH,ALL
ASEL,U,,,ALL

!建立筒体模型
A,53,56,656,653
A,41,44,644,641

AATT,1,3,2
ESIZE,1,0
MSHAPE,0,2D
MSHKEY,1
AMESH,ALL
ASEL,U,,,1

A,41,53,653,641
A,44,56,656,644

AATT,1,3,2
ESIZE,1,0
MSHAPE,0,2D
MSHKEY,1
AMESH,ALL
NUMMRG,ALL
NUMCMP,ALL

ALLSEL
SAVE
FINISH
!至此，模型已经全建立完毕
！进入求解模块
/SOLU
ANtYPE,STATIC
DELTIM,0.1,0.05,0.2
AUTOTS,ON
NLGEOM,ON
PRED,ON
LNSRCH,ON
NSEL,S,LOC,Z,0
D,ALL,ALL
ACEL,,,10

!开始施工过程的模拟

ESEL,ALL
EKILL,ALL
ESEL,ALL
NSEL,ALL
D,ALL,ALL

!首先修建柱子和梁，同事再修建楼板
！将柱、梁和楼板所属的单元激活，并删除其全部约束
ESEL,S,TYPE,,1
NSLE,S
NSEL,R,LOC,Z,0,3
ESLN,R,1

EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
ALLS
ESEL,S,MAT,,2
NSLE,S
NSEL,R,LOC,Z,0,3
ESLN,R,1
EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
ALLS
NSEL,R,LOC,Z,0
D,ALL,ALL
NSEL,ALL
ESEL,ALL
SOLVE
!激活第二至第十层的柱、梁和楼板单元
*DO,I,2,10,1
ESEL,S,TYPE,,1
NSLE,S
NSEL,R,LOC,Z,3*(I-1),3*I
ESLN,R,1
EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
ALLS
ESEL,S,MAT,,2
NSLE,S
NSEL,R,LOC,Z,3*(I-1),3*I
ESLN,R,1
EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
NSEL,ALL
ESEL,ALL
SOLVE
*ENDDO

!至此，梁、柱和楼板施工完毕
!以下进行外墙及隔墙的施工

ESEL,S,MAT,,3
NSLE,S
NSEL,R,LOC,Z,0,3
ESLN,R,1
EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
ALLS
ESEL,S,MAT,,1
NSLE,S
NSEL,R,LOC,Z,0,3
ESLN,R,1
EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
ALLS
NSEL,S,LOC,Z,0
DDELE,ALL,ALL
NSEL,ALL
ESEL,ALL
SOLVE

!激活第二至第十层的墙面单元
*DO,I,2,10,1
ESEL,S,MAT,,3
NSLE,S
NSEL,R,LOC,Z,3*(I-1),3*I
ESLN,R,1
EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
ALLS
ESEL,S,MAT,,1
NSLE,S
NSEL,R,LOC,Z,3*(I-1),3*I
ESLN,R,1
EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
NSEL,ALL
ESEL,ALL
SOLVE
*ENDDO
FINISH
!求解过程结束
!进入通用后处理器查看结果
/POST1
PLNSOL,U,X
PLNSOL,U,y
PLNSOL,U,z
PLNSOL,S,X
PLNSOL,S,Y
PLNSOL,S,Z

!查看柱子和梁的轴力，剪力和弯矩

ESEL,S,TYPE,,1 ETABLE,IF,SMISC,1 ETABLE,JF,SMISC,7 ETABLE,II,SMISC,2 ETABLE,JI,SMISC,8 ETABLE,IM,SMISC,6 ETABLE,JM,SMISC,12

PLETAB,IF,NOAV
PLETAB,JF,NOAV
PLETAB,II,NOAV
PLETAB,JI,NOAV
PLETAB,IM,NOAV
PLETAB,JM,NOAV
FINISH


!定义轴力单元列表 1
!定义轴力单元列表 2
!定义剪力单元列表 1
!定义剪力单元列表 2
!定义弯矩单元列表 1
!定义弯矩单元列表 2
!显示轴力单元列表 1
!显示轴力单元列表 2
!显示剪力单元列表 1
!显示剪力单元列表 2
!显示弯矩单元列表 1
!显示弯矩单元列表 2

2.斜拉桥三维仿真分析

/filname,cable-stayed bridge,1
keyw,pr_struc,1
/prep7
!定义单元类型
et,1,beam4

et,2,link10
!定义材料属性
mp,ex,1,3.5e10
mp,prxy,1,0.17
mp,dens,1,2500
mp,ex,2,10e15
mp,prxy,2,0
mp,dens,2,0
mp,ex,3,1.9e10
mp,prxy,3,0.25
mp,dens,3,1200
mp,damp,3,0.5
!定义实常数
!定义实常数
r,1,25.6,5.46,546.133,16,1.6
r,2,16,29.417,15.394,3.4,4.7
r,3,54,364.5,162,6,9
r,4,40,213.3,83.3,5,8
r,5,1,1/12,1/12,1,1
r,6,0.012,0.012 !索的
!创建节点和单元
!建立主梁节点
/view,1,1,1,1
/angle,1,270,xm,0
/replot

*do,i,1,59 x=-174*2+(i-1)*6 y1=-14
y2=14
n,3*(i-1)+1,x
n,3*(i-1)+2,x,y1
n,3*i,x,y2
*enddo

!此循环用于建立主梁的半跨节点
!最左端 x=174*2， x=0 左边的节点 x 坐标值，间距为 6
!桥面宽 28 米，故左边节点为-14
!桥面宽 28 米，故右边节点为-14
!建立主梁节点 3* (i-1) +1 为节点号
!以下两行建立桥面两边节点
!能想出这种建模命令的绝对是编程高手，哈哈
!完全可以先建立端部的三个节点，然后用这三个节点在 x

方向上复制 59 份，间距为 6
!建立主梁单元
type,1
real,1

mat,1
*do,i,1,58,1
j=3*(i-1)+1
e,j,j+3
*enddo
!建立鱼刺刚横梁
type,1
real,5
mat,2
*do,i,1,59,1
j=3*(i-1)+1


!以下循环建立建立桥面中线主梁单元

!以下循环用于建立桥面鱼刺横梁的节点

j1=3*(i-1)+2
j2=3*i
e,j,j1
e,j,j2
*enddo
!建立半跨主塔
i=59*3 !变量用于记录桥面的节点数， 即至此已经建立了 59*3 个节点了，
用于指导以后设定节点的编号
n,i+1,-174,-10,-30 !以下两行记录塔脚节点
n,i+2,-174,10,-30
n,i+3,-174,-15 !以下两行用于建立与桥面齐高的主塔节点
n,i+4,-174,15
*do,j,1,5,1 !以下循环用于建立索塔在桥面以上的节点
k=i+4+j
n,k,-174,0,60+(j-1)*18
*enddo
!建立下索塔单元
type,1
real,4
mat,1
e,i+1,i+3 !以下用于建立主塔在桥面以下的两根塔柱单元
e,i+2,i+4
!建立中索塔单元
type,1
real,3

mat,1
e,i+3,i+5
e,i+4,i+5
!建立上索塔单元
type,1
real,2
mat,1
*do,j,1,4,1
k=i+4+j
e,k,k+1
*enddo


!以下用于建立倒 Y 分叉点到桥面间的两根塔柱单元

!以下用于建立倒 Y 分叉点以上的塔柱单元

!建立与塔的倒 Y 分叉点链接的索单元
type,2
real,6
mat,3
e,i+5,89
e,i+5,90
!建立主塔倒 Y 分叉点以上第一个张拉点连接的索单元
*do,j,1,8,1
!此循环用于建立主塔倒 Y 分叉点以上第一个张拉点连接的所有索单元，共 32 个
e,i+6,89+3*j
e,i+6,89-3*j
e,i+6,90+3*j
e,i+6,90-3*j
*enddo

!建立与主塔的其他三个张拉点连接的单元
*do,k,1,3,1
*do,j,1,7,1

e,i+6+k,113+(k-1)*21+3*j
e,i+6+k,65-(k-1)*21-3*j
e,i+6+k,114+(k-1)*21+3*j
e,i+6+k,66-(k-1)*21-3*j
*enddo
*enddo
!生成全桥模型节点
i=i+9
nsym,x,i,all
esym,,i,all
nummrg,all
!建立索塔连接横梁单元
type,1
real,5
mat,2
j=i
i=i-9
n,1000,-174
e,1000,i+3
e,1000,i+4
n,2000,174
e,2000,i+3+j

!一共有 28 个索单元连接在每个张拉点上

!记录半跨的所有节点数
!用映射法直接建立另半跨节点
!用映射法直接建立另半跨单元
!合并所有节点和单元

e,2000,i+4+j
!施加主塔的四个脚上的全约束
nsel,s,loc,z,-30
d,all,all
allsel

!在左桥端施加 y,z 约束 nsel,s,loc,x,-348
nsel,r,loc,y,0


!仅给左端主梁施加约束

d,all,uy
d,all,uz
allsel

!在右桥端施加 y 约束 nsel,s,loc,x,348
nsel,r,loc,y,0


!仅给右端主梁施加约束

d,all,uy
allsel

numcmp,all
!施加重力场
acel,,,9.8

!耦合节点 ,耦合跨中由于对称而重复的单元节点以及两主塔上塔横梁和主梁的重合 节点，
cpintf,uy

cpintf,uz
cpintf,rotx
cpintf,rotz

!成桥状态的确定
!静力的初步计算
!直接进行静力计算
/solu
solve
finish
!得到最大位移为 1.288m，由于偏差太大需要重新计算(与事实不符) !修改实常数后重新计算：令 r,6,0.012,0.012,即给索以预应变 0.012
/solu !为了将计算应力用于下面的动力分析，这里打开预应力和集中质量设置
开关
lumpm,on
pstres,on
solve
finish
!求得最大位移为 0.0329m，说明已经达到成桥状态要求,(与事实不符)

!模态分析
!分析设置
/solu
antype,2
!MODOPT,LANB,20
!EQSLV,SPAR
MXPAND,10, , ,0
LUMPM,1
PSTRES,1

MODOPT,LANB,20,0,100, ,OFF
UPCOORD,1,ON !更新模型计算坐标， 目的是为了将预应力效应准确
的应用能够到模态分析上来
!分析设置完毕
solve
finish
!激励耦合分析
!恢复成桥静力计算结果
RESUME, cable-stayedbridge,db, !路径： utility menu>file>resume jobname.db
!生成地震激励向量， 将记事本格式的地震波数据调入到工作目录下,并执行以下命令 *DIM,aay1,ARRAY,2,50,1
*DIM,aaz1,ARRAY,2,50,1

*CREATE,ansuitmp
*VREAD,aay1(1,1),'tjx','txt',' ',50
(e9.3,e11.3)
*END

/INPUT,ansuitmp

*CREATE,ansuitmp
*VREAD,aaz1(1,1),'tjy','txt',' ',50
(e9.3,e11.3)
*END
/INPUT,ansuitmp

!输入车辆激励波， Z=1000*cos(10*t).车子跨越一个单元的时间为 0.1 秒(速度 216 公里/小时)，全桥共 696 米，历时 11.6 秒，纵向共 116 个单元， 117 个节点

n=117
*dim,fcar,array,n
*do,i,1,n,1
写入数据
tt=(i-1)*0.1
fcar(i)=1000*cos(10*tt)
*enddo

!定义向量维数
!定义车载荷向量
!以下循环为生成车载荷数据， 即向 “far”阵中

!输入风载荷激励 p=50*sin(1.5*t)。由于梁长 6 米，高 3 米，故每个节点的受力可表 示为：

*dim,fwind,array,n,n
***
*do,i,1,n,1
tt=(i-1)*0.1

!定义风荷载向量， ***注意要与车辆荷载步即时间相等

!按不同时刻，每一次都是给列赋值，即各节点赋予荷

载值，第 a 行第 b 列表示第b 时刻第 a 个节点的荷载值
p=50*sin(1.5*tt)
fwind(1,i)=3*3*p !两端的每根梁的面积为 6*3，两个节点上所受的力应为一半面
积与压力的乘积： 3*3*p
*do,j,1,n-2
fwind(j+1,i)=6*3*p !中间的每根梁的面积为 6*3，两个节点上所受的力应为两个相 邻梁各一半面积与压力的乘积： 6*3*p
*enddo
fwind(n,i)=3*3*p
*enddo

!加载步耦合
!由于瞬态分析是以荷载步为阶段来求解的，所以应该将 117 个载荷步分别写入荷载 文件中
finish !完成上一阶段的设置
/config,nres,2000 !设置最大子步数为 2000

/solu
antype,trans
trnopt,full
m=n


!定义分析类型为瞬态
!求解方式为完全法
!m 变量表示向量的需要存储的时间点数 m=n=117

*do,i,1,m,1
time,i*0.1
auto,off
deltim,0.05
kbc,0
outres,basic cnvtol,f,,0.01 !施加地震激励

!总循环开始，开始激励的耦合循环
!给每个载荷步文件设置结束点时间
!关闭自动时间步选项
!设置子步求解的时间间隔
!选择 ramped 选项
!选择基本结果输出
!选择力收敛准则，容差为 0.01

*if,i,le,50,then !如果 i 小于等于 50 (le 为 little equal to)
acel,,-aay1(2,i),10-aaz1(2,i) !其中最后一项为重力加速度与地震竖向加速度的差值
*else
acel,0,0,10
*endif
!施加移动载荷激励
*if,i,eq,1,then !如果力作用在最左端的主梁上时，单独定义施加的载荷，目的为 了以下删除载荷的语句可以使用
nsel,s,loc,x,-174*2
nsel,r,loc,y,0
nsel,r,loc,z,0
f,all,fz,-fcar(i)
allsel
*else !如果力不是作用在最左端时

nsel,s,loc,x,-174*2+i*6
nsel,r,loc,y,0
nsel,r,loc,z,0
f,all,fz,-fcar(i)
allsel
nsel,s,loc,x,-174*2+(i-1)*6
nsel,r,loc,y,0
nsel,r,loc,z,0
fdele,all,fz
allsel
*endif

!选择受力节点施加载荷

!删除前一节点上的载荷

!施加风载荷, (子循环)
*do,j,1,m,1
nsel,s,loc,y,-14 !选择有风作用一侧的节点
nsel,r,loc,x,-174*2+(j-1)*6
f,all,fy,fwind(j,i) !给所选节点施加载荷
allsel
*enddo
lswrite,i !写数据到载荷步文件 i
*enddo !总循环结束
lssolve,1,m,1

3.用ANSYS 对钢筋混凝土梁进行计算模拟

!分离式模型
FINISH

/CLEAR
/PREP7 !1.定义单元与材料属性 ET,1,SOLID65,,,,,,,1 ET,2,LINK8

MP,EX,1,13585
MP,PRXY,1,0.2

FC=14.3
FT=1.43

TB,CONCR,1 TBDATA,,0.5,0.95,FT,-1

TB,KINH,1,,11
TBPT,,0.0002,FC*0.19
TBPT,,0.0004,FC*0.36
TBPT,,0.0006,FC*0.51
TBPT,,0.0008,FC*0.64
TBPT,,0.0010,FC*0.75
TBPT,,0.0012,FC*0.84
TBPT,,0.0014,FC*0.91
TBPT,,0.0016,FC*0.96
TBPT,,0.0018,FC*0.99
TBPT,,0.002,FC
TBPT,,0.0033,FC*0.85

MP,EX,2,2.0E5
MP,PRXY,2,0.3
TB,BISO,2
TBDATA,,300,0


!混凝土材料的初始弹模以及泊松比

!混凝土单轴抗压强度和单轴抗拉强度

!定义混凝土材料及相关参数，关闭压碎 !定义混凝土应力应变曲线，用 KINH 模型

!钢筋材料的初始弹模以及泊松比

!钢筋的应力应变关系，用 BISO 模型

PI=ACOS(-1)
R,1
R,2,0.25*PI*22*22
R,3,0.25*PI*10*10
TBPLOT,MISO,1
TBPLOT,BISO,2

!2.创建几何模型
BLC4,,,150,300,2000
*DO,I,1,19
WPOFF,,,100
VSBW,ALL
*ENDDO

WPCSYS,-1
WPOFF,,,950
VSBW,ALL
WPOFF,,,100
VSBW,ALL

WPCSYS,-1
WPROTA,,-90
WPOFF,,,30
VSBW,ALL
WPOFF,,,240
VSBW,ALL

WPCSYS,-1
WPOFF,30
WPROTA,,,90
VSBW,ALL
WPOFF,,,45
VSBW,ALL
WPOFF,,,45
VSBW,ALL
WPCSYS,-1
!3.划分钢筋网格
ELEMSIZ=50
LSEL,S,LOC,X,30
LSEL,R,LOC,Y,30
LATT,2,2,2
LESIZE,ALL,ELEMSIZ
LMESH,ALL


!混凝土材料的数据表绘图
!钢筋材料的数据表绘图

!切出箍筋位置

!切出拟加载面

!网格尺寸变量，设置为 50mm

LSEL,S,LOC,X,75
LSEL,R,LOC,Y,30
LATT,2,2,2
LESIZE,ALL,ELEMSIZ
LMESH,ALL

LSEL,S,LOC,X,120
LSEL,R,LOC,Y,30
LATT,2,2,2
LESIZE,ALL,ELEMSIZ
LMESH,ALL

LSEL,S,LOC,X,30
LSEL,R,LOC,Y,270
LATT,2,3,2
LESIZE,ALL,ELEMSIZ
LMESH,ALL

LSEL,S,LOC,X,120
LSEL,R,LOC,Y,270
LATT,2,3,2
LESIZE,ALL,ELEMSIZ
LMESH,ALL

LSEL,S,TAN1,Z
LSEL,R,LOC,Y,30,270
LSEL,R,LOC,X,30,120
LSEL,U,LOC,X,75
LSEL,U,LOC,Z,0
LSEL,U,LOC,Z,2000
LSEL,U,LOC,Z,950
LSEL,U,LOC,Z,1050
LATT,2,3,2
LESIZE,ALL,ELEMSIZ
LMESH,ALL

LSEL,ALL
!4.划分混凝土网格
VATT,1,1,1
MSHKEY,1
ESIZE,ELEMSIZ
VMESH,ALL

ALLSEL,ALL

!5.施加荷载和约束
LSEL,S,LOC,Y,0
LSEL,R,LOC,Z,100
DL,ALL,,UY
LSEL,S,LOC,Y,0
LSEL,R,LOC,Z,1900
DL,ALL,,UY

DK,KP(0,0,100),UX,,,,UZ
DK,KP(0,0,1900),UX

P0=180000
Q0=P0/150/100
ASEL,S,LOC,Z,950,1050
ASEL,R,LOC,Y,300
SFA,ALL,1,PRES,Q0

ALLSEL,ALL
FINISH
!6.求解控制设置与求解
/SOLU
ANTYPE,0
NSUBST,60
OUTRES,ALL,ALL
AUTOS,ON
NEQIT,20
CNVTOL,U,,0.15
ALLSEL,ALL
SOLVE
FINISH
!7.进入 POST1 查看结果 /POST1
SET,LAST
PLDISP,1

ESEL,S,TYPE,,2
ETABLE,SAXL,LS,1
PLLS,SAXL,SAXL


!设定最后荷载步，查看变形

!绘制钢筋应力图

ESEL,S,TYPE,,1
/DEVICE,VECTOR,ON
PLCRACK,1,1


!绘制裂缝和压碎图

4.排水管道仿真分析
/PREP7

!定义单元类型，实常数和材料参数
ET,1,SOLID45
ET,2,LINK8
R,1,137.4E-6,2.5E-4
MP,EX,1,4E10
MP,PRXY,1,0.17
MP,DENS,1,2500
MP,EX,2,2.0E11
MP,PRXY,2,0.3
MP,DENS,2,7800

!创建模型
CSYS,1
N,1,0.58
N,2,0.63
N,3,0.7
NGEN,2,3,1,3,1,0,18,0,1
NGEN,2,6,1,6,1,0,0,0.05,1

!创建管道单元
TYPE,1
MAT,1
E,1,2,5,4,7,8,11,10
E,2,3,6,5,8,9,12,11
EPLOT
EGEN,20,12,1,2,1,,,,,,,18,,
NUMMRG,ALL
NUMCMP,ALL
EGEN,120,120,1,40,1,,,,,,,,0.05
NUMMRG,ALL
NUMCMP,ALL

!创建钢筋单元
TYPE,2
MAT,2
!创建 19 根箍筋
*DIM,NONE,ARRAY,21
*DO,I,1,19
*DO,J,1,21
NSEL,S,LOC,X,0.63
NSEL,R,LOC,Z,0.3*I
NSEL,R,LOC,Y,18*(J-1)
*GET,NONE(J),NODE,,NUM,MAX
*ENDDO
*DO,J,1,20
E,NONE(J),NONE(J+1)
*ENDDO
*ENDDO
!创建 10 根纵筋单元
*DIM,NTWO,ARRAY,121
*DO,I,1,10
*DO,J,1,121
NSEL,S,LOC,X,0.63
NSEL,R,LOC,Y,36*(I-1)
NSEL,R,LOC,Z,0.05*(J-1)
*GET,NTWO(J),NODE,,NUM,MAX
*ENDDO
*DO,J,1,120
E,NTWO(J),NTWO(J+1)
*ENDDO
*ENDDO
NUMMRG,ALL
NUMCMP,ALL

/SOLU
ANTYPE,0
NSEL,S,LOC,Z,0,0.3
NSEL,R,LOC,X,0.7
D,ALL,ALL
NSEL,S,LOC,Z,5.7,6
NSEL,R,LOC,X,0.7
D,ALL,ALL
ALLSEL

NSEL,S,LOC,X,0.58
SF,ALL,PRES,1E5
ALLSEL
NSEL,S,LOC,Y,0,180
NSEL,R,LOC,X,0.7
SF,ALL,PRES,3E5
CSYS,0
ACEL,,10
ALLSEL
SOLVE

/POST1
CSYS,1
PLNSOL,U,X

PLNSOL,U,Y
PLNSOL,U,Z
PLNSOL,U,SUM

PLNSOL,S,X
PLNSOL,S,Y
PLNSOL,S,1
PLNSOL,S,3

5.钢筋混凝土板

!在四个角点增加钢性垫片，在中心点处直接施加-2mm，压碎关闭计算完毕
/clear,nostart
/config,nres,5000
/prep7
!定义单元类型
et,1,solid65 !混凝土单元
et,2,solid45
keyopt,1,6,3
keyopt,1,5,1
*afun,deg
R,1,2,0.01,90,0,2,0.01,0,0 !整体模型，实常数
!定义材料属性
mp,ex,1,24000 !混凝土材料属性
mp,prxy,1,0.2
tb,conc,1,1,9
tbdata,,0.35,1,3.1125,-1
mp,ex,2,2e5 !纵向受拉钢筋，钢支座材料属性
mp,prxy,2,0.25
tb,bkin,2,1,2,1
tbdata,,360,20000
blc4,,,1000,1000,100
blc4,,,100,100,-50

wpoff,1000,0,0
blc4,,,-100,100,-50
wpoff,0,1000,0
blc4,,,-100,-100,-50
wpoff,-1000,0,0
blc4,,,100,-100,-50
/view,1,1,1,1
/ANG,1,240,ZS,1
/replot
/REP,FAST
lsel,s,loc,z,50
lesize,all,,,2
lsel,all
lsel,s,loc,x,500
lesize,all,,,20
lsel,all
lsel,s,loc,y,500
lesize,all,,,20
vsel,s,loc,z,50
vmesh,all
type,2
mat,2
lsel,s,loc,z,-25
lesize,all,,,1
lsel,all
lsel,s,loc,z,-50
lesize,all,,,2
vsel,s,loc,z,-25
vmesh,all
nummrg,all
numcmp,all
fini
/solu
asel,s,loc,z,