双重非对称转子支承系统实体有限元动力分析_马威猛_王建军

有限元动力学分析知识点

B 网格划分 1、网格划分三步骤: 定义单元属性、指定网格的控制参数、生成网格 2、单元属性（单元类型 (TYPE)、实常数 (REAL)、材料特性 (MAT)） 3、单元类型单元类型是一个重要选项，它决定如下单元特性：自由度(DOF)设置、单元形状、维数、假设的位移形函数。 1）线单元（梁单元、杆单元、弹簧单元） 2）壳用来模拟平面或曲面。 3）二维实体用于模拟实体截面 4）三维实体 ?用于几何属性，材料属性，荷载或分析要求考虑细节，而无法采用更简单的单元进行建模的结构。 ?也用于从三维CAD系统转化而来的几何模型，而这些几何模型转化成二维模型或壳体会花费大量的时间和精力 4、单元阶次与形函数 ?单元阶次是指单元形函数的多项式阶次。 ?什么是形函数? –形函数是指给出单元内结果形态的数值函数。因为FEA 的解答只是节点自由度值，需要通过形函数用节点自由度的值来描述单元内任一点的值。 –形函数根据给定的单元特性给出。

多体动力学软件和有限元软件的区别(优.选)

有限元软件与多体动力学软件数值分析技术与传统力学的结合在结构力学领域取得了辉煌的成就，出现了以ANSYS 、NASTRAN 等为代表的应用极为广泛的结构有限元分析软件。计算机技术在机构的静力学分析、运动学分析、动力学分析以及控制系统分析上的应用，则在二十世纪八十年代形成了计算多体系统动力学，并产生了以ADAMS 和DADS 为代表的动力学分析软件。两者共同构成计算机辅助工程（CAE ）技术的重要内容。商业通用软件的广泛应用给我们工程师带来了极大的便利，很多时候我们不需要精通工程问题中的力学原理，依然可以通过商业软件来解决问题，不过理论基础的缺失还是会给我们带来不少的困扰。随着动力有限元与柔性多体系统分析方法的成熟，有时候正确区分两者并不是很容易。机械领域应用比较广泛的有两类软件，一类是有限元软件，代表的有：ANSYS, NASTRAN, ABAQUS, LS-DYNA, Dytran 等；另一类是多体动力学软件，代表的有ADAMS, Recurdyn ， Simpack 等。在使用时，如何选用这两类软件并不难，但是如果深究这两类软件根本区别并不容易。例如，有限元软件可以分析静力学问题，也可以分析“动力学”问题，这里的“动力学”与多体动力学软件里面的动力学一样吗？有限元软件在分析动力学问题时，可以模拟物体的运动，它与多体动力学软件中模拟物体运动相同吗？多体动力学软件也可以分析柔性体的应力、应变等，这与有限元软件分析等价吗？ 1 有限元软件有限单元法是一种数学方法，不仅可以计算力学问题，还可以计算声学，热，磁等多种问题，我们这里只探讨有限元法在机械领域的应用。计算结构应力、应变等的力学基础是弹性力学，弹性力学亦称为弹性理论，主要研究弹性体在外力作用或温度变化等外界因素下所产生的应力、应变和位移，从而为工程结构或构件的强度、刚度设计提供理论依据和计算方法。也就是说用有限元软件分析力学问题时，是用有限元法计算依据弹性力学列出的方程。考虑下面这个问题，在()0t ，时间内给一个结构施加一个随时间变化的载荷()P t ，我们希望得到结构的应力分布，在刚刚施加载荷的时候，结构中的应力会有波动，应力场是变化的，但很久以后，应力场趋于稳定。如果我们想得到载荷施加很久以后，稳定的应力场分布，那么应该用静力学分析方法分析

第10章(非线性有限元) (1)分解

