CN113158400A - 一种利用改进耦合能预测封闭箱体结构振动特性方法 - Google Patents

一种利用改进耦合能预测封闭箱体结构振动特性方法 Download PDF

Info

Publication number
CN113158400A
CN113158400A CN202011470175.0A CN202011470175A CN113158400A CN 113158400 A CN113158400 A CN 113158400A CN 202011470175 A CN202011470175 A CN 202011470175A CN 113158400 A CN113158400 A CN 113158400A
Authority
CN
China
Prior art keywords
energy
box structure
closed box
coupling
vibration
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
Application number
CN202011470175.0A
Other languages
English (en)
Other versions
CN113158400B (zh
Inventor
石双霞
刘孔超
王雨轩
肖斌
王靖宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeast Electric Power University
Original Assignee
Northeast Dianli University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeast Dianli University filed Critical Northeast Dianli University
Priority to CN202011470175.0A priority Critical patent/CN113158400B/zh
Publication of CN113158400A publication Critical patent/CN113158400A/zh
Application granted granted Critical
Publication of CN113158400B publication Critical patent/CN113158400B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141Discrete Fourier transforms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Discrete Mathematics (AREA)
  • Operations Research (AREA)
  • Buildings Adapted To Withstand Abnormal External Influences (AREA)

Abstract

本发明提供了一种利用改进耦合能预测封闭箱体结构振动特性方法,该方法包括以下步骤:由六个各向同性的弹性板构成封闭箱体结构;引用二维改进傅里叶级数方法给出弯曲振动位移与面内振动位移的表达形式;建立封闭箱体结构的能量函数;封闭箱体结构的能量由箱体结构系统总动能、总势能所组成,由建立不同坐标系下的耦合能公式得出耦合能;将结构弯曲振动位移及面内振动位移代入所构建的能量函数中,并采用瑞利里兹法对封闭箱体结构能量函数中未知傅里叶系数取极值获得封闭箱体结构控制方程;求解控制方程的特征值及特征向量获得封闭箱体结构的固有频率及振型。

Description

一种利用改进耦合能预测封闭箱体结构振动特性方法
技术领域
本发明涉及机械工程技术领域,具体涉及一种应用耦合能模型预测箱体振动特性的方法。
背景技术
封闭箱体结构是航空、土木和海洋工程等领域中应用较为广泛的载体结构。由于结构的振动特性对箱体内设备的正常运行有重要影响,因此预测并控制箱体结构振动是保证箱体结构内设备稳定运行的前提。由组合结构而构成的船舶舱室往往处于低频噪声环境,对长期处于舱室内工作人员精神及身体会产生一定危害。因此了解结构振动特性、提高舱室舒适性是船舶舱室振动噪声控制中急需解决问题。现有的借助有限元软件通过仿真分析得到箱体振动特性的步骤较为繁琐、耗时较长,计算精度受到网格细密程度的影响,而且受样本点数量限制,无法对箱体振动特性进行准确全面的预测。现有技术所给出的耦合结构的耦合能公式只适用于坐标系为同方向的组合结构,当组合结构中涉及到不同方向坐标系时,此耦合能公式不符合需求。
发明内容
有鉴于此,本发明针对现有技术存在的限制,提供了一种利用改进耦合能预测封闭箱体结构振动特性方法。本发明的技术方案为:
一种利用改进耦合能预测封闭箱体结构振动特性方法,该方法包括以下步骤:
步骤1:由六个各向同性的弹性板构成封闭箱体结构;
步骤2:引用二维改进傅里叶级数方法给出弯曲振动位移与面内振动位移的表达形式;
步骤3:建立封闭箱体结构的能量函数;封闭箱体结构的能量由箱体结构系统总动能、总势能所组成;
步骤4:将结构弯曲振动位移及面内振动位移代入所构建的能量函数中,并采用瑞利里兹法对封闭箱体结构能量函数中未知傅里叶系数取极值获得封闭箱体结构控制方程;
步骤5:求解控制方程的特征值及特征向量获得封闭箱体结构的固有频率及振型。
优选的是,步骤2的具体过程为:弹性约束边界条件下矩形板结构的面内振动位移和横向弯曲振动位移采用二维改进傅立叶级数方法进行求解;所采用的二维改进傅里叶级数表达式由标准二维傅里叶余弦级数和补充函数所构成,其中,补充函数由辅助函数与一维傅立叶余弦级数的乘积而得。封闭箱体结构中各个板的弯曲振动位移wi(xi,yi)、x方向面内振动位移ui(xi,yi)及y方向面内振动位移vi(xi,yi)的表达形式如下:
弯曲振动位移:
Figure RE-GDA0003086438090000021
式中Lxi为第i个板结构沿x方向长度,Lyi为第i个板结构沿y方向长度。 m,n为正整数。(xi,yi)为第i个板的位置坐标。λLxim=mπ/Lxi,λLyin=nπ/Lyi。Aimn, aim,bim,cim,dim,ein,fin,gin,hin为第i个板结构弯曲振动位移内未知的傅里叶系数。补充函数中辅助函数的具体表形式如下:
Figure RE-GDA0003086438090000022
Figure RE-GDA0003086438090000023
Figure RE-GDA0003086438090000031
Figure RE-GDA0003086438090000032
Figure RE-GDA0003086438090000033
Figure RE-GDA0003086438090000034
Figure RE-GDA0003086438090000035
Figure RE-GDA0003086438090000036
x方向面内振动位移:
Figure RE-GDA0003086438090000037
y方向面内振动位移:
Figure RE-GDA0003086438090000038
面内振动位移表达式(式(2)-(3))中辅助函数ξ1Lxi(xi)和ξ2Lxi(xi)为:
Figure RE-GDA0003086438090000039
Figure RE-GDA00030864380900000310
式(4)中
Figure RE-GDA00030864380900000311
通过将辅助函数ξ1Lxi(xi)和ξ2Lxi(xi)中x与Lxi替换为y和Lyi可得面内振动位移表达式(式(2)-(3))中辅助函数ξ1Lyi(yi)和ξ2Lyi(yi)。x方向面内振动位移中Bimn,ain in,bin in,cin im,din im为第i个板内未知的傅里叶系数。y方向面内振动位移中Cimn,ein in,fin in,gin im,hin im为第i个板内未知的傅里叶系数。
优选的是,封闭箱体结构总势能包含6个板结构的弯曲振动势能、6个板结构的面内振动势能及12处结构之间的耦合能,封闭箱体结构总动能应包含6个板结构弯曲振动动能及6个板结构面内振动动能;
采用能量原理构建箱体结构的能量函数L:
L=V-T (5)
式(5)中V代表封闭箱体结构的总势能,T代表封闭箱体结构的总动能,
总势能:
Figure RE-GDA0003086438090000041
总动能:
Figure RE-GDA0003086438090000042
其中i代表箱体结构中6个板的序号,j代表箱体结构中12处耦合能的序号,Vi bending代表第i个板的弯曲振动势能,Vi in-panel代表第i个板的面内振动势能, Vj coupling代表第j处结构之间的耦合能,Ti bending代表第i个板的弯曲振动动能, Ti in-panel代表第i个板的面内振动动能。
优选的是,建立不同坐标系下的耦合能公式;建立的多种方向坐标系,根据所建立的坐标系给出12处耦合能表达形式。
优选的是,将结构弯曲振动位移表达式(1)、结构面内振动位移表达式(2) 与(3)代入至封闭箱体结构系统的拉格朗日函数表达形式,采用Rayleigh-Ritz 方法对拉格朗日函数中每个未知Fourier系数取极值可得到一系列线性方程组,其矩阵表达形式为:
(K-ω2M)E=0 (8)
矩阵表达形式中ω为封闭箱体结构角频率,K为封闭箱体结构系统的刚度矩阵,M为箱体结构的质量矩阵,E为箱体结构位移表达式的未知系数。对箱体结构系统控制方程(式(8))进行特征值求解,可获得封闭箱体结构系统特征值(ω)及特征向量(E)。封闭箱体结构系统固有频率利用特征值可得 (ω/(2π)),其模态振型可通过将所对应特征向量代入至位移表达式中获得。
与现有技术相比,本发明可以获得包括以下技术效果:
本发明采用由六个各向同性的弹性板所构成的封闭箱体模型,采用改进的傅里叶级数形式进行结构位移的表达,主要采用能量原理构建箱体结构能量函数,将结构弯曲振动位移及面内振动位移代入所构建的能量函数中,并采用瑞利里兹法对封闭箱体结构能量函数中未知傅里叶系数取极值获得封闭箱体结构控制方程,通过求解控制方程的特征值及特征向量进而获得封闭箱体结构的固有频率及振型,过程简便、效率高、计算精度高,通过商业软件有限元方法进行模态比较,验证其封闭箱体结构固有频率和振型,良好吻合性验证理论模型预测的准确性。
为了解决组合板结构的多处耦合问题,应将原有耦合能公式进行修改以适应不同坐标系下耦合结构振动问题分析。由于结构之间耦合需要,封闭箱体结构模型中所建立的多种方向坐标系使得封闭箱体结构具有多种形式耦合能。对原有耦合能公式进行修改建立了不同坐标系下的耦合能公式,根据封闭箱体结构模型所建立的坐标系给出12处耦合能表达形式,以供利用改进耦合能预测封闭箱体结构振动特性。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1是封闭箱体结构系统示意图;
图2(a)现有技术中的相同方向坐标系构成的耦合结构;
图2(b)本发明中的不同方向坐标系构成的耦合结构;
图3是封闭箱体结构7-13阶模态振型对比;
所示附图中1:1板,2:2板,3:3板,4:4板,5:5板,6;6板。
图3中a1-f1为本发明得到的封闭箱体结构模态振型,a2-f2为有限元方法得到的封闭箱体结构模态振型。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式,熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面结合附图和具体的实施例对本发明做进一步详细说明,所述是对本发明的解释而不是限定。
一种利用改进耦合能预测封闭箱体结构振动特性方法,该方法包括以下步骤:
步骤1:构建封闭箱体结构模型,构建的封闭箱体结构模型如图1所示,该封闭箱体结构模型由六个各向同性的弹性板所构成;
步骤2:引用二维改进傅里叶级数方法给出弯曲振动位移与面内振动位移的表达形式;
步骤2的具体过程为:弹性约束边界条件下矩形板结构的面内振动位移和横向弯曲振动位移采用二维改进傅立叶级数方法进行求解;所采用的二维改进傅里叶级数表达式由标准二维傅里叶余弦级数和补充函数所构成,其中,补充函数由辅助函数与单维傅立叶余弦级数的乘积而得,引进的补充函数可克服二维傅里叶余弦级数在每一边界上面内或弯曲振动位移一阶求导的不连续性。封闭箱体结构中各个板的弯曲振动位移wi(xi,yi)、x方向面内振动位移ui(xi,yi)及 y方向面内振动位移vi(xi,yi)的表达形式如下:
弯曲振动位移:
Figure RE-GDA0003086438090000071
式中Lxi为第i个板结构沿x方向长度,Lyi为第i个板结构沿y方向长度。 m,n为正整数。(xi,yi)为第i个板的位置坐标。λLxim=mπ/Lxi,λLyin=nπ/Lyi。Aimn, aim,bim,cim,dim,ein,fin,gin,hin为第i个板结构弯曲振动位移内未知的傅里叶系数。补充函数中辅助函数的具体表达形式如下:
Figure RE-GDA0003086438090000072
Figure RE-GDA0003086438090000073
Figure RE-GDA0003086438090000074
Figure RE-GDA0003086438090000075
Figure RE-GDA0003086438090000076
Figure RE-GDA0003086438090000077
Figure RE-GDA0003086438090000078
Figure RE-GDA0003086438090000079
x方向面内振动位移:
Figure RE-GDA00030864380900000710
y方向面内振动位移:
Figure RE-GDA0003086438090000081
面内振动位移表达式(式(2)-(3))中辅助函数ξ1Lxi(xi)和ξ2Lxi(xi)为:
Figure RE-GDA0003086438090000084
Figure RE-GDA0003086438090000085
式(4a)、(4b)中
Figure RE-GDA0003086438090000086
通过将辅助函数ξ1Lxi(xi)和ξ2Lxi(xi)中x与Lxi替换为y和Lyi可得面内振动位移表达式(式(2)-(3))中辅助函数ξ1Lyi(yi)和ξ2Lyi(yi)。x 方向面内振动位移中Bimn,ain in,bin in,cin im,din im为第i个板内未知的傅里叶系数。 y方向面内振动位移中Cimn,ein in,fin in,gin im,hin im为第i个板内未知的傅里叶系数。
步骤3:采用能量原理构建箱体结构的能量函数:
L=V-T (5)
V代表封闭箱体结构总势能,T代表封闭箱体结构总动能。封闭箱体结构总势能包含6个板结构的弯曲振动势能、6个板结构的面内振动势能及12处结构之间的耦合能,封闭箱体结构总动能应包含6个板结构弯曲振动动能及6个板结构面内振动动能;
总势能:
Figure RE-GDA0003086438090000082
总动能:
Figure RE-GDA0003086438090000083
其中i代表箱体结构中6个板的序号,j代表箱体结构中12处耦合能的序号,Vi bending代表第i个板的弯曲振动势能,Vi in-panel代表第i个板的面内振动势能, Vj coupling代表第j处结构之间的耦合能,Ti bending代表第i个板的弯曲振动动能, Ti in-panel代表第i个板的面内振动动能。
其中箱体结构中第i个板的弯曲振动势能Vi bending如式(8)所示,箱体结构中第i个板的弯曲振动动能Ti bending如式(9)所示,箱体结构中第i个板的面内振动势能Vi in-panel如式(10)所示,箱体结构中第i个板的面内振动动能Ti in-panel如式 (11)所示。
Figure RE-GDA0003086438090000091
Figure RE-GDA0003086438090000092
Figure RE-GDA0003086438090000093
Figure RE-GDA0003086438090000094
其中ω为圆频率;ρ为质量密度;t为时间;h为板结构厚度;D为板结构弯曲刚度;G为板结构拉伸刚度;μ为板结构材料泊松比;i为板的序号;j为(0,1) 表示坐标的起点和终点;Kbxij为x方向弯曲振动分量的旋转约束边界弹簧的(单位长度上)刚度系数;Kbyij为y方向弯曲振动分量的旋转约束边界弹簧的(单位长度上)刚度系数;kbxij为x方向弯曲振动分量的横向支撑边界弹簧的(单位长度上)刚度系数;kbyij为y方向弯曲振动分量的横向支撑边界弹簧的(单位长度上)刚度系数;knxij为x方向面内振动在法向的边界约束弹簧的(单位长度上)刚度系数;knyij为y方向面内振动在法向的边界约束弹簧的(单位长度上)刚度系数;kpxij为x方向面内振动在切向的边界约束弹簧的(单位长度上)刚度系数;kpyij为y方向面内振动在切向的边界约束弹簧的(单位长度上) 刚度系数;
图2(a)为任意边界与连接条件下耦合板结构分析模型中具有相同方向坐标的耦合结构。图2(b)为任意边界与连接条件下耦合板结构分析模型中具有不同方向坐标的耦合结构。图2中的耦合结构由公共边界连接的板1和板2组合而成,其中,板1位于x1-y1平面。弯曲振动分量边界条件通常可以用沿边界均匀分布的横向支撑kbx11和旋转约束弹簧Kbx11来模拟。类似地,面内振动分量的任意边界条件可以用在x1-y1平面中沿边界均匀分布的另外两种约束弹簧 (knx11和kpx11)来模拟。图2(a)中耦合板结构由公共边界(x1=0或x2=Lx2)连接的板 1和板2组合而成,坐标系的选取如图2(a)所示。图2(b)耦合板结构由公共边界(y1=0或x2=Lx2)连接的板1和板2组合而成,坐标系的选取如图2(b)。图2(a)或图2(b)中各板长宽分别为Lxi和Lyi(i=1或2)。以x1-y1-z1为全局坐标时,两板之间的相对位置可以通过定义耦合夹角θ表达出来。在结构连接处,通过沿耦合边界引入四类耦合弹簧(分别用Kc、kc1、kc2和kc3来表示),进而全面考虑了横向弯矩、面外剪切、面内纵向作用(垂直于耦合边界)和面内剪切(相切于耦合边界)四种耦合效应。
图2(a)与图2(b)内耦合结构中四类耦合弹簧所存储的弹性势能分别为式(12)与式(13):
Figure RE-GDA0003086438090000101
Figure RE-GDA0003086438090000111
式(12)-(13)中wi为沿zi轴方向的弯曲振动位移,ui和vi分别为沿x轴和y 轴方向的面内振动位移。Vcoupling1由四部分构成,第一项由考虑横向弯矩弹簧Kc 与耦合处角度差
Figure RE-GDA0003086438090000112
所提供的耦合势能,第二项由考虑面外剪切弹簧 kc1与耦合处横向振动位移差(w1-(w2cosθ-u2sinθ))所提供的耦合势能,第三项由考虑面内纵向作用弹簧kc2与耦合处x方向面内振动位移差(u1-(u2cosθ+w2sinθ)) 所提供的耦合势能,第四项由考虑面内剪切弹簧与耦合处y方向面内振动位移差(v1-v2)所提供的耦合势能。
式(12)所提出的耦合能表达方式仅适用于具有同方向坐标系的耦合结构 (如图2(a))。若耦合结构内坐标系(如图2(b))具有不同方向,式(12) 所给的耦合能不再适用,因此本文对已有耦合能(式(12))进行修改,给出了具有不同坐标方向的耦合结构(图2(b))的耦合能,其表达式如式(13)所示。式(13)所给出的耦合能Vcoupling2由四部分构成,第一项由考虑横向弯矩弹簧 Kc与耦合处角度差
Figure RE-GDA0003086438090000113
所提供的耦合能,第二项由考虑面外剪切弹簧 kc1与耦合处横向振动位移差(w1-(w2cosθ-u2sinθ))所提供的耦合能,第三项由考虑面内纵向作用弹簧kc2与耦合处x方向面内振动位移差(v1-(u2cosθ+w2sinθ)) 所提供的耦合能,第四项由考虑面内剪切弹簧y方向面内振动位移差(u1-(-v2)) 所提供的耦合能。根据式(13)的构建方式,本文给出了封闭箱体结构(图1) 结构内耦合能的表达形式。
发明者根据图1的封闭箱体结构所建立的坐标系给出12处耦合能表达式。封闭箱体结构图1中1板与5板的耦合能如式(14)所示,2板与5板的耦合能如式(15)所示,3板与5板的耦合能如式(16)所示,4板与5板的耦合能如式(17)所示,1板与6板的耦合能如式(18)所示,2板与6板的耦合能如式 (19)所示,3板与6板的耦合能如式(20)所示,4板与6板的耦合能如式(21) 所示,1板与2板的耦合能如式(22)所示,2板与3板的耦合能如式(23)所示,3板与4板的耦合能如式(24)所示,4板与1板的耦合能如式(25)所示。
Figure RE-GDA0003086438090000121
Figure RE-GDA0003086438090000122
Figure RE-GDA0003086438090000123
Figure RE-GDA0003086438090000131
Figure RE-GDA0003086438090000132
Figure RE-GDA0003086438090000133
Figure RE-GDA0003086438090000134
Figure RE-GDA0003086438090000141
Figure RE-GDA0003086438090000142
Figure RE-GDA0003086438090000143
Figure RE-GDA0003086438090000144
Figure RE-GDA0003086438090000151
式(14)-(25)中Lx、Ly和Lz分别代表图1中封闭箱体结构的几何尺寸。xi、 yi及zi代表封闭箱体结构中第i个板结构的位置坐标。θ代表图1封闭箱体结构内12组耦合板结构之间的夹角,为90度。耦合弹簧分别用Kc、kc1、kc2和kc3表示。
步骤4:将结构弯曲振动位移及面内振动位移代入所构建的能量函数中,并采用瑞利里兹法对封闭箱体结构能量函数中未知系数取极值获得箱体结构控制方程;将结构位移表达式(式(1)-(3))代入至箱体结构系统的拉格朗日函数表达式(式(5)),采用Rayleigh-Ritz方法对拉格朗日函数中每个未知Fourier 系数取极值可得到一系列线性方程组,其矩阵表达形式为:
(K-ω2M)E=0 (26)
式(26)中K为封闭箱体结构系统的刚度矩阵,M为封闭箱体结构的质量矩阵,E为封闭箱体结构位移表达式中未知傅里叶系数。对封闭箱体结构系统控制方程(式(26))进行特征值求解,可获得封闭箱体结构系统特征值(ω)及特征向量(E)。箱体结构系统固有频率利用特征值可得(ω/(2π)),其模态振型可通过将所对应特征向量代入至位移表达式中获得。
使用商业软件有限元模型进行数值验证
以图1所示的由6个板耦合而成的封闭箱体结构系统示意图为例。因各个弹性板物理参数(密度,厚度,杨氏模量,泊松比)相同,组合结构中每个弹性板的势能与动能的计算方式不变(式(8)-(11))。根据所设计的6种坐标系组成的封闭箱体结构存在12处结构耦合,其耦合能为式(14)-(25)。假定算例中板具有同样厚度,即h=0.008m。几何参数Lx=1.4m,Ly=1.2m,Lz=1.0 m。板材质均为钢,参数如下:杨氏模量E=2.16×1011N/m2,质量密度ρ=7800 kg/m3,泊松比μ=0.28。将其几何参数及物理参数代入至封闭箱体结构拉格朗日函数中,并通过瑞利里兹法获得封闭箱体结构特征方程,最终获得该封闭箱体的固有频率及模态振型。
表1封闭箱体结构前15阶固有频率(Hz)
Figure RE-GDA0003086438090000161
表1为封闭箱体结构前十五阶频率,通过比较可以看出截断级数为12时,封闭箱体结构频率已开始收敛,并与有限元方法所获得固有频率结果吻合较好。因此通过该方法所得到的数值可以表明该方法能够快速收敛并具有高的计算精度。图3为封闭箱体结构系统第7阶-第15阶模态振型,通过与有限元方法所获得模态振型相比较,不难发现本发明所获得模态振型与有限元方法所获得模态振型吻合很好。通过商业软件有限元方法进行模态比较,验证其箱体结构固有频率和振型,良好吻合性验证理论模型预测的准确性。
本发明采用由六个各向同性的弹性板所构成的箱体模型,采用改进的傅里叶级数形式进行位移表达,主要采用能量原理构建箱体结构能量函数,将结构弯曲振动位移及面内振动位移代入所构建的能量函数中,并采用瑞利里兹法对箱体结构能量函数中未知系数取极值获得箱体结构控制方程,通过求解控制方程的特征值及特征向量进而获得箱体结构的固有频率及振型,过程简便、效率高、计算精度高,通过商业软件有限元方法进行模态比较,验证其箱体结构固有频率和振型,良好吻合性验证理论模型预测的准确性。
如在通篇说明书及权利要求当中所提及的“包含”为一开放式用语,故应解释成“包含但不限定于”。“大致”是指在可接收的误差范围内,本领域技术人员能够在一定误差范围内解决所述技术问题,基本达到所述技术效果。说明书后续描述为实施本申请的较佳实施方式,然所述描述乃以说明本申请的一般原则为目的,并非用以限定本申请的范围。本申请的保护范围当视所附权利要求所界定者为准。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的商品或者系统不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种商品或者系统所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的商品或者系统中还存在另外的相同要素。
上述说明示出并描述了发明的若干优选实施例,但如前所述,应当理解发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述发明构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离发明的精神和范围,则都应在发明所附权利要求的保护范围内。

Claims (4)

1.一种利用改进耦合能预测封闭箱体结构振动特性方法,其特征在于,该方法包括以下步骤:
步骤1:由六个各向同性的弹性板构成封闭箱体结构;
步骤2:引用二维改进傅里叶级数方法给出弯曲振动位移与面内振动位移的表达形式;
步骤3:建立封闭箱体结构的能量函数;封闭箱体结构的能量由箱体结构系统总动能、总势能所组成,其中箱体结构总势能包含6个板结构的弯曲振动势能、面内振动势能及12处结构之间的耦合能,所述耦合能由建立不同坐标系下的耦合能公式得出;
步骤4:将结构弯曲振动位移及面内振动位移代入所构建的能量函数中,并采用瑞利里兹法对封闭箱体结构能量函数中未知傅里叶系数取极值获得封闭箱体结构控制方程;
步骤5:求解控制方程的特征值及特征向量获得封闭箱体结构的固有频率及振型。
2.如权利要求1所述的一种利用改进耦合能预测封闭箱体结构振动特性方法,其特征在于,步骤2的具体过程为:弹性约束边界条件下矩形板结构的面内振动位移和横向弯曲振动位移采用二维改进傅立叶级数方法进行求解;所采用的二维改进傅里叶级数表达式由标准二维傅里叶余弦级数和补充函数所构成,其中,补充函数由辅助函数与单傅立叶余弦级数的乘积而得。
3.如权利要求1所述的一种利用改进耦合能预测封闭箱体结构振动特性方法,其特征在于,箱体结构总势能包含6个板结构的弯曲振动势能、面内振动势能及12处结构之间的耦合能,箱体结构总动能应包含6个板结构弯曲振动动能及面内振动动能;
采用能量原理构建箱体结构的能量函数:
L=V-T (5)
V代表封闭箱体结构总势能,T代表封闭箱体结构总动能。
总势能:
Figure FDA0002833321870000021
总动能:
Figure FDA0002833321870000022
其中:i-箱体结构中6个板的序号,j-箱体结构中12处耦合能的序号。
4.如权利要求3所述的一种利用改进耦合能预测封闭箱体结构振动特性方法,其特征在于,根据所建立的坐标系给出12处耦合能表达形式,所给出的耦合结构的耦合能公式如下:
Figure FDA0002833321870000023
Figure FDA0002833321870000024
Figure FDA0002833321870000031
Figure FDA0002833321870000032
Figure FDA0002833321870000033
Figure FDA0002833321870000034
Figure FDA0002833321870000041
Figure FDA0002833321870000042
Figure FDA0002833321870000043
Figure FDA0002833321870000044
Figure FDA0002833321870000051
Figure FDA0002833321870000052
CN202011470175.0A 2020-12-14 2020-12-14 一种利用改进耦合能预测封闭箱体结构振动特性方法 Active CN113158400B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011470175.0A CN113158400B (zh) 2020-12-14 2020-12-14 一种利用改进耦合能预测封闭箱体结构振动特性方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011470175.0A CN113158400B (zh) 2020-12-14 2020-12-14 一种利用改进耦合能预测封闭箱体结构振动特性方法

Publications (2)

Publication Number Publication Date
CN113158400A true CN113158400A (zh) 2021-07-23
CN113158400B CN113158400B (zh) 2023-03-24

Family

ID=76882536

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011470175.0A Active CN113158400B (zh) 2020-12-14 2020-12-14 一种利用改进耦合能预测封闭箱体结构振动特性方法

Country Status (1)

Country Link
CN (1) CN113158400B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117171864A (zh) * 2023-11-02 2023-12-05 鲁东大学 一种梁结构线性振动预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101887474A (zh) * 2010-06-25 2010-11-17 哈尔滨工程大学 基于有限元法与广义傅里叶级数法的结构振动分析方法
CN110309579A (zh) * 2019-06-27 2019-10-08 复旦大学 一种针对弹性飞机阵风响应的仿真分析方法和系统
CN111382503A (zh) * 2020-02-27 2020-07-07 中南大学 弹性支承下旋转的柔性圆环的振动分析方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101887474A (zh) * 2010-06-25 2010-11-17 哈尔滨工程大学 基于有限元法与广义傅里叶级数法的结构振动分析方法
CN110309579A (zh) * 2019-06-27 2019-10-08 复旦大学 一种针对弹性飞机阵风响应的仿真分析方法和系统
CN111382503A (zh) * 2020-02-27 2020-07-07 中南大学 弹性支承下旋转的柔性圆环的振动分析方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DU J T 等: "Free vibration of two elastically coupled rectangular plates with uniform elastic boundary restraints", 《JOURNAL OF SOUND AND VIBRATION》 *
郭旗: "L型板—声腔系统动力学特性分析及实验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117171864A (zh) * 2023-11-02 2023-12-05 鲁东大学 一种梁结构线性振动预测方法
CN117171864B (zh) * 2023-11-02 2024-02-20 鲁东大学 一种梁结构线性振动预测方法

Also Published As

Publication number Publication date
CN113158400B (zh) 2023-03-24

Similar Documents

Publication Publication Date Title
Ferreira et al. Buckling and vibration analysis of isotropic and laminated plates by radial basis functions
Singha et al. Free vibration of rotating pretwisted FG-GRC sandwich conical shells in thermal environment using HSDT
Ma et al. Modeling complex crack problems using the numerical manifold method
Sadamoto et al. An effective meshfree reproducing kernel method for buckling analysis of cylindrical shells with and without cutouts
Wu Thermomechanical properties of fiber-reinforced plastic (FRP) bars
CN107220461A (zh) 一种变刚度复合材料板壳结构高效优化方法
Zhuang et al. Modal and aeroelastic analysis of trapezoidal corrugated-core sandwich panels in supersonic flow
Zheng et al. A bond‐based peridynamic model considering effects of particle rotation and shear influence coefficient
Zhao et al. Piecewise shear deformation theory and finite element formulation for vibration analysis of laminated composite and sandwich plates in thermal environments
Rezaiee-Pajand et al. Analyzing FG shells with large deformations and finite rotations
Sharma et al. Numerical analysis of acoustic radiation responses of shear deformable laminated composite shell panel in hygrothermal environment
Lal et al. Stochastic fracture analysis of laminated composite plate with arbitrary cracks using X-FEM
Swaminathan et al. Vibration and stability characteristics of functionally graded sandwich plates with/without porosity subjected to localized edge loadings
CN113158400B (zh) 一种利用改进耦合能预测封闭箱体结构振动特性方法
Ni et al. Free vibration of stiffened laminated shells of revolution with a free-form meridian and general boundary conditions
Chhipa et al. Removing fiber orientation uncertainty from the finite element model of a composite lamina with direct updating algorithm
Pal et al. Free vibration analysis of laminated composite plate with cut-out carrying concentrated and distributed mass
Ma et al. Nonlinear dynamic response of a stiffened plate with four edges clamped under primary resonance excitation
Wang et al. Nitsche-based isogeometric approach for free vibration analysis of laminated plate with multiple stiffeners and cutouts
Zinno et al. A three‐dimensional layer‐wise constant shear element for general anisotropic shell‐type structures
Zhou et al. Three-dimensional vibration analysis of a torus with circular cross section
Sylvia et al. A dynamic model of a reinforced thin plate with ribs of finite width
Aktas et al. Experimental confirmation for the validity of Ritz method in structural dynamic analysis
Anish et al. Rhombic laminates with mass variations under dual-axis compression
Yusa et al. Analysis of a many-hole problem using coupling-matrix-free iterative s-version FEM with multiple local meshes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant