某大学FLACFLACD在岩土工程中的应用讲义

FLAC/FLAC3D在岩土工程中的应用主要内容FLAC3D软件简介基本原理前后处理FISH语言简介流-固耦合分析初始应力的生成接触单元与应用完全非线性的动力分析自定义本构模型结构单元及应用FLAC(2D)GIIC操作接触面设置流体计算FLAC3D软件简介基本原理前后处理FISH语言简介流-固耦合分析初始应力的生成接触单元与应用完全非线性的动力分析自定义本构模型的基本方法结构单元及应用软件介绍FastLagrangianAnalysisofContinua美国Itasca咨询公司开发2D程序(1986)1990年代初引入中国有限差分法(FDM)DOS版→2.0→2.1→3.0→3.1Itasca其他软件软件介绍应用：岩土力学分析，例矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等岩土工程、采矿工程、水利工程、地质工程特色：大应变模拟完全动态运动方程使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍动力分析功能——地震工程衬砌功能——地下工程可开发功能FLAC3D软件简介基本原理前后处理FISH语言简介流-固耦合分析初始应力的生成接触单元与应用完全非线性的动力分析自定义本构模型的基本方法结构单元及应用有限差分法古老的方法(上世纪40年代)用差分格式转化控制方程中的微商格式流体力学；土工渗流问题；固结FDM&FEM的混合求解FDM的新进展Lagrangian网格源自流体力学中的拉格朗日法跟踪流体质点的运动状态跟踪固体力学中结点，按时步用Lagrangian法研究网格节点的运动节点和单元随材料移动，边界和接触面与单元的边缘一致固体力学大变形理论法国数学家、物理学家拉格朗日

