杆梁壳体单元连接
fini$/clear$/prep7
l1=0.8$l2=2$l3=1
h1=3$h2=1.2$b=1.5
et,1,solid95$et,2,link8
et,3,beam4
mp,ex,1,3.3e10
mp,prxy,1,0.2
mp,ex,2,2.1e11
mp,prxy,2,0.3
*afun,rad
r,1,0.25*acos(-1)*(0.16*0.16-0.144*0.144)
r,2,0.00908,2.14e-3,1.45e-2,0.3,0.2
rmore,,6.52e-5
r,3,0.00908,1.45e-2,2.14e-3,0.2,0.3
rmore,,6.52e-5
/view,,1,1,1
/ang,,-120,zs,1
blc4,,,l1,b,h1
wpoff,,,h1-h2
vsbw,all
wpcsys,-1
k,20,l1+l2,,h1
k,21,l1+l2+l3,,h1
k,22,l1+l2+2*l3,,h1
k,23,l1+l2+3*l3,,h1
kgen,2,20,23,1,,b,,4
k,28,l1+l2+3*l3,b/2,h1
k,29,l1+l2+3*l3,b/2,0
l,6,20$l,12,20
*do,i,20,22
l,i,i+1
*enddo
l,7,24
l,11,24
*do,i,24,26
l,i,i+1
*enddo
l,6,24
*do,i,20,22,1
l,i,i+4
*enddo
l,23,28
l,28,27
l,28,29
dk,29,all
da,1,all
lplot
ksel,s,,,20,23
fk,all,fz,-3e4
ksel,s,,,24,27
fk,all,fz,-2e4
fk,23,fx,10e4
fk,27,my,1e3
esize,0.2
mshape,0
mshkey,1
vatt,1,,1
vmesh,all
lsel,s,loc,x,l1+0.1,l1+l2
lsel,a,loc,x,l1+l2+l3
lsel,a,loc,x,l1+l2+2*l3
cm,lk,line
latt,2,1,2
lesize,all,,,1
lmesh,all
lsel,s,loc,x,l1+0.1,l1+l2+3*l3
cmsel,u,lk$lsel,u,loc,z,0,h1/2+0.1
latt,2,2,3
lesize,all,,,4
lmesh,all
lsel,s,loc,x,l1+l2+3*l3
lsel,r,loc,z,0,h1/2+0.1
latt,2,3,3
lesize,all,,,4
lmesh,all
allsel,all
fini
/solu
allsel,all
solve
fini
/post1
pldisp,1
fini$/clear$/prep7
l1=2$l2=4$b=1.6$h=2
et,1,solid95$et,2,beam4
mp,ex,1,3.3e10
mp,prxy,1,0.2
mp,ex,2,2.1e11
mp,prxy,2,0.3
r,1,0.00908,2.14e-3,1.45e-2,0.3,0.2
rmore,,6.52e-5
/view,,1,1,1
/ang,,-120,zs,1
blc4,,,l1,b,h$wpoff,,,h/2
vsbw,all
wpoff,,b/2$wprota,,90
vsbw,all
wpcsys,-1
k,50,l1+l2/2,b/2,h/2
k,51,l1+l2,b/2,h/2
l,16,50
l,50,51
asel,s,loc,x,0
da,all,all
dk,51,all
fk,50,fy,2e4
fk,50,fx,1e4
fk,50,mx,1.5e4
vsel,all
vatt,1,,1
esize,0.2
mshpe,0$mshkey,1
vmesh,all
lsel,s,loc,x,l1+0.1,l1+l2
latt,2,1,2
lesize,all,,,4
lmesh,all
allsel,all
fini
/solu
allsel
solve
fini
/post1
pldisp,1
etable,fxi,smisc,1$etable,fxj,smisc,7
etable,fyi,smisc,2$etable,fyj,smisc,8
etable,fzi,smisc,3$etable,fzj,smisc,9
etable,mxi,smisc,4$etable,mxj,smisc,10
etable,myi,smisc,5$etable,myj,smisc,11
etable,mzi,smisc,6$etable,mzj,smisc,12
plls,myi,myj$plls,mzi,mzj
fini$/clear$/prep7
l1=0.5$l2=2$h1=0.6$h2=0.02$t=0.016
q0=3000$q=q0*t
a0=t*h2$i0=t*h2*h2*h2/12
et,1,beam3$et,2,plane42,,,3
mp,ex,1,2.1e11$mp,prxy,1,0.3
r,1,t$r,2,a0,i0,h2
blc4,,,l1,h1
wpoff,,h1/2
wprota,,90
asbw,all
wpcsys,-1
k,50,l1+l2,h1/2
ksel,s,loc,x,l1
ksel,r,loc,y,h1/2
*get,kp0,kp,,num,max
ksel,all
l,kp0,50
lsel,s,loc,y,h1/2
lsel,r,loc,x,l1,l1+l2
cm,lk,line
latt,1,2,1
lesize,all,,,10
lmesh,all
asel,all
aatt,1,1,2
esize,0.1
mshpe,0
mshkey,1
amesh,all
/eshape,1
lsel,s,loc,x,0
dl,all,,all
cmsel,s,lk
esll,s
sfbeam,all,1,pres,q
ksel,s,,,kp0
nslk,s
*get,n1,node,,num,max
nsel,s,loc,x,l1$nsel,r,loc,y,h1/2,h1
n2=nnear(n1)
nsel,s,loc,x,l1$nsel,r,loc,y,0,h1/2
n3=nnear(n1)
ce,1,0,n2,ux,1,n3,ux,-1,n1,rotz,ny(n2)-ny(n3)
eplot
allsel,all
fini
/solu
solve
/post1
pldisp,1
fini$/clear$/prep7
l1=0.5$l2=2$h1=0.6$h2=0.02$t=0.016
q0=3000$q=q0*t
a0=t*h2$i0=t*h2*h2*h2/12
et,1,beam3$et,2,plane42,,,3
mp,ex,1,2.1e11$mp,prxy,1,0.3
r,1,t$r,2,a0,i0,h2$r,3,a0,1e4*i0,h2
blc4,,,l1,h1
wpoff,,h1/2
wprota,,90
asbw,all
wpcsys,-1
k,50,l1+l2,h1/2
ksel,s,loc,x,l1
ksel,r,loc,y,h1/2
*get,kp0,kp,,num,max
ksel,all
l,kp0,50
lsel,s,loc,y,h1/2
lsel,r,loc,x,l1,l1+l2
cm,lk,line
latt,1,2,1
lesize,all,,,10
lmesh,all
asel,all
aatt,1,1,2
esize,0.1
mshpe,0
mshkey,1
amesh,all
/eshape,1
type,1
real,3
e,24,1
lsel,s,loc,x,0
dl,all,,all
cmsel,s,lk
esll,s
sfbeam,all,1,pres,q
allsel,all
fini
/solu
solve
/post1
pldisp,1
fini$/clear$/prep7
l1=1.4$t=0.02$l2=1$r=0.1$p=20000$q=3000
et,1,shell63,,,2$et,2,beam189
mp,ex,1,2.1e11$mp,prxy,1,0.3
r,1,t
sectype,1,beam,csolid
secdata,r
blc5,,,l1,l1
wprota,,90
asbw,all
wprota,,,90
asbw,all
wpcsys,-1
k,50,,,l2+t/2
kp0=kp(0,0,0)
l,50,kp0
lsel,s,loc,z,0.1,10
latt,1,,2,,,,1$lesize,all,,,4
lmesh,all
aatt,1,1,1
esize,0.1
mshpe,0$mshkey,1
amesh,all
dk,50,all
fk,1,fy,p
fk,4,fx,p
fk,3,fy,-p
fk,2,fx,-p
fk,4,fz,-q
allsel,all
/pbc,f,,2
/pbc,u,,2
ftran
dtran
/eshape,1
ce,1,0,143,ux,1,23,ux,-1,2,rotz,ny(143)-ny(23)
ce,2,0,92,uy,1,30,uy,-1,2,rotz,-(nx(92)-nx(30))
fini
/solu
allsel
solve
fini
/post1
pldisp,1
fini$/clear$/prep7
l=500$b1=200$b2=12$t=2$h=15
q1=0.01$q2=10
et,1,shell63$et,2,beam188
mp,ex,1,2e5$mp,prxy,1,0.3
r,1,t
sectype,1,beam,rect
secdata,b2,h,3,3
secoffset,user,0,h/2+t/2
wprota,,90
rectng,-b1/2,b1/2,0,l
wprota,,,90
asbw,all
wpcsys,-1
lsel,s,tan1,x
lesize,all,,,50
lsel,s,tan1,z
lesize,all,,,10
aatt,1,1,1
amesh,all
k,100,,100,100
lsel,s,loc,x,0
latt,1,,2,,100,,1
lmesh,all
/eshape,1
lsel,s,loc,z,0
dl,all,,all
lsel,s,loc,x,b1/2
sfl,all,pres,q2
sfa,all,1,pres,q1
allsel,all
fini
/solu
allsel
solve
fini
/post1
pldisp,1
fini$/clear$/prep7
b0=0.1$h0=0.09$t0=0.01$l0=3
et,1,shell63,,,2$et,2,beam188$et,3,mpc184,1
mp,ex,1,2.1e11$mp,prxy,1,0.3
r,1,t0
sectype,1,beam,i
secdata,b0,b0,h0,t0,t0,t0
wprota,,90$blc4,,,b0,l0/2
agen,2,all,,,,h0-t0
wpoff,b0/2$wprota,,,90$blc4,,,h0-t0,l0/2
aglue,all
wpcsys,-1
wpoff,,(h0-t0)/2
wprota,,90
asbw,all
esize,b0/6
mshpe,0$mshkey,1
aatt,1,1,1$amesh,all
kp0=kp(b0/2,(h0-t0)/2,l0/2)
k,100,kx(kp0),ky(kp0),l0
k,200,kx(kp0),l0,l0
lsel,none
l,kp0,100
lesize,all,,,10
latt,1,,2,,,200,,1
lmesh,all
lsel,s,loc,z,l0/2
latt,1,,3
lmesh,all
lsel,s,loc,z,0
dl,all,,all
asel,s,loc,y,h0-t0
sfa,all,1,pres,,5000
asel,s,loc,x,b0/2
sfa,all,1,pres,3000
lsel,s,loc,z,l0/2+0.1,l0
esll,s
sfbeam,all,1,pres,5000*b0
sfbeam,all,2,pres,-3000*(h0-t0)
allsel,all
/eshape,1
fini
/solu
allsel
solve
/post1
pldisp,1
fini$/clear$/prep7
b=10$h=20$l1=100$l2=150
a1=b*h$i1=b*h**3/12$i2=h*b**3/12
et,1,solid95$et,2,beam4
mp,ex,1,3e5$mp,prxy,1,0.3$r,1,a1,i1,i2,b,h
blc5,,,b,h,l1$wprota,0,90
vsbw,all
wprota,,,90
vsbw,all
wpcsys,
-1
kp0=kp(0,0,l1)
k,100,,,l1+l2
l,kp0,100
lsel,s,loc,z,l1+1,l1+l2
latt,1,1,2
lesize,all,,,10
lmesh,all
vatt,1,,1
esize,2.5
mshpe,0
mshkey,1
lsel,s,loc,z,1,l1-1
lesize,all,b/2
vmesh,all
asel,s,loc,z,0
da,all,all
fk,100,fx,100
fk,100,fy,-200
nsel,s,loc,z,l1
cerig,1,all
allsel,all
fini
/solu
solve
fini
/post1
pldisp,1
fini$/clear$/prep7
b=10$h=20$l1=100$l2=150
a1=b*h$i1=b*h**3/12$i2=h*b**3/12
et,1,solid95$et,2,beam4
mp,ex,1,3e5$mp,prxy,1,0.3$r,1,a1,i1,i2,b,h
blc5,,,b,h,l1$wprota,0,90
vsbw,all
wprota,,,90
vsbw,all
k,99,,,l1
k,100,,,l1+l2
l,99,100
lsel,s,loc,z,l1+1,l1+l2
latt,1,1,2
lesize,all,,,10
lmesh,all
vatt,1,,1
esize,2.5
mshpe,0
mshkey,1
lsel,s,loc,z,1,l1-1
lesize,all,b/2
vmesh,all
asel,s,loc,z,0
da,all,all
fk,100,fx,100
fk,100,fy,-200
nsel,s,loc,z,l1
cerig,1,all
allsel,all
fini
/solu
solve
fini
/post1
pldisp,1
fini$/clear$/prep7
h=2.4$b=2$l1=1.8$l2=3.0$t=0.3
q1=10000$q2=50000
et,1,solid95$et,2,shell93
mp,ex,1,3e10$mp,prxy,1,0.2
r,1,t
blc4,,,b,h,l1
wpoff,,h-t/2,l1
wprota,,90
vsbw,all
blc4,,,b,l2
nummrg,all
vatt,1,,1
esize,t/2
vmesh,all
asel,s,loc,z,l1+0.1,l1+l2
aatt,1,1,2
amesh,all
sfa,all,1,pres,q1
lsla,s
lsel,r,loc,x,b
sfl,all,pres,q2*t
asel,s,loc,y,0
da,all,all
allsel,all
nsel,s,loc,z,l1
nsel,r,loc,y,h-t,h
cm,ne,node
*do,i,1,29
cmsel,s,ne
nsel,r,loc,x,nx(i)
cerig,i,all,ux,uz
*enddo
allsel,all
fini
/solu
solve
fini
/post1
pldisp,1
fini$/clear$/prep7
b=300$h=30$l=2000$q1=0.1$q2=2
et,1,solid95$et,2,shell93
et,3,beam189
mp,ex,1,2.1e5$Mp,prxy,1,0.3
r,1,h
sectype,1,beam,rect
secdata,b,h
blc4,,,b,h,l/4
wpoff,,h/2
wprota,,90
blc4,,l/4,b,l/4
vsbw,all
nummrg,all
wpoff,b/2
wprota,,,90
vsbw,all
asbw,all
k,100,b/2,h/2,l$k,200,b/2,l/2,l/2
l,kp(b/2,h/2,l/2),100
esize,h/2
vatt,1,,1
vmesh,all
asel,s,loc,z,l/4+0.1,l/2
aatt,1,1,2
amesh,all
lsel,s,loc,z,l/2+1,l
latt,1,,3,,200,,1
lmesh,all
asel,s,loc,y,h
sfa,all,1,pres,q1
lsel,s,loc,x,0
lsel,r,loc,z,l/4+1,l/2
sfl,all,pres,q2*h
lsel,s,loc,z,l/2+1,l
esll,s
sfbeam,all,1,pres,q1*b
sfbeam,all,2,pres,_q2*h
dk,100,all
asel,s,loc,z,0
da,all,all
allsel,all
nsel,s,loc,z,l/2
cerig,7986,all
nsel,s,loc,z,l/4
cm,n1cm,node
nsel,r,loc,y,h/2
*get,nodtol,node,,count
*dim,nodenum,,nodtol
*get,nodenum(1),node,,num,min
*do,i,2,nodtol
nodenum(i)=ndnext(nodenum(i-1))
*enddo
*do,i,1,nodtol
cmsel,s,n1cm
nsel,r,loc,x,nx(nodenum(i))
cerig,nodenum(i),all,ux,uz
*enddo
allsel,all
fini
/solu
solve
fini
/post1
pldisp,1