有限元动力学分析知识点汇总
材料力学有限元分析知识点总结
材料力学有限元分析知识点总结材料力学是研究物质力学性质和行为的学科,而有限元分析是一种利用计算机数值模拟方法对工程问题进行分析和计算的技术。
本文将从理论基础、有限元建模、求解方法和误差分析等方面总结材料力学有限元分析的关键知识点。
一、理论基础1. 材料力学基本原理:包括应力、应变、变形和弹性模量等基本概念,以及胡克定律和应力应变关系等基本理论。
2. 有限元法基本原理:包括将实际结构离散为有限个单元,建立节点和单元之间的关系,以及应用物理原理和数值方法求解得到数值解的基本思想。
3. 有限元离散方法:包括将连续问题离散化为有限个子问题,建立单元刚度矩阵和全局刚度矩阵,以及应用有限元法进行力学问题分析的基本步骤。
二、有限元建模1. 几何建模:将实际工程结构进行几何建模,通常使用CAD软件进行建模,包括建立节点和单元等。
2. 材料建模:根据实际材料的物理性质和力学行为,选择适当的材料模型,如线性弹性模型或非线性材料模型。
3. 网格划分:将结构离散为有限个单元,通常使用三角形单元或四边形单元进行网格划分,确保离散后的单元足够小且保证几何形状的准确性。
三、求解方法1. 单元应力应变计算:通过数值方法计算每个单元的应力和应变,可采用解析解、数值积分或有限元法求解。
2. 节点位移计算:根据应力应变关系和单元的几何形状,计算每个节点的位移,从而得到结构的变形情况。
3. 刚度矩阵的建立:根据单元的几何形状、材料性质和节点位移等信息,建立单元刚度矩阵和全局刚度矩阵,用于力学方程的求解。
4. 边界条件的施加:根据实际工程问题,施加适当的边界条件,如固支约束和荷载条件等,从而得到合理的求解结果。
四、误差分析1. 收敛性分析:通过逐步增加单元数目或减小网格大小,观察求解结果是否趋近于稳定值,从而判断数值解的收敛性。
2. 精度分析:通过与解析解或实验结果进行比较,评估数值解的精度,包括位移误差、应力误差和能量误差等指标。
3. 稳定性分析:判断数值解的稳定性和可靠性,防止数值发散或出现明显的计算错误。
动力学有限元
6.2结构动力有限元法理论与模型一、基本原理在实际问题的求解中,应用最广的是基于位移的有限元素法。
此法的基本思想是把本来为连续的工程结构分割成在结点上相联的单元组合体。
取这些结点的位移为基本未知量,并假定每个单元中的位移用单元位移函数来描述,这实质上是假定了单元的模态。
在此基础上,利用能量变分原理进行单元分析的全结构分析,得到全结构的振动平衡方程,从而把连续体的动力学问题化为多自由度系统的振动问题。
有限元动力分析的基本过程是首先将工程结构离散化,通过选择合理的单元确定出分析模型,在此基础上选择位移函数,进行单元分析,确定单元的刚度、质量、阻尼、载荷矩阵,再经过坐标变换,通过能量变分原理,进行全结构分析,建立系统的振动平衡方程。
最后运用有限元数值方法进行方程的求解。
结构动力有限元法采用的单元位移函数与静力分析相同,基本原理和求解过程也与静力分析相同,不同之处仅在分析模型的确定与运动方程的建立方面。
二、动态分析模型的确定由于结构动态分析中除考虑弹性力外,还要考虑惯性力和阻尼力,其运动方程是常微分方程组,所以动态分析的复杂程度高,计算工作量大,有限元分析模型要尽量精炼、简单。
1.模型确定的基本原则•分析模型应与分析的目的相适应。
动力分析的目的各不相同,有的是为了提供固有特性计算动态响应或供控制系统用;有的是为了舱内提供振动环境。
不同的目的,通常要求不同的模态数与计算精度。
显然,用于估算基本固有频率的模型应当比计算冲击响应的模型简单。
用于设计计算的模型应当比用于校核计算的模型简单。
•分析模型要与选用的计算工具与计算条件相适应。
计算机软件种类日益丰富,选择分析模型要与所用程序、所用计算机容量相适应。
如对于容量大的计算机,可选用较为复杂的有限元模型,而对于容量小的计算机则在能反映结构动态性能的前提下尽量简化模型,使求解规模尽量小。
对于大模型,可选用子结构模型,采用模态综合方法求解。
应注意, 不一定模型愈精细精度就愈高。
有限元动力学分析知识点
有限元动力学分析知识点复习目录一、模型输入、建模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)三维实体✓用于几何属性,材料属性,荷载或分析要求考虑细节,而无法采用更简单的单元进行建模的结构。
✓也用于从三维CAD系统转化而来的几何模型,而这些几何模型转化成二维模型或壳体会花费大量的时间和精力4、单元阶次与形函数•单元阶次是指单元形函数的多项式阶次。
•什么是形函数?–形函数是指给出单元内结果形态的数值函数。
有限元动力学分析知识点汇总
复习目录一、模型输入、建模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)三维实体✓用于几何属性,材料属性,荷载或分析要求考虑细节,而无法采用更简单的单元进行建模的结构。
✓也用于从三维CAD系统转化而来的几何模型,而这些几何模型转化成二维模型或壳体会花费大量的时间和精力4、单元阶次与形函数•单元阶次是指单元形函数的多项式阶次。
•什么是形函数?–形函数是指给出单元内结果形态的数值函数。
有限元静力学及动力学分析
1.2静力学分析步骤
基于ANSYS进行静力分析的基本步骤与ANSYS典型分析的 过程相同,一般包括建模、加载求解和检查分析结果等3个基本 步骤。
1.建模 2.加载求解 3.检查分析结果
2.1动力学有限元分析原理
1) 动力学分析的原因 2) 动力学有限元分析引例 3) 动力学分析的定义和目的 4) 动力学分析类型
1) 动力学分析的原因
静力分析也许能确保一个结
构可以承受稳定载荷的条件,但 这些还远远不够,尤其在载荷随 时间变化时更是如此。
umax= 位移幅值 =位移函数的相位角 u1 = 实部, umaxcos u2 = 虚部, umaxsin
3、 ANSYS谐响应分析步骤
频率 影响
要输 密度
前处理:定类型,画模型,设属性,分网格。 求 解:添约束,加载荷,查错误,求结果。 后处理:列结果,绘图形,显动画,下结论。
读结果 列固有频率 绘振型 制作振型动画
i j
(t) (t)
Fi (t)e [K ]e{ (t)}e
(3)整体分析
(4)通用运动方程
1.结构离散 与静力分析相同,选用适当的单元类型将连
续的弹性体离散成有限多个单元和节点。 2.单元分析
从离散的弹性体中任 意取出一个单元。利用给 定的位移插值方式表示单 元 内 任 一 点 的 位 移 {δ(t)}e, 进而确定节点的速度和加 速度。 3.整体分析 利用各节点处的变形协调条件和动力平衡条件即达朗贝尔原理,建立整体刚 度方程;
1. 基本方程和术语
通用运动方程: 假定为自由振动并忽略阻尼:
有限元基本知识归纳
有限元知识点归纳1.、有限元解的特点、原因?答:有限元解一般偏小,即位移解下限性原因:单元原是连续体的一部分,具有无限多个自由度。
在假定了单元的位移函数后,自由度限制为只有以节点位移表示的有限自由度,即位移函数对单元的变形进行了约束和限制,使单元的刚度较实际连续体加强了,因此,连续体的整体刚度随之增加,离散后的刚度较实际的刚度K为大,因此求得的位移近似解总体上将小于精确解。
2、形函数收敛准则(写出某种单元的形函数,并讨论收敛性)P49(1)在节点i处N i=1,其它节点N i=0;(2)在单元之间,必须使由其定义的未知量连续;(3)应包含完全一次多项式;(4)应满足∑Ni=1以上条件是使单元满足收敛条件所必须得。
可以推证,由满足以上条件的形函数所建单元是完备协调的单元,所以一定是收敛的。
4、等参元的概念、特点、用时注意什么?(王勖成P131)答:等参元—为了将局部坐标中几何形状规则的单元转换成总体(笛卡尔)坐标中的几何形状扭曲的单元,以满足对一般形状求解域进行离散化的需要,必须建立一个坐标变换。
即:为建立上述的变换,最方便的方法是将上式表示成插值函数的形式,即:其中m是用以进行坐标变换的单元节点数,xi,yi,zi是这些结点在总体(笛卡尔)坐标内的坐标值,Ni’称为形状函数,实际上它也是局部坐标表示的插值函数。
称前者为母单元,后者为子单元。
还可以看到坐标变换关系式和函数插值表示式:在形式上是相同的。
如果坐标变换和函数插值采用相同的结点,并且采用相同的插值函数,即m=n,Ni’=Ni,则称这种变换为等参变换。
5、单元离散?P42答:离散化既是将连续体用假想的线或面分割成有限个部分,各部分之间用有限个点相连。
每个部分称为一个单元,连接点称为结点。
对于平面问题,最简单、最常用的离散方式是将其分解成有限个三角形单元,单元之间在三角形顶点上相连。
这种单元称为常应变三角形单元。
常用的单元离散有三节点三角形单元、六节点三角形单元、四节点四边形单元、八节点四边形单元以及等参元。
有限元动力学分析知识点
复习目录一、模型输入、建模A 输入几何模型1、两种方法:No defeaturing 和 defeaturing(Merge合并选项、Solid实体选项、Small选项)2、产品接口。
输入IGES 文件的方法虽然很好,但是双重转换过程CAD > IGES > ANSYS 在很多情况下并不能实现100%的转换.ANSYS 的产品接口直接读入“原始”的CAD 文件,解决了上面提到的问题.3、输入有限元模型。
除了实体几何模型外, ANSYS 也可输入由某些软件包生成的有限元单元模型数据(节点和单元)。
B 实体建模1、定义实体建模:建立实体模型的过程。
(两种途径)1)自上而下建模:首先建立体(或面),对这些体或面按一定规则组合得到最终需要的形状.开始建立的体或面称为图元.工作平面用来定位并帮助生成图元.对原始体组合形成最终形状的过程称为布尔运算1第1页总体直角坐标系 [csys,0] 总体柱坐标系[csys,1]总体球坐标系[csys,2] 工作平面 [csys,4]2)自下而上建模:按照从点到线,从线到面,从面到体的顺序建立模型。
B 网格划分1、网格划分三步骤:定义单元属性、指定网格的控制参数、生成网格2、单元属性(单元类型 (TYPE)、实常数 (REAL)、材料特性 (MAT))3、单元类型单元类型是一个重要选项,它决定如下单元特性:自由度(DOF)设置、单元形状、维数、假设的位移形函数。
1)线单元(梁单元、杆单元、弹簧单元)2)壳用来模拟平面或曲面。
3)二维实体用于模拟实体截面4)三维实体2第2页用于几何属性,材料属性,荷载或分析要求考虑细节,而无法采用更简单的单元进行建模的结构。
也用于从三维CAD系统转化而来的几何模型,而这些几何模型转化成二维模型或壳体会花费大量的时间和精力4、单元阶次与形函数•单元阶次是指单元形函数的多项式阶次。
•什么是形函数–形函数是指给出单元内结果形态的数值函数。
有限元-动力学分析
1.电路分析的模态分析?
2. 什么是一阶电路、二阶电路?对应于机械的振动的一阶和 阶?
机翼的模态分析
模态分析
一般而言,模态分析就是分析器件的谐振频率。模态分析 是谐响应分析、瞬态动力学分析、谱分析的起点。
任何物体都有自身的固有频率,也称特征频率,用系统方 程描述后就是矩阵的特征值。很多工程问题都要涉及系统特 征频率问题,一个目的是防止共振、自激振荡之类的事故发 生,历史上有名的事件就是,步兵按统一步伐过大桥,结果 把大桥震塌了。
瞬态分析
瞬态动力学分析比静力学分析更复杂,因为按“工程”时间计算, 瞬态动力学分析通常要占用更多的计算机资源和更多的人力。
瞬态动力学分析可采用三种方法:完全(Full)法、缩减 (Reduced)法及模态叠加法。 完全法采用完整的系统矩阵计算瞬态响应(没有矩阵缩减)。它是 三种方法中功能最强的,允许包括各类非线性特性(塑性、大变形、大 应变等)。 缩减法通过采用主自由度及缩减矩阵压缩问题规模。在主自由度处 的位移被计算出来后,ANSYS可将解扩展到原有的完整自由度集上。 缩减法的 优点 是比完全法快且开销小。缩减法的 缺点是初始解只计算 主自由度的位移,第二步进行扩展计算,得到完整空间上的位移、应力 等。 模态叠加法通过对模态分析得到的振型(特征值)乘上因子并求和 来计算结构的响应
并非所有类型的计算都可以分为这几种状态,比如热分析就 有就没有模态分析。
静态分析
结构静力分析是有限元方法中最常用的一个应用领域。在 相当长的一段时间内,机械结构的设计,主要采用经验设计 计算模型非常简单、粗糙,有的还根本无法计算。
有限元第六讲 动力学分析
5.1.2谐响应分析
谐响应分析是用于确定线性结构在承受随时间按正弦(简谐) 规律变化的载荷时的稳态响应的一种技术。分析的目的是计 算结构在几种频率下的响应并得到一些响应值(通常是位移) 对频率的曲线,从这些曲线上可找到“峰值”响应并进一步 查看峰值频率对应的应力。
这种分析技术只计算结构的稳态受追振动,发生在激励开 始时的瞬态振动不在谐响应分析中考虑。作为一种线性分析, 该分析忽略任何即使己定义的非线性特性,如塑性和接触 (间隙)单元。但可以包含非对称矩阵,如分析在流体一结构 相互作用问题。谐响应分析也可用于分析有预应力的结构, 如小提琴的弦(假定简谐应力比预加的拉伸应力小得多)
MassMatrix Formulation[LLIMPMIL]:使用该选项可以选 定采用默认的质量矩阵形成方式(和单元类型有关)或集中质 量阵近似方式,建议在大多数情况下应采用默认形成方式。
PrestressEffectsca/culation[PSTRES]:选用该选项可以计 算有预应力结构的模态。默认的分析过程不包括预应力,即 结构是处于无应力状态的。
求解结构的前几阶模态,以了解结构如何响应的情形。该方法采用集中 质量阵(LUMPM,ON); Reduced(Householder)method:使用减缩的系统矩阵求解,速度快。但 由于减缩质量矩阵识近似矩阵,所以相应精度较低; Unsymmetric method:用于系统矩阵为非对称矩阵的问题,例如流体一 结构相同作用; Damped method:用于阻尼不可忽略的问题; QR Damped method:采用减缩的阻尼阵计算复杂阻尼问题,所以比 Damped method方法有更快的计算速度和更好的计算效率。
有限元讲义 6-动力分析有限元
第6章 结构动力分析有限元法此前述及的问题属于静力分析问题,即作用在结构上的荷载是与时间无关的静力。
由此求得的位移、应力等均与时间无关。
实际工程中的大部分都可简化成静力问题。
但当动载与静载相比不容忽略时,一般应进行动力分析。
如地震作用下的房屋建筑,风荷载作用下的高层建筑等,都应计算动荷载作用下的动力反应。
研究课题中以动力问题为主。
解决动力问题有两大工作要做:一是动荷载的模拟和计算,二是结构反应分析。
本章将讨论如何用有限元来解决动力计算问题。
6.1 结构动力方程一.单元的位移、速度和加速度函数设单元的位移函数为;}{[]}{ef N d = 6—1—1式中:单元位移函数列阵}{f 、结点位移函数列阵}{ed 均是时间t 的函数。
由6-1-1可求得单元的速度、加速度函数:}{[]}{e f N d = 6—1—2 }{[]}{ef N d = 6—1—3二.单元的受力分析设图示三角形单元,当它处于运动状态时,其上的荷载一般应包括:单元上的荷载;单元对结点的作用力,}{[]}{(,eeix iy F F F K d ⋅⋅⋅=结点力)单元内部单位体积的:惯性力:}{}{[]}{em F f N dρρ=-=- 6—1—4阻尼力(设正比于运动速度):}{}{[]}{ecF f N d αραρ=-=- 6—1—5干扰力(已知的条件):}{p F根据达朗贝尔原理,上述四力将构成一瞬时平衡力系,使单元处于动平衡状态。
为此寻求四者之间的关系;三.结点力与结点位移、速度和加速度之间的关系用虚功原理推导:令单元结点发生任意可能的虚位移}{*d,它满足单元所定义的位移场,即虚位移场}{[]}{**f N d =成立。
作用在单元上的外力所作的外力虚功:}{}{}{}{}{}{}{}{****TTTTPcmvvvT dF f F dv f F dv f F dv =+++⎰⎰⎰单元内部应力在由于虚位移所引起的虚应变上所做的内力虚功:}{}{[]}{[][]}{**TTvW dv B d D B d dv εσ==⎰()根据虚功原理(T=W ),若将惯性力}{m F ,阻尼力}{c F 用上面的6—1—4,6—1—5代替,得:}{}{[]}{}{[]}{[]}{[]}{[]}{[]}{[][]}{*****TPvvTvVd F N d F dv N d N d dv N d N d dv B d D B d dvαρρ+--=⎰⎰⎰⎰TTT ()()()()由于虚位移的任意性,可从等式两边各项中消去}{*d T,得:}{[][][]}{[][]}{[][]}{[]}{TTpvvvvF B D B dv d N N dv d N N dv d N F dv αρρ=++-⎰⎰⎰⎰TT简写为:}{[]}{[]}{[]}{}{eF k d c d m d R =++- 6—1—6式中:[][][][]Tv k B D B dv =⎰ 单刚(第一项为弹性恢复力) [][][]v c N N dv αρ=⎰T单元阻尼矩阵(第二项为阻尼力) [][][]v m N N dv ρ=⎰T 质量矩阵(第三项为惯性力)[][][]R e P v N F dv =⎰T 包括由作用在单元上的干扰力转化成的等效结点荷载6—1—6即为单元结点力之间的关系式。
有限元-第9讲-动力学问题有限单元法
a1 ae a2
... an
ui(t) ai vi(t)
wi(t)
(i 1,2,...n,)
(3)形成系统的求解方程
••
•
M a(t)C a(t)K(ta )Q (t)
(1.8)
其中
••
•
a(t)和a(t)
分别是系统的结点加速度向量和结点速度向量,
M,C,K和Q(t)分别是系统的质量、阻尼、刚度和结点载荷向量。9
•
at
1 2t
att att
中心差分法的递推公式
(3.1) (3.2)
1 t2 M 2 1 tC a t t Q t K 2 t2 M a t 1 t2 M 2 1 tC a t t(3.3)
上式是求解各个离散时间点解的递推公式,这种数值积分方法又 称为逐步积分法。
动力分析的计算工作量很大,因此提高效率,节省计算工作量的 数值方案和方法是动力分析研究工作中的重要组成部分。目前两 种普遍应用的减缩自由度的方法是减缩法和动力子结构法。
11
第2节 质量矩阵和阻尼矩阵
一、协调质量矩阵和集中质量矩阵
单元质量矩阵
Me NTNdV称为协调质量矩阵。 Ve
集中质量矩阵假定单元的质量集中在结点上,这样得到的质量矩 阵是对角线矩阵。以下分实体单元和结构单元进行讨论。
16
第2节 质量矩阵和阻尼矩阵
按第二种方法计算,得到集中质量矩阵与第一种方法结果一样。
注:对于8结点矩形单元,两种方法得到的集中质量矩阵不同。
在实际分析中,更多的是推荐用第二种方法来计算集中质量矩阵。 2.结构单元
2结点经典梁单元、协调质量矩阵和集中质量矩阵如下所示: (1)协调质量矩阵
位移插值函数是 N N 1 N 2 N 3N 4(2.7)
PROE5.0动力学与有限元分析笔记
第一篇机构动力学分析机构动力学分析工作模式:FEM模式:是对模型进行网格划分、边界约束、载荷、理想处理等前处理,后需第三方软件进行求解。
集成模式:运行于PROE野火平台之上,操作界面与PROE野火相同,能直接使用PROE野火的参数进行分析和优化。
应用程序——Mechanica——确定。
独立模式:不需要PROE野火平台支持,能独立运行,可导入第三方软件模型。
应用程序——机构:包括运动仿真和动态分析运动仿真:1、定义运动副、伺服电机以实现运动模拟。
2、观察记录分析。
3、测量位置、速度、加速度等运动特征。
4、图形显示这些测量值。
5、创建轨迹和运动包络,用物理方法描述运动。
动态分析:1、使用机械动态功能在机构上定义重力、力和力矩、弹簧、阻尼等特征。
2、可以对机构设置材料、密度等基本属性特征,使其更加接近现实中的机构。
1、若不涉及质量、重力等基本属性参数,只需使用机械设计分析就能实现运动分析。
2、若受到重力、外力和力矩、阻尼等参数影响,必须使用机械设计进行静态分析、机械动态进行分析。
运动学分析流程:机构运动学仅讨论与刚体本身有关的因素,不讨论引起这些运动的因素(如重力、外力和摩擦力等)。
因此,运动学属空间和时间等基本概念及其导致的速度和加速度。
运动仿真就是机构运动学分析,它是不考虑作用于机构系统上的力的情况下分析机构运动,并对主体位置、速度和加速度进行测量。
运动仿真流程:创建模型——检查模型——添加模型化要素——准备进行分析——分析模型——获取分析结果2、检查模型:在装配模型中,拖动可以移动的零部件,观察装配连接情况。
3、添加模型化要素:在机构中添加伺服电动机等运动分析要素。
4、准备进行分析:定义初始位置,建立测量方式。
5、创建分析模型:对所创建的机构模型进行运动学分析6、获取结果:回放分析结果、零件之间的干涉检查、获取轨迹曲线和运动包络线。
动力学分析流程:机构动力学是运动学和力学的统称。
力学是处理作用在物体上的力。
有限元仿真分析动力学-explicit总结
动力学-abaqus/explict总结动力学分为: 线性动力学和非线性动力学。
Standard适合模拟与模型的振动频率相比响应周期较长的问题;explicit:适合于模拟高速动力学问题。
线性动力学在abaqus/standard中求解,是基于模态的分析方法。
应用有:模态动力学:在时域内计算结构的线性动力学响应;可以使用直接积分稳态动力学: 计算由谐波激励引起的动态响应,可以使用直接积分。
响应谱分析:计算运动过程中的峰值响应;随即响应分析:计算随即连续激励的响应,如地震波。
非线性动力学:需要对运动方程进行直接积分;abaqus/standard中使用newmark积分方法,是隐式非线性直接积分法(无条件稳定,可以使用任意的时间增量,并且解仍然是有界的)。
Abaqus/explicit使用二阶精度的中心差分法(该方法是条件稳定的,只有在时间增量小于一定的临界值时才能给出有界的解)。
下面对explicit使用过程中的一些细节作简要的总结。
1.Abaqus/explicit:提供两种方案定义接触:1.1 General contact: 通用接触。
一般在模型中存在多个部件或复杂的拓扑结构情况下使用,该功能强大,不需像在abaqus/standard一样定义相互作用的接触对,在abaqus/explicit里会自动搜索相互作用的接触。
ExamplesThe following input specifies that the contact domain is based on self-contact of an all-inclusive, automatically generated surface but that contact (including self-contact in any overlap regions) should be ignored between the all-inclusive, automatically generated surface and surface_2:*CONTACT*CONTACT INCLUSIONS, ALL EXTERIOR 或ALL ELEMENT BASED*CONTACT PROPERTY ASSIGNMENT,,prop_1 (以全局的方式重新制定属性)*alum_surf,steel_surf,prop_2 (局部修改)*alum_surf,alum_surf,prop_3 (局部修改)*CONTACT EXCLUSIONS (不包括surface_2), surface_2Either of the following methods can be used to exclude self-contact for surface_1 fromthe contact domain:*CONTACT EXCLUSIONSsurface_1,or*CONTACT EXCLUSIONSsurface_1, surface_11.2.接触问题中调整初始节点位置Abaqus/explicit不允许接触表面的初始过盈。
有限元知识点汇总
有限元知识点汇总第一章1、何为有限元法?其基本思想是什么?》有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法。
》基本思想:化整为零,化零为整2、为什么说有限元法是近似的方法,体现在哪里?》有限元法的基本思想是几何离散和分片插值;》用离散单元的组合来逼近原始结构,体现了几何上的近似;用近似函数逼近未知量在单元内的真实解,体现了数学上的近似;利用及问题的等效的变分原理建立有限元基本方程,又体现了明确的物理背景。
3、单元、节点的概念?》单元:把参数单元划分成网格,这些网格就称为单元。
》节点:网格间相互连接的点称为节点。
4、有限元法分析过程可归纳为几个步骤?》3大步骤;——结构离散化;——单元分析;——整体分析。
5、有限元方法分几种?本课程讲授的是哪一种?》有限元方法分3种;——位移法、力法、混合法。
》本课程讲授的:位移法6、弹性力学的基本变量是什么?何为几何方程、物理方程及虚功方程?弹性矩阵的特点?》弹性力学的基本变量是——{外力、应力、应变、位移}》几何方程——{描述弹性体应变分量及位移分量之间关系的方程}》物理方程——{描述应力分量及应变分量之间的关系}》虚功方程——{描述内力和外力的关系的方程}》弹性矩阵特点——{ }7、何为平面应力问题和平面应变问题?》平面应力问题——{满足(1)几何条件——所研究的是一根很薄的等厚度薄板,即一个方向上的几何尺寸远远小于其余两个面上的几何尺寸;(2)载荷条件——作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布,而在两板面上无外力作用}》平面应变问题——{满足(1)几何条件——所研究的是长柱体,即长度方向的尺寸远远大于横截面的尺寸,且横截面沿长度方向不变;(2)载荷条件——作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力}第二章1 何为结构的离散化?离散化的目的?何为有限元模型?答:⑴所谓离散化,是用假想的线或面将连续物体分割成由有限个单元组成的集合体。
有限元分析的力学基础
2.1 变形体的描述与变量定义
(1) 变形体
变形体:即物体内任意两点之间可发生相对移动。 有限元方法所处理的对象:任意变形体
38
(2) 基本变量的定义
可以用以下各类变量作为任意变形体的描述
量
因此,在材料确定的情况下,基本的力学变量应该有:
位移、应变、应力
39
目的:对弹性体中的位移、应力、应变进行 定义和表达,进而建立平衡方程、几何方程 和材料物理方程
+
∂τ yz
∂y
dy dxdz
dy 2
+τ
yz dxdz
dy 2
−
τ
zy
+
∂τ zy
∂z
dz dxdy
dz 2
−τ zydxdy
dz 2
=
0
全式除以dxdydz,合并相同的项,得
τ
yz
+
1 2
∂τ yz
∂y
dy
−τ
zy
−
1 2∂τ zy∂z Nhomakorabeadz
=
0
略去微量项,得 τ yz = τ zy
∑ MY = 0 τ zx = τ xz
各个方向上具有相同特性;
(4) 线性弹性假定:物体的变形与外来作用的关系是线性的, 外力去除后,物体可恢复原状;
(5) 小变形假定:物体变形远小于物体的几何尺寸,在建立方 程时,可以高阶小量(二阶以上)。
以上基本假定将作为问题简化的出发点。
41
2.3 基本变量的指标表达
指标记法的约定:
自由指标:在每项中只有一个下标出现,如 σ ij ,i,
∑MZ =0
τ xy = τ yx
有限元第五章 有限元动力学基本原理
第五章 有限元动力学分析基本原理
在前面的介绍中,我们均假设作用在弹性体(或结 构)上的载荷与时间无关,与此相应的,位移、应力 及应变等也都和时间无关,即前面介绍的全部内容皆 称结构静力学有限元方法。但工程实际中还存在着另 外一类载荷与时间有关的动载荷作用于结构或弹性体, 此时,相应的位移、应力、应变等都与时间有关,而 且必须考虑惯性力和加速度等因素,这类分析或问题, 成为动力学分析。 对于质点—弹簧系统的振动,大家比较熟悉,例如 一个自由度为n的质点—弹簧振系,其动平衡方程为
停止迭代 此时为低阶特性
2
1
( i 1)
(i 1)
三、机械结构固有频率与振型
2.矩阵迭代法
例题:已知一振动系统的质量矩阵、刚度矩阵用迭 代法计算其最高阶固有频率和振型。
1 0 0 3 2 0 M 0 2 0 K 2 5 3 0 0 3 0 3 3 1 1 1 解: 1 1 1.5 1.5 K 1 1.5 11 / 6
& & & M C K P
第五章 有限元动力学分析基本原理
上式中每一项的含义不同
& & M C 为阻尼力
K 为弹性力
对于单元体而言,可以得到类似的上述方程
e T N N dV V
于是,令e T V来自m N N dV
一、单元质量矩阵的计算
1.一致质量矩阵
e
m 的计算式是通式,并因为计算质量矩阵和刚度矩
阵使用的形状函数一致,因此被称为一致质量阵。
有限元分析理论基础-大全-超详细
应力的单元平均或节点平均处理方法
最简单的处理应力结果的方法是取相邻单元或围绕节点各单元应力的平均值。
• 1.取相邻单元应力的平均值 这种方法最常用于 3 节点三角形单元中。这种最简单而又相当实用的单元得
到的应力解在单元内是常数。可以将其看作是单元内应力的平均值,或是单元 形心处的应力。由于应力近似解总是在精确解上下振荡,可以取相邻单元应力
们的平均值作为该节点的最后应力值 ,即 i
i
1
m
m
e i
e 1
其中,1~m 是围绕在 i 节点周围的全部单元。取平均值时也可进行面积加权。
有限元法求解问题的基本步骤
1.结构离散化
对整个结构进行离散化,将其分割成若干个单元,单元间彼此通过节点相连;
2.求出各单元的刚度矩阵[K](e)
虚应力原理的力学意义:如果位移是协调的,则虚应力和虚边界约束反力在他们 上面所作的功的总和为零。反之,如果上述虚力系在他们上面所作的功的和为零,则 它们一定是满足协调的。所以,虚应力原理表述了位移协调的必要而充分条件。
虚应力原理可以应用于线弹性以及非线性弹性等不同的力学问题。但是必须指 出,无论是虚位移原理还是虚应力原理,他们所依赖的几何方程和平衡方程都是基于 小变形理论的,他们不能直接应用于基于大变形理论的力学问题。
虚位移原理是平衡方程和力的边界条件的等效积分的“弱”形式; 虚应力原理是几何方程和位移边界条件的等效积分“弱”形式。 虚位移原理的力学意义:如果力系是平衡的,则它们在虚位移和虚应变上所作的 功的总和为零。反之,如果力系在虚位移(及虚应变)上所作的功的和等于零,则它 们一定满足平衡方程。所以,虚位移原理表述了力系平衡的必要而充分条件。一般而 言,虚位移原理不仅可以适用于线弹性问题,而且可以用于非线性弹性及弹塑性等非 线性问题。
有限元分析-动力学分析
谐响应分析
谐响应分析主要用于分析持续的周期载荷在结构中产生的 持续周期响应,以及确定线性结构承受随时间按正弦规律变 化的载荷时的稳态响应。 建立模型时必须指定弹性模量和材料密度,谐响应分析假 定所施加的所有载荷随时间按简谐规律变化,一个简谐载荷 需要三条信息:幅值、相位角、载荷频率范围。
谐响应分析实例
模态分析
一般而言,模态分析就是分析器件的谐振频率。模态分析 是谐响应分析、瞬态动力学分析、谱分析的起点。 任何物体都有自身的固有频率,也称特征频率,用系统方 程描述后就是矩阵的特征值。很多工程问题都要涉及系统特 征频率问题,一个目的是防止共振、自激振荡之类的事故发 生,历史上有名的事件就是,步兵按统一步伐过大桥,结果 把大桥震塌了。 飞机翅膀的频率分析的用途? 飞机飞行时更要注意频率问题,避免与气流共振,风洞试 验就是测试这种力学结构问题。模态分析的目的是想办法提 高结构的特征频率,现在的手段就是改变、优化设计尺寸和设 法减小结构的质量。
瞬态分析
阻力和阻尼的区别?
阻力和阻尼的区别
1.阻尼的外延大,阻力的外延小 2. 系统的能量的减小——阻尼振动不都是因“阻力”引起 的,就机械振动而言,一种是因摩擦阻力生热,使系统的 机械能减小,转化为内能,这种阻尼叫摩擦阻尼;另一种 是系统引起周围质点的震动,使系统的能量逐渐向四周辐 射出去,变为波的能量,这种阻尼叫辐射阻尼。 3. 阻尼、临界阻尼、欠阻尼、无阻尼 4. 力学中的阻尼是什么?电学中的阻尼是什么?
静态分析
结构静力分析是有限元方法中最常用的一个应用领域。在 相当长的一段时间内,机械结构的设计,主要采用经验设计, 计算模型非常简单、粗糙,有的还根本无法计算。 我们大部分的计算分析目的主要是 1) 结构的最优方案设计:根据计算结果的分析和比较,按强 度、刚度和稳定性要求,对原方案进行修改补充,从而保证 合理的应力。 2) 分析结构损害原因:当结构件在工作中发生故障如裂纹、 断裂、磨损过大等缺陷时,可应用有限元分析,研究损害原 因,找出危险区域和部位,直到找到合理结构。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
复习目录一、模型输入、建模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)三维实体✓用于几何属性,材料属性,荷载或分析要求考虑细节,而无法采用更简单的单元进行建模的结构。
✓也用于从三维CAD系统转化而来的几何模型,而这些几何模型转化成二维模型或壳体会花费大量的时间和精力4、单元阶次与形函数•单元阶次是指单元形函数的多项式阶次。
•什么是形函数?–形函数是指给出单元内结果形态的数值函数。
因为FEA 的解答只是节点自由度值,需要通过形函数用节点自由度的值来描述单元内任一点的值。
–形函数根据给定的单元特性给出。
–每一个单元的形函数反映单元真实特性的程度,直接影响求解精度,这一点将在下边说明。
5、网格密度•有限单元法的基本原则是:单元数(网格密度)越多,所得的解越逼近真实值。
•然而,随单元数目增加,求解时间和所需计算机资源急剧增加。
6、单元属性种类1)实常数和截面特性。
实常数用于描述那些用单元几何形状不能完全确定的几何参数。
2)材料特性每一分析都需要输入一些材料特性:结构单元所需的杨氏模量,热单元所需的热传导系数KXX等。
7、控制网格密度•ANSYS 提供了多种控制网格密度的工具, 既可以是总体控制也可以是局部控制:1)总体控制(智能网格划分、总体单元、尺寸缺省尺寸)2)局部控制(关键点尺寸、线尺寸、面尺寸)8、有两种主要的网格划分方法: 自由划分和映射划分.+2•自由划分–无单元形状限制.–网格无固定的模式.–适用于复杂形状的面和体.•映射划分–面的单元形状限制为四边形,体的单元限制为六面体 (方块).–通常有规则的形式,单元明显成行.–仅适用于“规则的” 面和体, 如矩形和方块.•过渡网格划分这一选择是在六面体单元和四面体单元间的过渡区生成金字塔形单元,(“集两家之长.”将四面体和六面体网格很好地结合起来而不破坏网格的整体性)•扫掠划分9、网格拖拉当把一个面拖拉成一个体时, 您可以将面上的网格随同它一起拖拉, 得到一个已网格化的体. 这称为网格拖拉.第4章 ANSYS建模基本方法(耿老师)1、ANSYS建模方法•直接建模•实体建模•输入在计算机辅助设计系统中创建的实体模型•输入在计算机辅助设计系统创建的有限元模型2、直接建模•直接创建节点和单元,模型中没有实体•优点–适用于小型或简单的模型–可实现对每个结点和单元的编号完全控•缺点–需要人工处理的数据量大,效率低–不能使用自适应网格划分功能–不适合进行优化设计–容易出错3、实体建模–先创建由关键点、线段、面和体构成的几何模型,然后用网格划分,生成节点和单元–优点–适合于复杂模型,尤其适合3D实体建模–需要人工处理的数据量小,效率高–允许对节点和单元实施不同的几何操作–支持布尔操作–支持ANSYS优化设计功能–可以进行自适应网格划分–可以进行进行局部网格细化–便于修改和改进•缺点–有时需要大量CPU处理时间–对小型、简单的模型有时比较繁琐–在特定条件下可能会失败4、工作平面工作平面—是一个可以移动的二维参考平面用于定位和确定体素的方向。
5、ANSYS中坐标系分类•整体坐标系和局部坐标系(定位几何体作用)•节点坐标系✓定义节点自由度的方向✓定义节点结果数据方向•单元坐标系✓规定正交材料特性的方向✓规定所施加面力的方向✓规定单元结构数据的方向•显示坐标系✓定义几何体被列表后显示•结果坐标系✓用来列表、显示或在通用后处理中节点或者单元结果转换到一个特定坐标系。
6、网格划分方法:自由划分、映射划分、延伸划分、自适应划分。
第5章加载、求解和后处理1、选择命令Selecting 功能可以将模型的一部分从整体中分离出来,为下一步工作做准备。
操作一般包括3步:•选择子集•对其所选择的图元执行操作•重新激活整个模型2、组元(Components):作为选择功能的一个延伸,通过给选中的一组图元命名,即可创建组元,组元可保存在数据库中。
集合(Assembly):集合可以由一个或者多个集合和一个或多个其他组元组成。
3、静力分析与动力分析的区别•静力分析假设只有刚度力有效。
•动力分析考虑所有三种类型的力。
•如果施加的荷载随时间快速变化,则惯性力和阻尼力通常很重要。
•因此,可以通过判断载荷是否与时间相关,选择静力分析还是动力分析。
•如果在相对较长的时间内载荷是一个常数,选择静力分析。
•否则,选择动力分析。
•总之,如果激励频率小于结构最低价固有频率的 1/3,则可以进行静力分析。
4、线性分析和非线性分析的区别•线性分析假设忽略荷载对结构刚度变化的影响。
典型特征是:–小变形–应力、应变在线性弹性范围内–没有诸如两物体接触或分离时的刚度突变•如果加载引起结构刚度显著变化,必须进行非线性分析。
引起结构刚度显著变化的典型原因:–应变超过弹性范围 (塑性)–大变形,例如承载的钓鱼竿–两物体之间的接触5、载荷分类•自由度DOF -定义节点的自由度( DOF )值 (结构分析_位移) •集中载荷-点载荷 (结构分析_力)•面载荷-作用在表面的分布载荷 (结构分析_压力)•体积载荷-作用在体积或场域内 (热分析_体积膨胀)•惯性载荷-结构质量或惯性引起的载荷 (重力、角速度等) 6、载荷的施加直接在实体模型加载的优点:+几何模型加载独立于有限元网格。
重新划分网格或局部网格修改不影响载荷。
+加载的操作更加容易,尤其是在图形中直接拾取时。
*无论采取何种加载方式,ANSYS求解前都将载荷转化到有限元模型。
因此,加载到实体的载荷将自动转化到其所属的节点或单元上。
7、载荷步及时间选项•一个载荷步是指边界条件和载荷选项的一次设置,用户可对此进行一次或多次求解。
•一个分析过程可以包括:–单一载荷步(常常这是足够的)–多重载荷步•有三种方法可以用来定义并求解多载荷步–多次求解方法–载荷步文件方法–向量参数方法8、求解器选择及求解求解器的功能是求解关于结构自由度的联立线性方程组。
三种求解器:波前求解器、power求解器和稀疏矩阵求解器9、求解前的模型检查•统一的单位单元类型和选项•材料性质参数–考虑惯性时应输入材料密度–热应力分析时应输入材料的热膨胀系数•实常数 (单元特性)•单元实常数和材料类型的设置•实体模型的质量特性 (Preprocessor > Operate > Calc Geom Items)•模型中不应存在的缝隙•壳单元的法向•节点坐标系•集中、体积载荷•面力方向求解失败原因:•约束不够! (通常出现的问题)。
•当模型中有非线性单元,整体或部分结构出现崩溃或“松脱”。
•材料性质参数有负值•屈曲 - 当应力刚化效应为负(压)时,在载荷作用下整个结构刚度弱化。
如果刚度减小到零或更小时,求解存在奇异性,因为整个结构已发生屈曲。
10、ANSYS的两个后处理器•通用后处理器(即“POST1”) 只能观看整个模型在某一时刻的结果。
•时间历程后处理器 (即“POST26”) 可观看模型在不同时间的结果。
但此后处理器只能用于处理瞬态和/或动力分析结果。
11、结构分析常见的单元性能单元选择的基本准则:➢在结构分析中,结构的应力状态决定单元类型的选择。
➢选择维数最低的单元去获得预期的结果 (尽量做到能选择点而不选择线,能选择线而不选择平面,能选择平面而不选择壳,能选择壳而不选择三维实体).➢对于复杂结构,应当考虑建立两个或者更多的不同复杂程度的模型。
你可以建立简单模型,对结构承载状态或采用不同分析选项作实验性探讨。
1)线单元–Beam (梁)单元是用于螺栓(杆),薄壁管件,C形截面构件,角钢或者狭长薄膜构件(只有膜应力和弯应力的情况)等模型。
–Spar(杆)单元是用于弹簧,螺杆,预应力螺杆和薄膜桁架等模型。
–Spring单元是用于弹簧,螺杆,或细长构件,或通过刚度等效替代复杂结构等模型。
2)平面单元–X-Y 平面单元:单元定义在整体笛卡尔X-Y平面内(有限元模型必须建在此面内),分平面应力、平面应变或轴对称几种受力状态。
3)板壳单元–Shell (壳)单元用于薄面板或曲面模型。
壳单元分析应用的基本原则是每块面板的主尺寸不低于其厚度的 5~10 倍。
4)实体单元三维实体单元:–用于那些由于几何、材料、载荷或分析结果要求考虑的细节等原因造成无法采用更简单单元进行建模的结构。
–四面体模型使用CAD建模往往比使用专业的FEA分析建模更容易,也偶尔得到使用。
第6章动力学分析1①动力学绪论1、动力学分析定义动力学分析是用来确定惯性(质量效应)和阻尼起着重要作用时结构或构件动力学特性的技术。
2、动力学特性:振动特性、随时间变化载荷的效应、周期或者随机载荷的效应3、动力学分析的类型(1)模态分析来确定结构的振动特性(2)瞬态动力学分析来计算结构对随时间变化载荷的响应(3)谐分析来确定结构对稳态简谐载荷的响应(4)进行谱分析来确定结构对地震载荷的影响(5)随机振动分析来确定结构对随机震动的影响4、运动方程其中:[M] = 结构质量矩阵 [C] = 结构阻尼矩阵 [K] = 结构刚度矩阵{F} = 随时间变化的载荷函数 {u} = 节点位移矢量{} = 节点速度矢量{ü}= 节点加速度矢量5、不同分析类型是对这个方程的不同形式进行求解–模态分析:设定F(t)为零,而矩阵 [C] 通常被忽略;–谐响应分析:假设F(t)和 u(t)都为谐函数,例如 Xsin(wt),其中,X 是振幅, w 是单位为弧度/秒的频率;–瞬间动态分析:方程保持上述的形式。