flac地震响应分析....docx

FLAC3D 边坡地震动力分析?? 2011-05-12 18:48:00|??分类：? HYPERLINK /blog/ \l m=0t=1c=fks_084071085084084067083080083095085095088071085094082069083 \o 默认分类 默认分类|举报|字号?订阅 DOCTYPE html PUBLIC -WCDTD XHTML TransitionalEN httpwwwworgTRxhtmlDTDxhtml-transitionaldtd 最近做的一个边坡地震动力分析，鉴于目前这一块资料实在缺少，拿出来与大家一块探讨一下！共同进步！（我也刚开始做这一块，不能保证完全正确，请多提意见!） ;静力计算部分 new config?dyn set?dyn?off gen?zone?brick?p0?0?0?0?p1?100?0?0?p2?0?2?0?p3?0?0?40?size?50?1?10 gen?zone?brick? p0?40?0?40?p1?100?0?40?p2?40?2?40?p3?60?0?60? p4?100?2?40?p5?60?2?60?p6?100?0?60?p7?100?2?60? size?30?1?10 plot?surf ;初始应力场生成 model?ela pro?bulk?1e9?shear?2e9?dens?2500 fix?y fix?x?range?x?-.1?.1 fix?x?range?x?99.9?100.1 fix?x?y?z?range?z?-.1?.1 set?gravity?0?0?-10 set?mech?ratio?1e-3 solve nbsp; ini?xvel?0?yvel?0?zvel?0 ini?xdis?0?ydis?0?zdis?0 nbsp; model?mohr pro?bulk?8.3e7?shear?3.8e7?coh?4.2e4?dila?20?ten?0.8e6?fric?17?dens?2500 nbsp; save?dizhen1.sav nbsp; restore?dizhen1.sav ;动力计算部分 def?iniwave freq?=5.0 end Iniwave ;定义动荷载（以速度形式） def?wave vv=9.36e-2*sin(2.0*pi*freq*dytime) ?;wave=0.5*(1.0-cos(2*pi*dytime/per)) ?;wave=sin(2*pi*dytime/per) if?dytime0.2 ?wave=dytime/0.2*vv ?else ??if?dytime2.0 ????wave=vv ???else? ?????if?dytime?=3.0 ???????wave=(3.0-dytime)/1.0*vv ?????endif ???endif endif ?if?dytime3.0 ?wave=0.0 endif end nbsp; set?dyn?on free?x?y?z?ran?z?-.1?.1??;去除静力边界 apply?nquiet?squiet?dquiet?ran?z?-.1?.1??;在底部施加静态边界 apply?xvel?1.0?hist?wave?range?z?-.1?.1??;施加动荷载 apply?ff;施加自由场边界 nbsp; pause?k hist?id?=1?gp?xvel?15?0?0 hist?id?=2?gp?xvel?15?0?10 hist?id?=3?gp?xvel?15?0?20 hist?id?=4?gp?xvel?15?0?30 hist?id=5?wave hist?id=6?dytime plot?create?hhh plot?show set?dyn?time=0 set?dyn?damp?local?.314??;局部阻尼 solve?age?4 plot?hist?5?vs?6? pause?K plot??1?2?3?4?vs?6 nbsp; ? Fig.1水平位移云图（最大水平位移53cm） nbsp; nbsp; ? Fig.2?剪切速率云图及速度矢量图 ? ? Fig.3地震波速度时程图（5vs6） ? Fig.4斜坡内各点速度时程图（1 2 3 4 vs6）