solidworks进行有限元分析的一般步骤

1.软件形式：

㈠. SolidWorks的内置形式：

◆COSMOSXpress——只有对一些具有简单载荷和支撑类型的零件的静态分析。

㈡. SolidWorks的插件形式：

◆COSMOSWorks Designer——对零件或装配体的静态分析。

◆COSMOSWorks Professional——对零件或装配体的静态、热传导、扭曲、频率、掉落测试、优化、疲劳分析。

◆COSMOSWorks Advanced Professional——在COSMOSWorks Professional的所有功能上增加了非线性和高级动力学分析。

㈢. 单独发行形式：

◆COSMOS DesignSTAR——功能与COSMOSWorks Advanced Professional相同。

2.使用FEA的一般步骤：

FEA=Finite Element Analysis——是一种工程数值分析工具，但不是唯一的数值分析工具！其它的数值分析工具还有：有限差分法、边界元法、有限体积法…

①建立数学模型——有时，需要修改CAD几何模型以满足网格划分的需要，

（即从CAD几何体→FEA几何体），共有下列三法：

▲特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。▲理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理（注：如果选中了“使用中面的壳网格”做为“网格类型”，COSMOSWorks会自动地创建曲面几何体）。

▲清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。如模型中的细长面、多重实体、移动实体及其它质量问题会造成网格划分的困难甚至无法划分网格—这时我们可以使用CAD质量检查工具（即SW菜单: Tools→Check…）来检验问题所在，另外含有非常短的边或面、小的特征也必须清除掉（小特征是指其特征尺寸相对于整个模型尺寸非常小！但如果分析的目的是找出圆角附近的应力分布，那么此时非常小的内部圆角应该被保留）。

②建立有限元模型——即FEA的预处理部分，包括五个步骤：

▲选择网格种类及定义分析类型（共有静态、热传导、频率…等八种类别）——这时将产生一个FEA算例，左侧浏览器中之算例名称之后的括号里是配置名称；

▲添加材料属性: 材料属性通常从材料库中选择，它不并考虑缺陷和表面条件等因素，与几何模型相比，它有更多的不确定性。

◇右键单击“实体文件夹”并选择“应用材料到所有”——所有零部件将被赋予相同的材料属性。

◇右键单击“实体文件夹”下的某个具体零件文件夹并选择“应用材料到所有实体”——某

个零件的所有实体(多实体)将被赋予指定的材料属性。

◇右键单击“实体文件夹”下具体零件的某个“Body”并选择“应用材料到实体”——只有该“Body”被赋予指定的材料属性。

▲施加约束：定义约束是最容易产生误差的地方。通常的误差来自于过约束模型，其后果是：结构过于刚硬并低估了实际变形量和应力值。对装配体而言，还要定义“接触/间隙”这种特殊的“约束”。约束的目的是禁止模型的刚体位移。

在COSMOSWorks中共有十种约束（不包括“接触/间隙”）。它也意味着处于指定的“点、线、面”上的全部这些节点所受到的约束。

约束符号中的箭头表示“平移”约束，而圆盘则表示“回转”约束（实体单元的每个节点仅有3个移动自由度，而壳单元有6个自由度）。

对“Solid mesh”而言，因为节点无转动自由度，所以选择“固定”和“不可移动”的效果是完全一样的。定义完约束之后，模型的空间位置就被固定下来了。此时，模型不可能再发生除弹形变形之外的位移（在FEA的静态分析中，可能存在的也只能是弹形位移），称之为“模型没有刚体位移”。

▲定义载荷：在现实中，只能大概地知道载荷的大小、分布、时间依赖关系。所以，必须在FEA分析中通过简化的假设做出近似的估计。因此，定义载荷会产生较大的建模误差（理想化误差）。

注：前面的四项统称为FEA分析的“预处理”，它们的不确定性程度从高到低依次为：约束、载荷、材料、几何模型。

▲网格划分：

a). COSMOSWorks中只有两类单元：一阶单元（草稿品质单元）和二阶单元（高品质单元）。或：实体四面体单元和三角形壳单元。这样，COSMOSWorks共有四种单元类型：一阶实体四面体单元（只有4个角节点，1个高斯点）、二阶实体四面体单元（有4个角节点和6个中间节点，共计10个节点，4个高斯点）、一阶三角形壳单元（只有3个角节点，1个高斯点）、二阶三角形壳单元（有3个角节点和3个中间节点，共计6个节点，3个高斯点）——这里的四面体不一定是正四面体，而三角形也不一定是正三角形。此外，二阶单元的边和面都可以是曲线形状，以模拟单元因加载而变形的实际情形。

b).单元的品质可通过SW菜单: COSMOSWorks→Options…→选Mesh标签…

c).一般FEA中拥有最少节点的单元是横梁单元，它只有2个节点（即梁的两个端点），但每个节点处均有6个自由度（即三个平移分量加三个转动位移分量）。

d). 二阶实体四面体单元和二阶三角形壳单元适用于曲线形的几何体。

e).某些类型的形状既可以使用实体单元也可以使用壳单元，具体选用什么类型的单元取决于分析的目的。然而，通常情况下，几何体的天然形状决定了所使用的单元类型，比如，一些铸件只能用实体网格划分，而一张金属板材最好使用壳单元。

f).有限单元网格中的自由度是指单元节点的自由度。实体单元的每个节点有三个自由度（三个平移分量），壳单元的每个节点有六个自由度（三个平移分量加三个转动位移分量）。节点的位移即为这些分量的几何合成矢量。

g).在进行网格划分时，单元在匹配几何体的过程中会经历变形扭曲，但过度的扭曲会导致单元的恶化，从而导致计算量徒增和计算精度大大地降低，甚至会无法计算。为此，需要通过控制默认单元的大小（即SW菜单: COSMOSWorks→Mesh→Create…，其中：Coarse对应大，Fine对应小）或应用局部网格控制（即SW菜单: COSMOSWorks→Mesh→Apply Control…）来避免单元的过度扭曲。

h).网格质量保证：包括长宽比检查和Jacobian检查, 这些检查由程序自动执行。

长宽比检查：正四面体的长宽比通常被用做计算其它单元的长宽比。一个单元的长宽比定义

为：四面体的最长边的长度值/四面体的顶点到其相对面的法向距离的最小长度值。这里，顶点的相对面需用正四面体正则化，并假定四面体的4个角点之间用直线相连。非常小的正四面体单元的长宽比可近似地认为是1.0。作为长宽比检查的一部分，COSMOSWorks还自动执行边长检查、内切圆和外接圆检查，以及法向长度检查。

Jacobian检查：即检查雅可比行列式的值，用于判断单元的弯曲程度。一个极端扭曲单元的雅可比行列式是负值，而负的雅可比行列式会导致FEA程序的终止。Jacobian检查是基于一系列点（高斯点或节点），这些点位于每个单元中。通常情况下，雅可比率小于或等于40是可以接受的。COSMOSWorks会自动调整扭曲单元的中节点位置，以确保所有的单元能通过雅可比检查。在二次单元中，单元边界上的中节点放置在真实的几何体上；但在尖劈和弯曲边界，将中节点放置在真实几何体上会导致产生边缘下相互重叠的扭曲单元。对正四面体而言，所有中节点均精确地定位在直边中点，其雅可比率为1.0，随着边缘曲率的增加，其雅可比率也增大。Jacobian检查设置可通过COSMOSWorks→options…→Mesh标签来实现。

i).局部网格控制：由三个参数来控制——所选实体的单元尺寸、层与层之间的单元尺寸比、受局部优化影响的单元层数。它们的缺省值分别为2.2、1.5、3。网格控制可用在点（顶点）、线（边界）、面（表面）、及装配体组件上。三个控制参数可通过命令：COSMOSWorks →Mesh→Apply Control…来实现。为了找出仍在工作的最大单元，可勾选COSMOSWorks→options…→Mesh标签中的Automatic Looping选项，“自动为实体循环”功能要求网格划分程序利用更小的全局单元尺寸网格对模型进行重新划分，用户可以控制：循环实验的最大次数、全局单元尺寸每次减少的幅度、公差。

对于应用于组件的网格控制由“Component significance (零件有效数)”来定义，对于不同的Slide位置，指示网格划分程序选用不同的单元尺寸来对每个选定的组件进行网格划分。但如果“use same element size”已勾选，那么所有组件均按“网格控制”窗口中指定相同单元尺寸来进行划分。

j).实际的网格划分过程，共分三个步骤：

第一步，评估几何模型——检查CAD几何体有无缺陷；

第二步，处理边界——即先将节点置于边界上，这一步被称做表面划分；

第三步，创建网格——用四面体单元来填充实体体积。

k).如果第一步失败，则最有可能的是几何模型错误，为了验证几何模型是否错误，以IGES 输出模型，观察是否出现错误信息“处理修整的表面实体失败”。

l). 如果第二步失败，分两种情况：i.在进度指示条到达最右端之前出现错误，则说明至少在一个面上的划分出现错误，此时，右键单击网格，选择“失败诊断”，以找出有问题的表面，再有分割线或网格控制来帮助划分该表面；ii.在进度指示条到达最右端之后且在第三步开始之前出现错误，此时，需要将公差从5%（默认）到10%对单元尺寸进行增加后重新划分网格，但如果公差为10%时仍旧失败，则可以继续增加公差，但最大不要超过25%。设置命令为：COSMOSWorks→Mesh→Create…→…

m). 如果第三步失败，则表明错误发生在体积填充阶段。此时，可将单元尺寸公差从5%减少到1%，如果仍然失败，则可以25%的幅度减少单元尺寸，并设公差为1%.

n). “失败诊断”工具只对实体单元有效，对壳单元不起作用。

o). 从2008版开始，COSMOSWorks实现了自动“局部网格控制”，因而“网格划分”完全不再需要人工干预。

③求解有限元模型——在结构分析中，FEA首先计算的是网格中每个节点的位移（矢量），再在此基础上计算应变和应力等其它物理量；在热分析中，FEA首先计算的是网格中每个节点的温度（标量），再在此基础上计算温度梯度和热流等其它物理量.

一般如果模型可划分网格，那么它就可以求解，但如果没有定义材料或载荷，则求解会

终止。解算器也可检查出由于约束不足而引起的刚体运动。但刚体运动可用解算器选项来处理，比如，使用软弹簧来稳定模型，或使用平面内作用、惯性卸除。影响选择合适的解算器的五个因素：

1).问题的大小——通常，FFEPlus在处理自由度(DOF)超过100,000时，速度比较快。FFEPlus随着问题的变大会变得更有效率。

2).计算机资源——在计算机可用的内存足够多时，Direct Sparse解算器的速度比较快。

3).分析选项；

4).单元类型；

5).材料属性——当模型中使用的材料弹性模量差异较大时（比如钢和尼龙），FFEPlus（迭代法）求解比Direct Sparse（直接法）求解的精度低。

如果不能确定选择哪个解算器是分析的最佳选择时，可将解算器的类型设为“自动”。

选择求解器的命令为：COSMOSWorks→Options…→选Results标签.

④结果分析——对结果的正确解释需要我们熟悉理解：i).各种假设，如在静态分析中的材料线性假设、小变形假设、静态载荷假设；ii).简化约定；iii).前面三步中产生的误差，如建模误差（也称理想化误差）、离散误差（划分网格时产生的误差）、数值误差（求解过程中产生的误差）。在这三种误差当中，只有离散化误差是FEA特有的，故只有这个误差能够在使用FEA时被控制——网格单元越小，离散误差越低；影响数学几何模型的建模误差，是在FEA之前引入的，故只能通过正确的建模技术来控制；数值误差（求解误差）是在计算过程中产生的，难于控制，但它们通常比较小。

执行“COSMOSWorks→Options…→Results标签→Automatic Results Plots按钮”，可确定要在程序界面中显示的计算项目的结果。

结果中的波节应力（Node values）是指单元节点上的应力，而单元应力（Element values）则是指单元高斯点上的应力。

单元应力和波节应力一般是不同的，但若两者相差过大，则说明网格划分不够精细。

解析解（用数学公式求出的解）只有在平面应力假设下，板的厚度非常薄时才有效——因为它不考虑应力沿板厚方向的分布（梯度分布：中间最大，两边缘最小），认为板厚方向的截面上的应力处处相等。所以，FEA解能够比较真实地反映应力的实际状况。

abaqus有限元分析过程

一、有限单元法的基本原理 有限单元法（The Finite Element Method）简称有限元（FEM），它是利用电子计算机进行的一种数值分析方法。它在工程技术领域中的应用十分广泛，几乎所有的弹塑性结构静力学和动力学问题都可用它求得满意的数值结果。 有限元方法的基本思路是：化整为零，积零为整。即应用有限元法求解任意连续体时，应把连续的求解区域分割成有限个单元，并在每个单元上指定有限个结点，假设一个简单的函数（称插值函数）近似地表示其位移分布规律，再利用弹塑性理论中的变分原理或其他方法，建立单元结点的力和位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程组，从而求解结点的位移分量. 进而利用插值函数确定单元集合体上的场函数。由位移求出应变, 由应变求出应力 二、ABAQUS有限元分析过程 有限元分析过程可以分为以下几个阶段 1.建模阶段: 建模阶段是根据结构实际形状和实际工况条件建立有限元分析的计算模型――有限元模型，从而为有限元数值计算提供必要的输入数据。有限元建模的中心任务是结构离散，即划分网格。但是还是要处理许多与之相关的工作：如结构形式处理、集合模型建立、单元特性定义、单元质量检查、编号顺序以及模型边界条件的定义等。

