FLAC3D基础介绍-FLAC3D在岩土工程中的应用.ppt

《FLAC3D基础介绍-FLAC3D在岩土工程中的应用.ppt》由会员分享，可在线阅读，更多相关《FLAC3D基础介绍-FLAC3D在岩土工程中的应用.ppt（73页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、GeoHohaiGeoHohai主要内容主要内容nFLAC3D软件简介软件简介nFLAC3D的基本原理的基本原理nFLAC3D的前后处理的前后处理n流流-固耦合分析固耦合分析n接触单元与应用接触单元与应用n完全非线性的动力分析完全非线性的动力分析n自定义本构模型的基本方法自定义本构模型的基本方法n结构单元及应用结构单元及应用1GeoHohaiGeoHohai主要内容主要内容nFLAC3D软件简介软件简介nFLAC3D的基本原理的基本原理nFLAC3D的前后处理的前后处理n流流-固耦合分析固耦合分析n接触单元与应用接触单元与应用n完全非线性的动力分析完全非线性的动力分析n自定义本构模型的基本方法

2、自定义本构模型的基本方法n结构单元及应用结构单元及应用2GeoHohaiGeoHohaiFLAC3D简介简介nFast Lagrangian Analysis of Continuan美国美国Itasca咨询公司开发咨询公司开发2D程序程序(1986)n1990年代初引入中国年代初引入中国n有限差分法有限差分法(FDM)nDOS版版2.0 2.1 3.0nItasca其他软件其他软件3GeoHohaiGeoHohaiFLAC3D简介简介n应用：应用：岩土力学分析，例矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等岩土工程、采矿工程、水利工程、地质工程 n特色：特色：大应

3、变模拟完全动态运动方程使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍 显示求解具有较快的非线性求解速度4GeoHohaiGeoHohai主要内容主要内容nFLAC3D软件简介软件简介nFLAC3D的基本原理的基本原理nFLAC3D的前后处理的前后处理n流流-固耦合分析固耦合分析n接触单元与应用接触单元与应用n完全非线性的动力分析完全非线性的动力分析n自定义本构模型的基本方法自定义本构模型的基本方法n结构单元及应用结构单元及应用5GeoHohaiGeoHohai基本原理基本原理n有限差分法有限差分法nLagrangian网格网格n空间混合离散技术空间混合离散技术nLagrangian

4、格式动量平衡方程格式动量平衡方程nFLAC3D的求解过程的求解过程nFLAC3D的本构模型的本构模型6GeoHohaiGeoHohai有限差分法有限差分法n古老的方法古老的方法(上世纪上世纪40年代年代)n用差分格式转化控制方程中的微商格式用差分格式转化控制方程中的微商格式n流体力学；土工渗流问题；固结流体力学；土工渗流问题；固结nFDM&FEM的混合求解的混合求解nFDM的新进展的新进展7GeoHohaiGeoHohaiLagrangian网格网格n源自流体力学中的拉格朗日法源自流体力学中的拉格朗日法跟踪流体质点的运动状态跟踪固体力学中结点，按时步用Lagrangian法研究网格节点的运动n

5、节点和单元随材料移动，边界和接节点和单元随材料移动，边界和接触面与单元的边缘一致触面与单元的边缘一致n固体力学大变形理论固体力学大变形理论法国数学家、物理学家拉格朗日 8GeoHohaiGeoHohai空间混合离散技术空间混合离散技术n结构域离散为可由四面体单元组合形成的五面体结构域离散为可由四面体单元组合形成的五面体或六面体等单元或六面体等单元n以以 为基本单元为基本单元(常应力、常应变常应力、常应变)n体应变的计算：体应变的计算：n偏应变的计算：偏应变的计算：9GeoHohaiGeoHohai空间混合离散技术空间混合离散技术+/2=10GeoHohaiGeoHohaiLagrangian格

