/TITLE, shield tunnel in Metro /NOPR
/COM, Structural
/PREP7
ET,1,BEAM3 !管片
ET,2,link10,,,1, !地基弹簧1
ET,3, COMBIN14,0,0, ! 地基弹簧2
ET,4,COMBIN14,0,6, !环间抗弯弹簧1外侧受拉 ET,5,COMBIN14,0,6, !环间抗弯弹簧2内侧受拉 ET,6, COMBIN14,0,1,0 !1环纵向抗剪弹簧径向 ET,7, COMBIN14,0,2,0 ! 1环纵向剪切弹簧切向 ET,8, COMBIN14,0,1,0 !2环纵向抗剪弹簧径向 ET,9, COMBIN14,0,2,0 ! 2环纵向剪切弹簧切向 ET,10,link10,,,1, !地基弹簧3
R,1,0.3, 0.00225,0.3, , , , R,2,0.0497, , , R,3,1e6, , , R,4,3e7, R,5,5e7, R,6, 2.79e10, R,7, 2.79e10, R,8, 2.79e10, R,9, 2.79e10, R,10,0.0497, , ,
MPTEMP,1,0 MPDATA,EX,1,,50e9 MPDATA,PRXY,1,,0.2 MPTEMP,1,0
MPDATA,DENS,1,,2500 MPDATA,EX,2,,170e6 MPDATA,EX,10,,170e6 SAVE
K,1000,0,0,0, csys,1
K,1,2.85,99,0 K,2,2.85,115,0 K,3,2.85,115,0 K,4, 2.85,180,0 K,5, 2.85,180,0 K,6, 2.85,252,0
K,7, 2.85,252,0 K,8, 2.85,324,0 K,9, 2.85,324,0 K,10, 2.85,36,0 K,11, 2.85,36,0 K,12,2.85,99,0 SAVE
Larc,1,2,1000,2.85 Larc,3,4,1000,2.85 Larc,5,6,1000,2.85 Larc,7,8,1000,2.85 Larc,9,10,1000,2.85 Larc,11,12,1000,2.85
Lgen,2,1,6,1,,-36,1.5,2000,1,0 Lgen,2,1,6,1,,- 36,-1.5,4000, 1,0
SAVE csys,1
lsel,all
LESIZE,all, ,1 ,, ,1 , , , TYPE,1 MAT,1 REAL,1 lmesh,1,6,
lsel,all
numstr,node,2001 lmesh,7,12, lsel,all
numstr,node,4001 lmesh,13,18,
/PNUM,ELEM,1 SAVE
!生成地弹簧1 csys,1
Nsel,s,loc,x,2.85
ngen,2,1000, all, ,,1 !生成节点 TYPE,2 MAT,2
REAL,2
Nsel,s,loc,x,3.85 Nsel,r,loc,y,-45,235 *get,na,node,0,count *get,nmin,node,,num,min bnum=nmin *do,m,1,na,
E,bnum,bnum-1000 bnum=ndnext(bnum) *ENDDO
TYPE,10 MAT,10 REAL,10
Nsel,s,loc,x,3.85 Nsel,r,loc,y,235,315 *get,na,node,0,count *get,nmin,node,,num,min bnum=nmin *do,m,1,na,
E,bnum,bnum-1000 bnum=ndnext(bnum) *ENDDO
!生成地弹簧2 Nsel,s,loc,x,2.85
ngen,2,1500,all, ,,0.043965844,10 TYPE,3 REAL,3
NSEL,s,node,,1501,1870 *get,na,node,0,count *get,nmin,node,,num,min bnum=nmin *do,m,1,na,
E,bnum,bnum-1500 bnum=ndnext(bnum) *ENDDO
NSEL,s,node,,3501,3870 *get,na,node,0,count *get,nmin,node,,num,min bnum=nmin
*do,m,1,na,
E,bnum,bnum-1500 bnum=ndnext(bnum) *ENDDO
NSEL,s,node,,5501,5870 *get,na,node,0,count *get,nmin,node,,num,min bnum=nmin *do,m,1,na,
E,bnum,bnum-1500 bnum=ndnext(bnum) *ENDDO
!第一个圆 csys,1
!生成环向耦合 Nsel,s,loc,x,2.85 Nsel,r,loc,z,0 cpintf, ux, 0.01, cpintf, uy, 0.01,
TYPE,4 REAL,4
NSEL,S,LOC,x, 2.85 NSEL,r, LOC,y,35.9,36.1 NSEL,r, LOC,z,0,
*get,i,node,0,num,max *get,j,node,0,num,min E,i,j
NSEL,S,LOC,x, 2.85
NSEL,r, LOC,y,179.9,180.1 NSEL,r, LOC,z,0,
*get,i,node,0,num,max *get,j,node,0,num,min E,i,j
TYPE,5 REAL,5
NSEL,S,LOC,x, 2.85 NSEL,r, LOC,y,98.9,99.1 NSEL,r, LOC,z,0,
*get,i,node,0,num,max *get,j,node,0,num,min E,i,j
NSEL,S,LOC,x, 2.85
NSEL,r, LOC,y,114.9,115.1 NSEL,r, LOC,z,0,
*get,i,node,0,num,max *get,j,node,0,num,min E,i,j
NSEL,S,LOC,x, 2.85
NSEL,r, LOC,y,-107.9,-108.1 NSEL,r, LOC,z,0,
*get,i,node,0,num,max *get,j,node,0,num,min E,i,j
NSEL,S,LOC,x, 2.85
NSEL,r, LOC,y,-35.9,-36.1 NSEL,r, LOC,z,0,
*get,i,node,0,num,max *get,j,node,0,num,min E,i,j
!第二个圆 csys,1
Nsel,s,loc,x,2.85 Nsel,r,loc,z,1.5 cpintf, ux, 0.01, cpintf, uy, 0.01,
TYPE,4 REAL,4
NSEL,S,LOC,x, 2.85
NSEL,r, LOC,y,-36+35.9, -36+36.1 NSEL,r, LOC,z,1.5,
*get,i,node,0,num,max *get,j,node,0,num,min E,i,j
NSEL,S,LOC,x, 2.85
NSEL,r, LOC,y, -36+179.9, -36+180.1
6盾构管片三维计算命令流-径向弹簧-不用循环语句加载耦合- 重合节点模拟纵向螺栓
