钢筋拉拔实验模拟

建模命令：
FINISH
/CLEAR,NOSTART
/TITLE,Analysis of la ba test
/PREP7
!申明，单位m,n,Pa,kg!
!定义混凝土的特性
! 单元属性
ET,1,SOLID65
KEYOPT,1,1,1
KEYOPT,1,3,0
KEYOPT,1,5,0
KEYOPT,1,6,3
KEYOPT,1,7,1
KEYOPT,1,8,1
!实参数
! 不考虑纵向,只有横向钢筋
R,1,
!材料属性
!混凝土基本材料属性
MP,EX,1,1.51E+10 !c30 3.0e10
mp,prxy,1,0.1667
!屈服准则
TB,MISO,1,1,19,
TBPT,,0,0
TBPT,,0.0002,3.02E+06
TBPT,,0.0004,5.93E+06
TBPT,,0.0006,8.63E+06
TBPT,,0.0008,1.10E+07
TBPT,,0.0010,1.31E+07
TBPT,,0.0012,1.48E+07
TBPT,,0.0014,1.60E+07
TBPT,,0.0016,1.69E+07
TBPT,,0.0018,1.74E+07
TBPT,,0.0020,1.75E+07
TBPT,,0.0022,1.74E+07
TBPT,,0.0024,1.70E+07
TBPT,,0.0026,1.65E+07
TBPT,,0.0028,1.58E+07
TBPT,,0.0030,1.50E+07
TBPT,,0.0032,1.41E+07
TBPT,,0.0034,13214557.14
TBPT,,0.0036,12290307.14
!!破坏准则
TB,CONC,1,1,9,
TBTEMP,0
TBDATA,1,0.7,1.0,2.0E6,-1,,
! 定义粘结单元类型
!*
ET,2,COMBIN39

KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,3
KEYOPT,2,4,0
KEYOPT,2,6,0
!*
!定义弹簧的位移和力的关系
R,2,0,0,0.0001,10,0.0002,20,
RMORE,0.0003,30.5,0.0004,40,0.0005,40.2,
RMORE,,,,,,,
RMORE,,,,,,,
RMORE,,,,,,,
RMORE,,,,,,,
RMORE,,,,,
!混凝土建模
!BLC4,-0.075,-0.075,0.15,0.15
CYL4,0,0,0.02
!asba,1,2
allsel,all
!lesize,all,,,4
k,100,
k,101,,,0.15
l,100,101
lesize,5,,,4
CYL4,0,0,0.02, ,0.05
lesize,5,,,4
vdrag,1,2,,,,,5
FLST,2,1,8
FITEM,2,0,0,0.2
WPAVE,P51X
wprota,,90
*do,i,1,4,1
wprota,,,45
vsbw,all
*enddo
lsel,s,loc,z
lesize,all,,,2
vsweep,all
WPCSYS,-1,0
allsel,all
!!!!!!!!!!!!!!!!!!!!!!!!!!
!allsel,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!弹簧
real,2
mat,2
type,2
cpnum=0
cm,n-zhong,node !需要耦合的节点
!建立一个耦合小程序
*GET,n_num,NODE,,COUNT, , , , !节点总数
*do,i,1,n_num !
cmsel,s,n-zhong
*GET,n_lowest,NODE,,NUM,MIN, , , , !号码最小的节点
*GET,n_x,NODE,n_lowest,LOC,X !该节点坐标
*GET,n_y,NODE,n_lowest,LOC,Y
*GET,n_z,NODE,n_lowest,LOC,Z

NSEL,s,LOC,X,n_x-0.001,n_x+0.001 !寻找与该节点位置相同的节点
NSEL,R,LOC,Y,n_y-0.001,n_y+0.001
NSEL,R,LOC,z,n_z-0.001,n_z+0.001
*GET,n_number,NODE,,NUM,Max, , , , !得到该点的编号

cm,n_cp_cp,node !位置相同的节点形成一个组
cmsel,s,n-zhong
cmsel,u,n_cp_cp
cm,n-zhong,node !取消这些点后剩余的点形成组
*GET,n_num_1,NODE,,COUNT, , , , !节点总数
*if,n_num_1,lt,2,exit !如果节点数小于二则退出
cmsel,s,n_cp_cp
*GET,n_num,NODE,,COUNT, , , ,
*if,n_num,gt,1,then
e,n_lowest,n_number
!CP,cpnum+1,uz,all
CP,cpnum+2,uX,all
CP,cpnum+3,uy,all

!CP,cpnum+1,ROTz,all
!CP,cpnum+2,ROTX,all
!CP,cpnum+3,ROTy,all
cpnum=cpnum+3
*else

*endif
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
计算：


/solu
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
FLST,2,32,1,ORDE,32
FITEM,2,91
FITEM,2,93
FITEM,2,113
FITEM,2,118
FITEM,2,123
FITEM,2,128
FITEM,2,133
FITEM,2,141
FITEM,2,177
FITEM,2,-178
FITEM,2,183
FITEM,2,189
FITEM,2,-190
FITEM,2,195
FITEM,2,203
FITEM,2,-204
FITEM,2,210
FITEM,2,-211
FITEM,2,217
FITEM,2,-218
FITEM,2,224
FITEM,2,-225
FITEM,2,231
FITEM,2,233
FITEM,2,282
FITEM,2,284
FITEM,2,292
FITEM,2,294
FITEM,2,308
FITEM,2,317
FITEM,2,326
FITEM,2,335
!*
/GO
D,P51X, , , , , ,ALL, , , , ,
FLST,2,17,1,ORDE,13
FITEM,2,1
FITEM,2,-6
FITEM,2,13
FITEM,2,-14
FITEM,2,252
FITEM,2,-253
FITEM,2,262
FITEM,2,268
FITEM,2,274
FITEM,2,300
FITEM,2,-301
FITEM,2,303
FITEM,2,-304
!*
/GO
F,P51X,FZ,0.1
!SAVE
allsel,all
!solve