CN107808048B - 一种均匀截面梁结构的振动分析方法 - Google Patents
一种均匀截面梁结构的振动分析方法 Download PDFInfo
- Publication number
- CN107808048B CN107808048B CN201711016927.4A CN201711016927A CN107808048B CN 107808048 B CN107808048 B CN 107808048B CN 201711016927 A CN201711016927 A CN 201711016927A CN 107808048 B CN107808048 B CN 107808048B
- Authority
- CN
- China
- Prior art keywords
- beam structure
- matrix
- displacement
- section
- core
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Geometry (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Data Mining & Analysis (AREA)
- Civil Engineering (AREA)
- Evolutionary Computation (AREA)
- Architecture (AREA)
- Structural Engineering (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Bridges Or Land Bridges (AREA)
Abstract
一种均匀截面梁结构的振动分析方法,包括以下步骤:提取梁结构的几何尺寸特性,将梁结构的位移场分离为截面面内位移和轴向位移;对梁结构截面面内位移进行拟合;对梁结构轴向位移进行延展;结合梁结构截面面内位移和轴向位移,计算得到梁结构的整体位移;计算弹性板结构的应变向量和应力向量;计算梁结构的应变能、动能方程、边界能;建立结构拉格朗日能量泛函得到梁结构核心质量矩阵和刚度矩阵。通过迭代循环核心矩阵求得整体的质量矩阵、刚度矩阵以及总体质量矩阵,进而得到结构的特征方程;求解梁结构的特征方程矩阵并计算得到固有频率和振型。本发明方法适用多边界条件、多截面形状梁结构,且精度高、收敛快、计算成本低和计算方法简单。
Description
技术领域
本发明属于结构动力学领域,具体涉及一种均匀截面梁结构的振动分析方法。
背景技术
细长型结构广泛应用于船舶与海洋工程、建筑工程和航空航天等工程设备中,例如桥梁、杆系、轴系以及飞机机翼等。深入研究这类结构的振动特性及其参数影响规律在设备早期设计阶段振动噪声水平预估和实现低噪声设计具有重要的理论与实践指导意义。工程应用中常采用梁理论来处理这种细长型结构,以达到兼顾求解效率和计算精度的目的。
经典的梁理论主要有欧拉-伯努利和铁木申科两种理论,然而,欧拉-伯努利梁理论忽略了结构横向拉伸和剪切变形以及转动惯量的影响,仅适用于求解特定的细长梁结构;铁木申科梁理论虽然考虑了横向剪切变形和转动惯量的影响,但仍忽略了横向拉伸变形以及扭转时的翘曲变形,仅适用于短粗梁结构的横向弯曲问题。当梁结构的高跨比较大或者为薄壁结构时,这两种理论计算结果往往偏差较大甚至无法用于计算分析,此时必须采用新的三维梁理论。卡诺基于三维弹性理论提出了卡诺高阶截取拟合技术(相关文献为Carrera E.A class of two-dimensional theories for anisotropic multilayeredplates analysis.[J].Mem.accad.sci.torino Cl.sci.fis.mat.natur,1995:49-87.Carrera E.Theories and finite elements for multilayered,anisotropic,composite plates and shells[J].Archives of Computational Methods inEngineering,2002,9(2):87-140.),并将该技术成功地引入到梁结构的研究中(相关文献为Petrolo M,Zappino E,Carrera E.Refined free vibration analysis of one-dimensional structures with compact and bridge-like cross-sections[J].Thin-Walled Structures,2012,56(2):49-61.)。该技术可以通过改变位移中未知变量的个数而衍生出很多的梁理论模型,而这只需对其核心部分稍加修改便可得到。基于卡诺高阶截取拟合技术,梁理论由低阶变成高阶会变得十分容易,这一系列的梁理论都可以通过同一个矩阵核心迭代而来,变化的只是迭代的次数。另一方面,该技术中位移里未知变量的个数可以根据问题的需求而确定,针对特定的工程梁问题,可以通过选择合适的未知变量的数量而达到相应的计算精度。如薄壁梁、弯曲、扭转、剪切局部载荷,各种静态和动态问题。
然而,已有的卡诺高阶截取拟合技术大部分是通过嵌入有限元方法来求解梁结构动力学问题,改变梁截面形状或增加单元数量往往需要重新建模分析,另外有限元方法普遍具有计算量大,计算精度不高,边界条件施加较为繁琐等缺点。因此,研究和建立一种能够适用多截面形状、多边界条件梁结构的振动分析方法具有十分重要的意义。
发明内容
本发明的目的在于提供一种适用多边界条件、多截面形状梁结构的,且精度高、收敛快、计算成本低和计算方法简单的均匀截面梁结构的振动分析方法。
本发明的目的是这样实现的,包括以下步骤:
(1)提取梁结构的几何尺寸特性,将梁结构的位移场分离为截面面内位移和轴向位移,具体表达式为:
Uk(x,y,z)=Φk(x,z)Fk(y)
其中x和z为结构空间沿着梁截面方向的坐标,y为梁轴向坐标,Uk(x,y,z)为梁结构的整体位移,k=1,2,3分别表示x,y和z三个方向上的分量,Φk(x,z)为截面面内位移,Fk(y)为轴向位移。
(2)利用卡诺高阶截取技术对梁结构截面面内位移进行拟合,拟合形式如下:
(3)采用改进傅里叶级数对梁结构轴向位移进行延展,具体形式如下:
其中,M为改进傅里叶级数的截断项数,m=-2,-1,...,M;bkm代表相应改进傅里叶展开项的系数,L为梁结构的总长度;
(4)结合梁结构截面面内位移和轴向位移,计算得到梁结构的整体位移,具体表达式如下:
其中,U1(x,y,z),U2(x,y,z)和U3(x,y,z)分别对应空间坐标x,y和z三个方向上的位移分量,Aτm,Bτm和Cτm为位移分量中相应项的系数。
(5)计算弹性板结构的应变向量和应力向量;
所涉及弹性板结构的应变向量的表达式为:
ε=[εx,εy,εz,γxy,γyz,γxz]T
其中,ε表示弹性板结构的应变向量;上标T表示转置;εx,εy和εz为正应变分量;γxy,γyz和γxz为切应变分量,且有
所涉及应力向量的表达式为:
σ=Dε
其中,σ表示弹性板结构的应力向量,D为结构材料系数矩阵。
(6)计算梁结构的应变能和动能方程,设置虚拟弹簧边界从而获取边界能,具体表达式如下:
其中,Vs,Tp和Vp分别为梁结构的应变能、动能和边界能方程;t表示时间,ρ为材料的密度。kx0,ky0和kz0为梁结构y=0端所设的虚拟弹簧边界的刚度系数,kxL,kyL和kzL为梁结构y=L端所设的虚拟弹簧边界的刚度系数。
(7)建立结构拉格朗日能量泛函Ω=Vs+Vp-Tp,然后对系数Aτm,Bτm和Cτm求偏导并令其结果为零,即可得到梁结构的3×3阶核心质量矩阵和刚度矩阵。核心矩阵中的元素如下:
其中,Kmnij为核心刚度矩阵,Mmnij为核心质量矩阵;上角标i,j=1,…,Γ;m,n=-2,-1,0,…,M;下角标a,b和c表示核心矩阵中元素的标号。和表示分别对x和z求一阶偏导,和表示分别对x和z求一阶偏导,Ψm,y和Ψn,y分别表示Ψm和Ψn对y求一阶导数。D11,…,D66为结构材料系数矩阵D中的元素。
(8)通过迭代循环核心矩阵求得整体的质量矩阵、刚度矩阵以及总体质量矩阵M,进而得到结构的特征方程;
所述质量矩阵和刚度矩阵的求解方法为:指针i,j由1取到N+1循环核心刚度矩阵Kmnpqij得到子子矩阵Kmnpq,指针p,q由1取到P+3循环子子矩阵Kmnpq得到子矩阵Kmn,指针m,n由1取到M+3循环子矩阵Kmn得到总体刚度矩阵K;
所述结构的特征方程表达式为:
(K-ω2M)A=0
其中ω为圆频率,A为对应ω的特征向量;
(9)求解梁结构的特征方程矩阵并计算得到固有频率和振型。
本发明的优势在于:本发明实现了对梁理论阶次和轴向拟合基函数截断级数的参数化分析,改变梁理论阶次或增加截断项数,只需简单地增加核心矩阵的循环迭代次数;本发明中梁截面位移函数以二维泰勒多项式进行拟合,并无特殊针对性,理论上该方法适用于任意截面形状的均匀梁结构;本发明中梁轴向位移函数以傅里叶级数进行全求解域拟合,与传统有限元方法相比,具有收敛速度快计算精度高等特点;本发明只需要通过控制边界弹簧的刚度来满足结构的不同边界条件要求,而不需要对程序做任何修改。总的来说本发明的方法具有适用任意边界条件和任意截面形状、精度高、收敛快、计算成本低等特点。
附图说明
图1为本发明的流程图;
图2为方形空心梁结构示意图;
图3为刚度矩阵的装配示意图。
具体实施方式
为使本发明解决的技术问题、采用的技术方案和达到的技术效果更加清楚,下面结合附图对本发明做进一步描述。本发明方法的执行步骤如图1所示。
考虑一个方形空心截面梁结构,如图2所示,截面尺寸b1=h1=0,b2=h2=0.2m,梁材料为各向同性,杨氏模量E=75GPa,密度ρ=7800kg/m3,泊松比μ=0.3。梁的长度L为2m或20m。对应的高跨比h2/L分别为0.1和0.01,分别代表着短粗梁和细长梁结构。梁两端简支无外加载荷。利用本发明方法对其进行求解,二维泰勒展开阶次N=2。具体步骤如下:
(1)提取空心截面梁的几何尺寸,建立空间直角坐标系,根据梁结构长度方向尺寸远大于另外两个方向这一特性,将梁结构的位移场分离为截面面内位移和轴向位移,具体表达式为:
Uk(x,y,z)=Φk(x,z)Fk(y)
其中x和z为结构空间沿着梁截面方向的坐标,y为梁轴向坐标,Uk(x,y,z)为梁结构的整体位移,k=1,2,3分别表示x,y和z三个方向上的分量,Φk(x,z)为截面面内位移,Fk(y)为轴向位移。
(2)利用卡诺高阶截取技术对梁结构截面面内位移进行拟合,由于该梁结构截面较为复杂,故选取二维泰勒展开阶次N=2进行计算,泰勒展开式项数Γ=(N+1)(N+2)/2=6,具体拟合形式如下:
其中,为二维泰勒展开的第τ项,ak0,ak1,ak2,…,akΓ为相应泰勒展开项的系数,akτ为二维泰勒展开第τ项的系数,b1,b2,h1和h2为梁截面尺寸参数,N为二维泰勒展开的阶次,Γ为二维泰勒展开的总项数。
(3)采用改进傅里叶级数对梁结构轴向位移进行延展,初步选取截断项数M=10,具体形式如下:
其中bkm代表相应改进傅里叶展开项的系数,L为梁结构的总长度,两个sine项(m=-2和m=-1),的引入是为了消除结构位移展开成传统傅里叶cosine级数时其本身及一阶导数在两端的不连续性,从而加快方法的收敛速度。
(4)结合梁结构截面面内位移和轴向位移,并将系数进行整合,即可得到梁结构的整体位移,具体表达式如下:
其中,U1(x,y,z),U2(x,y,z)和U3(x,y,z)分别对应空间坐标x,y和z三个方向上的位移分量,Aτm,Bτm和Cτm为梁结构的整体位移中相应项的系数。
(5)计算梁结构的应变向量和应力向量,振动应变具体形式如下:
ε=[εx,εy,εz,γxy,γyz,γxz]T
其中,εx,εy和εz为正应变分量,γxy,γyz和γxz为切应变分量;ε表示梁结构的应变向量,T表示转置。梁材料为各向同性,杨氏模量E=75GPa,泊松比μ=0.3,通过弹性力学中的物理方程计算振动应力具体形式如下:
σ=Dε
其中,σ表示梁结构的应力向量,D为结构材料系数矩阵。
(6)利用梁结构的应变向量和应力向量计算梁结构的应变能和动能,同时,对于两端简支边界,设置虚拟弹簧刚度以获取边界能:kx0=1010,ky0=0,kz0=1010;kxL=1010,kyL=0,kzL=1010。具体表达式如下:
其中,Vs,Tp和Vp分别为梁结构的应变能,动能和边界能方程。t表示时间,ρ为材料的密度。kx0,ky0和kz0为梁结构y=0端所设的虚拟弹簧边界的刚度系数,kxL,kyL和kzL为梁结构y=L端所设的虚拟弹簧边界的刚度系数。
(7)建立结构拉格朗日能量泛函Ω=Vs+Vp-Tp,然后对其中的系数Aτm,Bτm和Cτm求偏导并令其结果为零,即可得到梁结构的3×3阶核心质量矩阵和刚度矩阵。
求得核心矩阵中的元素如下:
其中Kmnij为核心刚度矩阵,Mmnij为核心质量矩阵;上角标i,j=1,…,Γ;m,n=-2,-1,0,…,M;下角标a,b和c用以表示核心矩阵中的各个元素,例如ab表示矩阵中第a行第b列的元素。和表示分别对x和z求一阶偏导,和表示分别对x和z求一阶偏导,同理Ψm,y和Ψn,y分别表示Ψm和Ψn对y求一阶导数。D11,…,D66为结构材料系数矩阵D中的元素。
(8)如图3所示,通过迭代循环核心矩阵求得整体的质量矩阵和刚度矩阵:指针i,j由1取到Γ=6循环核心矩阵Kmnij得到子矩阵Kmn,指针m,n由-2取到M=10循环子矩阵Kmn得到总体刚度矩阵K。通过相同的方法循环核心质量矩阵Mmnij得到总体质量矩阵M,进而得到结构的特征方程:
(K-ω2M)A=0
其中ω为固有频率,A为对应ω的特征向量。
(9)应用Arnoldi算法建立MATLAB求解器求解梁结构的特征方程矩阵输出其固有频率和振型。并进行收敛性判断:给予初始的改进傅里叶截断项数M1,算得结果并进行存储,增加截断项数至M2,通过前后两次结果的对比,判断计算结果是否收敛。
计算结果如表1所示,表1中给出了不同截断项数M下本方法计算出的结构整体的前五阶弯曲频率。细长梁指高跨比h2/L为0.01,短粗梁指高跨比h2/L为0.1,以传统的卡诺高阶梁理论计算结果作为参考。从表中可以看出,本发明方法无论是在处理细长梁还是短粗梁均具有较好的收敛速度,当截断项数为10时,前五阶弯曲模态已基本收敛。另一方面,通过与参考值对比可以发现本研究方法具有很好的计算精度。
表1梁结构前五阶弯曲模态收敛性分析(Hz)
注:参考值源自论文Carrera E,Petrolo M,Nali P.Unified formulationapplied to free vibrations finite element analysis of beams with arbitrarysection[J].Shock&Vibration,2015,18(3):485–502.
Claims (1)
1.一种均匀截面梁结构的振动分析方法,其特征在于,包括以下步骤:
(1)提取梁结构的几何尺寸特性,将梁结构的位移场分离为截面面内位移和轴向位移,具体表达式为:
Uk(x,y,z)=Φk(x,z)Fk(y)
其中x和z为结构空间沿着梁截面方向的坐标,y为梁轴向坐标,Uk(x,y,z)为梁结构的整体位移,k=1,2,3分别表示x,y和z三个方向上的分量,Φk(x,z)为截面面内位移,Fk(y)为轴向位移;
(2)利用卡诺高阶截取技术对梁结构截面面内位移进行拟合,拟合形式如下:
(3)采用改进傅里叶级数对梁结构轴向位移进行延展,具体形式如下:
其中,M为改进傅里叶级数的截断项数,m=-2,-1,...,M;bkm代表相应改进傅里叶展开项的系数,L为梁结构的总长度;
(4)结合梁结构截面面内位移和轴向位移,计算得到梁结构的整体位移,具体表达式如下:
其中,U1(x,y,z),U2(x,y,z)和U3(x,y,z)分别对应空间坐标x,y和z三个方向上的位移分量,Aτm,Bτm和Cτm为位移分量中相应项的系数;
(5)计算梁结构的应变向量和应力向量;
所涉及梁结构的应变向量的表达式为:
ε=[εx,εy,εz,γxy,γyz,γxz]T
其中,ε表示梁结构的应变向量;上标T表示转置;εx,εy和εz为正应变分量;γxy,γyz和γxz为切应变分量,且有
所涉及应力向量的表达式为:
σ=Dε
其中,σ表示梁结构的应力向量,D为结构材料系数矩阵;
(6)计算梁结构的应变能和动能方程,设置虚拟弹簧边界从而获取边界能,具体表达式如下:
其中,Vs,Tp和Vp分别为梁结构的应变能、动能和边界能方程;t表示时间,ρ为材料的密度;kx0,ky0和kz0为梁结构y=0端所设的虚拟弹簧边界的刚度系数,kxL,kyL和kzL为梁结构y=L端所设的虚拟弹簧边界的刚度系数;
(7)建立结构拉格朗日能量泛函Ω=Vs+Vp-Tp,然后对系数Aτm,Bτm和Cτm求偏导并令其结果为零,即可得到梁结构的3×3阶核心质量矩阵和刚度矩阵;核心矩阵中的元素如下:
其中,Kmnij为核心刚度矩阵,Mmnij为核心质量矩阵;上角标i,j=1,…,Γ;m,n=-2,-1,0,…,M;下角标a,b和c表示核心矩阵中元素的标号;和表示分别对x和z求一阶偏导,和表示分别对x和z求一阶偏导,Ψm,y和Ψn,y分别表示Ψm和Ψn对y求一阶导数;D11,…,D66为结构材料系数矩阵D中的元素;
(8)通过迭代循环核心矩阵求得总体质量矩阵M、总体刚度矩阵,进而得到结构的特征方程;
所述质量矩阵和刚度矩阵的求解方法为:指针i,j由1取到Γ循环核心矩阵Kmnij得到子矩阵Kmn,指针m,n由-2取到M循环子矩阵Kmn得到总体刚度矩阵K,通过相同的方法循环核心质量矩阵Mmnij得到总体质量矩阵M;
所述结构的特征方程表达式为:
(K-ω2M)A=0
其中ω为圆频率,A为对应ω的特征向量;
(9)求解梁结构的特征方程矩阵并计算得到固有频率和振型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711016927.4A CN107808048B (zh) | 2017-10-26 | 2017-10-26 | 一种均匀截面梁结构的振动分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711016927.4A CN107808048B (zh) | 2017-10-26 | 2017-10-26 | 一种均匀截面梁结构的振动分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107808048A CN107808048A (zh) | 2018-03-16 |
CN107808048B true CN107808048B (zh) | 2020-12-25 |
Family
ID=61582269
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711016927.4A Active CN107808048B (zh) | 2017-10-26 | 2017-10-26 | 一种均匀截面梁结构的振动分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107808048B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109580193B (zh) * | 2018-11-14 | 2020-06-23 | 中北大学 | 一种在爆炸荷载作用下等截面简支梁荷载系数的计算方法 |
CN111259470B (zh) * | 2020-01-13 | 2022-11-04 | 河海大学 | 任意变截面梁式构件的固有频率检测方法 |
CN111310321B (zh) * | 2020-02-07 | 2023-05-30 | 大连海事大学 | 一种基于Pasternak地基模型的层状土中单桩水平振动分析方法 |
CN111611693B (zh) * | 2020-04-27 | 2024-03-19 | 苏州科技大学 | 一种多段连续梁固有频率的计算方法 |
CN111881606B (zh) * | 2020-07-29 | 2022-09-23 | 河海大学常州校区 | 一种考虑截面变形的伸缩臂臂架动力学建模方法 |
CN112052524B (zh) * | 2020-09-25 | 2022-09-06 | 中国直升机设计研究所 | 一种直升机吊挂柔性机身建模方法 |
CN113051718B (zh) * | 2021-03-07 | 2023-01-03 | 天津大学 | 一种有延展假设分组拓扑径向受载圆环静力学分析方法 |
CN113591233B (zh) * | 2021-06-03 | 2024-04-30 | 南京林业大学 | 一种用于求解柔性梁系统动力学方程的分裂迭代方法 |
CN113486512B (zh) * | 2021-07-05 | 2022-05-10 | 哈尔滨工程大学 | 一种功能梯度变厚度叶片模型的颤振分析方法 |
CN113901544B (zh) * | 2021-09-14 | 2023-05-23 | 浙江工业大学 | 梁结构中波的空间分离与约束控制方法 |
CN116384162B (zh) * | 2023-06-05 | 2023-08-08 | 华东交通大学 | 基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887474A (zh) * | 2010-06-25 | 2010-11-17 | 哈尔滨工程大学 | 基于有限元法与广义傅里叶级数法的结构振动分析方法 |
CN102096736A (zh) * | 2011-01-28 | 2011-06-15 | 重庆大学 | 一种基于渐近变分法的复合材料层合板仿真及优化方法 |
CN103902504A (zh) * | 2014-04-11 | 2014-07-02 | 北京工业大学 | 应用改进微分变换法计算欧拉-伯努利梁固有频率的方法 |
CN105184060A (zh) * | 2015-08-24 | 2015-12-23 | 哈尔滨工程大学 | 一种基于空间配面与改进傅里叶级数的层合结构振动分析方法 |
CN105183958A (zh) * | 2015-08-24 | 2015-12-23 | 哈尔滨工程大学 | 一种复合材料层合结构三维振动分析方法 |
WO2017070187A1 (en) * | 2015-10-19 | 2017-04-27 | Massachusetts Institute Of Technology | Micro electromechanical system (mems) energy harvester with residual stress induced instability |
-
2017
- 2017-10-26 CN CN201711016927.4A patent/CN107808048B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887474A (zh) * | 2010-06-25 | 2010-11-17 | 哈尔滨工程大学 | 基于有限元法与广义傅里叶级数法的结构振动分析方法 |
CN102096736A (zh) * | 2011-01-28 | 2011-06-15 | 重庆大学 | 一种基于渐近变分法的复合材料层合板仿真及优化方法 |
CN103902504A (zh) * | 2014-04-11 | 2014-07-02 | 北京工业大学 | 应用改进微分变换法计算欧拉-伯努利梁固有频率的方法 |
CN105184060A (zh) * | 2015-08-24 | 2015-12-23 | 哈尔滨工程大学 | 一种基于空间配面与改进傅里叶级数的层合结构振动分析方法 |
CN105183958A (zh) * | 2015-08-24 | 2015-12-23 | 哈尔滨工程大学 | 一种复合材料层合结构三维振动分析方法 |
WO2017070187A1 (en) * | 2015-10-19 | 2017-04-27 | Massachusetts Institute Of Technology | Micro electromechanical system (mems) energy harvester with residual stress induced instability |
Non-Patent Citations (6)
Title |
---|
"A series solution for the vibrations of composite laminated deep curved beams with general boundaries";Tiangui Ye et al.;《Composite Structures》;20150314;第450-465页 * |
"A unified approach for the vibration analysis of moderately thick composite laminated cylindrical shells with arbitrary boundary conditions";GuoyongJin et al.;《International Journal of Mechanical Sciences》;20130820;第357-376页 * |
"A unified formulation for vibration analysis of open shells with arbitrary boundary conditions";Tiangui Ye et al.;《International Journal of Mechanical Sciences》;20140208;第42-59页 * |
"Elasticity solution for vibration of generally laminated beams by a modified Fourier expansion-based sampling surface method";Tiangui Ye et al.;《Computers and Structures》;20160302;第115-130页 * |
"任意边界及耦合条件下的多跨梁结构振动特性";郑超凡 等;《中国舰船研究》;20170831;第12卷(第4期);第95-101页 * |
"矩形厚梁结构弯曲自由振动精确化方程新形式";胡超 等;《中国科学》;20141231;第44卷(第4期);第441-448页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107808048A (zh) | 2018-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107808048B (zh) | 一种均匀截面梁结构的振动分析方法 | |
CN107818209B (zh) | 一种弹性板结构的振动分析方法 | |
Dell'Isola et al. | Elastic pantographic 2D lattices: a numerical analysis on the static response and wave propagation | |
Yu | An introduction to micromechanics | |
Lee et al. | A refined higher order shear and normal deformation theory for E-, P-, and S-FGM plates on Pasternak elastic foundation | |
CN105975645B (zh) | 一种基于多步的含激波区域飞行器流场快速计算方法 | |
Damanpack et al. | High-order free vibration analysis of sandwich beams with a flexible core using dynamic stiffness method | |
Han et al. | Accuracy of the adaptive GRP scheme and the simulation of 2-D Riemann problems for compressible Euler equations | |
Rezaiee-Pajand et al. | An efficient curved beam element for thermo-mechanical nonlinear analysis of functionally graded porous beams | |
CN103942381B (zh) | 用于飞机铝合金结构性能预测的状态近场动力学方法 | |
Mahran et al. | A comparison between different finite elements for elastic and aero-elastic analyses | |
CN105184060B (zh) | 一种基于空间配面与改进傅里叶级数的层合结构振动分析方法 | |
Gupta et al. | Extensions of the two-scale generalized finite element method to nonlinear fracture problems | |
D’Amore et al. | A Parallel Three‐dimensional Variational Data Assimilation Scheme | |
Rango et al. | Static and dynamic analysis of thick laminated plates using enriched macroelements | |
Liu et al. | A smoothed finite element method (S-FEM) for large-deformation elastoplastic analysis | |
Devarajan | Vibration Analysis of Timoshenko Beams using Isogeometric Analysis | |
CN107748821B (zh) | 一种三维耦合结构的振动分析方法 | |
Pagani et al. | Flutter analysis by refined 1D dynamic stiffness elements and doublet lattice method | |
Yeliseyev et al. | Two-dimensional (shell-type) and three-dimensional models for elastic thin-walled cylinder | |
Wu | Study on the self-propulsion of the rigid-flexible composite plate | |
Kianian et al. | Buckling analysis of functionally graded plates based on two-variable refined plate theory using the bubble finite strip method | |
Pachenari et al. | Analysis of tapered thin plates using basic displacement functions | |
Blanc et al. | Analysis and compression of time-accurate turbomachinery simulations using proper orthogonal decomposition | |
Das et al. | Nonlinear vibration analysis of clamped skew plates by a variational method |
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 |