2.计算阶段:计算阶段的任务是完成有限元方法有关的数值计算。 由于这一步运算量非常大，所以这部分工作由有限元分析软件控制并在计算机上自动完成 3.后处理阶段: 它的任务是对计算输出的结果惊醒必要的处理， 并按一定方式显示或打印出来，以便对结构性能的好坏或设计的合理性进行评估，并作为相应的改进或优化，这是惊醒结构有限元分析的目的所在。 下列的功能模块在ABAQUS/CAE操作整个过程中常常见到，这个表简明地描述了建立模型过程中要调用的每个功能模块。 “Part（部件） 用户在Part模块里生成单个部件，可以直接在ABAQUS/CAE环境下用图形工具生成部件的几何形状，也可以从其它的图形软件输入部件。 Property（特性） 截面（Section）的定义包括了部件特性或部件区域类信息，如区域的相关材料定义和横截面形状信息。在Property模块中，用户生成截面和材料定义，并把它们赋于(Assign)部件。 Assembly（装配件） 所生成的部件存在于自己的坐标系里，独立于模型中的其它部件。用户可使用Assembly模块生成部件的副本（instance），并且在整体坐标里把各部件的副本相互定位，从而生成一个装配件。 一个ABAQUS模型只包含一个装配件。

精讲solidworks有限元分析步骤

2013-08-29 17:31 by:有限元来源:广州有道有限元 1. 软件形式： ㈠. SolidWorks的内置形式： ◆COSMOSXpress——只有对一些具有简单载荷和支撑类型的零件的静态分析。 ㈡. SolidWorks的插件形式： ◆COSMOSWorks Designer——对零件或装配体的静态分析。 ◆COSMOSWorks Professional——对零件或装配体的静态、热传导、扭曲、频率、掉落测试、优化、疲劳分析。 ◆COSMOSWorks Advanced Professional——在COSMOSWorks Professional的所有功能上增加了非线性和高级动力学分析。 ㈢. 单独发行形式： ◆COSMOS DesignSTAR——功能与COSMOSWorks Advanced Professional相同。 2. 使用FEA的一般步骤： FEA=Finite Element Analysis——是一种工程数值分析工具，但不是唯一的数值分析工具！其它的数值分析工具还有：有限差分法、边界元法、有限体积法… ①建立数学模型——有时，需要修改CAD几何模型以满足网格划分的需要， （即从CAD几何体→FEA几何体），共有下列三法： ▲特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。▲理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理（注：如果选中了“使用中面的壳网格”做为“网格类型”，COSMOSWorks会自动地创建曲面几何体）。 ▲清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。如模型中的细长面、多重实体、移动实体及其它质量问题会造成网格划分的困难甚至无法划分网格—这时我们可以使用CAD质量检查工具（即SW菜单: Tools→Check…）来检验问题所在，另外含有非常短的边或面、小的特征也必须清除掉（小特征是指其特征尺寸相对于整个模型尺寸非常小！但如果分析的目的是找出圆角附近的应力分布，那么此时非常小的内部圆角应该被保留）。 ②建立有限元模型——即FEA的预处理部分，包括五个步骤： ▲选择网格种类及定义分析类型（共有静态、热传导、频率…等八种类别）——这时将产生一个FEA算例，左侧浏览器中之算例名称之后的括号里是配置名称； ▲添加材料属性: 材料属性通常从材料库中选择，它不并考虑缺陷和表面条件等因素，与几何模型相比，它有更多的不确定性。

有限元法的基本思想及计算 步骤

有限元法的基本思想及计算步骤 有限元法是把要分析的连续体假想地分割成有限个单元所组成的组合体，简称离散化。这些单元仅在顶角处相互联接，称这些联接点为结点。离散化的组合体与真实弹性体的区别在于：组合体中单元与单元之间的联接除了结点之外再无任何关联。但是这种联接要满足变形协调条件，即不能出现裂缝，也不允许发生重叠。显然，单元之间只能通过结点来传递内力。通过结点来传递的内力称为结点力，作用在结点上的荷载称为结点荷载。当连续体受到外力作用发生变形时，组成它的各个单元也将发生变形，因而各个结点要产生不同程度的位移，这种位移称为结点位移。在有限元中，常以结点位移作为基本未知量。并对每个单元根据分块近似的思想，假设一个简单的函数近似地表示单元内位移的分布规律，再利用力学理论中的变分原理或其他方法，建立结点力与位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程，从而求解结点的位移分量。然后利用插值函数确定单元集合体上的场函数。显然，如果单元满足问题的收敛性要求，那么随着缩小单元的尺寸，增加求解区域内单元的数目，解的近似程度将不断改进，近似解最终将收敛于精确解。 用有限元法求解问题的计算步骤比较繁多，其中最主要的计算步骤为： 1）连续体离散化。首先，应根据连续体的形状选择最能完满地描述连续体形状的单元。常见的单元有：杆单元，梁单元，三角形单元，矩形单元，四边形单元，曲边四边形单元，四面体单元，六面体单元以及曲面六面体单元等等。其次，进行单元划分，单元划分完毕后，要将全部单元和结点按一定顺序编号，每个单元所受的荷载均按静力等效原理移植到结点上，并在位移受约束的结点上根据实际情况设置约束条件。 2）单元分析。所谓单元分析，就是建立各个单元的结点位移和结点力之间的关系式。现以三角形单元为例说明单元分析的过程。如图1所示，三角形有三个结点i，j，m。在平面问题中每个结点有两个位移分量u，v和两个结点力分量F x，F y。三个结点共六个结点位移分量可用列

solidworks进行有限元分析的一般步骤