6、式动量平衡方程格式动量平衡方程F(t)m牛顿运动定律对于连续体在静力平衡条件下，加速度项为0，方程变为平衡方程11GeoHohaiGeoHohaiCase-1自由落体的模拟自由落体的模拟G=mgS=1/2gt2=20m命令流：命令流：config dyngen zon bri size 1 1 1ini x mul 0.1 y m 0.1 z m 0.1model elasprop bulk 3e8 shear 1e8ini dens 1000set grav 0 0-10solve age 212GeoHohaiGeoHohaiCase-1自由落体的模拟自由落体的模拟(movie)13Geo

7、HohaiGeoHohaiFLAC3D的求解过程的求解过程平衡方程(动量方程)应力应变关系(本构模型)Gauss定律单元积分应变率速度节点力新的应力对所有的网格节点对所有单元14GeoHohaiGeoHohaiFLAC3D中的本构模型中的本构模型n开挖模型开挖模型null n3个弹性模型个弹性模型各向同性弹性横观各向同性弹性正交各向同性弹性 n8个塑性模型（个塑性模型（Drucker-Prager模型、模型、Morh-Coulomb模型、应变硬化模型、应变硬化/软化模型、遍布节理软化模型、遍布节理模型、双线性应变硬化模型、双线性应变硬化/软化遍布节理模型、修正软化遍布节理模型、修正剑桥模型和剑

8、桥模型和胡克布朗模型胡克布朗模型）15GeoHohaiGeoHohaiFLAC3D中的本构模型中的本构模型16GeoHohaiGeoHohai主要内容主要内容nFLAC3D软件简介软件简介nFLAC3D的基本原理的基本原理nFLAC3D的前后处理的前后处理n流流-固耦合分析固耦合分析n接触单元与应用接触单元与应用n完全非线性的动力分析完全非线性的动力分析n自定义本构模型的基本方法自定义本构模型的基本方法n结构单元及应用结构单元及应用17GeoHohaiGeoHohaiFLAC3D的的前后处理前后处理n命令驱动命令驱动(推荐推荐)程序控制图形界面接口计算模型输出指定本构模型及参数指定初始条件及边

9、界条件，指定结构单元指定接触面指定自定义变量及函数(FISH)求解过程的变量跟踪进行求解模型输出18GeoHohaiGeoHohai菜单驱动菜单驱动(计算模式计算模式)命令栏19GeoHohaiGeoHohai菜单驱动菜单驱动(Plot)20GeoHohaiGeoHohaiCase-2 一个最简单的例子一个最简单的例子gen zon bri size 3 3 3 ;建立网格建立网格model elas ;材料参数材料参数prop bulk 3e8 shear 1e8ini dens 2000 ;初始条件初始条件fix z ran z-.1.1 ;边界条件边界条件fix x ran x-.1.1

10、fix x ran x 2.9 3.1fix y ran y-.1.1fix y ran y 2.9 3.1set grav 0 0-10solve ;求解求解app nstr-10e4 ran z 3 x 1 2 y 1 2solveRUN FLAC3D21GeoHohaiGeoHohai前后处理功能的优点前后处理功能的优点n多种多种zone类型类型n后处理快捷、方便、丰富后处理快捷、方便、丰富n计算过程中的计算过程中的hist变量动态显示变量动态显示nFISH可进行参数化模型设计可进行参数化模型设计n单元状态的可编程单元状态的可编程n计算暂停时的后处理与可保存计算暂停时的后处理与可保存22

11、GeoHohaiGeoHohai前后处理功能的缺点前后处理功能的缺点n复杂模型的建模功能不强复杂模型的建模功能不强可以编程导入其他软件形成的网格(比如：Ansys、Adina、GeoCAD)n无等值线的后处理功能无等值线的后处理功能(3D)可编程将.sav文件写入TecPlot等其他后处理软件n全命令操作，学习困难全命令操作，学习困难n鼠标功能单一鼠标功能单一(双击取击点坐标双击取击点坐标)23GeoHohaiGeoHohai主要内容主要内容nFLAC3D软件简介软件简介nFLAC3D的基本原理的基本原理nFLAC3D的前后处理的前后处理n流流-固耦合分析固耦合分析n接触单元与应用接触单元与应

