临时支撑CD法.doc

!修改网格及参数（Ⅴ级围岩） finish /clear *afun,deg x_zs=-50 y_zs=30 x_yx=50 y_yx=-40 num=24 th_s=0 th_f=180-th_s th=(th_f-th_s)/num yx_1_x=0 yx_1_y=0 yx_2_x=0 yx_2_y=11.37 yx_3_x=3.34 yx_3_y=0 yx_4_x=-3.34 yx_4_y=0 r_1=5.8 r_2=14.4 r_3=2.67 r_4=2.67 r_jg_1=8.8 r_jg_2=17.4 r_jg_3=5.67 r_jg_4=5.67 !beam3单元厚度（考虑钢纤维、二衬提高0.5） t1=0.25*1.5 !link1单元半径 r1=12.5e-3 !plane42单元 e_r=0.2e9 u_r=0.35 dens_r=1950 c_r=0.18e6 fi_r=24 !beam3单元 e_c=23e9 u_c=0.2 dens_c=2200 !link1单元 e_a=210e9 u_a=0.25 dens_a=7700 /prep7 !单元类型 et,1,beam3 et,2,plane42,,,2 keyopt,2,3,2 et,3,link1 !实常数 r,1,t1,t1*t1*t1/12,t1 r,2,3.14*r1*r1 !材料属性 !左上支护 mp,ex,1,e_c mp,prxy,1,u_c mp,dens,1,dens_c !右上支护 mp,ex,2,e_c mp,prxy,2,u_c mp,dens,2,dens_c !左下支护 mp,ex,3,e_c mp,prxy,3,u_c mp,dens,3,dens_c !右下支护 mp,ex,4,e_c mp,prxy,4,u_c mp,dens,4,dens_c !左上开挖土体 mp,ex,5,e_r mp,prxy,5,u_r mp,dens,5,dens_r tb,dp,5 tbdata,1,c_r,fi_r !右上开挖土体 mp,ex,6,e_r mp,prxy,6,u_r mp,dens,6,dens_r tb,dp,6 tbdata,1,c_r,fi_r !左下开挖土体 mp,ex,7,e_r mp,prxy,7,u_r mp,dens,7,dens_r tb,dp,7 tbdata,1,c_r,fi_r !右下开挖土体 mp,ex,8,e_r mp,prxy,8,u_r mp,dens,8,dens_r tb,dp,8 tbdata,1,c_r,fi_r !左上锚杆 mp,ex,9,e_a mp,prxy,9,u_a mp,dens,9,dens_a !右上锚杆 mp,ex,10,e_a mp,prxy,10,u_a mp,dens,10,dens_a !左临时支撑 mp,ex,11,e_c mp,prxy,11,u_c mp,dens,11,dens_c !上临时支撑 mp,ex,12,e_c mp,prxy,12,u_c mp,dens,12,dens_c !右临时支撑 mp,ex,13,e_c mp,prxy,13,u_c mp,dens,13,dens_c !下临时支撑 mp,ex,14,e_c mp,prxy,14,u_c mp,dens,14,dens_c !围岩 mp,ex,15,e_r mp,prxy,15,u_r mp,dens,15,dens_r tb,dp,15 tbdata,1,c_r,fi_r !实体建模 k,1,yx_1_x,yx_1_y k,2,yx_2_x,yx_2_y k,3,yx_3_x,yx_3_y k,4,yx_4_x,yx_4_y k,101,yx_1_x,yx_1_y,10 k,102,yx_2_x,yx_2_y,10 k,103,yx_3_x,yx_3_y,10 k,104,yx_4_x,yx_4_y,10 circle,1,r_1,101 circle,2,r_2,102 circle,3,r_3,103 circle,4,r_4,104 lcsl,all lsel,s,,,17,22,1 lsel,a,,,46,48,1 lsel,a,,,27,28,1 lsel,a,,,43 lsel,inve ldel,all lsel,all lcomb,17,18 lcomb,19,21 lcomb,43,27 lcomb,28,46 cm,zh,line circle,1,r_jg_1,101 circle,2,r_jg_2,102 circle,3,r_jg_3,103 circle,4,r_jg_4,104 cmsel,s,zh,line lsel,inve lcsl,a