CN106327578A - 一种基于不同物理介质的三维温度场插值方法 - Google Patents
一种基于不同物理介质的三维温度场插值方法 Download PDFInfo
- Publication number
- CN106327578A CN106327578A CN201510364238.7A CN201510364238A CN106327578A CN 106327578 A CN106327578 A CN 106327578A CN 201510364238 A CN201510364238 A CN 201510364238A CN 106327578 A CN106327578 A CN 106327578A
- Authority
- CN
- China
- Prior art keywords
- grid node
- model
- sensor points
- dielectric
- areas
- 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 21
- 238000013507 mapping Methods 0.000 claims abstract description 5
- 238000012544 monitoring process Methods 0.000 claims abstract description 4
- 238000012937 correction Methods 0.000 claims description 15
- 229910052757 nitrogen Inorganic materials 0.000 claims description 5
- 229910052698 phosphorus Inorganic materials 0.000 claims description 5
- 238000009877 rendering Methods 0.000 claims description 5
- 238000013480 data collection Methods 0.000 claims description 3
- 238000012800 visualization Methods 0.000 claims description 3
- 229910052739 hydrogen Inorganic materials 0.000 description 8
- 239000001257 hydrogen Substances 0.000 description 8
- 239000011229 interlayer Substances 0.000 description 7
- 239000007788 liquid Substances 0.000 description 5
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 4
- 150000002431 hydrogen Chemical class 0.000 description 4
- -1 interior Tank Substances 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
Landscapes
- Testing Or Calibration Of Command Recording Devices (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于不同物理介质的三维温度场插值方法,属于三维模型映射技术领域,该方法包括:绘制模型网格,获取网格节点坐标,采集模型中的传感器数据,判断传感器与模型节点的位置关系,依据位置关系设置空间插值算法的权重系数,对于每个网格节点,判断传感器点是否与该网格节点处于同一介质模型内,依据传感器点与网格节点所处的介质状态,确定插值权重系数,计算插值结果,得到每个模型的实时三维温度场表征图像,用于对模型进行实时的温度场监控;本发明通过传感器点与模型网格节点的位置关系结合三维模型的物理属性,修正空间插值权重系数,有效改善了模型的三维表征效果和实用性;以上方法不仅局限于三维温度场显示领域,还可以通用于其他物理场三维显示领域。
Description
技术领域
本发明涉及一种基于不同物理介质的三维温度场插值方法,属于三维模型映射技术领域,用于解决不同物理介质模型中的三维显示问题。
背景技术
空间插值方法常常用于GIS空间地理,电力系统,地质勘探等方面,在不同物理介质实体模型中应用尚属首次,由于GIS、空气、地质、电力系统的环境为同一介质,以往的空间插值方法能够适用于这些领域,但是实体模型往往由很多不同材料、不同成分的介质组成,比如用于储存液氢液氧的真空球罐,它由外罐、真空夹层、内罐构成,为了监测整个球罐的温度状态,传感器会分布在球罐的不同位置上,比如外罐内壁面、内罐外壁面、内罐内壁面以及球罐的中间支架上,若球罐中装入液氢,那么球罐的介质就由外罐、真空夹层、内罐、液氢以及氢气组成,若需要计算外罐的温度场分布,用以往的空间插值方法,内罐上的传感器值将映射到外罐上,实际上内罐外罐之间为真空夹层,内罐上的温度不会完全传递到外罐上,而是受到传热系数的影响,因此,本发明基于以往的空间插值方法,引入插值修正系数,通过传感器点与模型网格节点的位置关系结合三维模型的物理属性,修正空间插值权重系数,有效改善了实体模型的三维表征效果和实用性。
发明内容
1、一种基于不同物理介质的三维温度场插值方法,其特征在于,该方法包括以下步骤:
步骤1)绘制模型网格,网格节点的坐标为(xi,yi,zi),i=1,2,...P,P为自然数,将该模型按照物理属性划分为不同介质区域并编号(1,2,...,n),对应的传热系数为(k1,k2,...kn);
步骤2)遍历每个网格节点,判断传感器点与网格节点的位置关系,以第i个网格节点(xi,yi,zi)为研究对象,假设该网格节点所属的介质区域为a,a∈[1,n],遍历每个传感器点(xj′,yj′,zj′),j=1,2,...Q,Q为自然数,其中(xj′,yj′,zj′)是第j个传感器点的位置坐标,依次判断传感器点与网格节点是否重合,若是,则该网格节点的值为传感器点的值,若否,判断传感器点与网格节点的热传导方向上共涉及哪几层介质区域,若传感器点与网格节点处于同一介质区域内,即只涉及该介质区域,则插值修正系数为若传感器点所属的介质区域为a的相邻区域,设为b,b∈[1,n],即涉及a,b介质区域,则插值修正系数为若传感器点所属的介质区域为u,u∈[1,n],而网格节点与传感器点的热传导方向上还包括w个介质区域,即涉及a,b...,u介质区域,则插值修正系数为w为中间介质区域个数,依此类推...
步骤3)在一个数据采集周期开始时,实时采集传感器数据mj′,j=1,2,...N,N为自然数,以第i个网格节点(xi,yi,zi)为研究对象,则第i个网格节点的值mi为
mi=mj′ (dj(x,y,z)=O)
其中dj(x,y,z)是网格节点位置(xi,yi,zi)到传感器点位置(xj′,yj′,zj′)的距离:
[dj(x,y,z)]2=(xi-xj′)2+(yi-yj′)2+(zi-zj′)2
cj为第j个传感器点的插值修正系数;
步骤4)使用三维模型映射技术对网格节点数据mi(xi,yi,zi)进行颜色渲染,得到模型的三维可视化云图,用于对模型进行实时的温度场监控;
步骤5)下一个采集周期一到,返回步骤3)。
附图说明
图1为本发明的一种基于不同物理介质的三维温度场表征方法的实现步骤图。
图2为本发明的以低温真空球罐作为实施例的截面示意图。
图3为本发明的以低温真空球罐作为实施例的三维云图显示效果。
具体实施方式
本发明提出的一种基于不同物理介质的三维温度场插值方法,以低温真空球罐为例,其实施方式如下:
步骤1)绘制低温真空球罐模型网格,网格点的坐标为(xi,yi,zi),i=1,2,...P,P为自然数,将该低温真空球罐按照物理属性划分为外罐、真空夹层、内罐、液氢、氢气介质区域,相应编号为(1,2,3,4,5),对应的传热系数为(k1,k2,k3,k4,k5);
步骤2)遍历每个网格节点,判断传感器点与网格节点的位置关系,以外罐上的第i个网格节点(xi,yi,zi)为研究对象,那么该网格节点所属的介质区域为1,遍历每个传感器点(xj′,yj′,zj′),j=1,2,...N,N为自然数,其中(xj′,yj′,zj′)是第j个传感器点的位置坐标,依次判断传感器点与网格节点是否重合,若是,则该网格节点的值为传感器点的值,若否,判断传感器点与网格节点的热传导方向上共涉及哪几层介质区域,若传感器点位于外罐上,即只涉及1介质区域,则插值修正系数为若传感器点位于内罐外壁上,即涉及1、2介质区域,则插值修正系数为若传感器点位于内罐内壁上,而网格节点与传感器点的热传导方向上还包括真空夹层,即2介质区域,即涉及1、2、3介质区域,则插值修正系数为若传感器点位于球罐的真空支架上且处于液氢中,而网格节点与传感器点的热传导方向上还包括内罐和真空夹层,即2、3介质区域,即涉及1、2、3、4介质区域,则插值修正系数为若传感器点位于球罐的真空支架上且处于氢气中,而网格节点与传感器点的热传导方向上还包括内罐和真空夹层,即2、3介质区域,即涉及1、2、3、5介质区域,则插值修正系数为
步骤3)在一个数据采集周期开始时,实时采集传感器数据mj′,j=1,2,...N,N为自然数,以第i个网格节点(xi,yi,zi)为研究对象,则第i个网格节点的值mi为
mi=mj′ (dj(x,y,z)=0)
其中dj(x,y,z)是网格节点位置(xi,yi,zi)到传感器点位置(xj′,yj′,zj′)的距离:
[dj(x,y,z)]2=(xi-xj′)2+(yi-yj′)2+(zi-zj′)2
cj为第j个传感器点的插值修正系数;
步骤4)使用三维模型映射技术对低温真空球罐网格节点数据mi(xi,yi,zi)进行颜色渲染,得到模型的三维可视化云图,用于对低温真空球罐进行实时的温度场监控;
步骤5)下一个采集周期一到,返回步骤3)。
Claims (1)
1.一种基于不同物理介质的三维温度场插值方法,其特征在于,该方法包括以下步骤:
步骤1)绘制模型网格,网格节点的坐标为(xi,yi,zi),i=1,2,...P,P为自然数,将该模型按照物理属性划分为不同介质区域并编号(1,2,...,n),对应的传热系数为(k1,k2,...kn);
步骤2)遍历每个网格节点,判断传感器点与网格节点的位置关系,以第i个网格节点(xi,yi,zi)为研究对象,假设该网格节点所属的介质区域为a,a∈[1,n],遍历每个传感器点(xj′,yj′,zj′),j=1,2,...Q,Q为自然数,其中(xj′,yj′,zj′)是第j个传感器点的位置坐标,依次判断传感器点与网格节点是否重合,若是,则该网格节点的值为传感器点的值,若否,判断传感器点与网格节点的热传导方向上共涉及哪几层介质区域,若传感器点与网格节点处于同一介质区域内,即只涉及该介质区域,则插值修正系数为若传感器点所属的介质区域为a的相邻区域,设为b,b∈[1,n],即涉及a,b介质区域,则插值修正系数为若传感器点所属的介质区域为u,u∈[1,n],而网格节点与传感器点的热传导方向上还包括w个介质区域,即涉及a,b...,u介质区域,则插值修正系数为w为中间介质区域个数,依此类推...
步骤3)在一个数据采集周期开始时,实时采集传感器数据mj′,j=1,2,...N,N为自然数,以第i个网格节点(xi,yi,zi)为研究对象,则第i个网格节点的值mi为
mi=mj′ (dj(x,y,z)=0)
其中dj(x,y,z)是网格节点位置(xi,yi,zi)到传感器点位置(xj′,yj′,zj′)的距离:
[dj(x,y,z)]2=(xi-xj′)2+(yi-yj′)2+(zi-zj′)2
cj为第j个传感器点的插值修正系数;
步骤4)使用三维模型映射技术对网格节点数据mi(xi,yi,zi)进行颜色渲染,得到模型的三维可视化云图,用于对模型进行实时的温度场监控;
步骤5)下一个采集周期一到,返回步骤3)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510364238.7A CN106327578B (zh) | 2015-06-29 | 2015-06-29 | 一种基于不同物理介质的三维温度场插值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510364238.7A CN106327578B (zh) | 2015-06-29 | 2015-06-29 | 一种基于不同物理介质的三维温度场插值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106327578A true CN106327578A (zh) | 2017-01-11 |
CN106327578B CN106327578B (zh) | 2019-05-17 |
Family
ID=57722448
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510364238.7A Expired - Fee Related CN106327578B (zh) | 2015-06-29 | 2015-06-29 | 一种基于不同物理介质的三维温度场插值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106327578B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109580020A (zh) * | 2018-11-26 | 2019-04-05 | 武汉理工大学 | 船用低速机燃烧室零部件温度场可视化在线测试系统 |
CN109751044A (zh) * | 2019-01-29 | 2019-05-14 | 河南理工大学 | 一种油区热储层顶面温度等值线图编制方法 |
CN112417731A (zh) * | 2020-11-24 | 2021-02-26 | 江苏华路信息科技有限公司 | 基于cae分析结果的温度场插值方法 |
CN112560224A (zh) * | 2020-11-19 | 2021-03-26 | 北京航天试验技术研究所 | 一种火箭发动机试验工艺系统三维表征方法 |
CN113496082A (zh) * | 2020-04-08 | 2021-10-12 | 巴鲁夫公司 | 用于确定至少一个物理值的方法和系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0588759A (ja) * | 1991-09-26 | 1993-04-09 | Nippon Steel Corp | 移動する被加熱材の温度分布解析方法 |
CN103279586A (zh) * | 2013-01-25 | 2013-09-04 | 北京空间飞行器总体设计部 | 一种有限元温度场的二次映射方法 |
CN103614498A (zh) * | 2013-12-18 | 2014-03-05 | 中南大学 | 一种高炉炉壁三维温度场重建方法及计算机监控系统 |
CN103617367A (zh) * | 2013-12-06 | 2014-03-05 | 三峡大学 | 电磁场-流场-温度场耦合计算中的异型网格映射方法 |
-
2015
- 2015-06-29 CN CN201510364238.7A patent/CN106327578B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0588759A (ja) * | 1991-09-26 | 1993-04-09 | Nippon Steel Corp | 移動する被加熱材の温度分布解析方法 |
CN103279586A (zh) * | 2013-01-25 | 2013-09-04 | 北京空间飞行器总体设计部 | 一种有限元温度场的二次映射方法 |
CN103617367A (zh) * | 2013-12-06 | 2014-03-05 | 三峡大学 | 电磁场-流场-温度场耦合计算中的异型网格映射方法 |
CN103614498A (zh) * | 2013-12-18 | 2014-03-05 | 中南大学 | 一种高炉炉壁三维温度场重建方法及计算机监控系统 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109580020A (zh) * | 2018-11-26 | 2019-04-05 | 武汉理工大学 | 船用低速机燃烧室零部件温度场可视化在线测试系统 |
CN109751044A (zh) * | 2019-01-29 | 2019-05-14 | 河南理工大学 | 一种油区热储层顶面温度等值线图编制方法 |
CN113496082A (zh) * | 2020-04-08 | 2021-10-12 | 巴鲁夫公司 | 用于确定至少一个物理值的方法和系统 |
CN112560224A (zh) * | 2020-11-19 | 2021-03-26 | 北京航天试验技术研究所 | 一种火箭发动机试验工艺系统三维表征方法 |
CN112417731A (zh) * | 2020-11-24 | 2021-02-26 | 江苏华路信息科技有限公司 | 基于cae分析结果的温度场插值方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106327578B (zh) | 2019-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106327578A (zh) | 一种基于不同物理介质的三维温度场插值方法 | |
CN105513131B (zh) | 一种顾及褶皱的三维地质结构模型角点网格剖分方法 | |
CN104077134B (zh) | 获取数据中心机房平面温度云图的方法和系统 | |
CN103544728B (zh) | 一种在实景地图中基于法向量的被探面过滤方法 | |
Vergés et al. | 3D pore analysis of sedimentary rocks | |
CN111612903B (zh) | 一种基于混合数据模型的地质数据可视化方法 | |
CN108241777B (zh) | 基于非结构网格有限元法计算水合物沉积物中渗流速度场的方法 | |
CN105869202A (zh) | 电子地图的数字地面模型数据生成方法和装置 | |
CN104007479B (zh) | 一种基于多尺度剖分的电离层层析和电离层延迟改正方法 | |
CN107833282A (zh) | 一种地形建模和网格生成方法及装置 | |
CN109376433B (zh) | 基于土壤非饱和水和地下水耦合的区域水流运动模拟方法 | |
CN101051395A (zh) | 一种基于地球物理场数据的地质体三维可视化方法 | |
CN106682262A (zh) | 一种获取飞行器流场的数值模拟方法 | |
CN104966317A (zh) | 一种基于矿体轮廓线的三维自动建模方法 | |
CN104134203A (zh) | 一种近景摄影测量的快速密集匹配法 | |
CN107766516A (zh) | 一种基于边界保留的地理地图变形可视化方法 | |
CN109636912A (zh) | 应用于三维声呐图像重构的四面体剖分有限元插值方法 | |
Browne et al. | Fast three dimensional r-adaptive mesh redistribution | |
CN102117494A (zh) | 一种三维建模中的数据分层分块方法及装置 | |
CN109544687B (zh) | 一种定向分层分户图扩展方法 | |
CN104574513B (zh) | 基于几何形态的准确刻画夹层三维分布的表征方法 | |
CN104199949B (zh) | 一种实时环境下利用网格实现大规模动态目标可视化方法 | |
CN102520406A (zh) | 基于影像坐标映射的sar立体模型构建与量测方法 | |
CN105653881A (zh) | 基于多密度层次的流场可视化方法 | |
CN101533102A (zh) | 二维复杂结构三角网射线追踪全局方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190517 |