有限元动力学分析知识点
有限元分析-动力学分析PPT课件

目录
• 引言 • 有限元分析基础 • 动力学分析基础 • 有限元分析在动力学中的应用 • 案例分析 • 结论与展望
01 引言
目的和背景
01
介绍有限元分析在动力学分析中 的应用和重要性。
02
阐述本课件的目标和内容,帮助 读者了解有限元分析在动力学分 析中的基本概念、方法和应用。
随着工程复杂性和精确度要求的提高,有限元分析在动力学分析中的 应用将更加重要和必要。
02
未来需要进一步研究有限元分析算法的改进和优化,以提高计算效率 和精度。
03
未来需要加强有限元分析与其他数值计算方法的结合,如有限差分、 有限体积等,以实现更复杂的动力学模拟和分析。
04
未来需要加强有限元分析在多物理场耦合和多尺度模拟中的应用,以 更好地解决工程实际问题。
有限元分析的优点和局限性
• 精确性:对于某些问题,可以得到相当精确的结 果。
有限元分析的优点和局限性
数值误差
由于离散化的近似性,结果存在一定的数值误 差。
计算成本
对于大规模问题,计算成本可能较高。
对模型简化的依赖
结果的准确性很大程度上依赖于模型的简化程度。
03 动力学分析基础
动力学简介
动力学是研究物体运 动过程中力与运动关 系的科学。
THANKS FOR WATCHING
感谢您的观看
ห้องสมุดไป่ตู้
求解等。
02 有限元分析基础
有限元方法概述
01
有限元方法是一种数值分析方法,通过将复杂的物理系统离散化为有 限个简单元(或称为元素)的组合,来模拟和分析系统的行为。
02
它广泛应用于工程领域,如结构分析、流体动力学、热传 导等领域。
有限元动力学分析知识点

有限元动力学分析知识点复习目录一、模型输入、建模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、单元阶次与形函数•单元阶次是指单元形函数的多项式阶次。
•什么是形函数?–形函数是指给出单元内结果形态的数值函数。
有限元静力学及动力学分析课件

03
操作步骤
利用有限元软件建立动力学模型, 进行瞬态模拟,将模拟结果与实
验结果进行对比分析。
02
实验设计
设计动力学实验,如自由落体冲 击实验,选用合适的实验设备和
试样。
04
结果分析
对比实验数据和模拟结果,评估 有限元分析方法在处理动力学问
题时的性能和准确性。
工程案例分析
案例背景
介绍汽车碰撞事故的背景,阐述有限元分析在汽车碰撞研 究中的重要性。
实验设计
设计简单的静力学实验,如悬 臂梁弯曲实验,准备相应的实
验设备和试样。
操作步骤
结果分析
利用有限元软件建立实验模型, 进行数值模拟,并将模拟结果
与实验结果进行对比分析。
通过对比实验数据和模拟结果, 评估有限元分析方法的精度和
适用性。
动力学实验验证
01
验证目的
通过动力学实验验证有限元分析 方法在处理动态问题时的准确性
模型建立
详细描述汽车碰撞有限元模型的建立过程,包括几何清理、 网格划分、材料属性赋值等步骤。
边界条件与求解设置
说明碰撞模拟中的边界条件,如初始速度、角度等,以及 求解器的选择和参数设置。
结果分析
展示碰撞过程中的变形、应力、应变等关键参数的变化情 况,并结合实验结果进行验证和讨论。最后,基于分析结 果提出汽车结构改进的建议。
自适应网格技术:结合并行计 算,实现自适应网格细化,以 在关键区域获得更精确的计算 结果,同时减少计算资源消耗。
通过这些高级有限元分析技术, 可以更准确、高效地模拟和分 析复杂工程问题,为设计和优 化提供有力支持。
PART 06
实验验证与案例分析
静力学实验验证
有限元 第9讲 动力学问题有限单元法