12、用n完全非线性的动力分析完全非线性的动力分析n自定义本构模型的基本方法自定义本构模型的基本方法n结构单元及应用结构单元及应用24GeoHohaiGeoHohai流流-固耦合分析固耦合分析(单相流单相流)n基本功能基本功能n理论框架理论框架n计算模式计算模式n渗流边界条件，初始条件渗流边界条件，初始条件n单渗流计算及渗流耦合计算单渗流计算及渗流耦合计算 25GeoHohaiGeoHohai基本功能基本功能n渗流各向同性、各向异性渗流各向同性、各向异性n不同的渗流模型和属性不同的渗流模型和属性n流体流体压压力，涌入量，渗漏量和不渗水力，涌入量，渗漏量和不渗水边边界界n抽水井、点源、体抽水井、点源、

13、体积积源源n饱饱和渗流可采用和渗流可采用显显式差分法、式差分法、隐隐式差分法式差分法 n非非饱饱和渗流和渗流采用采用显显式差分法式差分法n渗流渗流-固体固体-热热的耦合的耦合n流流体体和和固固体体的的耦耦合合程程度度依依赖赖于于土土体体颗颗粒粒(骨骨架架)的的压压缩缩程度，用程度，用Biot系数表示系数表示颗颗粒的可粒的可压缩压缩程度。程度。n循环荷载引起的动水压力变化和循环荷载引起的动水压力变化和土体液化土体液化。26GeoHohaiGeoHohai理论框架理论框架n准静态准静态Biot理论理论 n多孔介质中遵循多孔介质中遵循Darcy定律的单相渗流定律的单相渗流n描述多孔介质中流体渗流的变

14、量描述多孔介质中流体渗流的变量孔隙水压力，饱和度，特定排水向量的三个分量 n质量守恒定律质量守恒定律n达西定律达西定律n本构定律本构定律考虑流体响应孔隙水压力改变，饱和度改变，体积应变改变和温度改变 27GeoHohaiGeoHohai流流-固耦合的计算模式固耦合的计算模式n无渗流模式无渗流模式孔压计算n有渗流模式有渗流模式瞬态渗流分析流-固耦合计算28GeoHohaiGeoHohai无渗流模式无渗流模式n不设置不设置CONFIG Fluidn孔压不改变孔压不改变nINITIAL ppnWATER tableWATER densitySET gravityWATER table facen手动

15、设置手动设置干湿密度干湿密度n设置设置CONFIG fluidn瞬态渗流分析瞬态渗流分析n有效应力计算有效应力计算n不排水计算不排水计算 n设置土体设置土体干密度干密度 n渗流模型渗流模型MODEL fl_isotropic MODEL fl_anisotropic MODEL fl_null 渗流模式渗流模式29GeoHohaiGeoHohai渗流边界条件，初始条件渗流边界条件，初始条件n默认的边界条件是不透水边界默认的边界条件是不透水边界 n孔隙压力自由孔隙压力自由(不透水边界不透水边界)n固定孔隙水压力固定孔隙水压力(透水边界透水边界)如：井 n孔隙压力，孔隙率，饱和度和流体属性的初始分

16、孔隙压力，孔隙率，饱和度和流体属性的初始分布可以用布可以用INITIAL命令或者命令或者PROPERTY命令定义。命令定义。30GeoHohaiGeoHohai单渗流计算及渗流耦合计算单渗流计算及渗流耦合计算 n时间比例时间比例 n完全耦合分析方法完全耦合分析方法 n孔压固定分析孔压固定分析(有效应力分析有效应力分析)n单渗流得到孔压分布单渗流得到孔压分布 n无渗流计算无渗流计算孔压的力学响应孔压的力学响应 n流流-固耦合计算固耦合计算 31GeoHohaiGeoHohai时间比例时间比例(scale)n力学过程的特征时间力学过程的特征时间 n流体扩散过程的特征时间流体扩散过程的特征时间 32

