带孔平板的线性静力分析
带孔平板的线性静力分析
本示例将对一个给定的带孔平板几何模型创建有限元模型、施加边界条件、进行有限元分析并在HyperView中观察受载平板的变形和应力结果。
本示例包括以下步骤：
?在HyperMesh中建立有限元模型
?施加载荷和边界条件
?求解
?观察结果
1．在HyperMesh中建立有限元模型
（1）载入OptiStruct用户界面并打开模型文件
1）启动HyperMesh。
2）在User Profile对话框中选择OptiStruct，点击OK。
这就加载了OptiStruct用户界面，它包括OptiStruct模板、宏菜单等。简化了与OptiStruct 使用相关的HyperMesh功能。
User Profiles…可以从下拉式菜单中的Preferences中进入。
3）在工具条选择按钮。
弹出Open file…窗口。
4）选择plate_hole.hm文件，模型位于
5）点击Open。
plate_hole.hm的数据被载入当前的HyperMesh中，替代了原有的数据。数据仅包含几何。
注意此时plate_hole.hm的路径显示在file:文本框中。
6）点击Return。
（2）定义材料属性、单元属性卡片及component
1）点击定义材料。
2）在面板左边选择create子面板。
３）点击name =并输入steel。
４）点击card image =并从弹出菜单中选择MAT1
５）点击create/edit。
弹出MAT1 的卡片信息。
如果括号中的量下面没有值，表示其处于关闭状态。要改变该状态，点击括号中的量，
其下会出现输入框。点击输入框输入数值即可。
对选项E输入2e5，对NU输入0.3，对Rho输入7.9e-09。
６）点击 return。
这就创建了一个新材料属性steel，该材料属性属于各向同性材料模型MAT1。该材料的杨氏弹性模量（Young's Modulus）为2E+05，泊松比（Poisson's Ratio）为0.3。仅进行静力分析不需要密度值，但其他求解过程需要用到材料密度值。
可以通过模型树的鼠标右键功能编辑修改。
７）点击创建一个属性卡片。
８）点击 name = 并输入 pshell。
９）点击card image =并从弹出菜单中选择PSHELL。
1０）点击material =并选择steel。
1１）点击create/edit。
弹出PSHELL 的卡片属性card image。
1２）点击Ｔ，输入厚度10。
1３）点击return。
这就创建了一个新的壳单元属性卡片，材料属性为steel中，单元厚度值为10mm，并指定了颜色。可以通过模型树的鼠标右键功能编辑修改。
14）点击return返回主菜单。
15) 点击创建component.
16) 点击 name = 并输入 shells。
17) 设置颜色，点击property,选择属性pshell。
（3）网格划分
利用automesh面板创建四边形单元，名为shells的component刚刚被创建，它是当前的component，生成的单元将放在这个component中。
1）在2D页面上选择automesh面板。
2）如果surfs 不是默认选项，点击entity selection开关并选择surfs。
3）点击surfs并从弹出菜单中选择displayed。
4）点击mesh。
5）为elem size =输入40。
6）点击recalc all。
注意在图形区边界的节点是如何调整的。
7）在面板左边选择mesh style子面板。
8）点击elem type:开关并从弹出菜单中选择quads。
9）点击set all（在elem type:中）。
此时该表面将划分为四边形网格。
10）点击mesh method:并从弹出菜单中选择autodecide。
11）在mesh method:中点击set all。
网格将按最优的运算法则自动调整。
12）点击mesh。
13）网格的创建取决于曲面边界的节点密度、运算法则和所选的网格类型。生成的网格如图4-1所示。
图4-1在automesh 面板中创建的四边形网格（单元大小为40mm）14）点击return。
这一步将网格保存在当前的component（shells）中。
15）点击return返回主菜单。
2．施加载荷和边界条件
在这一节中，模型将施加以下载荷和边界条件：四条外边界中的两相对边不能移动，另外两边没有约束。在孔的边界沿Z轴的正方向施加共1000N的力。
（1）创建载荷集（spcs和forces）
1）点击创建载荷集。
2）在面板左边选择create子面板。
3）点击 name = 并输入spcs。
4）点击color并从调色板中选择一种颜色。
6）点击向下箭头开关并从弹出菜单中选择no card image。
7）点击create。
创建了新的载荷集spcs。
8）点击name = 并输入forces。
9）点击color并选择另一种颜色。
10）点击create。
创建了新的载荷集forces。
11）点击return返回主菜单。
（2）创建约束
1）点击右下角条形按钮。打开了一个新窗口，通过它设置当前的load collectors。或者通过模型树上鼠标右键功能选择。
2）与1同样的窗口下，点击loadcol并从载荷集中选择spcs。
3）从Analysis页面选择constraints面板。
4）在面板左边选择create子面板。
5）点击entity selection开关并从弹出菜单中选择nodes。
6）点击nodes并从弹出的entity selection菜单中选择by window。
7）在图形区域拖拽一个窗口包围所选节点，如图4-2所示。
窗口是多边形，鼠标每点击一次创建一个顶点。
图4-2 需施加约束的节点
8）点击窗口内部（interior）并点击select entities。
9）约束dof1, dof2, dof3, dof4, dof5, 和dof6。
被激活的Dofs表示该自由度被约束，没有激活表示该自由度保持自由。
Dofs 1，2和3表示沿x, y, z方向的平动自由度。
Dofs4，5和6表示绕x, y, z方向的旋转自由度。
10）点击create。
对所选节点施加了约束。
11）点击return返回主菜单。
（3）对孔上的节点施加载荷。
1）点击右下角条形按钮。打开了一个新窗口，将forces设置成当前载荷集。或者通过
模型树上鼠标右键功能选择。
2）在Analysis页面中选择forces面板。
3）在面板左边选择create子面板。
4）点击黄色开关前向下箭头并从弹出菜单中选择nodes。
5）点击nodes黄色开关并从弹出的entity selection菜单中选择by window。
6）在图形区域拖拽一个窗口包围所选节点。
这个窗口是多边形，鼠标每点击一次创建一个顶点。
图4-3显示了应选节点。
图4-3 在孔上施加了载荷的节点
7）点击窗口内部并点击select entities。
8）点击nodes并从entity selection菜单中选择save。
9）点击return退出forces面板。
10）在Tool页面上选择count面板。
11）在面板左边选择model子面板。
12）点击黄色开关前向下箭头并从弹出菜单中选择nodes。
13）点击nodes并从弹出的entity selection菜单中选择retrieve。
14）重新得到了刚保存在forces面板中的节点。
15）点击selected。
这样就计算了孔上节点的数目。
16）记录节点数目。
17）点击return退出count面板。
18）在Analysis页面中选择forces面板。
19）点击nodes并从弹出的entity selection菜单中选择retrieve。
20）设置坐标系为global system）。
21）点击vector definition开关并选择constant vector。
22）点击magnitude =并输入数值使得孔上节点的总载荷为1000，即1000除以圆孔一周的节点数。
23）在magnitude =下面点击direction definition开关，并从弹出菜单中选择z-axis。
24）点击create。
对孔上的节点创建了一系列节点力并给定了Z方向的大小。
25）点击 return 返回主菜单。
（4）创建OptiStruct的载荷工况（即载荷步）
1）在Analysis页面中选择loadstep面板。
2）点击name =并输入lateral force。
3）点击type:开关，选择linear static。
4）选中SPC前面的复选框。
SPC右面出现等号按钮
6）点击等号按钮并从load collectors中选择SPC。
7）选中LOAD前面的复选框。
LOAD右面出现等号按钮。
8）点击等号按钮并从load collectors中选择forces。
9）点击create。
创建了OptiStruct一个工况，该工况引用了load collectors中的约束SPCS和载荷forces。
10）点击 return 返回主菜单。
3．求解
（1）运行OptiStruct
1）在Analysis页面中选择 OptiStruct 面板。
2）在input file:后面点击save as…。
弹出save files…窗口。
3）为OptiStruct模型文件选择路径并在File name:中输入模型文件名plate_hole.fem。
OptiStruct支持以.fem为扩展名的文件。
4）点击Save。
注意plate_hole.fem文件的名称和路径在input file: 文本框中显示。
5）设置memory options:为memory default。
6）点击run options:开关并选择analysis。
7）设置export options:为all。
8）点击OptiStruct。
这就运行了OptiStruct，如果运行成功，在plate_hole.fem所在路径下将出现新的结果文件。plate_hole.out文件可以用于查找错误信息，如果出现错误，这些信息将有助于调试。
4．察看结果与后处理
OptiStruct默认输出线性静力分析的位移（Displacement）和应力（Stress）结果。
本部分描述如何在HyperView中察看这些结果。
HyperView是一个完整的后处理软件和可视化环境，用于处理有限元分析（FEA）、多体系统仿真、视频和工程数据。
（1）察看应力云图
1）当命令窗口信息进程显示成功完成时，点击HyperView（位于OptiStruct面板）。
HyperView启动并加载结果，弹出信息窗口通知模型和结果文件成功载入HyperView 中。
2）点击Close关闭信息窗口。
3）点击工具条按钮Contour 。
4）在Result type:下选择Element Stresses [2D & 3D] (t)。
5）在Result type:下选择vonMises。
6）在Averaging method:）下面选择None。
7）点击Apply。
出现von Mises应力的云图。模型中的每一个单元分配一种颜色，表示该单元在所施加的载荷和边界条件下产生的应力值。
8）在右下方的view controls（视角控件）中点击Top以察看模型，如图4-4所示：
图4-4 von Mises 应力云图
最大应力值是多少？
最大应力出现在模型的什么位置？
将边界条件作用于模型有没有意义？
（2）察看位移云图
1）在Result type:下选择Displacement (v)。
2）在Result type:下选择Mag。
3）点击Apply。
显示在所施加的载荷和边界条件下的位移云图。
最大位移值是多少？
最大位移出现在模型的什么位置？
将边界条件作用于模型有没有意义？
（3）察看变形图
1）在右下方的view controls中点击Iso，显示模型的等角视图。
2）点击Deformed工具条按钮。
3）设置Result type:为Displacement(v)，Scale:为Scale factor，Type:为Uniform。4）将value设为500。
这意味着最大的位移为500个模型单位，其他所有的位移都相应的进行缩放。5）将Show:设为Wireframe。
6）点击Apply。
出现包含位移云图的变形图，覆盖在原有未变形的网格上。图4-5为其等角视图。
图4-5 变形图的等角视图
非线性有限元方法及实例分析
非线性有限元方法及实例分析 梁军 河海大学水利水电工程学院，南京（210098） 摘 要：对在地下工程稳定性分析中常用的非线性方程组的求解方法进行研究，讨论了非线性计算的迭代收敛准则，并利用非线性有限元方法分析了一个钢棒单轴拉伸的实例。 关键词：非线性有限元，方程组求解，实例分析 1引 言 有限单元法已成为一种强有力的数值解法来解决工程中遇到的大量问题，其应用范围从固体到流体，从静力到动力，从力学问题到非力学问题。有限元的线性分析已经设计工具被广泛采用。但对于绝大多数水利工程中遇到的实际问题如地下洞室等，将其作为非线性问题加以考虑更符合实际情况。根据产生非线性的原因，非线性问题主要有3种类型[1]： 1．材料非线性问题（简称材料非线性或物理非线性） 2．几何非线性问题 3．接触非线性问题（简称接触非线性或边界非线性） 2 非线性方程组的求解 在非线性力学中，无论是哪一类非线性问题，经过有限元离散后，它们都归结为求解一个非线性代数方程组[2]： ()()()00 021212211=… …==n n n n δδδψδδδψδδδψΛΛΛ （1.1） 其中n δδδ,,,21Λ是未知量，n ψψψ,,,21Λ是n δδδ,,,21Λ的非线性函数，引用矢量记 号 []T n δδδδΛ21＝ （1.2） []T n ψψψψΛ21= （1.3） 上述方程组（1.1）可表示为 ()0=δψ （1.4） 可以将它改写为 ()()()0=?≡?≡R K R F δδδδψ （1.5） 其中()δK 是一个的矩阵，其元素 是矢量的函数，n n ×ij k R 为已知矢量。在位移有限 元中，δ代表未知的结点位移，()δF 是等效结点力，R 为等效结点荷载，方程()0=δψ表示结点平衡方程。 在线弹性有限元中，线性方程组
非线性有限元分析
轨道结构的非线性有限元分析 姜建华 练松良 摘 要 实际轨道结构受载时的力学行为,属于典型的非线性力学问题。钢轨垫层刚度、钢轨抗扭刚度和扣件扣压力的大小是影响轨距扩大的主要因素。根据非线性有限元接触理论,建立了能准确反映扣件、钢轨与垫层的拧紧接触,以及受载车轮与钢轨侧向滑动接触的力学计算模型;并研究计算了不同扣件压力下,由于受载车轮与钢轨侧向滑动接触引起的轨距扩大问题。 关键词 轮轨关系,扣件压力,非线性弹性力学,有限元分析 1 引言 实际工程中常见的非线性问题一般可以归纳为三类:材料非线性、几何非线性以及边界条件非线性。材料非线性问题是由于材料的非线性本构关系所引起的,例如材料的弹塑性变形,材料的屈服和硬化等;几何非线性问题是由于结构的位移或变形相当大,以至必须按照变形后的几何位置来建立平衡方程;边界条件非线性问题是指边界条件随位移变化所引起的非线性问题。通常情况下,我们所遇到的非线性问题多数是上述三类非线性问题的组合[1,2]。 实际轨道结构受载时的力学行为,属于典型的非线性力学问题。比如基于轮轨接触的材料非线性、几何非线性及边界条件非线性问题,以及扣件、钢轨、垫层三者间相互作用时所表现的边界条件非线性行为等。所以,机车车辆在轨道结构上行驶时引起的力学现象是相当复杂的。以往在研究轨道各部分应力应变分布规律时,通常采用连续弹性基础梁理论或连续点支承,偶尔简单考虑扣件的作用和弹性垫层的使用。不管用哪一种支承方式建立模型,都由于这样那样的假设而带有一定程度的近似性。所以,如何利用现代力学理论的最新成果以及日益发展的计算机技术,根据轨道结构的具体情况,建立更为完整更为准确的轨道结构计算模型,为轨道设计部门提供更加可靠的设计依据或研究思路,已十分必要。 本文提出了用非线性有限元理论研究轮轨系统和轨道结构的思路。作为算例之一,本文将根据非线性有限元理论,建立能准确反映扣件、钢轨与垫层的拧紧接触,以及受载车轮与钢轨侧向滑动接触的力学计算模型。 2 轨道结构的有限元接触模型 对于非线性问题,不管是材料非线性、几何非线性,还是边界条件非线性,总是最终归结为求解一组非线性平衡方程及其控制方程。例如用位移作为未知数进行有限元分析时,最后可得到一组平衡方程及其控制方程为 : 图1 轮轨系统的对称性模型简图 [K(u)]{u}={R}(1) (u)= (u)(2)其中:{u}为节点位移列阵;{R}为节点载荷列阵; [K(u)]为总体刚度矩阵; (u)为边界条件。它们 36 姜建华:同济大学工程力学系,副教授、博士,上海200092
有限元非线性计算特点
有限元非线性计算特点 文章通过几个典型的工程计算模型，分析比较有限元线性与非线性计算结果，阐释了有限元非线性计算的特点及优点。 标签：工程计算；线性；非线性 1 引言 有限元单元法已成为强有力的数值解法来解决工程中遇到的大量问题，其应用范围从固体到流体，从静力到动力，从力学问题到非力学问题，有限元的线性分析已被广泛采用。但对于许多航空工程中遇到的问题，如进气道等，仅仅采用线性求解是不真实的，而采用非线性计算将更符号实际情况。本文借助MSC/NASTRAN有限元分析程序，对于典型的工程计算模型分析比较线性与非线性计算结果，从而给出非线性计算相对于线性计算的优点及特点。 2 有限元非线性计算的特点及优点 为了明确有限元非线性计算结果与线性计算结果的差异，更好的展现有限元非线性计算的特点，本节将借助于有限元分析软件MSC/NASTRAN，对一受外载的矩形薄板根据不同的边界条件，进行非线性及线性静力分析，通过分析比较计算结果，说明有限元非线性静力计算中的一些特点。 2.1 非线性与线性计算结果随载荷的变化 首先，给出薄板尺寸、载荷。 模型尺寸：薄板尺寸为500×500×1.5mm。 载荷：受法向气动压力（pressure），气动压力由小到大变化依次为0.01MPa、0.02MPa、0.04MPa、0.08MPa、0.16MPa。 取薄板中央节点位移、应力及薄板边缘中部节点位移，比较线性计算结果和非线性计算结果。在分别进行有限元线性及非线性分析后，给出位移、应力及支反力结果随载荷的变化曲线。图1、图3、图5分别为采用限元线性计算得到的参考点的位移、应力及支反力变化曲线；图2、图4、图6分别为采用有限元非线性计算得到的参考点的位移、应力及支反力变化曲线。 由圖1、3、5可见，采用线性静力分析后，参考点位移、应力、支反力均随载荷增加而线性增大，位移、应力、支反力与载荷呈明显的线性关系，这是线性静力分析的特点。对于本例，可以预言，在其它条件不变的情况下，计算出一套载荷下的结果，就可以按照线性关系求出压力载荷下的位移、应力及支反力结果。
带孔平板拉伸作业
带孔平板有限元分析 本文采用有限元法，对带圆孔的矩形平板进行了弹塑性受力分析，分析了圆孔处的应力集中现象，为其设计和应用提供了参考依据。 1. 研究问题概述 本文研究带圆孔矩形平板在轴对称拉力作用下的平面应力问题。平板开孔的应力问题是弹塑性力学平面中的一个经典的问题,也是实际工程中常见的问题。平板长200mm ，宽50mm ，厚8mm ，具体几何参数及受力见图1。 图1 平板几何参数及受力 2.弹性力学方法解答 由弹性力学知识知，在距圆孔圆心()r ρρ>处的径向正应力、环向正应力、切应力分别为： 222222 1c o s 211322p r p r r ρσψρρρ?????? =-+-- ? ????????? 22221cos 21322p r p r ?σψρρ????=+-+ ? ???? ? 2222sin 21132p r r ρψψρ ττψρρ???? ==--+ ?????? ? 沿着y 轴，90ψ=。，环向正应力为： 242413122r r p ?σρρ?? =++ ???
max 3q ?σ=由上表可知： ()max = 3K q ψ σ=故应力集中因子: 可见孔边最大应力比无孔时提高了3倍，应力集中系数k=3，如图2所示。 图2 孔边应力集中 3.有限元分析 3.1模型建立 图3 有限元模型 3.2边界条件和载荷 为避免在计算时平板产生移动引发计算问题，必须对试件的外部边界条件进行限定。对平板左侧进行铰接约束，示意图如下
图4 平板约束示意图 由于我们只关注孔附近的应力分布情况，根据圣维南原理，载荷的具体分布只影响载荷作用区附近的应力分布。故我们用均布力代替集中力施加在平板右侧的作用面上，其大小为225P MPa ，为负值。 图5 平板载荷示意图 3.3材料 平板的弹性模量为200GPa ，泊松比为0.3。其塑性的应力应变参数见下图 图6 塑性应力应变参数 3.4有限元网格划分 网格划分是非常重要的过程，它会对计算速度、精度、可靠性产生重要影响。网格划分主要包括两方面：尺寸、单元类型。
线性静力学分析实例
学号：p1******* 姓名：朱四海 线性静力学分析实例 1.1 问题的描述 一部件结构如图1-1所示，一端面受固定约束，另一端A、B两点受相反方向切向力，求受载后的Mises应力、位移分布。 ν 材料性质：弹性模量E=2e5，泊松比3.0 = 图1-1 部件模型 1.2 启动ABAQUS 启动ABAQUS有两种方法，用户可以任选一种。 （1）在Windows操作系统中单击“开始”--“程序”--ABAQUS 6.10 -- ABAQUS/CAE。 （2）在操作系统的DOS窗口中输入命令：abaqus cae。 启动ABAQUS/CAE后，在出现的Start Section（开始任务）对话框中选择Create Model Database。 1.3 创建部件 在ABAQUS/CAE顶部的环境栏中，可以看到模块列表：Module：Part，这表示当前处在Part（部件）模块，在这个模块中可以定义模型各部分的几何形体。可以参照下面步骤创建部件的几何模型。 （1）创建部件。对于如上图1-1所示的部件模型，可以先画出二维截面，再通过拉伸得到。步骤如下：
单击左侧工具区中的（Create Part）按钮，或者在主菜单里面选择Part--Create，弹出如图1-2所示的Create Part对话框。 图1-2 Create Part对话框 在Name（部件名称）后面输入ep2，Modeling Space（模型所在空间）设为3D，Shape选择Solid（实体），Type采用默认的Extrusion，在Approximate size里面输入600。单击Continue...按钮。 （2）绘制部件二维截面图。ABAQUS/CAE自动进入绘图环境，左侧的工具区显示出绘图工具按钮，视图区内显示栅格，视图区正中两条相互垂直的点划线即当前二维区域的X轴和Y轴。二者相交于坐标原点。 选择绘图工具箱中的工具，窗口提示区显示Pick a center point for the circle--or enter X,Y（选择一个中心点的圆，或输入X,Y的坐标），如图1-3所示。 图1-3 输入圆心坐标 输入圆上任意点坐标为（0,50），回车，第一个圆形就画出来了。继续画第二个圆，圆心坐标为（0,0），圆上任意一点（0,40）。
第9章 非线性问题的有限单元法
第9章非线性问题的有限单元法 9.1 非线性问题概述 前面章节讨论的都是线性问题，但在很多实际问题中，线弹性力学中的基本方程已不能满足，需要用非线性有限单元法。非线性问题的基本特征是变化的结构刚度，它可以分为三大类：材料非线性、几何非线性、状态非线性。 1. 材料非线性(塑性, 超弹性, 蠕变) 材料非线性指的是材料的物理定律是非线性的。它又可分为非线性弹性问题和非线性弹塑性问题两大类。例如在结构的形状有不连续变化（如缺口、裂纹等）的部位存在应力集中，当外载荷到达一定数值时该部位首先进入塑性，这时在该部位线弹性的应力应变关系不再适用，虽然结构的其他大部分区域仍保持弹性。 2. 几何非线性(大应变, 大挠度, 应力刚化) 几何非线性是有结构变形的大位移引起的。例如钓鱼杆，在轻微的垂向载荷作用下，会产生很大的变形。随着垂向载荷的增加，杆不断的弯曲，以至于动力臂明显减少，结构刚度增加。 3. 状态非线性(接触, 单元死活) 状态非线性是一种与状态相关的非线性行为。例如，只承受张力的电缆的松弛与张紧；轴承与轴承套的接触与脱开；冻土的冻结与融化。这些系统的刚度随着它们状态的变化而发生显著变化。 9.2 非线性有限元问题的求解方法 对于线性方程组，由于刚度方程是常数矩阵，可以直接求解，但对于非线性方程组，由于刚度方程是某个未知量的函数则不能直接求解。以下将简要介绍借助于重复求解线性方程组以得到非线性方程组解答的一些常用方法。 1.迭代法 迭代法与直接法不同，它不是求方程组的直接解，而是用某一近似值代人，逐步迭代，使近似值逐渐逼近，当达到允许的规定误差时，就取这些近似值为方程组的解。 与直接法相比，迭代法的计算程序较简单，但迭代法耗用的机时较直接法长。它不必存贮带宽以内的零元素，因此存贮量大大减少，且计算中舍入误差的积累也较小。以平面问题 为例，迭代法的存贮量一般只需直接法的14左右。在求解非线性方程组时，一般采用迭代 法。 2. 牛顿—拉斐逊方法 ANSYS程序的方程求解器计算一系列的联立线性方程来预测工程系统的响应。然而，非线性结构的行为不能直接用这样一系列的线性方程表示。需要一系列的带校正的线性近似来求解非线性问题。 一种近似的非线性救求解是将载荷分成一系列的载荷增量，即逐步递增载荷和平衡迭代。它可以在几个载荷步内或者在一个载荷步的几个子步内施加载荷增量。在每一个增量的
结构静力分析
第一章结构静力分析 1．1 结构分析概述 结构分析的定义：结构分析是有限元分析方法最常用的一个应用领域。结构这个术语是一个广义的概念，它包括土木工程结构，如桥梁和建筑物；汽车结构，如车身骨架；海洋结构，如船舶结构；航空结构，如飞机机身等；同时还包括机械零部件，如活塞，传动轴等等。 在ANSYS产品家族中有七种结构分析的类型。结构分析中计算得出的基本未知量（节点自由度）是位移，其他的一些未知量，如应变，应力，和反力可通过节点位移导出。 静力分析---用于求解静力载荷作用下结构的位移和应力等。静力分析包括线性和非线性分析。而非线性分析涉及塑性，应力刚化，大变形，大应变，超弹性，接触面和蠕变。 模态分析---用于计算结构的固有频率和模态。 谐波分析---用于确定结构在随时间正弦变化的载荷作用下的响应。 瞬态动力分析---用于计算结构在随时间任意变化的载荷作用下的响应，并且可计及上述提到的静力分析中所有的非线性性质。 谱分析---是模态分析的应用拓广，用于计算由于响应谱或PSD输入（随机振动）引起的应力和应变。 曲屈分析---用于计算曲屈载荷和确定曲屈模态。ANSYS可进行线性（特征值）和非线性曲屈分析。 显式动力分析---ANSYS/LS-DYNA可用于计算高度非线性动力学和复杂的接触问题。 此外，前面提到的七种分析类型还有如下特殊的分析应用： ●断裂力学 ●复合材料 ●疲劳分析 ●p-Method 结构分析所用的单元：绝大多数的ANSYS单元类型可用于结构分析，单元型 从简单的杆单元和梁单元一直到较为复杂的层合壳单元和大应变实体单元。 1．2 结构线性静力分析 静力分析的定义 静力分析计算在固定不变的载荷作用下结构的效应，它不考虑惯性和阻尼的影响，如结构受随时间变化载荷的情况。可是，静力分析可以计算那些固定不变的惯性载荷对结构的影响（如重力和离心力），以及那些可以近似为等价静力作用的随时间变化载荷（如通常在许多建筑规范中所定义的等价静力风载和地震载荷）。 静力分析中的载荷 静力分析用于计算由那些不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位移，应力，应变和力。固定不变的载荷和响应是一种假定；即假定载荷和结构的响应随时间的变化非常缓慢。静力分析所施加的载荷包括： ●外部施加的作用力和压力 ●稳态的惯性力（如中力和离心力） ●位移载荷 ●温度载荷 线性静力分析和非线性静力分析 静力分析既可以是线性的也可以是非线性的。非线性静力分析包括所有的非线性类型：大变形，塑性，蠕变，应力刚化，接触（间隙）单元，超弹性单元等。本节主要讨论线性静力分析，非线性静力分析在下一节中介绍。
第三章 非线性分析
第三章非线性分析 在工程结构实际中，常常会遇到许多不符合小变形假设的问题，例如板和壳等薄壁结构在一定载荷作用F，尽管应变很小，甚至未超过弹性极限，但是位移较大，材料微单元会有较大的刚体转动位移。这时平衡条件应如实地建立在变形后的位形上，以考虑变形对平衡的影响。同时应变表达式也应包括位移的二次项。这样，结构的几何形变关系将是非线性的。这种由于大位移和大转动引起的非线性问题称为几何非线性问题。在涉及几何非线性问题的有限元方法中，可以采用两种不同的表达格式来建立有限元方程。一种格式是所有静力学和运动学变量总是参考于初始位形的完全拉格朗日格式，即在整个分析过程中参考位形保持不变。而另一种格式中，所有静力学和运动学的变量参考于每一载荷步增量或时间步长开始的位形，即在分析过程中参考位形是不断被更新的，这种格式就称为更新的拉格朗日格式。下面将分别具体讨论大变形情况下应变和应力度量，几何非线性有限元方程的建立以及系数矩阵的形成。 在涉及几何非线性问题的有限元方法中，可以采用两种不同的表达格式来建立有限元方程。一种格式是所有静力学和运动学变量总是参考于初始位形的完全拉格朗日格式，即在整个分析过程中参考位形保持不变。而另一种格式中，所有静力学和运动学的变量参考于每一载荷步增量或时间步长开始的位形，即在分析过程中参考位形是不断被更新的，这种格式就称为更新的拉格朗日格式。下面将分别具体讨论大变形情况下应变和应力度量，几何非线性有限元方程的建立以及系数矩阵的形成。 第三章非线性分析的数值计算方法 3．1概述 非线性问题一般包括三类：材料非线性、几何非线性和边界非线性；而在许多实际的结构中，常常是三种非线性问题的融合，因此其解析方法能够得到的解答是十分有限的。对于非线性问题的求解，可以采用有限元分析的方法，因此非线性方程组的解法也就成为非线性问题有限元分析涉及的基本问题，也就是通常所说的非线性分析的数值计算方法I”。常用的有Newton—Raphson法(简称N-R)和弧长法。本文将详细介绍Newton-Raphson法和弧长法，且依据不同的约束方程形式介绍各种不同形式的弧长法并比较其准确性和可靠性，这在非线性分析计算中是非常有意义的。 3．2牛顿一拉夫森法
带孔平板的线性静力分析
带孔平板的线性静力分析 本示例将对一个给定的带孔平板几何模型创建有限元模型、施加边界条件、进行有限元分析并在HyperView中观察受载平板的变形和应力结果。 本示例包括以下步骤： ?在HyperMesh中建立有限元模型 ?施加载荷和边界条件 ?求解 ?观察结果 1．在HyperMesh中建立有限元模型 （1）载入OptiStruct用户界面并打开模型文件 1）启动HyperMesh。 2）在User Profile对话框中选择OptiStruct，点击OK。 这就加载了OptiStruct用户界面，它包括OptiStruct模板、宏菜单等。简化了与OptiStruct 使用相关的HyperMesh功能。 User Profiles…可以从下拉式菜单中的Preferences中进入。 3）在工具条选择按钮。 弹出Open file…窗口。 4）选择plate_hole.hm文件，模型位于
非线性有限元分析
非线性有限元分析 1 概述 在科学技术领域，对于许多力学问题和物理问题，人们已经得到了它们所应遵循的基本方程（常微分方程或偏微分方程）和相应的定解条件（边界条件）。但能够用解析方法求出精确解的只是少数方程性质比较简单，并且几何形状相当规则的问题。对于大多数工程实际问题，由于方程的某些特征的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析的答案。这类问题的解决通常有两种途径。一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。但是这种方法只是在有限的情况下是可行的，因为过多的简化可能导致误差很大甚至是错误的解答。因此人们多年来一直在致力于寻找和发展另一种求解途径和方法——数值解法。特别是五十多年来，随着电子计算机的飞速发展和广泛应用，数值分析方法已成为求解科学技术问题的主要工具。 已经发展的数值分析方法可以分为两大类。一类以有限差分法为代表，主要特点是直接求解基本方程和相应定解条件的近似解。其具体解法是将求解区域划分为网格，然后在网格的结点上用差分方程来近似微分方程，当采用较多结点时，近似解的精度可以得到改善。但是当用于求解几何形状复杂的问题时，有限差分法的精度将降低，甚至发生困难。 另一类数值分析方法是首先建立和原问题基本方程及相应定解条件相等效的积分提法，然后再建立近似解法并求解。如果原问题的方程具有某些特定的性质，则它的等效积分提法可以归结为某个泛函的变分，相应的近似解法实际上就是求解泛函的驻值问题。诸如里兹法，配点法，最小二乘法，伽辽金法，力矩法等都属于这一类方法。但此类方法也只能局限于几何形状规则的问题，原因在于它们都是在整个求解区域上假设近似函数，因此，对于几何形状复杂的问题，不可能建立合乎要求的近似函数。 1960年，R.W.CLOUGH发表了有限单元法的第一篇文献“The Finite Element Method in Plane Stress Analysis”，这同时也标志着有限单元法（FEM）的问世。有限单元法的基本思想是将连续的求解区域离散为一组有限个，且按一定方式相互联接在一起的单元的组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模型化几何形状复杂的求解域。并且可以利用在每一个单元假设的近似函数来分片地表示全求解域上待求的未知场函数，从而使一个连续的无限自由度问题变成离散的有限自由度问题。 现已证明，有限单元法是基于变分原理的里兹法的另一种形式，从而使里兹法分析的所有理论基础都适用于有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法。利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，而且事先不要求满足任何边界条件，因此可以用来处理很复杂的连续介质问题。 在短短四十余年的时间里，有限单元的分析方法已经迅速地发展为适合于使用各种类型计算机解决复杂工程问题的一种相当普及的方法。如今，有限元广泛地应用于各个学科门类，已经成为工程师和科研人员用于解决实际工程问题，进行科学研究不可或缺的有力工具。有限单元法的应用围已由弹性力学平面问题扩展到空间问题，板壳问题，由静力平衡问题扩展到稳定问题，动力问题和波动问题。分析的对象从弹性材料扩展到塑性，粘弹性，粘塑性和复合材料等，从固体
ABAQUS线性静力学分析实例
荿蚇衿膃蒂葿螅膂膁 线性静力学分析实例 线性静力学问题就是简单且常见得有限元分析类型,不涉及任何非线性(材料非线性、几何非线性、接触等),也不考虑惯性及时间相关得材料属性。在ABAQUS 中,该类问题通常采用静态通用(Static,General)分析步或静态线性摄动(Static,Linear perturbation)分析步进行分析。线性静力学问题很容易求解,往往用户更关系得就是计算效率与求解效率,希望在获得较高精度得前提下尽量缩短计算时间,特别就是大型模型。这主要取决于网格得划分,包括种子得设置、网格控制与单元类型得选取。在一般得分析中,应尽量选用精度与效率都较高得二次四边形/六面体单元,在主要得分析部位设置较密得种子;若主要分析部位得网格没有大得扭曲,使用非协调单元(如CPS4I 、C3D8I)得性价比很高。对于复杂模型,可以采用分割模型得方法划分二次四边形/六面体单元;有时分割过程过于繁琐,用户可以采用精度较高得二次三角形/四面体单元进行网格划分。一 悬臂梁得线性静力学分析 1、1 问题得描述 一悬臂梁左端受固定约束,右端自由,结构尺寸如图1-1所示,求梁受载后得Mises 应力、位移分布。 材料性质:弹性模量32e E =,泊松比3.0=ν 均布载荷:Mpa p 6.0= 图1-1 悬臂梁受均布载荷图 1、2 启动ABAQUS 启动ABAQUS 有两种方法,用户可以任选一种。
在Windows操作系统中单击“开始”--“程序”--ABAQUS 6、10 -- ABAQUS/CAE。 （2）在操作系统得DOS窗口中输入命令:abaqus cae。 启动ABAQUS/CAE后,在出现得Start Section(开始任务)对话框中选择Create Model Database。 1、3 创建部件 在ABAQUS/CAE顶部得环境栏中,可以瞧到模块列表:Module:Part,这表示当前处在Part(部件)模块,在这个模块中可以定义模型各部分得几何形体。可以参照下面步骤创建悬臂梁得几何模型。 （1）创建部件。对于如图1-1所示得悬臂梁模型,可以先画出梁结构得二维截面(矩形),再通过拉伸得到。 单击左侧工具区中得(Create Part)按钮,或者在主菜单里面选择Part--Create,弹出如图1-2所示得Create Part对话框。 图1-2 Create Part对话框
第八章几何非线性问题的有限元法
第八章 几何非线性问题的有限元法 引言 前面各章所讨论的问题都是在小变形假设的前提下进行的，即假定物体所发生的位移远小于物体自身的几何尺寸，应变远小于1。在此前提下，建立物体或微元体的平衡条件时可以不考虑物体的位置和形状（简称位形）的变化，因此在分析中不必区别变形前后位形的差别，且应变可用一阶无穷小的线性应变表达。 实际上，上述假设有时是不成立的。即使实际应变可能是小的，且不超过材料的弹性极限，但如果需要精确地确定位移，就必须考虑几何非线性，即平衡方程应该相对于变形后的位置得出，而几何关系应该计及二次项。例如平板大挠度理论中，由于考虑了中面内的薄膜应力，求得的挠度比小挠度理论的结果有显著的减低。再如在结构稳定性问题中，当载荷达到一定数值后，挠度比线性解答予示的结果更剧烈地增加，并且确实存在承载能力随继续变形而减低的现象。在冷却塔、薄壁结构及其它比较细长的结构中，几何非线性分析都显得十分重要。 几何非线性问题可以分为以下几种类型： （1）大位移小应变问题。一般工程结构所遇到的几何非线性问题大多属于这一类。例如高层建筑或高耸构筑物以及大跨度网壳等结构的分析常需要考虑到结构大位移的影响。 （2）大位移大应变问题，如金属压力加工中所遇到的问题就属于这一类型。 （3）结构的变形引起外载荷大小、方向或边界支承条件的变化等。 结构的平衡实际上是在结构发生变形之后达到的，对于几何非线性问题来说，平衡方程必须建立在结构变形之后的状态上。为了描述结构的变形需要设置一定的参考系统。一种做法是让单元的局部坐标系始终固定在结构发生变形之前的位置，以结构变形前的原始位形作为基本的参考位形，这种分析方法称作总体的拉格朗日（Ｌagrange ）列式法；另一种做法是让单元的局部坐标系跟随结构一起发生变位，分析过程中参考位形是不断被更新的，这种分析方法称作更新的拉格朗日列式法。 本章首先对几何非线性问题作一般性讨论，从中导出经典的线性屈曲问题的公式；然后建立平板大挠度问题和壳体的大位移（及大转动）分析的有限方法公式；接着还给出了大应变及大位移的一般公式，最后还详细讨论了杆系结构几何非线性问题的有关公式。在讨论中我们采用总体的拉格朗日列式法，但对杆系结构，为应用方便我们给出了两种列式法的公式。 & 一般性讨论 理论基础 无论是对于何种几何非线性问题，虚功原理总是成立的。由虚功原理，单元的虚功方程可以写成如下的形式 {}{}{}{}0=-???**v e eT e eT F dv δσε （） 其中{}F 为单元节点力向量，{}e *ε为单元的虚应变，{}e *δ为节点虚位移向量。 增量形式的应变一位移关系可表示为 {}[] {}e e d B d δε= （）
带孔板的建模及有限元分析Word版
基于SolidWorks带孔板的建模及有限元分析 李军 摘要：利用SolidWorks对带孔矩形板进行虚拟建模，通过赋予板材材质、载荷后进行网格划分，进而进行有限元分析，得出其应力、应变和位移的分布图，并对结果进行分析研究对板材安全性的影响。 关键词：SolidWorks；带孔板；建模；有限元分析 0 SolidWorks简介 Solidworks是一款优秀的三维设计软件，具有十分强大的零件设计功能及装配模块，同时也拥有丰富的后置处理模块。由于其功能强大，新手上手快，应用领域广，所以成为了主流的三维造型软件。经过17年的发展，在全球已经拥有30多万的客户，最新版本为SolidWorks 2011版。在中国SolidWorks在计算机辅助设计、计算机辅助工程、计算机辅助制造、计算机辅助工艺、数据管理等方面为企业提供了强大的动力，使企业在管理、设计和制造方面有了很大的提升。 1 带孔板的模型建立 矩形板材的尺寸为300*180*10mm，孔位于中心，直径为50mm，模型如图1。 图1 带孔矩形板模型 2前置处理 2.1在Command Manager中点击SIMULATION选项，建立新算例，名称默认，确认。 2.2赋予板材材料属性 材料为AISI304，材料属性如表1
表1 材料的属性 模型参考属性零部件 名称:AISI 304 模型类型:线性弹性同向性 默认失败准则:最大von Mises 应力屈服强度: 2.06807e+008 N/m^2 张力强度: 5.17017e+008 N/m^2 弹性模量: 1.9e+011 N/m^2 泊松比:0.29 质量密度:8000 kg/m^3 抗剪模量:7.5e+010 N/m^2 热扩张系数: 1.8e-005 /Kelvin SolidBody 1(凸台-拉伸1)(aisi304带孔矩形钢板静力分析) 曲线数据:N/A 2.3网格生成 在SIMULATION选项中选择“运行”中的“生成网格”，使用默认网格划分。网格 信息如表2，网格信息细节如表3，网格划分后的模型如图2。 表2 网格信息 网格类型实体网格 所用网格器: 基于曲率的网格 雅可比点 4 点 最大单元大小7.44196 mm 最小单元大小7.44196 mm 网格品质高 表3 网格信息细节 节点总数23523 单元总数13612 最大高宽比例 3.9347 单元(%)，其高宽比例< 3 99.7 单元(%)，其高宽比例> 10 0 扭曲单元(雅可比)的% 0 完成网格的时间(时;分;秒): 00:00:03 计算机名: PC-201009062016
非线性有限元读书报告
非线性有限元读书报告 1 对张量分析及其在结构非线性力学和有限元中应用的认识。 张量分析将张量中的每一项的计算过程（这些计算过程可能是相似的）用张量之间的计算方法代替，并统一，提供了一种更加方便的计算方法。张量标记中，指标不出现，因此张量标记的表示是独立于坐标系统的，并且可以用于笛卡尔坐标系，柱坐标系，曲线坐标等。如同矩阵一般，或者说张量就是矩阵的更一般的形式，而矩阵是二维的张量。此外，在张量标记中的方程非常容易记忆。 2 对各种应变、应力张量和材料本构关系的认识。 本段回答除本课程讲义外大部分参考文献[1] 2.1 应变张量 为了使刚体运动时，特别是刚体转动时，应变度量为0，工程应变张量被抛弃了。根据采用连续介质力学中的欧拉（Eulerian ）方法或拉格朗日方法（Lagrangian ），得到的应变分别为阿尔曼西（Almansi ）应变（或称欧拉应变） 张量（E ij ε）与格林（Green ）应变（或称拉格朗日应变）张量（G ij ε），它们与工 程应变（或称柯西小应变张量）（此处记为ij ε）均不相同。 三者的几何方程表达式为： 阿尔曼西（Almansi ）应变（或欧拉应变）张量： ()11,,,,11 ()22E ij ij ki kj i j j i k i k j F F u u u u εδ--= -=+-； 格林（Green ）应变（或拉格朗日应变）张量： (),,,,11 ()22G ij ki kj ij i j j i k i k j F F u u u u εδ= -=++； 工程应变（柯西小应变）张量： ,,1()2 ij i j j i u u ε=+。 其中u 为位移，ij F 为变形梯度矩阵，即i ij j x F X ??? =????? 另外，阿尔曼西（Almansi ）应变（或欧拉应变）张量（E ij ε）和格林（Green ）应变（或拉格朗日应变）张量（G ij ε）之间的转换关系为：
ABAQUS线性静力学分析实例
线性静力学分析实例 线性静力学问题是简单且常见的有限元分析类型，不涉及任何非线性（材料非线性、几何非线性、接触等），也不考虑惯性及时间相关的材料属性。在ABAQUS 中，该类问题通常采用静态通用（Static ，General ）分析步或静态线性摄动（Static ，Linear perturbation ）分析步进行分析。 线性静力学问题很容易求解，往往用户更关系的是计算效率和求解效率，希望在获得较高精度的前提下尽量缩短计算时间，特别是大型模型。这主要取决于网格的划分，包括种子的设置、网格控制和单元类型的选取。在一般的分析中，应尽量选用精度和效率都较高的二次四边形/六面体单元，在主要的分析部位设置较密的种子；若主要分析部位的网格没有大的扭曲，使用非协调单元（如CPS4I 、C3D8I ）的性价比很高。对于复杂模型，可以采用分割模型的方法划分二次四边形/六面体单元；有时分割过程过于繁琐，用户可以采用精度较高的二次三角形/四面体单元进行网格划分。 一 悬臂梁的线性静力学分析 问题的描述 一悬臂梁左端受固定约束，右端自由，结构尺寸如图1-1所示，求梁受载后的Mises 应力、位移分布。 材料性质：弹性模量32e E =，泊松比3.0=ν 均布载荷：Mpa p 6.0= 图1-1 悬臂梁受均布载荷图
启动ABAQUS 启动ABAQUS有两种方法，用户可以任选一种。 （1）在Windows操作系统中单击“开始”--“程序”--ABAQUS -- ABAQUS/CAE。 （2）在操作系统的DOS窗口中输入命令：abaqus cae。 启动ABAQUS/CAE后，在出现的Start Section（开始任务）对话框中选择Create Model Database。 创建部件 在ABAQUS/CAE顶部的环境栏中，可以看到模块列表：Module：Part，这表示当前处在Part（部件）模块，在这个模块中可以定义模型各部分的几何形体。可以参照下面步骤创建悬臂梁的几何模型。 （1）创建部件。对于如图1-1所示的悬臂梁模型，可以先画出梁结构的二维截面（矩形），再通过拉伸得到。 单击左侧工具区中的（Create Part）按钮，或者在主菜单里面选择Part--Create，弹出如图1-2所示的Create Part对话框。
研究生课程论文《非线性有限元分析》
1 引言 在科学技术领域内，对于许多力学问题和物理问题，人们已经得到了它们所应遵循的基本方程（常微分方程或偏微分方程）和相应的定解条件（边界条件）。但能够用解析方法求出精确解的只是少数方程性质比较简单，并且几何形状相当规则的问题。对于大多数工程实际问题，由于方程的某些特征的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析的答案。这类问题的解决通常有两种途径。一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。但是这种方法只是在有限的情况下是可行的，因为过多的简化可能导致误差很大甚至是错误的解答。因此人们多年来一直在致力于寻找和发展另一种求解途径和方法——数值解法。特别是五十多年来，随着电子计算机的飞速发展和广泛应用，数值分析方法已成为求解科学技术问题的主要工具。 已经发展的数值分析方法可以分为两大类。一类以有限差分法为代表，主要特点是直接求解基本方程和相应定解条件的近似解。其具体解法是将求解区域划分为网格，然后在网格的结点上用差分方程来近似微分方程，当采用较多结点时，近似解的精度可以得到改善。但是当用于求解几何形状复杂的问题时，有限差分法的精度将降低，甚至发生困难。 另一类数值分析方法是首先建立和原问题基本方程及相应定解条件相等效的积分提法，然后再建立近似解法并求解。如果原问题的方程具有某些特定的性质，则它的等效积分提法可以归结为某个泛函的变分，相应的近似解法实际上就是求解泛函的驻值问题。诸如里兹法，配点法，最小二乘法，伽辽金法，力矩法等都属于这一类方法。但此类方法也只能局限于几何形状规则的问题，原因在于它们都是在整个求解区域上假设近似函数，因此，对于几何形状复杂的问题，不可能建立合乎要求的近似函数。 1960年，R.W.CLOUGH发表了有限单元法的第一篇文献“The Finite Element Method in Plane Stress Analysis”，这同时也标志着有限单元法（FEM）的问世。有限单元法的基本思想是将连续的求解区域离散为一组有限个，且按一定方式相互联接在一起的单元的组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模型化几何形状复杂的求解域。并且可以利用在每一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数，从而使一个连续的无限自由度问题变成离散的有限自由度问题。 现已证明，有限单元法是基于变分原理的里兹法的另一种形式，从而使里兹法分析的所有理论基础都适用于有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法。利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，而且事先不要求满足任何边界条件，因此可以用来处理很复杂的连续介质问题。 在短短四十余年的时间里，有限单元的分析方法已经迅速地发展为适合于使用各种类型计算机解决复杂工程问题的一种相当普及的方法。如今，有限元广泛地应用于各个学科门类，已经成为工程师和科研人员用于解决实际工程问题，进行科学研究不可或缺的有力工具。有限单元法的应用范围已由弹性力学平面问题扩展到空间问题，板壳问题，由静力平衡问题扩展到稳定问题，动力问题和波动问题。分析的对象从弹性材料扩展到塑性，粘弹性，粘塑性和复合材料等，从固体力学扩展到流体力学，传热学等连续介质力学领域。在工程分析中的作用已从分析和校核扩展到优化设计并和计算机辅助设计技术相结合。各种各样商业化的大型通用有限元软件层出不穷，不断推陈出新。可以预见，随着现代力学，计算数学，计算机技术等学科的发展，有限单元法作为一个具有巩固理论基础和广泛应用范围的数值分析工具，必将得到进一步的完善和发展。
过盈配合应力的接触非线性有限元分析
过盈配合应力的接触非线性有限元分析 作者：许小强赵洪伦 摘要基于非线性有限元软件MARC，提出过盈配合应力的动态和静态两种有限元分析方法，并以铁道车辆某高速轮对组装的过盈装配为例进行了有限元仿真计算，比较了两种方法的计算结果，分析了过盈量、摩擦系数、形状误差对装配应力的影响，结果对于确定合理过盈量和改进加工工艺具有参考意义。 关键词过盈配合接触非线性接触应力 0引言 在机械工程实际中普遍采用过盈配合来传递扭矩和轴向力，例如轴承配合、轴瓦配合、铁道车辆的轮轴、制动盘等。它是利用过盈量产生半径方向的接触面压力，并依靠由该面压力产生的摩擦力来传递扭矩和轴向力。由于过盈配合两个相配合的接触面上不能粘贴应变片，因此难以对其应力状态进行测定，对整个组装过程的应力状态更难以进行跟踪研究，而且这种配合方式往往承受着交变载荷的作用，配合面间可能发生相对滑动，这一滑动是随着应力变化而变化的，因而配合面边缘的接触状态和应力状态也随着应力的交变而变化，表现出复杂的状态，因此一般只能凭经验确定采用的过盈量。从力学角度看，这类问题属于接触非线性问题，传统的弹性接触解法已难以处理，可采用光弹性模拟实验进行研究，但只能反映应力分布趋势。近年来，随着非线性理论的不断完善和计算机技术的飞速发展，利用非线性有限元法来分析这类问题已日趋成熟。 铁道车辆随着向高速、重载不断发展，对轮轴的安全性要求也越来越高。研究表明，轮轴配合部位的应力状态对车轴的疲劳强度具有重要的影响，因此对轮对配合部位的宏观接触应力状态进行研究将有助于指导轮对制造标准的制定、高速重载轮对的设计和加工工艺的改进，以提高轮对的抗疲劳性能。 本文利用著名非线性有限元软件MARC，针对过盈配合的压力压装法和温差组装法对这类问题提出动态和静态两种仿真计算方法，并以铁道车辆某高速轮对的配合为例进行了计算，对比了两种计算方法的结果，分析了过盈量、摩擦系数、形状误差等因素对装配应力的影响。
基于PANDA框架的非线性静力学有限元
基于PANDA框架的非线性静力学有限元 论文导读：基于PANDA框架。能够分析千万自由度规模的弹塑性静力学问题。非线性求解策略。形成了面向对象有限元并行计算框架PANDA。并行计算，基于PANDA框架的非线性静力学有限元。关键词：PANDA，静力学，非线性，有限元，并行计算 1 引言 特种武器结构复杂，在整个库存到靶序列（Stockpile to TargetSequence，STS）全寿命周期内要经历复杂严酷的载荷和环境条件，结构响应呈现出高度的材料非线性、边界非线性和几何非线性。为提高特种武器的设计、试验和库存维护水平，对武器结构在各种条件下响应的精细建模和分析至关重要，需要充分考虑结构的几何细节和物理内涵，所建立的有限元模型可达上千万自由度规模乃至更高，而传统的商用有限元程序由于国外对我国的出口限制，非线性有限元模型的分析规模被限制在几百万自由度以下，且计算周期较长，无法快速响应设计和维护的需要。 为了提升特种武器的工程数值模拟能力，适应不断提高的武器工程数值模拟需求，迎接和加速由现阶段小规模低效率计算向大规模高效并行计算的转变，2007年中国工程物理研究院启动了院预研重大项目“武器工程大规模并行计算框架研究及基础平台开发”。该项目在已有源码程序的基础上，通过在有限元并行计算方法方面开展研究与软件开发，初步形成了面向对象有限元并行计算框架PANDA，并基于PANDA框架初步开发了可应用于部分静力、振动、冲击和传热武器工
程问题求解的大规模有限元并行计算模拟程序。 针对特种武器研制中的非线性静力学有限元大规模精细分析需求，充分消化吸收开放源代码的程序设计思想和技巧，基于PANDA框架，开发非线性静力学有限元分析所需的单元类型、材料模型、非线性并行求解策略，集成大规模线性方程组并行求解算法，初步形成了可求解小应变、有限应变线弹性和弹塑性静力学问题的非线性静力学程序。悬臂梁弹塑性有限元分析模型达到了千万自由度规模，并行求解时间低于一小时。本文介绍了基于PANDA框架的单元类型、材料模型、非线性求解策略设计，并初步验证了非线性静力学有限元并行计算程序的计算精度和千万自由度规模分析能力。 2 基于PANDA框架的非线性静力学有限元并行计算程序设计 通过中国工程物理研究院的预研重大项目，采用面向对象、层次化、组件化的设计思想，对工程结构非结构网格有限元分析程序的基本数据结构、并行通信、求解控制等方面的共性和可重用部分进行抽象和程序实现，并集成了区域分割、解法器等服务组件，形成了面向对象有限元并行计算框架PANDA，提供经过系统规划设计的应用程序开发接口，以提供服务的形式引导应用程序的设计和实现，初步建立了结构分析有限元并行计算应用程序的集成开发环境。科技论文，并行计算。 基于PANDA框架，结构分析有限元并行计算应用程序的开发工作变得较为简单和高效，程序开发工作量大为减少。在PANDA框架既设的应用软件架构下，应用程序开发者可以将精力集中到本应用程序独