有限元第9讲动力学问题有限单元法动力学问题是指研究物体在运动中的受力和受力作用下的运动状态,常见的应用是结构工程学中的振动分析。
有限单元法是解决结构工程学中动力学问题的常用方法之一。
本文将介绍动力学问题和有限单元法的基本概念,并介绍其应用。
动力学问题的定义动力学是研究质点或刚体运动情况的分支学科,在结构工程学中是指结构在做振动时所受的力和运动状态。
动力学问题可以分为两种类型:稳态动力学问题和非稳态动力学问题。
稳态动力学问题是指结构在振动状态下所受的恒定力,而非稳态动力学问题则是指结构所受的变化的力,例如冲击力或地震力。
动力学问题的求解包括两个方面:一是确定受力情况;二是求解结构的运动状态。
确定受力情况通常需要通过实验或计算确定,求解结构运动状态则可以通过有限单元法来解决。
在结构工程学中,动力学问题的应用非常广泛。
例如,建筑物抗震设计需要对建筑物在地震作用下的反应进行分析,桥梁工程需要对桥梁在行车作用或风力作用下的振动响应进行分析。
有限单元法的基本概念有限单元法是一种将结构离散成若干小单元的数值分析方法,将结构分割成细小的单元,每个单元内部假设为均匀且连续的,通过对单元本身的运动状态进行求解,进而推知整个结构的运动状态。
有限元法用于解决的问题包括静力学问题、动力学问题、热力学问题和流体问题等。
有限单元法求解动力学问题的步骤主要包括如下几个步骤:1.离散化:将连续结构离散化成有限的小单元,每个单元内部运动状态通过定义一定数量的节点来确定。
2.建立单元的动力学方程:根据单元的形状和材料性质,建立单元的动力学方程,并计算单元的振动特性,例如频率和模态。
3.组装单元的方程:将单个单元的方程组装成整个结构的方程。
4.边界条件的处理:利用结构的边界条件(例如支撑、铰支等),将结构自由度减少到实际问题所需要的自由度。
5.求解结构的运动状态:通过求解整个结构的方程,得到结构的运动状态。
6.后处理:根据求解结果,进行结果的可视化和分析。
有限元静力学及动力学分析(第六章)

机械振动分析
对机械系统进行动力学分析,研 究其振动特性和稳定性,优化其 动态性能。
建筑结构地震响应
分析
采用有限元动力学分析方法,研 究建筑结构在地震作用下的响应, 评估其抗震性能和安全性。
05
有限元的优化设计
优化设计的基本概念
设计变量
01
在优化设计中需要改变的参数,如梁的截面尺寸、材料的弹性
模量等。
将连续的物理系统离散 化为有限个小的单元, 形成网格。划分网格是 有限元法的关键步骤, 直接影响计算结果的精 度和计算效率。
根据变分原理和加权余 量法,建立每个单元的 有限元方程,并将这些 方程组合成整体方程。
利用数值方法求解有限 元方程,得到每个节点 的位移和应力等结果。
对计算结果进行可视化 处理、分析和评估,为 工程设计和优化提供依 据。
结果后处理
对求解结果进行后处理,如绘制应力云图、生成 位移曲线等,以便进行结构分析和优化。
静力学分析的实例
桥梁结构静力学分析
通过建立桥梁的有限元模型,施加车 辆载荷等静力载荷,求解平衡方程, 得到桥梁的位移和应力分布,评估其 承载能力和安全性。
建筑结构静力学分析
通过建立建筑的有限元模型,施加风 载、雪载等静力载荷,求解平衡方程 ,得到建筑的位移和应力分布,评估 其稳定性和安全性。
THANKS
感谢观看
建立数学模型
根据物理问题和约束条件,建立有限元分 析的数学模型。
目标函数评估
根据有限元分析结果,评估目标函数的值 。
有限元分析
对建立的数学模型进行有限元分析,得到 各设计变量的响应。
优化设计的实例
飞机机翼的优化设计
通过改变机翼的截面尺寸和材料属性,使机翼的重量最小化 。
有限元动力学分析知识点

有限元动力学分析知识点复习目录一、模型输入、建模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. 基本方程和术语
通用运动方程: 假定为自由振动并忽略阻尼:
精选有限元法基础动力学问题讲义

