计算物理 ›› 2015, Vol. 32 ›› Issue (4): 386-394.

• 论文 • 上一篇    下一篇

三维光滑接触界面上的分配参数滑移面算法

姜玉曦, 周海兵, 熊俊, 刘文韬   

  1. 北京应用物理与计算数学研究所, 北京 100094
  • 收稿日期:2014-07-01 修回日期:2014-12-03 出版日期:2015-07-25 发布日期:2015-07-25
  • 作者简介:姜玉曦(1975-),男,江西上饶,副研究员,博士,从事计算流体力学和接触力学研究,E-mail:jiang_yuxi@iapcm.ac.cn
  • 基金资助:
    中国工程物理研究院基金(2012B0201028,2012B0201027,2012A0201011,2013A0101004,2013A0201009)资助项目

A Distributed Parameter Contact Algorithm for Sliding Surfaces on a Three-dimensional Smoothed Contact Surface

JIANG Yuxi, ZHOU Haibing, XIONG Jun, LIU Wentao   

  1. Institute of Applied Physics and Computational Mathematics, Beijing 100094, China
  • Received:2014-07-01 Revised:2014-12-03 Online:2015-07-25 Published:2015-07-25

摘要: 针对三维滑移界面上的分配参数滑移算法,提出一种三维滑移面光滑目的 ,仅根据单元面顶点的坐标和法向等局部信息,构造出曲面片边界G1连续(切平面连续),曲面片内C1连续的三维光滑曲面,此目的 为实际的滑移面提供一个准确的几何表达.在此光滑曲面上构建了具有C1连续的法向压力场和质量面密度场.然后在光滑滑移界面上通过Newton-Raphson迭代计算接触点,依据分配参数滑移算法施加接触约束和接触计算.数值算例显示光滑的接触滑移面可以减少滑移节点的振动,提高计算的收敛性.

关键词: 接触算法, 分配参数滑移算法, 三维滑移面光滑, 拉氏方法

Abstract: In distributed parameter contact algorithm on three-dimensional contact surface,a method for smoothing three-dimensional contact surface is developed. In the method,a smooth contact surface is interpolated with bi-cubic parametric patches. A bi-cubic parametric patch is defined by using local information such as coordinates and normals of vertexes of patch. The smoothing technique provides an accurate representation of actual contact surface which is C1 continuous in patches and G1continuous(tangent plane continuous) between adjacent patches. Then C1 continuous smooth fields of normal pressure and density on smoothed surface are interpolated. And contact points on smoothed contact surface are calculated with Newton-Raphson iteration. Finally,contact constraints are computed with distributed parameter contact algorithm. Numerical examples demonstrate that smoothed contact surface alleviates‘chatter’of nodes and improves convergence behavior.

Key words: contact algorithm, distributed parameter contact method, three-dimensional contact smoothing method, Lagrangian method

中图分类号: