CN113673127B - 一种复合材料激光烧蚀形貌快速获取方法 - Google Patents
一种复合材料激光烧蚀形貌快速获取方法 Download PDFInfo
- Publication number
- CN113673127B CN113673127B CN202110955858.3A CN202110955858A CN113673127B CN 113673127 B CN113673127 B CN 113673127B CN 202110955858 A CN202110955858 A CN 202110955858A CN 113673127 B CN113673127 B CN 113673127B
- Authority
- CN
- China
- Prior art keywords
- ablation
- composite material
- grid
- time
- laser
- 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
- 239000002131 composite material Substances 0.000 title claims abstract description 79
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000000608 laser ablation Methods 0.000 title claims abstract description 29
- 238000002679 ablation Methods 0.000 claims abstract description 100
- 238000004364 calculation method Methods 0.000 claims abstract description 41
- 230000003068 static effect Effects 0.000 claims abstract description 15
- 238000002474 experimental method Methods 0.000 claims abstract description 7
- 239000000463 material Substances 0.000 claims description 33
- 238000009826 distribution Methods 0.000 claims description 18
- 230000004044 response Effects 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 4
- 238000005530 etching Methods 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 description 10
- 238000011156 evaluation Methods 0.000 description 9
- 238000004088 simulation Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000007123 defense Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 231100000225 lethality Toxicity 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Laser Beam Processing (AREA)
Abstract
本发明公开了一种复合材料激光烧蚀形貌快速获取方法,该方法的主要实施步骤为:1、分析得到参数空间;2、通过常规实验方法获取多组复合材料线烧蚀速率和烧蚀起始时间数据;3、获取复合材料线烧蚀速率和烧蚀起始时间拟合公式;4、求取拟合公式中的拟合参数;5、构建复合材料三维空间表面网格模型,确定出网格模型的表面激光功率密度、表面气流速度和气流静压;6:启动复合材料激光烧蚀形貌快速获取程序的计算流程,按照时间步迭代循环计算,并重新绘制三维空间网格。使用该方法能过快速的、准确的获取复合材料激光烧蚀之后的形貌特性。
Description
技术领域
本发明涉及一种复合材料激光烧蚀形貌快速获取方法,特别涉及一种复合材料层合板及复合材料三明治夹层结构材料的激光烧蚀形貌快速获取方法。
背景技术
复合材料强度高、重量轻,在工业制造、航空航天、国防装备等各领域被大量使用。由于激光加工、激光防卫等需求,激光与复合材料相互作用过程一直受到人们的关注。在作用过程中,为了解激光辐照复合材料的烧蚀响应,科研人员采用了实验研究、数值模拟等方法,针对材料温度、烧蚀热、烧蚀失重、烧蚀穿孔时间、烧蚀孔尺寸/面积等烧蚀响应量开展研究。
例如,美国2007年版的《High energy laser(HEL)lethality data collectionstandards》中明确规定了激光毁伤研究中,穿孔时间为应采集的目标响应量之一。在一些需要进行快速评估的应用场景中,人们常常对材料烧蚀穿孔时间、烧蚀孔尺寸/面积等典型效应表征量开展研究,获得工程经验公式,为快速评估提供输入。然而,对多种快速评估应用场景而言,烧蚀穿孔时间、烧蚀孔尺寸/面积等效应表征量提供的信息有限,将其用于快速评估时误差较大。要提高烧蚀后效分析评估的精度,复合材料制备的结构/器件本身的实际烧蚀形貌(包括烧蚀孔实际形态)是必要的输入条件。
数值模拟方法因其经济性、便捷性、产生数据的丰富性等因素被广泛用于计算材料烧蚀形貌。在采用数值模拟方法进行计算时,科研人员通常考虑加入对各种复杂过程的模拟,以提高模型保真度和精度,丰富输出数据类型。但是在实际使用中,数值模拟方法却存在着以下问题:计算过程需耗费大量计算资源,且耗时较长,特别是考虑了复杂物理过程的计算模型,根据计算模型的规模一个算例在主流工作站上计算耗时长达数十分钟甚至几个小时。这对需要快速评估的应用场景而言是远远不能接受的。
而通过常规实验方法预先获得不同参数条件下激光辐照复合材料的烧蚀形貌,实验次数多,导致不仅费用高,而且获得的结果很难插值得到当前激光参数、气流参数条件的实际烧蚀形貌。
因而,目前复合材料激光烧蚀形貌尚缺乏快速有效的评估手段。
发明内容
为了解决采用数值仿真的方式获取复合材料激光烧蚀形貌时存在计算量大,耗时长,无法满足快速评估的问题,以及常规实验方法费用高,并难以插值得到当前激光参数、气流参数条件的实际烧蚀形貌的问题,本发明提供了一种复合材料激光烧蚀形貌快速获取方法。
本发明的具体技术方案是:
一种复合材料激光烧蚀形貌快速获取方法,其实施步骤如下:
步骤1:根据激光辐照复合材料的实际应用场景,分析得到参数空间;其中,参数空间包括激光功率密度I的取值范围、复合材料表面气流速度v的取值范围和复合材料表面气流静压p的取值范围;
步骤2:根据步骤1中获取的参数空间要求,通过常规实验方法获取多组复合材料线烧蚀速率和烧蚀起始时间数据;
步骤3:通过对步骤2中多组复合材料线烧蚀速率和烧蚀起始时间数据进行分析,得到了复合材料线烧蚀速率拟合公式和烧蚀起始时间拟合公式;
步骤4:在参数空间范围内,根据输入的材料线烧蚀速率、烧蚀起始时间通过常用响应面拟合方法求取线烧蚀速率拟合公式、烧蚀起始时间拟合公式中的拟合参数;
步骤5:在计算机内建立复合材料激光烧蚀形貌快速获取程序,在获取程序中构建复合材料三维空间表面网格模型,确定出网格模型的表面激光功率密度以及确定网格模型的表面气流速度和气流静压;
步骤6:启动复合材料激光烧蚀形貌快速获取程序的计算流程,进入材料烧蚀时间步迭代循环计算过程;
根据步骤5中确定的网格模型的表面激光功率密度以及网格模型的表面气流速度和气流静压,计算各网格节点的起始烧蚀时间;
若当前计算时刻未超过网格节点的起始烧蚀时间,则略过该网格节点;
否则根据网格模型的表面激光功率密度以及网格模型的表面气流速度和气流静压,以及步骤3中获得的材料线烧蚀速率拟合公式计算得到该时间步中材料各网格节点处的烧蚀速率,并利用该烧蚀速率和时间步间隔计算出各网格节点的烧蚀深度,之后更新网格节点坐标并重新绘制三维空间网格;
若三维空间网格中前表面网格节点移动到背网格面上时,在对应的节点状态表中标记该节点区域为“穿孔”状态,并将该网格节点绘制为背景色;
步骤7:每次时间步迭代计算结束后,统计处于“穿孔”状态的节点信息,当首次出现该状态节点时,将对应的时刻记为该复合材料的烧蚀穿孔时间;通过空间几何拓扑关系解算,获得烧蚀孔尺寸/面积;其它未处于“穿孔”状态的网格节点空间位置则还原出了材料的烧蚀形貌。
进一步地,上述复合材料线烧蚀速率vablate的拟合公式具体为:
vablate=a1+b1I+c1v+d1p+e1Iv+f1Ip+g1vp
所述复合材料烧蚀起始时间tstart的拟合公式为具体为:
tstart=a2+b2I+c2v+d2p+e2Iv+f2Ip+g2vp
其中,a1、b1、c1、d1、e1、f1、g1、a2、b2、c2、d2、e2、f2、g2为拟合参数。
进一步地,上述网格模型的表面激光功率密度的获取方法为:
根据激光入射角度、功率-时间曲线、光斑分布形态,结合复合材料表面网格节点法向矢量计算,确定出网格模型的表面激光功率密度。
进一步地,上述网格模型的表面激光功率密度的获取方法为:
按照复合材料激光烧蚀形貌快速获取程序接口规范导入实测激光光斑分布文件,结合材料表面网格节点法向矢量计算,确定网格模型的表面激光功率密度,该方式中由于导入了实测光斑,使得激光烧蚀形貌快速计算更具有实际意义。
进一步地,上述网格模型的表面气流速度和气流静压是按照程序接口规范导入材料表面气流状态数据确定的。材料表面气流状态数据可通过数值模拟、实验测量或采信其他可信外来数据等方式获得。
本发明具有以下有益效果:
1、本发明明确了烧蚀形貌计算依据的响应表征量及其主要影响因素,确定了烧蚀响应表征量的拟合公式,设计了快速计算算法,可与确立的烧蚀响应表征量和实验数据紧密结合,快速计算给出材料的烧蚀形貌,给出烧蚀穿孔时间、烧蚀孔尺寸/面积等信息,计算耗时可缩短一至两个数量级。
2、本发明在多种烧蚀响应表征量中,确立了线烧蚀速率、烧蚀起始时间为快速计算所依据的响应表征量,并确定了影响响应表征量的主要因素为激光功率密度、材料表面气流速度和静压,不仅减小了计算量和计算耗时,同时极大的方便了后续实验研究的开展。
3、本发明明确了适用的材料线烧蚀速率、烧蚀起始时间的拟合公式,方便后续实验研究的使用。
4、本发明通过网格节点颜色区分烧蚀状态、仅构建材料表面网格节点、结合烧蚀响应表征量减少不必要的计算、优化网格尺寸参数等多种措施联合使用,在满足计算精度较高的条件时实现了烧蚀形貌快速计算,整体计算过程相对数值模拟计算过程耗时缩短了一至两个数量级,为复合材料制备结构/器件的激光烧蚀后效快速评估提供了输入。
5、本发明的复合材料激光烧蚀形貌快速获取方法,计算过程采用小时间步迭代计算方法,将整个激光辐照时间划分为许多时间间隔小且连续的时间步,在上一个时间步计算结果的基础上计算下一个时间步内的材料烧蚀响应。
6、本发明通过建立复合材料三维空间表面网格模型,并合理的的调试优化三维空间网格尺寸,从而能够匹配计算精度和耗时均在预定范围内。
附图说明
图1为复合材料激光烧蚀形貌快速计算方法流程。其中,虚线框内为时间步迭代循环计算过程。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
在多种烧蚀响应表征量中,通过分析确立了线烧蚀速率和烧蚀起始时间作为复合材料烧蚀形貌的计算依据量;并在影响激光辐照复合材料线烧蚀速率、烧蚀起始时间的诸多因素中明确了激光功率密度、复合材料表面气流速度和气流静压为其主要影响因素,并通过以上的主要影响因素建立了复合材料激光烧蚀形貌快速获取方法。
参见图1,一种复合材料激光烧蚀形貌快速获取方法,主要实施步骤为:
步骤(1):根据激光辐照复合材料的实际应用场景,分析得到参数空间;其中,参数空间包括激光功率密度I的取值范围、复合材料表面气流速度v的取值范围和复合材料表面气流静压p的取值范围;
步骤(2):通过步骤(1)中获取的参数空间,通过常规实验方法获取多组复合材料线烧蚀速率和烧蚀起始时间数据;
步骤(3):通过对步骤(2)中多组复合材料线烧蚀速率和烧蚀起始时间数据进行分析,得到了复合材料线烧蚀速率拟合公式和烧蚀起始时间拟合公式;
复合材料线烧蚀速率vablate的拟合公式具体为:
vablate=a1+b1I+c1v+d1p+e1Iv+f1Ip+g1vp
复合材料烧蚀起始时间tstart的拟合公式为具体为:
tstart=a2+b2I+c2v+d2p+e2Iv+f2Ip+g2vp
上述公式中,a1、b1、c1、d1、e1、f1、g1、a2、b2、c2、d2、e2、f2、g2为拟合参数;
步骤(4):在参数空间范围内通过常用响应面拟合方法求取线烧蚀速率、烧蚀起始时间拟合公式中的拟合参数a1、b1、c1、d1、e1、f1、g1、a2、b2、c2、d2、e2、f2、g2;
步骤(5):在计算机内建立复合材料激光烧蚀形貌快速获取程序,在获取程序中构建复合材料三维空间表面网格模型,确定出网格模型的表面激光时空分布以及确定网格模型的表面气流分布;激光时空分布即为激光功率密度;气流分布即为复合材料表面气流速度和复合材料表面气流静压;
需要说明的一点是:
网格模型的表面激光时空分布的确定方式根据实际需求可以分为以下两种:
A:根据激光入射角度、功率-时间曲线、光斑分布形态等输入参数,结合复合材料表面网格节点法向矢量计算,确定出复合材料表面激光时空分布;
B:按照程序接口规范导入实测激光光斑分布文件,结合材料表面网格节点法向矢量计算,确定材料表面激光时空分布;
网格模型的表面气流分布是按照程序接口规范导入材料表面气流状态数据进行确定的。
步骤(6):启动程序计算流程,进入材料烧蚀时间步迭代循环计算过程;
设定激光辐照总时长为T,将总时长T划分为N个时间步,即t1,t2,…,ti,…,tN,每个时间步的时间间隔为Δt;其中i∈N;
根据步骤(5)中确定的网格模型的表面激光时空分布以及网格模型的表面气流分布,计算各网格节点的起始烧蚀时间;
若第i个时间步对应的计算时刻ti的未超过网格节点的起始烧蚀时间tstart,则略过该网格节点;
若第i个时间步对应的计算时刻ti的超过网格节点的起始烧蚀时间,则利用网格模型的表面激光时空分布以及网格模型的表面气流分布,以及步骤(3)中获得的材料线烧蚀速率拟合公式计算得到该时间步中材料各网格节点处的烧蚀速率,并利用烧蚀速率和时间步间隔Δt计算出各网格节点的烧蚀深度,之后更新网格节点坐标并重新绘制三维空间网格;
若三维空间网格中前表面网格节点移动到背网格面上时,在对应的节点状态表中标记该节点区域为“穿孔”状态,并将该网格节点绘制为背景色;
步骤(7):每次时间步迭代计算结束后,统计处于“穿孔”状态的节点信息,当首次出现该状态节点时,将对应的时刻记为该复合材料的烧蚀穿孔时间;通过空间几何拓扑关系解算,获得烧蚀孔尺寸/面积;其它未处于“穿孔”状态的网格节点空间位置则还原出了材料的烧蚀形貌。
Claims (6)
1.一种复合材料激光烧蚀形貌快速获取方法,其特征在于,包括以下实施步骤:
步骤1:根据激光辐照复合材料的实际应用场景,分析得到参数空间;其中,参数空间包括激光功率密度I的取值范围、复合材料表面气流速度v的取值范围和复合材料表面气流静压p的取值范围;
步骤2:根据步骤1中获取的参数空间要求,通过常规实验方法获取多组复合材料线烧蚀速率和烧蚀起始时间数据;
步骤3:通过对步骤2中多组复合材料线烧蚀速率和烧蚀起始时间数据进行分析,得到了复合材料线烧蚀速率拟合公式和烧蚀起始时间拟合公式;
步骤4:在参数空间范围内,根据输入的材料线烧蚀速率、烧蚀起始时间通过常用响应面拟合方法求取线烧蚀速率拟合公式、烧蚀起始时间拟合公式中的拟合参数;
步骤5:在计算机内建立复合材料激光烧蚀形貌快速获取程序,在获取程序中构建复合材料三维空间表面网格模型,确定出网格模型的表面激光功率密度以及确定网格模型的表面气流速度和气流静压;
步骤6:启动复合材料激光烧蚀形貌快速获取程序的计算流程,进入材料烧蚀时间步迭代循环计算过程;
根据步骤5中确定的网格模型的表面激光功率密度以及网格模型的表面气流速度和气流静压,计算各网格节点的起始烧蚀时间;
若当前计算时刻未超过网格节点的起始烧蚀时间,则略过该网格节点;
否则根据网格模型的表面激光功率密度以及网格模型的表面气流速度和气流静压,以及步骤3中获得的材料线烧蚀速率拟合公式计算得到该时间步中材料各网格节点处的烧蚀速率,并利用该烧蚀速率和时间步间隔计算出各网格节点的烧蚀深度,之后更新网格节点坐标并重新绘制三维空间网格;
若三维空间网格中前表面网格节点移动到背网格面上时,在对应的节点状态表中标记该节点区域为“穿孔”状态,并将该网格节点绘制为背景色;
步骤7:每次时间步迭代计算结束后,统计处于“穿孔”状态的节点信息,当首次出现该状态节点时,将对应的时刻记为该复合材料的烧蚀穿孔时间;通过空间几何拓扑关系解算,获得烧蚀孔尺寸/面积;其它未处于“穿孔”状态的网格节点空间位置则还原出了材料的烧蚀形貌。
2.根据权利要求1所述的复合材料激光烧蚀形貌快速获取方法,其特征在于:
所述复合材料线烧蚀速率vablate的拟合公式具体为:
vablate=a1+b1I+c1v+d1p+e1Iv+f1Ip+g1vp
所述复合材料烧蚀起始时间tstart的拟合公式为具体为:
tstart=a2+b2I+c2v+d2p+e2Iv+f2Ip+g2vp
其中,a1、b1、c1、d1、e1、f1、g1、a2、b2、c2、d2、e2、f2、g2为拟合参数。
3.根据权利要求1所述的复合材料激光烧蚀形貌快速获取方法,其特征在于:网格模型的表面激光功率密度的获取方法为:
根据激光入射角度、功率-时间曲线、光斑分布形态,结合复合材料表面网格节点法向矢量计算,确定出网格模型的表面激光功率密度。
4.根据权利要求1所述的复合材料激光烧蚀形貌快速获取方法,其特征在于:网格模型的表面激光功率密度的获取方法为:
按照复合材料激光烧蚀形貌快速获取程序接口规范导入实测激光光斑分布文件,结合材料表面网格节点法向矢量计算,确定网格模型的表面激光功率密度。
5.根据权利要求1所述的复合材料激光烧蚀形貌快速获取方法,其特征在于:网格模型的表面气流速度和气流静压是按照程序接口规范导入材料表面气流状态数据确定的。
6.根据权利要求1所述的复合材料激光烧蚀形貌快速获取方法,其特征在于:所述复合材料三维空间表面网格模型的网格节点颜色设置为材料材质颜色,且材料材质颜色与模型背景色有明显区分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110955858.3A CN113673127B (zh) | 2021-08-19 | 2021-08-19 | 一种复合材料激光烧蚀形貌快速获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110955858.3A CN113673127B (zh) | 2021-08-19 | 2021-08-19 | 一种复合材料激光烧蚀形貌快速获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113673127A CN113673127A (zh) | 2021-11-19 |
CN113673127B true CN113673127B (zh) | 2023-10-20 |
Family
ID=78544102
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110955858.3A Active CN113673127B (zh) | 2021-08-19 | 2021-08-19 | 一种复合材料激光烧蚀形貌快速获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113673127B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116305765B (zh) * | 2022-12-29 | 2024-06-14 | 中国航天三江集团有限公司 | 高能激光辐照树脂基纤维增强复合材料的仿真方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109623138A (zh) * | 2018-11-28 | 2019-04-16 | 深圳供电局有限公司 | 复合材料表面的处理方法和装置 |
CN112084607A (zh) * | 2020-10-23 | 2020-12-15 | 中国人民解放军国防科技大学 | 激光推进系统中聚合物推进剂烧蚀的二维模拟方法 |
CN112149321A (zh) * | 2020-08-13 | 2020-12-29 | 江苏大学 | 一种脉冲激光除漆的数值模拟方法及其应用 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8987632B2 (en) * | 2009-10-09 | 2015-03-24 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Modification of surface energy via direct laser ablative surface patterning |
-
2021
- 2021-08-19 CN CN202110955858.3A patent/CN113673127B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109623138A (zh) * | 2018-11-28 | 2019-04-16 | 深圳供电局有限公司 | 复合材料表面的处理方法和装置 |
CN112149321A (zh) * | 2020-08-13 | 2020-12-29 | 江苏大学 | 一种脉冲激光除漆的数值模拟方法及其应用 |
CN112084607A (zh) * | 2020-10-23 | 2020-12-15 | 中国人民解放军国防科技大学 | 激光推进系统中聚合物推进剂烧蚀的二维模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113673127A (zh) | 2021-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113673127B (zh) | 一种复合材料激光烧蚀形貌快速获取方法 | |
CN107392875A (zh) | 一种基于k近邻域划分的点云数据去噪方法 | |
CN102452020A (zh) | 一种数控机床刀具温度场和热变形定量分析方法 | |
TW201237647A (en) | Method and system for identifying rare-event failure rates | |
CN104485665A (zh) | 计及风速预测误差时空相关性的动态概率潮流计算方法 | |
CN104102836B (zh) | 一种电力系统快速抗差状态估计方法 | |
CN105554873A (zh) | 一种基于pso-ga-rbf-hop的无线传感器网络定位算法 | |
CN102903145A (zh) | 植物群体形态结构三维重建方法 | |
Zhou et al. | A new acoustic emission source location method using tri-variate kernel density estimator | |
JPH03278183A (ja) | 物理量の表示方法及びその装置 | |
CN117933095A (zh) | 一种基于机器学习的地表发射率实时反演及同化方法 | |
CN105956318A (zh) | 基于改进分裂 h-k 聚类方法的风电场机群划分方法 | |
CN104376026B (zh) | 基于网格和多维树混合结构的表格查找方法 | |
CN107480325B (zh) | 空间变异的非平稳非高斯地震动时程模拟方法 | |
CN109902350A (zh) | 对变截面梁的截面惯性矩进行模型修正中克服模态交换的方法 | |
CN106355641B (zh) | 三维粗糙目标雷达散射截面的仿真方法 | |
CN105335581A (zh) | 一种飞机线束质量分布数据优化方法 | |
CN112512267A (zh) | 基于移动机器人的温度监控方法、系统、设备和存储介质 | |
CN105631066B (zh) | 一种基于模态试验结果的有限元三维模型模态调整方法 | |
Sansregret et al. | Excalibbem: Tool development to calibrate building energy models | |
CN103258144B (zh) | 基于故障录波器数据的在线静态负荷建模方法 | |
CN111008453A (zh) | 一种金属材料二维晶粒组织图的生成方法 | |
CN112217199B (zh) | 一种风电场风速概率分布拟合方法 | |
CN111159857A (zh) | 一种音速喷嘴管壁二维瞬态温度场重构方法 | |
CN109509234A (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 |