!双层两跨框架结构!建模finish/clear/title,nanhudadaozhan!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!几何参数设置（根据工程修改）!顶板厚度d1=0.8!楼板厚度d2=0.4!底板厚度d3=0.9!两侧墙厚度d4=0.7!支柱等效厚度zhuchang=1.0!!!!柱长zhukuan=0.7!!!!柱宽zhuju=8!!!!柱距d5=(zhuchang*zhukuan**3/zhuju)**(1/3)!围护结构等效厚度d6=0.8!跨度w1=10!!!支柱离左侧墙距离nw1=10!!!划分数w2=10!!!支柱离右侧墙距离nw2=10!!!划分数!中板距顶板h1=5nh1=5!!!划分数!中板距底板h2=5nh2=5!!!划分数!底板距连续墙底h3=7nh3=7!!!划分数!顶板距连续墙顶h4=2nh4=2!!!划分数!围护结构与侧墙距离（侧墙单元别小于nn，它用来选择单元来用的）nn=0.9!!!!!!!!!!!!!!!!荷载参数（根据工程修改）!顶板水土压力加超载p1=80*1000!楼板荷载，恒载加活载p2=6*1000!底板水压p3=130*1000!围护结构顶水平土压qt1=10*1000!围护结构底水平土压,可以将土从中板处分层两个线性荷载来施加qt2=64*1000!围护结构底水平土压qt3=100*1000!侧墙顶水压qw1=30*1000!侧墙底水压qw2=p3!!!!!!!!!!!!!!!!!!物理参数!c30!衬砌容重r1=25e3*1!衬砌弹性模量e1=30e9!衬砌泊松比u1=0.2!c40!衬砌容重r2=25e3*1!衬砌弹性模量e2=32.5e9!衬砌泊松比u2=0.2!围岩弹性抗力系数,和单元划分细密有关，尽量将单元划分为1米长k1=10e6!!!!底板竖向基床系数k3=12e6!!!!水平基床系数!链杆单元弹性模量，按C30取k2=1e13!!!!!可以取个大数吧????!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义单元类型、实常数、材料属性。/prep7!定义梁单元et,1,beam3!定义链杆单元et,2,link10keyopt,2,3,1!设为只受压!定义弹簧单元et,3,combin14!定义实常数!定义梁单元的面积、惯性矩和梁高r,1,d1,d1*d1*d1/12,d1r,2,d2,d2*d2*d2/12,d2r,3,d3,d3*d3*d3/12,d3r,4,d4,d4*d4*d4/12,d4r,5,d5,d5*d5*d5/12,d5r,6,d6,d6*d6*d6/12,d6!!!围护结构若考虑刚度折减，则在此惯性矩可乘个系数!定义弹簧单元的弹性系数r,7,k1!!!!底板竖向基床系数r,9,k3!!!!水平基床系数!定义链杆单元的实常数(面积)r,8,1!定义材料属性!衬砌材料属性,C30mp,ex,1,e1mp,prxy,1,u1mp,dens,1,r1/10!衬砌材料属性,C40mp,ex,2,e2mp,prxy,2,u2mp,dens,2,r2/10!链杆单元属性mp,ex,3,k2!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建立几何模型。!创建隧道支护控制关键点k,1,-w1,0k,2,0,0k,3,w2,0k,4,-w1,h2k,5,0,h2k,6,w2,h2k,7,-w1,h2+h1k,8,0,h2+h1k,9,w2,h2+h1K,10,w2+nn,-h3K,11,w2+nn,0K,12,w2+nn,h2K,13,w2+nn,h2+h1K,14,w2+nn,h2+h1+h4K,15,-w1-nn,-h3K,16,-w1-nn,0K,17,-w1-nn,h2K,18,-w1-nn,h2+h1K,19,-w1-nn,h2+h1+h4!生成线l,1,2l,2,3l,4,5l,5,6l,7,8l,8,9l,1,4l,4,7l,3,6l,6,9l,2,5l,5,8l,10,11l,11,12l,12,13l,13,14l,15,16l,16,17l,17,18l,18,19!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!单元网格划分。!底板lsel,s,line,,1lesize,all,,,nw1latt,2,3,1lmesh,alllsel,s,line,,2lesize,all,,,nw2latt,2,3,1lmesh,all!中板lsel,s,line,,3lesize,all,,,nw1latt,2,2,1lmesh,alllsel,s,line,,4lesize,all,,,nw2latt,2,2,1lmesh,all!顶板lsel,s,line,,5lesize,all,,,nw1latt,2,1,1lmesh,alllsel,s,line,,6lesize,all,,,nw2latt,2,1,1lmesh,all!侧墙lsel,s,line,,7,9,2lesize,all,,,nh2latt,2,4,1lmesh,all!!lsel,s,line,,8,10,2lesize,all,,,nh1latt,2,4,1lmesh,all!中柱lsel,s,line,,11lesize,all,,,nh2latt,2,5,1lmesh,alllsel,s,line,,12lesize,all,,,nh1latt,2,5,1lmesh,all!围护lsel,s,line,,13,17,4lesize,all,,,nh3latt,1,6,1lmesh,alllsel,s,line,,14,18,4lesize,all,,,nh2latt,1,6,1lmesh,alllsel,s,line,,15,19,4lesize,all,,,nh1latt,1,6,1lmesh,alllsel,s,line,,16,20,4lesize,all,,,nh4latt,1,6,1lmesh,all!!!!!!!!!!!!!!!!!!!创建弹簧单元!显示节点nplot!显示节点编号/pnum,node,1/replot!确定单元类型type,3!赋予弹簧单元实常数real,7!创建底板弹簧单元!选择底板节点nsel,s,loc,y,0nsel,r,loc,x,-w1,w2ngen,2,300,all,,,,-0.9*do,i,1,nw1+nw2+1,1e,i,...
