CN115081148B - 一种基于势能理论的加筋板等效参数确定方法 - Google Patents

一种基于势能理论的加筋板等效参数确定方法 Download PDF

Info

Publication number
CN115081148B
CN115081148B CN202210850620.9A CN202210850620A CN115081148B CN 115081148 B CN115081148 B CN 115081148B CN 202210850620 A CN202210850620 A CN 202210850620A CN 115081148 B CN115081148 B CN 115081148B
Authority
CN
China
Prior art keywords
stiffened plate
plate
stiffened
potential energy
formula
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
Application number
CN202210850620.9A
Other languages
English (en)
Other versions
CN115081148A (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.)
Shanghai Suochen Information Technology Co ltd
Original Assignee
Shanghai Suochen Information Technology Co ltd
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 Shanghai Suochen Information Technology Co ltd filed Critical Shanghai Suochen Information Technology Co ltd
Priority to CN202210850620.9A priority Critical patent/CN115081148B/zh
Publication of CN115081148A publication Critical patent/CN115081148A/zh
Application granted granted Critical
Publication of CN115081148B publication Critical patent/CN115081148B/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/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • 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
    • 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
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/10Noise analysis or noise optimisation
    • 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)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Operations Research (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种基于势能理论的加筋板等效参数确定方法,包括以下步骤:基于势能理论构建加筋板弯曲振动最大势能公式和加筋板动能公式;将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式,得到平衡方程,并基于所述平衡方程计算加筋板的共振频率;令所述加筋板的共振频率与正交各向异性板的共振频率相等,得到所述加筋板等效为正交各向异性板的各个等效参数。本发明整个过程无需通过试验测试等手段,在减少了资源浪费的情况下能够保证结果准确度和计算效率。

Description

一种基于势能理论的加筋板等效参数确定方法
技术领域
本发明涉及振动噪声仿真技术领域,特别是涉及一种基于势能理论的加筋板等效参数确定方法。
背景技术
加筋板是航空航天、造船、自动化等工业领域常用的结构形式。在求解加筋板结构的振动问题时,加筋板的固有频率是一个重要的振动参数。目前针对加筋板的固有频率并没有普适通用的理论解析解,不同的专家学者通过试验方法、有限元数值计算方法进行计算,需要耗费巨大的人力财力和计算时间。
发明内容
本发明所要解决的技术问题是提供一种基于势能理论的加筋板等效参数确定方法,能够提高计算效率。
本发明解决其技术问题所采用的技术方案是:提供一种基于势能理论的加筋板等效参数确定方法,包括以下步骤:
基于势能理论构建加筋板弯曲振动最大势能公式和加筋板动能公式;
将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式,得到平衡方程,并基于所述平衡方程计算加筋板的共振频率;
令所述加筋板的共振频率与正交各向异性板的共振频率相等,得到所述加筋板等效为正交各向异性板的各个等效参数。
所述加筋板弯曲振动最大势能公式为:
Figure 712216DEST_PATH_IMAGE001
其中,
Figure 503455DEST_PATH_IMAGE002
为加筋板弯曲振动最大势能,
Figure 644586DEST_PATH_IMAGE003
表示平板的弯曲刚度,
Figure 204880DEST_PATH_IMAGE004
表示加筋板纵向的长度,
Figure 337047DEST_PATH_IMAGE005
表示加筋板横向的长度,
Figure 615581DEST_PATH_IMAGE006
分别表示加筋板纵向和横向坐标,
Figure 560403DEST_PATH_IMAGE007
表示原频率,
Figure 476669DEST_PATH_IMAGE008
表示加筋板材料的泊松比,
Figure 481534DEST_PATH_IMAGE009
表示加强筋材料的杨氏模量,
Figure 981786DEST_PATH_IMAGE010
表示加强筋材料的剪切模量,
Figure 464720DEST_PATH_IMAGE011
表示纵向加强筋的数量,
Figure 501071DEST_PATH_IMAGE012
表示横向加强筋的数量,
Figure 473575DEST_PATH_IMAGE013
分别表示纵向加强筋和横向加强筋的弯曲惯性矩,
Figure 461123DEST_PATH_IMAGE014
分别表示纵向加强筋和横向加强筋的扭转常数。
所述加筋板动能公式为:
Figure 249212DEST_PATH_IMAGE015
,其中,
Figure 638605DEST_PATH_IMAGE016
为加筋板动能,
Figure 47590DEST_PATH_IMAGE017
为加筋板的共振频率,
Figure 23898DEST_PATH_IMAGE018
分别为加筋板纵向和横向的模态阶数,
Figure 848635DEST_PATH_IMAGE019
为加筋板的面密度,
Figure 92534DEST_PATH_IMAGE004
表示加筋板纵向的长度,
Figure 610103DEST_PATH_IMAGE005
表示加筋板横向的长度,
Figure 73708DEST_PATH_IMAGE020
表示位移,
Figure 764452DEST_PATH_IMAGE006
分别表示加筋板纵向和横向坐标,
Figure 862858DEST_PATH_IMAGE011
表示纵向加强筋的数量,
Figure 52794DEST_PATH_IMAGE012
表示横向加强筋的数量,
Figure 502230DEST_PATH_IMAGE021
表示纵向加强筋的线密度,
Figure 668769DEST_PATH_IMAGE022
表示横向加强筋的线密度。
所述将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式前,需要将所述加筋板Z向位移根据边界条件展开为傅里叶级数形式。
有益效果
由于采用了上述的技术方案,本发明与现有技术相比,具有以下的优点和积极效果:本发明基于势能理论构建加筋板的弯曲振动最大势能公式和动能公式,再结合加筋板的Z向位移构建平衡方程得到加筋板的共振频率,再结合正交各向异性板的共振频率得到加筋板等效为正交各向异性平板的等效参数,整个过程无需通过试验测试等手段,在减少了资源浪费的情况下能够保证结果准确度和计算效率。
附图说明
图1是本发明实施方式的流程图;
图2是双向加筋平板结构示意图;
图3是纵向加筋平板结构示意图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
加筋板在较低频率时可等效为正交各向异性平板,平衡方程如下所示:
Figure 887260DEST_PATH_IMAGE024
求解加筋板等效参数,即求解
Figure 236378DEST_PATH_IMAGE026
等参数。
本发明的实施方式涉及一种基于势能理论的加筋板等效参数确定方法,如图1所示,包括以下步骤:基于势能理论构建加筋板弯曲振动最大势能公式和加筋板动能公式;将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式,得到平衡方程,并基于所述平衡方程计算加筋板的共振频率;令所述加筋板的共振频率与正交各向异性板的共振频率相等,得到所述加筋板等效为正交各向异性板的各个等效参数。
下面以双向加筋平板结构为例进一步说明本发明,其中,双向加筋平板结构如图2所示。
对于双向加筋的加筋板,根据势能理论,加筋板弯曲振动最大势能公式为:
Figure 907531DEST_PATH_IMAGE001
其中,
Figure 143340DEST_PATH_IMAGE002
为加筋板弯曲振动最大势能,
Figure 216338DEST_PATH_IMAGE003
表示平板的弯曲刚度,
Figure 810393DEST_PATH_IMAGE004
表示加筋板纵向的长度,
Figure 657257DEST_PATH_IMAGE005
表示加筋板横向的长度,
Figure 165599DEST_PATH_IMAGE006
分别表示加筋板纵向和横向坐标,
Figure 358683DEST_PATH_IMAGE007
表示原频率,
Figure 61322DEST_PATH_IMAGE008
表示加筋板材料的泊松比,
Figure 769384DEST_PATH_IMAGE009
表示加强筋材料的杨氏模量,
Figure 645198DEST_PATH_IMAGE010
表示加强筋材料的剪切模量,
Figure 692789DEST_PATH_IMAGE011
表示纵向加强筋的数量,
Figure 127181DEST_PATH_IMAGE012
表示横向加强筋的数量,
Figure 761687DEST_PATH_IMAGE013
分别表示纵向加强筋和横向加强筋的弯曲惯性矩,
Figure 877411DEST_PATH_IMAGE014
分别表示纵向加强筋和横向加强筋的扭转常数。上式中,等式右边三项中的第一项为无加筋平板的势能,第二项和第三项表示两个方向加强筋的势能。
加筋板的动能公式为:
Figure 779508DEST_PATH_IMAGE015
。其中,
Figure 322484DEST_PATH_IMAGE016
为加筋板动能,
Figure 178707DEST_PATH_IMAGE017
为加筋板的共振频率,
Figure 98122DEST_PATH_IMAGE018
分别为加筋板纵向和横向的模态阶数,
Figure 854725DEST_PATH_IMAGE019
为加筋板的面密度,
Figure 132385DEST_PATH_IMAGE021
表示纵向加强筋的线密度,
Figure 240018DEST_PATH_IMAGE022
表示横向加强筋的线密度,
Figure 697544DEST_PATH_IMAGE020
表示位移。
加筋板Z向位移根据边界条件可展开为傅里叶级数形式:
Figure 308654DEST_PATH_IMAGE028
,其中,
Figure 972196DEST_PATH_IMAGE030
表示位移幅值,
Figure 301547DEST_PATH_IMAGE032
表示加筋板纵向位移振型函数,
Figure 297184DEST_PATH_IMAGE034
表示加筋板横向位移振型函数。
将位移带入到加筋板弯曲振动最大势能公式和动能公式,并令二者相等得到平衡方程,通过求解平衡方程得到加筋板的共振频率:
Figure 28380DEST_PATH_IMAGE036
其中,
Figure 913422DEST_PATH_IMAGE038
表示纵向主振型函数的一阶和二阶导数,
Figure 402172DEST_PATH_IMAGE040
表示横向主振型函数的一阶和二阶导数。
正交各向异性板的共振频率也可以基于势能理论构建的势能公式和动能公式得到,具体为:
Figure 529397DEST_PATH_IMAGE042
Figure 616564DEST_PATH_IMAGE044
,即可得到加筋板等效为正交各向异性板的各个等效参数:
Figure 843146DEST_PATH_IMAGE046
对于单向加筋的加筋板,如图3所示,以只有平行X轴的加强筋为例进行说明。
Figure DEST_PATH_IMAGE048
两条边为自由边界,
Figure DEST_PATH_IMAGE050
可以得到:
Figure DEST_PATH_IMAGE052
Figure 835504DEST_PATH_IMAGE054
两条边简支约束的,并且梁间距相等,截面相同的情况下,
Figure 438523DEST_PATH_IMAGE056
对于n=1的模态,可以得到:
Figure 878732DEST_PATH_IMAGE058
不难发现,本发明基于势能理论构建加筋板的弯曲振动最大势能公式和动能公式,再结合加筋板的Z向位移构建平衡方程得到加筋板的共振频率,再结合正交各向异性板的共振频率得到加筋板等效为正交各向异性平板的等效参数,整个过程无需通过试验测试等手段,在减少了资源浪费的情况下能够保证结果准确度和计算效率。

Claims (2)

1.一种基于势能理论的加筋板等效参数确定方法,其特征在于,包括以下步骤:
基于势能理论构建加筋板弯曲振动最大势能公式和加筋板动能公式;所述加筋板弯曲振动最大势能公式为:
Figure 610897DEST_PATH_IMAGE001
,其中,
Figure 590355DEST_PATH_IMAGE002
为加筋板弯曲振动最大势能,
Figure 665365DEST_PATH_IMAGE003
表示平板的弯曲刚度,
Figure 388470DEST_PATH_IMAGE004
表示加筋板纵向的长度,
Figure 265159DEST_PATH_IMAGE005
表示加筋板横向的长度,
Figure 884359DEST_PATH_IMAGE006
分别表示加筋板纵向和横向坐标,
Figure 980753DEST_PATH_IMAGE007
表示原频率,
Figure 241971DEST_PATH_IMAGE008
表示加筋板材料的泊松比,
Figure 973166DEST_PATH_IMAGE009
表示加强筋材料的杨氏模量,
Figure 28847DEST_PATH_IMAGE010
表示加强筋材料的剪切模量,
Figure 78449DEST_PATH_IMAGE011
表示纵向加强筋的数量,
Figure 143357DEST_PATH_IMAGE012
表示横向加强筋的数量,
Figure 729059DEST_PATH_IMAGE013
分别表示纵向加强筋和横向加强筋的弯曲惯性矩,
Figure 955641DEST_PATH_IMAGE014
分别表示纵向加强筋和横向加强筋的扭转常数;所述加筋板动能公式为:
Figure 26628DEST_PATH_IMAGE015
,其中,
Figure 364068DEST_PATH_IMAGE016
为加筋板动能,
Figure 69856DEST_PATH_IMAGE017
为加筋板的共振频率,
Figure 467339DEST_PATH_IMAGE018
分别为加筋板纵向和横向的模态阶数,
Figure 491534DEST_PATH_IMAGE019
为加筋板的面密度,
Figure 898244DEST_PATH_IMAGE020
表示位移,
Figure 255276DEST_PATH_IMAGE021
表示纵向加强筋的线密度,
Figure 856284DEST_PATH_IMAGE022
表示横向加强筋的线密度;
将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式,得到平衡方程,并基于所述平衡方程计算加筋板的共振频率;
令所述加筋板的共振频率与正交各向异性板的共振频率相等,得到所述加筋板等效为正交各向异性板的各个等效参数。
2.根据权利要求1所述的基于势能理论的加筋板等效参数确定方法,其特征在于,所述将加筋板Z向位移带入到所述加筋板弯曲振动最大势能公式和所述加筋板动能公式前,需要将所述加筋板Z向位移根据边界条件展开为傅里叶级数形式。
CN202210850620.9A 2022-07-20 2022-07-20 一种基于势能理论的加筋板等效参数确定方法 Active CN115081148B (zh)

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 CN115081148A (zh) 2022-09-20
CN115081148B true 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)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115310188B (zh) * 2022-10-12 2023-01-10 上海索辰信息科技股份有限公司 一种基于最小势能理论确定加筋板等效参数的方法
CN115618163B (zh) * 2022-11-11 2023-05-12 上海索辰信息科技股份有限公司 一种加筋结构体固有频率计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013101598A (ja) * 2011-10-17 2013-05-23 Mizuho Information & Research Institute Inc 挙動解析システム、挙動解析方法及び挙動解析プログラム
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 中国航空工业集团公司西安飞机设计研究所 一种单向轴压矩形板后屈曲平衡路径计算方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013065764A1 (ja) * 2011-10-31 2013-05-10 国立大学法人大阪大学 解析装置、解析方法及び解析プログラム
CN108256214B (zh) * 2018-01-16 2021-07-02 滨州学院 一种航空正交加筋板刚度计算方法及装置
CN112378735B (zh) * 2020-10-28 2021-11-09 华南理工大学 一种考虑正交异性板残余应力效应的组合板梁单元分析法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013101598A (ja) * 2011-10-17 2013-05-23 Mizuho Information & Research Institute Inc 挙動解析システム、挙動解析方法及び挙動解析プログラム
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 中国航空工业集团公司西安飞机设计研究所 一种单向轴压矩形板后屈曲平衡路径计算方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
典型边界条件下加筋矩形板的横向振动特性分析;李国荣等;《振动与冲击》;20200825(第16期);第266-271页 *
基于FEM-MODENA的加筋板声-固耦合分析;于士甲等;《振动工程学报》;20170415(第02期);第100-107页 *
李凯 ; 何书韬 ; 邱永康 ; 吴国民 ; 郭文杰.附加多个集中质量加筋板的自由振动分析.《中国舰船研究》.2015, *
杨健生 ; 曾治平 ; 韦冬炎 ; 彭林欣.基于无网格法的非均匀弹性地基上变厚度加筋板弯曲与固有频率分析.《计算力学学报》.2020, *

