CN115310188A - 一种基于最小势能理论确定加筋板等效参数的方法 - Google Patents
一种基于最小势能理论确定加筋板等效参数的方法 Download PDFInfo
- Publication number
- CN115310188A CN115310188A CN202211243855.8A CN202211243855A CN115310188A CN 115310188 A CN115310188 A CN 115310188A CN 202211243855 A CN202211243855 A CN 202211243855A CN 115310188 A CN115310188 A CN 115310188A
- Authority
- CN
- China
- Prior art keywords
- represent
- plate
- stiffened
- displacement
- stiffened plate
- 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
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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Structural Engineering (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Laminated Bodies (AREA)
Abstract
本发明涉及一种基于最小势能理论确定加筋板等效参数的方法,包括:步骤(1):计算加筋板的应变能;步骤(2):计算加筋板的惯性功;步骤(3):根据所述加筋板的应变能和惯性功,基于最小势能理论得到加筋板的平衡方程;步骤(4):获取正交各向异性板的平衡方程,并将所述正交各向异性板的平衡方程和所述加筋板的平衡方程进行等效,得到加筋板的平衡参数。本发明建立了一个加强筋和正交各向异性板参数转化的桥梁,方便分析加强筋的振动方程。
Description
技术领域
本发明涉及加筋板结构计算技术领域,特别是涉及一种基于最小势能理论确定加筋板等效参数的方法。
背景技术
加筋板是指薄板和各种骨材的组合结构,骨材的存在能很好的影响和改变薄板结构的声振性能。并且骨材对平板的约束作用将薄板结构分割成一系列的小板格区域,这样骨材成为薄板抵抗失稳的主要承力部分,增强了薄板的力学性能,约束了薄板的振动。此外,骨材的存在还提高了加筋板结构的整体辐射阻抗,从而降低了结构的振速和辐射声功率水平。
关于加筋板结构的振动和声辐射研究,大体上主要分为解析法、数值法和实验法。解析法作为传统的研究方法,主要通过公式,考虑骨材对板材的作用,进行振动控制方程的求解,然后根据声振耦合连续性条件求解声场。
解析法针对的对象主要是两类,一类是无限大结构,一类是有限大结构。对于无限大结构主要的思路是通过傅里叶变化,将空间域转换到波数域,然后考虑周期函数的性质,运用无穷级数展开等方法处理振动控制方程进行求解,主要有波数法、无穷级数法;对于有限大结构主要研究方法有能量法、拉格朗日法、海姆霍子法等。
目前暂无有效的方法能针对针对加筋板和正交各向异性板的振动方程参数进行等效求解。
发明内容
本发明所要解决的技术问题是提供一种基于最小势能理论确定加筋板等效参数的方法,建立了一个加强筋和正交各向异性板参数转化的桥梁,方便分析加强筋的振动方程。
本发明解决其技术问题所采用的技术方案是:提供一种基于最小势能理论确定加筋板等效参数的方法,包括:
步骤(1):计算加筋板的应变能;
步骤(2):计算加筋板的惯性功;
步骤(3):根据所述加筋板的应变能和惯性功,基于最小势能理论得到加筋板的平衡方程;
步骤(4):获取正交各向异性板的平衡方程,并将所述正交各向异性板的平衡方程和所述加筋板的平衡方程进行等效,得到加筋板的平衡参数。
所述步骤(1)包括:
步骤(11):计算无加强筋的各向同性平板应变能;
所述步骤(11)中无加强筋的各向同性平板应变能的公式为:,其中,表示杨氏模量,表示泊松比,表示平板的长度,表示平板的宽度,表示,为方向位移,表示,为方向位移,表示平板的弯曲刚度且,表示平板的厚度,表示,表示,表示,为方向位移。
所述步骤(12)中N个平行轴加强筋的总应变能的公式为:,其中,为轴方向上加强筋之间的间距,表示平板的长度,表示平板的宽度,表示方向加强筋的弹性模量,表示方向梁截面的横截面积,表示,为方向位移,表示方向加强筋形心在板局部坐标系的坐标,表示,,表示方向加强筋截面惯性矩,表示方向梁截面的扭转刚度,表示,为方向位移。
所述步骤(12)中N个平行轴加强筋的总应变能的公式为:,其中,为轴方向上加强筋之间的间距,表示平板的长度,表示平板的宽度,表示方向加强筋的弹性模量,表示方向梁截面的横截面积,表示,表示方向位移,表示方向加强筋形心在板局部坐标系的坐标,表示,,表示方向加强筋截面惯性矩,表示方向梁截面的扭转刚度,表示,为方向位移。
所述步骤(2)中加筋板的惯性功公式为:,其中,表示平板的长度,表示平板的宽度,表示角频率,表示方向位移,表示将加强筋抹平之后所得等效面密度且,和分别表示圆柱体、梁和圆环的质量密度,表示平板的厚度,表示方向梁截面的面积,表示轴方向上加强筋之间的间距,表示方向梁截面的横截面积,表示轴方向上加强筋之间的间距。
所述步骤(3)中加筋板的平衡方程公式为:
所述步骤(4)中的正交各向异性板的平衡方程公式为:
其中,表示方向应力条件下的杨氏模量,表示,表示沿方向位移,表示方向的泊松比,表示,表示,为方向位移,表示板材的密度,表示,表示时间,表示平面的剪切模量,表示,表示,表示方向应力条件下的杨氏模量,表示,表示方向的泊松比,表示,表示方向的弯曲刚度,表示,表示方向的弯曲刚度,,表示平板的厚度,表示,表示。
所述步骤(4)具体为:
步骤(41):获取正交各向异性板的平衡方程;
步骤(42):在纯弯曲振动下将加筋板的平衡方程与正交各向异性板的平衡方程进行等效,得到加筋板的平衡参数。
有益效果
由于采用了上述的技术方案,本发明与现有技术相比,具有以下的优点和积极效果:本发明建立一个加强筋和正交各向异性板参数转化的桥梁,方便分析加强筋的振动方程;本发明由平板的能量原理导出了动态平衡方程和边界条件,并将加筋板平衡方程和正交各向异性板平衡方程的对比进行等效,进而能够让加筋板用正交各向异性板来等效。
附图说明
图1是本发明实施方式的方法流程图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
本发明的实施方式涉及一种基于最小势能理论确定加筋板等效参数的方法,请参阅图1,包括:
步骤(1):计算加筋板的应变能;
步骤(2):计算加筋板的惯性功;
步骤(3):根据所述加筋板的应变能和惯性功,基于最小势能理论得到加筋板的平衡方程;
步骤(4):获取正交各向异性板的平衡方程,并将所述正交各向异性板的平衡方程和所述加筋板的平衡方程进行等效,得到加筋板的平衡参数。
以下对本实施方式进行详细介绍:
1)加强筋平板应变能的求解
无加强筋的各向同性平板应变能如下:
根据线性Donnell-type应变位移关系:
加筋板中截面应变定义如下:
将应变与位移关系带入不考虑加强筋的平板可得应变能为:
在近似计算平板和加强筋在附着点上的位移相等的基础上,推导出了加强筋的应变能如下:
其中,为轴方向上加强筋之间的间距,表示平板的长度,表示平板的厚度,表示方向加筋的弹性模量,为方向梁截面的横截面积,表示,为方向加强筋形心在板局部坐标系的坐标,表示,,为方向加强筋截面惯性矩,为方向梁截面的扭转刚度,表示。
2)惯性载荷的势能
3)平衡方程和边界条件
系统总势能为加强筋应变能、均质板应变能、惯性功之和:
4)加筋板的等效参数
根据资料可得正交各向异性板的平衡方程如下:
其中,表示方向应力条件下的杨氏模量,表示,表示沿方向位移,表示方向的泊松比,表示,表示,为方向位移,表示板材的密度,表示,表示时间,表示平面的剪切模量,表示,表示,表示方向应力条件下的杨氏模量,表示,表示方向的泊松比,表示,表示方向的弯曲刚度,表示,表示方向的弯曲刚度,,表示平板的厚度,表示,表示。
将正交各向异性板的平衡方程:
和上述求得的加筋板平衡方程:
进行等效,对比上述两式,可得:
可进行等效:
由此可见,本发明提供了一种求解加筋板和正交各向异性板振动方程等效参数的方法,由平板的能量原理导出了动态平衡方程和边界条件,并将加筋板平衡方程和正交各向异性板平衡方程的对比进行等效,进而能够让加筋板用正交各向异性板来等效。
Claims (9)
1.一种基于最小势能理论确定加筋板等效参数的方法,其特征在于,包括:
步骤(1):计算加筋板的应变能;
步骤(2):计算加筋板的惯性功;
步骤(3):根据所述加筋板的应变能和惯性功,基于最小势能理论得到加筋板的平衡方程;
步骤(4):获取正交各向异性板的平衡方程,并将所述正交各向异性板的平衡方程和所述加筋板的平衡方程进行等效,得到加筋板的平衡参数。
9.根据权利要求1所述的基于最小势能理论确定加筋板等效参数的方法,其特征在于,所述步骤(4)具体为:
步骤(41):获取正交各向异性板的平衡方程;
步骤(42):在纯弯曲振动下将加筋板的平衡方程与正交各向异性板的平衡方程进行等效,得到加筋板的平衡参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211243855.8A CN115310188B (zh) | 2022-10-12 | 2022-10-12 | 一种基于最小势能理论确定加筋板等效参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211243855.8A CN115310188B (zh) | 2022-10-12 | 2022-10-12 | 一种基于最小势能理论确定加筋板等效参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115310188A true CN115310188A (zh) | 2022-11-08 |
CN115310188B CN115310188B (zh) | 2023-01-10 |
Family
ID=83868084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211243855.8A Active CN115310188B (zh) | 2022-10-12 | 2022-10-12 | 一种基于最小势能理论确定加筋板等效参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115310188B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115618163A (zh) * | 2022-11-11 | 2023-01-17 | 上海索辰信息科技股份有限公司 | 一种加筋结构体固有频率计算方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120245862A1 (en) * | 2009-09-14 | 2012-09-27 | Airbus Operations (S.A.S.) | Method for the structural analysis of panels consisting of an isotropic material and stiffened by triangular pockets |
CN105844033A (zh) * | 2016-03-30 | 2016-08-10 | 中国舰船研究设计中心 | 船舶总纵极限强度简化逐步破坏分析方法 |
CN108256214A (zh) * | 2018-01-16 | 2018-07-06 | 滨州学院 | 一种航空正交加筋板刚度计算方法及装置 |
CN110852012A (zh) * | 2019-11-11 | 2020-02-28 | 武汉理工大学 | 轴压加筋板极限强度预测方法 |
CN113722819A (zh) * | 2021-08-12 | 2021-11-30 | 中国舰船研究设计中心 | 一种计算加筋板弯曲变形与应力的半解析方法 |
CN114912302A (zh) * | 2022-07-18 | 2022-08-16 | 上海索辰信息科技股份有限公司 | 一种基于模态空间采样算法的加筋板模态密度获取方法 |
CN115081148A (zh) * | 2022-07-20 | 2022-09-20 | 上海索辰信息科技股份有限公司 | 一种基于势能理论的加筋板等效参数确定方法 |
-
2022
- 2022-10-12 CN CN202211243855.8A patent/CN115310188B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120245862A1 (en) * | 2009-09-14 | 2012-09-27 | Airbus Operations (S.A.S.) | Method for the structural analysis of panels consisting of an isotropic material and stiffened by triangular pockets |
CN105844033A (zh) * | 2016-03-30 | 2016-08-10 | 中国舰船研究设计中心 | 船舶总纵极限强度简化逐步破坏分析方法 |
CN108256214A (zh) * | 2018-01-16 | 2018-07-06 | 滨州学院 | 一种航空正交加筋板刚度计算方法及装置 |
CN110852012A (zh) * | 2019-11-11 | 2020-02-28 | 武汉理工大学 | 轴压加筋板极限强度预测方法 |
CN113722819A (zh) * | 2021-08-12 | 2021-11-30 | 中国舰船研究设计中心 | 一种计算加筋板弯曲变形与应力的半解析方法 |
CN114912302A (zh) * | 2022-07-18 | 2022-08-16 | 上海索辰信息科技股份有限公司 | 一种基于模态空间采样算法的加筋板模态密度获取方法 |
CN115081148A (zh) * | 2022-07-20 | 2022-09-20 | 上海索辰信息科技股份有限公司 | 一种基于势能理论的加筋板等效参数确定方法 |
Non-Patent Citations (1)
Title |
---|
王小明: "求解加筋板的位移与应力的一种解析计算方法", 《武汉理工大学学报(交通科学与工程版)》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115618163A (zh) * | 2022-11-11 | 2023-01-17 | 上海索辰信息科技股份有限公司 | 一种加筋结构体固有频率计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115310188B (zh) | 2023-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Thai et al. | A simple refined theory for bending, buckling, and vibration of thick plates resting on elastic foundation | |
Kolahchi et al. | Visco-nonlocal-refined Zigzag theories for dynamic buckling of laminated nanoplates using differential cubature-Bolotin methods | |
Banerjee et al. | Dynamic stiffness method for inplane free vibration of rotating beams including Coriolis effects | |
Banerjee | Dynamic stiffness formulation and free vibration analysis of centrifugally stiffened Timoshenko beams | |
Hosseini-Hashemi et al. | A new exact analytical approach for free vibration of Reissner–Mindlin functionally graded rectangular plates | |
Khalili et al. | Free vibration analysis of homogeneous isotropic circular cylindrical shells based on a new three-dimensional refined higher-order theory | |
Wang et al. | Weak form quadrature element method and its applications in science and engineering: a state-of-the-art review | |
Liu et al. | The singular edge-based smoothed finite element method for stationary dynamic crack problems in 2D elastic solids | |
Nguyen-Thoi et al. | Free and forced vibration analysis using the n-sided polygonal cell-based smoothed finite element method (nCS-FEM) | |
Jung et al. | Four-variable refined plate theory for forced-vibration analysis of sigmoid functionally graded plates on elastic foundation | |
Aydogdu et al. | Vibration analysis of cross-ply laminated square plates with general boundary conditions | |
Zhang et al. | Extension of Hencky bar-net model for vibration analysis of rectangular plates with rectangular cutouts | |
Pan et al. | Modeling large amplitude vibration of pretwisted hybrid composite blades containing CNTRC layers and matrix cracked FRC layers | |
Senjanović et al. | Modified Mindlin plate theory and shear locking-free finite element formulation | |
Yu et al. | A stabilized discrete shear gap extended finite element for the analysis of cracked Reissner–Mindlin plate vibration problems involving distorted meshes | |
Thinh et al. | Vibration and buckling analysis of functionally graded plates using new eight-unknown higher order shear deformation theory | |
Azrar et al. | Non-linear forced vibrations of plates by an asymptotic–numerical method | |
Zghal et al. | Vibration characteristics of plates and shells with functionally graded pores imperfections using an enhanced finite shell element | |
Cho et al. | Free vibration analysis of stiffened panels with lumped mass and stiffness attachments | |
Hu et al. | Dynamic and static isogeometric analysis for laminated Timoshenko curved microbeams | |
CN115310188B (zh) | 一种基于最小势能理论确定加筋板等效参数的方法 | |
Cheung et al. | Vibration analysis of symmetrically laminated rectangular plates with intermediate line supports | |
Xing et al. | Characteristic equations and closed-form solutions for free vibrations of rectangular Mindlin plates | |
Kumar et al. | Exact modal analysis of multilayered FG-CNT plate assemblies using the dynamic stiffness method | |
Wang et al. | Free vibration of soft-core sandwich panels with general boundary conditions by harmonic quadrature element 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 |