《【2017年整理】ABAQUS非线性有限元方法在结构接触分析中的应用》由会员分享，可在线阅读，更多相关《【2017年整理】ABAQUS非线性有限元方法在结构接触分析中的应用（10页珍藏版）》请在金锄头文库上搜索。
1、1非线性有限元方法在结构接触分析中的应用ABAQUS 范瑞娟Application of Nonlinear Finite Element Method on Structure Contact AnalysisFAN Ruijuan摘要： 针对某型机缝翼结构中所遇到的接触问题，应用ABAQUS 有限元分析软件，对其进行了边界非线性有限元分析。利用 ABAQUS 的接触算法，真实的模拟了缝翼滑轨和滚柱之间的接触关系，进而得到真实的传力路经和应力分布。关键词： 边界非线性；有限元分析；接触1 引言在有限元结构力学模拟中，非线性的来源有三种:材料非线性、边界非线性和几何非线性。材料非线性是指高应
2、变时材料发生屈服，其应力-应变关系成为非线性；边界非线性是指边界条件随着载荷的变化而发生变化；几何非线性是指在模型分析中位移的大小影响到结构响应的情况。飞机结构设计中遇到的接触关系属于边界非线性问题。在有限元结构分析过程中，接触条件是一类特殊的不连续约束，它允许力从模型的一部分传递到另一部分。因为只有当两个表面发生接触时才会有约束产生，而当两个接触的面分开时，就不存在约束关系了，所以这种约束是不连续的。下面结合飞机前缘缝翼滑轨设计中所遇到的接触问题，就如何应用2ABAQUS 软件，对其进行边界非线性有限元分析予以介绍。2 缝翼滑轨设计中的接触问题飞机的前缘缝翼是民用客机中机翼常用的增升活动面，
3、常规缝翼是通过滑轨在滑轮组架中的运动来改变机翼的翼型，以达到增加升力的目的。缝翼滑轨在滑轮组架滚柱上的支持就是一个典型的接触问题，如图 1 所示。滑轮组架内在每根滑轨的安装位置沿滑轨法向和侧向各布置了两组滚柱。当缝翼翼面上的载荷传到滑轨上时，滑轨受力变形，其上下表面就会有滚柱与滑轨表面发生接触，从而限制滑轨的法向运动；其左右两侧也会有滚柱与滑轨腹板表面发生接触，从而限制滑轨的侧向运动。在结构受载过程中，究竟是哪一个或哪些滚柱与滑轨发生接触，从而为其提供边界约束呢？这就是边界非线性有限元分析所要考虑的主要问题。图 1 滑轨与滚柱之间的接触关系3 接触类问题边界非线性有限元分析方法ABAQUS 是
4、国际上先进的大型通用有限元计算分析软件之一。滚柱与滑轨法向接触滚柱与滑轨侧向接触3特别是它的非线性力学（几何、材料、接触）分析功能有世界领先水平，应用此软件，使得该类接触问题的简便求解成为可能。ABAQUS 软件应用 Newton-Raphson 算法迭代求解非线性问题。ABAQUS 中的接触计算过程如图 2 所示。图中 p 表示从属节点上的接触压力；h 表示从属节点侵入主控表面的距离。图 2 接触逻辑3.1 模型简化缝翼本体有限元模型主要为板杆结构，模型网格较粗，主要用来施加载荷。滑轨和滚柱简化为细网格四面体单元，目的是分析滑轨受力变形后与滚柱的接触状态及滑轨的应力。缝翼加强肋连接耳片与滑轨
5、之间用耦合元（coupling）连接，使得载荷由缝结束增量步开开始增量确定接触状态 撤消约束 加上约束执行迭代检查接触的变化检查平衡点脱开：迭代点接触：迭代h 0 0P 0没有变化收敛脱开 接触不收敛4翼本体传给滑轨。缝翼有限元网格见图 3，2-1 滑轨和滚柱有限元网格见图 4。图 3 缝翼有限元网格图 4 2-1 滑轨及滚柱有限元网格2-1 滑轨及滚轮2-2 滑轨及滚轮 2-3 滑轨及滚轮53.2 接触定义 在 ABAQUS 软件中，在两个结构之间定义接触首先是要创建表面（surfaces）,再创建接触相互作用，使两个可能发生互相接触的表面成队,然后定义控制发生接触表面行为的力学性能模型。A
6、BAQUS 软件中的接触对包括主面和从面，一般来说，材料较硬的构件定义为主面，材料较软的构件定义为从面。缝翼滑轨材料为钛合金，滚柱材料为钢。因此，滚柱定义为主面，滑轨定义为从面，滑轨与滚柱之间的接触关系为无磨擦的硬接触。模型中接触关系定义如图 5 所示。在原始数据文件中以下述格式给出：*Surface, type=ELEMENT, name=_PickedSurf190, internal_PickedSurf190, _S3, S3_PickedSurf190, _S2, S2_PickedSurf190, _S4, S4_PickedSurf190, _S1, S1*Surface, ty
7、pe=ELEMENT, name=_PickedSurf191, internal_PickedSurf191, _S3, S3_PickedSurf191, _S2, S2_PickedSurf191, _S4, S4_PickedSurf191, _S1, S1* Interaction: Int-1*Contact Pair, interaction=IntProp-1_PickedSurf191, _PickedSurf190* INTERACTION PROPERTIES*Surface Interaction, name=IntProp-1.,*Friction0.,表面定义表面定
8、义节点集节点集接触对定义（前者为主面，后者为从面）接触属性6图 5 滑轨与滚柱之间的接触关系定义 3.3 模型加载及约束载荷以集中载荷的形式施加在缝翼本体结构上。约束是将每对滚柱的端头简支，在齿条与操纵系统小齿轮啮合位置建立了柱坐标系，并在此位置约束了柱坐标系下沿滑轨的切向位移，模型加载和约束见图 6。从面主面7图 6 模型加载和约束3.4 有限元应力分析结果应用 ABAQUS 软件对上述模型进行了边界非线性有限元分析。分析结果见图 7 和图 8。图 7 为缝翼及滑轨变形图，图 8 为 2-1 滑轨变形后的应力云图。缝翼最大变形为 12.8mm，滑轨最大应力为421.8MPa。2-1 滑轨局部
9、放大加载位置约束位置812FPM图 缝翼变形图 图 8 缝翼 2-1 滑轨应力云图3.5 结果合理性分析对于缝翼及滑轨这样一个多支点、多传力的复杂结构，力图求解它的精确解几乎是不可能的。这里，根据工程梁理论，对滑轨受力进行定性地分析。缝翼滑轨可简单视为双支点的外伸梁，见图 9。缝翼传来的气动载荷 P 和力矩 M，由滚柱提供的支反力 F1和 F2 相平衡。当滑轨端头作用图 9 所示的载荷时，滑轮组架内就会有一对滚柱与滑轨上、下表面发生接触，由 F1和 F2的载荷方向可以判断，前排上滚轮将提供 F1方向的支反力，后排下滚轮将提供 F2方向的支反力。并且，滑轨最大应力将出现在 F1支点处。图 9 滑
10、轨受力分析简图最大应力 421.8MPa9由 3.4 节变形和应力云图可以清楚的看出，滑轨受力变形后，滑轨上下表面与滚柱之间的接触状态及滑轨应力云图中显示的最大应力位置与上述工程梁理论分析相符合。由次可见，应用ABAQUS 软件分析该类接触问题，其结果是合理、可靠的。3.6 与其他分析软件对比随着计算机技术的突飞猛进，先后出现了越来越多的非线性有限元分析软件，比如 MSC/NASTRAN、MSC/MARC、ABAQUS 等等。以 MSC/NASTRAN 软件为例，如果用它的非线性功能分析缝翼滑轨与滚珠之间的接触问题，就需要用户填写许多 CGAP 元素，否则滑轨表面的接触状态及表面接触压力就无法
11、分析。但是，模型中如果 CGAP 元素太多，又往往影响解的收敛。而 ABAQUS 软件只需要在 CAE 中定义可能发生接触的两个表面，软件自动判断哪些节点在分析过程中可能会发生接触。另外，ABAQUS 软件可以自动调整载荷增量步的大小，因此它可以便捷而有效地得到非线性问题收敛的解。还有 ABAQUS 软件首次采用树状 CAE 界面，使初学用户非常容易接受。4 结束语通过以上分析可以看出，应用 ABAQUS 软件边界非线性有限元分析方法，能够更真实的模拟缝翼滑轨在机翼前缘仓内滑轮架中的支持条件，能够根据结构的受力状态判断滑轨与滚柱之间的接触状态，从而体现真实的传力路经得到准确的应力分布。该方法为
12、一飞院引进 ABAQUS 软件后首次在结构接触分析中的应用，并得10到理想的计算结果。建议在其它类似结构的细节分析中，可优先考虑使用 ABAQUS 软件。它不但降低了对分析人员的经验要求，还可以得到局部结构更准确的受力特性，进而为结构设计和试验提供理论依据。参考文献1 ABAQUS 有限元分析软件入门指南 清华大学出版社2 飞机设计手册 国防工业出版社作者简介： 范瑞娟（1970-）女，陕西省宝鸡市人， 中航第一飞机设计研究院强度研究室工程师 1995 年毕业于西北工业大学飞机系，主要从事飞机机翼强度设计研究。地址：西安市 72 号信箱 303 分箱邮编：710089电话：029-86832903