第10章非线性动力有限元法 (1) 10.1 几何非线性问题的有限元法 (2) 10.1.1 几何非线性问题的牛顿迭代法 ........................................................................... 2 10.1.2 典型单元的切线刚度矩阵 ................................................................................. 4 10.2 材料非线性问题的有限元法 (8) 10.2.1 弹/粘塑性问题的基本表达式 .............................................................................. 8 10.2.2 粘塑性应变增量和应力增量 ............................................................................... 9 10.2.3 弹/粘塑性平衡方程 ............................................................................................ 10 10.3 材料非线性问题的动力有限元法 ................................................................................ 11 10.4 应用举例 (14) 10.4.1 粘弹粘塑性动力有限元分析举例 ................................................................... 14 习题.. (15) 第10章非线性动力有限元法当机械结构受到较大的外载荷，或受到持续时间较短的冲击载荷作用时，结构会产生过大的变形, 以至于必须考虑结构几何大变形对结构整体刚度及固有频率的影响，即所谓的几何非线性影响。另外, 对于多数非线性动力学问题，还需要考虑材料非线性、接触非线性等方面的影响。非线性动力学分析求解的基本方程有如下形式 0=-+P I u M (4.141) 式中，Ku u C I += 为粘性效应项，考虑阻尼、粘塑、粘弹等效应。P 为外部激励。对于考虑各种非线性效应的动力学问题求解，需要对动力学方程进行直接时间积分。即非线性动力有限元分析具有如下特点：（1）问题分析过程需要考虑时间积分效应，不必做模态分析，不必提取固有频率；（2）采用直接积分方法求解非线性动力学方程，需要对时间作积分计算，因此计算量远远大于线性模态动力学方法；（3）非线性动力学分析中可以施加不同类型的载荷，包括结点力、非零位移、单元载荷；（4）在每个时间步上，进行质量、阻尼、及刚度的集成，采用完整矩阵，不涉及质量矩阵的近似；（5）可以同时考虑几何、材料和接触等多种非线性效应。非线性动力有限元分析程序常采用隐式Hilber-Hughes-Taylor 法进行时间积分运算。这种方法适于模拟非线性结构的动态问题，对于冲击、地震等激发的结构动态响应以及一些由于塑性或粘性阻尼造成的能量耗散，隐式算法特别有效。隐式积分方法需要对刚度矩阵求逆计算，并通过多次迭代求解增量步平衡方程。隐式Hilber-Hughes-Taylor 时间积分算法为无条件稳定，对时间步长没有特别的限制。采用子空间法也可以对动力学平衡方程作时间积分运算。子空间法是提取模态分析得到的各阶特征模态，并采用与线性模态动力学分析方法相近的分析方式进行求解。对于带有微小非线性效应的问题，如材料小范围进行入屈服、结点转角不大的情况，子空间法效率比进接积分法要高。

有限元动力学分析知识点汇总

复习目录一、模型输入、建模 A 输入几何模型 1、两种方法：No defeaturing 和 defeaturing （Merge合并选项、Solid实体选项、Small选项） 2、产品接口。输入IGES 文件的方法虽然很好，但是双重转换过程CAD > IGES > ANSYS 在很多情况下并不能实现100%的转换.ANSYS 的产品接口直接读入“原始”的CAD 文件，解决了上面提到的问题. 3、输入有限元模型。除了实体几何模型外， ANSYS 也可输入由某些软件包生成的有限元单元模型数据（节点和单元）。 B 实体建模 1、定义实体建模：建立实体模型的过程。（两种途径） 1）自上而下建模：首先建立体（或面）,对这些体或面按一定规则组合得到最终需要的形状. ?开始建立的体或面称为图元. ?工作平面用来定位并帮助生成图元. ?对原始体组合形成最终形状的过程称为布尔运算 ?总体直角坐标系 [csys,0] 总体柱坐标系[csys,1] 总体球坐标系[csys,2] 工作平面 [csys,4] 2）自下而上建模：按照从点到线，从线到面，从面到体的顺序建立模型。 B 网格划分 1、网格划分三步骤: 定义单元属性、指定网格的控制参数、生成网格 2、单元属性（单元类型 (TYPE)、实常数 (REAL)、材料特性 (MAT)） 3、单元类型单元类型是一个重要选项，它决定如下单元特性：自由度(DOF)设置、单元形状、维数、假设的位移形函数。 1）线单元（梁单元、杆单元、弹簧单元） 2）壳用来模拟平面或曲面。 3）二维实体用于模拟实体截面 4）三维实体 ?用于几何属性，材料属性，荷载或分析要求考虑细节，而无法采用更简单的单元进行建模的结构。

非线性有限元分析

