12-0框架结构模型

finish
/clear
/prep7
/title,Jiulimingzhu 12f-0
!建立第一层梁柱
k,1
k,2,-8
k,3,-16
k,4,-24
l,1,2
l,2,3
l,3,4

lsel,s,,,1,3
lgen,2,all,,,,,-6
lsel,s,,,4,6
lgen,2,all,,,,,-7
lsel,s,,,7,9
lgen,2,all,,,,,-6

l,13,9
l,9,5
l,5,1
lsel,s,,,13,15,1
lgen,4,all,,,-8
nummrg,all
numcmp,all

lsel,all
lgen,2,all,,,,3

*do,i,1,16
l,i,i+16
*enddo

!复制另一半时在此复制

!建立第一层楼板及筏板
a,17,29,32,20
lsel,s,loc,y,3
asbl,1,all

k,,1,,1
k,,-25,,1
k,,-25,,-20
k,,1,,-20
a,33,36,35,34

lsel,s,loc,y,0
asbl,1,all

nummrg,all
numcmp,all

!定义材料参数
et,1,beam188
mp,ex,1,3.25e10
mp,prxy,1,0.2
mp,dens,1,2700

et,2,shell63
mp,ex,2,3.25e10
mp,prxy,2,0.2
mp,dens,2,2700
r,1,0.11
r,2,0.6

et,3,combin14

!定义梁柱板尺寸并划分
sectype,1,beam,rect
secoffset,user,-0.3
secdata,0.6,0.3

sectype,2,beam,rect
secoffset,user,0.3
secdata,0.6,0.3

sectype,3,beam,rect
secdata,0.6,0.3


l_lesize=0.4
z_lesize=0.5

lsel,s,,,25,36
latt,1,,1,,,,1
lesize,all,l_lesize
lmesh,all

lsel,s,,,37,48
latt,1,,1,,,,2
lesize,all,l_lesize
lmesh,all

lsel,s,,,1,24
latt,1,,1,,,,3
lesize,all,l_lesize
lmesh,all

sectype,4,beam,rect
secoffset,user,-0.15,0.15
secdata,0.6,0.6

lsel,s,,,49
latt,1,,1,,,,4
lesize,all,z_lesize
lmesh,all

sectype,5,beam,rect
secoffset,user,0,0.15
secdata,0.6,0.6

lsel,s,,,50,51
latt,1,,1,,,,5
lesize,all,z_lesize
lmesh,all

sectype,6,beam,rect
secoffset,user,0.15,0.15
secdata,0.6,0.6

lsel,s,,,52
latt,1,,1,,,,6
lesize,all,z_lesize
lmesh,all

sectype,7,beam,rect
secoffset,user,0.15,0
secdata,0.6,0.6

lsel,s,,,56,60,4
latt,1,,1,,,,7
lesize,all,z_lesize
lmesh,all

sectype,8,beam,rect
secoffset,user,0.15,-0.15
secdata,0.6,0.6

lsel,s,,,64
latt,1,,1,,,,8
lesize,all,z_lesize
lmesh,all

sectype,9,beam,rect
secoffset,user,0,-0.15
secdata,0.6,0.6

lsel,s,,,62,63
latt,1,,1,,,,9
lesize,all,z_lesize
lmesh,all

sectype,10,beam,rect
secoffset,user,-0.15,-0.15
secdata,0.6,0.6

lsel,s,,,61
latt,1,,1,,,,10
lesize,all,z_lesize
lmesh,all

sectype,11,beam,rect
secoffset,user,-0.15
secdata,0.6,0.6

lsel,s,,,53,57,4
latt,1,,1,,,,11
lesize,all,z_lesize
lmesh,all

sectype,12,beam,rect
secdata,0.6,0.6

lsel,s,,,54,55
lsel,a,,,58,59
latt,1,,1,,,,12
lesize,all,z_lesize
lmesh,all

asel,s,loc,y,3
aatt,2,1,2
amesh,all

asel,s,loc,y,0
esize,0.5
mshape,0
mshkey,0
aatt,2,2,2
amesh,all

!复制1到12层
lsel,s,loc,y,0.1,3
lgen,12,all,,,,3
asel,s,loc,y,3
agen,12,all,,,,3
nummrg,all
numcmp,all

!加弹簧

y1=0

nsel,s,loc,y,y1

*Get,nnod,NODE,0,COUNT !得到底部节点数量

*Dim,xyz,array,nnod,3
*Dim,nodes,array,nnod

*Get,nd,NODE,0,NUM,MIN
*Do,I,1,nnod,1
nodes(I)= nd
xyz(I,1) = Nx(nd)
xyz(I,2) = Ny(nd)
xyz(I,3) =nz(nd)
nd = Ndnext(nd)
*Enddo


r,4,28.47e6 !50cm砂石垫层
real,4
type,3
*Get,nd1,N

ODE,0,NUM,MIN
*Do,I,1,nnod,1
n,,xyz(i,1),xyz(I,2)-0.2,xyz(I,3)
*Get,nd2,NODE,0,NUM,max
e,nodes(i),nd2
*enddo

allsel,all

!加荷载
/solu
acel,,13.23
allsel,all
aplot
*do,i,0,12
lou_h=3*i

asel,s,loc,y,lou_h
sfa,all,1,pres,-4.825e3
*enddo


!fu曲率变形,圆心在边上
nsel,s,loc,y,y1-0.2

*Get,nnod,NODE,0,COUNT !得到底部节点号

*Dim,xyz,array,nnod,3
*Dim,nodes,array,nnod

*Get,nd,NODE,0,NUM,MIN
*Do,I,1,nnod,1
nodes(I)= nd
xyz(I,1) = Nx(nd)
xyz(I,2) = Ny(nd)
xyz(I,3) =nz(nd)
nd = Ndnext(nd)
*Enddo

*Get,nd1,NODE,0,NUM,MIN

*Do,I,1,nnod,1
m_L=26
quR=5000
!yyy=sqrt(quR**2-m_L**2)-sqrt(quR**2-m_L**2-xyz(i,1)**2+2*m_L*xyz(i,1))
!yyy=-1*yyy
D,nd1,uy,yyy
d,nd1,ux,0
d,nd1,uz,0
nd1 = Ndnext(nd1)
*Enddo

allsel,all


solve