Also Published As

Publication number Publication date
CN115081148A (zh) 2022-09-20

Similar Documents

Publication Publication Date Title
CN115081148B (zh) 一种基于势能理论的加筋板等效参数确定方法
Zhang et al. On the dynamic response of porous functionally graded microbeam under moving load
Tornabene et al. Higher order formulations for doubly-curved shell structures with a honeycomb core
García-Macías et al. Static and free vibration analysis of functionally graded carbon nanotube reinforced skew plates
Mohammadi et al. An equivalent model for trapezoidal corrugated cores based on homogenization method
CN104537259B (zh) 一种基于细观力学模型的横观各向同性纤维材料性能的预测方法
Ji et al. On the importance of work-conjugacy and objective stress rates in finite deformation incremental finite element analysis
Guo et al. Nonlinear dynamics of fluid-conveying composite pipes subjected to time-varying axial tension in sub-and super-critical regimes
CN104778377B (zh) 一种组合梁弯曲振动的固有频率分析方法
US20070075450A1 (en) Parametrized material and performance properties based on virtual testing
CN101887474B (zh) 基于有限元法与广义傅里叶级数法的结构振动分析方法
CN107368634A (zh) 一种复杂结构耦合损耗因子的预示方法
CN104316388A (zh) 一种对各向异性材料结构件进行疲劳寿命测定的方法
CN110688789B (zh) 一种碳纤维复合单向板偏载冲击的多尺度响应分析方法
Tian et al. Analysis of nonlinear aeroelastic characteristics of a trapezoidal wing in hypersonic flow
Zhuang et al. Modal and aeroelastic analysis of trapezoidal corrugated-core sandwich panels in supersonic flow
Jafari-Talookolaei et al. In-plane and out-of-plane vibration modes of laminated composite beams with arbitrary lay-ups
CN107268820A (zh) 无屈曲波形结构耗能构件及其设计方法
CN112784359B (zh) 薄壁梁约束扭转极限承载力迭代计算方法
Su et al. Modeling and parametric analysis of in-plane free vibration of a floating cable-stayed bridge with transfer matrix method
Zheng et al. New analytic buckling solutions of non-Lévy-type cylindrical panels within the symplectic framework
Li et al. On the non-linear dynamics of a forced plate with boundary conditions correction in subsonic flow
Branner et al. Torsional performance of wind turbine blades–Part II: Numerical validation
Chwalowski et al. Progress on Transonic Flutter and Shock Buffet Computations in Support of the Third Aeroelastic Prediction Workshop
CN104598693A (zh) 一种确定薄壁结构高刚度连接区载荷传递的方法

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