非线性有限元分析 1 概述在科学技术领域，对于许多力学问题和物理问题，人们已经得到了它们所应遵循的基本方程（常微分方程或偏微分方程）和相应的定解条件（边界条件）。但能够用解析方法求出精确解的只是少数方程性质比较简单，并且几何形状相当规则的问题。对于大多数工程实际问题，由于方程的某些特征的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析的答案。这类问题的解决通常有两种途径。一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。但是这种方法只是在有限的情况下是可行的，因为过多的简化可能导致误差很大甚至是错误的解答。因此人们多年来一直在致力于寻找和发展另一种求解途径和方法——数值解法。特别是五十多年来，随着电子计算机的飞速发展和广泛应用，数值分析方法已成为求解科学技术问题的主要工具。已经发展的数值分析方法可以分为两大类。一类以有限差分法为代表，主要特点是直接求解基本方程和相应定解条件的近似解。其具体解法是将求解区域划分为网格，然后在网格的结点上用差分方程来近似微分方程，当采用较多结点时，近似解的精度可以得到改善。但是当用于求解几何形状复杂的问题时，有限差分法的精度将降低，甚至发生困难。另一类数值分析方法是首先建立和原问题基本方程及相应定解条件相等效的积分提法，然后再建立近似解法并求解。如果原问题的方程具有某些特定的性质，则它的等效积分提法可以归结为某个泛函的变分，相应的近似解法实际上就是求解泛函的驻值问题。诸如里兹法，配点法，最小二乘法，伽辽金法，力矩法等都属于这一类方法。但此类方法也只能局限于几何形状规则的问题，原因在于它们都是在整个求解区域上假设近似函数，因此，对于几何形状复杂的问题，不可能建立合乎要求的近似函数。 1960年，R.W.CLOUGH发表了有限单元法的第一篇文献“The Finite Element Method in Plane Stress Analysis”，这同时也标志着有限单元法（FEM）的问世。有限单元法的基本思想是将连续的求解区域离散为一组有限个，且按一定方式相互联接在一起的单元的组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模型化几何形状复杂的求解域。并且可以利用在每一个单元假设的近似函数来分片地表示全求解域上待求的未知场函数，从而使一个连续的无限自由度问题变成离散的有限自由度问题。现已证明，有限单元法是基于变分原理的里兹法的另一种形式，从而使里兹法分析的所有理论基础都适用于有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法。利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，而且事先不要求满足任何边界条件，因此可以用来处理很复杂的连续介质问题。在短短四十余年的时间里，有限单元的分析方法已经迅速地发展为适合于使用各种类型计算机解决复杂工程问题的一种相当普及的方法。如今，有限元广泛地应用于各个学科门类，已经成为工程师和科研人员用于解决实际工程问题，进行科学研究不可或缺的有力工具。有限单元法的应用围已由弹性力学平面问题扩展到空间问题，板壳问题，由静力平衡问题扩展到稳定问题，动力问题和波动问题。分析的对象从弹性材料扩展到塑性，粘弹性，粘塑性和复合材料等，从固体

第9章非线性问题的有限单元法

第9章非线性问题的有限单元法 9.1 非线性问题概述前面章节讨论的都是线性问题，但在很多实际问题中，线弹性力学中的基本方程已不能满足，需要用非线性有限单元法。非线性问题的基本特征是变化的结构刚度，它可以分为三大类：材料非线性、几何非线性、状态非线性。 1. 材料非线性(塑性, 超弹性, 蠕变) 材料非线性指的是材料的物理定律是非线性的。它又可分为非线性弹性问题和非线性弹塑性问题两大类。例如在结构的形状有不连续变化（如缺口、裂纹等）的部位存在应力集中，当外载荷到达一定数值时该部位首先进入塑性，这时在该部位线弹性的应力应变关系不再适用，虽然结构的其他大部分区域仍保持弹性。 2. 几何非线性(大应变, 大挠度, 应力刚化) 几何非线性是有结构变形的大位移引起的。例如钓鱼杆，在轻微的垂向载荷作用下，会产生很大的变形。随着垂向载荷的增加，杆不断的弯曲，以至于动力臂明显减少，结构刚度增加。 3. 状态非线性(接触, 单元死活) 状态非线性是一种与状态相关的非线性行为。例如，只承受张力的电缆的松弛与张紧；轴承与轴承套的接触与脱开；冻土的冻结与融化。这些系统的刚度随着它们状态的变化而发生显著变化。 9.2 非线性有限元问题的求解方法对于线性方程组，由于刚度方程是常数矩阵，可以直接求解，但对于非线性方程组，由于刚度方程是某个未知量的函数则不能直接求解。以下将简要介绍借助于重复求解线性方程组以得到非线性方程组解答的一些常用方法。 1.迭代法迭代法与直接法不同，它不是求方程组的直接解，而是用某一近似值代人，逐步迭代，使近似值逐渐逼近，当达到允许的规定误差时，就取这些近似值为方程组的解。与直接法相比，迭代法的计算程序较简单，但迭代法耗用的机时较直接法长。它不必存贮带宽以内的零元素，因此存贮量大大减少，且计算中舍入误差的积累也较小。以平面问题为例，迭代法的存贮量一般只需直接法的14左右。在求解非线性方程组时，一般采用迭代法。 2. 牛顿—拉斐逊方法 ANSYS程序的方程求解器计算一系列的联立线性方程来预测工程系统的响应。然而，非线性结构的行为不能直接用这样一系列的线性方程表示。需要一系列的带校正的线性近似来求解非线性问题。一种近似的非线性救求解是将载荷分成一系列的载荷增量，即逐步递增载荷和平衡迭代。它可以在几个载荷步内或者在一个载荷步的几个子步内施加载荷增量。在每一个增量的

