CN111553019A - 一种回转体撞击船体板架横梁剩余速度与转角计算方法 - Google Patents
一种回转体撞击船体板架横梁剩余速度与转角计算方法 Download PDFInfo
- Publication number
- CN111553019A CN111553019A CN202010313602.8A CN202010313602A CN111553019A CN 111553019 A CN111553019 A CN 111553019A CN 202010313602 A CN202010313602 A CN 202010313602A CN 111553019 A CN111553019 A CN 111553019A
- Authority
- CN
- China
- Prior art keywords
- revolving body
- plate frame
- calculating
- motion
- corner
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000003116 impacting effect Effects 0.000 title claims description 12
- 238000006073 displacement reaction Methods 0.000 claims abstract description 19
- 238000005452 bending Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 3
- 238000013277 forecasting method Methods 0.000 abstract description 3
- 210000000988 bone and bone Anatomy 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction 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/15—Vehicle, aircraft or watercraft design
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供一种回转体撞击船体板架横梁剩余速度与转角计算方法,首先确定回转体撞击横梁之前的初始运动物理量与横梁力学参数;然后通过编程求解方程计算横梁破坏时刻塑性铰位置ξm;计算横梁破坏时刻tm和破坏前的最大位移Wm;最后计算横梁破坏后回转体转角αm、剩余速度Vq以及运动方向角βm;本发明可计算光滑头部回转体撞击船体板架横梁剩余速度、转角以及运动方向角,对于高速回转体撞击舰船板架结构的运动提供一种快速预报方法。
Description
技术领域
本发明涉及一种回转体撞击船体板架横梁剩余速度与转角计算方法,属于穿甲力学领域。
背景技术
穿甲力学是一门经典学科,对于弹丸撞击板结构或半无限厚介质有着丰富的研究。对于弹丸高速撞击船体运动规律的研究,由于其直接关系到船体防护能力,因此具有重要意义。与均质平板的撞击问题相比,船体板架结构的撞击运动具有特殊之处。船体板架结构是由薄板和横梁、纵骨等型材构成,其主要承力构建为横梁和纵骨,因此撞击运动规律不同于均质平板撞击,横梁和纵骨对弹丸运动速度和偏转会产生较大影响。
目前均质平板撞击剩余速度经验公式和半经验公式较多,但对于船体板架结构撞击运动特性研究较少。宋卫东[宋卫东,宁建国.刚性回转体撞击加筋靶板的力学模型[J].弹道学报,2007(04):47-50.]和姚熊亮[姚熊亮,吴子奇,王治,等.战斗部对舰船靶标撞击毁伤效能研究[J].哈尔滨工程大学学报,2019,40(1):141-145.]等人对于船体板架结构撞击剩余速度进行了理论模型建立与预报,但其工作存在不足之处,一是其理论在均质平板撞击模型基础上对板架结构修正,并未能真正考虑到横梁的变形模式;二是其理论仅能预报剩余速度,对弹丸的偏转没有预报方法,而船体板架横梁对于弹丸的偏转运动十分显著,对弹丸运动有着重要影响。为此,本专利将弹丸视为光滑头部回转体,将船体板架横梁视为刚塑性梁模型,给出回转体撞击船体板架横梁剩余速度与转角计算方法。
发明内容
本发明的目的是为了对光滑头部回转体撞击船体板架横梁后的运动情况提供一种理论预报方法而提供一种回转体撞击船体板架横梁剩余速度与转角计算方法。
本发明的目的是这样实现的:步骤如下:
步骤1:确定回转体撞击横梁之前的初始运动物理量与横梁力学参数;
步骤2、得到横梁破坏时刻塑性铰位置ξm;
步骤3、计算横梁破坏时刻tm和破坏前的最大位移Wm;
步骤4、计算横梁破坏后回转体转角αm、剩余速度Vq以及运动方向角βm,完成回转体撞击船体板架横梁后的运动情况的理论预报。
本发明还包括这样一些结构特征:
1.步骤1中所述的回转体撞击横梁之前的初始运动物理量包括回转体质量M、转动惯量Jc、头部到质心距离Lc、初始速度V0、初始运动方向与横梁高度方向夹角β0、回转体轴线与初始运动方向夹角α0;横梁力学参数包括单位长度质量m、极限弯矩M0。
2.步骤2中横梁破坏时刻塑性铰位置ξm计算方程为:
3.步骤3中横梁破坏时刻tm和破坏前的最大位移Wm计算公式为:
4.步骤4中横梁破坏后回转体转角αm、剩余速度Vq以及运动方向角βm的计算公式分别为:
αm=-AWm+AtmV0 cosβ0+α0
与现有技术相比,本发明的有益效果是:本发明提出的一种光滑头部回转体撞击船体板架横梁剩余速度与转角计算方法,根据回转体运动初始参数与横梁结构参数,可以理论计算出回转体撞击后剩余速度、转角与运动方向角,对回转体撞击板架结构运动特性进行快速预报。由于目前并无回转体撞击板架结构横梁的转角与运动方向角理论模型,因此本发明填补了这一领域的空白。通过本发明公式,可以对船体板架结构撞击问题进行快速和准确的预报,对于船体受撞击毁伤程度的评估与防护结构设计具有理论意义与工程实用价值。
附图说明
图1是本发明方法的流程图;
图2是光滑回转体撞击船体板架横梁简化模型图;
图3是横梁运动速度场图;
图4是回转体质心位移图;
图5是撞击结束后回转体质心速度图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述。
一种光滑头部回转体撞击船体板架横梁剩余速度与转角计算方法,包括以下步骤:
步骤1、确定回转体撞击横梁之前的初始运动物理量与横梁力学参数;
步骤2、通过编程求解方程计算横梁破坏时刻塑性铰位置ξm;
步骤3、计算横梁破坏时刻tm和破坏前的最大位移Wm;
步骤4、计算横梁破坏后回转体转角αm、剩余速度Vq以及运动方向角βm。
上述方案中,步骤1中所述的回转体撞击横梁之前的初始运动物理量包括回转体质量M、转动惯量Jc、头部到质心距离Lc、初始速度V0、初始运动方向与横梁高度方向夹角β0、回转体轴线与初始运动方向夹角α0;横梁力学参数包括单位长度质量m、极限弯矩M0。
上述方案中,步骤2中横梁破坏时刻塑性铰位置ξm计算方程由回转体运动方程与横梁变形方程推导而得。
首先针对此问题作出如下假设:1)回转体假设为刚体;2)船体板架横梁视为刚塑性梁模型;3)横梁与回转体运动面垂直,即运动为二维的。4)由于撞击过程短,横梁运动效应在撞击过程中不会抵达边界,因此将横梁视为无限长梁,而塑性铰在撞击过程中沿着梁移动。
通过以上假设将问题转化为回转体撞击刚塑性梁问题,简化后的模型如图2所示。由于对称性,回转体在x-y平面作二维平面运动。其中水平方向x为回转体速度方向,板架结构平面法线与水平方向角度为β0。回转体的质心为C。回转体转角α为回转体轴线与水平方向的角度。
在撞击过程中,由于横梁的位移相比板架结构尺寸是微小的,因此假设撞击过程中板架结构角度不变,而弹头受力FN沿靶面法线方向。横梁受力大小为FN,方向与回转体受力反向。定义板架结构法线方向为q方向,如图2所示,横梁中点沿q方向位移为W。横梁运动局部坐标系如图3所示。在回转体作用下,横梁塑性铰不断向边界移动,塑性铰坐标为ξ。假设横梁速度场为线性分布,如图3所示,沿全梁分布的速度场可表达为:
首先建立回转体运动方程。转动力矩与回转体转角关系为:
实际问题中,转角α随时间而改变,即:
α(t)=α0+Δα(t) (4)
其中,Δα(t)为转角随时间变化量。但由于撞击时间较短,Δα(t)相比于β0+α0很小,因此式简化为:
回转体质心C垂直于靶面位移为Wc,平行于靶面位移为Uc,如图4所示。回转体沿q方向运动方程为:
由刚体运动可知,回转体质心位移Wc与梁中点挠度W的关系为:
Wc-LcΔα=W (7)
由式和可得:
进一步代入式,可得:
由梁的受力平衡方程可得:
代入式可得:
在塑性铰处,仅有梁的极限弯矩M0作用,而没有剪力,因此由弯矩平衡方程可得:
对式和进行求解,可得梁挠度公式。文献[Normam.Jones著.结构冲击[M].第二版.许俊,蒋平译.北京:国防工业出版社,2018:72-79.]中对求解过程有详细论述,本文不做赘述,仅引用其结论。梁变形区域挠度公式为:
对于图3所示梁的运动,梁的塑性应变可表示为:
随着中点位移增大,塑性应变增加,直到塑性应变达到断裂应变εm,梁破坏,撞击过程结束。利用式可得梁破坏的条件为:
利用式和,可得到梁破坏时的塑性铰坐标ξm满足的方程:
上述方案中,步骤3中所述的横梁破坏时刻tm和破坏前的最大位移Wm的计算公式建立方法如下:
由文献[Normam.Jones著.结构冲击[M].第二版.许俊,蒋平译.北京:国防工业出版社,2018:72-79.],撞击过程中时间与塑性铰位置关系为:
利用式即可得到梁发生破坏的时间tm:
由式可得撞击结束时刻梁中点位移为:
由式可得撞击结束时刻梁中点速度为:
上述方案中,步骤4中所述的横梁破坏后回转体转角αm、剩余速度Vq以及运动方向角βm的计算公式建立方法如下:
利用式和,得到回转体转角α的方程为:
对式进行积分,并利用回转体运动初始条件:
W(0)=0,α(0)=α0 (26)
可得回转体转角的解为:
α(t)=-AW+AtV0cosβ0+α0 (28)
式中,
撞击结束时回转体转角为:
αm=-AWm+AtmV0 cosβ0+α0 (30)
如图4所示,回转体初始质心为C点,撞击结束后质心为C’点,回转体质心在q-z坐标系下的位移为:
Wcm=Wm+LcΔαm (31)
Ucm=tmV0 sinβ0 (32)
回转体的剩余速度为:
回转体质心C水平与垂直方向位移为:
Ucx=Ucm sinβ0+Wcm cosβ0 (34)
Ucy=Ucm cosβ0-Ucm sinβ0 (35)
下面结合具体数值给出本发明的实施例:
结合图1,本发明的一种光滑头部回转体撞击船体板架横梁剩余速度与转角计算方法,包括以下步骤:
步骤1,确定回转体撞击横梁之前的初始运动物理量与横梁力学参数;
步骤1中所述的回转体撞击横梁之前的初始运动物理量包括回转体质量M、转动惯量Jc、头部到质心距离Lc、初始速度V0、初始运动方向与横梁高度方向夹角β0、回转体轴线与初始运动方向夹角α0;横梁力学参数包括单位长度质量m、极限弯矩M0。
本实施例中回转体与横梁力学参数如下表1所示。
表1回转体与横梁力学参数
步骤2,通过编程求解方程计算横梁破坏时刻塑性铰位置ξm;
求解方程式,得到横梁破坏时刻塑性铰位置ξm=0.17m。
步骤3,计算横梁破坏时刻tm和破坏前的最大位移Wm;
利用式和计算横梁破坏时刻tm和破坏前的最大位移Wm为:
步骤4,计算横梁破坏后回转体转角αm、剩余速度Vq以及运动方向角βm;
利用式、、计算横梁破坏后回转体转角αm、剩余速度Vq以及运动方向角βm为:
αm=-AWm+AtmV0 cosβ0+α0=3.30°,
至此,已求解了光滑头部回转体撞击横梁后剩余速度与转角等运动参数。
综上,本发明提出了一种光滑头部回转体撞击船体板架横梁剩余速度与转角计算方法,首先确定回转体撞击横梁之前的初始运动物理量与横梁力学参数;然后通过编程求解方程计算横梁破坏时刻塑性铰位置ξm;计算横梁破坏时刻tm和破坏前的最大位移Wm;最后计算横梁破坏后回转体转角αm、剩余速度Vq以及运动方向角βm;本发明可计算光滑头部回转体撞击船体板架横梁剩余速度、转角以及运动方向角,对于高速回转体撞击舰船板架结构的运动提供一种快速预报方法。
Claims (5)
1.一种回转体撞击船体板架横梁剩余速度与转角计算方法,其特征在于:步骤如下:
步骤1:确定回转体撞击横梁之前的初始运动物理量与横梁力学参数;
步骤2、得到横梁破坏时刻塑性铰位置ξm;
步骤3、计算横梁破坏时刻tm和破坏前的最大位移Wm;
步骤4、计算横梁破坏后回转体转角αm、剩余速度Vq以及运动方向角βm,完成回转体撞击船体板架横梁后的运动情况的理论预报。
2.根据权利要求1所述的一种回转体撞击船体板架横梁剩余速度与转角计算方法,其特征在于:步骤1中所述的回转体撞击横梁之前的初始运动物理量包括回转体质量M、转动惯量Jc、头部到质心距离Lc、初始速度V0、初始运动方向与横梁高度方向夹角β0、回转体轴线与初始运动方向夹角α0;横梁力学参数包括单位长度质量m、极限弯矩M0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010313602.8A CN111553019B (zh) | 2020-04-20 | 2020-04-20 | 一种回转体撞击船体板架横梁剩余速度与转角计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010313602.8A CN111553019B (zh) | 2020-04-20 | 2020-04-20 | 一种回转体撞击船体板架横梁剩余速度与转角计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111553019A true CN111553019A (zh) | 2020-08-18 |
CN111553019B CN111553019B (zh) | 2023-01-03 |
Family
ID=71999204
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010313602.8A Active CN111553019B (zh) | 2020-04-20 | 2020-04-20 | 一种回转体撞击船体板架横梁剩余速度与转角计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111553019B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114297860A (zh) * | 2021-12-30 | 2022-04-08 | 中国人民解放军军事科学院国防工程研究院 | 延时引信弹药与跳弹式防护结构碰撞分析方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103308273A (zh) * | 2013-05-28 | 2013-09-18 | 哈尔滨工程大学 | 大型封闭结构受风载作用刚度强度测试方法 |
CN104008226A (zh) * | 2014-04-17 | 2014-08-27 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种尾钩与道面撞击后反弹产生的纵向阻尼力计算方法 |
EP2793043A1 (en) * | 2013-04-18 | 2014-10-22 | Airbus Defence and Space GmbH | Determination of weapon locations and projectile trajectories by using automatic and hybrid processing of acoustic and electromagnetic detections |
CN104657611A (zh) * | 2015-02-14 | 2015-05-27 | 中国人民解放军海军工程大学 | 水下爆炸作用下舰船整体弹塑性运动响应预测方法及系统 |
CN105844033A (zh) * | 2016-03-30 | 2016-08-10 | 中国舰船研究设计中心 | 船舶总纵极限强度简化逐步破坏分析方法 |
CN106596028A (zh) * | 2016-12-23 | 2017-04-26 | 中国船舶工业集团公司第七0八研究所 | 一种悬臂横梁支撑的甲板纵骨梁柱屈曲临界载荷的确定方法 |
CN106599506A (zh) * | 2016-12-23 | 2017-04-26 | 中国船舶工业集团公司第七0八研究所 | 一种悬臂横梁对甲板纵骨支撑刚度的确定方法 |
CN110532656A (zh) * | 2019-08-20 | 2019-12-03 | 哈尔滨工程大学 | 一种基于侵彻毁伤相似性的金属板架结构等效设计方法 |
CN110852012A (zh) * | 2019-11-11 | 2020-02-28 | 武汉理工大学 | 轴压加筋板极限强度预测方法 |
-
2020
- 2020-04-20 CN CN202010313602.8A patent/CN111553019B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2793043A1 (en) * | 2013-04-18 | 2014-10-22 | Airbus Defence and Space GmbH | Determination of weapon locations and projectile trajectories by using automatic and hybrid processing of acoustic and electromagnetic detections |
CN103308273A (zh) * | 2013-05-28 | 2013-09-18 | 哈尔滨工程大学 | 大型封闭结构受风载作用刚度强度测试方法 |
CN104008226A (zh) * | 2014-04-17 | 2014-08-27 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种尾钩与道面撞击后反弹产生的纵向阻尼力计算方法 |
CN104657611A (zh) * | 2015-02-14 | 2015-05-27 | 中国人民解放军海军工程大学 | 水下爆炸作用下舰船整体弹塑性运动响应预测方法及系统 |
CN105844033A (zh) * | 2016-03-30 | 2016-08-10 | 中国舰船研究设计中心 | 船舶总纵极限强度简化逐步破坏分析方法 |
CN106596028A (zh) * | 2016-12-23 | 2017-04-26 | 中国船舶工业集团公司第七0八研究所 | 一种悬臂横梁支撑的甲板纵骨梁柱屈曲临界载荷的确定方法 |
CN106599506A (zh) * | 2016-12-23 | 2017-04-26 | 中国船舶工业集团公司第七0八研究所 | 一种悬臂横梁对甲板纵骨支撑刚度的确定方法 |
CN110532656A (zh) * | 2019-08-20 | 2019-12-03 | 哈尔滨工程大学 | 一种基于侵彻毁伤相似性的金属板架结构等效设计方法 |
CN110852012A (zh) * | 2019-11-11 | 2020-02-28 | 武汉理工大学 | 轴压加筋板极限强度预测方法 |
Non-Patent Citations (3)
Title |
---|
R.E.OSHIRO: "Experimental Study of Collision in Scaled Naval Structures", 《IMPACT ENGINEERING》 * |
姚熊亮: "战斗部为舰船靶标侵彻毁伤效能研究", 《哈尔滨工程大学学报》 * |
张楷弘: "大外飘船舶非线性波浪载荷时域水弹性分析方法与应用研究", 《CNKI博士电子期刊》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114297860A (zh) * | 2021-12-30 | 2022-04-08 | 中国人民解放军军事科学院国防工程研究院 | 延时引信弹药与跳弹式防护结构碰撞分析方法 |
CN114297860B (zh) * | 2021-12-30 | 2022-11-22 | 中国人民解放军军事科学院国防工程研究院 | 延时引信弹药与跳弹式防护结构碰撞分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111553019B (zh) | 2023-01-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111553019B (zh) | 一种回转体撞击船体板架横梁剩余速度与转角计算方法 | |
CN109344452A (zh) | 一种目标易损性分析方法 | |
CN101289733A (zh) | 激光辅助预应力喷丸成形-强化复合方法 | |
Zhang et al. | Research on impact resistance of ceramic matrix composites | |
CN114136779A (zh) | 一种准脆性材料i-ii型断裂韧度试验求解方法 | |
Liu et al. | A study on the surface overpressure distribution and formation of a double curvature liner under a two-point initiation | |
Karmakar et al. | Transient dynamic response of delaminated composite rotating shallow shells subjected to impact | |
Yang et al. | Load and wear experiments on the impact hammer of a vertical shaft impact crusher | |
Bell | The dynamic buckling of rods at large plastic strain | |
Jena et al. | Simulation of bullet penetration using finite element method | |
Rosenberg et al. | More on the ricochet of eroding long rods—validating the analytical model with 3D simulations | |
Zhang et al. | Study on the dynamic fatigue performance of shot blasting machine’s blades based on modal analysis | |
Angel et al. | Numerical modeling of impact interaction between an active element and a deformable barrier | |
Steverding et al. | The propagation law of cleavage fracture | |
Jiang et al. | High rotating speed projectile penetrating into moving vehicle door at different incident angle | |
CN104166402B (zh) | 单滑块变质心控制飞行器平衡运动状态的近似解析方法 | |
Stépán et al. | Impact induced internal fatigue cracks | |
LU102013B1 (en) | The P-I Curve Evaluation Method of Reinforced Concrete Structural Components Based on Equivalent Single Degree of Freedom | |
Xing et al. | A preliminary study of using high-speed Digital Image Correlation (DIC) to characterize the penetration on geomaterial | |
Fu et al. | Laser forming of the panel with crossed reinforcing bars | |
Chun et al. | Behavior of uniform anisotropic beams of rectangular section under transverse impact of a mass | |
Gerasimov et al. | Numerical simulation of high-velocity projectile interactions with groups of spaced rods and plates | |
Rumyantsev | Kinetics of cumulative jet penetration into glass | |
Titov et al. | Structural and Imitation Model of a Concrete Mix Rotational Seal | |
McKee et al. | Design of masonry walls for blast loading |
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 |