weixin_70982324的博客

粘弹性土体具有阻尼特性，应力作用在模型后的变形有阻尼作用而产生的时间差，因此应变的发展之后于应力的发展。局部阻尼是FLAC3D静力计算所采用的阻尼形式，与频率无关，无需据估计被建模系统的固有频率，但是优点：进行简单模型震动波的动力分析；对于地震模拟，不推荐局部阻尼，因为这种类型的阻尼不能完全取代多次循环加载的能量损失；f2最小中心频率：最小频率与三倍最小频率包含了最多能量的区域的此时的最小频率；2.FLAC力学阻尼：瑞利阻尼、滞后阻尼、局部阻尼、人工粘度。Ⅱ对比瑞利阻尼，滞后阻尼不影响动力计算的时间步；

①对于震动荷载源于模型内部（ 如地铁、隧道中的列车振动、内部爆破等问题）， 可以将震动荷载直接施加在节点上，这种情况下使用静态边界可以有效减小人工上的反射现象， 并且不需要再施加自由场边界。5.自由场边界：通过在模型边界侧面四周生成二维网格，角落生成一维网格的方法来实现；注意：关于荷载施加过程中的计算 ：动力学分析过 程中应避免静力荷载的变化。1.边界：抵消或者吸收振动波的反射响应；②刚性基础一般选择速度或者加时程，自由场边界即可；②如果动力源在顶部或者底，则用自由场更能抵消反射对结果的影响；

计算模型中存在刚度差异较大、模型网格尺寸不均匀的情况时，flac3d可以采用“动态多步”的过程来有效减小计算所需要的时间。如图所示-----滤波以及基线校正后。3.动态多步（分区使用时间步计算）时程转换与地震波滤波。

注意：将间隔指定为0.0，则假定随后是不均匀间隔的数据对。，添加Table Chart显示时程曲线。2.table地震波的文件形式。4.fish动荷载输入。

答：渗流场很多流体相关的参数都是在节点上的，但是力学模型参数基本上都是在zone上，如何实现历遍后关联到zone上！①渗流时长Timestep={（zone尺寸^2+孔隙率）/（渗透系数k+比奥模量M）} * 步数。注：：自定义时间步的时间不长约≥1.25 倍渗流特征长度（时间步）时，软件会报错。问：渗流场状态的改变关联到力学参数上，最难解决的问题是什么？4、渗流场参数关联力学参数的自定义耦合方法。6、流固耦合的控制和章节总结。1、渗流时间步与渗流时长。3、流固耦合的选择依据。

zone gridpoint fix pore-pressure 孔压值，可以跟随一个固定的孔压值，表示所选范围内的所有节点拥有相同的孔压；③比奥系数：孔隙压力变化时，材料单元中流体体积的获得（或损失）与该单元体积变化的的比值。一般地，力学扰动都是瞬时的，而渗流往往需要更长的时间，比如土体中超孔隙水压力的消散与固结往往要发生数小时、数天的时间；在设置特定的孔压边界时，可以配合initial pore-pressure使用，表示有效应力不变，总应力变化；①比奥系数的默认值=1，孔隙率的默认值=0.5；

混合离散化程序使它们更准确地响应涉及等体积塑性变形的情况(例：在莫尔-库仑剪切破坏期间);③更精细的网格导致更准确的结果，因为它们提供了更好的高应力梯度表示。⑤如果需要不同的区域大小，那么从最小到最大的变化越缓慢，结果就越好。然而，在构建几何复杂的模型时，它们比四面体区域更难使用。混合离散化过程一般使其响应更精确(略高阶响应)。②计算速度随区域数量的增加而降低。④准确度增加，区域纵横比趋于统一。

zone interface 1/"b"（id或者“名字”） node property ……keyword ……注：可以换算成young's modulus杨氏模量、poisson's ratio泊松比。Tip ：用空间范围约定接触面建立的区域， 只能附在两个分组构件表面，⑤FLAC3D接触面是单面性质的地质界面（设置参数组合仅为一套）③接触面构件等于或者小于接触的面单元网格大小；①选小原则，大小不一接触表面时；④接触面接触构件（非模型表面）注意：加上x和y的约束范围。②选密原则，上下网格密度不一。

