广场 › [ANSYS] › 预应力混凝土连续梁的ANSYS模拟分析

szg1999 发表于 2011-8-23 16:13:10

预应力混凝土连续梁的ANSYS模拟分析

利用Solid65单元与Link8分析预应力混凝土结构

/prep7

ET,1,solid65
et,2,link8
et,3,solid45
keyopt,1,7,1

r,1,4,0.01744,0,0      
r,2,2,0.056,0,0         
r,3,4,0.01,90,0         
mp,ex,1,3.25e10
mp,nuxy,1,0.2
mp,dens,1,2500
tb,concr,1,1
tbtemp,0
tbdata,,0.5,0.95,3.5e6,-1
tb,miso,1,1,12         !
tbpt,,0.00040985,13.32e6
tbpt,,0.0006768,21.312e6
tbpt,,0.0010152,27.972e6
tbpt,,0.0013536,32e6
tbpt,,0.0015228,33e6
tbpt,,0.001692,33.3e6
tbpt,,0.0018612,32.86e6
tbpt,,0.0021996,30.2e6
tbpt,,0.002538,26.7e6
tbpt,,0.0028764,23.32e6
tbpt,,0.0032148,20.4e6
tbpt,,0.0037224,16.9e6


mp,EX,2,2.0e11
mp,PRXY,2,.27
tb,bkin,2
tbtemp,0
tbdata,,407e6,0


mp,ex,3,1.95e11
mp,prxy,3,0.27
mp,alpx,3,1e-5
tb,bkin,3
tbtemp,0
tbdata,,1465e6,6.633e10
r,4,139.98e-6
r,5
r,6,5,0.0251,0,0         


mp,EX,4,2.1e11
mp,PRXY,4,.27
tb,bkin,4
tbtemp,0
tbdata,,275e6,0



mp,EX,5,2.0e11
mp,PRXY,5,.27
tb,bkin,5
tbtemp,0
tbdata,,366e6,0









/prep7
local,11,1,1139.971e-3,7296.1857e-3
csys,11
*dim,a,,6
*do,i,1,6
a(i)=260.9365+7.9553/5*(i-1)
*enddo
/prep7
*do,i,1,6
k,i,7236.5395e-3,a(i)
*enddo
*do,i,1,5
l,i,i+1
*enddo

local,12,1,1104.1693e-3,5269.1667e-3
csys,12
*dim,b,,2
*do,i,1,2
b(i)=268.8542+2.75/2*i
*enddo
/prep7
*do,i,1,2
k,i+6,5209.2083e-3,b(i)
*enddo
*do,i,1,2
l,i+5,i+6
*enddo

local,13,1,1124.814e-3,4446.1271e-3
csys,13
*dim,c,,5
*do,i,1,5
c(i)=271.6356+13.2296/5*i
*enddo
/prep7
*do,i,1,5
k,i+8,4385.9141e-3,c(i)
*enddo
*do,i,1,5
l,i+7,i+8
*enddo

local,14,1,2500e-3,-723.4697e-3
csys,14
*dim,d,,5
*do,i,1,5
d(i)=105.0391-30.0783/5*i
*enddo
/prep7
*do,i,1,5
k,i+13,963.4697e-3,d(i)
*enddo
*do,i,1,5
l,i+12,i+13
*enddo

local,15,1,3875.1859e-3,4446.1271e-3
csys,15
*dim,e,,5
*do,i,1,5
e(i)=255.1348+13.2296/5*i
*enddo
/prep7
*do,i,1,5
k,i+18,4385.9141e-3,e(i)
*enddo
*do,i,1,5
l,i+17,i+18
*enddo

local,16,1,3895.8307e-3,5269.1667e-3
csys,16
*dim,f,,2
*do,i,1,2
f(i)=268.3958+2.75/2*i
*enddo
/prep7
*do,i,1,2
k,i+23,5209.2083e-3,f(i)
*enddo
*do,i,1,2
l,i+22,i+23
*enddo

