CN101887474A - 基于有限元法与广义傅里叶级数法的结构振动分析方法 - Google Patents
基于有限元法与广义傅里叶级数法的结构振动分析方法 Download PDFInfo
- Publication number
- CN101887474A CN101887474A CN 201010208748 CN201010208748A CN101887474A CN 101887474 A CN101887474 A CN 101887474A CN 201010208748 CN201010208748 CN 201010208748 CN 201010208748 A CN201010208748 A CN 201010208748A CN 101887474 A CN101887474 A CN 101887474A
- Authority
- CN
- China
- Prior art keywords
- finite element
- stiffness matrix
- fourier series
- displacement
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明的目的在于提供基于有限元法与广义傅里叶级数法的结构振动分析方法。分为以下步骤:将需要进行振动分析的结构区域划分,分别形成相应的有限元表述区域与广义傅里叶级数表述区域;对有限元表述区域进行有限元网格划分并形成相应的质量刚度矩阵,依据广义傅里叶级数表述区域特点选择相应的假设位移形式,形成质量刚度矩阵;之后在两个区域之间建立虚拟弹簧利用能量变分方法,将虚拟弹簧势能转化为总体耦合刚度矩阵;然后对形成的质量刚度矩阵根据位移进行排列,形成总体结构质量刚度矩阵;求解线性方程组得到相应的节点位移和级数展开中的未知系数。本发明对大型较复杂结构,既可以获得比有限元法更高的精度,又可以节省大量的计算成本。
Description
技术领域
本发明涉及的是一种应用于工程力学和振动工程领域的结构振动分析方法。
背景技术
作为一种数值计算方法,有限元法在结构振动领域有着广泛的应用。在处理不规则复杂结构,该方法有着解析方法无法比拟的优势。但有限元法仍然存在一些缺点,如在计算大型复杂结构特别是高频问题甚至中频问题,由于现在计算技术限制,计算所需网格过多,导致求解非常的困难甚至无法求解。
在结构振动分析中,解析法具有结果精确,使用计算机求解具有计算结果的频带宽、精确可信以及占有计算机资源小,计算速度快等优点。但是,解析法往往只适用于简单规则结构。在解析方法中,傅里叶级数解法是近年来受到重视的方法之一,它可以适用于各种边界条件、运算方便,并由于级数间具有正交性,可以使计算量大为减小并保证很高精度。美国韦恩州立大学李文龙提出一种广义傅里叶级数方法成功解决了任意边界条件下横梁的弹性振动(W.L.Li,Vibrationanalysis of rectangular plates with general elastic boundary supports,Journal of Sound andVibration 273(2004)619-635.),求解出的固有频率和振型的精度及级数收敛速度都达到了十分理想的效果,此解法还被应用于求解和分析多跨度桥梁接受运动负载时的振动问题(W.L.Li,M.Daniels,A Fourier series method for the vibrations of elasticallyrestrained plates arbitrarily loaded with springs and masses,Journal of Sound and Vibration 252(2002)768-781.)。文献Vibrations of rectangular plates with arbitrary non-uniform elastic edgerestraints(X.Zhang,Wen L.Li*,Journal of Sound and Vibration 326(2009)221-234)中,国内杜敬涛等人将傅里叶级数方法解决了任意边界弹性边界矩形板振动问题,不均匀边界问题乃至是板与板耦合振动分析问题。对于梁、板、圆柱壳等规则结构,这些结构的微分方程是四阶的,展开级数具有四阶(或更高阶)的逐项可导的性质,这些结构均可用广义傅里叶级数方法进行结构振动求解。对于不规则结构或者规则结构的布尔运算,如矩形板上有三角形孔,广义傅里叶级数方法还无法求解,原因是一般结构的高阶逐项可导的条件很难满足,所以,傅里叶级数解法解决问题的范围受到很大限制。
将有限元法解决复杂结构能力强的优点和傅里叶级数展开法计算精度高、计算资源消耗小且计算速度快的优点有机结合,可以解决大型较复杂结构振动问题。目前还没有将两种方法结合的技术出现以及相关文献报道。
发明内容
本发明的目的在于提供可用来解决大型较复杂结构振动特别是中高频振动计算难题的基于有限元法与广义傅里叶级数法的结构振动分析方法。
本发明的目的是这样实现的:
本发明基于有限元法与广义傅里叶级数法的结构振动分析方法,其特征是:
(1)将待分析的结构划分为规则结构与非规则结构两部分,两部分之间利用虚拟弹簧连接;
(2)规则结构部分设定位移函数,即设定傅里叶级数及其附加容许函数,依据能量原理将傅里叶级数及其容许函数中的系数作为未知量,规则结构部分转化为当量刚度矩阵与质量矩阵,非规则结构部分采用有限元法构造该区域的总体刚度矩阵与质量矩阵;
(3)规则结构部分与非规则结构部分之间的虚拟弹簧两端位移分别以有限元节点位移以及级数展开的系数表示,利用弹簧两端位移表示弹簧储存势能,利用能量方法中的变分方法得到有限元离散区域与级数展开区域之间的耦合质量矩阵与刚度矩阵;
(4)有限元质量刚度矩阵、广义傅里叶级数展开质量刚度矩阵和虚拟弹簧的耦合质量刚度矩阵依据位移排列组合得到结构总体质量刚度矩阵;
(5)由结构总体质量刚度矩阵得到线性方程组,求解线性方程组得到相应的节点位移和级数展开中的未知系数。
本发明的优势在于:与有限元法相比,傅里叶级数展开方法作为解析方法虽然只能用于一些规则结构但该方法无需任何网格且有着收敛快速,计算所需资源少的优点;而有限元方法作为一种成熟技术,在结构领域已经得到非常广泛的应用,适用于任意形状的结构。由于结构的大型化复杂化,加上计算机的计算能力限制,导致有限元法计算有一定误差以及计算效率低下,对于一些复杂结构的中高频甚至无法求解。通过本发明提出的将两者结合的连接方法,对大型较复杂结构,既可以获得比有限元法更高的精度,又可以节省大量的计算成本。
附图说明
图1为本发明的流程示意图;
图2是本发明实施方式1的结构示意图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1,本发明可以分为以下步骤:
1、建立求解域,并将之分为规则结构如梁、矩形板等与非规则结构两部分,二者之间利用虚拟无穷大弹簧连接;
2、非规则部分采用有限元法。将非规则结构离散化成有限元单元,即将该区域分解成节点和单元;
3、假设代表单元物理行为的形函数,即假设代表单元解的近似连续函数,并对单元建立方程;
4、将单元组合成总体的问题,构造该区域的总体刚度矩阵与质量矩阵;
5、规则区域根据其控制方程特点,设定位移函数,即设定傅里叶级数及其附加容许函数;
6、依据能量原理将傅里叶级数及其容许函数中的系数作为未知量,规则区域转化为当量刚度矩阵与质量矩阵;
7、应用边界条件到不同区域的质量与刚度矩阵之中;
8、规则区域与非规则区域之间的虚拟弹簧两端位移分别以有限元节点位移以及级数展开的系数表示之;
9、利用弹簧两端位移,表示弹簧储存势能,同样的利用能量方法中的变分方法对之进行处理,得到有限元离散区域与级数展开区域之间的耦合质量矩阵与刚度矩阵。
10、分别利用形成的有限元质量刚度矩阵、傅里叶级数展开质量刚度矩阵以及连接二者的虚拟弹簧的耦合矩阵形成总体结构的总体质量与刚度矩阵。
11、求解得到的线性方程组,得到相应的节点位移和级数展开中的未知系数。
12、通过计算结果得到其他信息,如频率、各阶模态,亦可进一步做结构响应计算。
实施方式1:
结合图2,悬臂梁支撑左侧部分由有限元fa表示,支撑右侧部分由傅里叶级数展开方法表述。梁弯曲问题的基本方程可表示如下:
几何关系
应力应变关系
平衡方程
边界条件
或
或
用有限元法分析梁弯曲问题时,采用Hermite多项式作为单元的插值函数。对于两个节点的一维单元,函数φ采用Hermite多项式的插值表达式可写成
因此单元内挠度函数的插值表示如下
对泛函
取变分可以得到有限元的刚度矩阵
同理得到单元的质量矩阵
由此,有限元形成结构自由振动方程如式(11)所示:
解决弹性支承问题通常采用修改刚度矩阵的方法。以上述连续梁为例,可在相应节点上主刚度加相应值。
依据广义傅里叶级数方法及梁的特性假设位移条件形式如下
因此梁的曲率可表达如下:
梁的势能动能表达式如下
一个边界上有两个弹簧:分别为线性弹簧kw0,kw1,扭簧Kw0,Kw1。这样的话,边界上的能量表达式可以写为如下形式:
因此,***中储存的总势能为V=U+Ub0+Ub1,直梁拉格朗日函数可以表示为如下形式
L=V-T (17)
将式(12)~(16)代入式17,采用Rayleigh-Ritz法使拉格朗日函数对每个未知Fourier系数取极值,我们可以得到8个线性方程组,进一步写为矩阵表达式形式:
其中
通过以上推导,分别得到有限元与广义傅里叶级数展开方法的质量刚度矩阵。推导过程可以发现其质量刚度矩阵都由有变分导出,因此对于左右梁之间的弹簧kc,Kc作一个相同的处理,从而得到其耦合刚度与质量矩阵。
耦合项能量,
在耦合能量项中,存在左端有限元梁的最后一个节点的两个位移以及右端的Fourier系数,同样的对每个未知系数取极值,可以得到5个线性方程组。通过这五个线性方程组,将左侧的有限元梁与右侧的广义傅里叶级数梁耦合在一起使之成为一体。这五个线性方程组与式(11)、(18)合在一起整理之,可以得到这样的总体质量刚度矩阵
其中
通过式(20)即可求解相应结构的振动特性。
根据以上推导,对一根直梁几种不同边界(S-S,C-C,F-F)进行了验算,其参数如下所示:
长度为L=1m;
横截面积=0.0003;
惯性矩I=7.0e-009;
杨氏模量E=2.1e11。
计算所得结果如下所示,耦合方法结果与有限元,广义傅里叶级数方法结果相差很小。
两边简支边界
两边固支边界
两边自由边界
双边自由边界条件第五十阶达到精确值0.001%精度时的计算时间比较对比结果如下:
Claims (1)
1.基于有限元法与广义傅里叶级数法的结构振动分析方法,其特征是:
(1)将待分析的结构划分为规则结构与非规则结构两部分,两部分之间利用虚拟弹簧连接;
(2)规则结构部分设定位移函数,即设定傅里叶级数及其附加容许函数,依据能量原理将傅里叶级数及其容许函数中的系数作为未知量,规则结构部分转化为当量刚度矩阵与质量矩阵,非规则结构部分采用有限元法构造该区域的总体刚度矩阵与质量矩阵;
(3)规则结构部分与非规则结构部分之间的虚拟弹簧两端位移分别以有限元节点位移以及级数展开的系数表示,利用弹簧两端位移表示弹簧储存势能,利用能量方法中的变分方法得到有限元离散区域与级数展开区域之间的耦合质量矩阵与刚度矩阵;
(4)有限元质量刚度矩阵、广义傅里叶级数展开质量刚度矩阵和虚拟弹簧的耦合质量刚度矩阵依据位移排列组合得到结构总体质量刚度矩阵;
(5)由结构总体质量刚度矩阵得到线性方程组,求解线性方程组得到相应的节点位移和级数展开中的未知系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102087482A CN101887474B (zh) | 2010-06-25 | 2010-06-25 | 基于有限元法与广义傅里叶级数法的结构振动分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102087482A CN101887474B (zh) | 2010-06-25 | 2010-06-25 | 基于有限元法与广义傅里叶级数法的结构振动分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101887474A true CN101887474A (zh) | 2010-11-17 |
CN101887474B CN101887474B (zh) | 2012-02-01 |
Family
ID=43073392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102087482A Expired - Fee Related CN101887474B (zh) | 2010-06-25 | 2010-06-25 | 基于有限元法与广义傅里叶级数法的结构振动分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101887474B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104636556A (zh) * | 2015-02-09 | 2015-05-20 | 武汉理工大学 | 成任意角度连接的有限尺寸板结构振动响应计算方法 |
CN104918184A (zh) * | 2015-04-14 | 2015-09-16 | 哈尔滨工程大学 | 一种耦合声场的声学耦合方法 |
CN105184060A (zh) * | 2015-08-24 | 2015-12-23 | 哈尔滨工程大学 | 一种基于空间配面与改进傅里叶级数的层合结构振动分析方法 |
CN105426343A (zh) * | 2015-11-02 | 2016-03-23 | 株洲时代新材料科技股份有限公司 | 一种基于傅里叶级数的复杂结构振动解析分析方法 |
CN106294975A (zh) * | 2016-08-05 | 2017-01-04 | 大连理工大学 | 一种基于降阶模型的梁式结构自由振动分析方法 |
CN106407582A (zh) * | 2016-09-28 | 2017-02-15 | 中国工程物理研究院总体工程研究所 | 一种用于柔性机构静力学和动力学分析的高效计算方法 |
CN106844875A (zh) * | 2016-12-28 | 2017-06-13 | 湖南大学 | 一种基于傅里叶级数的高速凸轮优化设计方法 |
CN107808048A (zh) * | 2017-10-26 | 2018-03-16 | 哈尔滨工程大学 | 一种均匀截面梁结构的振动分析方法 |
CN109948180A (zh) * | 2019-01-25 | 2019-06-28 | 北京航空航天大学 | 一种正交各向异性对边简支矩形薄板振动分析方法 |
CN113158400A (zh) * | 2020-12-14 | 2021-07-23 | 东北电力大学 | 一种利用改进耦合能预测封闭箱体结构振动特性方法 |
-
2010
- 2010-06-25 CN CN2010102087482A patent/CN101887474B/zh not_active Expired - Fee Related
Non-Patent Citations (2)
Title |
---|
《哈尔滨工业大学学报》 20001031 周传月等 旋转对称结构振动特性计算的新方法 24-28 1 第32卷, 第5期 2 * |
《振动与冲击》 20030228 吴国荣等 用于结构振动模态分析的改进的有限元法 5-8 1 第22卷, 第2期 2 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104636556B (zh) * | 2015-02-09 | 2017-12-29 | 武汉理工大学 | 成任意角度连接的有限尺寸板结构振动响应计算方法 |
CN104636556A (zh) * | 2015-02-09 | 2015-05-20 | 武汉理工大学 | 成任意角度连接的有限尺寸板结构振动响应计算方法 |
CN104918184A (zh) * | 2015-04-14 | 2015-09-16 | 哈尔滨工程大学 | 一种耦合声场的声学耦合方法 |
CN104918184B (zh) * | 2015-04-14 | 2018-12-25 | 哈尔滨工程大学 | 一种耦合声场的声学耦合方法 |
CN105184060A (zh) * | 2015-08-24 | 2015-12-23 | 哈尔滨工程大学 | 一种基于空间配面与改进傅里叶级数的层合结构振动分析方法 |
CN105184060B (zh) * | 2015-08-24 | 2018-12-07 | 哈尔滨工程大学 | 一种基于空间配面与改进傅里叶级数的层合结构振动分析方法 |
CN105426343A (zh) * | 2015-11-02 | 2016-03-23 | 株洲时代新材料科技股份有限公司 | 一种基于傅里叶级数的复杂结构振动解析分析方法 |
CN106294975A (zh) * | 2016-08-05 | 2017-01-04 | 大连理工大学 | 一种基于降阶模型的梁式结构自由振动分析方法 |
CN106294975B (zh) * | 2016-08-05 | 2019-01-01 | 大连理工大学 | 一种基于降阶模型的梁式结构自由振动分析方法 |
CN106407582A (zh) * | 2016-09-28 | 2017-02-15 | 中国工程物理研究院总体工程研究所 | 一种用于柔性机构静力学和动力学分析的高效计算方法 |
CN106407582B (zh) * | 2016-09-28 | 2019-10-01 | 中国工程物理研究院总体工程研究所 | 一种用于柔性机构静力学和动力学分析的高效计算方法 |
CN106844875A (zh) * | 2016-12-28 | 2017-06-13 | 湖南大学 | 一种基于傅里叶级数的高速凸轮优化设计方法 |
CN106844875B (zh) * | 2016-12-28 | 2020-02-18 | 湖南大学 | 一种基于傅里叶级数的高速凸轮优化设计方法 |
CN107808048A (zh) * | 2017-10-26 | 2018-03-16 | 哈尔滨工程大学 | 一种均匀截面梁结构的振动分析方法 |
CN107808048B (zh) * | 2017-10-26 | 2020-12-25 | 哈尔滨工程大学 | 一种均匀截面梁结构的振动分析方法 |
CN109948180A (zh) * | 2019-01-25 | 2019-06-28 | 北京航空航天大学 | 一种正交各向异性对边简支矩形薄板振动分析方法 |
CN113158400A (zh) * | 2020-12-14 | 2021-07-23 | 东北电力大学 | 一种利用改进耦合能预测封闭箱体结构振动特性方法 |
CN113158400B (zh) * | 2020-12-14 | 2023-03-24 | 东北电力大学 | 一种利用改进耦合能预测封闭箱体结构振动特性方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101887474B (zh) | 2012-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101887474B (zh) | 基于有限元法与广义傅里叶级数法的结构振动分析方法 | |
Chen et al. | A local adaptive sampling method for reliability-based design optimization using Kriging model | |
Martínez-Ferrer et al. | An efficient finite-volume method to study the interaction of two-phase fluid flows with elastic structures | |
CN103838913B (zh) | 曲线箱梁弯桥的有限单元法 | |
Sevim et al. | Water length and height effects on the earthquake behavior of arch dam-reservoir-foundation systems | |
Liu et al. | Random vibration study of functionally graded porous curved beams with elastically restrained ends | |
CN106295028A (zh) | 一种局部结构动力学建模方法及装置 | |
CN106354954A (zh) | 一种基于叠层基函数的三维力学模态仿真模拟方法 | |
CN112949065A (zh) | 模拟层状岩体力学行为的双尺度方法、装置、存储介质及设备 | |
CN115081148A (zh) | 一种基于势能理论的加筋板等效参数确定方法 | |
CN101261651A (zh) | 主梁断面气动自激力的全过程数值模拟方法 | |
Williams et al. | Towards deep and simple understanding of the transcendental eigenproblem of structural vibrations | |
JP2005050063A (ja) | 流体の数値計算方法及び装置及び制御プログラム及び記憶媒体 | |
Righi et al. | ROM-based uncertainties quantification of flutter speed prediction of the BSCW wing | |
Zhang et al. | A stochastic wavelet finite element method for 1d and 2d structures analysis | |
CN115906333A (zh) | 一种桁架结构的几何非线性等效板动力学建模与响应分析方法 | |
CN104239588A (zh) | 一种基于超单元构造方法的结构功能分析方法 | |
Yan et al. | Dynamic Characteristics of Vibrating Screen with Determinate Structure and Statically Indeterminate Structure | |
CN115310188A (zh) | 一种基于最小势能理论确定加筋板等效参数的方法 | |
Malekzadeh et al. | Improvement of Numerical Manifold Method using Nine-node Quadrilateral and Ten-node Triangular Elements along with Complex Fourier RBFs in Modeling Free and Forced Vibrations | |
Wang et al. | In-plane free vibration of inhomogeneous curved beam with variable curvature under elastic constraints | |
Jiang et al. | Reduced multiscale computation on adapted grid for the convection-diffusion Robin problem | |
CN101276381A (zh) | 一种参数化的静态超单元构造方法 | |
Hatami et al. | Nonlinear analysis of axially moving plates using FEM | |
Wang et al. | Stress characteristics and overload failure analysis of cemented sand and gravel dam in Naheng Reservoir |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120201 Termination date: 20170625 |
|
CF01 | Termination of patent right due to non-payment of annual fee |