1.软件形式： ㈠. SolidWorks的内置形式： ◆COSMOSXpress——只有对一些具有简单载荷和支撑类型的零件的静态分析。 ㈡. SolidWorks的插件形式： ◆COSMOSWorks Designer——对零件或装配体的静态分析。 ◆COSMOSWorks Professional——对零件或装配体的静态、热传导、扭曲、频率、掉落测试、优化、疲劳分析。 ◆COSMOSWorks Advanced Professional——在COSMOSWorks Professional的所有功能上增加了非线性和高级动力学分析。 ㈢. 单独发行形式： ◆COSMOS DesignSTAR——功能与COSMOSWorks Advanced Professional相同。 2.使用FEA的一般步骤： FEA=Finite Element Analysis——是一种工程数值分析工具，但不是唯一的数值分析工具！其它的数值分析工具还有：有限差分法、边界元法、有限体积法… ①建立数学模型——有时，需要修改CAD几何模型以满足网格划分的需要， （即从CAD几何体→FEA几何体），共有下列三法： ▲特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。▲理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理（注：如果选中了“使用中面的壳网格”做为“网格类型”，COSMOSWorks会自动地创建曲面几何体）。▲清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。如模型中的细长面、多重实体、移动实体及其它质量问题会造成网格划分的困难甚至无法划分网格—这时我们可以使用CAD质量检查工具（即SW菜单: Tools→Check…）来检验问题所在，另外含有非常短的边或面、小的特征也必须清除掉（小特征是指其特征尺寸相对于整个模型尺寸非常小！但如果分析的目的是找出圆角附近的应力分布，那么此时非常小的内部圆角应该被保留）。 ②建立有限元模型——即FEA的预处理部分，包括五个步骤： ▲选择网格种类及定义分析类型（共有静态、热传导、频率…等八种类别）——这时将产生一个FEA算例，左侧浏览器中之算例名称之后的括号里是配置名称； ▲添加材料属性: 材料属性通常从材料库中选择，它不并考虑缺陷和表面条件等因素，与几何模型相比，它有更多的不确定性。 ◇右键单击“实体文件夹”并选择“应用材料到所有”——所有零部件将被赋予相同的材料属性。 ◇右键单击“实体文件夹”下的某个具体零件文件夹并选择“应用材料到所有实体”——某个零件的所有实体(多实体)将被赋予指定的材料属性。 ◇右键单击“实体文件夹”下具体零件的某个“Body”并选择“应用材料到实体”——只有

ANSYS 有限元分析基本流程

第一章实体建模 第一节基本知识 建模在ANSYS系统中包括广义与狭义两层含义，广义模型包括实体模型和在载荷与边界条件下的有限元模型，狭义则仅仅指建立的实体模型与有限元模型。建模的最终目的是获得正确的有限元网格模型，保证网格具有合理的单元形状，单元大小密度分布合理，以便施加边界条件和载荷，保证变形后仍具有合理的单元形状，场量分布描述清晰等。 一、实体造型简介 1．建立实体模型的两种途径 ①利用ANSYS自带的实体建模功能创建实体建模： ②利用ANSYS与其他软件接口导入其他二维或三维软件所建立的实体模型。 2．实体建模的三种方式 (1)自底向上的实体建模 由建立最低图元对象的点到最高图元对象的体，即先定义实体各顶点的关键点，再通过关键点连成线，然后由线组合成面，最后由面组合成体。 (2)自顶向下的实体建模 直接建立最高图元对象，其对应的较低图元面、线和关键点同时被创建。 (3)混合法自底向上和自顶向下的实体建模 可根据个人习惯采用混合法建模，但应该考虑要获得什么样的有限元模型，即在网格划分时采用自由网格划分或映射网格划分。自由网格划分时，实体模型的建立比较1e单，只要所有的面或体能接合成一体就可以：映射网格划分时，平面结构一定要四边形或三边形的面相接而成。 二、ANSYS的坐标系 ANSYS为用户提供了以下几种坐标系，每种都有其特定的用途。 ①全局坐标系与局部坐标系：用于定位几何对象(如节点、关键点等)的空间位置。 ②显示坐标系：定义了列出或显示几何对象的系统。 ③节点坐标系：定义每个节点的自由度方向和节点结果数据的方向。 ④单元坐标系：确定材料特性主轴和单元结果数据的方向。 1．全局坐标系 全局坐标系和局部坐标系是用来定位几何体。在默认状态下，建模操作时使用的坐标系是全局坐标系即笛卡尔坐标系。总体坐标系是一个绝对的参考系。ANSYS提供了4种全局坐标系：笛卡尔坐标系、柱坐标系、球坐标系、Y-柱坐标系。4种全局坐标系有相同的原点，且遵循右手定则，它们的坐标系识别号分别为：0是笛卡尔坐标系(cartesian)，1是柱坐标系 (Cyliadrical)，2是球坐标系(Spherical)，5是Y-柱坐标系(Y-aylindrical)，如图2-1所示。

有限元法分析过程

有限元法分析过程 有限元法分析过程大体可分为：前处理、分析、后处理三大步骤。 对实际的连续体经过离散化后就建立了有限元分析模型，这一过程是有限元的前处理过程。在这一阶段，要构造计算对象的几何模型，要划分有限元网格，要生成有限元分析的输入数据，这一步是有限元分析的关键。 有限元分析过程主要包括：单元分析、整体分析、载荷移置、引入约束、求解约束方程等过程。这一过程是有限元分析的核心部分，有限元理论主要体现在这一过程中。 有限元法包括三类：有限元位移法、有限元力法、有限元混合法。 在有限元位移法中，选节点位移作为基本未知量； 在有限元力法中，选节点力作为未知量； 在有限元混合法中，选一部分基本未知量为节点位移，另一部分基本未知量为节点力。 有限元位移法计算过程的系统性、规律性强，特别适宜于编程求解。一般除板壳问题的有限元应用一定量的混合法外，其余全部采用有限元位移法。因此，一般不做特别声明，有限元法指的是有限元位移法。 有限元分析的后处理主要包括对计算结果的加工处理、编辑组织和图形表示三个方面。它可以把有限元分析得到的数据，进一步转换为设计人员直接需要的信息，如应力分布状态、结构变形状态等，并且绘成直观的图形，从而帮助设计人员迅速的评价和校核设计方案。 附：FELAC 2.0软件简介 FELAC 2.0采用自定义的有限元语言作为脚本代码语言，它可以使用户以一种类似于数学公式书写和推导的方式，非常自然和简单的表达待解问题的微分方程表达式和算法表达式，并由生成器解释产生完整的并行有限元计算C程序。 FELAC 2.0的目标是通过输入微分方程表达式和算法之后，就可以得到所有有限元计算的程序代码，包含串行程序和并行程序。该系统采用一种语言(有限元语言)和四种技术(对象技术、组件技术、公式库技术生成器技术)开发而成。并且基于FELAC 1.0的用户界面，新版本扩充了工作目录中右键编译功能、命令终端输入功能，并且丰

solidworks有限元分析范例

注意:本文件内容只是一个简短的分析报告样板，其内相关的分析条件、设置和结果不一定是正确的，您还是要按本书正文所教的自行来做。 一、范例名： (Gas Valve气压阀) 1 设计要求： （1）输入转速1500rpm。 （2）额定输出压力5Mpa，最大压力10Mpa。 2 分析零件 该气压泵装置中，推杆活塞、凸轮轴和箱体三个零件是主要的受力零件，因此对这三个零件进行结构分析。 3 分析目的 （1）验证零件在给定的载荷下静强度是否满足要求。 （2）分析凸轮轴零件和推杆活塞零件的模态，在工作过程中避开共振频率。 （3）计算凸轮轴零件的工作寿命。 4 分析结果 1.。推杆活塞零件 材料：普通碳钢。 在模型上直接测量得活塞推杆的受力面积S为：162mm2，由F=PS计算得该零件端面的力F为：1620N。所得结果包括： 1 静力计算： （1）应力。如图1-1所示，由应力云图可知，最大应力为21Mpa，静强度设计符合要求。 （2）位移。如图1-2所示，零件变形导致的最大静位移为2.2e-6m。 （3）应变。如图1-3所示，应变云图与应力云图的对应的，二者之间存在一转换关系。

