收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

弹性接触问题的边界面法分析及其低秩近似快速算法

郑兴帅  
【摘要】:在机械工程领域中,机械部件之间通过接触关系进行运动和载荷的传递,如齿轮啮合、滚动轴承等。接触过程中往往伴随着应力集中现象的产生,这也是机械部件出现早期裂纹、裂纹扩展甚至疲劳断裂的重要原因。因此,接触问题的研究具有重要的工程实际意义。随着CAE技术的发展,有限元法、边界元法等数值计算方法被用于分析工程中的接触问题。然而,无论是有限元法还是边界元法,它们的CAE模型中都只包含了网格的几何信息,CAE模型和CAD模型彼此相互独立。因此,从CAD实体模型到CAE网格模型,就不可避免地引入了离散几何误差。边界面法是一种基于边界积分方程和计算机图形学的新的边界类型数值计算方法。它将边界积分方程和CAD造型系统中的边界表征(B-Rep)数据结构有机地结合起来,在CAD模型的曲面参数空间中进行边界积分方程的离散和变量插值近似,积分单元上的几何数据直接通过曲面的参数空间获得,而不是根据网格近似计算,避免了离散几何误差。而且,在边界面法中位移和面力同时作为方程的未知量,可以通过位移和面力直接施加接触条件。本文基于边界面法理论实现了弹性接触问题的数值分析,解决了接触非线性计算过程中的一些关键问题。同时,对于边界面法系数矩阵的低秩近似,本文在研究自适应交叉近似算法的基础上,尝试提出基于节点几何信息的无迭代交叉近似算法。本文的主要研究工作和成果如下:(1)将边界面法扩展到多域弹性力学问题。接触问题中,接触条件的单边不等式约束特点,即便是线弹性材料的相互接触,也往往产生非线性。但是,从数值分析的角度来看,弹性接触问题是一种具有特殊界面条件的多域问题。因此,作为分析弹性接触问题的基础,本文首先将边界面法扩展到多域弹性力学问题中,设计独立的交界面C++数据结构,用于处理多域问题和接触问题中的交界面条件,也为下面实现基于曲面参数空间的接触搜索算法提供便利。(2)提出基于曲面参数空间的接触搜索算法。边界面法的数据结构中不仅包含了离散的网格信息,还包含了边界曲面的参数化信息。本文在边界面法的理论框架下,提出在接触曲面的参数空间中实现接触搜索。首先计算节点到整个接触面的投影过程,可以得到接触面上投影点的参数坐标;然后,在接触面的参数空间中引入四叉树结构,通过树结构的分级特性,可以快速搜索到投影点所在的单元(单元在参数空间中所对应的区域)。(3)采用点对面(NTS)接触离散模型,建立无摩擦和摩擦接触约束方程。在接触问题的数值算法中,接触离散模型定义了接触体之间运动和载荷的传递方式。早期的点对点(ntn)接触离散模型,要求在接触面上节点与节点一一对应,在节点对之间建立接触约束方程,网格划分困难;本文将有限元接触算法中广泛使用的点对面(nts)接触离散模型引入到边界面接触算法中,在接触边界上通过单元插值的方式,定义投影点与单元节点的接触约束关系,从而建立无摩擦和摩擦接触问题的接触约束方程。这样使得两个接触体可以独立划分网格,降低了划分网格难度。(4)提出一种接触迭代算法,实现边界面法分析无摩擦接触问题。由于不考虑摩擦影响,最终的接触状态与加载路径无关。因此,无摩擦接触问题的非线性来自于接触边界的事先未知。本文采用迭代过程处理无摩擦接触的非线性。首先假设初始的接触区,根据初始假设建立接触约束方程,组装整体系统方程并进行第一次计算。然后判断接触区节点是否满足接触条件(非穿透性条件和压应力条件)。如果不满足接触条件,则更新接触约束重新进行计算。每个迭代步中,需要不断修正接触区的大小,直至接触区上所有节点满足接触条件。(5)提出一种增量-迭代算法,实现边界面法分析摩擦接触问题。在摩擦接触问题中,由于摩擦过程是不可逆的,最终的接触状态与加载路径有关,需要采用增量形式的加载过程。由于接触边界和摩擦状态均事先未知,摩擦接触问题具有更强的非线性。本文借鉴有限元法中处理摩擦接触问题的增量技术,提出了一种增量-迭代算法。在法向接触中,根据非穿透接触条件和压应力条件对接触区进行迭代计算。在切向接触中,根据库伦摩擦定律判断节点的摩擦状态并建立相应的接触约束方程。当增量步中发生穿透时,根据节点的穿透量来细分该增量步,并重新开始该增量步的计算。在增量步之间,将上一增量步接触区节点的摩擦状态和滑移角作为下一增量步开始时的初始假设,然后进行迭代计算。(6)提出基于节点几何信息的无迭代交叉近似算法。无论是快速多极算法还是自适应交叉近似算法都是利用了边界积分方程中基本解的远场衰减特性,对其远场矩阵快进行低秩近似。在边界面法中采用高次单元时,有效的数值积分流程往往是:先在一个积分单元上计算得到数值积分所需的数据(如高斯积分点坐标、法向量、雅可比系数等),然后对所有源点循环。如果反过来,用一个节点对所有单元循环,则需要对单元上的积分数据进行重复计算,降低了计算效率。在自适应交叉近似算法的迭代过程中,需要对单元上的积分数据进行重复计算,而且每个迭代步中都要估计低秩近似的误差,误差估计占整个迭代过程计算量的主要部分。无需迭代的交叉近似算法显然能大大提高计算效率。因此,本文尝试提出了基于节点几何信息的无迭代交叉近似算法,可以在数值积分之前确定骨架点的位置,这样就可以继续采用一个单元对所有源点进行循环的计算方式,避免了对单元上的积分数据进行重复计算。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 赖宝全,邓贵杜;基于维度的系统边界面行为分析[J];系统辩证学学报;2005年02期
2 梁平英;陈梦成;;双材料自由边界面端部问题的数值分析[J];应用力学学报;2008年02期
3 钦亚洲;陈建龙;熊孝波;刘成;;土体各向异性边界面模型研发的几个问题探讨[J];力学与实践;2013年05期
4 麦文建,李崇银;论振荡型对流[J];中国科学(B辑 化学 生物学 农学 医学 地学);1982年08期
5 张福征;胡存;刘海笑;;饱和粘土循环动力分析边界面模型的积分算法[J];港工技术;2012年06期
6 王湛,甄永辉;混凝土本构关系边界面模型[J];工程力学;1999年03期
7 林国钧;石明奎;;线性二级规划的最优边界面算法[J];贵州大学学报(自然科学版);1995年01期
8 李兴照;王录民;范量;;软粘土的循环加载边界面模型[J];河南工业大学学报(自然科学版);2006年04期
9 刘艳秋;;边界面模型隐式积分算法的开发[J];中国水运(下半月);2013年09期
10 刘艳秋;胡存;刘海笑;;一种适用于饱和黏土循环动力分析边界面塑性模型的隐式积分算法[J];岩土力学;2013年12期
11 伊颖锋;;基于边界面理论的小应变土体本构模型[J];市政技术;2006年04期
12 B·赫雷拉;J·帕拉斯;吴承平;;流动中边界面的识别[J];应用数学和力学;2010年09期
13 胡存;刘海笑;;适用于饱和黏土循环动力分析的新型边界面塑性模型[J];水利学报;2011年10期
14 黄茂松;杨超;崔玉军;;循环荷载下非饱和结构性土的边界面模型[J];岩土工程学报;2009年06期
15 王刚;张建民;;边界面模型在MARC中的开发实现及应用[J];岩土力学;2006年09期
16 费康;刘汉龙;;边界面模型在ABAQUS的开发应用[J];解放军理工大学学报(自然科学版);2009年05期
17 张卫华;赵成刚;傅方;;饱和砂土相变状态边界面本构模型[J];岩土工程学报;2013年05期
18 孙吉主;施戈亮;;循环荷载作用下接触面的边界面模型研究[J];岩土力学;2007年02期
19 栾茹;;三维气流数值仿真中悬浮边界面处理技术的研究[J];系统仿真学报;2008年13期
20 徐舜华;郑刚;徐光黎;;循环荷载下砂土的剪切硬化边界面本构模型[J];岩土力学;2010年01期
中国重要会议论文全文数据库 前1条
1 魏星;黄茂松;;天然结构性粘土的各向异性边界面模型[A];第一届中国水利水电岩土力学与工程学术讨论会论文集(上册)[C];2006年
中国博士学位论文全文数据库 前6条
1 郑兴帅;弹性接触问题的边界面法分析及其低秩近似快速算法[D];湖南大学;2015年
2 李源;弹性动力学问题的时域边界面法[D];湖南大学;2015年
3 陆陈俊;基于边界面法的CAE与CAD一体化研究[D];湖南大学;2016年
4 王现辉;自适应快速多极子边界面算法分析大规模声学问题[D];湖南大学;2013年
5 覃先云;边界面法的单元实现及其在复杂结构分析中的应用[D];湖南大学;2012年
6 丰土根;饱和砂土不排水动力特性及多机构边界面塑性模型研究[D];河海大学;2002年
中国硕士学位论文全文数据库 前10条
1 王成;发动机喷油嘴弹性静力问题的边界面法分析[D];湖南大学;2014年
2 向科铭;三轴加载条件下的混凝土边界面模型的研究[D];北京交通大学;2008年
3 刘艳秋;适用于各向异性边界面塑性模型的隐式积分算法[D];天津大学;2012年
4 历永杰;偏平面上单调及循环荷载作用下混凝土的边界面模型研究[D];北京交通大学;2009年
5 余列祥;基于GPU加速的边界面法的研究[D];湖南大学;2013年
6 李湘贺;多域边界面法在稳态热传导和弹性静力学分析中的应用[D];湖南大学;2014年
7 费柏平;基于三维实体的弹性结构屈曲分析及其边界面法实现[D];湖南大学;2012年
8 王伟;二维弹塑性土层的波动数值模拟[D];中国地震局工程力学研究所;2005年
9 于雷;考虑应力水平影响的软土弹塑性模型研究[D];天津大学;2014年
10 向晶晶;循环荷载下基于薄层边界面模型的桩土界面特性分析[D];郑州大学;2014年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978