空间混合离散技术结构域离散为可由四面体单元组合形成的五面体或六面体等单元以为基本单元(常应力、常应变)体应变的计算：偏应变的计算：空间混合离散技术+/2=Lagrangian格式动量平衡方程F(t)m牛顿运动定律对于连续体在静力平衡条件下，加速度项为0，方程变为平衡方程Case-1自由落体的模拟G=mgS=1/2gt2=20m命令流：configdyngenzonbrisize111inixmul0.1ym0.1zm0.1modelelaspropbulk3e8shear1e8inidens1000setgrav00-10solveage2Case-1自由落体的模拟(movie)FLAC3D的求解过程平衡方程(动量方程)应力—应变关系(本构模型)Gauss定律单元积分应变率速度节点力新的应力对所有的网格节点对所有单元FLAC3D中的本构模型FLAC3D软件简介基本原理前后处理FISH语言简介流-固耦合分析初始应力的生成接触单元与应用完全非线性的动力分析自定义本构模型的基本方法结构单元及应用FLAC3D的前后处理命令驱动(推荐)程序控制图形界面接口计算模型输出指定本构模型及参数指定初始条件及边界条件，指定结构单元指定接触面指定自定义变量及函数(FISH)求解过程的变量跟踪进行求解模型输出菜单驱动(计算模式)命令栏菜单驱动(Plot)Case-2一个最简单的例子genzonbrisize333;建立网格modelelas;材料参数propbulk3e8shear1e8inidens2000;初始条件fixzranz-.1.1;边界条件fixxranx-.1.1fixxranx2.93.1fixyrany-.1.1fixyrany2.93.1setgrav00-10solve;求解appnstr-10e4ranz3x12y12solveRUNFLAC3D前后处理鞠功能的优船点多种zone类型后处理快塘捷、方便横、丰富计算过陕程中的his床t变量动态音显示FISH可进行参配数化模型猾设计单元状宗态的可敢编程计算暂违停时的芦后处理暂与可保防存FLAC3D软件简宜介基本原理前后处理FISH语言简介流-固耦合推分析初始应闻力的生郑成接触单拥元与应偿用完全非线来性的动力岩分析自定义本疲构模型的毅基本方法结构单元瞒及应用FISH语言简伪介软件自带写的编程语跪言语法简单xxx……end_令xxx注意事在项与FLA任C本身的关骗键字冲突保留字不池可缩写变量可渠不定义珠，因此义注意检单查程序prin扔tfi尘shtabl某e,e痰xtra等命令偶使用Case枣-3F芝ISH与建模两个圆形必隧道的连访接部分变直径蝇的隧道背部分主要内容软件简介基本原仪理前后处垄理FISH语言简脂介流-固耦合分餐析接触单元乐与应用完全非线仅性的动力励分析自定义竞本构模顷型的基怎本方法结构单元蹦及应用FLAC3D流-固耦合输分析(单相流)基本功贝能理论框份架计算模式渗流边界扇条件，初共始条件单渗流计次算及渗流劳耦合计算基本功能渗流各向贝同性、各瞒向异性不同的渗香流模型和匹属性流体压泡力，涌屈入量，护渗漏量惕和不渗批水边界抽水井给、点源厅、体积竟源饱和渗流拖可采用显占式差分法熟、隐式差积分法非饱和渗定流采用显式衰差分法渗流-固体-热的耦合流体和妥固体的抖耦合程暖度依赖袜于土体墓颗粒(骨架)的压缩程俩度，用Biot系数表详示颗粒潜的可压咬缩程度近。循环荷腹载引起弱的动水样压力变榜化和土体液化。地下水模伯拟方法有效应力喝计算不耦合孔压为了暗正确计算器有效应力渗流计薪算已得虑到孔压础分布饱和有水面线疫的部分饱挎和力学变形听产生孔压静力或匀者动力不排水指孔压或遮液化流固耦茫合计算僻模式时间比张例有效应力棵计算不设置CONF仇IGF计luid孔压不绿改变设置孔压防分布INI屈TIA辈Lp迅pWAT嘉ER姑tab宋leWATE辩Rde丹nsit本yWAT康ER搬tab辽le体fac风eSET枯gr蔑avi锯ty手动设置干湿密度设置CON酸FIG现fl骗uid设置土体干密度渗流模型MOD洪EL俘fl_汤iso书tro傅picMOD窜EL快fl_拒ani叼sot伴rop搅icMOD奏EL涨fl_哄nul驶lSET蝇fl沟uid游of连fset泪WATE留Rbu率lk=消0无渗流浩模式渗流模式不耦合骑计算CON俭FIG秒fl炮uidSETmec棒ho瞧ff正确的渗目透系数得到孔压拾分布和水净面线稳态渗流冒可以减小KfCONF染IGf味luidSETflu掩id恒off正确的流商体模量Kf不需要变渗透系仇数单渗流盆模式单力学猫模式流固耦合雨计算CON续FIG昆fl准uidSET借fl图uid昏on真实的撕流体模立量Kf和渗透香系数耦合方卡式DpDevDsDevDp计算模式手动调整嚼的STEP求解主从进缩慧程的SOLV但E求解自动STEP求解流固耦合渡的计算方合法手动调辅整的STE根P求解SET湖fl删uid至on遵me屑ch蜜offSTEPSET拢fl伪uid工of疼fm盒ech青onSTEP主从进程赌的SOL窗VE求解SET拘me休ch创for督ceSET判mech洗sub狮step香na晴uto辣(从进程)SET神mech已sub工step妥m碰(主进程)SOL喜VE侮age自动STEP求解STE瓣P渗流问题(CON溉FIGflu益id)分析步骤时间比欺例(ts,tc)稳态不排水状同态相当扰动类型力学扰归动孔压扰起动流固刚度香比Rk是否>>>骆1完全耦合务模式时间比统例相当宏；力学播扰动渗流边界坊条件，初纵始条件默认的处边界条仿件是不团透水边伶界孔隙压饭力自由(不透水洪边界)固定孔役隙水压度力(透水边界)如：井孔隙压力盒，孔隙率劲，饱和度胆和流体属迟性的初始鹊分布可以相用INIT抓IAL命令或者PROP错ERTY命令定株义。Cas饱e-4真空预宰压的简铸单模拟孔压边界畜条件ts>>tc长期分垮析(排水)Rk>>1骨架很软孔压扰动进行bio邻t_m洋od调整砂层软土层粘土层PVD2m8m10mDat摔af吹ile框:数值分猾析过程(mov汗ie)初始应浮力的生率成为什么扭要单独济列出？分析过程新中出现的扇很多问题熊都与初始卫应力是否惨合理有关手册中的使例子五花丘八门生成方盒法弹性求箩解更改强度拨参数的弹精塑性求解设置初川始应力暴的弹塑桶性求解存在水压撞力的初始从应力生成水下建筑门的初始应吐力生成弹性求旦解模型尺寸单元数量密度KGU1×1×2(m3)1×1×2200030MPa10MPa0.35gen稳zo年nb应ri拜siz罢e1摩1们2me批lasprop碗bul查k3e炸7sh扰ear猴1e7fix遮z餐ran粱z地0fix牵x鼓ran栋x挠0fix鞠x肢ran德x茎1fix都yra索ny摇0fix瓦y岔ran般y捆1ini秃de命ns逮200白0set蛋gr失av简00法-1段0solv兔eStep午=1员62sz=京-40e丑3sx=临-2选1.5须4e3更改强饶度参数着的弹塑泉性求解模型尺寸单元数量密度KGcfu1×1×2(m3)1×1×2200030MPa10MPa10kPa150.35gen久zon丛bri给size痛11姓2mode环lmo桶hrprop杆bul相k3e厨7sh精ear爸1e7c1授e10f1丹5ten较1e10fix则z山ran退z群0fix牛xra五nx亭0fix笼x浮ran饲x辽1fix匹y撕ran剃y丛0fix驱yra要ny嗽1ini会dens视200脖0set捆grav柿00抽-10solv筋eprop掘bul富k3e炭7sh堤ear步1e7c10蜡e3f1傍5ten裤0sol谷veSte乒p=弄16邻3sz=治-4蛙0e3sx=大-2甚1.5触4e3Or:或so厌lve国el碧ast约ic设置初始沈应力的弹瘦塑性求解模型尺寸单元数量密度KGcfu1×1×2(m3)1×1×2200030MPa10MPa10kPa150.35gen环zon默bri择size企11扁2mode煤lmo井hrprop严bul潮k3e句7sh截ear筒1e7c10顿e3f藏15误ten在0fix铲zra奸nz寇0fix游x标ran都x邪0fix抹x瘦ran粉x辅1fix滋y拿ran郊y做0fix贷y肠ran蚂y符1ini焦de所ns颠200江0ini衔szz谣-40e辞3gr睡ad0滩02皂0e3刮ran乳z0香2inisyy-20e奔3gr聚ad0台01啊0e3喉ran盒z0午2inisxx-20迷e3镰gra才d0嘉0悠10e呆3r姻an衰z0千2set巡寿grav呜00浪-10solv推eStep偷=0sz=家-40e穴3sx=-2直0e3存在水办压力的丹初始应邮力生成模型尺寸单元数量饱和密度KGcfu水位线孔隙率1×1×2(m3)1×1×2200030MPa10MPa10kPa150.351m0.5gen固zon枯bri诉size杜11秒2mod跑el浩mprop粪bul着k3e浪7sh个ear访1e7c10定e10兄f15载ten移1e1半0fix榆z耐ran新z微0fix叔x劲ran跪x架0fix颈xra莲nx姿1fix衫yra他ny堆0fix秒y泪ran踢y圈1ini蹄de旱ns饼200间0r所an羽z0忽1ini嫂de筐ns贯150行0r芹an随z1猴2ini艘sz泼z-找35e欣3院gr盛ad忆00劈燕20仿e3挨ra球nz读0览1ini窑syy义-17.舍5e3卵grad供00辛10e搁3r晕anz劫01ini托sxx晶-17.谢5e3肥grad辣00督10e胃3r冶anz岂01ini扎szz处-15e烦3什grad茅00狗15e辅3r酒anz赌12ini删sy陵y-养7.5冈e3砍gr找ad肌00茫7.关5e3鬼ra久nz百1轨2ini雕sxx序-7.5妻e3圈grad风00萌7.5还e3r堂anz情12ini霉pp1批0e3瘦grad省00鹊-10旨e3r辆anz摘01set材grav炕00逝-10Ste矿p=83sz=亚-33.治75e3sx=-韵23.亿21e矿3rd=rs–尼n×际s×rf水下建前筑的初轧始应力括生成模型尺寸单元数量饱和密度KGcfu水位线1×1×2(m3)1×1×2200030MPa10MPa10kPa150.353mgen棵zon辟bri喝size析11馅2mode华lmpro爆pb拳ulk衫3e己7s霞hea鄙r1菌e7席c1烂0e1扶0f起15溉te遭n1诱e10fix垂z贿ran淘z吐0fix承xra矛nx倘0fix厦x浑ran逐x缓1fix贝y狭ran吼y每0fix扁y辩ran爹y盐1ini请dens议200翻0ra模nz王02ini雀sz紫z-雀50e箭3g材rad盛0瞎02鸡0e3宋r罩an宗z0妖1ini崭syy撕-30e希3gr稿ad0曲01续0e3赤ran未z0羞1ini购sxx-30e查3grad咸00采10e仙3r说anz照01ini壮pp3编0e3辟grad猛00粪-10幕e3r偶anz曾02app错nstr恒ess吊-10e吐3ra碍nz厦2set思grav佛00凤-10solv捐eStep唐=111sz=纳-5钓0e3sx=谋-44.赖4e3初始应力皂检查办法初始应威力的计稳算时间供不会“骑非常长拣”经常检查盘模型的响馒应plo悼co幅ns衫zz谅(sy倡y,载sxx乱)应力场plo维con晃zdis茄(yd税is,俱xdis进)位移场plo辈blo霜sta屈服状态plo若gp符fix盘re俩ds瘦k速度约滨束条件plo零fap贺red样sk体力plo箱hi自st富(un肾bal削)不平衡杀力plo赖inte效rfac戴ens丸tres饭s(s彼stre拐ss)接触面单糟元主要内容软件简介基本原惭理前后处理FISH语言简怀介流-固耦合分吹析初始应计力的生污成接触单胜元与应宰用完全非线贱性的动力抱分析自定义本逝构模型的地基本方法结构单代元及应网用接触面单酸元的用途岩体介质尸中的解理庭、断层、肥岩层面地基与土漏体的接触箱、槽及呀其内充填微物的接触空间中吩无变形香的固定讽“障碍注”接触面倡的原理三角形单絮元(无厚度!)参数较布多三种工刷作模式粘结界竹面粘接滑移库伦滑动接触单元太模型的建爱立(1)关键要形成同攻一位置的皇两个节点(面)“移来移题去”(推荐)建两个分攀开的模型建立接老触单元通过INI*add使模型冤接触注意dist的含义NOm键erge钥!!!接触面dis蛋t1234接触单妨元模型湾的建立(2)“导来导去竭”利用exp爆gri脑d,逐imp上gri寺d命令进委行网格笨导出与甩导入配合DELE毁TE命令适于内部剖接触面的怜建立，或其他前逐处理工磨具建立慈的网格“导来导去创”具体方或法save姐1.s跃avdel储ran计grop锄2n每otInte腔rfac耀e1表facesav登e2驰.sa延vres且t1吊.sa墙v

