摩擦接触问题参变量二次规划模型.docxVIP

摩擦接触问题参变量二次规划模型 在机械、航空航天、土木工程等领域有许多接触现象。随着现代工程结构分析的精细化，接触问题的描述和解决成为一个回避的问题。然而，接触问题是最困难的非线性工程问题之一。从数学结构分析的结果来看，等式函数的接触限制是非线性的本质。因此，有效应用力学模型中的接触条件是成功解决接触问题的关键。近年来，该方法在优化领域中采用了不均匀边界法，如拉格兰柱子法、惩罚函数法、补偿法、数学规划法、非光滑半光滑方程法等。 我国钟万勰提出的求解摩擦接触问题的参变量变分原理 增广Lagrange乘子法是为了克服罚函数法的缺点而提出的一种施加约束的最优化方法, 它可以被视为一种约束优化问题的精确正则化方法, 该方法无需过大的罚因子就可以获得高精度的解.在1992年以后, 基于增广Lagrange乘子法构造的摩擦接触算法由于其多重优点而成为一类颇为流行的方法 1 单向应力应变关系的不稳定性解析 为便于比较, 沿用文献[3]的表述方式和假设条件.考虑二维问题, 考虑由两个物体B 并且在接触面两侧几乎处处满足 其中, n 接触边界S 其中:ε 在S (1) 法向单边接触定律 式中, 各式的意义分别为非穿透条件、法向不受拉条件和互补条件. (2) 切向Coulomb摩擦定律 式中:第一式为粘贴条件, 第二式为滑动条件;μ为摩擦因数.仿照理想刚塑性体单向应力应变关系的表述方式, 式 (3) 和式 (4) 的摩擦接触条件还可以统一表述为 式中:g为接触滑动势函数; 参变量变分原理求解摩擦接触问题时, 文献[3]的作法是首先引入惩罚函数, 将接触条件 (5) 转化为如下形式: 式中 本文给出摩擦接触条件的另一种表述形式: 不难证明, 只需E 由式 (7) 中的第一式可得, 式 (7) 的意义在于它提供了一个可实际操作的迭代格式: 联系数学规划中的增广Lagrange乘子法, 迭代格式 (8) 中关于 在实际计算过程中, 随着迭代过程的进行, 当||p 2 应用参变量变分原理求解接触问题 式 (8) 给出的接触条件, 同时适用于二维和三维摩擦接触问题.为了说明简单, 下面仅以平面弹性摩擦接触问题为例, 展示接触条件 (8) 在参变量变分原理解摩擦接触问题中的应用. 2.1 扩大和开发lagranga参与变量的最小势能原理 平面弹性接触力学边值问题在区域Ω=Ω 对于给定的 其中 上述参变量变分原理的策略就是作系统能量指标泛函Π 2.2 增广lagrange线性互补算法 对两接触体Ω=Ω 其中 式中: 其中 因此, 有限元离散后, 平面弹性接触问题的增广Lagrange参变量最小势能原理可转化为由式 (10) 和 (11) 组成的参数二次规划问题, 经过一系列的推导, 可将其进一步转化为下列线性互补问题: 若令M=U-CK 这是一个标准的线性互补问题, 但它包含了3类未知变量: 据此, 给出如下摩擦接触问题的增广Lagrange线性互补算法: (1) 初始化, 设定接触刚度矩阵 (2) 确定 (3) 由式 (14) 得p 3 使用示例 3.1 接触条件的验证 如图1所示, 一个截面为圆形的弹性柱体同刚性基础的接触, 这是一个经典的Hertz接触问题.该问题可以模型化为平面应变问题:力F=1.528, 8, k N, 弹性模量E=19.60, MPa, 泊松比ν=0.3, 截面半径R=8, cm, 并取1/4结构进行分析.有限元模型中1/4柱体剖分为225个平面四边形四节点等参元, 接触边界划分为15个平面点接触单元. 计算过程中, 设定算法的迭代终止条件为 首先考证算法的有效性, 设定接触刚度E 前文指出, 罚因子的选取对于算法的精度不会有太大影响, 但会影响算法的收敛速度.下面的计算经验表明了这一观点.以摩擦因数μ=0.4的问题为例, 分别选取E 从计算量的角度考虑, 本文算法要求多次求解线性互补问题 (13) , 其中迭代的第一步所求解的互补问题与文献[3]中所求解的互补问题相同.显然, 本文算法会增加计算量.为了考证计算量的增加幅度, 对摩擦因数μ=0.4, 惩罚因子选为E 3.2 上部梁的最大挠度 悬臂梁接触问题 结果三种工况均迭代2次即收敛. 结果比较:文献[13]中, 无摩擦时, 位于上部的梁的最大挠度6.717, mm, 本文算法结果为6.717, mm;文献[13]中, 有摩擦情况 (摩擦因数μ=0.2) 时, 位于上部的梁的最大挠度为6.550, mm, 本文算法结果为6.540, mm. 4 增广lagrange线性互补算法 本文通过将增广Lagrange乘子法的思想引入参变量变分原理中, 构造了求解摩擦接触问题的参变量二次规划分析的增广Lagrange线性互补算法.算法的显著优点是在不需要过大的惩罚因子就能获得