wpcsys,-1,0
wpoffs,-mod_w1,-(mod_w2+hn_sec_h)*0.5,mod_h-hn_sec_t2*0.5 asel,r,loc,z,
wpoffs,,,-hn_sec_h asel,a,loc,z, cm,hnsect2,area
AATT, 2, 5, 2, 0, !!!!梁翼缘
cmsel,s,gangjiegou,area cmsel,u,hnsect2,area wpcsys,-1,0
wpoffs,-mod_w1,-mod_w2*0.5,mod_h-hn_sec_t2*0.5 asel,r,loc,y,hw_sec_w*0.5+0.001,mod_w2 cm,hnsect1,area
AATT, 2, 3, 2, 0, !!!!!梁腹板
cmsel,s,gangjiegou,area cmsel,u,hnsect2,area cmsel,u,hnsect1,area cm,hwall,area asel,r,loc,x,
cm,hwsect1,area
AATT, 2, 2, 2, 0, !!!!!柱子腹板
cmsel,s,hwall,area cmsel,u,hwsect1,area
AATT, 2, 4, 2, 0, !!!!!柱子翼缘
!!!!钢结构网格
cmsel,s,gangjiegou,area ESIZE,gridding,0, AMESH,all
nummrg,all
!!!!!!!选钢筋切体 wpcsys,-1,0 allsel,all
KWPAVE, 33 wpoffs,,,-20
vsbw,all
wpcsys,-1,0
KWPAVE, 33 wpoffs,20, wpro,,,90 vsbw,all
*do,i,1,mod_w1/steelbar,1
wpoffs,,,steelbar !!!!!!改了 vsbw,all *enddo
!!!!!!!!!!!!!!!!!!!! !体模型save.amesh !!!!!!!!!!!!!!!!!!!! !选钢筋的线 wpcsys,-1,0 allsel,all
KWPAVE, 33 wpoffs,20,,-20
lsel,r,loc,z
cm,ltemp,line !选出那一层的线 lsel,r,loc,x,0 cm,ltemp1,line
cmsel,s,ltemp,line
*do,i,1,mod_w1/steelbar,1 lsel,r,loc,x,i*steelbar cmsel,a,ltemp1,line cm,ltemp1,line cmsel,s,ltemp,line *enddo
cmsel,s,ltemp1,line LATT,3,7,5, , , , ESIZE,gridding,0, LMESH,all
!!!!!!!!!!!!!!!!!体单元
allsel,all
ESIZE,gridding,0,
VATT, 1, 1, 1, 0 !体属性 vmesh,all
!**************************************************** !加弹簧、耦合 !上混泥土板
!*************************** type,3 !设定弹簧单元 real,6
wpcsys,-1,0
wpoffs,-mod_w1, nsel,s,loc,x
nsel,r,loc,y,0,-(mod_w2-hw_sec_w)*0.5 nsel,r,loc,z,mod_h-0.001,mod_h+0.001 cm,enod,node
*get,max0,node,0,count !max1=15 *dim,ojd0,,max0 *dim,jd0,,max0
*get,nod0,node,0,num,min ojd0(1)=nod0 *do,i,2,max0
ojd0(i)=ndnext(ojd0(i-1)) *enddo allsel,all nsel,all
cmsel,u,enod *do,i,1,max0 nod0=ojd0(i) j0=nnear(nod0) jd0(i)=j0 *enddo nsel,all
*do,i,1,max0 e,ojd0(i),jd0(i) *enddo allsel,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!node(x,y,z) !返回最近节点的编号
!!1400 !nnear(n) !返回最接近n的节点编号 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!耦合纵梁节点 wpcsys,-1,0
wpoffs,-mod_w1,
!nsel,r,loc,y,0,-2400 !!1400
nsel,r,loc,z,mod_h-0.001,mod_h+0.001 nsel,r,loc,x,-hn_sec_w*0.5,hn_sec_w*0.5 cm,cmljnod,node
*get,max1,node,0,count !max1=15 *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,cmljnod *do,i,1,max1 nod1=ojd(i) j=nnear(nod1) jd(i)=j *enddo nsel,all ji = 1
*do,i,1,max1
cp,ji,ux,ojd(i),jd(i) cp,ji+1,uz,ojd(i),jd(i) ji=ji+2 *enddo allsel,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!node(x,y,z) !返回最近节点的编号
!nnear(n) !返回最接近n的节点编号 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!耦合端部节点 wpcsys,-1,0
wpoffs,-mod_w1,
nsel,s,loc,x,-hn_sec_w*0.5,hn_sec_w*0.5
nsel,r,loc,y,-(mod_w2+hw_sec_w)*0.5 ,-(mod_w2-hw_sec_w)*0.5 !!1400
nsel,r,loc,z,mod_h-0.001,mod_h+0.001 cm,egnod,node
*get,max2,node,0,count !max1=15 *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,egnod *do,i,1,max2 nod2=ojd2(i) j2=nnear(nod2) jd2(i)=j2 *enddo nsel,all ji2 = ji + 1 *do,i,1,max2
cp,ji2,uy,ojd2(i),jd2(i) ji2=ji2+1 *enddo allsel,all
!!!!!!!!!!!载荷
wpcsys,-1,0 nsel,s,loc,z, D,all, , , , , ,ALL, , , , , !ACEL,0,0,0, nsel,s,loc,y,
D,all, , , , , ,UY, , , , , allsel,all
wpcsys,-1,0
wpoffs,-mod_w1,-(mod_w2-hw_sec_h)*0.5,mod_h-hn_sec_t2*0.5 nsel,r,loc,z, nsel,r,loc,x,
nsel,r,loc,y,0,mod_h-hn_sec_t2*0.5 cm,N_load,node allsel,all
wpcsys,-1,0
!cmsel,s,N_load,node !*get,Nnod,node,0,count !*Dim,Nodes,array,Nnod !*get,Nd,node,0,num,min !*do,I,1,Nnod,1 ! Nodes(I)=Nd
! F,Nd,FZ,-11000 !!!!!2000-20000 ! Nd=NDNEXT(Nd)
!*ENDDO !!SFTRAN !allsel,all !wpcsys,-1,0 ! !/sol
!cnvtol,f,,0.05,2 !定义收敛条件,使用缺省的VALUE !!nsubst,50 !定义子步数
!outres,all,all !输出每一子步的结果 !autots,1 !打开自动时间步控制 !lnsrch,1 !打开线性搜索
!ncnv,2 !如果不收敛时结束而不退出
!!neqit,50 !每一子步中方程的迭代次数限值 !pred,on !打开预测器 !
!ANTYPE,0 !NLGEOM,1
!NSUBST,100,0,0 !OUTRES,ERASE !OUTRES,ALL,LAST !AUTOTS,-1.0 !PSTRES,1 !TIME,20