17、GeoHohaiGeoHohai完全耦合分析方法完全耦合分析方法n时间比例时间比例短期行为(不排水)n ts(分析时间)tcn施加扰动的属性施加扰动的属性流体扰动：渗流可不与力学过程耦合力学扰动：耦合等级取决于流固刚度比n流固刚度比流固刚度比33GeoHohaiGeoHohai单渗流得到孔压分布单渗流得到孔压分布n用途：排水沟；抽水井；耦合计算用途：排水沟；抽水井；耦合计算n计算步骤计算步骤CONFIG fluid SET mech offSET fluid implicit on/offMODEL fl_;PROP STEP;SOLVE age;SET fluid ratioSET flui

18、d off mech onnPROP biot_c 0(or INI fmod 0)34GeoHohaiGeoHohai无渗流计算无渗流计算孔压的力学响应孔压的力学响应n不排水短期响应不排水短期响应n两种分析方法：干法和湿法两种分析方法：干法和湿法干法：Ku=K+a2Mn两种破坏形式WATER或INI获得常孔压，不排水的c,(孔压改变较小)=0,c=cu(MK+4/3G)湿法：耦合体系的短期行为n使用排水的K,c,n若SET fluid off,Biot_mod(fmod)真实35GeoHohaiGeoHohai流流-固耦合计算固耦合计算nCONFIG fluid;M(Kf);K(渗透系数渗透

19、系数)真实，则真实，则FLAC3D默默认耦合计算认耦合计算pvvpn预估流预估流/力特征时间力特征时间n耦合计算前先达到一个平衡状态耦合计算前先达到一个平衡状态SET fluid on mech off;SET fluid off mech on;STEPSET mech force;SET mech substep n auto;SET fluid substep m(=1)STEP:渗流步足够小36GeoHohaiGeoHohaiCase-3真空预压的简单模拟真空预压的简单模拟n孔压边界条件孔压边界条件ntstc长期分析(排水)nRk1骨架很软n孔压扰动孔压扰动进行biot_mod调整砂层

20、软土层粘土层PVD2m8m10mData file:Data file:37GeoHohaiGeoHohai数值分析过程数值分析过程(movie)38GeoHohaiGeoHohai主要内容主要内容nFLAC3D软件简介软件简介nFLAC3D的基本原理的基本原理nFLAC3D的前后处理的前后处理n流流-固耦合分析固耦合分析n接触单元与应用接触单元与应用n完全非线性的动力分析完全非线性的动力分析n自定义本构模型的基本方法自定义本构模型的基本方法n结构单元及应用结构单元及应用39GeoHohaiGeoHohai接触面单元的用途接触面单元的用途n岩体介质中的解理、断层、岩层面岩体介质中的解理、断层、

21、岩层面n地基与土体的接触地基与土体的接触n箱、槽及其内充填物的接触箱、槽及其内充填物的接触n空间中无变形的固定空间中无变形的固定“障碍障碍”40GeoHohaiGeoHohai接触面的原理接触面的原理n三角形单元三角形单元(无厚度无厚度!)n参数较多参数较多n三种工作模式三种工作模式粘结界面粘接滑移库伦滑动41GeoHohaiGeoHohai接触单元模型的建立接触单元模型的建立(1)n关键关键要形成同一位置的两个节要形成同一位置的两个节点点(面面)n“移来移去移来移去”(推荐推荐)建两个分开的模型建立接触单元通过INI*add使模型接触注意dist的含义接触面dist42GeoHohaiGeo

22、Hohai接触单元模型的建立接触单元模型的建立(2)n“导来导去导来导去”利用expgrid,impgrid命令进行网格导出与导入配合DELETE命令适于内部接触面的建立，或其他前处理工具建立的网格43GeoHohaiGeoHohai“导来导去导来导去”具体方法具体方法nsave 1.savndel ran grop 2 notnInterface 1 facensave 2.savnrest 1.savdel ran group 2expgrid 1.fac3dnrest 2.savimpgrid 1.flac3d44GeoHohaiGeoHohai接触面参数的确定接触面参数的确定n虚构的为