独塔斜拉桥非线性动力特性有限元分析

– 12 – 现代物业?新建设 2012年第11卷第5期（1）式中：[M]、[C]、[K]分别为桥梁结构的质量、阻尼、刚度矩阵；、和δ分别为桥梁结构的加速度、速度和位移向量；F为作用于桥梁空间梁单元的力向量。求桥梁自振特性时，阻尼影响不大，一般不考虑其影响。令[C]=0、{F}=0,则得到其无阻尼自振方程，即（2）式(2)具有非零解的条件为: （3）式中：[M]、[K]含义同式子（1），ω2为振型的特征值(自振频率)。式（2）的特征方程(频率方程)为: （4）式中：[M]、[K]含义同式子（1），ωn 2为第n 阶振型的特征值（自振频率）；δn 为第n 阶振型向量，即主振型（模态）。对于式（3）求解广义特征值问题求解方法比较多，常用的有Lanczos向量迭代法、逆迭代法、Rayleigh- Ritz 法、Jacobi（雅可比）法、Ritz向量迭代法、子空间迭代法等。０前言斜拉桥发展至今已有四十多年的历史，斜拉桥的主跨跨径也不断增大。目前，世界著名海峡正准备建造主跨1,300～3,000m的斜拉桥。斜拉桥具有外观轻巧、跨度大、结构性能好、施工简便等优点，在国内外得到迅速发展。斜拉桥的动力特性是其结构动力分析及设计的前提,结构的动力特性取决于结构的组成体系、刚度、质量和支承条件等。因此建立一种理想的大跨度斜拉桥动力分析模型，并进行动力分析，掌握其特性，具有十分重要的意义。 1 分析方法及原理桥梁结构的振动特性主要取决于其各阶自振频率和主振型等。自振频率首先是表征结构刚性的指标，同时也是判断结构在动力作用下是否会发生车桥共振的依据。桥梁的动力方程可写为：现代建设 Modern Construction 独塔斜拉桥非线性动力特性有限元分析何军拥（广东工贸职业技术学院，广州　广东　510510）摘要：自振频率是评价桥梁动力性能的重要依据。结合某斜拉桥工程设计实例，考虑结构的非线性，采用大型有限元分析程序ANSYS建立三维有限元模型。以等效弹性模量法模拟拉索，利用主从约束模拟塔梁的连接及边墩、辅助墩与主梁的连接；主塔与主梁之间设纵向阻尼器，模拟纵向漂浮体系，对模型进行动力分析，其分析结果为桥梁维修和养护提供参考意见。关键词：独塔斜拉桥；ANSYS；动力特性中图分类号：U448.41 文献标识码：A 文章编号：1671-8089（2012）05-0012-03 Finite Element Analysis of Nonlinear Dynamic Characteristics of Single Pylon Cable-Stayed Bridge HE Jun-yong (Guangdong College of Industry and Commerce, Guangzhou 510510) Abstract: Natural frequencies are critical bases to evaluate the dynamic performance of bridges. This paper takes a large cable-stayed bridge as an example and a ? nite element computing model is set up for this bridge based on software ANSYS, which is to carry out modal analysis and get relevant modal parameters. The equivalent elastic modulus method is applied to analog cable modulus. The master-slave constraint is used to analog tower and beam connections, the side piers, auxiliary pier and beam connection, and longitudinal damper located between the main beam and the main tower, to simulate vertical ? oat system. The results may be used as reference in bridge inspection etc. Keywords: Single pylon cable-stayed bridges; ANSYS; Dynamic behavior [作者简介] 何军拥（1971- ），男，湖南邵阳人，硕士，广东工贸职业技术学院，副教授，高工，主要研究方向：工程力学、材料工程。

UG有限元分析教程

