装配体有限元分析.docx

装配体有限元分析.docx

《装配体有限元分析.docx》由会员分享，可在线阅读，更多相关《装配体有限元分析.docx（10页珍藏版）》请在冰豆网上搜索。

装配体有限元分析

基于ANSYSWORKBENCH的装配体有限元分析

模拟装配体的本质就是设置零件与零件之间的接触问题。

装配体的仿真所面临的问题包括：

（1）模型的简化。

这一步包含的问题最多。

实际的装配体少的有十几个零件，多的有上百个零件。

这些零件有的很大，如车门板；有的体积很小，如圆柱销；有的很细长，如密封条；有的很薄且形状极不规则，如车身；有的上面钻满了孔，如连接板；有的上面有很多小突起，如玩具的外壳。

在对一个装配体进行分析时，所有的零件都应该包含进来吗？

或者我们只分析某几个零件？

对于每个零件，我们可以简化吗？

如果可以简化，该如何简化？

可以删除一些小倒角吗？

如果删除了，是否会出现应力集中？

是否可以删除小孔，如果删除，是否会刚好使得应力最大的地方被忽略？

我们可以用中面来表达板件吗？

如果可以，那么，各个中面之间如何连接？

在一个杆件板件混合的装配体中，我们可以对杆件进行抽象吗？

或者只是用实体模型？

如果我们做了简化，那么这种简化对于结果造成了多大的影响，我们可以得到一个大致的误差范围吗？

所有这些问题，都需要我们仔细考虑。

（2）零件之间的联接。

装配体的一个主要特征，就是零件多，而在零件之间发生了关系。

我们知道，如果零件之间不能发生相对运动，则直接可以使用绑定的方式来设置接触。

如果零件之间可以发生相对运动，则至少可以有两种选择，或者我们用运动副来建模，或者，使用接触来建模。

如果使用了运动副，那么这种建模方式对于零件的强度分析会造成多大的影响？

在运动副的附近，我们所计算的应力其精确度大概有多少？

什么时候需要使用接触呢？

又应该使用哪一种接触形式呢？

（3）材料属性的考虑。

在一个复杂的装配体中所有的零件，其材料属性多种多样。

我们在初次分析的时候，可以只考虑其线弹性属性。

但是对于高温，重载，高速情况下，材料的属性不再局限于线弹性属性。

此时我们恐怕需要了解其中的每一种材料，它是超弹性的吗？

是哪一种超弹性的？

它发生了塑性变形吗？

该使用哪一种塑性模型？

它是粘性的吗？

它是脆性的吗？

它的属性随着温度而改变吗？

它发生了蠕变吗？

是否存在应力钢化问题？

如此众多的零件，对于每一个零件，我们都需要考察其各种各样的力学属性，这真是一个丰富多彩的问题。

（4）有限元网格的划分。

我们知道，通过WORKBENCH，我们只需要按一个按钮，就可以得到一个粗糙的网格模型。

但是如果从HYPERMESH的角度来看，ANSYS自动划分的网格，很多都是不合理的，质量较差而不能使用。

那么对于装配体中的每个零件，我们该如何划分网格？

对于每一个零件，我们是否要对之进行切割形成规则的几何体后，然后尽量使用六面体网格？

如果我们这样做的话，那么单单划分网格这一项，就要消耗我们大量的时间。

而且，当这种网格划分完以后，我们还需要反复加密网格，反复计算，直到结果的收敛。

就如同减速器这样的一个装配体，稍微粗略的划分网格，都是10万多个节点，如果我们网格划分得细密一些，很容易上百万个节点。

这么大量的节点，一般的笔记本和台式机计算起来都很困难。

这给我们的仿真工作带来了极大的困扰。

这些问题都是前处理中出现的。

如何解决这些问题，恐怕要我们广大的CAE工程师和CAE研究人员共同努力，从各个侧面进行研究，得到一些个别的成果，然后在某些时候，再集成起来，得到具有普遍指导意义的方法和结论。

ANSYS WORKBENCH提供的六种接触类型

不少朋友提到了关于接触类型的问题，对于如何使用接触类型弄不清楚。

为了帮助刚入门的朋友们了解这些接触类型，笔者首先翻译了ANSYS关于接触类型的帮助，然后对之进行点评。

翻译的部分帮助如下：

ANSYSWORKBENCH提供了6种接触类型，这些接触类型大多只对面接触使适用。

（1）bonded.使用绑定以后，在接触面或者接触边之间不存在切向的相对滑动或者法向的相对分离。

这是缺省的接触类型，适用于所有的接触区域（实体接触，面接触，线接触）。

（2）noseparation.这与绑定类似。

在接触面或者接触线之间不允许发生法向的相对分离，但是允许发生少量的切向无摩擦滑动。

（3）frictionless:

用于模拟无摩擦的单边接触。

所谓单边接触，就是说，一旦两个物体之间出现了分离，则法向力就为零。

因此当外力发生改变时，接触面之间可能会分开，也可能会闭合。

这种情况下假设摩擦系数为零，即当发生切向相对滑动时，没有摩擦力。

（4）rough：

与无摩擦接触类型相似。

它模拟非常粗糙的接触，保证两个物体之间只是发生静摩擦，而不会发生切向的滑移，从而不会产生滑动摩擦。

它相当于在两个物体之间施加了无限大的摩擦系数。

（5）frictional:

有摩擦的接触。

这是最实际的情况，两个接触面之间既可以法向分离，也可以切向滑动。

当切向外力大于最大静摩擦力后，发生切向滑动。

一旦发生切向滑动后，会在接粗面之间出现滑动摩擦力，该滑动摩擦力要根据正压力和摩擦系数来计算。

此时需要用户输入摩擦系数。