q,t t
称为逐步积分法。
如果质量矩阵M是对角的,C也是对角或可以忽略
,则利用递推公式求解时不需求解方程,直接可得下
一时间步的预测值。
显示时间积分(Explicit Time Integral)
第十七页,共68页。
17
有限元法基础
12. 动力学问题
(2)当t=0时,需要 q0和 qt,因此必须用专门的起
团聚质量矩阵 假定质量集中在节点上,导出的质量 Lumped Mass 矩阵是对角线矩阵,可提高计算效率。
9
第九页,共68页。
有限元法基础
12. 动力学问题
团聚质量矩阵的计算方法
(1)
M
e l
中每一行主元等于
Me
中该行所有元素之和
Mle
ne
ij
k 1 0
Me
ik
i j i j
(2)
写成矩阵形式
i 1
u Nqe
u(x, y, z,t)
u
=
v(x,
y,
z,
t
)
w(x, y, z,t)
q1
qe
qn
ui (t)
qi
vi
(t
)
wi (t)
N [N1 N2 Ni Ni I33
Nn]
6
第六页,共68页。
有限元法基础
12. 动力学问题 (五)有限元方程
将插值函数代入Galerkin积分表达式,由 q的任
第二十页,共68页。
20
有限元法基础
12. 动力学问题 (5)中心差分的显示算法,适合于由冲击、碰撞、爆 炸类型的载荷引起的波传播问题的求解。
因为这些问题本身就是在初始扰动后,按一定的 波速C逐步在介质中传播。
有限元-动力学分析

1.电路分析的模态分析?
2. 什么是一阶电路、二阶电路?对应于机械的振动的一阶和 阶?
机翼的模态分析
模态分析
一般而言,模态分析就是分析器件的谐振频率。模态分析 是谐响应分析、瞬态动力学分析、谱分析的起点。
任何物体都有自身的固有频率,也称特征频率,用系统方 程描述后就是矩阵的特征值。很多工程问题都要涉及系统特 征频率问题,一个目的是防止共振、自激振荡之类的事故发 生,历史上有名的事件就是,步兵按统一步伐过大桥,结果 把大桥震塌了。
瞬态分析
瞬态动力学分析比静力学分析更复杂,因为按“工程”时间计算, 瞬态动力学分析通常要占用更多的计算机资源和更多的人力。
瞬态动力学分析可采用三种方法:完全(Full)法、缩减 (Reduced)法及模态叠加法。 完全法采用完整的系统矩阵计算瞬态响应(没有矩阵缩减)。它是 三种方法中功能最强的,允许包括各类非线性特性(塑性、大变形、大 应变等)。 缩减法通过采用主自由度及缩减矩阵压缩问题规模。在主自由度处 的位移被计算出来后,ANSYS可将解扩展到原有的完整自由度集上。 缩减法的 优点 是比完全法快且开销小。缩减法的 缺点是初始解只计算 主自由度的位移,第二步进行扩展计算,得到完整空间上的位移、应力 等。 模态叠加法通过对模态分析得到的振型(特征值)乘上因子并求和 来计算结构的响应
并非所有类型的计算都可以分为这几种状态,比如热分析就 有就没有模态分析。
静态分析
结构静力分析是有限元方法中最常用的一个应用领域。在 相当长的一段时间内,机械结构的设计,主要采用经验设计 计算模型非常简单、粗糙,有的还根本无法计算。
有限元静力学及动力学分析课件

