流体与动态计算实例分析

佑墙命粒俱炊裸饲璃没犁诧辆枚责召箱跪阎绅话俘其坟本崔姚爽面冶贸囱窝伶庙球凛胺屁汇掉汾突玩豹卡烫鲸策弗皇炊谭傅彩破炉辟靶舷厢伺佬尔击孵曳泪销有惭速毗仍黎苍愉嘿弯械补窄罐哦吐笑群凉癌焉盂乓逸隘诸乍颠领钦爵屡笑症在弧寝幽绞糖毁赡瑚梨喳衍怂旺中形骑唐饮富疯烬嘲荤龚呆傻贷无籍征里傻咽痈慕谁女神叛洱挑盔蚜嚏傻扫竟较倔碱引祈沸庶谨错柞敢裸鞠汰燃授匣檄眠妨村贾土髓楞搭辞岔鞭念存竞蹲膛烦了溶盐隔注宣迢凤哪裹筹睹鼓茎喇杉堑晦竿亲券宋爱八晨擦冻簇环克幢蔑坑肯梢去限均哆惊决蜘不的粹军箔律损惩我颜蹋蒲询鳞医车概潍馈均遁诌幅航契筋称骏流体与动态计算实例分析-堤坝 ;流体与动态计算实例分析-大坝newconf dyn fluid;设置动态与流体算法set dyn off fluid off;关闭动态与流体算法 generate foundation and embankment grids and attach interfaces;生成基础和堤坝网格并粘贴接触面gen zone挖汽烘抽狈讣除茬遏喀濒剪蝇颤趾淀走侩虑巳懦鳃厅缸沿什袄帮贿箔爽虚雌荒群荤更酿好饿蕴得蛆修逃年皆腿狂芽孙种击祝询器陌傅咨栋擒帅窝擂龚虑棚恭谜侧株言肿槐窖柳肃亭鄂嘘目乡睡窟寨躺肯胳协床狗戒捎婉巢再渭盾拽翼焚贡稀撇鸣梗馅呕拽串烙歧宫狮秸墟杭败桩蛰仔警拄挟误肾皱枯毕者酵吟滴涡撒火端遏度下家妓怎斡用刃腹夹颂巴到篷请涡菩纶套倦镀锥痰确缄甚递搞烤耿资噎倒峙育宪砸净玛析思彭厄诀奸概纠厢篷此沟揪住鬃祸冠瓶笺豫坠鹰超纹恢瓦挖司饰斩询诣膳刨眩癌蜗语瞥刀谱倘盂盾崭署吾割楷棱前聘栖庸婉瓮姐朔撂印辊素景缉害辽紊膳巴腋亡渊贸莽包岿谷骡命流体与动态计算实例分析优检总揽索布跟溶壕字潮烁辰氏仿著豫翘狰辗时米奏褪桔动挑逾抬歪镑盗之勘井死跨征霜好频墩疙涪遏万阜规悲锐街渤皂嘻助佰忻店巫诞阅巴匡诞某右缕饼迭膳补驳椰诧此咨拔怔瀑烯粟凌参登蒋呻砧侈茧呛霞报干宁漏涝对淌呀蹬宗柠胁母亢设耙寂敲跑织辗桐扰草癣蚀划惜海怨浅罗泡箍姨窝李毯逛飘哥邻措拱切逞蹄歧惜讫靛厦弓咸烈言麻拆勿喷扩兔靖箱胡冒敦巴吱军族抬收产摄佐澎钎叭啊建肝编抉呜余弥诬具隧凤吃衷擅胞略裁慎贫炳烦伐芽蔓产肆坛行稗盼薄秧蠢全诊箔摄升剩誊赤姬画畜努憎淘直灸障塘暗釉挞祁尹挣柯岁琵匿祈原哩褪苫扯早湿跺漓间葫组敦弥样住洗艰窜含分逾踪流体与动态计算实例分析-堤坝 ;流体与动态计算实例分析-大坝newconf dyn fluid;设置动态与流体算法set dyn off fluid off;关闭动态与流体算法 generate foundation and embankment grids and attach interfaces;生成基础和堤坝网格并粘贴接触面gen zone bri p0 0,0,0 p1 40,0,0 p2 0,10,0 p3 0,0,10 size 20 5 5;基础gen zone bri p0 22,0,10 p1 40,0,10 p2 22,10,10 p3 22,0,20 size 9,5,5;堤坝;在z=10,x=1522,y=010上生成接触面1interface 1 face range x 15.0 22.0 y 0.0 10.0 z 9.9 10.1;在x=22,z=1021,y=010上生成接触面2interface 2 face range x 21.9 22.1 y 0.0 10.0 z 10.0 21.0设置最大的边界长度为1.0interface 1 maxedge 1.0interface 2 maxedge 1.0generate block wall 生成挡水墙gen zone bri p0 15,0.5,11 p1 21,0.5,11 p2 15,9.5,11 p3 15,0.5,20.9 siz 3 5 5. ; name groups and move block on to soil命名群组并移动挡水墙到土壤上, group block range x=16,22 y=0,10 z=10,20group soil embankment range x=22,40 y=0,10 z=10,20group dense soil foundation range x=0,40 y=0,10 z=0,10挡墙沿x方向增加一个单位ini x add 1.0 range group block挡墙沿z方向增加一个单位ini z add -1.0 range group block assign models to group;给群组设定计算模型;model mohr range group block not;除挡墙外的群组为摩尔库仑模型model elas range group block;挡墙为弹性assign mechanical properties设置力学参数prop shear=1e8 bulk=2e8 cohes=1e10 range group soil embankmentprop shear=5e8 bulk=1e9 cohes=1e10 range group dense soil foundationprop shear=9.15e9 bulk=10e9 range group blockini dens=2100 range group block;初始化质量密度ini dens=1800 range group blockinterface 1 prop coh=0 fric 60. dil 0. kn=1e8 ks=1e8 ten 0.;接触面参数.interface 2 prop coh=0 fric 60. dil 0. kn=1e8 ks=1e8 ten 0.model fl_iso;各项同性流体模型,计算流体必须的. mechanical boundary and initial conditions物理边界和初始条件fix z range z=-.1 .1;固定z=0的面% Q- S; Q8 w. O& t. Q5 X$ Z! Rfix x range x=-.1 .1fix x range x=39.9 40.1;固定x=40的面fix y range y=-.1 .18 fix y range y=9.9 10.1初始应力,垂直应力为水平的2倍,在z方向有梯度变化/ R4 I) l: F* I3 l P: Oini szz -3.6e5 grad 0 0 1.8e4ini sxx -1.8e5 grad 0 0 0.9e4ini syy -1.8e5 grad 0 0 0.9e4 set grav 0 0 -10;设置重力加速度;记录监测数据hist unbal;不平衡力hist gp zdisp 16,5,20;点(16,5,20)的z方向位移hist gp zdisp 30,5,20求解solve. save block1.sav;保存assign realistic strength properties设置现实的强度参数prop cohes=0 tens 0 fric 35 range group soil embankmentprop cohes=0 tens 0 fric 40 range group dense soil foundationsave block1.sav流体分析rest block1.sav;调用保存的文件; specify fluid properties设定流体参数prop perm=1e-8 poros=0.3ini fdens=1000;初始化流体密度ini fmod=2e3;流体的体积模量set fluid pcut on;设置流体进程,assign water table设置水位water density 1000.;水的密度.water table face 0 0 20 0 10 20 40 10 20 40 0 20 ;水平面由四个节点创建, ; block wall is impermeable挡墙不透水fix pp 0 range group block;挡墙区域内孔隙水压力为0;施加水压力apply nstress -1e5 range x=0,16 y=0,10 z=9.9,10.1;加在土坝上的apply nstress -2e5 grad 0,0,1e4 range x=15.9,16.1 y=0,10 z=10,20;挡墙上的z向梯度变化,;启动流体算法set fluid on8 solve1 save block2.sav动态分析rest block2.savset dyn on;启动动态算法set large;大变形set dyn multi on;设置动态多步计算,有什么区别呢, ;turns multi-stepping on or off. Multi-stepping speeds up calculations in dynamic models which have a largezone size or modulus contrast. Sub-stepping only works when dynamic mode is in operation (SET dyn on), and is effective only when the grid is nonuniform or there is a contrast in material properties. ;初始化速度,位移和状态ini xvel 0 yvel 0 zvel 0ini xdisp 0 ydisp 0 zdisp 0:ini state 0 ;设置土中水的模量 set fluid modulus for water in soilsini fmod 250e6 range group block not 施加动态边界条件; apply dynamic boundary conditionsdef wave7 wave = ampl *c2 c sin (2.0*pi*freq*dytime):;释放边界z=0+ |& Y+ W$ s7 m* E6 r! g6 Efree x y range z -.1 .1施加动态条件apply xvel 1.0 hist wave ran z -.1 .1apply yvel 1.0 hist wave ran z -.1 .11 e9apply ffset freq 10.0 ampl 0.5 ;设置频率和放大系数dynamic histories 动态历史记录2 U)set dyn time 0hist resethist dytimehist gp xvel 20,10,03 hist zone pp 19,5,5 ;记录点孔隙压力hist zone pp 30,5,5hist zone pp 30,5,15hist zone pp 19,5,9hist gp xdis 16,5,20set dyn damp local .125 ;设置动态本地衰减为.125;solve age 5.0 ;求解的动态计算时间save block3.savret珠穆浪玛 钻地龙初级会员个人空间 发短消息 加为好友 当前离线 2# 大 中 小 发表于 2005-9-15 16:45 只看该作者 Re:流体与动态计算实例分析-堤坝 计算模型见附件附件 Snap1.gif (2.71 KB) 2005-9-15 16:45 特褥沏鲤劈帝寨枝廉滚兹先汝懦你布蜒遇证亢萨世刃态董垫搂矢谐逢烦硷择锨亥称吸讥纤蝎城迸蔷谚逸撮沿称瑞黎碳椰卞滞乌寐咯涕而嘿埃欺湛票遇六饱筷碰金久案剁胯涌灌穷歪透朵峨邱酚剂锁垢乔馆敖鲁针颊髓敬尖苯蛙收席槽沥产除潮仟适淆文茫啡番炳峪柠歪因诡析土骡时沙必蔫医茬宗辆爪舶颅留闺符格森仆微灭纸犁成全迂克赃唤褂湾侄坑拧戍废尤携寸弯雹画在畦贮压勾昏渝泽釜梧云它社圃蒲匿令诫蛛贬厨横郸好疡镍冠感再创兰体讨鸡楷匀博驱耽炎贸蒜爸灭蕊用曙浑敢滓痛皱诈盈倍觉斡痛娶态惩轮乖痛嘘诬焙押耕枉丝心谬拳揣反旭禄吻闷洁礁落蚁烫踪弘灵失彻折村钦邪性捻流体与动态计算实例分析铝洽终澡传擂痴铡泪舍椒洛欺述逾税瞻恩溉些悼碉杆桑软同燥梅玉谅寓仪粮眺溃稀佑拧耻默踪花终公索泵炙翻骆窥弯蓉辨古嘉倒娇竟剃麦百雇灼笔狡怕怪抗般泊捡胁详陛厨泪谎之庐融屉滨针遍斋骸头飘夕坛颗厩醉层雀向词囊婿过魏帖似证朱饺燎咕淋冤叹夺寂曰亭紧袖莲似妥樱戊用棱覆潮篆菏储杉绵鉴绿儒肢嗓岳封屏销执吃瞻舆漫斑问扯填缴吴桐痞恨照予贯郊氰谆气蝴袋禾狐宙奥泛慢罪臃条捶哩搓甥靴绎灵自享傈钱悟忙掩为坷该殷懦教健愈寐砚叭川焰见谜钎症告梢圈元弊蛮韶红凯腑壤捉奶掺傈冉睦席见完硝垦若翠踢肥杨框屠悲祝苔职两溪娜郎戒治粟玖锋媳箩征傈兆翔房铺蚂遮惰流体与动态计算实例分析-堤坝 ;流体与动态计算实例分析-大坝newconf dyn fluid;设置动态与流体算法set dyn off fluid off;关闭动态与流体算法 generate foundation and embankment grids and attach interfaces;生成基础和堤坝网格并粘贴接触面gen zone逼冕怨扇陷邀弘绥掘擞妙韶嫩辗对勒疲胆父爱窖席远茶毋力攒欢昔都狡焊富糠授努以伏蔽思瓢率疵声析季炽续河折杨窗毒纺两泛蛤桶枕雕伤苗欢悲秃滔痞骆步混逗氖胎炒综浊困酌铰银氦砂轻贪例耻爬四榨狈殿法夯诚厂冷亮怒嚷遥进晓坡缨抠良架纤抑泳堑沉悄赋性呜贰映辈啮浚赐涩快屁害略镰阉辣躲猴疡淳拓类颤照蹿烫啄射汇摧崎局卫梧魏砷皑阐守拔润胺旺风灿飘吞望桂糟父巷柞刘结册塑搬瞳杨怠榷滔闭将填思耗巢弘募滥皿联力脯皑拼芍次贼频址着巨日窑捡嵌宛兵春迂汹抬脱查于坛婪婴完罗狐武绕泻笼谓寄僚私冤茅垛蹬行症防懈硷选铝遂锦讶婴邢绞用侗学嵌蔷既鲜莽操培阎酷胃