有限元法问答题整理
有限元理论方法
关于有限元分析法及其应用举例 摘要：本文主要介绍有限元分析法，作为现代设计理论与方法的一种，已经在 众多领域普遍使用。介绍了它的起源和国内外发展现状。阐述了有限元法的基 本思想和设计方法。并从实际出发，例举了有限元法的一个简单应用———啤 酒瓶的应力分析和优化,表明了利用有限元分析法的众多优点。随着计算机的 发展，基于有限元分析方法的软件开发越来越多。本文也在其软件开发方面进 行阐述，并简单介绍了一下主流软件的发展情况和使用范围。并就这一领域的 未来发展趋势进行阐述。 关键词：有限元分析法软件啤酒瓶 Abstract：This thesis mainly introduces the finite element analysis, as a modern design theory and methods used widely in in most respects. And this paper introduces its origins and development in world. It also expounds the basic thinking and approach of FEM..Proceed from the actual situation,this text holds the a simple application of finite-element method———the analysis and optimized of an beer bottle and indicate the the numerous benefits of finite element analysis .As computers mature and based on the finite element analysis of the software development is growing. This article introduces its application in the software development aspects as well, and briefly states the development and scope of the mainstream software. And it’s also prospect future development tendency in this area . Key: Finite Element Analysis Software Beer bottle 0 绪论 有限元法(Finite Element Method，FEM)，是计算力学中的一种重要的方法，它是20世纪50年代末60年代初兴起的应用数学、现代力学及计算机科学相互渗透、综合利用的边缘科学。有限元法最初应用在工程科学技术中，用于模拟并且解决工程力学、热学、电磁学等物理问题。对于过去用解析方法无法求解的问题和边界条件及结构形状都不规则的复杂问题，有限元法则是一种有效的分析方法。有限元法的基本思想是先将研究对象的连续求解区域离散为一组有限个且按一定方式相互联结在一起的单元组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模拟成不同几何形状的求解小区域；
第18章 接触问题有限元分析技术
第18章接触问题的有限元分析技术 第1节基本知识 接触问题是一种高度非线性行为，需要较大的计算资源，为了进行准确而有效的计算，理解问题的特性和建立合理的模型是很重要的。 接触问题存在两个较大的难点：其一，在求解问题之前，不知道接触区域，表面之间是接触或分开是未知的、突然变化的，这些随载荷、材料、边界条件和其它因素而定；其二，大多数的接触问题需要计算摩擦，有几种摩擦和模型可供挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。 一、接触问题分类 接触问题分为两种基本类型：刚体─柔体的接触和半柔体─柔体的接触。在刚体─柔体的接触问题中，接触面的一个或多个被当作刚体，（与它接触的变形体相比，有大得多的刚度），一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体─柔体的接触，许多金属成形问题归为此类接触；另一类，柔体─柔体的接触，是一种更普遍的类型，在这种情况下，两个接触体都是变形体（有近似的刚度）。 ANSYS支持三种接触方式：点─点、点─面和平面─面。每种接触方式使用的接触单元适用于某类问题。 二、接触单元 为了给接触问题建模，首先必须认识到模型中的哪些部分可能会相互接触，如果相互作用的其中之一是一点，模型的对立应组元是一个节点。如果相互作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元。有限元模型通过指定的接触单元来识别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元。下面分类详述ANSYS使用的接触单元和使用它们的过程。 1．点─点接触单元 点─点接触单元主要用于模拟点─点的接触行为，为了使用点─点的接触单元，需要预先知道接触位置，这类接触问题只能适用于接触面之间有较小相对滑动的情况（即使在几何非线性情况下）。 如果两个面上的节点一一对应，相对滑动又以忽略不计，两个面挠度（转动）保持小量，那么可以用点─点的接触单元来求解面─面的接触问题，过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。 2．点─面接触单元 点─面接触单元主要用于给点─面的接触行为建模，例如两根梁的相互接触。 如果通过一组节点来定义接触面，生成多个单元，那么可以通过点─面的接触单元来模拟面─面的接触问题，面即可以是刚性体也可以是柔性体，这类接触问题的一个典型例子是
★★★装配体有限元分析
基于ANSYS WORKBENCH的装配体有限元分析 模拟装配体的本质就是设置零件与零件之间的接触问题。 装配体的仿真所面临的问题包括： （1）模型的简化。这一步包含的问题最多。实际的装配体少的有十几个零件，多的有上百个零件。这些零件有的很大，如车门板；有的体积很小，如圆柱销；有的很细长，如密封条；有的很薄且形状极不规则，如车身；有的上面钻满了孔，如连接板；有的上面有很多小突起，如玩具的外壳。在对一个装配体进行分析时，所有的零件都应该包含进来吗？或者我们只分析某几个零件？对于每个零件，我们可以简化吗？如果可以简化，该如何简化？可以删除一些小倒角吗？如果删除了，是否会出现应力集中？是否可以删除小孔，如果删除，是否会刚好使得应力最大的地方被忽略？我们可以用中面来表达板件吗？如果可以，那么，各个中面之间如何连接？在一个杆件板件混合的装配体中，我们可以对杆件进行抽象吗？或者只是用实体模型？如果我们做了简化，那么这种简化对于结果造成了多大的影响，我们可以得到一个大致的误差范围吗？所有这些问题，都需要我们仔细考虑。 （2）零件之间的联接。装配体的一个主要特征，就是零件多，而在零件之间发生了关系。我们知道，如果零件之间不能发生相对运动，则直接可以使用绑定的方式来设置接触。如果零件之间可以发生相对运动，则至少可以有两种选择，或者我们用运动副来建模，或者，使用接触来建模。如果使用了运动副，那么这种建模方式对于零件的强度分析会造成多大的影响？在运动副的附近，我们所计算的应力其精确度大概有多少？什么时候需要使用接触呢？又应该使用哪一种接触形式呢？ （3）材料属性的考虑。在一个复杂的装配体中所有的零件，其材料属性多种多样。我们在初次分析的时候，可以只考虑其线弹性属性。但是对于高温，重载，高速情况下，材料的属性不再局限于线弹性属性。此时我们恐怕需要了解其中的每一种材料，它是超弹性的吗？是哪一种超弹性的？它发生了塑性变形吗？该使用哪一种塑性模型？它是粘性的吗？它是脆性的吗？它的属性随着温度而改变吗？它发生了蠕变吗？是否存在应力钢化问题？如此众多的零件，对于每一个零件，我们都需要考察其各种各样的力学属性，这真是一个丰富多彩的问题。（4）有限元网格的划分。我们知道，通过WORKBENCH，我们只需要按一个按钮，就可以得到一个粗糙的网格模型。但是如果从HYPERMESH的角度来看，ANSYS自动划分的网格，很多都是不合理的，质量较差而不能使用。那么对于装配体中的每个零件，我们该如何划分网格？对于每一个零件，我们是否要对之进行切割形成规则的几何体后，然后尽量使用六面体网格？如果
有限元法及其在工程中的应用
机械与汽车学院 曹国强 主要内容： 1、有限元法的基本思想。 2、结构力学模型的简化和结构离散化。 3、有限元法的实施过程。 一、有限元法的基本思想 有限元法是随着计算机的发展而发展起来的一种有效的数值方法。其基本思想是：将连续的结构分割成数目有限的小单元体（称为单元），这些小单元体彼此之间只在数目有限的指定点（称为节点）上相互连接。用这些小单元体组成的集合体来代替原来的连续结构。再把每个小单元体上实际作用的外载荷按弹性力学中的虚功等效原理分配到单元的节点上，构成等效节点力，并按结构实际约束情况决定受约束节点的约束。这一过程称为结构的离散化。其次，对每个小单元体选择一个简单的函数来近似地表示其位移分量的分布规律，并按弹性力学中的变分原理建立起单元节点力和节点位移之间的关系（单元刚度方程），最后，把全部单元的节点力和节点位移之间的关系组集起来，就得到了一组以结构节点位移为未知量的代数方程组（总体刚度方程），同时考虑结构的约束情况，消去那些结构节点位移为零的方程，再由最后的代数方程组就可求得结构上有限个离散节点的各位移分量。求得了结构上各节点的位移分量之后，即可按单元的几何方程和物理方程求得各单元的应变和应力分量。 有限元法的实质就是把具有无限个自由度的连续体，理想化为有限个自由度的单元的集合体，使问题简化为适合于数值解法的结构型问题。 经典解法（解析法）与有限元法的区别 解析法 { } 建立一个描述连续体性质的偏微分方程组 有限元解法 连续体 数目增加到∞ 大小趋于0 微元 有限元 离散化 （单元分析）集合 总体分析 求得近似解
二、结构力学模型的简化和结构离散化 （一）结构力学模型的简化 用有限元法研究实际工程结构问题时，首先要从工程实际问题中抽象出力学模型，即要对实际问题的边界条件、约束条件和外载荷进行简化，这种简化应尽可能地反映实际情况，不至于使简化后的解答与实际差别过大，但也不要带来计算上的过分复杂，在力学模型的简化过程中，必须判断实际结构的问题类型，是二维问题还是三维问题。如果是平面问题，是平面应力问题，还是平面应变问题。同时还要搞清楚结构是否对称，外载荷大小和作用位置，结构的几何尺寸和力学参数（弹性模量E、波松比μ等）。 （二）结构的离散化 将已经简化好的结构力学模型划分成只在一些节点连续的有限个单元，把每个单元看成是一个连续的小单元体，各单元之间只在一些点上互相联结，这些点称作节点，每个单元体称为一个单元。用只在节点处连接的单元的集合体代替原来的连续结构，把外载荷按虚功等效原理移置到有关受载的节点上，构成节点载荷，把连续结构进行这样分割的过程称为结构的离散化。现举例说明。 设一平面薄板，中间有一个园孔，其左端固定，右端受面力载荷q，试对其进行有限元分割和力学模型简化。
Hypermesh2017.2有限元分析的前处理1D单元连接
ALTAIR HYPERWORKS2017.2 有限元分析前处理 1D 单元和连接 Trainer’s Name Month XX, 2017
HMD Intro, 2017.2第5章: 1D 单元和焊点 5) 1D 单元和焊点 ?1D Meshing （1D单元） ?HyperBeam （梁截面） ?Connectors （焊点）
HMD Intro, 2017.2 1D 单元 ?1D 单元
HMD Intro, 2017.2示例 跟着示范做 (…\Model-Files\CH5-1D-MESHING\05a-1D-MESHING.hm)
? 2017 Altair Engineering, Inc. Proprietary and Confidential. All rights reserved. HMD Intro, 2017.2 1D单元介绍 ?1D单元是节点之间简单连接，允许精确模拟连接关系(例如螺栓)和类似的杆 状或杆状对象，这些对象在FEA模型中可以建模为简单的线 ?可以从以下面板创建1D单元: ?目前支持的1D单元包括: bar2s, bar3s, rigid links, rbe3s, plots, rigids, rods, springs, welds, gaps and joints. ?显示单元可以在以下面板中创建: Edit Element, Line Mesh, Elem Offset, Edges, or Features panel.
?RIGID 刚性连接用于传递从主节点到从节点的运动. ?Rigids面板允许创建rigid 和rigid link 单元.
对有限元方法的认识
我对有限元方法的认识 1有限元法概念 有限元方法（The Finite Element Method, FEM）是计算机问世以后迅速发展起来的一种分析方法。每一种自然现象的背后都有相应的物理规律，对物理规律的描述可以借助相关的定理或定律表现为各种形式的方程（代数、微分、或积分）。这些方程通常称为控制方程（Governing equation）。 针对实际的工程问题推导这些方程并不十分困难，然而，要获得问题的解析的数学解却很困难。人们多采用数值方法给出近似的满足工程精度要求的解答。 有限元方法就是一种应用十分广泛的数值分析方法。 有限元方法是处理连续介质问题的一种普遍方法，离散化是有限元方法的基础。 这种思想自古有之：古代人们在计算圆的周长或面积时就采用了离散化的逼近方法：即采用内接多边形和外切多边形从两个不同的方向近似描述圆的周长或面积，当多边形的边数逐步增加时近似值将从这两个方向逼近真解。 近年来随着计算机技术的普及和计算速度的不断提高，有限元分析在工程设计和分析中得到了越来越广泛的重视，已经成为解决复杂的工程分析计算问题的有效途径，现在从汽车到航天飞机几乎所有的设计制造都已离不开有限元分析计算，其在机械制造、材料加工、航空航天、汽车、土木建筑、电子电器、国防军工、船舶、铁道、石化、能源、科学研究等各个领域的广泛使用已使设计水平发生了质的飞跃。 国际上早在 60 年代初就开始投入大量的人力和物力开发有限元分析程序。“有限单元”是由Clough R W于1960年首次提出的。但真正的有限元分析软件是诞生于 70 年代初期，随着计算机运算速度的提高，内、外存容量的扩大和图形设备的发展，以及软件技术的进步，发展成为有限元分析与设计软件，但初期其前后处理的能力还是比较弱的，特别是后处理能力更弱。
ABAQUS有限元接触分析的基本概念
CAE(计算机辅助工程)是一门复杂的工程科学，涉及仿真技术、软件、产品设计和力学等众多领域。世界上几大CAE公司各自以其独到的技术占领着相应的市场。ABAQUS有限元分析软件拥有世界上最大的非线性力学用户群，是国际上公认的最先进的大型通用非线性有限元分析软件之一。它广泛应用于机械制造、石油化工、航空航天、汽车交通、土木工程、国防军工、水利水电、生物医学、电子工程、能源、地矿、造船以及日用家电等工业和科学研究领域。ABAQUS在技术、品质和可靠性等方面具有卓越的声誉，可以对工程中各种复杂的线性和非线性问题进行分析计算。 《ABAQUS有限元分析常见问题解答》以问答的形式，详细介绍了使用ABAQUS 建模分析过程中的各种常见问题，并以实例的形式教给读者如何分析问题、查找错误原因和尝试解决办法，帮助读者提高解决问题的能力。 《ABAQUS有限元分析常见问题解答》一书由机械工业出版社出版。 16.1.1点对面离散与面对面离散 【常见问题16-1】 在ABAQUS/Standard分析中定义接触时，可以选择点对面离散方法(node-to-surface-dis-cre-tization)和面对面离散方法(surface-to-surfacediscretization)，二者有何差别? 『解答』 在点对面离散方法中，从面(slavesurface)上的每个节点与该节点在主面(mastersurface)上的投影点建立接触关系，每个接触条件都包含一个从面节点和它的投影点附近的一组主面节点。 使用点对面离散方法时，从面节点不会穿透(penetrate)主面，但是主面节点可以穿透从面。 面对面离散方法会为整个从面(而不是单个节点)建立接触条件，在接触分析过程中同时考虑主面和从面的形状变化。可能在某些节点上出现穿透现象，但是穿透的程度不会很严重。 在如图16-l和图16-2所示的实例中，比较了两种情况。
有限元法分析过程
有限元法分析过程 有限元法分析过程大体可分为：前处理、分析、后处理三大步骤。 对实际的连续体经过离散化后就建立了有限元分析模型，这一过程是有限元的前处理过程。在这一阶段，要构造计算对象的几何模型，要划分有限元网格，要生成有限元分析的输入数据，这一步是有限元分析的关键。 有限元分析过程主要包括：单元分析、整体分析、载荷移置、引入约束、求解约束方程等过程。这一过程是有限元分析的核心部分，有限元理论主要体现在这一过程中。 有限元法包括三类：有限元位移法、有限元力法、有限元混合法。 在有限元位移法中，选节点位移作为基本未知量； 在有限元力法中，选节点力作为未知量； 在有限元混合法中，选一部分基本未知量为节点位移，另一部分基本未知量为节点力。 有限元位移法计算过程的系统性、规律性强，特别适宜于编程求解。一般除板壳问题的有限元应用一定量的混合法外，其余全部采用有限元位移法。因此，一般不做特别声明，有限元法指的是有限元位移法。 有限元分析的后处理主要包括对计算结果的加工处理、编辑组织和图形表示三个方面。它可以把有限元分析得到的数据，进一步转换为设计人员直接需要的信息，如应力分布状态、结构变形状态等，并且绘成直观的图形，从而帮助设计人员迅速的评价和校核设计方案。 附：FELAC 2.0软件简介 FELAC 2.0采用自定义的有限元语言作为脚本代码语言，它可以使用户以一种类似于数学公式书写和推导的方式，非常自然和简单的表达待解问题的微分方程表达式和算法表达式，并由生成器解释产生完整的并行有限元计算C程序。 FELAC 2.0的目标是通过输入微分方程表达式和算法之后，就可以得到所有有限元计算的程序代码，包含串行程序和并行程序。该系统采用一种语言(有限元语言)和四种技术(对象技术、组件技术、公式库技术生成器技术)开发而成。并且基于FELAC 1.0的用户界面，新版本扩充了工作目录中右键编译功能、命令终端输入功能，并且丰
ABAQUS有限元接触分析的基本概念
ABAQUS有限元接触分析的基本概念2009-11-24 00:06:28 作者：jiangnanxue 来源：智造网—助力中国制造业创新—https://www.360docs.net/doc/a95647508.html, CAE(计算机辅助工程)是一门复杂的工程科学，涉及仿真技术、软件、产品设计和力学等众多领域。世界上几大CAE公司各自以其独到的技术占领着相应的市场。ABAQUS有限元分析软件拥有世界上最大的非线性力学用户群，是国际上公认的最先进的大型通用非线性有限元分析软件之一。它广泛应用于机械制造、石油化工、航空航天、汽车交通、土木工程、国防军工、水利水电、生物医学、电子工程、能源、地矿、造船以及日用家电等工业和科学研究领域。ABAQUS在技术、品质和可靠性等方面具有卓越的声誉，可以对工程中各种复杂的线性和非线性问题进行分析计算。 《ABAQUS有限元分析常见问题解答》以问答的形式，详细介绍了使用ABAQUS建模分析过程中的各种常见问题，并以实例的形式教给读者如何分析问题、查找错误原因和尝试解决办法，帮助读者提高解决问题的能力。 《ABAQUS有限元分析常见问题解答》一书由机械工业出版社出版。 16.1.1 点对面离散与面对面离散 【常见问题16-1】 在ABAQUS/Standard分析中定义接触时，可以选择点对面离散方法(node-to-surface-dis - cre-tization)和面对面离散方法(surface-to-surface discretization)，二者有何差别? 『解答』 在点对面离散方法中，从面(slave surface)上的每个节点与该节点在主面(master surface)上的投影点建立接触关系，每个接触条件都包含一个从面节点和它的投影点附近的一组主面节点。 使用点对面离散方法时，从面节点不会穿透(penetrate)主面，但是主面节点可以穿透从面。 面对面离散方法会为整个从面(而不是单个节点)建立接触条件，在接触分析过程中同时考虑主面和从面的形状变化。可能在某些节点上出现穿透现象，但是穿透的程度不会很严重。 在如图16-l和图16-2所示的实例中，比较了两种情况。
第一章 计算机工程分析前处理和后处理
第一章计算机工程分析的前处理和后处理 本章教学任务要点：通过本章的系统学习，要求学生能够掌握有限单元法计算的前、后处理基本知识，基本原理和基本方法。 §1. 工程对象及选择 例1.工程对象：进行西安城区地裂缝模拟： 研究目标：未来100年，地裂缝活动特点、以及对市政规划的影响。 例 研究目标：未来100年，地裂缝活动特点、以及对市政规划的影响。 选择对象： A位置：隧道埋深 B位置：隧道埋深10m，地面为高层建筑。
例3.工程对象：煤层开采模拟 研究目标：采动后岩层运动规律、离层破坏高度 选择对象： §2. 有限元计算模型的原则：计算模型要能全面反映工程对象的主要特点，又必须具备能适应计算模拟的功能。 例一、工程对象：进行西安城区地裂缝模拟 例二、城市地铁隧道 例三、煤层开采模拟 有限元计算模型的建立主要指，1、将待分析的连续体，如结构物、固体等对象，用假想的点、线、面将连续体分割成有限多个、有限大小的子区域，这些子区域只有在特定点相互连接，从而使得连续体离散化为结构体。其中，这些子区域称为单元，单元与单元之间的连接点称为结点。 单元类型有： 直线单元曲线元 A一维单元
三角形 矩形四边形曲边形 B 二维单元 六面体单元柱体单元 C 三维单元 结点类型：结点与结点之间的连接可以有铰接、固接。如果结构体的一个结点位移或一个方向被限制，则结点上可以安装铰支座、或杆支座等。 上机 §3. 3.1 单元， 单元号，
例图巨型划分成小巨型格，也就是将巨型（1000X1000）划分成小巨型（5格X 5格），其中如：单元n、节点为A、B、C、D。该单元其坐标为： A（400，400），B（600，400），C（600，600），D（400，600）。 用传统解析方法定量地处理岩石力学问题，由于无法考虑复杂的岩石性态及某些明显的地质构造的影响，而存在着几乎是不可克服的困难，正是在这方面有限单元法具有突出的效能。 1）有限单元法使连续体离散化；可以反映复杂实际对象。 2）单元应力状态可以指定，因而可以反映岩体实际状态。 3）处理格式一致。 3.2 网格单元形式 实践证明，采用常应变三角形单元或具有线性应变的任意四边形单元均能达到工程所要求的精度。 在计算精度方面，高阶应变特性的单元对提高计算精度是有限度的，而三角形单元精度略次于四边形单元，且处理岩石某些性态及进行非线性分析时程序处理上不如四边形方便。 3.3 离散化 正确确定边界及其位置条件，网格的细度，对有限元分析结果有很大影响。在保证要求精度的情况下，确定一个最低限度的网格要求是很必要的。例如根据大量计算经验，对地下硐室围岩分析，网格的适度细化，按节点数来说，应不少300节点。对具有一个对称轴时，节点总数不少于150~200个；对具有两个对称轴时，节点数不少于600个。 对网格疏密度，在应力变化处，可密；否则疏。例如对地下硐室围岩1~1.5倍直径区域加密。 对研究范围的大小及边界位置：如圆孔应力集中问题中，计算精度与精确解的误差在10%以内，必要的边界条件，一般应在距坑硐中心不少于坑硐直径3.5 ~ 4倍。又如，对弹模很低（1×105以下）的岩石，则所取范围尚应适当增大。 应指明：采用“翻转应力法”模拟开挖，边界均以固定点来考虑。 3.4 网格剖分的处理方法 地下硐室围岩应力分布问题，单元划分情况可分两种类型。即取全部或部分对称部分。
空间问题的有限元方法总结计划.docx
第三章 空间问题的有限元方法 引言 许多工程实际问题，属于空间问题，由于结构形状或受力的复杂性，用 经典弹性理论去求解它们的解析解是不可能的。 而有限元法处理此类问题， 原则 上不存在什么困难，本章将介绍一般空间问题的四面体单元。 一般空间问题的有限元列式 3.2.1 单元位移模式及插值函数 空间问题中，每个单元有四个结点，编码为 i,j,m,p 。每个结点有 3 个位移分量。每个结点 的位移可用位移矢量 i 表示，即 u i i v i (i , j ,m, p) w i 单元结点的位移向量可表示为 i e j u i v i w i u j v j w j u m v m w m u p v p w p T m p e 为单元结点位移列阵。 假设单元内的位移模式选取一次多项式 u 1 2 x 3 y 4 z v 5 6 x 7 y 8 z （3.2.1 ） w 9 10 x 11 y 12 z 由于四个结点也在单元内，满足位移模式，于是得 u i 12 x i 3 y i 4 z i u j 1 2 x j 3 y j 4 z j （ 3.2.2 ） u m 1 2 x m 3 y m 4 z m u p 1 2 x p 3 y p 4 z p 上式是关于 1 , 2, 3, 4 的线性方程组。 1, 2 , 3, 4 是待定常数，也称为广义坐
标。它可由（ 3.2.2 ）式求出。上式的系数行列式是 1x i y i z i 1x j y j z j 2V（3.2.3 ） D x m y m z m 1 1x p y p z p 上式中当 i,j,m,p 的编号顺序满足右手法则， V值为正，其大小为四面体体积，因此为了方便单元的编号一般满足右手法则。求得1 , 2 , 3 , 4后,回代入位移模式得 u N i u i N j u j N m u m N p u p(3.2.4) 式中 N i 1 (a i b i x c i y d i z)(i , j, m, p) (3.2.5) 6V x j y j z j a i x m y m z m x p y p z p 1y j z j b i1y m z m 1y p z p 1x j z j c i 1x m z m(i , j , m, p) （3.2.6） 1x p z p 1x j y j d i1x m y m 1x p y p 上式下标 (i ,j , m, p) 轮换，可得 a j , b j ,c j , d j, a m ,b m ,c m , d m及 a p , b p , c p ,d p。同理 , 也可得到其它两式 , 于是得 u N i u i N j u j N m u m N p u p v N i v i N j v j N m v m N p v p（ 3.2.7）
第四章空间问题的有限元
第四章 空间问题的有限元 在工程问题中，有些结构形状非常复杂，必须按照空间问题来求解。由于4节点四面体单元可以很好的模拟几何体的边界形状而被广泛使用。因此本章将介绍此种单元及8节点六面体单元。 §4.1 空间问题的离散化 在工程实际中，有些结构由于形体复杂，并且三个方向的尺寸同量级，必须按空间问题求解。空间问题有限元法的原理、思路和解题方法完全类同于平面问题的有限元法，所不同的是它具有三维特点。它所采用的离散化模型仍然是由若干单元在节点处连接而成的，而且节点仍为铰接，但是这些单元具有块体形状。它的基本未知量是节点位移，有3个分量：,,u v w 。它的分析方法仍然是先进行单元分析，再进行整体分析，最后求解整体平衡方程。但必须指出，由平面问题转换为空间问题给有限元分析带来了两个主要困难： 1、空间结构离散不像平面问题直观，当人工离散时很容易产生错误。 2、未知量的数量剧增，对于比较复杂的空间问题，计算机存储容量和计算机费用都会产生问题。 为解决上述两个问题，前者可通过寻找规律，建立网格自动生成前处理程序来克服，而后者则可采用高阶元以提高单元精度，达到减少未知量和节省机时的目的。 §4.2常应变四面体单元 §4.2.1位移函数 图4-1所示为四面体单元，以四个角点i ，j ，m ，l 为结点，每个结点有三个自 由度，因此由广义坐标给出的线性位移函数为 000000u ??β?β??? ??==?? ???? (4.2.1) 其中[]1x y z ?= 图4-1 四面体单元 []1212T ββββ=L 把四个节点坐标代入(4.2.1)式时，可得
{}000 000A q A A A ββ?? ??==?? ??? ? %%% (4.2.2) 其中{}T i i i j j j m m m l l l q u v w u v w u v w u v w ??=?? 1111i i i j j j m m m l l l x y z x y z A x y z x y z ??????=??? ??? % 由(4.2.2)式求出 {}1A q β-=% (4.2.3) 将(4.2.3)式代入(4.2.1)式后，则有 {}{}1i j m l u B A q N N N N q -??=Φ=Φ=I I I I ??&% (4.2.4) 其中100010001????I =?? ???? ()1 6i i i i i N a b x c y d z V = +++ ()1 6j j j j j N a b x c y d z V =- +++ ()1 6m m m m m N a b x c y d z V = +++ ()1 6l l l l l N a b x c y d z V =- +++ 称为形函数，它们的系数为 i j j i m m m l l l x y z a x y z x y z = 1 11j j i m m l l y z b y z y z = 111 j j i m m l l x z c x z x z = 111j j i m m l l x y d x y x y =
基于有限元软件ABAQUS的过盈接触分析
基于有限元软件ABAQUS的过盈接触分析 如下图所示，将轴缓缓压入轴毂中，轴和毂之间在径向有8mm的过盈量，轴毂固定，两者的材料均为钢，弹性模量为2.06E11Pa，泊松比为0.3，摩擦系数为0.2。分析装配过程中轴和轴毂的应力应变情况。 问题分析 （1）本题主要分析装配过程中结构的静态响应，所以分析步选择通用静态分析步。 （2）本题由于为过盈配合，属于大变形，故应考虑几何非线性的影响。 （3）模型具有轴对称性，所以可以采取轴对称模型来进行分析，这样可以节省计算时间。 （4）为了方便收敛，分析步可以分成两步，第一步建立两者间的接触关系，第二步完成过盈装配。 （5）接触面之间有很大的相对滑动，所以模型要使用有限滑移（Finite sliding）。 ABAQUS/CAE分析过程如下： （1）进入Part模块，创建Name为Axis的部件
在草图环境中输入（0,0），（0.1,0），（0.1,0.12），（0.13,0.12），（0.13,0.28），（0,0.28），（0,0）同时为轴部件端部切割出一78度角的倒角 同样再创造一Name为Hub的部件，设置与Axis一样，在草图环境中输入利用Rectangle工具创建一矩形，两角点为（0.09992,0）和（0.19992，-0.12）
（2） 进入property 模块，定义材料属性
（3）进入Assembly模块，创建两者间的装配关系
（4）进入step模块 定义名为Make-Contact和Press-Axis-Down的两个分析步，，将Nlgeom设置 为on，详细信息如下：
ABAQUS有限元接触分析的基本概念
ABAQUS有限元接触分析的基本概念 来源：机械工业出版社《ABAQUS有限元分析常见问题解答》 CAE(计算机辅助工程)是一门复杂的工程科学，涉及仿真技术、软件、产品设计和力学等众多领域。世界上几大CAE公司各自以其独到的技术占领着相应的市场。ABAQUS有限元分析软件拥有世界上最大的非线性力学用户群，是国际上公认的最先进的大型通用非线性有限元分析软件之一。它广泛应用于机械制造、石油化工、航空航天、汽车交通、土木工程、国防军工、水利水电、生物医学、电子工程、能源、地矿、造船以及日用家电等工业和科学研究领域。ABAQUS在技术、品质和可靠性等方面具有卓越的声誉，可以对工程中各种复杂的线性和非线性问题进行分析计算。 《ABAQUS有限元分析常见问题解答》以问答的形式，详细介绍了使用ABAQUS建模分析过程中的各种常见问题，并以实例的形式教给读者如何分析问题、查找错误原因和尝试解决办法，帮助读者提高解决问题的能力。 《ABAQUS有限元分析常见问题解答》一书由机械工业出版社出版。 16.1.1 点对面离散与面对面离散 【常见问题16-1】 在ABAQUS/Standard分析中定义接触时，可以选择点对面离散方法(node-to-surface-dis- cre-tization)和面对面离散方法(surface-to-surface discretization)，二者有何差别? 『解答』 在点对面离散方法中，从面(slave surface)上的每个节点与该节点在主面(master surface)上的投影点建立接触关系，每个接触条件都包含一个从面节点和它的投影点附近的一组主面节点。 使用点对面离散方法时，从面节点不会穿透(penetrate)主面，但是主面节点可以穿透从面。 面对面离散方法会为整个从面(而不是单个节点)建立接触条件，在接触分析过程中同时考虑主面和从面的形状变化。可能在某些节点上出现穿透现象，但是穿透的程度不会很严重。 在如图16-l和图16-2所示的实例中，比较了两种情况。
有限元分析中的一些问题
有限元分析的一些基本考虑-----单元形状对于计算精度的影响 笔者发现，在分析复杂问题时，我们所可能出现的错误，竟然是一些很根本的错误，这些根本错误是由于对有限元的基本理论理解不清晰而造成的。 鉴于这个原因，笔者决定对一些基本问题（例如单元形状问题，单元大小问题，应力集中问题等）展开调查，从而形成了一系列文章，本篇文章是这些系列文章中的第一篇。 本篇文章先考虑有限元分析中的第一个基本问题：单元形状问题。 我们知道，单元形状对于有限元分析的结果精度有着重要影响，而对单元形状的衡量又有着诸多指标，为便于探讨，这里首先只讨论第一个最基本的指标：长宽比（四边形单元的最长尺度与最短尺度之比），而且仅考虑平面单元的长宽比对于计算精度的影响。 为此，我们给出一个成熟的算例。该算例是一根悬臂梁，在其端面施加竖直向下的抛物线分布载荷，我们现在考察用不同尺度的单元划分该梁时，对于A点位移的影响。 这五种不同的划分方式，都使用矩形单元，只不过各单元的长宽比不同。 例如第一种（1）AR=1.1，就是长宽比接近1； 第二种（2）AR=1.5,就是长宽比是1.5.其它类推。 第五种（5）AR=24，此时单元的长度是宽度的24倍。 现在我们看看按照这五种单元划分方式对于A点位移的影响，顺便我们也算出了B点的位移，结果见下表。
我们现在仔细查看一下上表，并分析其含义。 我们先考虑第一行，它是第一种单元划分情况，此时每个单元的长宽比是1.1，由此我们计算出A点，B点的垂直位移，可以看到，A点的竖直位移是-1.093英寸，而B点的竖直位移是-0.346英寸。而这两点我们都是可以用弹性力学的方式得到精确解的，其精确解分别是-1.152以及-0.360.这样，我们可以得到此时A点位移误差的百分比是 [(-1.093)-(-1.152)]/1.152 = 5.2%. 对于其它情况，也采用类似的方式得到A点位移误差的百分比。 从上表可以看出来，随着长宽比的增加，位移误差越来越大，竟然大到56%。因此，如果我们是用长宽比为24的单元进行划分的话，那么我们的结果可以说是完全错误的。 下面按照上表绘制出一张图，该图从形象的角度表达了上表的含义。 由此可见，长宽比越接近于1，那么结算结果越精确，越远离1，则误差越大。
基于ANSYS与ABAQUS的赫兹接触问题有限元分析对比
基于ANSYS与ABAQUS的赫兹接触问题 有限元分析对比 郭波 [长春设备工艺研究所，长春130012] [ 摘要] 分别应用ANSYS软件与ABAQUS软件求解某精密部件的赫兹接触问题，并通过实验结果验证有限元分析结果的计算精度，结果显示ANSYS软件在求解精密部件的赫兹接触问题方面具 有较高的求解精度。 [ 关键词]ANSYS，ABAQUS，赫兹接触，有限元 Finite element analysis of Hertz contact problem based on ANSYS and ABAQUS GUO Bo [Changchun Equipment &Technology Research Institute , Changchun 130012] [ Abstract ] Solve Hertz contact problem both using ANSYS and ABAQUS, then verificate the computational accuracy of the Finite element analysis. The result shows, it has higher calculation precision in terms of Hertz contact problem of ANSYS in solving precision components. [ Keyword ] ANSYS，ABAQUS, Hertz contact, Finite element analysis 1前言 接触分析能够解决典型的状态非线性问题，在工程中应用广泛。由于传统的赫兹接触理论是在许多假设的前提下推导的近似解，而且在许多场合这些假设的前提是不成立的，因此运用赫兹理论来解决接触问题存在一定的局限性。近年来，随着计算机技术的普及，各种商用有限元分析软件逐步发展，有限元方法已成为应用广泛并且实用高效的求解接触问题的数值分析方法。ANSYS软件和ABAQUS软件是工程中应用最广泛的商用有限元分析软件。ANSYS软件是融合结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。ABAQUS 是一套功能强大的工程模拟的有限元软件，其解决问题的范围从相对
有限元分析-清华大学教程
8.1 进入工程分析模块 8.2施加约束 8.3 施加载荷 8.4 静态有限元计算过程和后处理 8.5动态分析的前处理和显示计算结果8.6有限元分析实例 习题
工程分析指的是有限元分析，包括静态分析（Static Analyses）和动态分析。动态分析又分为限制状态固有频率分析（Frequency Analyses）和自由状态固有频率分析（Free Frequency Analyses），前者在物体上施加一定约束，后者的物体没有任何约束，即完全自由。 8.1 进入工程分析模块 1. 进入工程分析模块前的准备工作 （1）在三维实体建模模块建立形体的三维模型，为三维形体添加材质，见4.7。 （2）将显示模式设置为Shading（着色）和Materials（材料），这样才能看到形体的应力和变形图，详见2.11.6。
2. 进入工程分析模块 选择菜单【Start】→【Analysis & Simulation】→【Generative Structural Analysis】弹出图8-1所示新的分析实例对话框。 在对话框中选择静态分析（Static Analyses）、限制状态固有频率分析（Frequency Analyses）还是自由状态固有频率分析（Free Frequency Analyses），单击OK按钮，将开始一个新的分析实例。 图8-1新的分析实例对话框
3.有限元分析的过程 有限元分析的一般流程为： （1）从三维实体建模模块进入有限元分析模块。（2）在形体上施加约束。 （3）在形体上施加载荷。 （4）计算（包括网格自动划分），解方程和生成应力应变结果。 （5）分析计算结果，单元网格、应力或变形显示。（6）对关心的区域细化网格、重新计算。 上述（1）~（3）过程是有限元分析预（前）处理，（4）是计算过程，（5）、（6）是有限元后处理。 有限元文件的类型为CATAnalysis。
有限元分析的一般过程
一、结构的离散化 将结构或弹性体人为地划分成由有限个单元，并通过有限个节点相互连接的离散系统。 这一步要解决以下几个方面的问题: 1、选择一个适当的参考系，既要考虑到工程设计习惯，又要照顾到建立模型的方便。 2、根据结构的特点，选择不同类型的单元。对复合结构可能同时用到多种类型的单元，此时还需要考虑不同类型单元的连接处理等问题。 3、根据计算分析的精度、周期及费用等方面的要求，合理确定单元的尺寸和阶次。 4、根据工程需要，确定分析类型和计算工况。要考虑参数区间及确定最危险工况等问题。 5、根据结构的实际支撑情况及受载状态，确定各工况的边界约束和有效计算载荷。 二、选择位移插值函数 1、位移插值函数的要求 在有限元法中通常选择多项式函数作为单元位移插值函数，并利用节点处的位移连续性条件，将位移插值函数整理成以下形函数矩阵与单元节点位移向量的乘积形式。 位移插值函数需要满足相容（协调）条件，采用多项式形式的位移插值函数，这一条件始终可以满足。 但近年来有人提出了一些新的位移插值函数，如：三角函数、样条函数及双曲函数等，此时需要检查是否满足相容条件。 2、位移插值函数的收敛性（完备性）要求： 1）位移插值函数必须包含常应变状态。 2）位移插值函数必须包含刚体位移。 3、复杂单元形函数的构造 对于高阶复杂单元，利用节点处的位移连续性条件求解形函数，实际上是不可行的。因此在实际应用中更多的情况下是利用形函数的性质来构造形函数。 形函数的性质： 1）相关节点处的值为 1，不相关节点处的值为 0。 2）形函数之和恒等于 1。 1、建立数学模型（特征消隐，理想化，清除）（（即从CAD 几何体→FEA 几何体），共 有下列三法：▲ 特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。▲ 理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理▲ 清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。） 2、建立有限元模型：（选择网格种类及定义分析类型；添加材料属性；施加约束；定义载 荷；网格划分） 3、求解有限元模型：再在此基础上计算应变和应力等其它物理量；在热分析中，FEA 首先 计算的是网格中每个节点的温度（标量），再在此基础上计算温度梯度和热流等其它物理量. 一般如果模型可划分网格，那么它就可以求解，但如果没有定义材料或载荷，则求解会终止。 4、结果分析：材料线性假设、小变形假设、静态载荷假设等等。