第1章高级仿真入门在本章中，将学习： ?高级仿真的功能。 ?由高级仿真使用的文件。 ?使用高级仿真的基本工作流程。 ?创建FEM和仿真文件。 ?用在仿真导航器中的文件。 ?在高级仿真中有限元分析工作的流程。 1.1综述 UG NX4高级仿真是一个综合性的有限元建模和结果可视化的产品，旨在满足设计工程师与分析师的需要。高级仿真包括一整套前处理和后处理工具，并支持广泛的产品性能评估解法。图1-1所示为一连杆分析实例。图1-1连杆分析实例高级仿真提供对许多业界标准解算器的无缝、透明支持，这样的解算器包括NX Nastran、MSC Nastran、ANSYS和ABAQUS。例如，如果结构仿真中创建网格或解法，则指定将要用于解算模型的解算器和要执行的分析类型。本软件使用该解算器的术语或“语言”及分析类型来展示所有网格划分、边界条件和解法选项。另外，还可以求解模型并直接在高级仿真中查看结果，不必首先导出解算器文件或导入结果。高级仿真提供基本设计仿真中需要的所有功能，并支持高级分析流程的众多其他功能。 ?高级仿真的数据结构很有特色，例如具有独立的仿真文件和FEM文件，这有利于在分布式工作环境中开发有限元（FE）模型。这些数据结构还允许分析师轻松地共享FE数据去执行多种类型分析。

UG NX4高级仿真培训教程 2 ?高级仿真提供世界级的网格划分功能。本软件旨在使用经济的单元计数来产生高质量网格。结构仿真支持完整的单元类型（1D、2D和3D）。另外，结构级仿真使分析师能够控制特定网格公差。例如，这些公差控制着软件如何对复杂几何体（例如圆角）划分网格。 ?高级仿真包括许多几何体简化工具，使分析师能够根据其分析需要来量身定制CAD几何体。例如，分析师可以使用这些工具提高其网格的整体质量，方法是消除有问题的几何体（例如微小的边）。 ?高级仿真中专门包含有新的NX传热解算器和NX流体解算器。 NX传热解算器是一种完全集成的有限差分解算器。它允许热工程师预测承受热载荷系统中的热流和温度。 NX流体解算器是一种计算流体动力学（CFD）解算器。它允许分析师执行稳态、不可压缩的流分析，并对系统中的流体运动预测流率和压力梯度，也可以使用NX传热和NX流体一起执行耦合传热/流体分析。 1.2仿真文件结构当向前通过高级仿真工作流时，将利用4个分离并关联的文件去存储信息。要在高级仿真中高效地工作，需要了解哪些数据存储在哪个文件中，以及在创建那些数据时哪个文件必须是激活的工作部件。这4个文件平行于仿真过程，如图1-2所示。图1-2仿真文件结构设计部件文件的理想化复制当一个理想化部件文件被建立时，默认有一.prt扩展名，fem#_i是对部件名的附加。例如，如果原部件是plate.prt，一个理想化部件被命名为plate_fem1_i.prt。一个理想化部件是原设计部件的一个相关复制，可以修改它。理想化工具让用户利用理想化部件对主模型的设计特征做改变。不修改主模型部件，

有限元动力学分析方程及解法

动力分析中平衡方程组的解法 1前言描述结构动力学特征的基本力学变量和方程与静力问题类似，但所有的变量都是时间的函数。基本变量三大类变量(,)i u t ξ、(,)ij t εξ和(,)ij t σξ是坐标位置(,,)x y z ξ和时间t 的函数，一般将其记为()()()i ij ij u t t t εσ。基本方程 (1) 平衡方程利用达朗贝尔原理将惯性力和阻尼力等效到静力平衡方程中，有 ,()()()()0ij j i i i t b t u t u t σρν+--=&&& (1) 其中ρ为密度，ν为阻尼系数。 (2) 几何方程 ,,1 ()(()())2ij i j j i t u t u t ε=+ (2) (3) 物理方程 ()()ij ijkl kl t D t σε= (3) 其中ijkl D 为弹性系数矩阵。 (4) 边界条件位移边界条件()BC u 为， ()()i i u t u t = 在u S 上 (4) 力的边界条件()BC p 为， ()()ij j i t n p t σ= 在p S 上 (5) 初始条件 0(,0)()i i u t u ξξ== (6) 0(,0)()i i u t u ξξ==&& (7)

