CN112651156B - 一种高速来流条件下激光毁伤效应的多场耦合分析方法 - Google Patents
一种高速来流条件下激光毁伤效应的多场耦合分析方法 Download PDFInfo
- Publication number
- CN112651156B CN112651156B CN202011518700.1A CN202011518700A CN112651156B CN 112651156 B CN112651156 B CN 112651156B CN 202011518700 A CN202011518700 A CN 202011518700A CN 112651156 B CN112651156 B CN 112651156B
- Authority
- CN
- China
- Prior art keywords
- structural
- coupling surface
- coupling
- ablation
- fluid
- 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/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Laser Beam Processing (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供一种高速来流条件下激光毁伤效应的多场耦合分析方法,包括如下步骤:预先建立包含激光的热力载荷下结构分析模型和高速来流作用下流体分析模型;根据结构瞬态热力耦合分析确定结构热力响应,获得耦合面的烧蚀退化量;通过网格重构更新结构分析模型和流体分析模型;基于耦合面温度场,在耦合界面处进行数据映射,得到固体耦合面气动力、气动热,然后更新结构分析模型的载荷条件;基于更新的结构分析模型进行下一时刻瞬态热力耦合分析;依此循环,最终得到高速来流条件下激光对飞行目标毁伤效应的全耦合数值计算结果。本发明提出考虑局部烧蚀效应的多场耦合数值分析策略,实现高速来流条件下激光对飞行目标毁伤效应的全耦合数值计算。
Description
技术领域
本发明涉及热力学结构领域,特别是涉及一种对高速来流条件下激光毁伤效应的多场耦合分析方法。
背景技术
高速飞行导弹具有飞行速度快、毁伤威力大、突防能力强等战略优势,现有动能拦截导弹防御系统难以百分之百奏效。高能激光武器具有即瞄即打、持续打击、转移迅速等特点,在反制高速飞行靶目标方面具有独特的优势。
激光对高速靶目标的拦截本质上是高速来流条件下的激光热力毁伤问题。然而,目前对该问题的研究还停留在原理性和概念性阶段,主要的原因在于缺乏定量化的分析与预报能力,导致对其中的热力破坏机理认知不明。
高速来流条件下的激光热力毁伤效应数值分析涉及到传热、高速流动、结构响应、热力烧蚀等多种效应之间的相互作用,是一类含有移动边界的复杂流-热-固多物理场耦合问题,相关的研究目前还很少见。现有技术中,气流的影响往往以等效、固定的常数描述,没有考虑激光诱导的烧蚀行为所带来的移动边界、壁面高温等效应对高速流动的反作用。这种情况下,无法对激光热力烧蚀行为给出准确的描述,更不论对其中的热力破坏机理进行全面深入的理解。
因此,必须发展适合于高速来流条件下激光毁伤效应的多场耦合分析方法,定量化地给出目标特性、打击距离、激光参数、毁伤模式、毁伤时间等众多因素之间的关系,从而明确激光对高速靶目标的毁伤机理,实现高效激光打击策略的制定。
发明内容
本发明的目的是提供一种精确且全面对高速来流条件下激光毁伤效应的多场耦合分析方法
具体地,本发明提供一种高速来流条件下激光毁伤效应的多场耦合分析方法,包括如下步骤:
步骤100,预先建立包含激光的热力载荷下结构分析模型和高速来流作用下流体分析模型;
步骤200,根据结构瞬态热力耦合分析确定结构整体温度、结构耦合面温度场和结构耦合面变形量的热力响应,同时基于结构整体温度、结构耦合面气动力和结构耦合面气动热,结合结构烧蚀模型进行结构烧蚀分析,获得结构耦合面的烧蚀退化量;
步骤300,基于结构耦合面变形量和烧蚀退化量,相加得到结构耦合面总变形量,然后进行结构网格重构,更新结构分析模型;同时基于结构耦合面总变形量在耦合界面处进行数据映射,获得流体耦合面总变形量,然后进行流体网格重构,更新流体分析模型;
步骤400,基于耦合面温度场,在耦合界面处进行数据映射,得到流体耦合面温度场,再次更新流体分析模型并进行流体稳态分析,获得流体耦合面气动力和流体耦合面气动热的流场信息;基于流体耦合面气动力和气动热,在耦合界面处进行数据映射,得到固体耦合面气动力、气动热,然后更新结构分析模型的载荷条件;
步骤500,基于更新的结构分析模型,进行瞬态热力耦合分析,时间间隔为Δt;判断是否满足预设的停止条件,如果是,则停止整个流程,否则,返回执行步骤200,最终得到高速来流条件下局部热流对快速飞行目标毁伤效应的全耦合数值计算结果。
本发明克服现有技术中对流-热-固多场耦合效应估计不足的问题,提出的多场耦合分析方法,通过构建局部烧蚀效应的多场耦合数值分析策略,从而实现高速来流条件下激光对快速飞行目标毁伤效应的全耦合数值计算,为明确高速来流条件下激光热力毁伤机理提供有效的技术支撑。
本发明能够定量化地给出目标特性、打击距离、激光参数、毁伤模式、毁伤时间等众多因素之间的关系,从而明确激光对高速靶目标的毁伤机理,实现高效激光打击策略的制定,为未来高技术战争提供基础性和技术性支撑。
附图说明
图1是本发明一个实施方式的方法流程图;
图2是本发明一个实施方式中预先建立的结构分析模型和流场分析模型;
图3是本发明实施例中不同时刻流体网格与固体网格示意图;
图4是本发明实施例中不同时刻结构对称面烧蚀形貌;
图5是本发明实施例中不同激光辐照时刻光斑中心烧蚀速率与流场参数对应关系;
图6是本发明实施例中激光辐照4s时结构对称面烧蚀坑型线试验与仿真结果对比。
具体实施方式
以下通过具体实施例和附图对本方案的具体结构和实施过程进行详细说明。
如图1所示,在本发明的一个实施方式中,公开一种高速来流条件下激光毁伤效应的多场耦合分析方法,包括如下步骤:
步骤100,预先建立包含激光的热力载荷下结构分析模型和高速来流作用下流体分析模型;
这里的热力载荷还包括气动力和气动热。其中激光以空间热流的形式作用在结构的表面。结构分析模型通过有限元分析方法可得到结构热力响应,流体分析模型通过有限体积方法可得到流场特性。
步骤200,根据结构瞬态热力耦合分析确定结构整体温度、结构耦合面温度场和结构耦合面变形量的热力响应,同时基于结构整体温度、结构耦合面气动力和结构耦合面气动热,结合结构烧蚀模型进行结构烧蚀分析,获得结构耦合面的烧蚀退化量;
在本步骤中,结构烧蚀模型可通过结构耦合面法向线烧蚀速率v进行表征。以下述的C/SiC复合材料薄板的烧蚀模型为例,其对应的表达式为:
v=vTC+vTM
其中,vTC为热化学烧蚀线速率,vTM为机械剥蚀线速率。基于结构耦合面法向线烧蚀速率v,对其在分析时间步内积分得到结构耦合面烧蚀退化量。
步骤300,基于结构耦合面变形量和烧蚀退化量,相加得到结构耦合面总变形量,然后进行结构网格重构,更新结构分析模型;同时基于结构耦合面总变形量在耦合界面处进行数据映射,获得流体耦合面总变形量,然后进行流体网格重构,更新流体分析模型;
在本步骤中,网格重构通是过插值方法实现。将结构耦合面总变形量作为待插值量,运用插值函数对结构耦合面网格节点的总变形量进行拟合,然后利用构造出来的插值函数将结构耦合面变形效应光滑地分散到整个结构网格区域的节点上,从而得到结构所有节点的总变形量。结合结构节点坐标和节点变形量,可得到结构节点的变形后坐标,从而实现结构网格重构。该网格重构方法能够保证耦合界面的高精度捕捉与重构。
本步骤中的所有的网格重构方式相同。
步骤400,基于耦合面温度场,在耦合界面处进行数据映射,得到流体耦合面温度场,再次更新流体分析模型并进行流体稳态分析,获得流体耦合面气动力和流体耦合面气动热的流场信息;基于流体耦合面气动力和气动热,在耦合界面处进行数据映射,得到固体耦合面气动力、气动热,然后更新结构分析模型的载荷条件。该步骤实现了流体和结构在耦合界面的信息传递。
基于结构耦合面总变形的数据映射是通过插值方法将总变形量从结构耦合面传递到流体耦合面。
基于流体耦合面气动力和气动热的数据映射是通过插值方法将气动力和气动热从流体耦合面传递到结构耦合面。
步骤500,基于更新的结构分析模型,进行瞬态热力耦合分析,时间间隔为Δt;判断是否满足预设的停止条件,如果是,则停止整个流程,否则,返回执行步骤200,最终得到高速来流条件下局部热流对快速飞行目标毁伤效应的全耦合数值计算结果。
时间间隔Δt根据实际需求预先设置;停止条件是以当前时间时刻为最后的计算时刻。
本实施方式克服现有技术中对流-热-固多场耦合效应估计不足的问题,提出的多场耦合分析方法,通过构建局部烧蚀效应的多场耦合数值分析策略,从而实现高速来流条件下激光对快速飞行目标毁伤效应的全耦合数值计算,为明确高速来流条件下激光热力毁伤机理提供有效的技术支撑。
本实施方式能够定量化地给出目标特性、打击距离、激光参数、毁伤模式、毁伤时间等众多因素之间的关系,从而明确激光对高速靶目标的毁伤机理,实现高效激光打击策略的制定,为未来高技术战争提供基础性和技术性支撑。
以下以具体实施例对本方法的工作过程作详细说明。
如图2所示,在Ma=1.8,h=25km的超声速气流条件下,对C/SiC复合材料薄板的激光毁伤效应进行数值仿真分析;来流条件中马赫为1.8,攻角0,高度25km;结构条件为,材料种类为C/SiC复合材料,薄板长、宽、高分别为50mm,50mm,3mm;激光载荷为,功率2kW,高斯分布,作用于薄板中心,半径为5mm,辐照时间为4s,时间间隔Δt取为0.1s。
在前述基础上利用本方法,对高速气流条件下C/SiC复合材料薄板的激光毁伤效应进行分析,得到了相应的毁伤结果,其中图3为不同时刻流体网格固体网格示意图,体现了本实施例中的高精度网格重构过程;图4为不同时刻结构对称面烧蚀形貌;图5为不同时刻光斑中心烧蚀速率与流场参数的关系;图6为激光辐照4s时,结构对称面烧蚀坑型线试验与仿真结果对比,从而验证了本方法的有效性。
至此,本领域技术人员应认识到,虽然本文已详尽示出和描述了本发明的多个示例性实施例,但是,在不脱离本发明精神和范围的情况下,仍可根据本发明公开的内容直接确定或推导出符合本发明原理的许多其他变型或修改。因此,本发明的范围应被理解和认定为覆盖了所有这些其他变型或修改。
Claims (4)
1.一种高速来流条件下激光毁伤效应的多场耦合分析方法,其特征在于,包括如下步骤:
步骤100,预先建立包含激光的热力载荷下结构分析模型和高速来流作用下流体分析模型,所述热力载荷还包括气动力和气动热,所述激光以空间热流的形式作用在结构的表面;
步骤200,根据结构瞬态热力耦合分析确定结构整体温度、结构耦合面温度场和结构耦合面变形量的热力响应,同时基于结构整体温度、结构耦合面气动力和结构耦合面气动热,结合结构烧蚀模型进行结构烧蚀分析,获得结构耦合面的烧蚀退化量;
步骤300,基于结构耦合面变形量和烧蚀退化量,相加得到结构耦合面总变形量,然后进行结构网格重构,更新结构分析模型;同时基于结构耦合面总变形量在耦合界面处进行数据映射,获得流体耦合面总变形量,然后进行流体网格重构,更新流体分析模型;
所述网格重构通是过插值方法实现,将结构耦合面总变形量作为待插值量,运用插值函数对结构耦合面网格节点的总变形量进行拟合,然后利用构造出来的插值函数将结构耦合面变形效应光滑地分散到整个结构网格区域的节点上,从而得到结构所有节点的总变形量;
结合结构节点坐标和节点变形量,可得到结构节点的变形后坐标,从而实现结构网格重构;
基于结构耦合面总变形的数据映射是通过插值方法将总变形量从结构耦合面传递到流体耦合面;
步骤400,基于耦合面温度场,在耦合界面处进行数据映射,得到流体耦合面温度场,再次更新流体分析模型并进行流体稳态分析,获得流体耦合面气动力和流体耦合面气动热的流场信息;基于流体耦合面气动力和气动热,在耦合界面处进行数据映射,得到固体耦合面气动力、气动热,然后更新结构分析模型的载荷条件;基于流体耦合面气动力和气动热的数据映射是通过插值方法将气动力和气动热从流体耦合面传递到结构耦合面;
步骤500,基于更新的结构分析模型,进行瞬态热力耦合分析,时间间隔为Δt;判断是否满足预设的停止条件,如果是,则停止整个流程,否则,返回执行步骤 200,最终得到高速来流条件下局部热流对快速飞行目标毁伤效应的流动-传热-结构-烧蚀全耦合数值计算结果。
2.根据权利要求1所述的多场耦合分析方法,其特征在于,
所述步骤300中,网格重构是在保证外形的情况下,网格节点数目不变,位置重新分布。
3.根据权利要求1所述的多场耦合分析方法,其特征在于,
所述步骤500中,时间间隔Δt根据实际需求预先设置。
4.根据权利要求1所述的多场耦合分析方法,其特征在于,
所述步骤500中,停止条件是以当前时间时刻为最后的计算时刻。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011518700.1A CN112651156B (zh) | 2020-12-21 | 2020-12-21 | 一种高速来流条件下激光毁伤效应的多场耦合分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011518700.1A CN112651156B (zh) | 2020-12-21 | 2020-12-21 | 一种高速来流条件下激光毁伤效应的多场耦合分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112651156A CN112651156A (zh) | 2021-04-13 |
CN112651156B true CN112651156B (zh) | 2022-03-04 |
Family
ID=75358543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011518700.1A Active CN112651156B (zh) | 2020-12-21 | 2020-12-21 | 一种高速来流条件下激光毁伤效应的多场耦合分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112651156B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114417747B (zh) * | 2022-01-19 | 2022-12-20 | 中山大学 | 一种水上飞机的水面起降状态预测方法及装置 |
CN114626269B (zh) * | 2022-03-21 | 2023-03-24 | 中国科学院力学研究所 | 一种分析多材料结构激光烧蚀效应的方法和装置 |
CN115825061B (zh) * | 2022-11-21 | 2023-07-07 | 中国科学院力学研究所 | 一种高速气流下激光诱导超高温烧蚀的原位观测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095603A (zh) * | 2015-09-09 | 2015-11-25 | 哈尔滨工业大学 | 一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法 |
CN105973563A (zh) * | 2016-05-12 | 2016-09-28 | 中国科学院力学研究所 | 一种强激光与超声速风洞联合破坏效应的试验系统及方法 |
CN107368661A (zh) * | 2017-07-27 | 2017-11-21 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种高超声速飞行器热气动弹性特性的耦合计算方法 |
CN111458366A (zh) * | 2020-04-17 | 2020-07-28 | 北京空天技术研究所 | 一种烧蚀热防护系统结构气动热/传热耦合分析方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109918808B (zh) * | 2019-03-13 | 2023-10-20 | 北京宇航系统工程研究所 | 一种气热弹三场耦合仿真分析方法 |
CN110991097B (zh) * | 2019-08-06 | 2023-05-23 | 电子科技大学 | 一种激光烧蚀碳纤维靶材温度场的计算方法 |
-
2020
- 2020-12-21 CN CN202011518700.1A patent/CN112651156B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095603A (zh) * | 2015-09-09 | 2015-11-25 | 哈尔滨工业大学 | 一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法 |
CN105973563A (zh) * | 2016-05-12 | 2016-09-28 | 中国科学院力学研究所 | 一种强激光与超声速风洞联合破坏效应的试验系统及方法 |
CN107368661A (zh) * | 2017-07-27 | 2017-11-21 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种高超声速飞行器热气动弹性特性的耦合计算方法 |
CN111458366A (zh) * | 2020-04-17 | 2020-07-28 | 北京空天技术研究所 | 一种烧蚀热防护系统结构气动热/传热耦合分析方法 |
Non-Patent Citations (1)
Title |
---|
强激光辐照典型目标多场耦合效应的实验与数值模拟研究;黄亿辉;《中国科学院力学研究所机构知识库》;20150828;第1-128页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112651156A (zh) | 2021-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112651156B (zh) | 一种高速来流条件下激光毁伤效应的多场耦合分析方法 | |
Hebert | Metallurgical aspects of powder bed metal additive manufacturing | |
Linke et al. | Performance of different tungsten grades under transient thermal loads | |
Deepanraj et al. | Theoretical analysis of gas turbine blade by finite element method | |
CN111627503B (zh) | 一种氧化铝陶瓷基复合材料涂层激光熔覆制造过程中应力场的预测方法 | |
CN112329201B (zh) | 一种集成化离子电推进仿真模型及方法 | |
CN101799336B (zh) | 涡轮叶片热障涂层应力的预测方法 | |
Yuan et al. | Impact point prediction for thrusting projectiles in the presence of wind | |
CN112329247A (zh) | 一种霍尔电推力器放电等离子体和自溅射的仿真方法 | |
CN114441590B (zh) | 一种针对梯度防热材料传热及力学特性的确定方法及系统 | |
CN115169017A (zh) | 一种用于激光-电磁复合推力器宏观性能参数的预测方法 | |
Yakinthos et al. | Modelling operation of system of recuperative heat exchangers for aero engine with combined use of porosity model and thermo-mechanical model | |
Narahari et al. | MHD natural convection flow past an impulsively started infinite vertical porous plate with Newtonian heating in the presence of radiation | |
CN105403502A (zh) | 一种高温合金涡轮叶片服役温度的实验评估方法 | |
Long et al. | Research on the heat insulation performance of advanced thermal barrier coatings based on dense vertical cracked structure | |
Zhou et al. | An axisymmetric hybrid model for the plasma transport in a helicon plasma thruster | |
Pavei et al. | Thermo-mechanical design of the Plasma Driver Plate for the MITICA ion source | |
Gonzalez et al. | Simulation of plume-plasma expansion with one-dimensional Particle-in-Cell | |
Sadowski et al. | Numerical and experimental analysis of foreign objects impact into the surface with TBC coating | |
CN112651150B (zh) | 一种Al2O3/316L不锈钢阻氚系统表面多重裂纹行为的扩展有限元模拟分析方法 | |
Reilly et al. | Thermal analysis of the 100-kW class X3 hall thruster | |
Azevedo et al. | Magnetic field enhancement of the quad confinement thruster (QCT): Design and early development of the QCT phoenix | |
Brezinsek et al. | Preparation of PFCs for the efficient use in ITER and DEMO plasma-wall interaction studies within the EUROfusion consortium | |
Karthik et al. | Material Selection and Structural Integration of Actively Cooled High Speed Combustion Chamber | |
Lee et al. | Numerical Study of Ablation Phenomena of Flame Deflector |
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 |