local,17,1,3860.029e-3,7296.1857e-3
csys,17
*dim,g,,5
*do,i,1,5
g(i)=271.1083+7.9553/5*i
*enddo
/prep7
*do,i,1,5
k,i+25,7236.5395e-3,g(i)
*enddo
*do,i,1,5
l,i+24,i+25
*enddo


lsel,s,,,1,29,1
latt,3,4,2
bfl,all,temp,-142000*0.8/(139.98*1.95)   


allsel,all
lsel,s,,,1,29,1,
lesize,all,,,1
lmesh,all


/prep7
csys,0
k,101,0,0,0.075
k,102,0,0,-0.075
k,103,0,0.06,-0.075
k,104,0,0.06,0.075
k,105,0,0.24,0.075
k,106,0,0.24,-0.075
k,107,0,0.3,-0.075
k,108,0,0.3,0.075


kgen,2,101,108,1,1,0,0,100
kgen,2,201,208,1,0.25,0,0,100
kgen,2,301,308,1,1,0,0,100
kgen,2,401,408,1,0.5,0,0,100
kgen,2,501,508,1,1,0,0,100
kgen,2,601,608,1,0.25,0,0,100
kgen,2,701,708,1,1,0,0,100
a,101,102,103,104
a,103,104,105,106
a,105,106,107,108
a,201,202,203,204
a,203,204,205,206
a,205,206,207,208
a,301,302,303,304
a,303,304,305,306
a,305,306,307,308
a,401,402,403,404
a,403,404,405,406
a,405,406,407,408
a,501,502,503,504
a,503,504,505,506
a,505,506,507,508
a,601,602,603,604
a,603,604,605,606
a,605,606,607,608
a,701,702,703,704
a,703,704,705,706
a,705,706,707,708
a,801,802,803,804
a,803,804,805,806
a,805,806,807,808
vext,1,3,1,-0.02,0,0
vext,1,3,1,1,0,0
vext,4,6,1,0.25,0,0
vext,7,9,1,1,0,0
vext,10,12,1,0.5,0,0
vext,13,15,1,1,0,0
vext,16,18,1,0.25,0,0
vext,19,21,1,1,0,0
vext,22,24,1,0.02,0,0




vsel,s,,,6,12,3
vatt,1,1,1
vsel,s,,,18,24,3
vatt,1,1,1



vsel,s,,,15
vatt,1,2,1

vsel,s,,,4,10,3
vatt,1,2,1

vsel,s,,,16,22,3
vatt,1,2,1







vsel,s,,,13
vatt,1,6,1

vsel,s,,,5,23,3
vatt,1,3,1




vsel,s,,,1,3,1
vsel,a,,,25,27,1
vatt,2,5,3





allsel,all
vsel,s,,,1,27,1
nsel,u,,,1,30,1


nummrg,all
numcmp,all









lsel,s,,,30,100,10
lsel,a,,,32,102,10
lsel,a,,,35,105,10
lsel,a,,,38,108,10
lesize,all,,,2
allsel,all
lsel,s,,,31,101,10
lsel,a,,,33,103,10
lsel,a,,,37,107,10
lsel,a,,,39,109,10
lesize,all,,,1
allsel,all
lsel,s,,,34,104,10
lsel,a,,,36,106,10
lesize,all,,,6
allsel,all
lsel,s,,,127,201,74
lesize,all,,,1
allsel,all
lsel,s,,,134,182,48
lesize,all,,,5
allsel,all
lsel,s,,,142,174,32
lesize,all,,,2
allsel,all
lsel,s,,,150,166,16
lesize,all,,,5
allsel,all
lsel,s,,,158
lesize,all,,,5

allsel,all
mshape,0,3d
mshkey,1
vmesh,all