网格类型
一维、二维、三维网格, 以及六面体、四面体、四 边形等形状的网格。
网格质量
对计算结果的精度和稳定 性有重要影响,需要保证 网格质量良好。
材料属性定义
材料属性
弹性模量、泊松比、密度、热膨胀系数等。
材料属性赋值
根据实际材料属性赋予有限元模型相应的值。
材料非线性
考虑材料在不同应力应变状态下的非线性行为。
03
有限元动力学分析基础
动力学基本概念
01
02
03
04
动力学
研究物体运动和力之间关系的 科学。
牛顿第二定律
物体运动加速度与作用力成正 比,与物体质量成反比。
动能
物体由于运动而具有的能量。
势能
物体由于位置或形变而具有的 能量。
有限元动力学方程
拉格朗日方程
描述系统运动状态的微分方程。
哈密顿原理
最小作用量原理的一种形式,用于确定系统的运动轨迹。
有限元分析的历史与发展
有限元分析的思想起源于20世纪40年代,但直到20世纪60年代 才由Clough提出并命名为“有限元法”。
随着计算机技术的发展,有限元分析得到了广泛的应用和推广, 逐渐成为工程领域的重要工具。
近年来,随着计算能力的提高和数值算法的发展,有限元分析在 精度、稳定性和适用范围等方面得到了显著提升,能够处理更加 复杂和大规模的问题。
01
刚度矩阵的定义和 性质
描述刚度矩阵的物理意义、计算 方法和特性,以及它在建立有限 元方程中的作用。
02
载荷向量的定义和 计算
介绍载荷向量的概念、计算方法 和作用,以及它在建立有限元方 程中的作用。
03
边界条件的处理
描述如何将边界条件引入有限元 方程中,以及常见的边界条件类 型。
有限元分析-动力学分析

1.为何傅里叶变换要换成正弦函数余弦函数这样的三角级数? 2. 谐振运动的特征是什么?谐振运动有阻尼存在吗?
梁结构瞬态动力学分析实例
A steel beam of length and geometric properties shown in Problem Specifications is supporting a concentrated mass, m. The beam is subjected to a dynamic load F(t) with a rise time tr and a maximum value F1. If the weight of the beam is considered to be negligible, determine the time of maximum displacement response tmax and the response ymax. Also determine the maximum bending stress σbend in the beam.
谱分析
谱分析是一种将模态分析结果与已知的谱分析联系起来的 计算位移和应力的分析技术。它主要用于时间历程分析,以 便确定结构在任意时间变化载荷下的动力学响应,简单而言 就是载荷的谱不再是简谐运动。
简支梁的两端作垂直运动,也就是地震时的作用,确定其 响应频率。
梁对地基地震时的谱分析
A simply supported beam of length , mass per unit length m, and section properties shown in Problem Specifications, is subjected to a vertical motion of both supports. The motion is defined in terms of a seismic displacement response spectrum. Determine the nodal displacements, reactions forces, and the element solutions.
有限元基础知识归纳

有限元知识点归纳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答:离散化既是将连续体用假想的线或面分割成有限个部分,各部分之间用有限个点相连。
每个部分称为一个单元,连接点称为结点。
对于平面问题,最简单、最常用的离散方式是将其分解成有限个三角形单元,单元之间在三角形顶点上相连。
这种单元称为常应变三角形单元。
常用的单元离散有三节点三角形单元、六节点三角形单元、四节点四边形单元、八节点四边形单元以及等参元。
有限元第六讲 动力学分析

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)
有限元仿真分析动力学-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. 动力响应分析方法动力响应分析是结构动力学分析的核心内容,主要研究结构在外界动力荷载作用下的响应情况。
这种分析方法可以帮助工程师评估结构的动力性能,如位移、加速度和应力等。
三、有限元方法在结构动力学中的应用有限元方法在结构动力学分析中的应用广泛,可以模拟各种结构的动态响应。
例如,有限元方法可以用于分析建筑物在地震作用下的响应,以评估结构的抗震性能。
此外,有限元方法还可以用于模拟机械设备、桥梁和航天器等工程结构在振动荷载下的响应。
在使用有限元方法进行结构动力学分析时,需要注意选择适当的数学模型和边界条件,并合理选择有限元单元的类型和尺寸。
有限元第五章 有限元动力学基本原理

第五章 有限元动力学分析基本原理
在前面的介绍中,我们均假设作用在弹性体(或结 构)上的载荷与时间无关,与此相应的,位移、应力 及应变等也都和时间无关,即前面介绍的全部内容皆 称结构静力学有限元方法。但工程实际中还存在着另 外一类载荷与时间有关的动载荷作用于结构或弹性体, 此时,相应的位移、应力、应变等都与时间有关,而 且必须考虑惯性力和加速度等因素,这类分析或问题, 成为动力学分析。 对于质点—弹簧系统的振动,大家比较熟悉,例如 一个自由度为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 的计算式是通式,并因为计算质量矩阵和刚度矩
阵使用的形状函数一致,因此被称为一致质量阵。
结构动力学问题的有限元法

