P01_Chap_03_界面单元9.docVIP

3-1 概要 界面单元是为了模拟相同材料上的凹凸面之间的接触或异质材料之间的接触而开发的单元，主要用于描述混凝土不均匀的离散式裂缝、钢筋与混凝土间的滑移、钢筋和混凝土的接触或砌体结构的接触面等接触形式。 界面单元虽然使用一般的有限元公式，但单元的厚度为0。 为了进行数值分析，使用罚函数方法将相对位移和界面力的关系定义为如下: (3.1.1) 二维时、、的向量如下，边界上的一个积分点的相对位移和界面力如图3.1.1所示。 ，， (3.1.2) (a) 相对位移 (b) 界面力 图3.1.1 二维分析时界面的相对位移和界面力 若是三维的情况，则各成分可表示为式（3.1.3）和图3.1.2。 ,， (3.1.3) (a) 相对位移 (b) 界面力 图3.1.2 三维分析时的界面相对位移和界面力 虚线表示边界面，是法向界面力，和表示切向界面力。且是法向相对位移, 和是切向相对位移。 3-2 坐标系和相对位移 在界面单元中是由界面单元上端（top）和下端（bottom）的单元的位移决定的。界面单元上端节点上的位移可表示如下: (3.2.1) 单元内任意坐标x、y、z和平动位移u、v、w可表示如下: ，， ，，, (3.2.2) 界面单元下端节点上的位移和坐标也可采用相同的方法表示。即在式(3.2.1)和(3.2.2)将上标“top”换成上标“bot”，见下式(3.2.3）和（3.2.4）。 (3.2.3) ，， ，， (3.2.4) 在结构分析中位移使用边界上单元的坐标系，但界面单元的刚度使用法向——切向组成的坐标系，界面单元不传递旋转刚度，因此式（3.2.5）是单元轴方向的相对位移。 (3.2.5) 其中， 将界面单元上的相对位移矩阵采用分块矩阵定义如下： (3.2.6) 其中， 3-3 点界面单元 图3.3.1 点界面单元 点界面单元如图3.3.1所示由边界面上的坐标轴和界面上下的节点—上节点和下节点组成。上节点和下节点可以有相同的坐标位置，用户可以自由设定边界面的坐标轴。为了方便说明在图3.3.1中将两节点分开表示。应注意的是在midas FEA中用户定义的单元坐标轴决定了界面单元的坐标轴。单元轴的x、y、z轴的顺序决定了界面单元的法向、 剪切向和切向坐标轴的方向，因此点对点界面单元如图3.3.1所示由下节点指向上节点的方向为该界面单元的x轴。 单元的任意节点在整体坐标系的坐标是通过形函数()确定的，如式（3.3.1）所示： (3.3.1) 某一点上的位移由等参形函数定义如下： (3.3.2) 上下节点的等参形函数如式（3.3.3）所示： (3.3.3) 因此相对位移与单元位移矩阵的关系可定义如下： (3.3.4) 3-4 线界面单元 图3.4.1 线界面单元 在midas FEA中线界面单元提供如图3.4.1所示的低阶和高阶两种单元。线界面单元常被用于模拟平面单元与平面单元或平面单元与线单元的接触。 刚度矩阵D与式(3.1.2)相同，且由于线界面单元是2D单元，因此与使用同样的刚度值。 线界面单元的任意点在整体坐标系中的坐标可采用如下式(3.4.1)所示的函数定义。 (3.4.1) 带括号的表示高阶单元。 某一点的整体位移如下式（3.4.2）所示： (3.4.2) 线界面单元的等参形函数可定义如下所示。 ，