（6）forcedfrictionalsliding:

该选项只对刚体动力学适用。

它与frictional类型类似，只是没有静摩擦阶段。

此时，系统会在每个接触点上施加一个切向的阻力。

该切向阻力正比于法向接触力。

到底使用哪种接触类型，取决于你需要解决的问题。

如果

（1）需要模拟两个物体之间轻微的分离

（2）要获得接接触面附近的应力，那么可以考虑下列三种接触类型：

frictionless,rough和frictional.它们可以模拟间隙，并能更精确的建模真实的接触区域。

不过使用这三种接触会导致更长的求解时间，也可能会导致收敛问题。

如果出现了收敛问题，那么可以对接触区域使用更细的网格。

笔者的点评如下：

装配体的分析中，如何对两个物体之间的连接关系进行建模是一个关键技术问题。

对于连接关系，总体考虑如下：

（1）如果两个相邻物体在分析中始终不会有相对运动，最好直接在DM中用多体部件来表达，这最省事。

（2）如果两个相邻物体在分析中存在相对运动，而我们并不关注其连接点附近的应力情况，那么用运动副来表述更简单。

（3）如果相邻两物体在分析中有相对运动，而且我们对这种相对运动的接触面及其附近点的应力情况感兴趣，那么使用接触。

关于接触类型的分类问题。

实际上，接触就是依据两个物体之间是否有切向和法向的相对分离来进行划分的。

在两个相互接触的物体之间，也只能发生这两种运动。

要么，在法线方向上可以分开；要么在切线方向上可以发生相对移动。

如果

（1）法线方向不可分开，切线方向也不可发生相对滑动，则使用boneded。

（2）法线方向不可分开，切线方向可以发生轻微的无摩擦滑动，则使用noseparation.

（3）法线方向可以分开，切线方向不可以发生相对滑动，则用rough.

（4）法线方向可以分开，切线方向可以发生相对滑动，且没有摩擦力。

则是frictionless。

（5）法线方向可以分开，切线方向可以发生相对滑动，存在摩擦力。

则是frictional。

基于Ansys Workbench的接触分析例子1

前面一篇基于Ansys经典界面的接触分析例子做完以后，不少朋友希望了解该例子在Workbench中是如何完成的。

我做了一下，与大家共享，不一定正确。

毕竟这种东西，教科书上也没有，我只是按照自己的理解在做，有错误的地方，恳请指正。

1．问题描述

一个钢销插在一个钢块中的光滑销孔中。

已知钢销的半径是0.5units,长是2.5units，而钢块的宽是4Units,长4Units,高为1Units,方块中的销孔半径为0.49units,是一个通孔。

钢块与钢销的弹性模量均为36e6，泊松比为0.3.

由于钢销的直径比销孔的直径要大，所以它们之间是过盈配合。

现在要对该问题进行两个载荷步的仿真。

（1）要得到过盈配合的应力。

（2）要求当把钢销从方块中拔出时，应力，接触压力及约束力。

2．问题分析

由于该问题关于两个坐标面对称，因此只需要取出四分之一进行分析即可。

进行该分析，需要两个载荷步：

第一个载荷步，过盈配合。

求解没有附加位移约束的问题，钢销由于它的几何尺寸被销孔所约束，由于有过盈配合，因而产生了应力。

第二个载荷步，拔出分析。

往外拉动钢销1.7units，对于耦合节点上使用位移条件。

打开自动时间步长以保证求解收敛。

在后处理中每10个载荷子步读一个结果。

本篇只谈第一个载荷步的计算。

3．生成几何体

上述问题是ANSYS自带的一个例子。

对于几何体，它已经编制了生成几何体的命令流文件。

所以，我们首先用经典界面打开该命令流文件，运行之以生成四分之一几何体；然后导出为一个IGS文件，再退出经典界面，接着再到WORKBENCH中，打开该IGS文件进行操作。

（3.1）首先打开ANSYSAPDL14.5.

（3.2）然后读入已经做好的几何体。

从【工具菜单】-->【File】-->【ReadInputFrom】打开导入文件对话框

找到ANSYS自带的文件

\ProgramFiles\AnsysInc\V145\ANSYS\data\models\block.inp

【OK】后四分之一几何模型被导入，结果如下图

（3.3）导出几何模型

从【工具菜单】】-->【File】-->【Export】打开导出文件对话框，在该对话框中设置如下

即把数据库中的几何体导出为一个block.igs文件。

【OK】以后该文件被导出。

（3.4）退出ANSYSAPDL14.5.

选择【OK】退出经典界面。

4．打开AnsysWorkBench，并新建一个静力学分析系统。

结果如下图

导入几何体模型。

在Geometry单元格中，选择ImportGeometry-->Browse，如下图

找到上一步所生成的block.igs文件。

则该静力学系统示意图更新如下。

可见，几何单元格后面已经打勾，说明文件已经关联。

5．浏览几何模型

双击Geometry单元格，打开几何体。

在弹出的长度单位对话框内，选择米（Meter）的单位。

然后按下工具栏中的Generate按钮如下图

则主窗口中模型如下图

可见，长方形的变长是2m，这与题目中给定的大小是一致的。

然后退出DesignModeler，则又重新回到WorkBench界面中。

6．定义材料属性

双击EngineeringData，则默认材料是钢材。

这里直接修改该钢材的属性即可。

只有线弹性材料属性：

弹性模量36E6和泊松比0.3

然后在工具栏中选择“ReturnToProject”以返回到WorkBench界面中。

7．创建接触

在主窗口中分别选择目标面，接触面如下

然后对该接触的细节面板设置如下

其中，

（1）说明接触类型是带摩擦的接触，摩擦系数是0.2，是非对称接触