CN110162881A - 一种弯剪扭组合下船舯结构极限承载能力的确定方法 - Google Patents

一种弯剪扭组合下船舯结构极限承载能力的确定方法 Download PDF

Info

Publication number
CN110162881A
CN110162881A CN201910431249.0A CN201910431249A CN110162881A CN 110162881 A CN110162881 A CN 110162881A CN 201910431249 A CN201910431249 A CN 201910431249A CN 110162881 A CN110162881 A CN 110162881A
Authority
CN
China
Prior art keywords
bearing capacity
bending
finite element
element analysis
torsion
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
Application number
CN201910431249.0A
Other languages
English (en)
Other versions
CN110162881B (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.)
Zhejiang University of Technology ZJUT
708th Research Institute of CSIC
Original Assignee
Zhejiang University of Technology ZJUT
708th Research Institute of CSIC
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 Zhejiang University of Technology ZJUT, 708th Research Institute of CSIC filed Critical Zhejiang University of Technology ZJUT
Priority to CN201910431249.0A priority Critical patent/CN110162881B/zh
Publication of CN110162881A publication Critical patent/CN110162881A/zh
Application granted granted Critical
Publication of CN110162881B publication Critical patent/CN110162881B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及船舶结构设计领域，尤其涉及一种弯剪扭组合下船舯结构极限承载能力的确定方法，包括：建立非线性有限元分析模型，依据有限元分析模型计算模型的约束扭转极限强度，依据有限元分析模型计算模型的自由扭转极限强度，计算约束扭转和自由扭转的极限承载能力之差，依据约束扭转极限承载能力和自由扭转极限承载能力的差计算模型在弯剪扭组合加载下弯矩和剪力的极限强度并获得变形云图和载荷‑位移曲线，依据有限元分析模型在弯剪扭组合下弯矩和剪力的极限承载能力计算模型在弯剪扭组合加载下的自由扭转极限强度并获得载荷‑位移曲线。本发明用于确定弯剪扭组合加载作用下船舯区域船体梁的极限承载能力。

Description