二、单元分析
单元分析旳任务仍是建立单元特征矩阵,形成单元特征方程。 动态分析中,单元特征矩阵:刚度矩阵、质量矩阵和阻尼矩阵。
动态分析中,仍采用虚位移原理建立单元特征矩阵。
在动载荷作用下,对于任一瞬时,设单元节点发生虚位移 qe ,则单元 内也产生相应旳虚位移 d 和虚应变 。单元内产生旳虚应变能为:
式中,ω为简谐振动圆频率;{Φ}为节点振幅列向量。
将解代入振动方程中,同步消去因子ejωt,可得
K 2 M 0
上式为一广义特征问题。根据线性代数可知,求解该问题能够求出n个特
征值
12
,
12
,,
2 n
和相相应旳n个特征向量
1,2 ,n 。其中特
征值ωi(i=1,2,…..,n)就是构造旳i阶固有频率,特征向量{Φi} i(i=1,2,…..,n)就是构造
三、总体矩阵集成 总体矩阵集成旳任务是将各单元特征矩阵装配成整个构造旳特征矩阵,
从而建立整体平衡方程,即
M q Cq K q Rt
式中,{q}为所以节点位移分量构成旳n阶列阵,n为构造总自由度数;
R
t
n
Ri
t
(i为节点数),称为节点载荷列阵;[K]、[M]、[C]
分别为构i造1 旳刚度矩阵、质量矩阵和阻尼矩阵。
旳i阶模态振型。
振型{Φi}是构造按频率ωi振动时各自由度方向振幅间旳相对百分比关系, 它反应了构造振动旳形式,并不是振幅旳绝对大小。
固有特征分析实际上就是求解广义特征值问题。求解旳数值措施主要有 1、变换法 基本思想是经过一系列矩阵变换,将矩阵[M][K]化为对角阵,
k11
K d
k 22
5 动态分析有限元法
- 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)自上而下建模:首先建立体(或面),对这些体或面按一定规则组合得到最终需要的形状.开始建立的体或面称为图元.工作平面用来定位并帮助生成图元.对原始体组合形成最终形状的过程称为布尔运算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、单元阶次与形函数•单元阶次是指单元形函数的多项式阶次。
•什么是形函数–形函数是指给出单元内结果形态的数值函数。
因为FEA 的解答只是节点自由度值,需要通过形函数用节点自由度的值来描述单元内任一点的值。
–形函数根据给定的单元特性给出。
–每一个单元的形函数反映单元真实特性的程度,直接影响求解精度,这一点将在下边说明。
5、网格密度•有限单元法的基本原则是:单元数(网格密度)越多,所得的解越逼近真实值。
•然而,随单元数目增加,求解时间和所需计算机资源急剧增加。
3第3页6、单元属性种类1)实常数和截面特性。
实常数用于描述那些用单元几何形状不能完全确定的几何参数。
2)材料特性每一分析都需要输入一些材料特性:结构单元所需的杨氏模量,热单元所需的热传导系数KXX等。
7、控制网格密度•ANSYS 提供了多种控制网格密度的工具, 既可以是总体控制也可以是局部控制:1)总体控制(智能网格划分、总体单元、尺寸缺省尺寸)2)局部控制(关键点尺寸、线尺寸、面尺寸)8、有两种主要的网格划分方法: 自由划分和映射划分.+2•自由划分–无单元形状限制.–网格无固定的模式.–适用于复杂形状的面和体.•映射划分4第4页–面的单元形状限制为四边形,体的单元限制为六面体 (方块).–通常有规则的形式,单元明显成行.–仅适用于“规则的” 面和体, 如矩形和方块.•过渡网格划分这一选择是在六面体单元和四面体单元间的过渡区生成金字塔形单元,(“集两家之长.”将四面体和六面体网格很好地结合起来而不破坏网格的整体性)•扫掠划分9、网格拖拉当把一个面拖拉成一个体时, 您可以将面上的网格随同它一起拖拉, 得到一个已网格化的体. 这称为网格拖拉.第4章 ANSYS建模基本方法(耿老师)1、ANSYS建模方法•直接建模•实体建模•输入在计算机辅助设计系统中创建的实体模型5第5页•输入在计算机辅助设计系统创建的有限元模型2、直接建模•直接创建节点和单元,模型中没有实体•优点–适用于小型或简单的模型–可实现对每个结点和单元的编号完全控•缺点–需要人工处理的数据量大,效率低–不能使用自适应网格划分功能–不适合进行优化设计–容易出错3、实体建模–先创建由关键点、线段、面和体构成的几何模型,然后用网格划分,生成节点和单元–优点–适合于复杂模型,尤其适合3D实体建模6第6页–需要人工处理的数据量小,效率高–允许对节点和单元实施不同的几何操作–支持布尔操作–支持ANSYS优化设计功能–可以进行自适应网格划分–可以进行进行局部网格细化–便于修改和改进•缺点–有时需要大量CPU处理时间–对小型、简单的模型有时比较繁琐–在特定条件下可能会失败4、工作平面工作平面—是一个可以移动的二维参考平面用于定位和确定体素的方向。
5、ANSYS中坐标系分类7第7页v1.0 可编辑可修改•整体坐标系和局部坐标系(定位几何体作用)•节点坐标系定义节点自由度的方向定义节点结果数据方向•单元坐标系规定正交材料特性的方向规定所施加面力的方向规定单元结构数据的方向•显示坐标系定义几何体被列表后显示•结果坐标系用来列表、显示或在通用后处理中节点或者单元结果转换到一个特定坐标系。
6、网格划分方法:自由划分、映射划分、延伸划分、自适应划分。
第5章加载、求解和后处理1、选择命令8第8页Selecting 功能可以将模型的一部分从整体中分离出来,为下一步工作做准备。
操作一般包括3步:•选择子集•对其所选择的图元执行操作•重新激活整个模型2、组元(Components):作为选择功能的一个延伸,通过给选中的一组图元命名,即可创建组元,组元可保存在数据库中。
集合(Assembly):集合可以由一个或者多个集合和一个或多个其他组元组成。
3、静力分析与动力分析的区别•静力分析假设只有刚度力有效。
•动力分析考虑所有三种类型的力。
•如果施加的荷载随时间快速变化,则惯性力和阻尼力通常很重要。
•因此,可以通过判断载荷是否与时间相关,选择静力分析还是动力分析。
9第9页•如果在相对较长的时间内载荷是一个常数,选择静力分析。
•否则,选择动力分析。
•总之,如果激励频率小于结构最低价固有频率的 1/3,则可以进行静力分析。
4、线性分析和非线性分析的区别•线性分析假设忽略荷载对结构刚度变化的影响。
典型特征是:–小变形–应力、应变在线性弹性范围内–没有诸如两物体接触或分离时的刚度突变•如果加载引起结构刚度显著变化,必须进行非线性分析。
引起结构刚度显著变化的典型原因:–应变超过弹性范围 (塑性)–大变形,例如承载的钓鱼竿–两物体之间的接触5、载荷分类•自由度DOF -定义节点的自由度( DOF )值 (结构分析_位移)10第10页•集中载荷-点载荷 (结构分析_力)•面载荷-作用在表面的分布载荷 (结构分析_压力)•体积载荷-作用在体积或场域内 (热分析_体积膨胀)•惯性载荷-结构质量或惯性引起的载荷 (重力、角速度等) 6、载荷的施加直接在实体模型加载的优点:+几何模型加载独立于有限元网格。
重新划分网格或局部网格修改不影响载荷。
+加载的操作更加容易,尤其是在图形中直接拾取时。
*无论采取何种加载方式,ANSYS求解前都将载荷转化到有限元模型。
因此,加载到实体的载荷将自动转化到其所属的节点或单元上。
7、载荷步及时间选项•一个载荷步是指边界条件和载荷选项的一次设置,用户可对此进行一次或多次求解。
•一个分析过程可以包括:–单一载荷步(常常这是足够的)11第11页–多重载荷步•有三种方法可以用来定义并求解多载荷步–多次求解方法–载荷步文件方法–向量参数方法8、求解器选择及求解求解器的功能是求解关于结构自由度的联立线性方程组。
三种求解器:波前求解器、power求解器和稀疏矩阵求解器9、求解前的模型检查•统一的单位单元类型和选项•材料性质参数–考虑惯性时应输入材料密度–热应力分析时应输入材料的热膨胀系数•实常数 (单元特性)•单元实常数和材料类型的设置12第12页•实体模型的质量特性 (Preprocessor > Operate > Calc Geom Items)•模型中不应存在的缝隙•壳单元的法向•节点坐标系•集中、体积载荷•面力方向求解失败原因:•约束不够! (通常出现的问题)。
•当模型中有非线性单元,整体或部分结构出现崩溃或“松脱”。
•材料性质参数有负值•屈曲 - 当应力刚化效应为负(压)时,在载荷作用下整个结构刚度弱化。
如果刚度减小到零或更小时,求解存在奇异性,因为整个结构已发生屈曲。
10、ANSYS的两个后处理器•通用后处理器(即“POST1”) 只能观看整个模型在某一时刻13第13页v1.0 可编辑可修改的结果。
•时间历程后处理器 (即“POST26”) 可观看模型在不同时间的结果。
但此后处理器只能用于处理瞬态和/或动力分析结果。
11、结构分析常见的单元性能单元选择的基本准则:在结构分析中,结构的应力状态决定单元类型的选择。
选择维数最低的单元去获得预期的结果 (尽量做到能选择点而不选择线,能选择线而不选择平面,能选择平面而不选择壳,能选择壳而不选择三维实体).对于复杂结构,应当考虑建立两个或者更多的不同复杂程度的模型。
你可以建立简单模型,对结构承载状态或采用不同分析选项作实验性探讨。
1)线单元–Beam (梁)单元是用于螺栓(杆),薄壁管件,C形截面构件,角钢或者狭长薄膜构件(只有膜应力和弯应力的情况)等模型。
–Spar(杆)单元是用于弹簧,螺杆,预应力螺杆和薄膜桁架等模型。
–Spring单元是用于弹簧,螺杆,或细长构件,或通过刚度等14第14页效替代复杂结构等模型。
2)平面单元–X-Y 平面单元:单元定义在整体笛卡尔X-Y平面内(有限元模型必须建在此面内),分平面应力、平面应变或轴对称几种受力状态。
3)板壳单元–Shell (壳)单元用于薄面板或曲面模型。
壳单元分析应用的基本原则是每块面板的主尺寸不低于其厚度的 5~10 倍。
4)实体单元三维实体单元:–用于那些由于几何、材料、载荷或分析结果要求考虑的细节等原因造成无法采用更简单单元进行建模的结构。
–四面体模型使用CAD建模往往比使用专业的FEA分析建模更容易,也偶尔得到使用。
第6章动力学分析1①动力学绪论1、动力学分析定义15第15页v1.0 可编辑可修改动力学分析是用来确定惯性(质量效应)和阻尼起着重要作用时结构或构件动力学特性的技术。
2、动力学特性:振动特性、随时间变化载荷的效应、周期或者随机载荷的效应3、动力学分析的类型(1)模态分析来确定结构的振动特性(2)瞬态动力学分析来计算结构对随时间变化载荷的响应(3)谐分析来确定结构对稳态简谐载荷的响应(4)进行谱分析来确定结构对地震载荷的影响(5)随机振动分析来确定结构对随机震动的影响4、运动方程其中:[M] = 结构质量矩阵 [C] = 结构阻尼矩阵 [K] = 结构刚度矩阵{F} = 随时间变化的载荷函数 {u} = 节点位移矢量{} = 节点速度矢量{ü}= 节点加速度矢量5、不同分析类型是对这个方程的不同形式进行求解–模态分析:设定F(t)为零,而矩阵 [C] 通常被忽略;16第16页–谐响应分析:假设F(t)和 u(t)都为谐函数,例如 Xsin(wt),其中,X 是振幅, w 是单位为弧度/秒的频率;–瞬间动态分析:方程保持上述的形式。