FLACFLAC3D在岩土与采矿的应用 PT.pptVIP

力学阻尼 瑞利(rayleigh)阻尼 假设阻尼与质量、刚度的线性关系 参数确定简单(等价平均应变=60%*emax) 中心频率(共振计算，地震平均频率) 临界阻尼比 计算速度慢，不推荐 局部(local)阻尼 FLAC3D的静力分析阻尼 参数简单 适合简单情况 滞回阻尼(Hysteretic Damping) 模拟岩土介质的动模量衰减曲线 initial damp hysteretic name sig3 (3参数) sig4 (4参数) Hardin(1参数) (哈丁模型) Default（2参数) 优点 直接采用模量降级曲线 阻尼比不会影响时间步 缺点 输出的曲线会不一致 滞回阻尼(Hysteretic Damping) 低循环应变下得到的阻尼比要小于试验结果，这会导致低级的噪声，尤其在高频情况下。可以在中心频率上增加一个小量的Rayleigh阻尼(~0.2%刚度比例)，这样也不会降低时步； 若初始应力不为0，剪应力-剪应变曲线可能不匹配。因此在生成初始应力时就要调用Hyst阻尼； Hyst阻尼不仅会增加能量损失，还会导致在大循环应变下的平均剪切模量的降低，在输入波的基频接近共振频率的时候，由于可能会导致动力反应幅值的增大； Hyst阻尼之前要做一次弹性无阻尼求解，以获得发生循环应变的最大水平，若循环应变过大导致剪切模量过多的降低，那么用Hyst阻尼是有问题的； 即使应变较小，使用屈服模型也会增大应变，因此若有广泛屈服的现象，则使用屈服模型，不用Hyst阻尼 参考了Itasca的中国培训资料? 地震波的调整 基线校正 对于地震分析的加速度时程，其积分得到的速度和位移应归0 美国地质调查研究所 Basic Strong-Motion Accelerogram Processing Software (BAP) 对网格施加一个固定速度从而使残余的位移变为0 动力荷载的频率与单元尺寸的双向调整 高频的输入要求单元尺寸很小 一定的单元尺寸对应输入的最大频率 一般进行滤波处理 滤掉低能量的高频 FFT.FIS Origin SeismoSignal 地震波的调整 El-Centro波 FFT 修正后的时程 修正后FFT 5Hz 动孔压的生成——液化 干沙剪应变循环加载试验 初始加载阶段，沙土通常先压实再膨胀。卸载时，沙土遵循与加载相似的路径，但在零应变时，有些残余体积应变存在。取决于初始孔隙率，这可能代表纯粹的压实 假定孔隙中充满水 对于常体积测试，有效应力降低，孔隙水压保持不变 对于常荷载测试，(例如，盒子上法向荷载固定)，孔隙水压增加，有效应力减小 有效应力为零时发生液化 动孔压的生成——液化 因此孔隙水压增加不是液化的基本原因 由于颗粒间(重组以后)的低接触力导致有效应力的减小 描述液化的模型 高级模型：BSHP (边界面低塑性本构模型, Wang et al. 1990) 简单模型：MC + 体积应变增量模型 Finn模型： Byrne模型： 主要内容 软件简介 基本原理 前后处理 FISH语言简介 流-固耦合分析 初始应力的生成 接触单元与应用 完全非线性的动力分析 自定义本构模型的基本方法 结构单元及应用 基本功能 渗流各向同性、各向异性 不同的渗流模型和属性 流体压力，涌入量，渗漏量和不渗水边界 抽水井、点源、体积源 饱和渗流可采用显式差分法、隐式差分法 非饱和渗流采用显式差分法 渗流-固体-热的耦合 流体和固体的耦合程度依赖于土体颗粒(骨架)的压缩程度，用Biot系数表示颗粒的可压缩程度。 循环荷载引起的动水压力变化和土体液化。 地下水模拟方法 有效应力计算 不耦合 孔压为了正确计算有效应力 渗流计算已得到孔压分布 饱和 有水面线的部分饱和 力学变形产生孔压 静力或者动力 不排水孔压或液化 流固耦合计算模式 时间比例 有效应力计算 不设置CONFIG Fluid 孔压不改变 设置孔压分布 INITIAL pp WATER table WATER density WATER table face SET gravity 手动设置干湿密度 设置CONFIG fluid 设置土体干密度 渗流模型 MODEL fl_isotropic MODEL fl_anisotropic MODEL fl_null SET fluid off set WATER bulk = 0 无渗流模式 渗流模式 不耦合计算 CONFIG fluid SET mech off 正确的渗透系数 得到孔压分布和水面线 稳态渗流可以减小Kf CONFIG fluid SET fluid off 正确的流体模量Kf 不需要渗透系数 单渗流模式 单力学模式 流固耦合计算 CONFIG fluid SET fluid on 真实的流体模量K