一种弯剪扭组合下船舯结构极限承载能力的确定方法
技术领域
本发明涉及船舶结构设计领域，尤其涉及一种弯剪扭组合下船舯结构极限承载能力的确定方法。
背景技术
随着计算技术的迅速发展，船舶结构极限强度的评估越来越多地采用有限元方法，而这一工作的首要问题就是如何建立有限元模型。垂向弯矩作用下船体梁的极限强度最受关注也最为成熟。对于具有大开口、扭转刚度较低的集装箱船，弯剪扭的组合作用尤其重要。船体的扭转分析表明船舯双力矩较小，翘曲应力也较小，以自由扭转为主，因此船舯区域主要承受自由扭转和弯剪组合荷载作用。然而，由于问题的复杂性，该位置极限强度的理论和试验研究并不多见。
发明内容
为了解决以上技术问题，本发明提供了一种弯剪扭组合下船舯结构极限承载能力的确定方法，通过寻找船舯区域扭转和弯剪组合下的边界条件和确定方法，用有限元分析来进行研究。
本发明所解决的技术问题可以采用以下技术方案来实现：
一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，包括：
步骤S1，取船舯船体梁横截面模型建立非线性有限元分析模型；
步骤S2，依据所述有限元分析模型计算约束扭转极限承载能力；
步骤S3，依据所述有限元分析模型计算自由扭转极限承载能力；
步骤S4，计算所述约束扭转极限承载能力和所述自由扭转极限承载能力的差；
步骤S5，依据所述约束扭转极限承载能力和所述自由扭转极限承载能力的差计算所述有限元分析模型在弯剪扭组合下弯矩和剪力的极限承载能力，并获得变形云图和载荷-位移曲线；
步骤S6，依据所述有限元分析模型在弯剪扭组合下弯矩和剪力的极限承载能力计算所述有限元分析模型在弯剪扭组合下的所述自由扭转极限承载能力，并获得载荷-位移曲线。
优选的，所述步骤S1中，所述有限元分析模型取船舯船体梁横截面模型，长度方向采用一跨范围，所述有限元分析模型中不含横向构件。
优选的，所述步骤S2中，所述有限元分析模型采用一端x、y、z方向线位移约束，另一端在形心位置作用扭矩，并通过所述有限元分析模型计算所述约束扭转极限值。
优选的，所述步骤S3中包括：
步骤S31,对所述有限元分析模型非加载端面的所有纵向构件节点进行y、z线位移约束。
优选的，所述步骤S3中还包括：
步骤S32，对加载端面的所有纵向构件节点与关联点的y、z线位移、x角位移进行耦合关联。
优选的，所述步骤S3中还包括：
步骤S33，在扭心处施加x方向的扭矩。
优选的，所述步骤S3中还包括：
步骤S34,对外底板中心线上的节点约束x、y、z线位移。
优选的，所述步骤S5中包括：
步骤S51，一跨所述有限元分析模型采用一端x、y、z方向线位移约束，另一端在形心位置作用所述步骤S4中计算得到的所述约束扭转与所述自由扭转极限承载能力之差。
优选的，所述步骤S5中还包括：
步骤S52，在所述步骤S51的基础上，在加载端的形心位置处作用所述扭矩、所述弯矩和所述剪力，通过所述有限元分析模型计算可得组合状态下的所述弯矩、所述剪力极限值，并获得载荷-位移曲线及变形云图。
优选的，所述步骤S6中，将所述步骤S5中计算得到的弯剪扭组合下所述弯矩和所述剪力的极限承载能力值减去所述约束扭转与所述自由扭转的扭矩值差，并获得载荷-位移曲线。
其有益效果在于：
一种弯剪扭组合下船舯结构极限承载能力的确定方法，可用以确定弯剪扭组合加载作用下船舯区域船体梁的极限承载能力，减轻计算问题的复杂性。
附图说明
图1为非线性有限元分析模型和非线性有限元分析模型截面；
图2为约束扭转的边界条件和加载示意；
图3为自由扭转的边界条件和加载示意；
图4为船舯区域弯剪扭组合作用下的边界条件和加载示意；
图5为船舯区域弯剪扭组合作用下的弯矩-转角曲线；
图6为船舯区域弯剪扭组合作用下的剪力-位移曲线；
图7为船舯区域弯剪扭组合作用下的扭矩-转角曲线；
图8为船舯区域弯剪扭极限状态下的变形云图；
图9为一种弯剪扭组合下船舯结构极限承载能力的确定方法步骤框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
下面结合附图和具体实施例对本发明作进一步说明，但不作为本发明的限定。
参照图9为本发明提供的一种弯剪扭组合下船舯结构极限承载能力的确定方法步骤框图，包括：
步骤S1，取船舯船体梁横截面模型建立非线性有限元分析模型；
进一步的，鉴于船体梁舱段模型能准确地模拟船舯区域结构在弯剪扭组合作用下极限承载能力，但缺点是模型大，非线性计算时间长且收敛性差，因此有限元分析模型取船舯船体梁横截面模型，长度方向采用1跨范围，即强框架间距范围，模型中不含横向构件。参照图1按照实际试验的模型尺寸建立有限元分析模型，材料属性和实际试验模型一致，如表1：
单位；mm
表1
有限元网格划分如下：水平板纵骨间划分18个单元，侧壁纵骨间划分27个单元，腹板划分3个单元，扶强材翼缘为梁单元。
进一步的，步骤S2，依据所述有限元分析模型计算约束扭转极限承载能力。
参照图2模型采用一端x、y、z方向线位移约束，一端在形心位置作用扭矩，通过有限元分析模型计算约束扭转极限值。
进一步的，步骤S3，依据所述有限元分析模型计算自由扭转极限承载能力。
参照图3对有限元分析模型非加载端面的所有纵向构件节点进行y、z线位移约束，加载端面的所有纵向构件节点与关联点的y、z线位移、x角位移进行耦合关联，并在扭心处施加x方向的扭矩。为防止刚体运动，对外底板中心线上的节点约束x、y、z线位移；亦或采用理论方法计算梁的自由扭转极限值。
进一步的，步骤S4，计算所述约束扭转极限承载能力和所述自由扭转极限承载能力的差；
进一步的，步骤S5，依据所述约束扭转极限承载能力和所述自由扭转极限承载能力的差计算所述有限元分析模型在弯剪扭组合下弯矩和剪力的极限承载能力，并获得变形云图和载荷-位移曲线。
参照图4自由扭转和弯剪组合下的极限承载能力计算通过对两个分析步作用荷载来实现。第一个分析步中，一跨有限元分析模型采用一端x、y、z方向线位移约束，另一端在形心位置作用步骤S4中计算得到的约束扭转与自由扭转极限承载能力之差；第二个分析步在第一步的基础上，在加载端的形心位置处作用扭矩、弯矩和剪力值；通过确定可得组合状态下的弯矩、剪力极限值，并获得载荷-位移曲线和变形云图。
本发明较佳的实施例中：
(1)极限强度计算结果汇总，如表2：
表2
(2)弯矩-转角曲线，参照图5；
(3)剪力-位移曲线，参照图6；
(4)极限状态下的变形云图，参照图8。
进一步的，步骤S6，依据所述有限元分析模型在弯剪扭组合下弯矩和剪力的极限承载能力计算所述有限元分析模型在弯剪扭组合下的所述自由扭转极限承载能力，并获得载荷-位移曲线，如表2所示，并获得载荷-位移曲线，参照图8。
以上所述仅为本发明较佳的实施例，并非因此限制本发明的实施方式及保护范围，对于本领域技术人员而言，应当能够意识到凡运用本发明说明书及图示内容所作出的等同替换和显而易见的变化所得到的方案，均应当包含在本发明的保护范围内。