23、了合并节点而设置的接触面虚构的为了合并节点而设置的接触面Kn=ks=10*n真实的刚性接触面真实的刚性接触面如料仓下料c,D,Tension重要，kn,ks不重要n真实的柔性接触面真实的柔性接触面断层；水力劈裂材料试验得到参数对于kn,ks：岩石断层10100MPa/m(粘土);100GPa(岩石)反分析方法：通过断层中岩石的变形与原岩的变形45GeoHohaiGeoHohai主要内容主要内容nFLAC3D软件简介软件简介nFLAC3D的基本原理的基本原理nFLAC3D的前后处理的前后处理n流流-固耦合分析固耦合分析n接触单元与应用接触单元与应用n完全非线性的动力分析完全非线性的动力分析n自定

24、义本构模型的基本方法自定义本构模型的基本方法n结构单元及应用结构单元及应用46GeoHohaiGeoHohai完全非线性的动力分析完全非线性的动力分析n特点特点n动力荷载动力荷载n动力边界条件动力边界条件n地震波的调整地震波的调整n动孔压的生成动孔压的生成47GeoHohaiGeoHohaiFLAC3D动力分析特点动力分析特点n完全非线性分析完全非线性分析n遵循任何指定的非线性本构关系遵循任何指定的非线性本构关系 n不同频率间会出现干涉和混合不同频率间会出现干涉和混合 n模拟不可恢复的位移和模拟不可恢复的位移和永久变形 n合适的塑性理论，塑性应变增量与应力有关合适的塑性理论，塑性应变增量与应力

25、有关 n易进行不同本构模型的对比分析易进行不同本构模型的对比分析 48GeoHohaiGeoHohai动力荷载动力荷载n动力输入的类型动力输入的类型加速度时程速度时程应力(压力)时程力时程nAPPLY INTERIOR(内部)TABLEFISH49GeoHohaiGeoHohai动力边界条件动力边界条件n静态静态(quiet,粘性粘性)边界边界Lysmer and Kuhlemeyer(1969)模型边界法向和切向设置独立的阻尼器 n自由场自由场(free field)边界边界Cundall et al.(1980)自由场网格与主体网格的耦合粘性阻尼器，自由场网格的不平衡力施加到主体网格边界上

26、设置条件n底部水平，重力方向为z向n侧面垂直，法向分别为x,y向n其他边界条件在APPLY ff之前50GeoHohaiGeoHohai力学阻尼力学阻尼n瑞利瑞利(rayleigh)阻尼阻尼假设阻尼与质量、刚度的线性关系参数确定简单计算速度慢，不推荐n局部局部(local)阻尼阻尼FLAC3D的静力分析阻尼参数简单适合简单情况51GeoHohaiGeoHohai滞回阻尼滞回阻尼(Hysteretic Damping)n模拟岩土介质的动模量衰减曲线模拟岩土介质的动模量衰减曲线ninitial damp hysteretic namesig3(三参数)sig4(四参数)Hardin n(哈丁模型)

27、defaultn计算速度快计算速度快推荐52GeoHohaiGeoHohai地震波的调整地震波的调整n基线校正基线校正对于地震分析的加速度时程，其积分得到的速度和位移应归0美国地质调查研究所nBasic Strong-Motion Accelerogram Processing Software(BAP)对网格施加一个固定速度从而使残余的位移变为053GeoHohaiGeoHohai地震波的调整地震波的调整n动力荷载的频率与单元尺寸的双向调整动力荷载的频率与单元尺寸的双向调整高频的输入要求单元尺寸很小一定的单元尺寸对应输入的最大频率一般进行滤波处理n滤掉低能量的高频nFFT.FISnOrigi

