二次等参数体元在接触问题中的应用.pdfVIP

维普资讯 72一 21- Vo12I No 2 Jun． 1997 国 燕山大学 齐齐齐略尔藕； 是计学台 ～ ．． ‘／ ·√ 摘要探讨了用二次等参数体元解算三维接触问题时的载荷移置接触f 蠡 理等问题 ．给出了载荷移置矩阵及其逆矩阵——光滑矩阵 ．用考题验证了上述两个矩阵，并 应用这两个矩阵解算了工程中的三维接触问题 ． · 关键词 自盘掌 !塑童堡丝 骂 望睦． 自 分类号 0343．3 J 0 引言 在使用有限元法分析工程结构中经常遇到实体之间的接触问鄹 寸(例如机床太件之间 的接触．万向节轴滚子轴承内圈与外圈之间通过滚子的接触．直升飞机自动倾斜器中导 筒与滑筒自 等。广泛应用二次等参牧体元 ．但由于二次(以及二次以上)等参数体元 的构成特性．致使在计算中出现一些特殊现象，给有限元接触计算的过程和计算结果的 整理造成困难 ．这包括： 1)把作用在体元表面的分布载荷向单元该表面节点移置时，得到的等效节点载荷与 分布载荷的符号相反； 2，在接触面上同为接触点的角节点和边中节点，其接触法向力的符号不同： 3)采用常规的接触判据和相应的程序计算接触问题时，只能单独使用角节点或边中 节点作为接触 ： 4)不能用常规的方法由计算出的接触节点力得到接触面压力分布 ． 为了解决 匕述问题，本文导出了二次等参数体元的载荷移置矩阵和它的逆矩阵一 光滑矩阵．载荷移置矩阵用于完成某些分布载荷的移置计算；光滑矩阵的作用是采用一 致载荷法由等效节点力反推分布载荷，以解决接被判断和对计算结果的整理等问题 ． 1 载荷移置矩阵和光滑矩阵 若单元某表面只作用着法向分布力，则载荷移置公式为 】=l_Jv]。PdS (1) 借用等参数单元的插值函数，把分布力P写成如下形式liJ： ‘ = ∑ =_Jv]) (2) l 式中， 为等参体元变载表面节点i处的分布力集度 ．把式(2)代入式 (1)，因为 {)的 元素 是节点值．则有 )=f_Jv][MdS ) (3) 将式 (3)积分．．有 1994年6月17日收到 维普资讯 、 第2期 王守信 等 二欢荨参数体元在接触问题中的应用 173 )=[GJ ) (4) 式 (4)中的 [G_即为—个单元的载荷移置矩阵 ．对于多个单元的表面．可按有限元方式 集合．形成整体的载荷移置矩阵 [ ． 由[GJ矩阵求逆，即可得到光滑矩阵 [GJ_。．于是．结构表面分布载荷可以由等效 节点载荷表示为 ]=[q 以j 2 考题 1)载荷移置 图1给出由3+--次等参数体元组成的结构 (只画出变载表面)，其表面作用着线性分布载荷 ．底部最大载 荷集度为 i35kPa． 对于由多个等参数体元组成的结构，式 (4)可写为 )=[G] )