图1-1 应力云图图1-2 位移云图 图1-3 应变云图图1-4 模态分析 2 模态分析： 图1-4的“列举模式”对话框中列出了“推杆活塞”零件在工作载荷下，其前三阶的模态的频率远远大于输入转速的频率，因此在启动及工作过程中，该零件不会发生共振情况。模态验证符合设计要求。 2。凸轮轴零件 材料：45钢，屈服强度355MPa。 根据活塞推杆的受力情况，换算至该零件上的扭矩约为10.5N·m。 1 静力分析： 如图1-5所示为“凸轮轴”零件的应力云图，零件上的最大应力为212Mpa，平均应力约为120MPa，零件的安全系数约为1.7，符合设计要求。 图1-5 应力云图图1-6 模态分析 2 模态分析

Matlab有限元分析操作基础共11页

Matlab有限元分析20140226 为了用Matlab进行有限元分析，首先要学会Matlab基本操作，还要学会使用Matlab进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵

2. Matlab有限元分析的基本操作 （1）单元划分（选择何种单元，分成多少个单元，标号）（2）构造单元刚度矩阵（列出…） （3）组装系统刚度矩阵（集成整体刚度矩阵） （4）引入边界条件（消除冗余方程） （5）解方程 （6）后处理（扩展计算）

3. Matlab有限元分析实战【实例1】

分析： 步骤一：单元划分

>>k1=SpringElementStiffness(100)

a) 分析SpringAssemble库函数 function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？ 今天的系统刚度矩阵是什么？ 因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ????-- ??？

有限元分析及其应用思考题附答案2012

有限元分析及其应用-2010 思考题： 1、有限元法的基本思想是什么？有限元法的基本步骤有那些？其中“离散”的含义是什 么？是如何将无限自由度问题转化为有限自由度问题的？ 答：基本思想：几何离散和分片插值。 基本步骤：结构离散、单元分析和整体分析。 离散的含义：用假想的线或面将连续物体分割成由有限个单元组成的集合，且单元之间仅在节点处连接，单元之间的作用仅由节点传递。当单元趋近无限小，节点无限多，则这种离散结构将趋近于实际的连续结构。 2、有限元法与经典的差分法、里兹法有何区别？ 区别：差分法：均匀离散求解域，差分代替微分，要求规则边界，几何形状复杂精度较低； 里兹法：根据描述问题的微分方程和相应的定解构造等价的泛函表达式，求得近似解； 有限元：基于变分法，采用分片近似进而逼近总体的求解微分方程的数值计算方法。 3、一根单位长度重量为q的悬挂直杆，上端固定，下端受垂直向下的外力P，试 1）建立其受拉伸的微分方程及边界条件； 2）构造其泛函形式； 3）基于有限元基本思想和泛函求极值构造其有限元的计算格式（即最小势能原理）。4、以简单实例为对象，分别按虚功原理和变分原理导出有限元法的基本格式（单元刚度矩 阵）。 5、什么是节点力和节点载荷？两者有何区别？ 答：节点力：单元与单元之间通过节点相互作用 节点载荷：作用于节点上的外载 6、单元刚度矩阵和整体刚度矩阵各有何特点？其中每个矩阵元素的物理意义是什么（按自 由度和节点解释）？ 答：单元刚度矩阵：对称性、奇异性、主对角线恒为正 整体刚度矩阵：对称性、奇异性、主对角线恒为正、稀疏性、带状性。 Kij，表示j节点产生单位位移、其他节点位移为零时作用i节点的力，节点力等于节点位移与单元刚度元素乘积之和。 7、单元的形函数具有什么特点？有哪些性质？ 答：形函数的特点：Ni为x，y的坐标函数，与位移函数有相同的阶次。 形函数Ni在i节点的值为1，而在其他节点上的值为0； 单元内任一点的形函数之和恒等于1； 形函数的值在0~1间变化。 8、描述弹性体的基本变量是什么？基本方程有哪些组成？ 答：基本变量：外力、应力、应变、位移 基本方程：平衡方程、几何方程、物理方程、几何条件 9、何谓应力、应变、位移的概念？应力与强度是什么关系？ 答：应力：lim△Q/△A=S △A→0 应变：物体形状的改变 位移：弹性体内质点位置的变化 10、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系？何谓“强形 式”？何谓“弱形式”，两者有何区别？建立弱形式的关键步骤是什么？

solidworks有限元分析的分析方法

solidworks有限元分析的分析方法 solidworks有限元分析可应用于机械、汽车、家电、电子产品、家具、建筑、医学骨科等产品设计及研发。其作用是：确保产品设计的安全合理性，同时采用优化设计，找出产品设计最佳方案，降低材料的消耗或成本; 在产品制造或工程施工前预先发现潜在的问题; 模拟各种试验方案，减少试验时间和经费; 是产品设计研发的核心技术， 学习solidworks有限元分析最重要的是随机应变而不是根据理论一成不变的做。这是看板网经过超过十几年的实践经验和培训经验中总结出来的教训。这也是初学者常常出现的问题，过度重视理论，行为上变现为按部就班，结果往往出现问题。按部就班对于刚开始学习的人是有好处的，但是，学习一段时间后就应该学会创新。 第一步，要知道零部件之间的接触关系。 一般来讲，有限元分析的前要有符合实际的模型，再有符合实际的约束条件，如果是装配体需要知道零部件之间的接触关系。 第二步，建立有限元模型 在SW的有限元分析中可以用非常简单的步骤做到，知道零部件之间的接触关系。首先建立有限元模型，为什么叫有限元模型？因为为了节约分析的时间，降低有些应力集中区域的产生，我们需要对模型简化，所以，一般来讲有限元模型和实际的模型有一点的区别。 第三步，启动有限元分析插件 模型建好后就可以启动有限元分析插件，具体：点插件---Simulation框内打勾，启动后就可以进行边界条件的设置，然后点运行就可以了，当然，如果不设定网格精度，软件会默认网格的大小。 第四步，后处理 关于后处理：前面是i过程，而后处理就是我们要的结果，这个需要你的一些力学上的知识支撑了，比如应力、应变、位移、安全系数、作用力、反作用力等等。具体的还需要对软件进一步的了解！ 看板网培训的主要Simulation有限元分析理论分类有静态分析、频率分析、扭曲分析、热分析、设计优化、掉落测试、疲劳理论、疲劳分析、振动分析。

有限元分析的一般过程

