|
利用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
|
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
显身卡