Claims (10)

1.一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，包括：
步骤S1，取船舯船体梁横截面模型建立非线性有限元分析模型；
步骤S2，依据所述有限元分析模型计算约束扭转极限承载能力；
步骤S3，依据所述有限元分析模型计算自由扭转极限承载能力；
步骤S4，计算所述约束扭转极限承载能力和所述自由扭转极限承载能力的差；
步骤S5，依据所述约束扭转极限承载能力和所述自由扭转极限承载能力的差计算所述有限元分析模型在弯剪扭组合下弯矩和剪力的极限承载能力，并获得变形云图和载荷-位移曲线；
步骤S6，依据所述有限元分析模型在弯剪扭组合下弯矩和剪力的极限承载能力计算所述有限元分析模型在弯剪扭组合下的所述自由扭转极限承载能力，并获得载荷-位移曲线。
2.根据权利要求1所述的一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，所述步骤S1中，所述有限元分析模型取船舯船体梁横截面模型，长度方向采用一跨范围，所述有限元分析模型中不含横向构件。
3.根据权利要求1所述的一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，所述步骤S2中，所述有限元分析模型采用一端x、y、z方向线位移约束，另一端在形心位置作用扭矩，并通过所述有限元分析模型计算所述约束扭转极限值。
4.根据权利要求1所述的一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，所述步骤S3中包括：
步骤S31,对所述有限元分析模型非加载端面的所有纵向构件节点进行y、z线位移约束。
5.根据权利要求1所述的一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，所述步骤S3中还包括：
步骤S32，对加载端面的所有纵向构件节点与关联点的y、z线位移、x角位移进行耦合关联。
6.根据权利要求1所述的一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，所述步骤S3中还包括：
步骤S33，在扭心处施加x方向的扭矩。
7.根据权利要求1所述的一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，所述步骤S3中还包括：
步骤S34,对外底板中心线上的节点约束x、y、z线位移。
8.根据权利要求4所述的一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，所述步骤S5中包括：
步骤S51，一跨所述有限元分析模型采用一端x、y、z方向线位移约束，另一端在形心位置作用所述步骤S4中计算得到的所述约束扭转与所述自由扭转极限承载能力之差。
9.根据权利要求5所述的一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，所述步骤S5中还包括：
步骤S52，在所述步骤S51的基础上，在加载端的形心位置处作用所述扭矩、所述弯矩和所述剪力，通过所述有限元分析模型计算可得组合状态下的所述弯矩、所述剪力极限值，并获得载荷-位移曲线及变形云图。
10.根据权利要求1所述的一种弯剪扭组合下船舯结构极限承载能力的确定方法，其特征在于，所述步骤S6中，将所述步骤S5中计算得到的弯剪扭组合下所述弯矩和所述剪力的极限承载能力值减去所述约束扭转与所述自由扭转的扭矩值差，并获得载荷-位移曲线。
CN201910431249.0A 2019-05-22 2019-05-22 一种弯剪扭组合下船舯结构极限承载能力的确定方法 Active CN110162881B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910431249.0A CN110162881B (zh) 2019-05-22 2019-05-22 一种弯剪扭组合下船舯结构极限承载能力的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910431249.0A CN110162881B (zh) 2019-05-22 2019-05-22 一种弯剪扭组合下船舯结构极限承载能力的确定方法

Publications (2)

Publication Number Publication Date
CN110162881A true CN110162881A (zh) 2019-08-23
CN110162881B CN110162881B (zh) 2023-05-16

Family

ID=67632010

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910431249.0A Active CN110162881B (zh) 2019-05-22 2019-05-22 一种弯剪扭组合下船舯结构极限承载能力的确定方法

Country Status (1)