虚功原理基于上述基本方程，可以写出平衡方程及力边界条件下的等效积分形式， ,() ()0p ij j i i i ij j i S u u b u d n p dA δσρνδσΩ∏=---+Ω+-=??&&& (8) 对该方程右端第一项进行分部积分，并应用高斯-格林公式，整理得， ()()0p ijkl ij kl i i i i i i i i S D u u u u d b u d p u dA εδερδνδδδΩΩ-++Ω-Ω+=???&&& (9) 有限元分析列式单元的节点位移列阵为， 111222()[(),(),(),(),(),()(),(),()]e t k k k U t u t v t w t u t v t w t u t v t w t =L (10) 单元内的插值函数为， (,)()()e t u t N U t ξξ= (11) 其中()N ξ为单元的形状函数矩阵，与相应的静力问题单元的形状函数矩阵完全相同，ξ为单元中的几何位置坐标。基于上面的几何方程和物理方程及(11)式，将相关的物理量表达为节点位移的关系，有， (,)[](,)[]()()()()e e t t t u t N U t B U t εξξξξ=?=?= (12) (,)()()()()e e t t t D DB U t S U t σξεξξ=== (13) (,)()()e t u t N U t ξξ=&& (14) (,)()()e t u t N U t ξξ=&&&& (15) 将(12)-(15)供稿到虚功方程(9)中，有， [()()()()]()0e e e e e e e T e t t t t t M U t C U t K U t R t U t δδ∏=++-=&&&g (16) 由于()e t U t δ具有任意性，消去该项并简写有， e e e e e t t t t U C U KU R ++=&&& (17) 其中， e e T M N Nd ρΩ= Ω? (18) e e T C N Nd νΩ=Ω? (19)

主流CAE有限元分析软件的比较

随着现代科学技术的发展，人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能，需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响，看看是否会发生破坏性事故；分析计算核反应堆的温度场，确定传热和冷却系统是否合理；分析涡轮机叶片内的流体动力学参数，以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式，这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析（FEA，Finite Element Analysis）方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中，有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃，主要表现在以下几个方面：增加设计功能，减少设计成本；缩短设计和分析的循环周期；增加产品和工程的可靠性；采用优化设计，降低材料的消耗或成本；在产品制造或工程施工前预先发现潜在的问题；模拟各种试验方案，减少试验时间和经费；进行机械事故分析，查找事故原因。在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分析计算，FEA 在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在，世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。以下对一些常用的软件进行一些比较分析： 1. LSTC公司的LS-DYNA系列软件 LS-DYNA是一个通用显式非线性动力分析有限元程序，最初是1976年在美国劳伦斯利弗莫尔国家实验室(Lawrence Livermore National Lab.)由J.O.Hallquist 主持开发完成的，主要目的是为核武器的弹头设计提供分析工具，后经多次扩充和改进，计算功能更为强大。此软件受到美国能源部的大力资助以及世界十余家著名数值模拟软件公司（如ANSYS、MSC.software、ETA等）的加盟，极大地加强了其的前后处理能力和通用性，在全世界范围内得到了广泛的使用。在软件的广告中声称可以求解各种三维非线性结构的高速碰撞、爆炸和金属成型等接触非线性、冲击载荷非线性和材料非线性问题。即使是这样一个被人们所称道的数值模拟软件，实际上仍在诸多不足，特别是在爆炸冲击方面，功能相对较弱，其欧拉混合单元中目前最多只能容许三种物质，边界处理很粗糙，在拉格朗日——欧拉结合方面不如DYTRAN灵活。虽然提供了十余种岩土介质模型，但每种模型都有不足，缺少基本材料数据和依据，让用户难于选择和使用。2. MSC.software公司的DYTRAN软件当前另一个可以计算侵彻与爆炸的商业通用软件是MSC.Software Corporation ( MSC公司) 的MSC.DYTR AN程序。该程序在是在LS-DYNA3D的框架下，在程序中增加荷兰PISCES INTERNATIONAL公司开发的PICSES的高级流体动力学和流体——结构相互作用功能，还在PISCES的欧拉模式算法基础上，开发了物质流动算法和流固耦合算法。在同类软件中，其高度非线性、流—固耦合方面有独特之处。MSC.DYTR AN的算法基本上可以概况为：MSC.DYTRAN采用基于Lagrange格式的有限单元方法（FEM）模拟结构的变形和应力，用基于纯Euler格式的有限体积方法（FVM）描述材料（包括气体和液体）流动，对通过流体与固体界面传递相互作用的流体—结构耦合分析，采用基于混合的Lagrange格式和纯Euler 格式的有限单元与有限体积技术，完成全耦合的流体-结构相互作用模拟。MSC.DYTRAN用有限体积法跟踪

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

