【基本建模Project1】2D问题：带孔平板的有限元分析

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。三个结点共六个结点位移分量可用列

ANSYS 有限元分析 平面薄板

《有限元基础教程》作业二：平面薄板的有限元分析 班级：机自101202班 姓名：韩晓峰 学号：201012030210 一．问题描述： P P h1mm R1mm 10m m 10mm 条件：上图所示为一个承受拉伸的正方形板，长度和宽度均为10mm ，厚度为h 为1mm ，中心圆的半径R 为1mm 。已知材料属性为弹性模量E=1MPa ，泊松比为0.3，拉伸的均布载荷 q ＝1N/mm 2。根据平板结构的对称性，只需分析其中的二分之一即可，简化模型如上右图所 示。 二．求解过程： 1 进入ANSYS 程序 →ANSYS 10.0→ANSYS Product Launcher →File management →input job name: ZY2→Run 2设置计算类型 ANSYS Main Menu: Preferences →select Structural → OK 3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK → Options… →select K3: Plane Strs w/thk →OK →Close 4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX: 1e6, PRXY:0.3 → OK 5定义实常数以及确定平面问题的厚度 A NSYS Main Menu: Preprocessor →Real Constants …→Add/Edit/Delete →Add →Type 1→OK →Real Constant Set No.1,THK:1→OK →Close 6生成几何模型 a 生成平面方板 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Rectangle →By 2 Corners →WP X:0,WP Y:0,Width:5,Height:5→OK b 生成圆孔平面 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Circle →Solid Circle →WPX=0,WPY=0,RADIUS=1→OK b 生成带孔板 ANSYS Main Menu: Preprocessor →Modeling →Operate →Booleans → Subtract →Areas →点击area1→OK →点击area2→OK 7 网格划分 A NSYS Main Menu: Preprocessor →Meshing →Mesh Tool →(Size Controls) Global: Set →SIZE: 0.5 →OK →iMesh →Pick All → Close

带孔平板拉伸作业

带孔平板有限元分析 本文采用有限元法，对带圆孔的矩形平板进行了弹塑性受力分析，分析了圆孔处的应力集中现象，为其设计和应用提供了参考依据。 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有限元网格划分 网格划分是非常重要的过程，它会对计算速度、精度、可靠性产生重要影响。网格划分主要包括两方面：尺寸、单元类型。

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所示。

带孔平板的线性静力分析

带孔平板的线性静力分析 本示例将对一个给定的带孔平板几何模型创建有限元模型、施加边界条件、进行有限元分析并在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文件，模型位于/tutorials/os/。 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 的卡片信息。 如果括号中的量下面没有值，表示其处于关闭状态。要改变该状态，点击括号中的量，

有限元法分析过程

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

带孔板的建模及有限元分析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

开孔板问题(仅限借鉴)

开孔板的问题（应力集中的问题） 一. 引言 应力集中即Stress concentration，是指受力构件由于外界因素或自身因素几何形状、外形尺寸发生突变而引起局部范围内应力显著增大的现象。 在弹性力学中，这是一类问题，应力在固体局部区域内显著增高的现象。多出现于尖角、孔洞、缺口、沟槽以及有刚性约束处及其邻域。应力集中会引起脆性材料断裂；使物体产生疲劳裂纹。在应力集中区域，应力的最大值（峰值应力）与物体的几何形状和加载方式等因素有关。局部增高的应力值随与峰值应力点的间距的增加而迅速衰减。由于峰值应力往往超过屈服极限而造成应力的重新分配，所以，实际的峰值应力常低于按弹性力学计算出的理论峰值应力。反映局部应力增高程度的参数称为应力集中系数k，它是峰值应力与不考虑应力集中时的应力的比值，恒大于1且与载荷大小无关。 二．产生应力集中的原因 构件中产生应力集中的原因主要有： (1) 截面的急剧变化。如：构件中的油孔、键槽、缺口、台阶等。 (2) 受集中力作用。如：齿轮轮齿之间的接触点，火车车轮与钢轨的接触点等。 (3) 材料本身的不连续性。如材料中的夹杂、气孔等。 (4) 构件中由于装配、焊接、冷加工、磨削等而产生的裂纹。 (5) 构件在制造或装配过程中，由于强拉伸、冷加工、热处理、焊接等而引起的残余应力。这些残余应力叠加上工作应力后，有可能出现较大的应力集中。(6) 构件在加工或运输中的意外碰伤和刮痕。 三．弹性力学中的应力集中 1.工程结构中常开设孔口最简单的为圆孔。 弹性力学研究‘小孔口问题’，应符合 (1)孔口尺寸＜＜弹性体尺寸， 孔口引起的应力扰动局限于小范围内。 (2)孔边距边界较远（＞1.5倍孔口尺寸） 孔口与边界不相互干扰。

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、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系？何谓“强形 式”？何谓“弱形式”，两者有何区别？建立弱形式的关键步骤是什么？

开孔板分析

Step
00
概要
? 二维线性静力分析
开孔板
? 模型
- ? 对称模型 - 单位 : N, mm - 各向同性弹性材料 - 平面应力单元
? 荷载和边界条件
- 线压力 - 约束(对称)
? 模型
? 输出结果
- 范.梅塞斯应力 - 主应力向量 - 结果标注
开孔板
Step
01
1 2 3 4 5 6
分析 > 分析控制 操作步骤
– 控制 表单
分析类型 : [选择2D和XY平面] 点击 [
] [单位系统]
1
长度 : [mm] 点击 [确认] 键 点击[确认] 键 在工作窗口点击鼠标右键并选择[隐藏基 2
准面，栅格，三角标]
5
6 3
4
________________________________________________________________________________________________________________________________
开孔板 开孔板
|
1

