第九讲结构动力学问题有限元法
合集下载
7.2 结构动力学的有限单元法
B：应变矩阵 ν：线性阻尼系数
三、结构的动力学方程
单元的刚度矩阵、质量矩阵和阻尼矩阵，要用来形成整体 的刚度矩阵、质量矩阵和阻尼矩阵。在不考虑体积力的条 件下，整个结构的动力学方程为：
& & Mδ& + Cδ + Kδ = f
当f = 0, c = 0时得到自由振动时的无阻尼动力方程 0时得到自由振动时的无阻尼动力方程
四、有限元法求解动力学问题实例
设第j个梁单元的单位长度质量为m 抗弯刚度为EI 设第j个梁单元的单位长度质量为mj, 抗弯刚度为EIj， 长度为L 两端的节点编号为j j+1，节点j 长度为Lj， 两端的节点编号为j和j+1，节点j的横向 位移为u 角位移为u 节点j+1的横向位移为u 位移为uj1, 角位移为uj2, 节点j+1的横向位移为uj3, 角 位移为u 位移为uj4。设梁单元的广义坐标向量为：
关于振型的正交性
δ K δm = δ ω M δm
2 m
^ T n
^
^ T n
^
2 δ K δ m = δ ωn M δ m
^ T n
^
^ T n
^
(ω − ω ) δ M δ m = 0
2 m 2 n
^ T n
^
ω ≠ω
2 m
^ T n ^
2 n
δ M δm = 0
^ T n
^
同样可以求得 δ K δ = 0 振型关于质量矩阵和刚度矩阵正交 m
ω T = [ω1 ω 2 ... ω N ]
四、方程的特征值及振型的正交性
对于齐次代数方程 ( K − ω M ) δ = 0 ，只能求得振型矢量 的各元素的相对值。
有限元 第9讲 动力学问题有限单元法
有限元第9讲动力学问题有限单元法动力学问题是指研究物体在运动中的受力和受力作用下的运动状态，常见的应用是结构工程学中的振动分析。
有限单元法是解决结构工程学中动力学问题的常用方法之一。
本文将介绍动力学问题和有限单元法的基本概念，并介绍其应用。
动力学问题的定义动力学是研究质点或刚体运动情况的分支学科，在结构工程学中是指结构在做振动时所受的力和运动状态。
动力学问题可以分为两种类型：稳态动力学问题和非稳态动力学问题。
稳态动力学问题是指结构在振动状态下所受的恒定力，而非稳态动力学问题则是指结构所受的变化的力，例如冲击力或地震力。
动力学问题的求解包括两个方面：一是确定受力情况；二是求解结构的运动状态。
确定受力情况通常需要通过实验或计算确定，求解结构运动状态则可以通过有限单元法来解决。
在结构工程学中，动力学问题的应用非常广泛。
例如，建筑物抗震设计需要对建筑物在地震作用下的反应进行分析，桥梁工程需要对桥梁在行车作用或风力作用下的振动响应进行分析。
有限单元法的基本概念有限单元法是一种将结构离散成若干小单元的数值分析方法，将结构分割成细小的单元，每个单元内部假设为均匀且连续的，通过对单元本身的运动状态进行求解，进而推知整个结构的运动状态。
有限元法用于解决的问题包括静力学问题、动力学问题、热力学问题和流体问题等。
有限单元法求解动力学问题的步骤主要包括如下几个步骤：1.离散化：将连续结构离散化成有限的小单元，每个单元内部运动状态通过定义一定数量的节点来确定。
2.建立单元的动力学方程：根据单元的形状和材料性质，建立单元的动力学方程，并计算单元的振动特性，例如频率和模态。
3.组装单元的方程：将单个单元的方程组装成整个结构的方程。
4.边界条件的处理：利用结构的边界条件（例如支撑、铰支等），将结构自由度减少到实际问题所需要的自由度。
5.求解结构的运动状态：通过求解整个结构的方程，得到结构的运动状态。
6.后处理：根据求解结果，进行结果的可视化和分析。
有限元法应用举例
核反应堆运行过程中涉及高温、 高压、高辐射等极端条件，热工 水力学分析是确保安全性的重要
环节。
有限元法可以对核反应堆的热工 水力学进行模拟，评估冷却剂流 动、热能传递、压力容器应力分
布等关键参数。
通过模拟分析，可以优化反应堆 设计，提高运行效率，降低事故
风险。
建筑物的能耗模拟与优化
建筑物的能耗是节能减排的重要领域，能耗模拟与优化有助于降低能源消耗和碳排 放。
况，为设备的电磁兼容性设计和优化提供依据。
通过有限元分析，可以评估设备的电磁辐射是否符合相关标准
03
和规定，以及优化设备的天线布局和结构设计等。
高压输电线路的电场分析
高压输电线路在运行过程中会 产生电场和磁场，其强度和分 布情况对环境和人类健康具有 一定影响。
有限元法可以用来分析高压输 电线路的电场分布情况，包括 电场强度的计算和分布规律的 分析等。
通过有限元分析，可以评估高 压输电线路对环境和人类健康 的影响，为线路的规划、设计 和优化提供依据。
07
有限元法应用举例：声学分析
消声室的声学设计
消声室是用于测试和测量声音的特殊 实验室，其内部环境需要极低的噪音 水平。
通过模拟和分析，可以确定最佳的吸 音材料和布局，以及最佳的隔音结构， 以达到最佳的消声效果。
有限元法应用举例
• 有限元法简介 • 有限元法应用领域 • 有限元法应用举例：结构分析 • 有限元法应用举例：流体动力学分析 • 有限元法应用举例：热传导分析 • 有限元法应用举例：电磁场分析 • 有限元法应用举例：声学分析
01
有限元法简介
定义与原理
定义
有限元法是一种数值分析方法，通过将复杂的物理系统离散 化为有限数量的简单单元（或称为元素），并建立数学模型 ，对每个单元进行单独分析，再综合所有单元的信息，得到 整个系统的行为。
结构动力学问题的有限元法
K Q
K Q
对于结构动力学问题，节点载荷阵还包括惯性力和阻尼力。
e e e K Q (M C ) e e 1 m
或改写为：
C K M Q
代入：
dV Q N u
T T T
M N N dV
dV N N
e T
e
e dV Q N u
e T T
N N dV C
其中：
M M C C
e
e
质量阵和阻尼阵的叠加方法与刚度阵的叠加方法相同，也 是对称稀疏阵。
三、动力方程的简化
M e N T N dV
称为一致质量矩阵，是稀疏带状阵。
如果将单元质量阵近似作为对角阵，则方程变成彼此独立，避免 联立，称为集中质量阵或团聚质量阵。 解耦 例如长度为L，截面积为A，密度为ρ的梁单元。 i
A，ρ
L
j
x
1 A L 0 集中质量阵： m 2 0 0
0 0 0 0
0 0 1 0
0 0 0 0
156 22L 22L 2 一致质量阵： 4 L AL m 13L 420 54 2 13 L 3 L
54 13L 13L 3L2 156 22L 2 22L 4 L
ˆ P K P K
T
在变换[K]和[M]的过程中，有时使用一次雅克比变换将一个 非对角线元素化为零以后，它在另一次变换中会重新变为非零 元素，但在素质上有所减小。这说明需要反复使用雅克比变换， 最终非对角线元素将趋于零。 在实际求解过程中，不必严格地把矩阵[K]和[M]所有的非对 角线元素变换为零，通常在完成一次变换后进行判断是否达到预 l 1 (l ) 设的精度：
动力学问题的有限元法
❖ 就结构的瞬态响应分析而言，典型的有结构在冲击载 荷下的响应问题。结构动力学中这类问题的特点是， 载荷作用前沿时间与构件的自振基频周期相近，远大 于应力波在构件中的传播时间。或者构件上长时间作 用随时间剧烈变化的载荷。
❖ 结构动力学问题在工程中具有普遍性。
3) 弹塑性动力学问题
❖ 这是连续介质变形体动力学问题的另一个重要领域。 涉及许多科学和工程领域，如高速碰撞，爆炸冲击， 人工地震勘探，无损探伤等。
力学问题。对等效系统应用虚功原理：
T
V
dV VuT(fu u)dV SuTTdS
• 将前面位移空间离散表达式和单元的几何方程、物理方 程代入上式虚功方程，并考虑到变分的任意性，得到离
散系统控制方程——结构有限元动力学方程：
Ma(t)Ca(t)Ka(t)Q(t)
方程中的系数矩阵分别为：系统质量矩阵，阻尼 矩阵，整体刚度矩阵。右端项为整体节点载荷向量。
u N ae
u(x, y, z,t) u来自v(x,
y,
z, t)
w( x, y, z, t)
a
e
a a
1 2
a n
ai
uvii
(t) (t)
(i
1,2,, n)
wi (t)
• 为建立有限元动力学响应控制方程，利用达朗倍尔原
理，在每个时刻 t，将连续介质中质点加上惯性力 u 和阻尼力 u ，则系统的动力学问题转化为等效静
• 如果忽略阻尼，则结构动力学方程简化为：
Ma(t) Ka(t) Q(t)
• 上式动力学方程的右端项为零时就得到结构自由振动 方程。
• 从动力学方程导出过程可以看出，动力学问题的有限元 分析中，由于平衡方程中出现了惯性力和阻尼力，从而 引入了质量矩阵和阻尼矩阵，运动方程是耦合的二阶常 微分方程组，而不是代数方程组。该方程又称为有限元 半离散方程，因为对空间是有限元离散的，对时间是连 续的。
❖ 结构动力学问题在工程中具有普遍性。
3) 弹塑性动力学问题
❖ 这是连续介质变形体动力学问题的另一个重要领域。 涉及许多科学和工程领域，如高速碰撞，爆炸冲击， 人工地震勘探，无损探伤等。
力学问题。对等效系统应用虚功原理：
T
V
dV VuT(fu u)dV SuTTdS
• 将前面位移空间离散表达式和单元的几何方程、物理方 程代入上式虚功方程，并考虑到变分的任意性，得到离
散系统控制方程——结构有限元动力学方程：
Ma(t)Ca(t)Ka(t)Q(t)
方程中的系数矩阵分别为：系统质量矩阵，阻尼 矩阵，整体刚度矩阵。右端项为整体节点载荷向量。
u N ae
u(x, y, z,t) u来自v(x,
y,
z, t)
w( x, y, z, t)
a
e
a a
1 2
a n
ai
uvii
(t) (t)
(i
1,2,, n)
wi (t)
• 为建立有限元动力学响应控制方程，利用达朗倍尔原
理，在每个时刻 t，将连续介质中质点加上惯性力 u 和阻尼力 u ，则系统的动力学问题转化为等效静
• 如果忽略阻尼，则结构动力学方程简化为：
Ma(t) Ka(t) Q(t)
• 上式动力学方程的右端项为零时就得到结构自由振动 方程。
• 从动力学方程导出过程可以看出，动力学问题的有限元 分析中，由于平衡方程中出现了惯性力和阻尼力，从而 引入了质量矩阵和阻尼矩阵，运动方程是耦合的二阶常 微分方程组，而不是代数方程组。该方程又称为有限元 半离散方程，因为对空间是有限元离散的，对时间是连 续的。
9 杆系结构的有限元法 -1
苏轼诗词的风格特点简析_苏轼的诗词代表作有哪些苏轼是北宋文坛革新的杰出领袖，也是宋代乃至中国历史上一个全能的文学天才。
每个作者都有自己的风格特点，那么苏轼诗词的风格特点是怎样的?下面就是小编给大家整理的苏轼诗词的风格特点，希望对你有用!苏轼诗词的风格特点一、气势恢弘、激情磅礴的豪放风格苏轼开创了豪放词风，是豪放词的代表，他能够借助瑰丽恢宏的意象来抒发慷慨豪情，将充沛激昂或悲壮苍凉的激情融入词中，并开拓了词的写作范围。
(一)气势磅礴、场景宏阔雄壮气势恢弘。
苏词的豪放词有气势豪迈飞动、场景宏阔雄壮的特点。
如：《念奴娇·赤壁怀古》大江东去，浪淘尽，千古风流人物。
故垒西边，人道是、三国周郎赤壁。
乱石穿空，惊涛拍岸，卷起千堆雪。
江山如画，一时多少豪杰。
遥想公瑾当年，小乔初嫁了，雄姿英发。
羽扇纶巾，谈笑间、强虏灰飞烟灭。
故国神游，多情应笑我，早生华发。
人生如梦，一樽还酹江月。
这是一首著名的豪放词，磅礴的气势、壮美的场景、雄浑的境界，正是豪放词风的体现。
词的上片写景，描写赤壁的景色，前三句，仿佛是描述远景，长江水浩浩荡荡，滔天波浪如大浪淘沙，送走了一代风流人物;次三句，描写近景，目光投向赤壁古战场;上片最后两句，是上片景物描写的总结。
下片怀古，词人抓住周瑜年轻有为的主要特征，塑造了他雄姿英发的英雄形象，并由此联想到自己华发早生，功业无成，不免产生人生如梦之感。
总体来看，这首词从江水的东流，感受到时光的逝去，进而把江山与人物合写，使“江山如画”与“风流人物”都得到形象的表现，写得雄浑豪放，气象恢弘，堪称历代咏史怀古诗词之绝唱，亦开后世豪放一派之先河，给人以撼魂荡魄的艺术力量。
这类词还如《水调歌头黄州快哉亭赠张握诠》(落日绣帘卷)等。
(二)直抒胸臆，自由豪放豪情奔放。
苏轼的词注重抒情言志，直抒胸臆，自由豪放。
如：《江城子·密州出猎》：“老夫聊发少年狂，左牵黄，右擎苍。
锦帽貂裘，千骑卷平冈。
欲报倾城随太守，亲射虎，看孙郎。
有限元法概述
但真正的应用实际问题是到1960年以后，随着电子数 值计算机的广泛应用和发展，有限单元法的发展速度才显 著加快。现代有限元法第一个成功的尝试，是将刚架位移 法推广应用于弹性力学平面问题，这是Turner，Clough 等人在分析飞机结构时于1956年得到的成果。他们第一 次给出了用三角形单元求得平面应力问题的正确解答。
（2）MSC/NASTRAN。 MSC/NASTRAN是在原NAST RAN基础上进行大量改进后的系统软件，主要包括MS C.Patran并行框架式有限元前后处理及分析系统、 MS C.GS-Mesher快速有限元网格、 MSC.MARC非线性有 限元软件等。其中MSC.MARC具有较强的结构分析能
.
5.在产品制造或工程施工前预先发现潜在的问题； 6. 模拟各种试验方案，减少试验时间和经费； 7. 进行机械事故分析，查找事故原因。
轴承强度分析
.
汽车碰撞实验
.
刹车制动时地盘的应力分析
.
钢板精轧机热轧制分析
.
三维椭圆封头开孔补强
.
水轮机叶轮的受力分析模拟
.
人体股骨端受力分析
.
半导体芯片温度场的数值仿真
知量时称为混合法。 位移法易于实现计算自动化，所以，在有限单元法
中位移法应用范围最广。
.
2、有限元法的发展
有限单元法基本思想的提出，可以追溯到Courantl在1 943年的工作，他第一次尝试应用定义在三角形区域上的 分片连续函数和最小位能原理相结合，来求解St·Venant 扭转问题。相继一些应用数学家、物理学家和工程师由于 各种原因都涉足过有限单元的概念。
.
4、有限元的特点
(1) 概念清楚,容易理解。可以在不同的专业背景和水平 上建立起对该方法的理解。从使用的观点来讲,每个人的 理论基础不同,理解的深度也可以不同,既可以通过直观的 物理意义来学习,也可以从严格的力学概念和数学概念推 导。
（2）MSC/NASTRAN。 MSC/NASTRAN是在原NAST RAN基础上进行大量改进后的系统软件，主要包括MS C.Patran并行框架式有限元前后处理及分析系统、 MS C.GS-Mesher快速有限元网格、 MSC.MARC非线性有 限元软件等。其中MSC.MARC具有较强的结构分析能
.
5.在产品制造或工程施工前预先发现潜在的问题； 6. 模拟各种试验方案，减少试验时间和经费； 7. 进行机械事故分析，查找事故原因。
轴承强度分析
.
汽车碰撞实验
.
刹车制动时地盘的应力分析
.
钢板精轧机热轧制分析
.
三维椭圆封头开孔补强
.
水轮机叶轮的受力分析模拟
.
人体股骨端受力分析
.
半导体芯片温度场的数值仿真
知量时称为混合法。 位移法易于实现计算自动化，所以，在有限单元法
中位移法应用范围最广。
.
2、有限元法的发展
有限单元法基本思想的提出，可以追溯到Courantl在1 943年的工作，他第一次尝试应用定义在三角形区域上的 分片连续函数和最小位能原理相结合，来求解St·Venant 扭转问题。相继一些应用数学家、物理学家和工程师由于 各种原因都涉足过有限单元的概念。
.
4、有限元的特点
(1) 概念清楚,容易理解。可以在不同的专业背景和水平 上建立起对该方法的理解。从使用的观点来讲,每个人的 理论基础不同,理解的深度也可以不同,既可以通过直观的 物理意义来学习,也可以从严格的力学概念和数学概念推 导。
有限元法PPT课件
和时间。
如何克服局限性
改进模型
通过更精确地描述实际 结构，减少模型简化带
来的误差。
优化网格生成
采用先进的网格生成技 术，提高网格质量，降
低计算误差。
采用高效算法
采用并行计算、稀疏矩 阵技术等高效算法，提
高计算效率。
误差分析和验证
对有限元法的结果进行误 差分析和验证，确保结果
的准确性和可靠性。
05 有限元法的应用实例
有限元法ppt课件
目 录
• 引言 • 有限元法的基本原理 • 有限元法的实现过程 • 有限元法的优势与局限性 • 有限元法的应用实例 • 有限元法的前沿技术与发展趋势 • 结论
01 引言
有限元法的定义
01
有限元法是一种数值分析方法， 通过将复杂的结构或系统离散化 为有限个简单元（或称为元素） 的组合，来模拟和分析其行为。
有限元法在流体动力学分析中能够处理复杂的流体流动和 压力分布。
详细描述
通过将流体域离散化为有限个小的单元，有限元法能够模 拟流体的流动、压力、速度等状态，广泛应用于航空、航 天、船舶等领域。
实例
分析飞机机翼在不同飞行状态下的气动性能，优化机翼设 计。
热传导分析
总结词
有限元法在热传导分析中能够处理复杂的热传递过程。
实例
分析复杂电磁设备的电磁干扰问题，优化设备性能。
06 有限元法的前沿技术与发 展趋势
多物理场耦合的有限元法
总结词
多物理场耦合的有限元法是当前有限元法的重要发展方向， 它能够模拟多个物理场之间的相互作用，为复杂工程问题提 供更精确的解决方案。
详细描述
多物理场耦合的有限元法涉及到流体力学、热力学、电磁学 等多个物理场的耦合，通过建立统一的数学模型，能够更准 确地模拟多物理场之间的相互作用。这种方法在航空航天、 能源、环境等领域具有广泛的应用前景。
如何克服局限性
改进模型
通过更精确地描述实际 结构，减少模型简化带
来的误差。
优化网格生成
采用先进的网格生成技 术，提高网格质量，降
低计算误差。
采用高效算法
采用并行计算、稀疏矩 阵技术等高效算法，提
高计算效率。
误差分析和验证
对有限元法的结果进行误 差分析和验证，确保结果
的准确性和可靠性。
05 有限元法的应用实例
有限元法ppt课件
目 录
• 引言 • 有限元法的基本原理 • 有限元法的实现过程 • 有限元法的优势与局限性 • 有限元法的应用实例 • 有限元法的前沿技术与发展趋势 • 结论
01 引言
有限元法的定义
01
有限元法是一种数值分析方法， 通过将复杂的结构或系统离散化 为有限个简单元（或称为元素） 的组合，来模拟和分析其行为。
有限元法在流体动力学分析中能够处理复杂的流体流动和 压力分布。
详细描述
通过将流体域离散化为有限个小的单元，有限元法能够模 拟流体的流动、压力、速度等状态，广泛应用于航空、航 天、船舶等领域。
实例
分析飞机机翼在不同飞行状态下的气动性能，优化机翼设 计。
热传导分析
总结词
有限元法在热传导分析中能够处理复杂的热传递过程。
实例
分析复杂电磁设备的电磁干扰问题，优化设备性能。
06 有限元法的前沿技术与发 展趋势
多物理场耦合的有限元法
总结词
多物理场耦合的有限元法是当前有限元法的重要发展方向， 它能够模拟多个物理场之间的相互作用，为复杂工程问题提 供更精确的解决方案。
详细描述
多物理场耦合的有限元法涉及到流体力学、热力学、电磁学 等多个物理场的耦合，通过建立统一的数学模型，能够更准 确地模拟多物理场之间的相互作用。这种方法在航空航天、 能源、环境等领域具有广泛的应用前景。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3、利用节点位移插值表示单元内任一点的位移
d=Nqe
一般仍采用与静力分析相同的形函数，[N]。当单元数量较多时，上述 插值可以得到较好的插值精度。 4、在线弹性条件下，单元内的应变和应力与节点位移的关系仍为
= B qe =DBqe
但这时的位移、应变和应力都是某一时刻的瞬时值，它们都是随时间t 变化的函数。
●三角形平面问题单元
2 0 1 0 1 0
2 0 1 0 1
惯性力是分布力，按分布力向节点等效的原则和 实施过程，有： 第九讲结构动力学问题有限元法
1.一致质量矩阵
R
e q
N
T
qdV
V
V
N
T
Nห้องสมุดไป่ตู้
2 t 2
e
dV
N T N dV e
V
于是，令
meN TNdV
第九讲结V构动力学问题有限元法
1.一致质量矩阵
me 的计算式是通式，并因为计算质量矩阵和刚度矩
第九讲结构动力学问题有限元法
第一节 动态分析有限元法的特点
一、载荷特点 结构所受的载荷是随时间变化的动载荷。 这是与静力分析的一个根
本区别。
二、位移特点 1、节点位移{q}不仅是坐标的函数，而且也是时间的函数。仍以节
点位移{q}作为基本未知量。
2、节点具有速度 q和q加速度。
第九讲结构动力学问题有限元法
第九讲结构动力学问题有限元法
式中 keBTDBdV
V
me NTNdV
V
ce NTNdV
V
分别称为单元的刚度矩阵、质量矩阵和阻尼矩阵，它们就是决定单元动态性能的 特性矩阵。
R t e N T P v d V N T P s d A N T P c
V
A
称为单元节点动载荷列阵，它是作用在单元上的体力、面力和集中力向单元节 点移置的结果。
第九讲结构动力学问题有限元法
由于 dN qe,B qe
且形函数仅为坐标x、y、z的函数，与时间无关，因此有
d N qe
d N qe
d N qe
B q e
根据虚位移原理，有
UW
代入经整理，可得单元运动方程为
m e q e c e q e k e q e R t e
动载荷（又称动力分析）
固有特性分析
响应分析
固 有 频 率
振 型
位 移 响 应
速 度 响 应
加 速 度 响 应
动 应 变
动 应 力
第九讲结构动力学问题有限元法
固有特性：是一组模态参数构成，它由结构本身（质量与刚度 分布）决定，而与外部载荷无关，但决定了结构对动载荷的响应； 响应分析：是计算结构对给定动载荷的各种响应特性。
在动态分析和静力分析中，单元的刚度矩阵是相同的，外部载荷的移置原理也一样。
第九讲结构动力学问题有限元法
单元质量矩阵根据其形成过程分为一致质量阵和 集中质量阵，各有自身的优点和缺点。
1.一致质量矩阵
在离散后的结构中，取出一个单元，根据达朗贝 尔原理，单位体积上作用的惯性力为：
q t2 2 t2 2N e N t2 2e
9 动态分析有限元法
第九讲结构动力学问题有限元法
工程中受动载荷的产品：受道路载荷的汽车；受风载的雷达； 受海浪冲击的海洋平台；受偏心离心力作用的旋转机械等。 动态分析的必要性：当产品受到随时间变化的动载荷时，需 要进行动态分析，以了解产品动态特性。
第九讲结构动力学问题有限元法
第九讲结构动力学问题有限元法
UTdV
V
单元除受动载荷外，还有加速度和速度引起的惯性力 ddV 和阻尼力 ddV ，其中ρ为材料密度，v是线性阻尼系数。外力所做的虚功为：
WdTPvdVdTPsdAdTPc
V
A
dTddVdTddV
V
V
式中，{Pv}、{Ps}、{Pc}分别为作用于单元上的动态体力、动态面力和动态 集中力；V为单元面积；A为单元面积。
第九讲结构动力学问题有限元法
二、单元分析
单元分析的任务仍是建立单元特性矩阵，形成单元特性方程。 动态分析中，单元特性矩阵：刚度矩阵、质量矩阵和阻尼矩阵。
动态分析中，仍采用虚位移原理建立单元特性矩阵。
第九讲结构动力学问题有限元法
在动载荷作用下，对于任一瞬时，设单元节点发生虚位移 qe ，则单元 内也产生相应的虚位移 d 和虚应变 。单元内产生的虚应变能为:
阵使用的形状函数一致，因此被称为一致质量阵。
2.集中质量矩阵
在工程实际中，为了求解方便，有人把单元质量 平均分到单元的各个节点上，如平面三角形单元的 质量可分配为：
mi mj mk 3 V dV
第九讲结构动力学问题有限元法
2.集中质量矩阵 单元质量矩阵为：
m e dm i im a i m jg m j m km k
第九讲结构动力学问题有限元法
5、由于节点具有速度和加速度，结构将受到阻尼和惯性力的作用。 根据达朗伯原理，引入惯性力和阻尼力之后结构仍处于平衡状态，因 此动态分析中仍可采用虚位移原理来建立单元特性方程，然后再集成。 整个结构的平衡方程为
M q C q K q R t
式又称运动方程，它不再是静力问题那样的线性方程，而是一个二阶 常微分方程组。
dx
4
Al
30
1
8
1 8 4 8 8 16 第九讲结构动力学问题有限元法
3.常用单元的一致质量矩阵
●三次梁单元
156 22l 54 13l
me Al 22l 4l2 13l
3l2
420 54 13l 156 22l
13l 3l2
22l
4l2
第九讲结构动力学问题有限元法
3.常用单元的一致质量矩阵
求解过程复杂，建立有限元模型时要特别注意控制模型规模。
第九讲结构动力学问题有限元法
第二节 动态分析有限元法的一般步骤
第九讲结构动力学问题有限元法
一、结构离散
该步骤与静力分析完全相同，只是应该分析内容不同，对网格形式的要求 有可能不一样。 静力分析：要求在应力集中部位加密网格； 动态分析：由于固有频率和振型主要与结构的质量和刚度分布有关，要求 整个结构采用尽可能均匀的网格形式。
3.常用单元的一致质量矩阵
●一次杆单元 N 1
me AlNTNdxAl121 2dx
A l
2 1
122dx
1第九2讲结构动力2学问题有限元法
l
6
2 1
1 2
3.常用单元的一致质量矩阵
●二次杆单元
me
A
NT
l
Ndx
A
l
222122
1 2
412
222122
1 2
412
T