一、结构的离散化 将结构或弹性体人为地划分成由有限个单元，并通过有限个节点相互连接的离散系统。 这一步要解决以下几个方面的问题: 1、选择一个适当的参考系，既要考虑到工程设计习惯，又要照顾到建立模型的方便。 2、根据结构的特点，选择不同类型的单元。对复合结构可能同时用到多种类型的单元，此时还需要考虑不同类型单元的连接处理等问题。 3、根据计算分析的精度、周期及费用等方面的要求，合理确定单元的尺寸和阶次。 4、根据工程需要，确定分析类型和计算工况。要考虑参数区间及确定最危险工况等问题。 5、根据结构的实际支撑情况及受载状态，确定各工况的边界约束和有效计算载荷。 二、选择位移插值函数 1、位移插值函数的要求 在有限元法中通常选择多项式函数作为单元位移插值函数，并利用节点处的位移连续性条件，将位移插值函数整理成以下形函数矩阵与单元节点位移向量的乘积形式。 位移插值函数需要满足相容（协调）条件，采用多项式形式的位移插值函数，这一条件始终可以满足。 但近年来有人提出了一些新的位移插值函数，如：三角函数、样条函数及双曲函数等，此时需要检查是否满足相容条件。 2、位移插值函数的收敛性（完备性）要求： 1）位移插值函数必须包含常应变状态。 2）位移插值函数必须包含刚体位移。 3、复杂单元形函数的构造 对于高阶复杂单元，利用节点处的位移连续性条件求解形函数，实际上是不可行的。因此在实际应用中更多的情况下是利用形函数的性质来构造形函数。 形函数的性质： 1）相关节点处的值为 1，不相关节点处的值为 0。 2）形函数之和恒等于 1。 1、建立数学模型（特征消隐，理想化，清除）（（即从CAD 几何体→FEA 几何体），共 有下列三法：▲ 特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。▲ 理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理▲ 清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。） 2、建立有限元模型：（选择网格种类及定义分析类型；添加材料属性；施加约束；定义载 荷；网格划分） 3、求解有限元模型：再在此基础上计算应变和应力等其它物理量；在热分析中，FEA 首先 计算的是网格中每个节点的温度（标量），再在此基础上计算温度梯度和热流等其它物理量. 一般如果模型可划分网格，那么它就可以求解，但如果没有定义材料或载荷，则求解会终止。 4、结果分析：材料线性假设、小变形假设、静态载荷假设等等。

有限元分析基础教程(ANSYS算例)(曾攀)

有限元分析基础教程Fundamentals of Finite Element Analysis (ANSYS算例) 曾攀 清华大学 2008-12

有限元分析基础教程曾攀 有限元分析基础教程 Fundamentals of Finite Element Analysis 曾攀 (清华大学) 内容简介 全教程包括两大部分，共分9章；第一部分为有限元分析基本原理，包括第1章至第5章，内容有：绪论、有限元分析过程的概要、杆梁结构分析的有限元方法、连续体结构分析的有限元方法、有限元分析中的若干问题讨论；第二部分为有限元分析的典型应用领域，包括第6章至第9章，内容有：静力结构的有限元分析、结构振动的有限元分析、传热过程的有限元分析、弹塑性材料的有限元分析。本书以基本变量、基本方程、求解原理、单元构建、典型例题、MATLAB程序及算例、ANSYS算例等一系列规范性方式来描述有限元分析的力学原理、程序编制以及实例应用；给出的典型实例都详细提供有完整的数学推演过程以及ANSYS实现过程。本教程的基本理论阐述简明扼要，重点突出，实例丰富，教程中的二部分内容相互衔接，也可独立使用，适合于具有大学高年级学生程度的人员作为培训教材，也适合于不同程度的读者进行自学；对于希望在MATLAB程序以及ANSYS平台进行建模分析的读者，本教程更值得参考。 本基础教程的读者对象：机械、力学、土木、水利、航空航天等专业的工程技术人员、科研工作者。

