CN106528935A - 一种用于捷联式海洋重力仪的热设计方法 - Google Patents
一种用于捷联式海洋重力仪的热设计方法 Download PDFInfo
- Publication number
- CN106528935A CN106528935A CN201610884717.6A CN201610884717A CN106528935A CN 106528935 A CN106528935 A CN 106528935A CN 201610884717 A CN201610884717 A CN 201610884717A CN 106528935 A CN106528935 A CN 106528935A
- Authority
- CN
- China
- Prior art keywords
- strapdown
- temperature
- fan
- gravimeter
- thermal design
- 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.)
- Pending
Links
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
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
一种用于捷联式海洋重力仪的热设计方法,为重力仪制定热设计方案,采用保温材料隔热降低自动传热、加温片加热和风道风扇***散热增大主动可控传热,为捷联式海洋重力仪提供一级温控。在有限元分析软件中利用参数目标驱动优化方法对捷联式海洋重力仪三维模型进行仿真分析,得到温度场响应曲面。在温箱中进行实物温度试验,根据规律趋势调整参数,得到热设计参数,实现对重力仪的热设计。本发明给出了热设计方案,并能够快速精确确定热设计参数,利用该方法对重力仪进行热设计,能够有效保证重力仪在海洋重力测量时温度恒定,提高测量精度。
Description
技术领域
本发明涉及一种用于捷联式海洋重力仪的热设计方法,属于重力测量产品设计领域。
背景技术
目前,重力测量在资源勘探、固体潮监测等民用和科学研究方面应用广泛。航空重力测量***的出现使得高效、大范围的重力测量成为可能。在资源勘探方面要求重力测量***的精度优于1mGal(1mGal=10-5m/s2=1μg)。
重力测量在作业实施时,无论载体是飞机、车辆还是船舶,基本上都没有温度调节装置,载体内部温度受海拔高度、太阳辐射的影响很大,特别是无人机环境中,尤为恶劣,而重力仪***所用的核心测试器件都对温度极其敏感,在恶劣的环境下重力测量精度降低,甚至测试数据无效。由于海洋重力测量的特殊性,每次测量作业成本很高,因此对重力仪进行热设计,使重力仪本体温度恒定,迫在眉睫。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供了一种用于捷联式海洋重力仪的热设计方法,给出了热设计方案,并能够快速精确确定热设计参数,利用该方法对重力仪进行热设计,能够有效保证重力仪在海洋重力测量时温度恒定,提高测量精度。
本发明的技术解决方案是:一种用于捷联式海洋重力仪的热设计方法,步骤如下:
(1)在捷联式海洋重力仪本体以及箱体内壁上沿周向均匀分布加温片,在箱体外表面敷设一层隔热材料,隔热材料与箱体侧壁之间设计有多条风道,在箱体上盖和下盖上分别布置多个风扇;
(2)对捷联式海洋重力仪进行三维建模;
(3)根据捷联式海洋重力仪中发热元器件的功率确定热源参数,设定外界环境温度、模型中不同材料的热参数和辐射系数,设计各个加温片的功率和风扇的风量参数,在有限元分析软件中利用目标驱动优化方法根据上述参数对捷联式海洋重力仪三维模型进行温度仿真分析,得到本体温度场响应曲面;
(4)根据本体温度场响应曲面确定加温片功率和风扇风量参数对捷联式海洋重力仪本体温度的影响规律趋势;根据捷联式海洋重力仪传感器的最佳工作温度范围和步骤(3)得到的本体温度场响应曲面以及海洋重力测量的外界环境温度范围,确定本体稳定工作温度和不同外界环境温度下每个加温片的功率和每个风扇的风量;
(5)构建捷联式海洋重力仪实物模型,其材料的热参数和辐射系数与步骤(3)相同,在高低温箱中实现所要求的外界环境温度,根据仿真分析结果对加温片功率和风扇风量进行设置,进行实物温度试验,得到本体温度;
(6)如果本体温度为步骤(4)确定的稳定工作温度,则根据外界环境温度按照步骤(4)确定的加温片功率和风扇风量对重力仪进行热设计;否则,根据加温片功率和风扇风量参数对重力仪本体温度的影响规律趋势对每个加温片的功率和每个风扇的风量进行微调,使本体温度达到步骤(4)确定的稳定工作温度,利用此时每个加温片的功率和每个风扇的风量对重力仪进行热设计。
所述步骤(3)中在有限元分析软件中对捷联式海洋重力仪三维模型进行仿真分析时,网格模型采用六面体和四面体的实体混合模型,网格大小控制在5mm-7mm。
所述步骤(4)中确定不同外界环境温度下每个加温片功率和每个风扇风量的方法为:
(3.1)根据本体稳定工作温度和本体温度场响应曲面确定不同环境温度下每个加温片功率范围和每个风扇风量范围;
(3.2)每个加温片的功率p按照下式计算:
其中p1、p2分别为该加温片功率范围的中间值;
每个风扇的风量f按照下式计算:
f1、f2为该风扇风量范围的中间值。
所述保温材料采用聚氨酯泡沫,导热系数0.02,厚度为8mm-15mm。
所述箱体上盖上的风扇沿上盖中心对称,且靠近上盖的边缘;箱体下盖上的风扇沿下盖中心对称,且靠近下盖的边缘,所述风扇放置方式满足风道下进风、上吸风。。
所述加温片均匀布置在本体和箱体上,保证本体和箱体温度快速均匀扩散。
本发明与现有技术相比的优点在于:
(1)本发明热设计方案采用加温片加热,保温材料隔热和设计风道利用风扇进行散热,保温材料采用导热系数0.02、厚度为8mm-15mm的聚氨酯泡沫,能够将重力仪与外界环境隔离成分开的两个部分,减小了重力仪与外界的自动传热,保证隔热效果;通过加温片加热以及设计风道利用风扇进行散热方式增大主动可控传热,在特定的环境下加快了与外界的主动可控传热。风道采用下进风,上吸风的设置,增大了空气对流散热。因此本发明的热设计方案既能够适应低温的外界环境,又能在高温的环境下保持相对稳定的散热。
(2)本发明采用目标驱动优化方法通过一次仿真分析即可得到响应曲面,能够快速确定热设计参数,与选择不同参数进行多次实物试验确定热设计参数相比,本发明降低了成本,且参数设计更准确,同时能够对各种设计参数对重力仪本体温度的影响趋势做出全面准确的判断。
(3)本发明在确定热设计参数后,还要进行实物温度试验,充分保证热设计参数可靠,确保重力仪在海洋重力测量时温度恒定,提高测量精度,降低实施风险。
附图说明
图1为本发明工作流程图;
图2为重力仪的热设计模型示意图,其中(a)为本体、隔热材料和风扇示意图,(b)为风道示意图;
图3为某组参数下重力仪热设计模型的温度场图;
图4为利用本发明方法设计的某重力仪在海洋重力测量时的温度变化曲线。
具体实施方式
本发明提出了一种用于捷联式海洋重力仪的热设计方法,采用加温片加热,保温材料隔热和设计风道利用风扇进行散热方式,既能够适应低温的外界环境,又能在高温的环境下保持相对稳定的散热,利用仿真分析方法采用参数目标驱动优化方法得到响应曲面,进而指导实物温度试验来确定热设计参数,如图1所示,其具体步骤如下所述:
(1)在捷联式海洋重力仪本体以及箱体内壁上沿周向均匀布置加温片,在箱体外表面敷设一层保温材料,保温材料与箱体侧壁之间设计有多条风道,在箱体上盖和下盖上分别布置多个风扇。如图2所示为重力仪的热设计模型图,其中(a)为本体、隔热材料和风扇示意图,(b)为风道示意图。其中加温片均匀布置在本体和箱体上,能够保证本体和箱体温度快速均匀扩散。
箱体上盖上的风扇沿上盖中心对称,且靠近上盖的边缘;箱体下盖上的风扇沿下盖中心对称,且靠近下盖的边缘。下盖的风扇和上盖的风扇均向上吹风,以使风道下进风、上吸风,从而增大空气对流,提高散热效果。
为了降低自动传热、实现有效隔热,本发明隔热材料采用聚氨酯泡沫,导热系数0.02,厚度为8mm-15mm,优选为10mm。本发明设计风道沿箱体外表面均匀分布,风道表面积与箱体体积比例为27%,高温环境下对流散热效果明显。
(2)对捷联式海洋重力仪进行三维建模;
(3)根据捷联式海洋重力仪中发热元器件的功率确定热源参数,设定模型中不同材料的热参数(比热容、传热系数等)和辐射系数,设计各个加温片的功率和风扇的风量参数以及不同的外界环境温度,在有限元分析软件中利用目标驱动优化方法对外界环境温度、各加温片功率、各风扇风量进行温度仿真分析,得到本体温度场响应曲面。图3为某组参数下重力仪模型的温度场图。
进行仿真分析时,网格模型采用六面体和四面体的实体混合模型,网格大小控制在5mm左右,网格划分完成后节点数48811,单元数27859。
(4)根据本体温度场响应曲面确定加温片功率和风扇风量参数对捷联式海洋重力仪本体温度的影响规律趋势;根据捷联式海洋重力仪传感器的最佳工作温度范围和步骤(3)得到的本体温度场响应曲面以及海洋重力测量的外界环境温度范围,确定本体稳定工作温度和不同环境温度下每个加温片的功率范围和每个风扇的风量范围。每个加温片的功率为其功率范围的中间值,每个风扇的风量为其风量范围的中间值。
(5)构建捷联式海洋重力仪实物模型,其材料的热参数和辐射系数与步骤(3)相同,在高低温箱中实现所要求的外界环境温度,根据仿真分析结果对加温片功率和风扇风量进行设置,进行实物温度试验,得到本体温度;
(6)如果本体温度为本体稳定工作温度,则按照步骤(4)确定的参数实现对重力仪的热设计;否则,根据加温片功率和风扇风量参数对重力仪本体温度的影响规律趋势对每个加温片的功率和每个风扇的风量进行微调,使本体温度达到本体稳定工作温度,利用此时每个加温片的功率和每个风扇的风量实现对重力仪的热设计。
在实际对重力仪进行热设计时,为了保证重力仪本体的温度稳定,可以进一步设计有控制器,事先将不同外界环境温度对应的每个加温片的功率和每个风扇的风量写入控制器,然后根据当前外界环境温度对重力仪进行热设计,当外界环境温度变化时,控制器可以根据事先写入的内容对加温片的功率和风扇的风量进行调节。
图4为利用本发明方法对某重力仪进行热设计后,进行海洋重力测量时的温度变化曲线。由图可知,30天的测量过程中,重力仪本体的温度基本稳定在58℃。从而有效保证了安装在本体上的核心传感器能够高精度的测量,提高了重力仪整体测量精度。
本发明未详细描述内容为本领域技术人员公知技术。
Claims (6)
1.一种用于捷联式海洋重力仪的热设计方法,其特征在于步骤如下:
(1)在捷联式海洋重力仪本体以及箱体内壁上沿周向均匀分布加温片,在箱体外表面敷设一层隔热材料,隔热材料与箱体侧壁之间设计有多条风道,在箱体上盖和下盖上分别布置多个风扇;
(2)对捷联式海洋重力仪进行三维建模;
(3)根据捷联式海洋重力仪中发热元器件的功率确定热源参数,设定外界环境温度、模型中不同材料的热参数和辐射系数,设计各个加温片的功率和风扇的风量参数,在有限元分析软件中利用目标驱动优化方法根据上述参数对捷联式海洋重力仪三维模型进行温度仿真分析,得到本体温度场响应曲面;
(4)根据本体温度场响应曲面确定加温片功率和风扇风量参数对捷联式海洋重力仪本体温度的影响规律趋势;根据捷联式海洋重力仪传感器的最佳工作温度范围和步骤(3)得到的本体温度场响应曲面以及海洋重力测量的外界环境温度范围,确定本体稳定工作温度和不同外界环境温度下每个加温片的功率和每个风扇的风量;
(5)构建捷联式海洋重力仪实物模型,其材料的热参数和辐射系数与步骤(3)相同,在高低温箱中实现所要求的外界环境温度,根据仿真分析结果对加温片功率和风扇风量进行设置,进行实物温度试验,得到本体温度;
(6)如果本体温度为步骤(4)确定的稳定工作温度,则根据外界环境温度按照步骤(4)确定的加温片功率和风扇风量对重力仪进行热设计;否则,根据加温片功率和风扇风量参数对重力仪本体温度的影响规律趋势对每个加温片的功率和每个风扇的风量进行微调,使本体温度达到步骤(4)确定的稳定工作温度,利用此时每个加温片的功率和每个风扇的风量对重力仪进行热设计。
2.根据权利要求1所述的一种用于捷联式海洋重力仪的热设计方法,其特征在于:所述步骤(3)中在有限元分析软件中对捷联式海洋重力仪三维模型进行仿真分析时,网格模型采用六面体和四面体的实体混合模型,网格大小控制在5mm-7mm。
3.根据权利要求1所述的一种用于捷联式海洋重力仪的热设计方法,其特征在于:所述步骤(4)中确定不同外界环境温度下每个加温片功率和每个风扇风量的方法为:
(3.1)根据本体稳定工作温度和本体温度场响应曲面确定不同环境温度下每个加温片功率范围和每个风扇风量范围;
(3.2)每个加温片的功率p按照下式计算:
其中p1、p2分别为该加温片功率范围的中间值;
每个风扇的风量f按照下式计算:
f1、f2为该风扇风量范围的中间值。
4.根据权利要求1所述的一种用于捷联式海洋重力仪的热设计方法,其特征在于:所述保温材料采用聚氨酯泡沫,导热系数0.02,厚度为8mm-15mm。
5.根据权利要求1所述的一种用于捷联式海洋重力仪的热设计方法,其特征在于:所述箱体上盖上的风扇沿上盖中心对称,且靠近上盖的边缘;箱体下盖上的风扇沿下盖中心对称,且靠近下盖的边缘,所述风扇放置方式满足风道下进风、上吸风。
6.根据权利要求5所述的一种用于捷联式海洋重力仪的热设计方法,其特征在于:所述加温片均匀布置在本体和箱体上,保证本体和箱体温度快速均匀扩散。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610884717.6A CN106528935A (zh) | 2016-10-10 | 2016-10-10 | 一种用于捷联式海洋重力仪的热设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610884717.6A CN106528935A (zh) | 2016-10-10 | 2016-10-10 | 一种用于捷联式海洋重力仪的热设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106528935A true CN106528935A (zh) | 2017-03-22 |
Family
ID=58333092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610884717.6A Pending CN106528935A (zh) | 2016-10-10 | 2016-10-10 | 一种用于捷联式海洋重力仪的热设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106528935A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526868A (zh) * | 2017-07-17 | 2017-12-29 | 安徽四创电子股份有限公司 | 一种用于雷达电子机柜***的热设计方法 |
CN108594318A (zh) * | 2018-07-12 | 2018-09-28 | 广州海洋地质调查局 | 一种用于水下的重力仪 |
CN108897058A (zh) * | 2018-07-12 | 2018-11-27 | 中国人民解放军国防科技大学 | 一种捷联式水下动态重力测量仪承压温控舱 |
CN113536515A (zh) * | 2020-04-17 | 2021-10-22 | 株洲中车时代电气股份有限公司 | 一种针对高海拔的热设计参数等效换算方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101424534A (zh) * | 2008-12-09 | 2009-05-06 | 东南大学 | 惯性/重力组合导航半实物模拟装置 |
CN104914716A (zh) * | 2015-04-08 | 2015-09-16 | 中国人民解放军海军工程大学 | 海洋航空重力测量平台自适应控制与容错保护***及方法 |
CN105204083A (zh) * | 2015-09-14 | 2015-12-30 | 武汉大学 | 基于绝对重力仪的垂直重力梯度提取方法及*** |
-
2016
- 2016-10-10 CN CN201610884717.6A patent/CN106528935A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101424534A (zh) * | 2008-12-09 | 2009-05-06 | 东南大学 | 惯性/重力组合导航半实物模拟装置 |
CN104914716A (zh) * | 2015-04-08 | 2015-09-16 | 中国人民解放军海军工程大学 | 海洋航空重力测量平台自适应控制与容错保护***及方法 |
CN105204083A (zh) * | 2015-09-14 | 2015-12-30 | 武汉大学 | 基于绝对重力仪的垂直重力梯度提取方法及*** |
Non-Patent Citations (1)
Title |
---|
郭纬川马杰楚贤李东明: "重力仪散热装置的热仿真分析", 《导航定位与授时》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526868A (zh) * | 2017-07-17 | 2017-12-29 | 安徽四创电子股份有限公司 | 一种用于雷达电子机柜***的热设计方法 |
CN107526868B (zh) * | 2017-07-17 | 2020-09-25 | 安徽四创电子股份有限公司 | 一种用于雷达电子机柜***的热设计方法 |
CN108594318A (zh) * | 2018-07-12 | 2018-09-28 | 广州海洋地质调查局 | 一种用于水下的重力仪 |
CN108897058A (zh) * | 2018-07-12 | 2018-11-27 | 中国人民解放军国防科技大学 | 一种捷联式水下动态重力测量仪承压温控舱 |
CN108897058B (zh) * | 2018-07-12 | 2019-11-29 | 中国人民解放军国防科技大学 | 一种捷联式水下动态重力测量仪承压温控舱 |
CN113536515A (zh) * | 2020-04-17 | 2021-10-22 | 株洲中车时代电气股份有限公司 | 一种针对高海拔的热设计参数等效换算方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106528935A (zh) | 一种用于捷联式海洋重力仪的热设计方法 | |
CN109883660B (zh) | 一种热模拟试验控制方法 | |
EP3845465B1 (en) | Aerostat icing characteristic numerical simulation and experimental verification system | |
Kuang | Weakly forced mock Walker cells | |
CN113549997B (zh) | 单晶生长的方法、装置及单晶体 | |
CN104730107B (zh) | 滨海电厂温排水数学模型表面散热系数实验室观测方法 | |
Zhang et al. | 3D Model‐based simulation analysis of energy consumption in hot air drying of corn kernels | |
Tekasakul et al. | Three-dimensional numerical modeling of heat and moisture transfer in natural rubber sheet drying process | |
CN103994652A (zh) | 薄层干燥实验设备及其干燥工艺 | |
CN104749347A (zh) | 研究地铁区间隧道土壤温湿度场演化规律的实验装置 | |
CN109684717A (zh) | 基于量纲分析的油炒烹饪中表面换热系数的预测方法 | |
Leng et al. | A parallel computational model for three-dimensional, thermo-mechanical Stokes flow simulations of glaciers and ice sheets | |
Madhiyanon et al. | High temperature spouted bed paddy drying with varied downcomer air flows and moisture contents: Effects on drying kinetics, critical moisture content, and milling quality | |
CN107991116A (zh) | 一种基于点阵热源的散热器性能评估与面积计算平台 | |
CN204374028U (zh) | 一种带加热功能的高温流体密度测试装置 | |
Mao et al. | A convective heat transfer coefficient algorithm for thermal analysis of machine tools considering a temperature change | |
CN112560251B (zh) | 一种伪装网表面温度数值计算的简化方法 | |
CN101968521A (zh) | 电器产品热稳定试验加热方法及其专用烘箱 | |
CN115329631A (zh) | 一种基于数字孪生的激光切割热误差补偿方法 | |
Ukaji et al. | A comparison of laboratory experiments and numerical simulations of steady baroclinic waves produced in a differentially heated rotating fluid annulus | |
CN105760574A (zh) | 一种模拟固态发酵过程中温度场分布的方法 | |
Fang et al. | Characteristics of thermal stress evolution during the cooling stage of multicrystalline silicon | |
Lutsko et al. | The influence of meridional gradients in insolation and longwave optical depth on the climate of a gray radiation GCM | |
CN204065088U (zh) | 纺织品干燥速率测试机 | |
CN117707257B (zh) | 一种绿豆芽培育温度控制方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170322 |