混凝土非线性分析应力应变曲线输出解读Word格式文档下载.docx
《混凝土非线性分析应力应变曲线输出解读Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《混凝土非线性分析应力应变曲线输出解读Word格式文档下载.docx(32页珍藏版)》请在冰豆网上搜索。
graphvariables.在对话框上"
thefirstvariable"
对应的交互框中输入2。
点击确定键,则预想的应力-应变曲线就显示在屏幕上。
ok!
试试看!
欢迎大家继续批评指正!
/prep7
et,1,solid65
et,2,link8
et,3,solid45
r,1
r,2,236e-6
r,3,75.4e-6
r,4,151.2e-6
!
混凝土本构关系
mp,ex,1,2.522e10
mp,prxy,1,0.2
mp,dens,1,2500
tb,concr,1,1
tbdata,,0.5,0.95,5e6,-1
tb,miso,1,1,7
tbtemp,0
tbpt,,0.00068,17.15e6
tbpt,,0.00136,28.03e6
tbpt,,0.00204,33.62e6
tbpt,,0.00272,35.96e6
tbpt,,0.0034,36.53e6
tbpt,,0.00408,36.53e6
纵筋的本构关系
mp,ex,2,2.0e11
mp,prxy,2,0.28
tb,bkin,2
tbdata,,362e6,2e9
!
箍筋的本构关系
mp,EX,3,2.1E11
mp,PRXY,3,.27
tb,bkin,3
tbtemp,0
tbdata,,272.5e6,2.1e9
建立模型
local,11,1,0,0,0,0,-90,0!
建立局部坐标系,采用柱坐标,绕y轴顺时针转90度
csys,11
*do,j,1,8
n,j,0.185,j*360/8,-0.6
*enddo
ngen,3,100,1,8,1,0,0,0.3!
复制节点,3次,编号增量,源节点1~8,0,0,0.3
ngen,4,100,201,208,1,0,0,0.15
ngen,7,100,501,508,1,0,0,0.3
type,2!
指定单元属性,建立基础箍筋单元
real,3
mat,3
*do,i,1,6
*do,j,1,7
e,(i-1)*100+j,(i-1)*100+j+1
*enddo
e,(i-1)*100+1,(i-1)*100+8
建立柱上箍筋单元
real,4
*do,i,1,5
e,(i+5)*100+j,(i+5)*100+j+1
e,(i+5)*100+1,(i+5)*100+8
建立纵向钢筋
real,2
mat,2
*do,i,1,11
e,(i-1)*100+j,i*100+j
allsel,all!
合并节点,压缩节点
nummrg,all
numcmp,all
wpcsys,1,11
cylind,0,0.185,-0.6,0,0,360!
建立圆柱,内外径i-or,高度z1-z2,起终角度
cylind,0,0.185,0,0.45,0,360
cylind,0.185,0.3,-0.6,0,0,360
cylind,0,0.185,0.45,2.25,0,360
cylind,0,0.185,2.25,2.55,0,360
cylind,0.185,0.20,2.25,2.55,0,360
cylind,0,0.185,2.55,2.565,0,360
wprota,,90!
绕y轴逆时针转90度
vsel,s,,,1,7,1!
选择体,选择类型,item,comb,最小编号,最大编号,增量
vsbw,all!
体相减
numcmp,all
wprota,,,90!
绕z轴逆时针转90度
vsel,s,,,1,14,1
vsbw,all
vsel,s,,,1,28,1
vglue,all
wprota,,,-90
wprota,,-90
allsel,all
vsel,s,,,1,16,1
vsel,a,,,21,24,1
vatt,1,1,1
vsel,s,,,17,20,1
vsel,a,,,25,28,1
vatt,2,1,3
allsel,all
lsel,s,,,1,8,1!
选择线
lsel,a,,,11,22,1!
补选(a)线
lsel,a,,,25,40,1
lsel,a,,,43,46,1
lesize,all,,,2!
线尺寸为2
lsel,s,,,47,50,1
lsel,a,,,53,56,1
lsel,a,,,58,73,1
lesize,all,,,2
lsel,s,,,98,107,1
lsel,a,,,110,112,2
lesize,all,,,1
lsel,s,,,23,24,1
lsel,a,,,9,10,1
lsel,s,,,78,82,1
lesize,all,,,3
lsel,s,,,83,87,1
lesize,all,,,6
lsel,s,,,41,42,1
lsel,a,,,88,92,1
lsel,s,,,93,97,1
mshkey,1
mshape,0,3d
vsel,s,,,1,28
vmesh,all
wprota,,90
csys,0
/solu!
进入加载、求解
nsel,s,loc,x,0.3!
选择半径在x=3处的节点,并施加约束
d,all,all
OUTRES,ALL,2!
输出获得的解
pred,on!
打开预应力选项
cnvtol,f,,0.05,2!
定义收敛条件,使用缺省的VALUE
autots,1!
打开自动时间步控制
lnsrch,1!
打开线性搜索
ncnv,2!
如果不收敛时结束而不退出
neqit,50!
每一子步中方程的迭代次数限值
save
time,1!
定义第1载荷步
nsubst,10
asel,s,,,29,32,1
sfa,all,,pres,5468699!
在已选择的面上施加压力荷载
acel,,9.8!
重力加速度y=9.8
kbc,1!
阶跃加载方式
lswrite,1
time,10!
定义第2荷载步
cp,1,ux,298,300,303!
将节点298,300,303ux向位移耦合
d,298,ux,0.015!
再借点298处施加ux位移0.015
KBC,0
NSUBST,20,1000000,10
lswrite,2
time,20!
定义第3载荷步
d,298,ux,-0.015
NSUBST,100,1000000,10
lswrite,3
time,30!
定义第4载荷步
d,298,ux,0.015
lswrite,4
time,40!
定义第5载荷步
d,298,ux,-0.015
lswrite,5
time,50!
定义第6载荷步
lswrite,6
time,60!
定义第7载荷步
lswrite,7
time,70!
定义第8载荷步
lswrite,8
time,80
lswrite,9
time,90
lswrite,10
time,100
lswrite,11
time,110
d,298,ux,0.03
NSUBST,200,1000000,10
lswrite,12
time,120
d,298,ux,-0.03
lswrite,13
time,130
lswrite,14
time,140
lswrite,15
time,150
lswrite,16
time,160
lswrite,17
time,170
lswrite,18
time,180
lswrite,19
time,190
lswrite,20
time,200
lswrite,21
time,210
d,298,ux,0.045
NSUBST,300,1000000,10
lswrite,22
lssolve,1,22,1
命令流:
1、单调水平荷载作用下的命令流:
建模:
monobrickwall
solid65-整体式建模
*
Finish
/clear
/units,si
/com,Structural
/config,nres,400000
/prep7
/title,mono-brick-wall
depth_sup=0.24!
弹性支垫240mm
单元属性
forbrick,mortarandlevelingcoat
et,1,solid65
KEYOPT,1,1,0
KEYOPT,1,5,0
KEYOPT,1,6,0
KEYOPT,1,7,1
forelasticsupport
et,2,solid45
*--------------------------------砌体材料--------------------------------------
MU10,M5砌体材料属性
mp,prxy,1,0.15
mp,dens,1,1700
mp,ex,1,0.24e10
砌体屈服准则
TB,MKIN,1
TBTEMP,,strain
TBdata,1,0.2E-3,1.0E-3,2.0E-3,3.0E-3,4.8E-3
TBTEMP,,,
TBDATA,1,0.48E6,0.82E6,1.32E6,1.5E6,1.32E6
TBDATA,1,0.48E6,1.31E6,2.11E6,2.4E6,2.11E6
/XRANGE,0,0.02
TBPLOT,MKIN,1
砌体破坏准则
hntrl=0.13E6
hntrl=0.21E6
tb,concr,1,,4,
tbdata,,0.2,0.9,hntrl,-1!
不考虑压碎
*--------------------------------弹性支撑--------------------------------------
弹性支撑暂取钢筋材料属性
mp,dens,2,7.600
mp,prxy,2,0.30
realconstants
*--------------------------------modeling--------------------------------------
geometrymodel
block,0,3.12,0,0.24,0,3.12
block,-0.24,3.12,0,0.24,3.12,3.36
Vsel,s,volu,,1,2,1$aslv,s$lsla,s
lesize,all,0.24
Vsel,s,volu,,1$vatt,1,1,1,!
brickwall
Vsel,s,volu,,2$vatt,2,1,2,!
vsel,all!
必须再次选择所有单元
VSWEEP,ALL
allsel
/device,vector,1
/eshape,1
*-----------------------------------耦合---------------------------------------
elasticsupportnodesset
esel,s,type,,2$nsle,s,all$nsel,r,loc,z,3.12
nsel,u,loc,x,-0.24$cm,slavenode,node
othernodeset
esel,s,type,,1$nsle,s,all$cm,masternode,node
/input,ucouple,mac
finish
耦合文件ucouple,mac:
*-----------------------------------开始---------------------------------------
allsel!
最好保留这句命令
*******将从属节点编号依次存入数组****************
cmsel,s,slavenode
*get,count1_node,node,0,count
*del,slave_node
*dim,slave_node,array,count1_node
*get,slave_node
(1),node,0,num,min
*do,i,2,count1_node
slave_node(i)=ndnext(slave_node(i-1))
*******将主节点编号依次存入数组****************
cmsel,s,masternode
*get,count2_node,node,0,count
*del,master_node
*dim,master_node,array,count2_node
*get,master_node
(1),node,0,num,min
*do,i,2,count2_node
master_node(i)=ndnext(master_node(i-1))
********将与从属节点耦合的节点数组初始化****************
*del,cp_node
*dim,cp_node,array,count1_node
*do,i,1,count1_node
cp_node(i)=0
*********开始选择程序****************
*do,i,1,count1_node
kk=1
k=1
*dowhile,kk
k=nnear(slave_node(i))
nsel,s,cp,,all
kk=nsel(k)+0.001
nsel,u,node,,k
cm,masternode,node
cp_node(i)=k
*******选择完毕****************
*******开始耦合****************
cp,next,UX,slave_node(i),cp_node(i)
cp,next,UY,slave_node(i),cp_node(i)
cp,next,UZ,slave_node(i),cp_node(i)
*******耦合完毕****************
求解方法NR法:
brickwallNRsolution
displacementload
用NR法
*---------------------------------solution-------------------------------------
/solu
antype,0
-----------constraint----------------
nsel,s,loc,z,0
esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24
d,all,uy
d,all,uz
lumpm,on
solcontrol,on
pred,on
lnsrch,on
打开自动时间步控制
如果不收敛时结束而不退出
compression=0.24E7
*------------------------------------loadsteps--------------------------------
*------------------------------------1st---------------------------------------
acel,,,9.8
nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.119
esln,s,0,all
sfe,all,6,pres