目录 [[[[[[\\\\\\ 【ANSYS算例】3.3.7(3) 三梁平面框架结构的有限元分析 1 【ANSYS算例】4.3.2(4) 三角形单元与矩形单元的精细网格的计算比较 3 【ANSYS算例】5.3(8) 平面问题斜支座的处理 6 【ANSYS算例】6.2(2) 受均匀载荷方形板的有限元分析9 【ANSYS算例】6.4.2(1) 8万吨模锻液压机主牌坊的分析(GUI) 15 【ANSYS算例】6.4.2(2) 8万吨模锻液压机主牌坊的参数化建模与分析(命令流) 17 【ANSYS算例】7.2(1) 汽车悬挂系统的振动模态分析(GUI) 20 【ANSYS算例】7.2(2) 汽车悬挂系统的振动模态分析(命令流) 23 【ANSYS算例】7.3(1) 带有张拉的绳索的振动模态分析(GUI) 24 【ANSYS算例】7.3(2) 带有张拉的绳索的振动模态分析(命令流) 27 【ANSYS算例】7.4(1) 机翼模型的振动模态分析(GUI) 28 【ANSYS算例】7.4(2) 机翼模型的振动模态分析(命令流) 30 【ANSYS算例】8.2(1) 2D矩形板的稳态热对流的自适应分析(GUI) 31 【ANSYS算例】8.2(2) 2D矩形板的稳态热对流的自适应分析(命令流) 33 【ANSYS算例】8.3(1) 金属材料凝固过程的瞬态传热分析(GUI) 34 【ANSYS算例】8.3(2) 金属材料凝固过程的瞬态传热分析(命令流) 38 【ANSYS算例】8.4(1) 升温条件下杆件支撑结构的热应力分析(GUI) 39 【ANSYS算例】8.4(2) 升温条件下杆件支撑结构的热应力分析(命令流) 42 【ANSYS算例】9.2(2) 三杆结构塑性卸载后的残余应力计算(命令流) 45 【ANSYS算例】9.3(1) 悬臂梁在循环加载作用下的弹塑性计算(GUI) 46 【ANSYS算例】9.3(2) 悬臂梁在循环加载作用下的弹塑性计算(命令流) 49 附录 B ANSYS软件的基本操作52 B.1 基于图形界面(GUI)的交互式操作(step by step) 53 B.2 log命令流文件的调入操作(可由GUI环境下生成log文件) 56 B.3 完全的直接命令输入方式操作56 B.4 APDL参数化编程的初步操作57

有限元分析步骤

有限元建模与分析 有限元分析（FEA）是一种预测结构的偏移与其它应力影响的过程，有限元建模（FEM）将这个结构分割成单元网格以形成实际结构的模型，每个单元具有简单形态（如正方形或三角形）。这样有限元程序就有了可写出在刚度矩阵结构中控制方程方面的信息。每个单元上的未知量就是在节点上的位移，这个点就是单元元的连接点。有限元程序将这些单个单元的刚度矩阵组合起来以形成整个模型的总刚度矩阵，并给予已知力和边界条件来求解该刚度矩阵以得出未知位移，从节点上位移的变化就可以计算出每个单元中的应力。 有限单元由假定的应变方程式导出，有些单元可假设其应变是常量，而另外一些可采用更高阶的函数。利用给定单元的这些方程和实际几何体，则可以写出外力和节点位移之间的平衡方程。对于单元的每个节点来说，每个自由度就有一个方程，这些方程被十分便利地写成矩阵的形式以用于计算机的演算中，这个系数的矩阵就变成了一个显示出力对位移的关系的刚度矩阵：{F}=[K]、{d} 尽管求知量处于离散的自由度，内部方程仍被写成表述为连续集的应变函数。这就意味着如果选择了正确单元的话，纵然这个有限元模型有一组离散的方程，只要用有限的节点和单元也可以收敛出正确的答案。 有限元模型是解决全部结构问题的完全理想的模型。这些问题包括节点的定位，单元，物理的和材料的特性，载荷和边界条件，根据分析类型的不同，如静态结构载荷，动态的或热力分析，这个模型就确定得不同。 一个有限元模型常常由不止一种单元类型来建立，有限元模型是以结构的偏移来建立成数学模型，而不只是在外观上象原结构。也许某个零件用梁单元最好，而另外的零件则可能用薄壳单元最理想。 对于给定的问题来讲，求解结果的准确性将取决于结构建模的好坏，负载和边界条件的确定，以及所用单元的精度。 一般来讲，如模型细分更小的单元，则求解将更准确。了解你在最终的求解结果上有充分收敛的唯一确信的方法是用更细网格的单元来建立更多的模型，以检查求解结果的收敛性。 新的有限元用户经常产生想象上的错误，即建立一个有限元模型的目的是建立一个看起来象这种结构的模型。有限元建模的目的是建立一个从数学意义是“相似”的模型，而不是一个外观相似的模型。一个有经验的使用者学会了怎样选择单元的正确类型，和在模型的不同区域中怎样来细分网格。 一个经常忽略的错误根源是在一个模型中的负载和边界条件上进行了错误的假设。同时也很轻易地相信一个有限元模型的每个十进位的结果。以及忘掉了在负载和边界条件上粗糙的假设。如果有一个关于怎样建立边界条件模型的问题的话，宁可用你的模型以不同的方法去测试其灵敏度，而不是仅遵循一种方法，得出一种答案，

Solidworks应力分析实例

基于Solidworks 软件的应力分析 Solidworks 中有限元分析插件CosMos/Works 分析零件的静力学性能，得出载荷分布情况，定性的分析极限载荷（这里指的是最大扭矩）下的应力，应变分布及其安全性能。 其分析流程如下： 1、建立一个简化的分析模型； 2、指定材料、元素和截面； 3、加约束和载荷； 4、设定网格； 5、执行分析； 6、结果显示； 7、生成研究报告。 分析对象 电机轴及啮合处的变速器输入轴，离合器花键轴及啮合处的离合器从动盘，电机轴和离合器花键轴之间的联接螺栓（M12x40，10.9级）。 材料 目前公司所用的变速器输入轴材料为20CrMnTi ，考虑其受力情况，材料不一致，其强度就会不一样，容易导致强度差的失效，因此根据目前情况，电机轴和离合器花键轴均选用20CrMnTi 。 20CrMnTi 用于制作渗碳零件，渗碳淬火后有良好的耐磨性和抗弯强度，有较高的低温冲击韧性，切削加工性能良好，承受高速、中载或重载以及冲击和摩擦的主要零件。 对于截面为15的样件，经过第一次淬火880℃，第二次淬火870℃，油冷；在经过回火200℃，水冷和空冷。得到的力学性能：抗拉强度MPa b 1080=σ，屈服强度MPa s 835=σ，伸长率（式样的标距等于5倍直径时的伸长率）%105=δ，断面收缩率%45=ψ，冲击韧度2/55cm J A kU =，硬度217HB 。

对于截面尺寸小于等于100的样件，经过调质处理，力学性能：抗拉强度 MPa b 615=σ，屈服强度MPa s 395=σ，伸长率%175=δ，断面收缩率%45=ψ， 冲击韧度2/47cm J A kU =。本分析还要使用到的参数：泊松比25.0=μ，抗剪模量G=7.938GPa ，弹性模量E=207GPa ，密度23/108.7m N ?=ρ。 螺栓联接受力分析 螺纹联接根据载荷性质不同，其失效形式也不同。受静载荷螺栓的失效形式多为螺纹部分的塑性变形或螺栓被拉断；受变向载荷螺栓的失效形式多为螺栓的疲劳断裂；对于受横向载荷的绞制孔用螺栓联接，其失效形式主要为螺栓杆被剪断，螺栓杆或连接孔接触面被挤压破坏。 对于10.9级M12的普通螺栓，屈服强度MPa s 900=σ，拧紧力矩T=120N.m 。 为了增强螺纹连接的刚性、防松能力及防止受载螺栓的滑动，装配时需要预紧。 其拧紧扳手力矩T 用于克服螺纹副的阻力矩T1及螺母与被连接件支撑面间的摩擦力矩T2，装配时可用力矩扳手法控制力矩。 公式： d * F *K =T2+T1=T 0 拧紧扳手力矩T=120N.m ，其中K 为拧紧力矩系数，0 F 为预紧力N ，d 为螺 纹公称直径12mm 。 摩擦表面状态 K 值 有润滑 无润滑 精加工表面 0.1 0.12 一般工表面 0.13-0.15 0.18-0.21 表面氧化 0.2 0.24 镀锌 0.18 0.22 粗加工表面 - 0.26-0.3

solidworks有限元分析使用方法

solidworks有限元分析使用方法 solidworks有限元分析应用于机械、汽车、家电、电子产品、家具、建筑、医学骨科等产品设计及研发。其作用是：确保产品设计的安全合理性，同时采用优化设计，找出产品设计最佳方案，降低材料的消耗或成本; 在产品制造或工程施工前预先发现潜在的问题; 模拟各种试验方案，减少试验时间和经费; 是产品设计研发的核心技术。看板网根据超过十年的项目经验和培训经验，提醒各位朋友，有限元分析，不同于绘图。以下是看板网总结的solidworks有限元分析使用方法，希望对大家有用。 一、软件形式： （一）solidworks的内置形式： SimulationXpress——只有对一些具有简单载荷和支撑类型的零件的静态分析。 （二）SolidWorks的插件形式： SimulationWorksDesigner——对零件或装配体的静态分析。 SimulationWorksProfessional——对零件或装配体的静态、热传导、扭曲、频率、掉落测试、优化、疲劳分析。 SimulationWorksAdvancedProfessional——在SimulationWorksProfessional的所有功能上增加了非线性和高级动力学分析。 （三）单独发行形式： Simulation DesignSTAR——功能与SimulationWorks Advanced Professional相同。 二、使用FEA的一般步骤： FEA=Finite Element Analysis——是一种工程数值分析工具，但不是唯一的数值分析工具！其它的数值分析工具还有：有限差分法、边界元法、有限体积法等等。 （一）建立数学模型 有时，需要修改CAD几何模型以满足网格划分的需要，（即从CAD几何体→FEA几何体），共有下列三法：

Matlab有限元分析操作基础

Matlab 有限元分析20140226 为了用Matlab 进行有限元分析，首先要学会Matlab 基本操作，还要学会使用Matlab 进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵 11221 2 1200k k k k k k k k -?? ??-????--+??

2. Matlab有限元分析的基本操作 （1）单元划分（选择何种单元，分成多少个单元，标号）（2）构造单元刚度矩阵（列出…） （3）组装系统刚度矩阵（集成整体刚度矩阵） （4）引入边界条件（消除冗余方程） （5）解方程 （6）后处理（扩展计算）

3. Matlab有限元分析实战【实例1】

分析： 步骤一：单元划分

步骤二：构造单元刚度矩阵 >>k1=SpringElementStiffness(100) >>…?

步骤三：构造系统刚度矩阵 a) 分析SpringAssemble库函数 function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？ 今天的系统刚度矩阵是什么？ 因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ?? ?? -- ?? ？

有限元分析过程

有限元分析过程 有限元分析过程可以分为以下三个阶段：1.建模阶段: 建模阶段是根据结构实际形状和实际工况条件建立有限元分析的计算模型——有限元模型，从而为有限元数值计算提供必要的输入数据。有限元建模的中心任务是结构离散，即划分网格。但是还是要处理许多与之相关的工作：如结构形式处理、集合模型建立、单元特性定义、单元质量检查、编号顺序以及模型边界条件的定义等。2.计算阶段: 计算阶段的任务是完成有限元方法有关的数值计算。由于这一步运算量非常大，所以这部分工作由有限元分析软件控制并在计算机上自动完成。3.后处理阶段: 它的任务是对计算输出的结果惊醒必要的处理，并按一定方式显示或打印出来，以便对结构性能的好坏或设计的合理性进行评估，并作为相应的改进或优化，这是惊醒结构有限元分析的目的所在。注意：在上述三个阶段中，建立有限元模型是整个有限分析过程的关键。首先，有限元模型为计算提供所以原始数据，这些输入数据的误差将直接决定计算结果的精度；其次，有限元模型的形式将对计算过程产生很大的影响，合理的模型既能保证计算结构的精度，又不致使计算量太大和对计算机存储容量的要求太高；再次，由于结构形状和工况条件的复杂性，要建立一个符合实际的有限元模型并非易事，它要考虑的综合因素很多，对分析人员提出了较高的要求；最后，建模所花费的时间在整个分析过程中占有相当大的比重，约占整个分析时间的70%，因此，把主要精力放在模型的建立上以及提高建模速度是缩短整个分析周期的关键。原始数据的计算模型，模型中一般包括以下三类数据：1.节点数据: 包括每个节点的编号、坐标值等；2.单元数据：a.单元编号和组成单元的节点编号；b.单元材料特性，如弹性模量、泊松比、密度等；c.单元物理特征值，如弹簧单元的刚度系数、单元厚度、曲率半径等；d.一维单元的截面特征值，如截面面积、惯性矩等；e.相关几何数据3.边界条件数据：a.位移约束数据；b.载荷条件数据；c.热边界条件数据；d.其他边界数据.建立有限元模型的一般过程：1.分析问题定义在进行有限元分析之前，首先应对结果的形状、尺寸、工况条件等进行仔细分析，只有正确掌握了分析结构的具体特征才能建立合理的几何模型。总的来说，要定义一个有限元分析问题时，应明确以下几点：a.结构类型； b.分析类型； c.分析内容； d.计算精度要求； e.模型规模； f.计算数据的大致规律2.几何模型建立几何模型是从结构实际形状中抽象出来的，并不是完全照搬结构的实际形状，而是需要根据结构的具体特征对结构进行必要的简化、变化和处理，以适应有限元分析的特点。3.单元类型选择划分网格前首先要确定采用哪种类型的单元，包括单元的形状和阶次。单元类型选择应根据结构的类型、形状特征、应力和变形特点、精度要求和硬件条件等因素综合进行考虑。4.单元特性定义有限元单元中的每一个单元除了表现出一定的外部形状外，还应具备一组计算所需的内部特征参数，这些参数用来定义结构材料的性能、描述单元本身的物理特征和其他辅助几何特征等.5.网格划分网格划分是建立有限元模型的中心工作，模型的合理性很大程度上可以通过所划分的网格形式反映出来。目前广泛采用自动或半自动网格划分方法，如在Ansys中采用的SmartSize网格划分方法就是自动划分方法。6.模型检查和处理一般来说，用自动或半自动网格划分方法划分出来的网格模型还不能立即应用于分析。由于结构和网格生成过程的复杂性，划分出来的网格或多或少存在一些问题，如网格形状较差，单元和节点编号顺序不合理等，这些都将影响有限元计算的计算精度和计算时间。7.边界条件定义在对结构进行网格划分后称为离散模型，它还不是有限元模型，只有在网格模型上定义了所需要的各类边界条件后，网格模型才能成为完整的有限元模型。计算机几何建模方法㈠.几何模型的形式1.线框模型：用组成结构的棱边表示结构形状和大小的模型称为线框模型，或线架模型。它是使用最早的几何模型，其特点是数据量少、数据结构简单、算法处理方便，模型输入可以通过定义线段端点坐标来实现。但是这种模型有很大的局限性，它的几何描述能力差，只能提供一个框架，对几何形状的理解很容易产生多义性，也不能计算结构的重量、

有限元单元法求解问题的的基本步骤

诚信·公平·开放·共赢 Loyalty Fair Opening Win-win 有限元单元法求解问题的的基本步骤 (1)建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发点。 (2) 区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作，这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外，还要表示节点的位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。 （3) 确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函数作为单元基函数。有限元方法中的基函数是在单元中选取的，由于各单元具有规则的几何形状，在选取基函数时可遵循一定的法则。 (4) 单元分析：将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近;再将近似函数代入积分方程，并对单元区域进行积分，可获得含有待定系数(即单元中各节点的参数值)的代数方程组，称为单元有限元方程。 (5) 总体合成：在得出单元有限元方程之后，将区域中所有单元有限元方程按一定法则进行累加，形成总体有限元方程。 (6) 边界条件的处理：一般边界条件有三种形式，分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。对于自然边界条件，一般在积分表达式中可自动得到满足。对于本质边界条件和混合边界条件，需按一定法则对总体有限元方程进行修正满足。 (7) 解有限元方程：根据边界条件修正的总体有限元方程组，是含所有待定未知量的封闭方程组，采用适当的数值计算方法求解，可求得各节点的函数值。