d趟el奔ran兔gr威oup恭2弦exp途gri郊d1磁.fa顷c3dres驶t2澡.sa泻v

i猴mpg盈rid胁1.航fla庸c3d接触面战参数的苦确定虚构的贞为了合少并节点瘦而设置椅的接触汽面Kn=k拍s=10木*真实的刚婚性接触面如料仓斩下料c,D杏,Te文nsi蹄on重要，kn,k哨s不重要真实的祝柔性接蔑触面断层；水力劈俯裂材料试验得留到参数对于kn,笋ks：岩石虾断层10~商100议MPa肉/m(粘土);1面00GP姐a(岩石)反分析方宏法：通过棚断层中疼岩石的浊变形与窃原岩的识变形主要内攀容软件简武介基本原迅理前后处糠理FIS汗H语言简贡介流-固耦合分诸析初始应森力的生吩成接触单元株与应用完全非线涛性的动力疤分析自定义本初构模型的筒基本方法结构单臭元及应筛用完全非挪线性的屯动力分石析与等效硬线性方从法的比盆较动力荷话载动力边界杏条件力学阻尼熟于滞回阻狂尼地震波的云调整动孔压的怒生成FLAC与等效线墨性方法等效线衔性方法述是岩土漠地震工族程中的喘常用方婚法动本构模片型等效线赤性模型抗剪强度息的降级曲父线阻尼比询随剪应古变的变甲化FLAC常规模型(MC)多种频率你成分的干凯涉和混合永久变形费计算弹塑性计射算动力荷载动力输脏入的类喊型加速度时识程速度时程应力(压力)时程力时程APPL霜YINT读ERI脸OR饮(内部)TAB婶LEFISHQuie垃t边界静态(qui淹et,粘性)边界Lysm键era去ndK垮uhle李meye纲r(19粘69)模型边界酿法向和切岗向设置独论立的阻尼乎器性能对于法向p波和s波能很糕好的吸膛收对于倾斜杯入射的波嚷和Ray眨lei刚gh波也有盘所吸收湖，但存枣在反射人工边浪界仍应货当足够菊远Quie脱t边界应勾用内部振动(如隧道逝中的列套车振动医问题)动力荷载秒直接施加柜在节点上使用Quie枕t边界减浅小人工轮边界上铁的反射不需要FF边界外部荷见载的底近部边界软土地希基上的隔地震荷平载不适鲁合用加有速度或族速度边六界条件使用应凳力条件t=-2Csrvs地震底部滋输入的侧猴向边界扭曲了入茫射波quie替tquie抄tquie也tFre知e-f替iel洗d边界Cund折all守eta冶l.(脆1980乌)自由场均网格与绑主体网艺格的耦穷合粘性峰阻尼器去，自由秤场网格池的不平户衡力施筒加到主挖体网格伟边界上设置条件底部水近平，重返力方向免为z向侧面垂直佛，法向分绸别为x,y向其他边颤界条件税在APP借LY梁ff之前Fre书e-f袋iel睬d边界APP吓LY旦ff将边界上短单元的属架性、条件佩和变量全首部转移ff单元上幅；设置以后条主体网格桶上的改动乒将不会被FF边界所束响应可存在弱任意的嘱本构模省型以及狱流体耦港合(仅竖向)FF边界进行煌小变形计境算，主体昼网格可大跨变形，FF边界上兼的变形央要相对城较小存在att碍ach的边界将粗不能设置FF边界边界上古的Int迫erf观ace将不能连装续动力边界匆设置需在FF边界设叮置之前力学阻尼瑞利(ray威leig句h)阻尼假设阻尼她与质量、蓬刚度的线代性关系参数确肺定简单(等价平均蝇应变=60级%*emax)中心频坚率(共振计算切，地震平地均频率)临界阻尼钥比计算速度马慢，不推茧荐局部(loc双al)阻尼FLA剑C3D的静力分澡析阻尼参数简单适合简单滨情况滞回阻尼(Hy言ste题ret丸ic笔Dam融pin电g)模拟岩土齐介质的动枕模量衰减妹曲线init菜ial闹damp脊hys释tere晃ticnamesig3唐(3参数)sig4她(4参数)Har肾din葬(1参数)(哈丁模庄型)Defa扭ult（2参数)优点直接采芦用模量尼降级曲妹线阻尼比妇不会影盼响时间丸步缺点输出的曲系线会不一遍致滞回阻尼(Hy喜ste哲ret练ic颂Dam其pin逝g)低循环赶应变下洋得到的搞阻尼比久要小于奇试验结浩果，这写会导致淋低级的鲁噪声，纯尤其在步高频情益况下。饼可以在蒜中心频离率上增沿加一个筋小量的Rayl葵eigh阻尼(~0悄.2%刚度比著例)，这样也捎不会降低玩时步；若初始应列力不为0，剪应碎力-剪应变探曲线可折能不匹蔬配。因六此在生哥成初始趟应力时虑就要调配用Hys嗽t阻尼；Hyst阻尼不摇仅会增威加能量科损失，蕉还会导帖致在大拍循环应演变下的裳平均剪缠切模量孤的降低蒙，在输梯入波的耕基频接绿近共振届频率的亦时候，荷由于可次能会导朱致动力梨反应幅袋值的增骨大；Hys泉t阻尼之分前要做级一次弹幸性无阻企尼求解肉，以获青得发生性循环应桂变的最救大水平仗，若循煤环应变丑过大导存致剪切继模量过搜多的降博低，那妙么用Hys蚀t阻尼是有况问题的；即使应变茎较小，使跳用屈服模怀型也会增涉大应变，脏因此若有猛广泛屈服亏的现象，晚则使用屈肯服模型，联不用Hyst阻尼参考了Ita怎sca的中国浊培训资阿料地震波的乒调整基线校固正对于地择震分析宰的加速建度时程判，其积描分得到轧的速度屠和位移袄应归0美国地质鬼调查研究菌所Basi廉cSt恒rong蜘-Mot旱ion敢Acce耗lero控gram境Pro抚cess段ing过Soft省ware做(BA膝P)对网格施删加一个固批定速度从流而使残余驰的位移变杨为0动力荷载循的频率与弦单元尺寸纯的双向调开整高频的输兵入要求单令元尺寸很度小一定的单敌元尺寸对乞应输入的贪最大频率一般进凉行滤波印处理滤掉低能窄量的高频FFT.弹FISOrig睁inSei沿smo卖Sig采nal地震波的冲调整El-Cen美tro波FFT修正后煎的时程修正后FFT5Hz动孔压缎的生成——液化干沙剪应瓜变循环加膊载试验初始加愧载阶段瞧，沙土着通常先劳压实再沉膨胀。价卸载时含，沙土昆遵循与教加载相新似的路掉径，但棉在零应于变时，始有些残戚余体积珠应变存凉在。取殿决于初独始孔隙饼率，这兄可能代敢表纯粹痕的压实假定孔隙田中充满水对于常体挠积测试，侦有效应力妥降低，孔鞠隙水压保吨持不变对于常翼荷载测浮试，(例如，县盒子上株法向荷肝载固定)，孔隙祖水压增零加，有妻效应力售减小有效应力信为零时发无生液化动孔压的嫁生成——液化因此孔六隙水压缺增加不交是液化暂的基本竹原因由于颗粒渴间(重组以后)的低接触跪力导致有撞效应力的溜减小描述液化岂的模型高级模型田：BSH深P(边界面低弯塑性本构贤模型,Wa巨nge朋tal剂.19概90)简单模柔型：MC+体积应变浅增量模型Finn模型：Byr时ne模型：主要内严容软件简魄介基本原理前后处理FIS丹H语言简介流-固耦合向分析初始应力师的生成接触单掘元与应挽用完全非酷线性的嚼动力分茶析自定义辟本构模个型的基音本方法结构单胡元及应级用自定义肯本构模灵型的基屠本方法必要性试验总结疲的本构模猜型特定条姨件下的蛛本构模屑型交叉学尊科的本摇构模型二次开未发环境自定义本侨构模型的相功能自定义悼本构模选型的基室本方法二次开发固环境FLA丸C3D采用面向区对象的语引言标准C++编写本构模型训都是以动辟态连接库稳文件(.D乏LL文件)的形式提菜供VC+雁+6.朵0(S但P4)或更高坚版本的敞开发环誉境优点自定义的越本构模型鼻和软件自慌带的本构触模型的执陶行效率处插在同一个万水平自定义本锋构模型(.D座LL文件)适用于零高版本赵的FLAC突(2D)、3DE摇C、UDEC等其他Itas呢ca软件中自定义本米构模型的颗功能主要功分能：对扛给出的拒应变增胶量得到寇新的应寨力辅助功能画：模型名称川、版本读写操作模型文件王的编写基类(cl莫ass乓Co跟nst夫itu闭tiv扭eM架ode检l)的描述成员函推数的描船述模型的注滋册模型与FLA棵C3D之间的信柳息交换模型状先态指示狂器的描愚述自定义本驰构模型的艇基本方法头文件(use接rmod嘱el.h贩)中进行新恶的本构模枝型派生类炮的声明修改模型图的ID(>螺100)、名称和洽版本修改派生映类的私有黄成员C++文件(us杀erm挨ode拼l.c况pp)中修改品模型结至构(Us池erM纺ode果l::某Use魔rMo尝del蹈(bo衔ol理bRe肝gis猾ter裕):邻Con精sti厕tut丧ive迅Mod脚el)cons继tch罪ar*售*Use煮rMod记el::征Prop迟erti养es()函数模型的参秩数名称字蹄符串con挤st稍cha雀r*愈*Us芝erM雕ode场l::蚊Sta扛tes犁()函数计算过社程中的溪状态指道示器自定义本岸构模型的就基本方法dou校ble拌Us屡erM捧ode拘l::隆Get裤Pro写per泄ty(嫩)和void焦Use丽rMod召el::持Set轻Prop企erty性()函数con皆st袋cha出r*送Us耽erM盆ode地l::姿Ini诵tia带liz挠e()函数参数和旧状态指企示器的惹初始化某，并对颈派生类珍声明中仓定义的压私有变愚量进行阶赋值con滩st耀cha仍r*斤Us央erM炒ode呼l::Run(涉)函数由应变靠增量计让算得到揭应力增雕量，从靠而获得猪新的应忘力con堵st婶cha双r*南Us鞋erM兽ode缸l::蔑Sav狭eRe尺sto诊re(节)函数对计算结范果进行保书存。程序的调醋试在VC+慨+的工程设年置中将FLA碎C3D软件中的EXE文件路价径加入克到程序白的调试命范围中倚，并将FLA屑C3D自带的DLL文件加入止到附加动杀态链接库(Add渐itio救nal竞DLLs村)中，然后剪在Init序iali饺ze()或Run协()函数中修设置断哥点，进秘行调试洋；在程序文猜件中加入retu泼rn()语句，搏这样可救以将希陷望得到棋的变量稳值以错树误提示抚的形式杜在FLA进C3D窗口中找得到。例1:非线性弹车性(Dun另can-症Chan烤g)模型例2:棵FLA夸C液化模蚊型的修蓬正Yu-眨min风CH氧EN,非Ha双n-l凑ong挣LI丹U.昂Cou插ple闷dh敢ydr杜aul德ic-咳mec钟han语ica仁la将nal丢ysi浸so碍fl扛arg饼ed唯efo布rma馋tio封ni叼ndu炕ced适by流po激st-可liq容uef百ied撒sa亦nd.吧Ge龟oPr爆oc2男006园,N仗anj连ing眠,C遭hin页a:贷700围-70葱5BeginInitialliquefaction?LiquefiedcalculationYESPost-liquefiedcalculationLiquefaction?NONormalFinnModelYESNODynamictimefinished?EndNOYES主要内容软件简行介基本原直理前后处理FIS伟H语言简介流-固耦合分落析接触单均元与应抗用完全非歇线性的裁动力分辫析自定义酿本构模矿型的基兔本方法结构单元融及应用FLAC3D中的结构并单元有限单但元梁(bea辆m)单元锚索(cab斗le)单元桩(pi妥le)单元锚杆:r巧ock导bol甲to建n壳(sh携ell膨)单元格栅(ge纺ogr斯id)单元土工织引物；土拉工格栅初衬(li地ner莫)单元beamcablepileshellgeogridlinerCase似-5结构的祸动力响称应con壤fig侧dy抄nsel驼pi梁le粥id=歌1b家eg侮00捆0厚end牛0沈01sel钢pile思pro继pde栋ns2窃400敞&Emo曾d1喷.0e及10轮Nu河0.3效XC丙Are难a0糊.3外&XCJ坚0.16纵375眼XCIy贴0.0酱0625浮XCI掩z0.课0157必5&Per储2.创8化CS_迷sK云1.3企e11宗CS脸_nK逼1.蹦3e1话1&CS_n搬Gap腰offsel主node菌fix减xy愤zx冠ryr馅zr覆ran乞id=1sel些se秆td蔬amp松co似mbi滚neddeff1whil帜este仆ppin剧gf0=恢100办00*肢sin袜(10耽*dy豆tim架e)np=植nd_汤headloop冈whi牙len攀p#踏nullifn板d_po客s(np螺,1,3渔)=1nd_a跟pply总(np,浇1)=f各0end梢ifnp呢=n移d_n究ext丢(np诵)endl坏oopendsolv叶eag纹e1pile1000看0*si跌n(10针*t)Case瓣-4结构的压动力响旗应结构单锤元的应浙用土与结珠构的相涉互作