28、n54GeoHohaiGeoHohai地震波的调整地震波的调整El-Centro波FFT修正后的时程修正后FFT5Hz55GeoHohaiGeoHohai动孔压的生成动孔压的生成液化液化n干沙剪应变循环加载试验干沙剪应变循环加载试验初始加载阶段，沙土通常先压实再膨胀。卸载时，沙土遵循与加载相似的路径，但在零应变时，有些残余体积应变存在。取决于初始孔隙率，这可能代表纯粹的压实n假定孔隙中充满水假定孔隙中充满水对于常体积测试，有效应力降低，孔隙水压保持不变对于常荷载测试，(例如，盒子上法向荷载固定)，孔隙水压增加，有效应力减小n有效应力为零时发生液化有效应力为零时发生液化56GeoHohaiGeo

29、Hohai动孔压的生成动孔压的生成液化液化n因此孔隙水压增加不是液化的基本原因因此孔隙水压增加不是液化的基本原因n由于颗粒间由于颗粒间(重组以后重组以后)的低接触力导致有效应力的低接触力导致有效应力的减小的减小n描述液化的模型描述液化的模型高级模型：BSHP(边界面低塑性本构模型,Wang et al.1990)简单模型：MC+体积应变增量模型nFinn模型：nByrne模型：57GeoHohaiGeoHohai主要内容主要内容nFLAC3D软件简介软件简介nFLAC3D的基本原理的基本原理nFLAC3D的前后处理的前后处理n流流-固耦合分析固耦合分析n接触单元与应用接触单元与应用n完全非线性

30、的动力分析完全非线性的动力分析n自定义本构模型的基本方法自定义本构模型的基本方法n结构单元及应用结构单元及应用58GeoHohaiGeoHohai自定义本构模型的基本方法自定义本构模型的基本方法n必要性必要性试验总结的本构模型特定条件下的本构模型交叉学科的本构模型n二次开发环境二次开发环境n自定义本构模型的功能自定义本构模型的功能n自定义本构模型的基本方法自定义本构模型的基本方法59GeoHohaiGeoHohai二次开发环境二次开发环境nFLAC3D采用面向对象的语言标准采用面向对象的语言标准C+编写编写n本构模型都是以动态连接库文件本构模型都是以动态连接库文件(.DLL文件文件)的形的形式

31、提供式提供nVC+6.0(SP4)或更高版本的开发环境或更高版本的开发环境n优点优点自定义的本构模型和软件自带的本构模型的执行效率处在同一个水平 自定义本构模型(.DLL文件)适用于高版本的FLAC(2D)、3DEC、UDEC等其他Itasca软件中60GeoHohaiGeoHohai自定义本构模型的功能自定义本构模型的功能n主要功能：对给出的应变增量得到新的应力主要功能：对给出的应变增量得到新的应力n辅助功能：辅助功能：模型名称、版本读写操作n模型文件的编写模型文件的编写基类(class Constitutive Model)的描述成员函数的描述模型的注册模型与FLAC3D之间的信息交换模型

32、状态指示器的描述 61GeoHohaiGeoHohai自定义本构模型的基本方法自定义本构模型的基本方法n头文件头文件(usermodel.h)中进行新的本构模型派生类中进行新的本构模型派生类的声明的声明修改模型的ID(100)、名称和版本修改派生类的私有成员nC+文件文件(usermodel.cpp)中修改模型结构中修改模型结构(UserModel:UserModel(bool bRegister):ConstitutiveModel)nconst char*UserModel:Properties()函数函数模型的参数名称字符串nconst char*UserModel:States()函数

33、函数计算过程中的状态指示器62GeoHohaiGeoHohai自定义本构模型的基本方法自定义本构模型的基本方法ndouble UserModel:GetProperty()和和void UserModel:SetProperty()函数函数nconst char*UserModel:Initialize()函数函数参数和状态指示器的初始化，并对派生类声明中定义的私有变量进行赋值nconst char*UserModel:Run()函数函数由应变增量计算得到应力增量，从而获得新的应力nconst char*UserModel:SaveRestore()函数函数对计算结果进行保存。n程序的调试程序