结构非线性动力分析方法综述_周文峰

·自然科学研究· 结构非线性动力分析方法综述周文峰　郭　剑 (攀枝花学院土木工程学院,四川攀枝花　617000) 摘　要　时程分析法是一种计算机模拟分析方法,其优势在于能模拟出结构进入非弹性阶段的受力性能。该方法主要包括结构分析模型、单元模型和恢复力模型三个重要方面。本文从这三个方面简单介绍了结构非线性动力反应分析方法。关键词　非线性;动力分析;模型结构抗震设计方法经历了静力阶段、反应谱阶段和动力阶段。从本质上说,前二者所采用的方法均为静力法,且只能进行弹性分析。动力阶段的形成建立在计算机的普及和数值分析方法的出现基础之上,其分析方法称为时程分析法。时程分析法本质上是一种计算机模拟分析方法,能够计算出结构地震反应的全过程,该方法的突出优势在于能模拟出结构进入非弹性阶段的受力性能。时程分析法的出现促进了结构非线性地震反应分析的发展。它主要包括结构分析模型、单元模型和恢复力模型三个重要方面,下面从这三个方面进行简单介绍。 1　结构分析模型结构的模型化是非线性动力反应分析的第一步,结构模型的模拟应着重于其动力特性的模拟。因此体系恢复力、质量、阻尼模型的准确性是模拟精度的前提。目前的结构分析模型可分为以下几类: 1.1　层间模型考虑到框架结构质量的分布规律,很容易形成以楼层为单元的多质点体系的思路,故将这种模型称之为层间模型。在研究框架结构动力反应时,层间模型中采用得最多的是层间剪切型模型。该模型假定框架结构层间变形以剪切变形为主,忽略其它形式变形的影响,故而比较适用于高跨比不大、层数不多的框架。为了进一步拓宽此模型的适用范围,在此模型基础上又发展了层间剪弯型模型,使之能适用于层数较多和高跨比较大的框架。但是层间模型在实际使用中却存在比较大的困难,这主要反映在如何具体确定层间的剪切刚度及弯曲刚度的问题上,而且这二者之间又是耦合在一起的。这一问题层间模型自身是无法解决的。目前,层间模型只是对于常见的层数不多且平面布置十分简单、规则、对称并且能简化为平面结构的框架有一定的实用性,也就是说对于这类框架通常能根据经验进行适当的假设后进行简单推导得到层间单元刚度。 1.2　杆系模型杆系模型是将整体结构离散为梁、柱单元进行分析。杆系分析模型的出现不仅解决了层间模型所面临的层间刚度无法确定的困难,而且它还解决了层间模型所固有的另外两个缺陷。其一,如果说层间模型从宏观(层单元)角度展示了结构总体动力反应规律,那么由于框架各杆进入非弹性阶段的先后次序不同所造成的整个框架动力反应规律的不同,则是层间模型所不能解释、反映的。其二,无论从抗震研究还是设计角度来看,框架结构的梁、柱构件在地震作用下的反应规律到底如何也是人们所关心的,因为结构的设计最终要落实到构件的设计。如柱端弯矩增大系数应如何取值等,这些问题采用层间模型是无法回答的,从这个角度看也必须将框架结构细化到至少是构件层次才有可能解决这些问题。杆系分析模型分为两大类,平面杆系分析模型与空间杆系分析模型。目前,平面杆系分析模型的研究相对较为成熟,国内外已开始将注意力转向空间杆系分析模型的研究上。 2　单元模型对于杆系分析模型,目前用于模拟单元滞回性能的模型已有很多,这些单元分析模型可采取分类的方式加以比较考察。这些模型大致可分为两大类若干小类。 2.1　集中塑性铰模型单分量模型是集中塑性铰模型中最简单的一类,该模型将杆单元的非弹性性能用非线性弹簧反映,而不对非弹性变· 109·第23卷第4期攀枝花学院学报 2006年8月V o l .23.N o .4 J o u r n a l o f P a n z h i h u a U n i v e r s i t y A u g .2006

非线性有限元分析(学习总结报告)

非线性有限元博士研究生专业课课程报告

