第七章 结构振动的有限元分析
探索土木工程中的结构分析方法
探索土木工程中的结构分析方法提要:土木工程是一门关于设计、建造和维护基础设施的学科,其
中结构分析是其中重要的一部分。本文将探索土木工程中常用的结构
分析方法,包括有限元分析、结构振动分析和结构可靠性分析。这些
方法能够帮助工程师评估和改进建筑物和桥梁等结构的性能和可靠性。
一、有限元分析
有限元分析是土木工程中最常用的结构分析方法之一。它通过把结
构划分为有限数量的单元,然后利用数学模型来计算每个单元的应力
和应变。这些单元的力学行为可以用微分方程和矩阵运算表示。有限
元分析可以帮助工程师了解结构在不同载荷下的行为及其强度。此外,有限元分析还可以进行结构优化,并提供改进设计的指导。
二、结构振动分析
结构振动分析是研究结构在外部激励下的振动响应的方法。土木工
程中的结构通常会受到地震、风力和交通等外部激励的影响,因此对
结构的振动特性进行分析十分重要。结构振动分析可以帮助工程师判
断结构的自然频率、共振情况以及其对外界激励的响应。这些信息对
于设计抗震性能良好的建筑和桥梁至关重要。
三、结构可靠性分析
结构可靠性分析是评估结构在服役期内的可靠性和安全性的方法。
土木工程中的结构一般会处于不断变化的环境条件下,如荷载、温度
和湿度等。通过结构可靠性分析,工程师可以计算出结构在特定使用
寿命内的可靠性水平,并评估所需的维护和修复工作。这有助于确保
结构在使用过程中不会发生失效或损坏,保证人们的生命和财产安全。
结论:在土木工程中,结构分析方法的应用至关重要。有限元分析
能够帮助工程师理解和优化结构的力学行为。结构振动分析则可以帮
航空航天工程中的结构振动与噪声分析与优化
航空航天工程中的结构振动与噪声分析与优
化
航空航天工程是一个高度复杂和要求严苛的领域,结构振动与噪声控制是其中非常重要的一个方面。航空航天器的结构在运行过程中会受到各种力的作用,这些力引起的振动会产生噪声,进而影响到乘客的舒适感和航空设备的性能。因此,在航空航天工程中,进行结构振动与噪声分析与优化显得尤为重要。
一、结构振动分析
1. 自由振动分析
自由振动是指结构在没有外部强迫力作用下的振动。自由振动分析能够帮助工程师了解结构的天然频率、振型和阻尼等特性。通过获取这些信息,可以为后续的设计和优化提供基础。
2. 强迫振动分析
强迫振动是指结构在受到外界激励力的作用下产生的振动。强迫振动分析能够帮助工程师了解结构的响应情况,判断是否会产生共振,从而避免结构的失效或破坏。
二、噪声分析
1. 噪声源的识别
通过对飞机各个部位进行噪声分析,可以确定哪些部位是主要的噪声源。例如,发动机、起落架、风扇等都可能产生噪声,并且不同的噪声源会受到不同的因素影响,如转速、气流速度等。
2. 噪声传播路径分析
噪声的传播路径是指噪声从噪声源传播到乘客耳朵的路径。在噪声传播路径分析中,需要考虑噪声的传播特性,例如传播的距离、传播的方式(空气传播或结构传播)等。
三、优化措施
1. 结构优化
针对结构振动与噪声问题,可以通过优化结构设计来减小振动和噪声的影响。例如,采用合适的材料、调整结构的几何形状、增加结构的刚度等方法可以减小结构的振幅和共振现象,从而减少噪声产生。
2. 声学优化
声学优化是指通过合理的声学设计和控制方法来减少噪声的产生和传播。例如,在发动机设计中,可以通过优化燃烧室结构、改进排气系统等方式来减少发动机噪声。
机械振动第7章-弹性体振动
C3 C1
Y ' (0) 0 Y ' (l) 0
C2 C4 0 C4 C2
(chl cosl)C1 (shl sin l)C2 0
(shl sin l)C1 (chl cosl)C2 0
特征方程
coslchl 1
1l 4.730
2l 7.853
由于弦线上传播x 的横波在 x 方向无振动,
所以作用在微分段 ds 上的张力的 x分量应
该为零,即 T2 cos2 T1 cos 1 0
又根据牛顿第二定律,在y 方向微分段的运 动方程为:
T2 sin 2
T1 sin
1
ds
2 y t 2
对于小的振动,可取
ds dx
而 1, 2 都很小,所以
2u E 2u
t 2 z 2
2u t 2
VP2
2u z 2
(1)
u u(C1cosωnt C2sinωnt)
VP-桩身材料的[纵波]波速(m/s)
C1、C2由桩顶桩端的边界条件(见下)确定:
摩擦桩:
u
|z0
0,
u z
| z l
0
端承桩:
u z
|z0
0,
u z
| z l
0
2.考虑桩周土的影响 同上述,但结果很复杂。
第七章结构振动的有限元分析
第七章结构振动的有限元分析
第一节引言
结构振动是指结构在外力的作用下发生的同步振动。它在工程结构的
设计和分析中具有重要的意义。传统的结构振动分析方法主要有模态分析
法和频域分析法。近年来,随着计算机技术的发展,有限元方法在结构振
动分析中的应用越来越广泛。
第二节有限元方法概述
有限元方法是一种通过将连续结构离散化为有限个单元,并在每个单
元内进行局部计算,然后将单元组装起来进行全局分析的一种方法。有限
元方法的基本思想是将连续体分解为有限个离散单元,通过求解每个单元
的位移和应变,进而得到整个结构的力学行为和响应。在结构振动分析中,有限元方法可以更准确地描述结构的边界条件和模态特性。
第三节有限元建模
有限元建模是有限元分析的关键步骤之一、在有限元建模过程中,需
要根据实际情况选择适当的单元类型和单元尺寸,并确定边界条件。有限
元建模的准确性直接影响到振动分析的有效性和准确性。
第四节模态分析
模态分析是结构振动分析的常用方法之一、它可以通过求解结构的本
征频率和本征振型,对结构进行全面的振动特性分析。在有限元分析中,
模态分析主要通过求解结构的特征值问题来实现。
第五节动力分析
动力分析是结构振动分析的另一种常用方法。与模态分析相比,动力
分析能够更真实地反映结构在外力作用下的振动响应。在有限元分析中,
动力分析主要通过求解结构的动力方程来实现。
第六节振动问题的求解技巧与注意事项
在进行结构振动的有限元分析时,需要注意一些技巧和问题。首先,
应正确选择结构的边界条件和单元类型,以保证分析结果的准确性。其次,应注意振动问题的约束条件和模态解耦技巧,以简化计算过程。此外,在
结构有限元分析
第七章 结构有限元分析
引 言
求解具体结构工程中的问题是有限元素法的最终目的,而实际工程结构是复杂多样的,要很好的运用有限元素法还得解决好像坐标变换、对称边界条件运用以及复杂结构连接等问题。本章即为解决有限元方法应用于工程结构中实际问题的算法。
一、杆系或梁系的刚度坐标变换
1、向量的坐标变换公式
i) 一维向量的平面分解
θθsin cos v u u +=
⎭
⎬⎫
⎩⎨⎧=v u ]sin [cos θθ
ii) 一维向量的三维空间分解
⎪⎭
⎪
⎬⎫⎪⎩⎪⎨⎧=w v u u ]cos cos [cos γβα
iii) 平面向量的坐标变换:⎭
⎩⎥⎦⎤⎢
⎣⎡-=⎭⎬⎫⎩⎨⎧v v u θθ
θθcos sin sin cos
⎪⎭
⎪
⎬⎫⎪⎩⎪⎨⎧⎥⎦⎤⎢⎣⎡=⎭⎬⎫⎩⎨⎧w v u n m l n m l v u 22
2
111
2、杆元局部系下刚阵与整体系下刚阵的变换 i )局部系下的单元平衡方程:
⎪⎭
⎪
⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎦⎤⎢⎣⎡--j i j i p p u u l EA 1111 []{}{}p K =δ
由坐标变换(对节点力)
x
x
⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎥⎦
⎤
⎢⎢⎢⎢⎣⎡=
⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧j i jy jx iy ix p p p p p p θθθθsin 0cos 0
0sin 0cos =>{}{}
P P T T ⎥⎦
⎤⎢⎣⎡=λλ00 由局部坐标系下的平衡方程
{}[]{}
δλλK P T T
⎥⎦
⎤
⎢
⎣⎡=00 由位移(节点)的坐标变换
[]⎭
⎬⎫
⎩⎨⎧=i i i v u u θθs i n c o s
机械结构的振动模态分析与优化设计
机械结构的振动模态分析与优化设计引言
机械结构的振动问题一直是工程领域的研究热点之一。振动问题主要影响结构的安全性、可靠性和性能。因此,在机械结构设计过程中,振动模态分析和优化设计显得尤为重要。本文将探讨机械结构振动模态分析和优化设计的方法与实践。
一、振动模态分析
振动模态分析是研究结构振动特性的一种方法。通过振动模态分析,可以获取结构的固有频率、振型形态和模态的阻尼特性等信息。振动模态分析的目的是为了了解结构的振动特性,为优化设计提供依据。
在进行振动模态分析时,首先需要使用有限元分析(FEA)的方法建立结构的有限元模型。然后,通过求解结构的特征值问题,可以得到结构的固有频率和模态形态。振动模态分析的结果可以通过模态分析软件进行可视化展示,更加直观地观察结构的振动行为。
二、振动模态的影响因素
振动模态的特性受到多个因素的影响。首先,结构的几何形状和材料性质是影响振动模态的主要因素。例如,结构的尺寸和形状会影响固有频率和振型形态。材料的弹性模量和密度也会影响结构的固有频率。
其次,结构的边界条件和约束条件也会影响振动模态的特性。边界条件是指结构与周围环境的约束关系,例如结构的支座条件。约束条件是指结构内部各部件之间的连接关系。边界条件和约束条件会影响结构的振动自由度,并且改变结构的固有频率和振型。
最后,结构的质量分布和强度分布也会对振动模态产生影响。质量分布是指结构各部件的质量分布情况,不同的质量分布会导致结构的固有频率不同。强度分布
是指结构各部件的强度和刚度分布情况,不同的强度分布会导致结构的振动行为不同。
《弹性力学与有限元》第7章空间问题的有限元分析
11 1 1
11 1 1
Vi
=
1 6
x y
xj yj
xm ym
xp yp
,Vj
=
1 6
xi yi
x y
xm ym
xp , yp
z z j zm zp
zi z zm z p
1 111
11 11
Vm
=
1 6
xi yi
xj yj
x y
xp yp
,Vp
Li
=
Vi V
, Lj
=
Vj V
, Lm
=
Vm V
, Lp
=
Vp V
为体积坐标。其中 V 是四面体 ijmp 的体积, Vi、Vj、Vm、Vp 分别是四面体 ojmp、oimp、 oijp 和 oijm 的体积。
热动3同学修改
王正伟 13601363209
11 1 1 V = 1 xi x j xm xp ,
热动3同学修改
王正伟 13601363209
形状函数
而且满足
∑8
u = N1u1 + N2u2 +" + N8u8 = Niui
⎫ ⎪
i =1
⎪
∑8
v = N1v1 + N2v2 +" + N8v8 = Nivi
机械工程中的结构振动分析与优化设计
机械工程中的结构振动分析与优化设计
引言
机械工程是一门综合性的学科,多个领域交叉与融合,其中结构振动分析与优化设计是其重要的一部分。随着科技的发展和工业的进步,机械结构的性能要求越来越高,因此对结构振动的研究与优化设计显得尤为重要。本文将从结构振动的基本概念入手,介绍结构振动的分析方法和优化设计的原理,并探讨现代工程中的一些应用案例。
一、结构振动的基本概念
结构振动是指机械结构在受到外界激励作用后,由于其自身的刚度和阻尼特性导致的一种以周期性变化为特征的运动。振动现象既存在于自然界中,也广泛应用于各个领域的工程设计中。结构振动的研究主要包括振动的起因、振动的传播和振动的控制等方面。
二、结构振动的分析方法
1.模态分析
模态分析是结构振动研究的重要方法之一。通过求解结构的本征值问题,可以得到结构的固有频率和振型。固有频率代表了结构振动的固有特性,而振型则代表了结构在不同时刻的变形情况。模态分析主要通过有限元分析来实现,在计算软件的支持下,可以准确地得到结构的模态参数。
2.频域分析
频域分析是一种将时域信号转换为频域特性的方法。在结构振动分析中,频域分析可以用来研究振动信号的频率成分和谐波分析等问题。常用的频域分析方法包括傅里叶变换、功率谱密度分析等。
3.时域分析
时域分析是对结构振动信号的直接观测和分析。通过时域分析可以得到振动信号的时间响应和幅值特性。时域分析方法包括时序分析、自相关分析和互相关分析等。
三、结构振动的优化设计原理
结构振动的优化设计是通过对结构参数的调整和优化,最大限度地改善结构的振动性能。在优化设计中,需要确定优化目标函数和约束条件,通过数学模型和计算方法得到最优解。常用的结构优化设计方法包括数值优化方法、遗传算法、粒子群算法等。
动力学问题的有限元法(PDF)
第七章 动力学问题的有限元法
结构动力学是研究动载荷作用下结构动力反应规律的学科,讨论结构在动力荷载作用下反应的分析方法,寻找结构固有动力特性、动力荷载和结构反应三者间的相互关系。研究结构在动力荷载作用下的反应规律,能够为结构的动力可靠性(安全、舒适)设计提供依据。
前面介绍的静力学问题的研究对象是受不随时间变化的载荷作用。而动力学问题的对象受随时间而变的载荷的作用,从而使在结构中产生的位移、速度、应力和应变都随时间而变。当结构受随时间变化的载荷作用,且这种载荷的作用对结构的变形和应力的产生起主要作用,以致影响设备的安全性,或舒适性。这时就要进行动力学分析,充分认识其规律性,从设计阶段就抑制这种不利状况的发生。例如,有时虽然动载荷不大,但结构在交变力的作用下,其某些固有频率与激励力的作用频率相接近时,就会引起很大的振动、变形或应力,这时,就必须对结构作动力学分析。又如,要利用结构在周期性作用力驱动下的定向振动,例如利用这种运动输送产品,这时,就必须巧妙地设计结构,使其具有某些与激励频率一致的固有频率,并且使结构对激励具有适当的响应能力。总之,不管是利用振动,还是抑制振动,都需要进行结构动力学分析。
当前结构动力学的研究内容有三类。第一类问题:反应分析(结构动力计算),第二类问题:参数(或称系统)识别,第三类问题:荷载识别。第一类问题是已知系统动态特性和动载荷作用部位及大小,求出系统的响应——随时间变化的位移,速度,加速度和应力等。第二类问题是已知系统的输入输出特性,分析系统固有的动态特性,结构模态分析就属于这一类问题。第三类问题是在已知系统动态特性的条件下, 通过测量系统的响应,或由响应准则预先给出响应要求, 以此识别对响应的外载荷。三类结构动力学研究内容的载荷、结构和响应之间的关系如图7-1所示。动载荷种类大致分类如图7-2所示。
有限元分析第七章第一部分
第七章 四阶问题(板的弯曲)
在常见的工程结构中,板或板梁结构较为普遍。在有限元分析中板可分为薄板和中厚板,在对薄板分析时采用了克希霍夫(Kirchhoff )假设。板中面上任一点(x, y )允许有三个位移分量,其中面内位移u 、v 构成一平面应力问题(二阶问题)。横向位移则构成一个四阶问题(弯曲问题)。对于线性问题(小挠度),这两个问题之间没有耦合。可以分别进行研究,再将结果迭加。平面应力问题在本章以前进行了研究。本章着重讨论板的弯曲问题。在对中厚板进行分析时,重点将介绍位移和转角各自独立插值的板单元,这种板单元考虑了板的剪切变形。而且,弯曲问题可以降阶为二阶问题来描述。(对插值函数要求C O 连续)
§7-1薄板小挠度弯曲的基本方程
设板中面的横向位移(挠度)为w (x, y )
1、几何关系
中面法线绕x, y 轴的转角
曲率
其中 又称为扭率
2、弹性关系 弯矩
其中M xy 又称为扭矩。t 为板的厚度
图7-1
图7-2 ⎪⎪⎩
⎪⎪⎨
⎧∂∂-=∂∂=x w
y w y x θθ{}⎪⎪⎪
⎭⎪
⎪
⎪⎬⎫⎪⎪
⎪⎩⎪⎪⎪⎨⎧∂∂∂-∂∂-∂∂-=y x w y w x w 222222χ(7-1-1) (7-1-2) y
x w
∂∂∂-22{}[]{}χννννD y x w y w x w Et M M M M xy y x =⎪
⎪⎪⎭⎪⎪⎪⎬⎫
⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂∂-∂∂-∂∂-⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=222222
32)1(000101)1(12(7-1-3)
3、平衡方程
有限元分析在工程结构优化中的应用研究
有限元分析在工程结构优化中的应用研究
简介:
工程结构的优化设计是现代工程领域中的核心问题之一。在传统设计方法中,
需要进行大量的试验和实际测试来确定最佳设计方案,这样不仅费时费力,还无法保证最终设计结果的准确性。有限元分析作为一种数字仿真技术,为工程结构优化提供了一种高效、精确和可靠的方法。本文将探讨有限元分析在工程结构优化中的应用及其研究进展。
一、有限元分析的基本原理
有限元分析是基于数值计算方法的一种结构分析技术。其基本原理是将复杂的
结构分割成有限数量的单元,通过求解单元的力学行为来推导整个结构的力学行为。有限元分析基于强大的计算能力,可以模拟各种力的作用下,结构的变形和应力分布情况。通过有限元分析,工程师可以对结构进行优化设计、评估其性能,并预测其在各种工况下的响应。
二、有限元分析的应用领域
1. 结构强度分析
有限元分析可用于评估结构的强度和刚度,帮助设计师选择合适的材料和几何
形状,以满足设计要求。通过有限元分析,可以确定结构在各种载荷条件下的变形和应力分布,进而优化结构的设计。
2. 结构振动分析
振动问题在工程中具有重要意义,如航空航天工程中的飞行器结构振动、建筑
工程中的地震响应等。有限元分析可以帮助工程师预测结构的固有频率和振型,并
评估其对外界激励的响应。通过优化结构的材料和几何参数,可以降低振动干扰和振动噪声。
3. 热力学分析
有限元分析还可以用于分析结构的温度分布、热应力和热传导等问题。在工程中,如锅炉、发动机等高温设备、电子器件等都需要进行热力学分析。利用有限元分析,可以预测结构在不同温度下的变形和损伤情况,并优化设计,以保证其稳定运行。
第七章 结构振动的有限元分析
/POST1 !进入一般的后处理 /VIEW, 1 ,1,1,1 !设置视角 /ANG, 1 /REP,FAST SET,FIRST !调出第1阶的模态结果 SET,NEXT !调出下一阶的模态结果 SET,NEXT !调出下一阶的模态结果(实际上,这时为 第3阶模态) PLDI, , !显示变形后的结构图 ANMODE,10,0.5, ,0 !进行动画显示,设置10帧,每 帧显示0.5秒
静力学情形(static case) 静力学情形 由于与时间无关,则方程(7-19)退化为
Βιβλιοθήκη Baidu
无阻尼情形(undamped system) 无阻尼情形 此时 ,则方程(7-19)退化为
无阻尼自由振动情形(free vibration of undamped system) 无阻尼自由振动情形 则 , ,方程(7-19)退化为
单元内的位移插值函数为
其中N(ξ) 为单元的形状函数矩阵,与相对应的静力问题单 元的形状函数矩阵完全相同,ξ为单元中的几何位置坐标。
基于上面的几何方程和物理方程以及式(7-11),将相关的 物理量(应变和应力)表达为节点位移的关系,有
单元有限元方程
将单元的各个矩阵进行组装,可形成系统的整体有限元方 程,即
第7章 结构振动的有限元分析 章 7.1 结构振动分析的基本原理 结构的振动分析将涉及到 模态分析(modal analysis)、 瞬态动力学分析(transient dynamics analysis)、 简谐响应分析(harmonic response analysis)、 随机谱分析(spectrum analysis)等方面, 其中结构的模态分析(固有频率与振型)将是所有振动分析的 基础,下面将就结构的模态分析进行阐述。
结构动力学第七章
概述
• 振型叠加法是求解线弹性多自由度体系动力反应的行之有效方法,在 定前若干阶振型和自振频率之后,任何线性结构动力反应的近似解都 容易求得。 • 我们实际所面临的结构范围十分广泛,从只有几个自由度的高度简化 数学模型、只需要考虑一、二阶模态就能求得动力反应的近似解,一 到包含几百甚至数万个自由度的高度复杂的有限单元模型,其中可能 50或100阶模态对结构动力反应有重要影响。 • 要求解大型结构至要求阶数的振型和频率,完全利用行列式方程的解 是困难的。从数学观点来看,求解各类结构的振型和自振频率属于矩 特征值问题,自然地,可以利用矩阵特征值的求解技术来处理结构振 和自振频率的求解问题。
B
(1) ij
(2) ij
=B
(1) ji
(2) ji
= ∫ m( x)φi ( x)φ j ( x)dx
0
l
B
=B
= ∑ mkφi ( xk )φ j ( xk )
k =1
n
∑(A
j =1
m
ij
− ω 2 Bij )ai = 0
(i = 1, 2,L, m)
令上式的系数行列式等于零,即得频率方程为
ω2 =
∫ ∫
l 0
l
0
EI [U ′′( x)]2 dx
n 2 i =1
m( x)U ( x)dx + ∑ mU 2 ( xi ) i
有限元分析基础
i i
uj
j
T j
结点力向量为
F i eU i V i M ieT F jeU j V j M je T
单元e结点力列向量为
Fe F F ije e U i V i M i U j V j M jeT
31
(1)结点: ① 铰结点;② 刚结点;③ 混合结点。 (2)支座: ① 活动铰支座;② 固定铰支座 ;
③ 固定支座 ;④ 定向支座
15
第二章 结构几何构造分析
2.2.2 结构的分类与基本特征
(1) 按结构在空间的位置分 结构可分为平面结构和空间结构两大类
(2) 按结构元件的几何特征分 ① 杆系结构: 梁、拱、桁架、刚架、桁构结构等 。 ② 板壳结构 ③ 实体结构实体结构的长、宽、高三个尺寸都很 大,具有同一量级。 ④ 混合结构
对称结构在正对称载荷下,对称轴截面上只能产生 正对称的位移,反对称的位移为零;对称结构在反对称 载荷下,对称轴截面上只有反对称的位移,正对称的位 移为零。 (1) 具有奇数跨的刚架
① 正对称载荷作用
(a) 对称刚架
(b) 变形状态分析
(c) 对称性利用
图2-22对称性利用示意图
19
第二章 结构几何构造分析
单元之间的位移协调性。
由单元结点位移,确定待定系数项
当 x 0时, u ui
用于振动分析的有限元方法资料
求解特征值: 乘以l/2EI得: 令系数矩阵的行列式等于0得: 方程的根即梁的自然频率:
结果可以和精确解比较:
12.7——一致、集中质量矩阵
12.3节中推出的质量矩阵是一致质量矩阵,因为用于推导刚度 矩阵的位移模型也用于推导质量矩阵。
一些动态问题可以用形式简单的质量矩阵求解。最简单的质量 矩阵——集中质量矩阵,可以通过将质点指定到节点上。
▪ 前文中,节点没有固定,结构在节点力的作用下会发生刚体 位移。 ▪ 也就是说,矩阵[K]是奇异矩阵。通常情况下,我们希望结构 的位移为零。 ▪ 因此,我们需要添加边界条件对矩阵[M]、[K]和向量F进行约 束。
例1:杆件分析
如图:均质;长0.5m;断面截面积5e-4m^2;杨氏模量 200GPa;密度7850Kg/m^3;左端固定。
(19)
梁单元的动能、势能、虚功表达式分别为:
式中I是 横截面的
惯性矩
上式中:
(20) (21) (22)
通过上式,可以得到梁单元的质量、刚度矩阵,等效节点力:
二 单元矩阵的坐标变换
局部坐标系:以各个单元本身的轴线为基准所设立的坐标系。便于计算节点位移。 缺点:如果整个系统里各个单元取向各异,各个节点位移方向不一致,如下图。如 何使汇交于一个节点的各个杆件的节点位移真正相等? 解决方法:进行坐标变换
a.节点2处施加1000N静态轴向外力u2,求应力 b.求系统固有频率
有限元分析—模态分析
有限元分析—模态分析
有限元分析是一种结构力学领域的分析方法,可以对结构进行数值求解,以获得其固有频率和振型。模态分析是其中的一种应用,用于研究结
构在固有频率下的振动情况。本文将介绍有限元分析的基本原理、模态分
析的步骤和应用,并讨论其在实际工程中的重要性。
有限元分析是一种利用数值方法对结构进行力学分析的技术。它将结
构离散化为有限数量的单元,通过单元之间的相互作用来模拟整个结构的
力学行为。在进行模态分析时,通常采用线性弹性模型,即假设结构在固
有频率下是线性弹性振动的。
模态分析的主要目标是确定结构的固有频率和振型。固有频率是结构
自由振动的频率,与结构的几何形状、材料性质和边界条件等相关。振型
则描述了结构在不同频率下的振动模式。通过模态分析,可以了解结构在
特定频率下的振动情况,为结构设计和改进提供依据。
模态分析的步骤主要包括:建模、网格划分、边界条件的定义、求解
和结果分析。建模是指将实际结构抽象为数学模型,在计算机上进行仿真。网格划分是将结构划分为有限数量的单元,以便进行数值求解。边界条件
的定义是指确定结构的受力和支撑情况,包括约束、荷载等。求解是指通
过数值计算方法求解结构的固有频率和振型。结果分析是对求解结果进行
解释和评价,了解结构的振动特性。
模态分析在工程中具有广泛的应用。首先,它可以用于优化结构设计。通过模态分析,可以评估结构在不同固有频率下的振动情况,从而优化结
构的设计参数,使其在工作频率下保持稳定。其次,模态分析可以用于故
障诊断。结构的振动特性在受到损伤或故障时会发生变化,通过模态分析
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
LSTR,1,2 !连接点1和2生成直线1 LSTR,5,1 !连接点5和1生成直线2 BSPLIN,2,3,4,5,,,-0.025,0,0,-0.025,-0.00625,0, !采用B样条,连 接点2,3,4,5生成曲线3 AL,1,2,3 !由线1,2,3围成一个面 ESIZE,0.00625 !在单元划分前,定义单元的边的尺度为 0.00625 MSHAPE,0,2D !设置单元划分的类型为2D四边形(key=0) MSHKEY,0 !设置网格的自由划分(0) AMESH,all !对所有的面进行网格划分(无设置时,则默认为 采用第1号类型单元)
静力学情形(static case) 静力学情形 由于与时间无关,则方程(7-19)退化为
无阻尼情形(undamped system) 无阻尼情形 此时 ,则方程(7-19)退化为
无阻尼自由振动情形(free vibration of undamped system) 无阻尼自由振动情形 则 , ,方程(7-19)退化为
TYPE, 2 !设置单元类型2(空间单元) EXTOPT,ESIZE,10,0, !设置体单元扩展为10段 VEXT,all,,,0,0,0.25,,,, !对所有的面进行z方向的体(包括单 元)扩展,每次扩展的z方向增量为0.25 ESEL,U,TYPE,,1 !除单元类型1外,选择所有的单元(实际 上就是体单元) NSEL,S,LOC,Z,0 !选择z=0的节点 D,all, , , , , ,ALL, , , , , !对所选择的节点施加全部的固定约 束 NSEL,ALL !选择所有的节点 FINISH !前处理结束
/SOLU !进入求解模块 ANTYPE,2 !设置模态分析(2) MODOPT,LANB,5 !设定LANB方法,提取5阶模态 MXPAND,5, , ,0 !设定模态扩展数为5 MODOPT,LANB,5,0,0, ,OFF !设定LANB方法,计算5 阶模态 SOLVE !求解 FINISH !退出求解模块
杆单元的质量矩阵 质量矩阵分为两种,即一致质量矩阵和集中质量矩阵。 (1)一致质量矩阵 对于二节点杆单元,在局部坐标内有节点位移列阵和形状 函数矩阵
相应的质量矩阵为
所谓一致质量矩阵(consistent mass matrix)是指推导质量 矩阵时与推导刚度矩阵时所使用的形状函数矩阵相“一 致”。
(2)集中质量矩阵 将该二节点杆单元的质量直接对半平分,集中到二个节点 上,就可以得到集中质量矩阵(lumped mass matrix)为
可以看出,集中质量矩阵的系数都集中在矩阵的对角线 上,也就是说对应于各个自由度的质量系数相互独立, 相互之间无耦合;而一致质量矩阵的系数则有相互耦合。
平面三节点三角形单元的质量矩阵 (1)一致质量矩阵
7.1.3 常用单元的质量矩阵
结构振动分析将涉及到结构的刚度矩阵、质量矩阵和阻尼 矩阵,由(7-17)式可知,动力学问题中的刚度矩阵与静力 问题的刚度矩阵完全相同,而质量矩阵则通过(7-15)式来 进行计算,对于一种单元,只要得到它的形状函数矩阵, 就可以容易地计算出质量矩阵;由阻尼矩阵的计算公式 (7-16)可知,它的计算与质量矩阵相同,只是有关的系数 不同而已。下面给出常见单元的质量矩阵。
第7章 结构振动的有限元分析 章 7.1 结构振动分析的基本原理 结构的振动分析将涉及到 模态分析(modal analysis)、 瞬态动力学分析(transient dynamics Biblioteka Baidunalysis)、 简谐响应分析(harmonic response analysis)、 随机谱分析(spectrum analysis)等方面, 其中结构的模态分析(固有频率与振型)将是所有振动分析的 基础,下面将就结构的模态分析进行阐述。
/POST1 !进入一般的后处理 /VIEW, 1 ,1,1,1 !设置视角 /ANG, 1 /REP,FAST SET,FIRST !调出第1阶的模态结果 SET,NEXT !调出下一阶的模态结果 SET,NEXT !调出下一阶的模态结果(实际上,这时为 第3阶模态) PLDI, , !显示变形后的结构图 ANMODE,10,0.5, ,0 !进行动画显示,设置10帧,每 帧显示0.5秒
几何方程
物理方程
边界/初始条件 BC/IC 位移边界条件BC(u)
力边界条件BC(p)
初始条件IC(initial condition):
7.1.2 结构振动的有限元分析列式 用于动力学问题分析的单元构造与前面静力问题时相同,不 同之处是所有基于节点的基本力学变量也都是时间的函数 单元的节点位移列阵为
(2) 集中质量矩阵
机翼模型的振动模态分析 一个简化的飞机机翼模型如图所示,该机翼沿延翼方向为 等厚度。有关的几何尺寸见图,机翼材料的常数为:弹性 模量 ,泊松比 ,密度 ; 对该结构进行振动模态的分析。
/PREP7 ! 进入前处理 ET,1,PLANE42 !选取单元类型1(平面单元) ET,2,SOLID45 !选取单元类型2(空间单元) MP,EX,1,0.26e9 !定义材料的弹性模量(1号材料) MP,DENS,1,886 !定义材料的密度(1号材料) MP,PRXY,1,0.3 !定义材料的泊松比(1号材料) K, ,,,, !生成几何点1,坐标(0,0,0) K, ,0.05,,, !生成几何点2,坐标(0.05,0,0) K, ,0.0575,0.005,, !生成几何点3,坐标(0.0575,0.005,0) K, ,0.0475,0.0125,, !生成几何点4,坐标(0.0475,0.0125,0) K, ,0.025,0.00625,, !生成几何点5,坐标(0.025,0.00625,0)
7.1.1 结构振动分析的基本方程 描述结构动力学特征的基本力学变量和方程与前面的静 力问题类似,但增加了惯性力项和阻尼力项,且所有的 变量都将随时间而变化 结构振动的三大类变量 位移 应变 应力 是坐标位置ξ(x,y,z)和时间t的函数。
结构振动的三大类方程及边界/初始条件 结构振动的三大类方程及边界 初始条件 平衡方程(考虑惯性力和阻尼力)
其振动形式叫做自由振动(free vibration),该方程解的形式为 这是简谐振动的形式,其中ω为常数;将其代入(7-23)中, 有
该方程有非零解的条件是
这就是特征方程(eigen equation),ω为自然圆频率 (natural circular frequency)(rad/sec),也叫圆频率,对 应的频率为f=ω/2π(Hz)。求得自然圆频率ω后,再将其 代入方程(7-26)中,可求出对应的特征向量(eigen vector) ˆq ,这就是对应于振动频率ω的振型(mode)。
单元内的位移插值函数为
其中N(ξ) 为单元的形状函数矩阵,与相对应的静力问题单 元的形状函数矩阵完全相同,ξ为单元中的几何位置坐标。
基于上面的几何方程和物理方程以及式(7-11),将相关的 物理量(应变和应力)表达为节点位移的关系,有
单元有限元方程
将单元的各个矩阵进行组装,可形成系统的整体有限元方 程,即