16m后张法简支梁ansys计算命令流,绝对原创
fini
/clea
/filn,recaculate
/prep7
*afun,deg
local,11,,,,,,,1.15
k,1,
k,2,0,,-0.12
k,3,0.46,,-0.25
k,4,0.46,,-0.8
k,5,1.71,,-0.8
k,6,1.71
l,1,2
*repeat,5,1,1
l,6,1
l1=lsinqr(,14)
k,7,0.46+0.25,,-0.405
k,8,0.46+0.25,,-0.21
k,9,0.46+0.25+0.15,,-0.21+0.08
k,10,1.3,,-0.13
k,11,1.45,,-0.21
k,12,1.45,,-0.405
l,7,8
*repeat,5,1,1
clocal,12,,1.45,,-0.405,,,
lsel,s,,,7,lsinqr(,14)
lsymm,z,all
nummrg,kp,0.1,0.1
alls
numcmp,all
al,all
l2=lsinqr(,14)
!
csys,11
k1=kpinqr(,14)
k,k1+1,1.71,
k,k1+2,1.71,,-0.8
k,k1+3,1.71+1.26,,-0.8
k,k1+4,1.71+1.26,,
l,k1+1,k1+2
*repe,3,1,1
l,k1+1,k1+4
lgen,2,l1+1,l2,1,1.26
l3=lsinqr(,14)
lsel,s,,,l2+1,l3
al,all
agen,8,2,,,1.26
clocal,12,,13.5/2,,,
arsym,x,1
ksel,s,loc,z
nummrg,kp,1e-3,1e-3
alls
numcmp,all
!拖拉成体(上部面以处理)
a1=arinqr(,14)
voffst,1,16,
*repeat,a1,1
v1=vlinqr(,14)
vsymm,y,all
vsel,s,,,1,v1
vdele,all,,,1
alls
numcmp,all
k1=kpinqr(,14)
l1=lsinqr(,14)
!预应力筋
csys,11
cloc,12,,0.46+0.125,8,-0.8
k,k1+1,,0,0.085
k,k1+2,,5.5,0.085
k,k1+3,,6,0.086
k,k1+4,,6.5,0.096
k,k1+5,,7,0.112
k,k1+6,,7.5,0.13
k,k1+7,,7.8,0.14
l,k1+1,k1+2
*repe,6,1,1
l2=lsinqr(,14)
lsymm,y,l1+1,l2
k1=kpinqr(,14)
l3=lsinqr(,14)
k,k1+1,,0,0.195
k,k1+2,,3,0.195
k,k1+3,,3.5,0.196
k,k1+4,,4,0.206
k,k1+5,,4.5,0.229
k,k1+6,,5,0.264
k,k1+7,,5.5,0.312
k,k1+8,,6,0.373
k,k1+9,,6.5,0.446
k,k1+10,,7,0.532
k,k1+11,,7.5,0.62
k,k1+12,,7.728,0.66
l,k1+1,k1+2
*repe,11,1,1
l4=lsinqr(,14)
lsymm,y,l3+1,l4
lsel,s,,,l1+1,lsinqr(,14)
ksll
nummrg,kp,1e-3,1e-3
lgen,2,all,,,12.33
cm,ls1,line
!平弯筋
k1=kpinqr(,14)
l1=lsinqr(,14)
csys,11
cloc,12,,0.46+0.125+0.95,8,-0.8
k,k1+1,,0,0.085
k,k1+2,,5.5,0.085
k,k1+3,,6,0.086
k,k1+4,,6.5,0.096
k,k1+5,0.02,7,0.112
k,k1+6,.04,7.5,0.13
k,k1+7,0.04,7.8,0.14
l,k1+1,k1+2
*repe,6,1,1
l2=lsinqr(,14)
lsymm,y,l1+1,l2
k1=kpinqr(,14)
l3=lsinqr(,14)
k,k1+1,,0,0.195
k,k1+2,,3,0.195
k,k1+3,,3.5,0.196
k,k1+4,,4,0.206
k,k1+5,,4.5,0.229
k,k1+6,,5,0.264
k,k1+7,,5.5,0.312
k,k1+8,,6,0.373
k,k1+9,,6.5,0.446
k,k1+10,0.02,7,0.532
k,k1+11,0.04,7.5,0.62
k,k1+12,0.04,7.728,0.66
l,k1+1,k1+2
*repe,11,1,1
l4=lsinqr(,14)
lsymm,y,l3+1,l4
lsel,s,,,l1+1,lsinqr(,14)
ksll
nummrg,kp,1e-3,1e-3
cm,ls2,line
csys,11
clocal,12,,1.71
lsymm,x,ls2
cm,ls2,line
lgen,9,ls2,,,1.26
cmsel,a,ls1
cm,ls1,line
alls
numcmp,all
!单元属性以及单元类型
et,1,65
keyopt,1,1,1
keyopt,1,5,2
keyopt,1,7,1
tb,concr,
tbdata,1,0.5,1,2.07e6,-1
mp,ex,1,3.45e10
mp,prxy,1,0.2
mp,dens,1,2.6e3
mp,ex,2,2.1e11
mp,prxy,2,0.28
mp,dens,2,7850
mp,ex,3,1.95e11
mp,dens,3,7850
et,3,63
r,1,0.001
asel,s,loc,y,0
aatt,1,1,3
esize,0.08
amesh,all
r,11,2,0.005,90,,2,0.005
rmore,90,,2,0.005,91.5
r,2,0.000181,0.00357692308
et,2,link8
!划分边板
csys,11
alls
vatt,1,11,1
alls
esiz,0.2
vswe,all
!预应力筋
cmsel,s,ls1
latt,3,2,2
esize,0.2
lmesh,all
!耦合
nsll
,,1
n1=ndinqr(,13)
*dim,nd,,n1
n2=1
*do,i,1,n1
*get,nd(i),node,n2,nxth
n2=nd(i)
*enddo
alls
numcmp,all
nslv,s,1
*do,i,1,n1
nsel,a,,,nd(i)
n2=nnear(nd(i))
cp,next,all,nd(i),n2
nsel,u,,,n2
nsel,u,,,nd(i)
*enddo
sfcum,pres,add
!???
asel,s,loc,z,
nsla,s,1
nsel,r,loc,x,0,1
sf,all,pres,100000/ndinqr(,13)
asel,s,loc,z,
nsla,s,1
nsel,r,loc,x,13.5,13
sf,all,pres,80000/ndinqr(,13)
asel,s,loc,z,
nsla,s,1
sf,all,pres,2.6e4*13.5*16*0.15/ndinqr(,13)
!约束
alls
csys,11
numcmp,all
asel,s,loc,z,-0.8
nsla,s,1
nsel,r,loc,x,0.46,13.5-0.46
nsel,r,loc,y,0,0.15
d,all,uz,0,,,,,uy,ux
asel,s,loc,z,-0.8
nsla,s,1
nsel,r,loc,x,0.46,13.5-0.46
nsel,r,loc,y,16,16-0.15
d,all,uz,0,,,,
csys,11
asel,s,loc,z,
nsla,s,1
*dim,nod1,,10
*dim,dis,,10,2
n1=1
*do,i,2,6,4
*do,j,0,16,4
nod1(n1)=node(i,j,0)
n1=n1+1
*enddo
*enddo
alls
csys,0
acel,,,9.8
lswrite,1
csys,11
asel,s,loc,z,
nsla,s,1
nsel,r,loc,x,3,3.3
nsel,r,loc,y,8,8.4
sf,all,pres,1.608e5/2/ndinqr(,13)
asel,s,loc,z,
nsla,s,1
nsel,r,loc,x,3+2.3,3.3+2.3
nsel,r,loc,y,8,8.4
sf,all,pres,1.608e5/2/ndinqr(,13)
asel,s,loc,z,
nsla,s,1
nsel,r,loc,x,3,3.3
nsel,r,loc,y,8+4,8.4+4
sf,all,pres,6.12e4/2/ndinqr(,13)
asel,s,loc,z,
nsla,s,1
nsel,r,loc,x,3+2.3,3.3+2.3
nsel,r,loc,y,8+4,8.4+4
sf,all,pres,6.12e4/2/ndinqr(,13)
alls
lswrite,2
save
fini
/sol
nropt,full
solcontrol,on,1
eqslv,pcg
cnvtol,u,,0.01,0
!autots,on
!deltim,0.02,0,1
lssolve,1,2,1
save
fini
/post
alls