k,1001,0,-0.02,0.075
k,1002,0,-0.02,-0.075
k,1003,-0.02,-0.02,-0.075
k,1004,-0.02,-0.02,0.075
a,1001,1002,1003,1004
k,1005,2.45,-0.02,0.075
k,1006,2.45,-0.02,-0.075
k,1007,2.55,-0.02,-0.075
k,1008,2.55,-0.02,0.075
a,1005,1006,1007,1008
k,1009,5,-0.02,0.075
k,1010,5,-0.02,-0.075
k,1011,5.02,-0.02,-0.075
k,1012,5.02,-0.02,0.075
a,1009,1010,1011,1012
vext,121,123,1,0,0.02,0
vsel,s,,,28,30,1
vatt,2,5,3
lsel,s,,,202,212,2
lesize,all,,,2
allsel,all
lsel,s,,,203,205,2
lesize,all,,,1
allsel,all
lsel,s,,,207,209,2
lesize,all,,,2
allsel,all
lsel,s,,,211,213,2
lesize,all,,,1
allsel,all
lsel,s,,,218,229,11
lsel,a,,,237
lesize,all,,,1
allsel,all
vsel,s,,,28,30,1
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,1,30,1
nsel,u,,,1,30,1

nummrg,all
numcmp,all


allsel,all   
csys,11
lsel,s,,,1,5
nsll,s,1   
cm,cm1jnod,node   
*get,max1,node,0,count
*dim,ojd,,max1
*dim,jd,,max1   
*get,nod1,node,0,num,min
ojd(1)=nod1
*do,i,2,max1
ojd(i)=ndnext(ojd(i-1))
*enddo
allsel,all
nsel,all
cmsel,u,cm1jnod   
*do,i,1,max1
nod1=ojd(i)
j=nnear(nod1)
jd(i)=j
*enddo
nsel,all
ji=1
*do,i,1,max1
csys,11
nrotat,1,6,1
nrotat,58,90,32
nrotat,114
nrotat,186,187,1
nrotat,189
cp,ji,ux,ojd(i),jd(i)   
cp,ji+1,uz,ojd(i),jd(i)   
ji=ji+2
*enddo
allsel,all


allsel,all
csys,12
lsel,s,,,7
nsll,s,1   
cm,cm1jnod2,node   
*get,max2,node,0,count
*dim,ojd2,,max2
*dim,jd2,,max2   
*get,nod2,node,0,num,min
ojd2(1)=nod2
*do,i,2,max2
ojd2(i)=ndnext(ojd2(i-1))
*enddo
allsel,all
nsel,all
cmsel,u,cm1jnod2   
*do,i,1,max2
nod2=ojd2(i)
j=nnear(nod2)
jd2(i)=j
*enddo
nsel,all
ji=13
*do,i,1,max2
csys,12
nrotat,7,8,1
nrotat,225,231,6
cp,ji,ux,ojd2(i),jd2(i)      
cp,ji+1,uz,ojd2(i),jd2(i)   
ji=ji+2
*enddo
allsel,all


allsel,all   
csys,13
lsel,s,,,9,12,1
nsll,s,1   
cm,cm1jnod3,node   
*get,max3,node,0,count
*dim,ojd3,,max3
*dim,jd3,,max3   
*get,nod3,node,0,num,min
ojd3(1)=nod3
*do,i,2,max3
ojd3(i)=ndnext(ojd3(i-1))
*enddo
allsel,all
nsel,all
cmsel,u,cm1jnod3   
*do,i,1,max3
nod3=ojd3(i)
j=nnear(nod3)
jd3(i)=j
*enddo
nsel,all
ji=17
*do,i,1,max3
csys,13
nrotat,9,13,1
nrotat,300,321,21
nrotat,375,385,5
cp,ji,ux,ojd3(i),jd3(i)      
cp,ji+1,uz,ojd3(i),jd3(i)   
ji=ji+2
*enddo
allsel,all


allsel,all
csys,14
lsel,s,,,14,17,1
nsll,s,1   
cm,cm1jnod4,node   
*get,max4,node,0,count
*dim,ojd4,,max4
*dim,jd4,,max4   
*get,nod4,node,0,num,min
ojd4(1)=nod4
*do,i,2,max4
ojd4(i)=ndnext(ojd4(i-1))
*enddo
allsel,all
nsel,all
cmsel,u,cm1jnod4   
*do,i,1,max4
nod4=ojd4(i)
j=nnear(nod4)
jd4(i)=j
*enddo
nsel,all
ji=27
*do,i,1,max4
csys,14
nrotat,14,18,1
nrotat,485,488,1
nrotat,456
cp,ji,ux,ojd4(i),jd4(i)   
cp,ji+1,uz,ojd4(i),jd4(i)   
ji=ji+2
*enddo
allsel,all