Step
02
1 2 3 4 5
几何 > 曲线 > 在工作平面上创建 > 操作步骤
2D 矩形(线框)… , 圆…
3
位置:“(0) , <10, 30>”回车 点击[取消] 键 点击[显示全部] 中心:“(0)”, 半径:“5”回车 点击[取消] 键
?
2
?
1
1
2
4
?
1
( ): “坐标 x, y”, 输入一个点回车 < >: “相对 dx, dy”，输入一个点回车 (0) 等于 (0, 0) 5
?
2
[Esc] 是命令[取消]的快捷键.
开孔板
Step
03
1 2 3 4 5 6
几何 > 曲线 > 操作步骤
交叉线…
点击[
] 选择屏显对象
?
2
?
1
点击[适用] 键 点击[取消] 键
2
3
选择有标记[ O ]的3条线 按键盘上的 [删除] 键 点击[确认] 键 6
4
?
1
“Ctrl+A” 是“选择屏显示对象”快 捷键.
?
2
[Enter] 是[适用]命令的快捷键.
________________________________________________________________________________________________________________________________
2
| 开孔板
开孔板

带孔平板的应力集中分析

有限元方法 Finite Element Method ——基于ANSYS的有限元建模与分析 姓名吴威 学号20100142 班级10级土木茅以升班2班 西南交通大学 2014年4月

综合练习——带孔平板的应力分布及应力集中系数的计算一、问题重述 计算带孔平板的应力分布及应力集中系数。 二、模型的建立与计算 在ANSYS中建立模型，材料的设置属性如下 分析类型为结构（structural），材料为线弹性（Linear Elastic），各向同性（Isotropic）。弹性模量、泊松比的设定均按照题目要求设定，以N、cm为标准单位，实常数设置中设板厚为1。

采用solid 4 node 42板单元，Element Behavior设置为Plane strs w/thk。 建立模型时先建立完整模型，分别用单元尺度为5cm左右的粗网格和单元尺度为2cm左右的细网格计算。 然后取四分之一模型计算比较精度，为了使粗细网格单元数与完整模型接近，四分之一模型分别用单元尺度为2.5cm左右的粗网格和单元尺度为1cm左右的细网格计算。 (1) 完整模型的计算 ①粗网格

单元网格的划分及约束荷载的施加如图（单元尺度为5cm） 约束施加时在模型左侧边界所有节点上只施加x方向的约束，即令U X=0，在左下角节点上施加x、y两个方向的约束，即U X=0、U Y=0。荷载施加在右侧边界上，大小为100。 对模型进行分析求解得到： 节点应力云图（最大值222.112）

单元应力云图（最大值256.408） 可看出在孔周围有应力集中现象，其余地方应力分布较为均匀，孔上部出现最大应力。 ②细网格 单元网格的划分及约束荷载的施加如图（单元尺度为2cm）

有限元分析的一般过程

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

孔板结构原理特点

孔板结构原理特点 中文名称：孔板 英文名称：orifice plate 定义：设在管道或泄水隧洞中，用以量测流量或进行消能的开孔隔板。 孔板orifice plate 安装在封闭管道中，按节流装置的原理，测量液体、气体和蒸汽流量的检出元件。 标准孔板是一块具有圆形开孔的金属薄板,圆孔壁与孔板前端面成直角,安装时孔板轴心与管道轴线同心。 孔板已有国际标准(ISO 5167)。测出孔板两端压差，按此标准即可算出具有一定精确度的流量值。孔板取压方式在国际标准中规定为径距取压、法兰取压和角接取压（取压孔紧靠孔板）3种。当测量含有少量固体的液体或含有少量液体的气体时，为便于少量固体或液体通过，孔板的开孔可制成扇形的,或制成与管道的轴线是偏心的. 孔板是测量流量的差压发生装置，配合各种差压计或差压变送器可测量管道中各种流体的流量。节流装置包括环室孔板，喷嘴等。节流装置与差压变送器配套使用，可测量液体、蒸汽、气体的流量，它广泛应用于石油、化工、冶金、电力、轻工等部门。 充满管道的流体，当它们流经管道内的节流装置时，流束将在节流装置的节流件处形成局部收缩，从而使流速增加，静压力低，于是在节流件前后便产生了压力降，即压差，介质流动的流量越大，在节流件前后产生的压差就越大，所以可以通过测量压差来衡量流体流量的大小。这种测量方法是以能量守衡定律和流动连续性定律为基准的。 智能节流装置（流量计）是集流量、温度、压力检测功能于一体，并能进行温度、压力自动补偿的新一代流量计，该仪表采用先进的微机技术与微功耗新技术，功能强，结构紧凑，操作简单，使用方便。 特点 ▲节流装置结构易于复制，简单、牢固，性能稳定可靠，使用期限长，价格低廉。 ▲孔板计算采用国际标准与加工 ▲应用范围广，全部单相流皆可测量，部分混相流亦可应用。 ▲标准型节流装置无须实流校准，即可投用。 ▲一体型孔板安装更简单，无须引压管，可直接接差压变送器和压力变送器。 智能型特点 ▲采用进口单晶硅智能差压传感器 ▲高精度，完善的自诊断功能 ▲智能孔板流量计其量程可自编程调整。 ▲可同时显示累计流量、瞬时流量、压力、温度。 ▲具有在线、动态全补偿功能外，还具有自诊断、自行设定量程。 ▲配有多种通讯接口

有限元分析基础教程(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