Country Link
CN (1) CN110162881B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111241727A (zh) * 2020-01-03 2020-06-05 浙江工业大学 一种利用单跨有限元模型进行豪华邮轮极限强度计算的方法
CN112818480A (zh) * 2021-01-18 2021-05-18 浙江工业大学 压扭组合作用下横撑极限强度计算方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3956549A (en) * 1973-01-09 1976-05-11 Helmut Stoeberl Boat hull connector and method
CN1582240A (zh) * 2001-10-16 2005-02-16 马太及肖特工程师公司 混合结构的船
US20060288255A1 (en) * 2003-03-17 2006-12-21 Noboru Kikuchi Process of estimating relationship between element distortion and analysis error
CN101900648A (zh) * 2010-06-30 2010-12-01 哈尔滨工业大学深圳研究生院 多种受力情况下的钢管混凝土承载力确定方法及应用
CN206634167U (zh) * 2017-01-19 2017-11-14 中山市宏鸿船舶修造有限公司 一种开体驳船
CN109141820A (zh) * 2018-10-08 2019-01-04 中国船舶科学研究中心（中国船舶重工集团公司第七0二研究所） 船体模型扭矩及剪力测量方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3956549A (en) * 1973-01-09 1976-05-11 Helmut Stoeberl Boat hull connector and method
CN1582240A (zh) * 2001-10-16 2005-02-16 马太及肖特工程师公司 混合结构的船
US20060288255A1 (en) * 2003-03-17 2006-12-21 Noboru Kikuchi Process of estimating relationship between element distortion and analysis error
CN101900648A (zh) * 2010-06-30 2010-12-01 哈尔滨工业大学深圳研究生院 多种受力情况下的钢管混凝土承载力确定方法及应用
CN206634167U (zh) * 2017-01-19 2017-11-14 中山市宏鸿船舶修造有限公司 一种开体驳船
CN109141820A (zh) * 2018-10-08 2019-01-04 中国船舶科学研究中心（中国船舶重工集团公司第七0二研究所） 船体模型扭矩及剪力测量方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHENFENG LI 等: "Ultimate Bearing Capacity Assessment of Hull Girder With Asymmetric Cross Section", 《J. OFFSHORE MECH. ARCT. ENG.》 *
付敬杰: "海上风机一体化运输安装船船体结构强度分析及优化", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111241727A (zh) * 2020-01-03 2020-06-05 浙江工业大学 一种利用单跨有限元模型进行豪华邮轮极限强度计算的方法
CN111241727B (zh) * 2020-01-03 2024-05-07 浙江工业大学 一种利用单跨有限元模型进行豪华邮轮极限强度计算的方法
CN112818480A (zh) * 2021-01-18 2021-05-18 浙江工业大学 压扭组合作用下横撑极限强度计算方法

Also Published As

Publication number Publication date
CN110162881B (zh) 2023-05-16

Similar Documents

Publication Publication Date Title
Byklum et al. A semi-analytical model for global buckling and postbuckling analysis of stiffened panels
CN110162881A (zh) 一种弯剪扭组合下船舯结构极限承载能力的确定方法
CN105844033B (zh) 船舶总纵极限强度简化逐步破坏分析方法
Thill et al. Aerodynamic study of corrugated skins for morphing wing applications
CN110516293B (zh) 船舯结构弯剪扭组合作用下极限承载能力的有限元计算方法
CN108487542A (zh) 钢管混凝土束组合构件稳定承载力计算方法
Asadi et al. Buckling behaviour of thin-walled laminated composite beams having open and closed sections subjected to axial and end moment loading
CN110362940A (zh) 复杂载荷作用下海洋工程结构极限承载能力计算方法
Tiomkin et al. On the stability of two-dimensional membrane wings
Kwon et al. Analyses of cracks in thick stiffened plates repaired with single-sided composite patch
CN106599506B (zh) 一种悬臂横梁对甲板纵骨支撑刚度的确定方法
CN111753357B (zh) 一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法
CN106940736A (zh) 一种考虑侧向载荷作用的纵骨多跨失稳的极限载荷确定方法
Chapman Reduction of profile drag at supersonic velocities by the use of airfoil sections having a blunt trailing edge
Daskovsky The hydrofoil in surface proximity, theory and experiment
CN107066658A (zh) 一种悬臂横梁支撑的甲板纵骨梁柱屈曲的载荷‑端缩曲线确定方法
Arena et al. Nonlinear aeroelastic formulation for flexible high-aspect ratio wings via geometrically exact approach
Kianian et al. Buckling analysis of functionally graded plates based on two-variable refined plate theory using the bubble finite strip method
Oguaghamba et al. Buckling and Postbuckling Loads Characteristics of All Edges Clamped Thin Rectangular Plate
Ozioko et al. Algorithm for exact solution of thick anisotropic plates
Caponnetto et al. Sailing yacht design using advanced numerical flow techniques
Jiao et al. CFD-FEM simulation of water entry of aluminium flat stiffened plate structure considering the effects of hydroelasticity
Datta et al. Modelling a spade rudder as a hollow two-way tapered Kirchhoff’s plate: free dry and wet vibration study with numerical verification
CN220932590U (zh) 复合材料工型梁弯剪稳定性试验加载装置
Senjanović et al. INFLUENCE OF SHEAR ON THE TORSION OF THIN-WALLED GIRDERS.

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