CN115081148A - 一种基于势能理论的加筋板等效参数确定方法 - Google Patents
一种基于势能理论的加筋板等效参数确定方法 Download PDFInfo
- Publication number
- CN115081148A CN115081148A CN202210850620.9A CN202210850620A CN115081148A CN 115081148 A CN115081148 A CN 115081148A CN 202210850620 A CN202210850620 A CN 202210850620A CN 115081148 A CN115081148 A CN 115081148A
- Authority
- CN
- China
- Prior art keywords
- stiffened plate
- plate
- stiffened
- potential energy
- longitudinal
- 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
- 238000005381 potential energy Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 14
- 238000005452 bending Methods 0.000 claims abstract description 25
- 238000006073 displacement reaction Methods 0.000 claims abstract description 18
- 230000003014 reinforcing effect Effects 0.000 claims description 20
- 239000011324 bead Substances 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 9
- 239000003351 stiffener Substances 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 5
- 239000002699 waste material Substances 0.000 abstract description 2
- 230000002457 bidirectional effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
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/17—Mechanical parametric or variational design
-
- 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
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/10—Noise analysis or noise optimisation
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种基于势能理论的加筋板等效参数确定方法,包括以下步骤:基于势能理论构建加筋板弯曲振动最大势能公式和加筋板动能公式;将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式,得到平衡方程,并基于所述平衡方程计算加筋板的共振频率;令所述加筋板的共振频率与正交各向异性板的共振频率相等,得到所述加筋板等效为正交各向异性板的各个等效参数。本发明整个过程无需通过试验测试等手段,在减少了资源浪费的情况下能够保证结果准确度和计算效率。
Description
技术领域
本发明涉及振动噪声仿真技术领域,特别是涉及一种基于势能理论的加筋板等效参数确定方法。
背景技术
加筋板是航空航天、造船、自动化等工业领域常用的结构形式。在求解加筋板结构的振动问题时,加筋板的固有频率是一个重要的振动参数。目前针对加筋板的固有频率并没有普适通用的理论解析解,不同的专家学者通过试验方法、有限元数值计算方法进行计算,需要耗费巨大的人力财力和计算时间。
发明内容
本发明所要解决的技术问题是提供一种基于势能理论的加筋板等效参数确定方法,能够提高计算效率。
本发明解决其技术问题所采用的技术方案是:提供一种基于势能理论的加筋板等效参数确定方法,包括以下步骤:
基于势能理论构建加筋板弯曲振动最大势能公式和加筋板动能公式;
将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式,得到平衡方程,并基于所述平衡方程计算加筋板的共振频率;
令所述加筋板的共振频率与正交各向异性板的共振频率相等,得到所述加筋板等效为正交各向异性板的各个等效参数。
所述加筋板弯曲振动最大势能公式为:
其中,为加筋板弯曲振动最大势能,表示平板的弯曲刚度,表示加筋板纵向的长度,表示加筋板横向的长度,分别表示加筋板纵向和横向坐标,表示原频率,表示加筋板材料的泊松比,表示加强筋材料的杨氏模量,表示加强筋材料的剪切模量,表示纵向加强筋的数量,表示横向加强筋的数量,分别表示纵向加强筋和横向加强筋的弯曲惯性矩,分别表示纵向加强筋和横向加强筋的扭转常数。
所述加筋板动能公式为:,其中,为加筋板动能,为加筋板的共振频率,分别为加筋板纵向和横向的模态阶数,为加筋板的面密度,表示加筋板纵向的长度,表示加筋板横向的长度,表示位移,分别表示加筋板纵向和横向坐标,表示纵向加强筋的数量,表示横向加强筋的数量,表示纵向加强筋的线密度,表示横向加强筋的线密度。
所述将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式前,需要将所述加筋板Z向位移根据边界条件展开为傅里叶级数形式。
有益效果
由于采用了上述的技术方案,本发明与现有技术相比,具有以下的优点和积极效果:本发明基于势能理论构建加筋板的弯曲振动最大势能公式和动能公式,再结合加筋板的Z向位移构建平衡方程得到加筋板的共振频率,再结合正交各向异性板的共振频率得到加筋板等效为正交各向异性平板的等效参数,整个过程无需通过试验测试等手段,在减少了资源浪费的情况下能够保证结果准确度和计算效率。
附图说明
图1是本发明实施方式的流程图;
图2是双向加筋平板结构示意图;
图3是纵向加筋平板结构示意图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
加筋板在较低频率时可等效为正交各向异性平板,平衡方程如下所示:
本发明的实施方式涉及一种基于势能理论的加筋板等效参数确定方法,如图1所示,包括以下步骤:基于势能理论构建加筋板弯曲振动最大势能公式和加筋板动能公式;将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式,得到平衡方程,并基于所述平衡方程计算加筋板的共振频率;令所述加筋板的共振频率与正交各向异性板的共振频率相等,得到所述加筋板等效为正交各向异性板的各个等效参数。
下面以双向加筋平板结构为例进一步说明本发明,其中,双向加筋平板结构如图2所示。
对于双向加筋的加筋板,根据势能理论,加筋板弯曲振动最大势能公式为:
其中,为加筋板弯曲振动最大势能,表示平板的弯曲刚度,表示加筋板纵向的长度,表示加筋板横向的长度,分别表示加筋板纵向和横向坐标,表示原频率,表示加筋板材料的泊松比,表示加强筋材料的杨氏模量,表示加强筋材料的剪切模量,表示纵向加强筋的数量,表示横向加强筋的数量,分别表示纵向加强筋和横向加强筋的弯曲惯性矩,分别表示纵向加强筋和横向加强筋的扭转常数。上式中,等式右边三项中的第一项为无加筋平板的势能,第二项和第三项表示两个方向加强筋的势能。
将位移带入到加筋板弯曲振动最大势能公式和动能公式,并令二者相等得到平衡方程,通过求解平衡方程得到加筋板的共振频率:
对于单向加筋的加筋板,如图3所示,以只有平行X轴的加强筋为例进行说明。
不难发现,本发明基于势能理论构建加筋板的弯曲振动最大势能公式和动能公式,再结合加筋板的Z向位移构建平衡方程得到加筋板的共振频率,再结合正交各向异性板的共振频率得到加筋板等效为正交各向异性平板的等效参数,整个过程无需通过试验测试等手段,在减少了资源浪费的情况下能够保证结果准确度和计算效率。
Claims (4)
1.一种基于势能理论的加筋板等效参数确定方法,其特征在于,包括以下步骤:
基于势能理论构建加筋板弯曲振动最大势能公式和加筋板动能公式;
将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式,得到平衡方程,并基于所述平衡方程计算加筋板的共振频率;
令所述加筋板的共振频率与正交各向异性板的共振频率相等,得到所述加筋板等效为正交各向异性板的各个等效参数。
4.根据权利要求1所述的基于势能理论的加筋板等效参数确定方法,其特征在于,所述将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式前,需要将所述加筋板Z向位移根据边界条件展开为傅里叶级数形式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210850620.9A CN115081148B (zh) | 2022-07-20 | 2022-07-20 | 一种基于势能理论的加筋板等效参数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210850620.9A CN115081148B (zh) | 2022-07-20 | 2022-07-20 | 一种基于势能理论的加筋板等效参数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115081148A true CN115081148A (zh) | 2022-09-20 |
CN115081148B CN115081148B (zh) | 2022-11-15 |
Family
ID=83260464
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210850620.9A Active CN115081148B (zh) | 2022-07-20 | 2022-07-20 | 一种基于势能理论的加筋板等效参数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115081148B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115310188A (zh) * | 2022-10-12 | 2022-11-08 | 上海索辰信息科技股份有限公司 | 一种基于最小势能理论确定加筋板等效参数的方法 |
CN115618163A (zh) * | 2022-11-11 | 2023-01-17 | 上海索辰信息科技股份有限公司 | 一种加筋结构体固有频率计算方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013101598A (ja) * | 2011-10-17 | 2013-05-23 | Mizuho Information & Research Institute Inc | 挙動解析システム、挙動解析方法及び挙動解析プログラム |
US20140379314A1 (en) * | 2011-10-31 | 2014-12-25 | Osaka University | Analyzer, analysis method, and analysis program |
CN108256214A (zh) * | 2018-01-16 | 2018-07-06 | 滨州学院 | 一种航空正交加筋板刚度计算方法及装置 |
CN108491591A (zh) * | 2018-03-06 | 2018-09-04 | 东南大学 | 一种高温环境下曲线加筋板有限元分析方法 |
CN109145369A (zh) * | 2018-07-11 | 2019-01-04 | 东南大学 | 一种计及非共振传输的中高频局部动响应预示方法 |
CN113722819A (zh) * | 2021-08-12 | 2021-11-30 | 中国舰船研究设计中心 | 一种计算加筋板弯曲变形与应力的半解析方法 |
CN114239153A (zh) * | 2021-12-31 | 2022-03-25 | 中国航空工业集团公司西安飞机设计研究所 | 一种单向轴压矩形板后屈曲平衡路径计算方法 |
US20220127802A1 (en) * | 2020-10-28 | 2022-04-28 | South China University Of Technology | Combined Plate-Beam Unit Analysis Method Considering Residual Stress Effect of Orthotropic Plate |
-
2022
- 2022-07-20 CN CN202210850620.9A patent/CN115081148B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013101598A (ja) * | 2011-10-17 | 2013-05-23 | Mizuho Information & Research Institute Inc | 挙動解析システム、挙動解析方法及び挙動解析プログラム |
US20140379314A1 (en) * | 2011-10-31 | 2014-12-25 | Osaka University | Analyzer, analysis method, and analysis program |
CN108256214A (zh) * | 2018-01-16 | 2018-07-06 | 滨州学院 | 一种航空正交加筋板刚度计算方法及装置 |
CN108491591A (zh) * | 2018-03-06 | 2018-09-04 | 东南大学 | 一种高温环境下曲线加筋板有限元分析方法 |
CN109145369A (zh) * | 2018-07-11 | 2019-01-04 | 东南大学 | 一种计及非共振传输的中高频局部动响应预示方法 |
US20220127802A1 (en) * | 2020-10-28 | 2022-04-28 | South China University Of Technology | Combined Plate-Beam Unit Analysis Method Considering Residual Stress Effect of Orthotropic Plate |
CN113722819A (zh) * | 2021-08-12 | 2021-11-30 | 中国舰船研究设计中心 | 一种计算加筋板弯曲变形与应力的半解析方法 |
CN114239153A (zh) * | 2021-12-31 | 2022-03-25 | 中国航空工业集团公司西安飞机设计研究所 | 一种单向轴压矩形板后屈曲平衡路径计算方法 |
Non-Patent Citations (4)
Title |
---|
于士甲等: "基于FEM-MODENA的加筋板声-固耦合分析", 《振动工程学报》 * |
李凯;何书韬;邱永康;吴国民;郭文杰: "附加多个集中质量加筋板的自由振动分析", 《中国舰船研究》 * |
李国荣等: "典型边界条件下加筋矩形板的横向振动特性分析", 《振动与冲击》 * |
杨健生; 曾治平; 韦冬炎; 彭林欣: "基于无网格法的非均匀弹性地基上变厚度加筋板弯曲与固有频率分析", 《计算力学学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115310188A (zh) * | 2022-10-12 | 2022-11-08 | 上海索辰信息科技股份有限公司 | 一种基于最小势能理论确定加筋板等效参数的方法 |
CN115310188B (zh) * | 2022-10-12 | 2023-01-10 | 上海索辰信息科技股份有限公司 | 一种基于最小势能理论确定加筋板等效参数的方法 |
CN115618163A (zh) * | 2022-11-11 | 2023-01-17 | 上海索辰信息科技股份有限公司 | 一种加筋结构体固有频率计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115081148B (zh) | 2022-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115081148B (zh) | 一种基于势能理论的加筋板等效参数确定方法 | |
Mohammadi et al. | An equivalent model for trapezoidal corrugated cores based on homogenization method | |
Guo et al. | Nonlinear dynamics of fluid-conveying composite pipes subjected to time-varying axial tension in sub-and super-critical regimes | |
Ji et al. | On the importance of work-conjugacy and objective stress rates in finite deformation incremental finite element analysis | |
Tian et al. | Analysis of nonlinear aeroelastic characteristics of a trapezoidal wing in hypersonic flow | |
US20070075450A1 (en) | Parametrized material and performance properties based on virtual testing | |
Zhuang et al. | Modal and aeroelastic analysis of trapezoidal corrugated-core sandwich panels in supersonic flow | |
CN103745066B (zh) | 一种大展弦比机翼结构刚度指标的确定方法 | |
CN104316388A (zh) | 一种对各向异性材料结构件进行疲劳寿命测定的方法 | |
Jafari-Talookolaei et al. | In-plane and out-of-plane vibration modes of laminated composite beams with arbitrary lay-ups | |
CN107368634A (zh) | 一种复杂结构耦合损耗因子的预示方法 | |
CN101887474A (zh) | 基于有限元法与广义傅里叶级数法的结构振动分析方法 | |
Su et al. | Modeling and parametric analysis of in-plane free vibration of a floating cable-stayed bridge with transfer matrix method | |
CN105956368A (zh) | 一种复合材料层合结构的非概率动力可靠性评估方法 | |
Liu et al. | Random vibration study of functionally graded porous curved beams with elastically restrained ends | |
CN105787167A (zh) | 聚合物基复合材料有效应力松弛系数的预测方法 | |
CN110852012A (zh) | 轴压加筋板极限强度预测方法 | |
Li et al. | On the non-linear dynamics of a forced plate with boundary conditions correction in subsonic flow | |
CN107268820A (zh) | 无屈曲波形结构耗能构件及其设计方法 | |
CN104794332B (zh) | 一种高层建筑风致响应分析模型的不确定性分析方法 | |
Yan et al. | Structure/material concurrent optimization of lattice materials based on extended multiscale finite element method | |
Chwalowski et al. | Progress on Transonic Flutter and Shock Buffet Computations in Support of the Third Aeroelastic Prediction Workshop | |
Yi | Finite element analysis of thermoelastic damping in contour-mode vibrations of micro-and nanoscale ring, disk, and elliptical plate resonators | |
Zhang et al. | Aeroelastic instability of an inverted cantilevered plate with cracks in axial subsonic airflow | |
Nallim et al. | An analytical–numerical approach to simulate the dynamic behaviour of arbitrarily laminated composite plates |
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 |