变截面连续刚构实体建模.docVIP

变截面连续刚构实体建模 !----------------- !author:KooFong !email:kooyukfong@ !----------------- finish$/clear$/prep7 et,1,solid95 !定义混凝土为solid95单元，参考号为1 et,2,link8 !定义钢束为link8单元，参考号为2 ec=3.45e10 !混凝土弹性模量值 es=1.95e11 !钢束弹性模量值 mp,ex,1,ec !定义混凝土弹性模量 mp,prxy,1,0.1667 !定义混凝土泊松比 mp,dens,1,2600 !定义混凝土密度 mp,ex,2,es !定义钢筋弹性模量 mp,prxy,2,0.3 !定义钢筋泊松比 mp,dens,2,7850 !定义钢筋密度 area_s=1.39e-4 !钢束面积 Fs=193905 !钢束张拉力 L=0.5 !钢束单元长度 strain=Fs/es/area_s !钢束初应变 r,2,area_s,strain !钢束实常数 *dim,z,array,26 *dim,h,array,26 *do,loop,1,26,1 *if,loop,le,12,then z(loop)=0+(loop-1)*3.5 !左端抛物线方程 h(loop)=-8/5929*z(loop)**2-2.5 *elseif,loop,eq,13,then z(loop)=39.5 !桥墩处 h(loop)=-4.5 *elseif,loop,eq,14,then z(loop)=40.5 !桥墩处 h(loop)=-4.5 *else z(loop)=41.5+(loop-15)*3.5 !左端抛物线方程 h(loop)=-8/5929*(z(loop)-80)**2-2.5 *endif *enddo *do,loop,1,26,1 k,1+(loop-1)*20,-6,-0.35,z(loop) !循环建立节点 k,2+(loop-1)*20,-6,0,z(loop) k,3+(loop-1)*20,-3.25,0,z(loop) k,4+(loop-1)*20,-2.45,0,z(loop) k,5+(loop-1)*20,2.45,0,z(loop) k,6+(loop-1)*20,3.25,0,z(loop) k,7+(loop-1)*20,6,0,z(loop) k,8+(loop-1)*20,6,-0.35,z(loop) k,9+(loop-1)*20,3.25,-0.45,z(loop) k,10+(loop-1)*20,3.25,0.45+h(loop),z(loop) k,11+(loop-1)*20,3.25,h(loop),z(loop) k,12+(loop-1)*20,2.45,h(loop),z(loop) k,13+(loop-1)*20,-2.45,h(loop),z(loop) k,14+(loop-1)*20,-3.25,h(loop),z(loop) k,15+(loop-1)*20,-3.25,0.45+h(loop),z(loop) k,16+(loop-1)*20,-3.25,-0.45,z