目录第一章绪言 (1) 1.1 非固体力学非线性问题的分类[1] (1) 1.2 非线性问题的分析过程[1] (2) 1.3 非线性有限元分析的基本原理 (2) 1.4 钢筋混凝土非线性分析的特点、现状及趋势 (3) 第二章非线性方程组的数值解法 (4) 2.1逐步增量法[3,4,5] (4) 2.2迭代法[3,4,5] (6) 2.3收敛标准 (8) 2.3.1.位移收敛准则 (8) 2.3.2.不平衡力收敛准则 (8) 2.3.3.能量收敛准则 (9) 2.4结构负刚度的处理[4,5] (9) 第三章材料的本构关系 (13) 3.1 钢筋的本构关系 (13) 3.1.1 单向加载下的应力应变关系 (13) 3.1.2 反复加载下的应力应变关系 (14) 3.2 混凝土的本构关系 (14) 3.2.1 单向加载下的应力应变关系 (14) 3.2.2 重复加载下的应力应变关系 (14) 3.2.3 反复加载下的应力应变关系 (14) 3.3 恢复力模型的分类 (14) 3.4 恢复力的获得方法 (15) 第四章非线性有限元在结构倒塌反应中的应用 (17) 4.1 钢筋混凝土结构倒塌反应研究现状 (17) 4.2 钢筋混凝土的有限元模型 (17) 4.2.1分离式模型 (18) 4.2.2组合式模型 (19) 4.2.3整体式模型 (20) 4.3 倒塌反应中RC结构有限元分析方法的选择 (20) 4.3.1隐式有限单元法 (21) 4.3.2显式有限单元法 (22) 4.4 钢筋混凝土框架结构的倒塌反应分析 (22) 4.4.1基于隐式有限单元法的倒塌分析 (22) 4.4.2 基于显式有限单元法的倒塌分析 (23) 4.5显式有限法在倒塌反应分析中的问题 (24)

第10章(非线性有限元)分解

公式号、图号等第十章非线性动力有限元法当机械结构受到较大的外载荷，或受到持续时间较短的冲击载荷作用时，结构会产生过大的变形, 以至于必须考虑结构几何大变形对结构整体刚度及固有频率的影响，即所谓的几何非线性影响。另外, 对于多数非线性动力学问题，还需要考虑材料非线性、接触非线性等方面的影响。非线性动力学分析求解的基本方程有如下形式 0=-+P I u M (4.141) 式中，Ku u C I += 为粘性效应项，考虑阻尼、粘塑、粘弹等效应。P 为外部激励。对于考虑各种非线性效应的动力学问题求解，需要对动力学方程进行直接时间积分。即非线性动力有限元分析具有如下特点：（1）问题分析过程需要考虑时间积分效应，不必做模态分析，不必提取固有频率；（2）采用直接积分方法求解非线性动力学方程，需要对时间作积分计算，因此计算量远远大于线性模态动力学方法；（3）非线性动力学分析中可以施加不同类型的载荷，包括结点力、非零位移、单元载荷；（4）在每个时间步上，进行质量、阻尼、及刚度的集成，采用完整矩阵，不涉及质量矩阵的近似；（5）可以同时考虑几何、材料和接触等多种非线性效应。非线性动力有限元分析程序常采用隐式Hilber-Hughes-Taylor 法进行时间积分运算。这种方法适于模拟非线性结构的动态问题，对于冲击、地震等激发的结构动态响应以及一些由于塑性或粘性阻尼造成的能量耗散，隐式算法特别有效。隐式积分方法需要对刚度矩阵求逆计算，并通过多次迭代求解增量步平衡方程。隐式Hilber-Hughes-Taylor 时间积分算法为无条件稳定，对时间步长没有特别的限制。采用子空间法也可以对动力学平衡方程作时间积分运算。子空间法是提取模态分析得到的各阶特征模态，并采用与线性模态动力学分析方法相近的分析方式进行求解。对于带有微小非线性效应的问题，如材料小范围进行入屈服、结点转角不大的情况，子空间法效率比进接积分法要高。此外，非线性动力有限元分析还可以采用显式动态算法，如中心差分法。显式时间积分算法为有条件稳定，其临界稳定时间步长限制了时间步长的大小，与有限元模型最小单元尺寸、材料应力波速等有关。显式时间积分法适于模拟高速冲击、接触等问题。上述方法的选择需要综合考虑计算量、分析问题的规模、单元限制等多方面因素，需要丰富的有限元模拟的理论、经验和实践知识。以下以几何非线性问题和材料非线性问题为例介绍非线性有限元法，其中粘弹粘塑性非线性材料问题的分析是典型的非线性动力有限元的求解思想。 9.1 几何非线性问题的有限元法几何非线性问题一般是指物体经历大的刚体位移和转动，但固连于物体坐标系中的应变分量仍假设为小量, 即大位移小应变情况。