mabizhucenan1 mabizhucenan1
关注数: 4 粉丝数: 63 发帖数: 423 关注贴吧数: 24
本人在做边坡动力液化分析,希望大家看看命令流有问题吗? 初始分析: impgrid deceng.flac3d config dyn fluid model mohr prop bulk 3.39e7 shear 1.56e7 coh 1.4e4 fri 25 ran gro soil prop bulk 8e8 shear 5.3e8 coh 1.5e5 fri 40 ran gro rock ini dens 1430 ran gro soil;计算干密度 ini dens 2307 ran gro rock model fl_iso;流体模型 prop poro 0.53 perm 7.8e-11 ran gro soil prop poro 0.128 perm 1.9e-12 ran gro rock ini fmod 2e8 fdens 1000 ini fmod 0 impgrid shuiwei.flac3d;建立水位面 gen separate water_below interface 2 wrap water_below water_above set grav 0 0 -9.81 water den 1000 def water_table p_i=i_head p_ie=i_elem_head(p_i) loop while p_ie # null p_gp1=ie_vert(p_ie,1) p_gp2=ie_vert(p_ie,2) p_gp3=ie_vert(p_ie,3) x1=in_pos(p_gp1,1) y1=in_pos(p_gp1,2) z1=in_pos(p_gp1,3) x2=in_pos(p_gp2,1) y2=in_pos(p_gp2,2) z2=in_pos(p_gp2,3) x3=in_pos(p_gp3,1) y3=in_pos(p_gp3,2) z3=in_pos(p_gp3,3) command water table face x1,y1,z1 x2,y2,z2 x3,y3,z3 endcommand p_ie=ie_next(p_ie) endloop end model null range group water_below model null range group water_above water_table int 2 dele dele range group water_below dele range group water_above gen separate soil;建立滑动面接触面 int 1 wrap soil rock int 1 prop ks 3.4e8 kn 3.4e8 fric 24 coh 1.2e4 fix x y z ran z -0.1 0.1;边界条件 fix y ran y -0.1 0.1 fix y ran y 689.9 690.1 fix x ran x -0.1 0.1 fix x ran x 0.9 1.1 set grav 0 0 -9.81 set fluid off dyn off hist id=1 unbal solve elas save chushi.sav 动力分析: restore chushi.sav set dyn on fluid on set dyn multi on;动态多步 set large;大变形 ini fmod 2e8 set fluid pcut on model finn ran gro soil;finn模型 prop bulk 3.39e7 shear 1.56e7 coh=1.4e4 fric=25 ran gro soil;finn参数 ini dens 1430 ran gro soil prop ff_latency=50 ff_switch = 0 ff_c1=0.8 ff_c2=0.79 ff_c3=0.45 ff_c4=0.73 ran gro soil ini xvel 0 yvel 0 zvel 0;位移速度设为0 ini xdis 0 ydis 0 zdis 0 ini state 0 fix x y z ran z -0.1 0.1 fix y ran y -0.1 0.1 fix y ran y 689.9 690.1 fix x ran x -0.1 0.1 fix x ran x 0.9 1.1 free y ran z -0.1 0.1 table 1 read 1.txt;y方向施加地震波 apply yacc=9.8 hist table 1 ran z -0.1 0.1 set dyn damp local 0.224;局部阻尼 apply ff;自由场边界 group ff_corner group ff_side ran x 0 1 group ff_side ran y 0 690 group main_grid ran x 0 1 y 0 690 config zextra 20;监测动孔压数据 def get_old_stress p_z = zone_head loop while p_z # null z_extra(p_z,5) = z_pp(p_z) ;初始孔压 p_z = z_next(p_z) endloop end get_old_stress def get_PPR whilestepping p_z = zone_head loop while p_z # null z_extra(p_z,7) = z_szz(p_z) + z_pp(p_z) ;竖向有效应力 z_extra(p_z,4) = z_szz(p_z) ;竖向应力 z_extra(p_z,6) = z_pp(p_z) - z_extra(p_z,5) ;超静孔隙水压力 z_extra(p_z,9) = z_szz(p_z) z_extra(p_z,2) = (-z_extra(p_z,6) ) / (z_extra(p_z,4)+ z_extra(p_z,5)) ;超孔压比 z_extra(p_z,3) = 1.0 - z_extra(p_z,7) / (z_extra(p_z,4)+ z_extra(p_z,5)) ;超孔压比 p_z = z_next(p_z) endloop end get_PPR def get_ppr_zone whilestepping p_z1 = z_near( 1,3.115e2,8.102e1 );监测点 ppr_11 = z_extra(p_z1,2) ppr_12 = z_extra(p_z1,3) ppr_13 = z_extra(p_z1,6) ppr_14 = z_extra(p_z1,7) end hist id 2 dytime hist id 3 ppr_11 hist id 4 ppr_12 hist id 5 ppr_13 hist id 6 ppr_14 set dyn time=0 solve age 10
1 下一页