allsel,all   
csys,15
lsel,s,,,19,22,1
nsll,s,1   
cm,cm1jnod5,node   
*get,max5,node,0,count
*dim,ojd5,,max5
*dim,jd5,,max5   
*get,nod5,node,0,num,min
ojd5(1)=nod5
*do,i,2,max5
ojd5(i)=ndnext(ojd5(i-1))
*enddo
allsel,all
nsel,all
cmsel,u,cm1jnod5   
*do,i,1,max5
nod5=ojd5(i)
j=nnear(nod5)
jd5(i)=j
*enddo
nsel,all
ji=37
*do,i,1,max5
csys,15
nrotat,19,23,1
nrotat,646,652,3
nrotat,549
nrotat,573
cp,ji,ux,ojd5(i),jd5(i)      
cp,ji+1,uz,ojd5(i),jd5(i)   
ji=ji+2
*enddo
allsel,all


allsel,all
csys,16
lsel,s,,,24
nsll,s,1   
cm,cm1jnod6,node   
*get,max6,node,0,count
*dim,ojd6,,max6
*dim,jd6,,max6   
*get,nod6,node,0,num,min
ojd6(1)=nod6
*do,i,2,max6
ojd6(i)=ndnext(ojd6(i-1))
*enddo
allsel,all
nsel,all
cmsel,u,cm1jnod6   
*do,i,1,max6
nod6=ojd6(i)
j=nnear(nod6)
jd6(i)=j
*enddo
nsel,all
ji=47
*do,i,1,max6
csys,16
nrotat,24,25,1
nrotat,684,690,6
cp,ji,ux,ojd6(i),jd6(i)   
cp,ji+1,uz,ojd6(i),jd6(i)   
ji=ji+2
*enddo
allsel,all


allsel,all
csys,17
lsel,s,,,26,29,1   
nsll,s,1   
cm,cm1jnod7,node   
*get,max7,node,0,count
*dim,ojd7,,max7
*dim,jd7,,max7   
*get,nod7,node,0,num,min
ojd7(1)=nod7
*do,i,2,max7
ojd7(i)=ndnext(ojd7(i-1))
*enddo
allsel,all
nsel,all
cmsel,u,cm1jnod7   
*do,i,1,max7
nod7=ojd7(i)
j=nnear(nod7)
jd7(i)=j
*enddo
nsel,all
ji=51
*do,i,1,max7
csys,17
nrotat,26,30,1
nrotat,834,835,1
nrotat,840
nrotat,759
nrotat,778
cp,ji,ux,ojd7(i),jd7(i)   
cp,ji+1,uz,ojd7(i),jd7(i)   
ji=ji+2
*enddo
allsel,all


cp,61,uy,30,778
csys,11
cp,62,uy,1,58
allsel,all
csys,0

/solu

nsel,s,loc,y,-0.02
nsel,r,loc,x,0
d,all,all
nsel,s,loc,y,-0.02
nsel,r,loc,x,2.5
d,all,uy
d,all,uz

nsel,s,,,913,915,1
d,all,uy
d,all,uz
allsel,all   


pred,on
cnvtol,f,,0.05,2
autots,1
lnsrch,1
ncnv,2
neqit,50
save


time,1e-6
outres,all,all
nsubst,1

acel,,9.8
lswrite,1


time,50
outres,all,all
nsubst,100

cp,63,uy,271,272,273
cp,64,uy,730,731,732
d,271,uy,-0.015
d,730,uy,-0.015
allsel,all
KBC,0
allsel,all
lswrite,2
allsel,all
lswrite,2




lssolve,1,2,1

hewenlong2008 发表于 2011-8-23 16:42:43

建议楼主命令后面加注释      

syky86 发表于 2011-10-6 12:12:49

确实是好东西，呵呵，谢谢啦

zhuece_tumu 发表于 2011-10-8 18:17:21

有些看不懂，实力有限啊

查看完整版本: 预应力混凝土连续梁的ANSYS模拟分析