面面接触算法研究.docxVIP

面面接触算法研究 1 接触算法的计算方法 接触问题是工程中常见的力学问题。接触问题的特点和难点是接触限制和接触力的未知。由于接触界面的区域大小、位置以及接触状态都是未知的、随时间变化的，因此接触问题表现出显著的非线性特征。接触问题的非线性决定了接触分析过程中需要经常插入接触界面的搜寻判定。有限元法为分析接触问题提供了有效的工具，在进行有限元求解接触问题时，接触体被离散为空间单元集，接触界面的相互作用被转化为离散的单元表面之间的作用或者离散的节点与单元表面的作用。接触问题的动力分析主要包括选取时间积分方案和接触算法，其中接触算法包含接触搜索算法和接触力算法。 接触搜索的目的是确定整个系统中有哪些部位发生了接触或者有哪些原已接触的部分发生了滑移或脱离。在接触动力分析中，接触搜索占较大的计算量，因此接触搜索算法的计算效率至关重要。由于接触搜索是接触力计算的基础，所以其计算精度也非常关键。在一般情况下接触搜索需要先进行全局搜索，再进行局部搜索。通过全局搜索先粗略找到围绕特殊点所有潜在可能的接触单元面。全局算法有主从面算法 局部搜索目的是计算节点到单元面的精确距离，找出接触投影点所在单元面的相对位置关系，从而确定接触节点与单元面间的接触状态。局部搜索算法主要有点面算法 在完成接触搜索之后，接触力求解通常采用拉格朗日乘子法或罚函数法 强震作用下的高混凝土坝系统灾变过程伴随着材料和接触复合非线性问题，对其进行全过程精细化数值模拟和全面深入的抗震安全评价，需要求解未知量高达百万甚至千万级，尽管各种商业软件在求解一般常见问题时可以显示出所谓的强大功能，但在计算方法和计算实施方案上很难适用于解决高混凝土坝复杂工况下地震响应分析所遇到的复杂而又特殊的问题，因此，高性能并行计算是解决该问题的必然选择。 本文拟在研究现有接触算法的基础上，提出多体接触问题快捷高效的面面接触算法。这种算法既可解决搜索盲区问题，又能回避求解非线性投影方程。本文算法采用拉格朗日乘子法动力接触方程的增量形式求解接触力，另外，为了便于并行计算及其程序编制，本文还根据接触面分布特点进行预分区处理，并基于提出的多体接触问题面面接触算法，开发实现其算法的FORTRAN源代码程序，为后续实现高混凝土坝和岩体高边坡静动力分析的并行计算提供支持。 2 搜索法的分析 2.1 全局搜索过程的实体面求解 本文将主从面算法与位码算法相结合，实现接触面的全局搜索。利用主从面算法 本文构建了描述“接触面对”的主面单元和从面单元的共享实常数组：主面单元用正整数表示，从面单元用其对应的负整数表示，当绝对值相同时主面单元和从面单元构成“接触面对”；不同的“接触面对”通过不同的实常数定义；一组实常数可以对应多个边界面。 以下是对一组“接触面对”的全局搜索步骤：（1）搜索预处理：统计“接触面对”单元面上的主节点和从节点个数，计算各单元面外法向量及各节点外法向量。其中节点外法向量由包含该节点的单元法向量取平均得到。（2）根据“接触面对”的位置关系，构造包含接触单元面的最小长方体。在接触问题的有限元计算过程中，全局搜索是相对耗时较大的计算环节，为了避免在每个时间步都进行全局搜索，适当扩大该长方体的区域，确定接触搜索范围。（3）确定立方格的网格尺寸，立方格尺寸与接触面中平均的单元尺寸相接近或可以稍大一些。（4）建立典型的立方格结构，其三个坐标（x,y,z）方向的网格尺寸一致。按先x方向，再y方向，最后z方向的顺序对立方格进行编号。（5）对当前“接触面对”的所有主从节点循环，根据节点坐标确定其所在的立方格的编号 基于以上全局搜索流程编制的程序框图如图1所示。 2.2 局部搜索流程 本文结合点面算法与内外算法完成接触局部搜索。由全局搜索获得“接触面对”的点面关系，采用内外算法判断“接触面对”主面上的节点落入哪些从面接触面单元内，为了避免搜索盲区，本文采用了相关单元节点法向“平均向量”。如图2所示，在接触单元中，节点B的法向量n取单元AB的法向量n 下面给出一组“点面接触测试对”的局部搜索步骤：（1）对测试对中的面单元数进行循环，根据“点”法向量与单元面法向量，进一步判断是否为潜在接触关系。（2）针对潜在接触单元，计算“点”到单元面投影点位置及其贯入量。（3）由内外算法 图4给出了基于局部搜索流程编制的程序框架。 综上所述，本文采用主从面算法与位码算法相结合完成接触全局搜索，点面算法与内外算法相结合完成接触局部搜索。同时在执行接触搜索前预先考虑了接触区域分块划分，使得本文算法可以方便地进行高混凝土坝和岩体高边坡静动力稳定分析的并行计算处理。 3 拉格朗日乘子法 根据可能接触边界的分布特点，将计算域分解成不同的子区域。考虑由N个子区域Ω 式中：ρ为质量密度；c为阻尼系数；u? 不同子区域间的可能接触边界条件为分离状态、黏着状