混凝土非线性分析应力应变曲线输出解读Word格式文档下载.docx

混凝土非线性分析应力应变曲线输出解读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