P_Chap__界面非线性(DOC可编).doc

3-1 概要 在midas FEA中提供图3.1.1所示的点、线、面三种类型界面单元。另外线和面单元既可以使用低阶单元也可以使用高阶单元。 (a) 二维四节点节点界面单元 (b) 二维六节点界面单元 (c) 三维八节点界面单元 (d) 三维十六节点界面单元 图3.1.1 midas FEA的界面单元 本手册中为了便于理解公式以二维分析为主进行说明，另外下标表示法线方向，和表示切线方向。 界面单元上发生的界面力向量(traction vector)分为法线和切线方向(见式(3.1.1))。此时界面力的单位为与应力单位相同。 (3.1.1) 界面单元上发生的相对位移(relative displacement)如下： (3.1.2) 在此，界面力向量和相对位移之间的关系为线性时，可以表达如下： (3.1.3) 其中： : 法线方向罚刚度(penalty stiffness) (单位：) : 切线方向罚刚度(单位：) 因为界面单元上的非线性特性反应在节点上，使用高斯积分法不能正确描述节点上的非线性响应，其结果也会在界面处发生摆动。为了克服这样的缺点，在midas FEA中使用了牛顿-科特斯(Newton-Cotes)积分公式。 将上面式(3.1.3)用增量形式表现如下： (3.1.4) 其中： : 界面力向量 : 相对位移增量 此时，本构矩阵表达式如下： (3.1.5) 3-2 离散裂缝 图3.2.1 离散裂缝的相对位移与界面力的关系 离散裂缝的本构关系使用全量理论用总变形来计算界面力。用全量理论可以定义法线和切线方向的非线性特性，表达式如下： (3.2.1) 上面的公式可以使用切线刚度定义如下： (3.2.2) 一般来说法线方向界面力会引发I型裂缝，具有受拉软化特性。在midas FEA中提供下面三种受拉软化本构模型。 1. 脆性裂缝模型(brittle cracking model) 2. 线性受拉软化模型(linear tension softening model) 3. 非线性受拉软化模型(nonlinear tension softening model) 3-2-1 脆性裂缝模型 脆性裂缝模型使用了完全脆性模型，达到最大应力以后的关系曲线如下图所示。 图3.2.2 脆性裂缝模型 (3.2.3) 3-2-2 线性受拉软化模型 线性软化模型中法线方向的应力与位移的关系如下图所示。 图3.2.3 线性受拉软化模型 (3.2.4) 其中： : : I型裂缝的断裂能 : 抗拉强度 在线性软化模型中提供了割线逼近法和弹性接近法两种卸载方法。割线逼近法是指卸载时使用与原点相连的直线的刚度卸载，超过原点时重新使用初始刚度的方法；弹性接近法是指卸载时使用初始刚度的方法。 3-2-3 非线性软化模型 Hordijk、Cornelissen、Reinhardt将混凝土的受拉软化定义为如下图指数变化方式。 图3.2.4 非线性受拉软化模型 (3.2.5) 其中： : : : 除了割线和弹性卸载方法以外，还可以选择如下的滞回模型(hysteresis model)卸载。 图3.2.5 滞回模型 3-2-4 剪力滞留法(Shear Retention) I型裂缝在发生裂缝后剪切界面力将减小，其公式如下： (3.2.6) 其中： : 折减的剪切刚度(reduced shear stiffness) 的范围为。发生I型裂缝后假设截面非常光滑则取，如果像混凝土那样不光滑则取，一般的取值范围为0.1～0.3。3-3 膨胀裂缝 在界面上，当相对位移大于时将发生裂缝。另外像混凝土那样裂缝表面不是很光滑的材料，在粗糙的裂缝表面骨料会互相咬合，所以当切线方向发生相对位移时会引起法向的膨胀(dilatancy)。即法线方向的位移与切线方向的位移不是相