后面补上~~

我感觉地应力平衡的效果不太好，按照ABAQUS的地应力平衡标准，位移结果因为得很小~我也不太懂啥原因！第一次把coh和ten最大，第二次把coh和ten改成正常值；ini赋值，可以加速平衡，效果一样；，后面加关键词ratio设置；适用于浅埋的地表工程。人为设置侧压力系数，

1-①.瞬时开挖sinstantaneous excavation。逐步开挖gradual excavation。1.模型建立&表面分组。

注：①Flac3D对外表面的位移约束首先Face Apply更加符合本地计算系统逻辑，并试图在多个约束之间自动进行协调，结果更可靠；2.随时间逐渐增加边界应力大小：zone face apply stress-xx -1e5。1.边界条件（真实边界&人工边界）：应力边界、位移边界（位移设定&速度设定）注：面上的力为1e10，四周的力为5e9，四个角点的力为2.5e9。1.解除应力边界（例如动荷载计算）：zone face。使用FISH或者Table对计算步进行改变应力边界，1.表面转角处理，注意二者不同！

(2) Anisotropic (Transversely) Elastic Model各向异性模型。(3) Orthotropic Elastic Model正交各向异性模型。(1) Isotropic Elastic Model各向同性模型。①剪切膨胀：摩尔库伦模型、节理模型、应变软化节理模型。Tip：每个六面体都可以划分为五个四面体，②剪切硬化/软化：应变软化模型、节理模型。④抗拉软化：应变软化模型、节理模型。③体积硬化/软化：修正剑桥模型。Null Model空模型。

注意与zone delete range group 'CC' 的区别;③范围设置可用分组名字，空间或者两者结合。Tip：①一个计算执行之前，后命令会覆盖前命令。②不设置范围默认全局施加。Modeifier 修饰。

比率：模型中所有网格点的最大不平衡力大小除以所有网格点的平均外加力大小；当比率在计算过程中低于值时，计算将停止。指定倾角和倾角方向时，FLAC3D假设:正Y为北(0度方位角)，正X为东(+90度方位角)，和正Z为上。Plugin：若不打开此本构关系，FISH插件可以加载，但是循环是禁用的（目前还不太理解 hhhhh）5.孔隙水压力：孔隙压力-流体孔隙压力在压缩过程中是正的；①阻尼的意义：动力学模拟静力学防止一直震动，实现虚拟能力的消散；3.正应变：张拉为正，压缩为负；1.正应力：拉为正，压为负；

B8--砖形 W6--楔形 P5--金字塔形 B7--退化砖形 T4--四面体形。Griddle：非流形独立闭合空间默认分组或者不同批次的gs网格即可分组；1.节点数据表头：节点标识；4.FACE面表头（节点ID）：面标识符；3.ZONE GROUPS分组表头（单元ID）：分组及槽子名称。2.单元表头：单元识别符；'新分组名字' <slot '名字'>5.FACE GROUPS（面ID）group '旧组名'

dip-direction F（度）：指定了全局x0y平面中从正y轴顺时针测量的平面的倾角方向，默认=0；dif F（度）：指定了全局x0y平面以负Z方向测量的平面的倾角，默认=0。ratio：控制网格的密度分布，下一个单元尺寸是上一个单元尺寸的尺寸比。能控制模型空间大小尺寸的只有point和edge和dimension。Size：控制网格的密度（XYZ等方向单元的数量）镜像的两种表达方法：①法向和某一点；edge：模型的边长（20*20*20）

是指求解过程中始终采用塑性模型，但为防止在计算过程中出现屈服区域，将粘聚力和抗拉强度设为大值，计算至平衡后，再将粘聚力和抗拉强度改为分析所采用的值计算至最终平衡状态。

命令自动初始化一个由重力引起的应力场。这个命令不考虑几何形状中的任何不规则性。这个命令默认了垂直应力和水平应力相等，即ratio=1.0. 因此在这种情况下所有方向的应力是相同的。立即删除开挖单元：zone cmodel assign null range group "Space"或者zone delete~逐步开挖：zone relax excavate range group "Space"