CN112329146B - 一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法 - Google Patents
一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法 Download PDFInfo
- Publication number
- CN112329146B CN112329146B CN202011223141.1A CN202011223141A CN112329146B CN 112329146 B CN112329146 B CN 112329146B CN 202011223141 A CN202011223141 A CN 202011223141A CN 112329146 B CN112329146 B CN 112329146B
- Authority
- CN
- China
- Prior art keywords
- load
- full
- equivalent
- spread
- deck
- 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/15—Vehicle, aircraft or watercraft design
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
Abstract
本发明涉及一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法,属于船舶结构设计技术领域。本发明在甲板板格尺寸、材料力学特性、轮印载荷大小与尺寸等相关参数信息确定下,利用有限元软件分别计算得到不同位置处多轮印载荷和满铺均布载荷作用下甲板板格的塑性变形值,利用板格等价变形媒介参数,得到多轮印载荷与满铺均布载荷之间的对应关系,从而得到多次轮印载荷与满铺均布载荷之间的等效比例因子计算公式。有效解决不同轮印载荷及尺寸下、不同板格尺寸下的甲板多轮印载荷等效问题;本发明设计方法完整、适用,可以结合甲板的实际情况进行计算,易于实现。
Description
技术领域
本发明涉及一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法,属于船舶结构设计技术领域。
背景技术
现代舰船甲板遭受轮印载荷作用的情况越来越多,例如,直升机甲板遭受直升机轮印载荷作用,车辆甲板遭受重载车辆轮印载荷等。在甲板使用寿命期间,轮印载荷可能作用在甲板板格不同位置上,由此甲板板格会产生累积的塑性变形,这将会严重影响甲板的工作性能和安全性能。因此,不同位置处多轮印载荷作用下甲板板格的塑性变形计算方法研究显得至关重要。
由于甲板遭受轮印载荷的大小和位置都是随机的,而且是反复多次的,极大增加了轮印载荷作用下甲板板格塑性变形计算问题的困难性和复杂性。目前针对不同位置处多轮印载荷作用下甲板板格塑性变形的计算方法比较缺乏,提出多次轮印载荷与满铺均布载荷等效方法显得尤为重要。
发明内容
本发明的目的是为解决针对不同位置处多轮印载荷作用下甲板板格塑性变形的计算方法技术问题。
为达到解决上述问题的目的,本发明所采取的技术方案是提供一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法,包括以下步骤:
步骤1:确定甲板板格尺寸、材料力学特性和轮印载荷大小的相关参数信息,利用有限元软件计算出板格上不同位置处多次轮印载荷作用下甲板板格塑性变形值大小,并进行转换得到板格的等价变形媒介值;
步骤2:利用有限元软件计算出对应的相同甲板板格尺寸下一系列满铺均布载荷下的板格永久变形值,并转换为板格的等价变形媒介值,从而得到满铺均布载荷与板格的等价变形媒介值之间的关系曲线;
步骤3:根据多次轮印载荷作用下得到的等价变形媒介参数值,从满铺均布载荷与板格等价变形媒介值之间的关系曲线中插值确定等效满铺均布载荷值;
步骤4:根据多次轮印载荷和等效满铺均布载荷参数之间的对应关系,得到多次轮印载荷与满铺载荷之间的等效比例因子,从而给出多次轮印载荷与满铺载荷等效方法。
优选地,所述步骤1、步骤2和步骤3中板格等价变形媒介可以描述为式(1):
优选地,所述步骤1中轮印载荷可以描述为式(2):
式中Qp表示无量纲轮印载荷参数,P表示为轮印集中力,单位N,E表示为甲板板格材料的杨氏模量,单位MPa,s表示为甲板纵骨间距,单位mm,σy表示为甲板板格材料的屈服强度,单位MPa。
优选地,所述步骤1中对确定尺寸、材料的甲板板格,建立双线性材料参数单板模型,施加四周刚性固定的边界条件。采用ABAQUS商业软件,在轮印区域施加局部均布载荷,在进行不同位置处多次轮印载荷作用下甲板板格塑性变形值计算时,采用多载荷步计算分析。
优选地,所述步骤2中满铺均布载荷可以描述为式(3):
式中Qe表示无量纲满铺均布载荷参数,其中Pe表示为满铺均布载荷,单位MPa。
优选地,所述步骤2中甲板板格与步骤1中甲板板格计算模型相同,在整个甲板板格内施加满铺均布载荷,对满铺均布载荷作用下甲板板格塑性变形值进行计算。
优选地,所述步骤4中多次轮印载荷与满铺均布载荷之间的等效比例因子可以描述为式(4):
式中r表示为多次轮印载荷与满铺均布载荷之间的等效比例因子,λ表示为载荷集中参数,α表示甲板板格的长宽比。
优选地,所述步骤4中的载荷集中参数λ可以描述为式(5):
式中a表示轮印长度,单位mm,b表示轮印宽度,单位mm。
相比现有技术,本发明具有如下有益效果:
本发明针对现有技术存在的不足,提供了一种基于甲板塑性变形关系下多次轮印载荷与满铺均布载荷等效方法,能有效解决不同轮印尺寸下、不同轮印载荷作用位置下以及不同甲板板格尺寸的轮印载荷等效问题。
本发明在甲板板格尺寸、材料力学特性、轮印载荷大小与尺寸等相关参数信息确定下,利用有限元软件分别计算得到不同位置处多轮印载荷和满铺均布载荷作用下甲板板格的塑性变形值,利用板格等价变形媒介参数,得到多轮印载荷与满铺均布载荷之间的对应关系,从而得到多次轮印载荷与满铺均布载荷之间的等效比例因子计算公式。它能有效解决不同轮印载荷及尺寸下、不同板格尺寸下的甲板多轮印载荷等效问题,根据基于甲板塑性变形关系下多次轮印载荷与满铺均布载荷等效方法,可以计算出多轮印载荷所对应的满铺均布载荷,并解决了不同位置处多轮印载荷作用下甲板板格的塑性变形计算问题,简单方便,准确高效。本发明所提供的设计方法完整、适用,可以结合甲板的实际情况进行计算,易于实现。
附图说明
图1为单个位置集中载荷的几何特征图
图2为甲板板格双线性材料模型
图4为l/s=3.428下等效因子r与1/β的关系图
具体实施方式
为使本发明更明显易懂,兹以优选实施例,并配合附图作详细说明如下:如图1-4所示,本发明提供一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法,包括以下步骤:
步骤1:确定甲板板格尺寸、材料力学特性和轮印载荷大小的相关参数信息,利用有限元软件计算出板格上不同位置处多次轮印载荷作用下甲板板格塑性变形值大小,并进行转换得到板格的等价变形媒介值;
步骤2:利用有限元软件计算出对应的相同甲板板格尺寸下一系列满铺均布载荷下的板格永久变形值,并转换为板格的等价变形媒介值,从而得到满铺均布载荷与板格的等价变形媒介值之间的关系曲线;
步骤3:根据多次轮印载荷作用下得到的等价变形媒介参数值,从满铺均布载荷与板格等价变形媒介值之间的关系曲线中插值确定等效满铺均布载荷值;
步骤4:根据多次轮印载荷和等效满铺均布载荷参数之间的对应关系,得到多次轮印载荷与满铺载荷之间的等效比例因子,从而给出多次轮印载荷与满铺载荷等效方法。
所述步骤1、步骤2和步骤3中板格等价变形媒介可以描述为式(1):
所述步骤1中轮印载荷可以描述为式(2):
式中Qp表示无量纲轮印载荷参数,P表示为轮印集中力,单位N,E表示为甲板板格材料的杨氏模量,单位MPa,s表示为甲板纵骨间距,单位mm,σy表示为甲板板格材料的屈服强度,单位MPa。
所述步骤1中对确定尺寸、材料的甲板板格,建立双线性材料参数单板模型,施加四周刚性固定的边界条件。采用ABAQUS商业软件,在轮印区域施加局部均布载荷,在进行不同位置处多次轮印载荷作用下甲板板格塑性变形值计算时,采用多载荷步计算分析。
所述步骤2中满铺均布载荷可以描述为式(3):
式中Qe表示无量纲满铺均布载荷参数,其中Pe表示为满铺均布载荷,单位MPa。
优选地,所述步骤2中甲板板格与步骤1中甲板板格计算模型相同,在整个甲板板格内施加满铺均布载荷,对满铺均布载荷作用下甲板板格塑性变形值进行计算。
所述步骤4中多次轮印载荷与满铺均布载荷之间的等效比例因子可以描述为式(4):
式中r表示为多次轮印载荷与满铺均布载荷之间的等效比例因子,λ表示为载荷集中参数,α表示甲板板格的长宽比。
所述步骤4中的载荷集中参数λ可以描述为式(5):
式中a表示轮印长度,单位mm,b表示轮印宽度,单位mm。
为了进一步说明上述基于甲板塑性变形关系下多次轮印载荷与满铺均布载荷等效方法,上述步骤1中的所包含甲板板格、轮印载荷等信息的计算模型如图1所示,图1中l表示甲板板格长度,s表示甲板板格宽度或者甲板纵骨间距,a和b分别表示轮印长度和轮印宽度,其中轮印载荷位置可在板格任意位置上出现,此外,甲板板格材料模型如图2所示,选定材料的相关具体力学参数如表1所示。
表1
所述步骤1中在确定甲板板格尺寸、材料力学特性和轮印载荷大小的相关参数信息之后,需利用有限元软件计算出不同板格尺寸、不同轮印尺寸以及不同位置处多次轮印载荷作用下甲板板格塑性变形值,图3中给出了特定板格长宽比l/s=2和特定板厚t=14mm下不同轮印尺寸以及不同板格等价变形媒介值下所对应的无量纲轮印载荷Qp值。
式中t表示甲板板厚。
上述的步骤2在确定甲板板格尺寸、材料力学特性等相关参数信息之后,需利用有限元软件计算出不同板格尺寸、不同轮印尺寸以及不同满铺均布载荷作用下甲板板格塑性变形值,图3给出了特定板格长宽比l/s=2和特定板厚t=14mm下不同板格等价变形媒介值下所对应的无量纲轮印载荷Qe值。
所述的步骤3在板格等效变形媒介参数确定下,可以从满铺均布载荷与板格等价变形媒介值之间的关系曲线中插值确定多轮印载荷Qe所对应的等效满铺均布载荷Qp值,图3给出了特定板格长宽比l/s=2和特定板厚t=14mm下相同板格等效变形媒介参数下多轮印载荷Qe与等效满铺均布载荷Qp值。
r=Qe/QP 式(7)
根据步骤3中的多次轮印载荷和等效满铺均布载荷参数之间的对应关系,可以计算多次轮印载荷与满铺载荷之间的等效比例因子r。图4给出了特定板格长宽比l/s=3.428下的两种载荷集中参数λ下的多次轮印载荷与满铺载荷之间的等效比例因子r与板格的柔度系数β之间的关系,可以发现,r与1/β成正比关系,为此,通过计算不同板格尺寸和不同轮印尺寸下等效比例因子r与板格的柔度系数β之间的关系,可以拟合得到等效比例因子r的计算公式,如式(4)所示。式(4)中的板格长宽比参数α可以描述为式(8):
α=l/s 式(8)。
综上所述,本发明在甲板板格尺寸、材料力学特性、轮印载荷大小与尺寸等相关参数信息确定下,利用有限元软件分别计算得到不同位置处多轮印载荷和满铺均布载荷作用下甲板板格的塑性变形值,利用板格等价变形媒介参数,得到多轮印载荷与满铺均布载荷之间的对应关系,从而得到多次轮印载荷与满铺均布载荷之间的等效比例因子计算公式。它能有效解决不同轮印载荷及尺寸下、不同板格尺寸下的甲板多轮印载荷等效问题,根据基于甲板塑性变形关系下多次轮印载荷与满铺均布载荷等效方法,可以计算出多轮印载荷所对应的满铺均布载荷,并解决了不同位置处多轮印载荷作用下甲板板格的塑性变形计算问题,简单方便,准确高效。本发明所提供的设计方法完整、适用,可以结合甲板的实际情况进行计算,易于实现。
以上所述,仅为本发明的较佳实施例,并非对本发明任何形式上和实质上的限制,应当指出,对于本技术领域的普通技术人员,在不脱离本发明的前提下,还将可以做出若干改进和补充,这些改进和补充也应视为本发明的保护范围。凡熟悉本专业的技术人员,在不脱离本发明的精神和范围的情况下,当可利用以上所揭示的技术内容而做出的些许更动、修饰与演变的等同变化,均为本发明的等效实施例;同时,凡依据本发明的实质技术对上述实施例所作的任何等同变化的更动、修饰与演变,均仍属于本发明的技术方案的范围内。
Claims (3)
1.一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法,其特征在于:包括以下步骤:
步骤1:确定甲板板格尺寸、材料力学特性和轮印载荷大小的相关参数信息,利用有限元软件计算出板格上不同位置处多次轮印载荷作用下甲板板格塑性变形值大小,并进行转换得到板格的等价变形媒介值;
板格等价变形媒介可以描述为式(1):
轮印载荷可以描述为式(2):
式中Qp表示无量纲轮印载荷参数,P表示为轮印集中力,单位N,E表示为甲板板格材料的杨氏模量,单位MPa,s表示为甲板纵骨间距,单位mm,σy表示为甲板板格材料的屈服强度,单位MPa;
步骤2:利用有限元软件计算出对应的相同甲板板格尺寸下一系列满铺均布载荷下的板格永久变形值,并转换为板格的等价变形媒介值,从而得到满铺均布载荷与板格的等价变形媒介值之间的关系曲线;
满铺均布载荷可以描述为式(3):
式中Qe表示无量纲满铺均布载荷参数,其中Pe表示为满铺均布载荷,单位MPa;
步骤3:根据多次轮印载荷作用下得到的等价变形媒介参数值,从满铺均布载荷与板格等价变形媒介值之间的关系曲线中插值确定等效满铺均布载荷值;
步骤4:根据多次轮印载荷和等效满铺均布载荷参数之间的对应关系,得到多次轮印载荷与满铺载荷之间的等效比例因子,从而给出多次轮印载荷与满铺载荷等效方法;
式中r表示为多次轮印载荷与满铺均布载荷之间的等效比例因子,λ表示为载荷集中参数,α表示甲板板格的长宽比;
载荷集中参数λ可以描述为式(5):
式中a表示轮印长度,单位mm,b表示轮印宽度,单位mm。
2.如权利要求1所述的一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法,其特征在于:所述步骤1中对确定尺寸、材料的甲板板格,建立双线性材料参数单板模型,施加四周刚性固定的边界条件;采用ABAQUS商业软件,在轮印区域施加局部均布载荷,在进行不同位置处多次轮印载荷作用下甲板板格塑性变形值计算时,采用多载荷步计算分析。
3.如权利要求1所述的一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法,其特征在于:所述步骤2中甲板板格与步骤1中甲板板格计算模型相同,在整个甲板板格内施加满铺均布载荷,对满铺均布载荷作用下甲板板格塑性变形值进行计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011223141.1A CN112329146B (zh) | 2020-11-05 | 2020-11-05 | 一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011223141.1A CN112329146B (zh) | 2020-11-05 | 2020-11-05 | 一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112329146A CN112329146A (zh) | 2021-02-05 |
CN112329146B true CN112329146B (zh) | 2023-03-28 |
Family
ID=74317139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011223141.1A Active CN112329146B (zh) | 2020-11-05 | 2020-11-05 | 一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112329146B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984885A (zh) * | 2018-07-05 | 2018-12-11 | 中国船舶工业集团公司第七0八研究所 | 一种基于许可永久变形下的装载甲板板设计方法 |
CN111824356A (zh) * | 2020-07-29 | 2020-10-27 | 中国船舶工业集团公司第七0八研究所 | 一种轮印载荷作用下的i型夹层结构设计方法 |
-
2020
- 2020-11-05 CN CN202011223141.1A patent/CN112329146B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984885A (zh) * | 2018-07-05 | 2018-12-11 | 中国船舶工业集团公司第七0八研究所 | 一种基于许可永久变形下的装载甲板板设计方法 |
CN111824356A (zh) * | 2020-07-29 | 2020-10-27 | 中国船舶工业集团公司第七0八研究所 | 一种轮印载荷作用下的i型夹层结构设计方法 |
Non-Patent Citations (2)
Title |
---|
Design Method for Steel Deck Plates under Quasi-static Patch Loads with Allowable Plastic Deformations;Zhu L等;《 Marine Structures》;20200511;全文 * |
轮印载荷下车辆甲板板厚设计方法研究;王维舟;《硕士论文工程科技Ⅱ辑》;20200315;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112329146A (zh) | 2021-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Eberhardt | Some developments in'distance sampling' | |
CN106951617A (zh) | 一种氯离子在混凝土中反常扩散动态数据重构的分形导数模拟方法 | |
CN105975711B (zh) | 一种基于材料损伤累积水平的沥青路面寿命评估方法 | |
CN112329146B (zh) | 一种基于甲板塑性变形关系下多轮印载荷与满铺载荷等效方法 | |
CN111709085A (zh) | 一种约束阻尼薄板结构拓扑优化设计方法 | |
CN112784359B (zh) | 薄壁梁约束扭转极限承载力迭代计算方法 | |
Perfect et al. | Multifractal model for soil aggregate fragmentation | |
CN106940736A (zh) | 一种考虑侧向载荷作用的纵骨多跨失稳的极限载荷确定方法 | |
CN105117577A (zh) | 输电铁塔用八地脚螺栓刚性塔座板规格计算方法 | |
CN107704657A (zh) | 一种由残余应力释放引起的板类零件翘曲变形的计算方法 | |
CN110044731B (zh) | 非均布法向荷载结构面直剪试验方法及辅助加载装置 | |
CN110516293A (zh) | 船舯结构弯剪扭组合作用下极限承载能力的有限元计算方法 | |
CN109446658A (zh) | 一种钛合金锥形压痕约束因子和代表应变的确定方法 | |
CN114912302B (zh) | 一种基于模态空间采样算法的加筋板模态密度获取方法 | |
CN110489916B (zh) | 基于损伤状态倾角影响线曲率的等截面梁损伤识别方法 | |
CN107330894A (zh) | 面向平面舱壁类型船舱的点云数据自动分割方法 | |
CN110489795A (zh) | 基于正则化广义逆的结构破损诊断方法 | |
CN103033467A (zh) | 二维多孔材料共面动态缓冲性能的测定方法 | |
CN108984885B (zh) | 一种基于许可永久变形下的装载甲板板设计方法 | |
CN113378354B (zh) | 一种考虑剪力作用的薄壁梁弯曲极限强度计算方法 | |
Dokainish et al. | Vibrations of orthotropic parallelogramic plates with variable thickness | |
CN110852013A (zh) | 基于结构力学法的加高重力坝新老混凝土分载计算方法 | |
CN107609231B (zh) | 一种蜗杆磨削齿轮表面微观形貌仿真方法和系统 | |
McCutcheon et al. | Simplified design procedure for glued-laminated bridge decks | |
CN105608291A (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 |