34、的调试在VC+的工程设置中将FLAC3D软件中的EXE文件路径加入到程序的调试范围中，并将FLAC3D自带的DLL文件加入到附加动态链接库(Additional DLLs)中，然后在Initialize()或Run()函数中设置断点，进行调试；在程序文件中加入return()语句，这样可以将希望得到的变量值以错误提示的形式在FLAC3D窗口中得到。63GeoHohaiGeoHohai一个例子一个例子(Duncan-Chang)64GeoHohaiGeoHohai主要内容主要内容nFLAC3D软件简介软件简介nFLAC3D的基本原理的基本原理nFLAC3D的前后处理的前后处理n流流-固耦合分析固

35、耦合分析n接触单元与应用接触单元与应用n完全非线性的动力分析完全非线性的动力分析n自定义本构模型的基本方法自定义本构模型的基本方法n结构单元及应用结构单元及应用65GeoHohaiGeoHohaiFLAC3D中的结构单元中的结构单元n有限单元有限单元n梁梁(beam)单元单元n锚索锚索(cable)单元单元n桩桩(pile)单元单元锚杆:rockbolt onn壳壳(shell)单元单元n格栅格栅(geogrid)单元单元土工织物；土工格栅n初衬初衬(liner)单元单元beamcablepileshellgeogridliner66GeoHohaiGeoHohaiCase-4结构的动力响应结

36、构的动力响应config dynsel pile id=1 beg 0 0 0 end 0 0 1sel pile prop dens 2400&Emod 1.0e10 Nu 0.3 XCArea 0.3&XCJ 0.16375 XCIy 0.00625 XCIz 0.01575&Per 2.8 CS_sK 1.3e11 CS_nK 1.3e11&CS_nGap off sel node fix x y z xr yr zr ran id=1sel set damp combined def f1 whilestepping f0=10000*sin(10*dytime)np=nd_head

37、loop while np#null if nd_pos(np,1,3)=1 nd_apply(np,1)=f0 endif np=nd_next(np)endloopendsolve age 1 pile10000*sin(10*t)67GeoHohaiGeoHohaiCase-4结构的动力响应结构的动力响应68GeoHohaiGeoHohai结构单元的应用结构单元的应用n土与结构的相互作用土与结构的相互作用桩基；基坑；边坡锚固地下硐室的支撑结构；采矿；盾构土工织物；土工合成材料n结构不宜复杂结构不宜复杂岩土工程软件，不宜单纯的结构分析复杂结构的模拟很困难结构单元仍不完善nplot显示n双向

38、接触结构(挡土墙)n结构单元的厚度69GeoHohaiGeoHohaiFLAC3D的使用心得的使用心得n 很好的岩土工程专业软件本构模型丰富强大的动力分析功能对内存的要求不高可开发性好nFISHnC+n 时步受网格尺寸的影响很大长期动力计算的时间很长长期渗流计算的时间很长前处理功能需提高结构单元仍需完善手册说明仍需完善70GeoHohaiGeoHohaiFLAC3D的学习经验的学习经验nManuals for Details.(dynamax)n中间时步表现出一些不合实际的结果，需要足够中间时步表现出一些不合实际的结果，需要足够的专业和数学知识进行判断与解释。的专业和数学知识进行判断与解释。(qinjianshe)n充分了解你写的每一条命令的含义。充分了解你写的每一条命令的含义。(qinjianshe)n少量单元的数值试验帮助理解软件的功能少量单元的数值试验帮助理解软件的功能n关键变量的参数化编程设计关键变量的参数化编程设计n参数参数 模型模型 方法方法(软件软件)n由简到繁，循序渐进由简到繁，循序渐进n-SimWe-仿真论坛仿真论坛 A17:FLAC3D/FLAC2D71GeoHohaiGeoHohai欢迎大家访问欢迎大家访问72GeoHohaiGeoHohai谢谢大家!敬请大家批评指正!73