CN114067079A - 复杂曲面电磁波矢量动态可视化方法 - Google Patents
复杂曲面电磁波矢量动态可视化方法 Download PDFInfo
- Publication number
- CN114067079A CN114067079A CN202111401043.7A CN202111401043A CN114067079A CN 114067079 A CN114067079 A CN 114067079A CN 202111401043 A CN202111401043 A CN 202111401043A CN 114067079 A CN114067079 A CN 114067079A
- Authority
- CN
- China
- Prior art keywords
- electromagnetic wave
- vector
- tbl
- data table
- data
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
Abstract
本发明公开了一种复杂曲面电磁波矢量动态可视化方法,包括以下步骤:S1.从nastran网格文件中读取复杂曲面模型的点列表[P]N和点连接列表形成三角面元列表[Tri]M;初始化一张维度为N×18的数据表格[Tbl]N,18;S2.对三角面元列表中的每一个面元计算单位法向矢量ni,并计算与ni垂直的两正交单位向量bi与ti,并根据计算结果对数据表格[Tbl]N,18的部分数据进行填充;S3.对步骤S2中得到的数据表格进行处理,计算数据表格[Tbl]N,18中的另一部分数据,对数据表格进行填充,得到完整的数据表格;S4.根据完整的数据表格,对复杂曲面电磁波矢量进行动态可视化处理。本发明能够将将电磁波的波动耦合动态地进行显示,有效实现了复杂曲面上的电磁波的矢量波动情况可视化。
Description
技术领域
本发明涉及电磁波领域,特别是涉及一种复杂曲面电磁波矢量动态可视化方法。
背景技术
在电磁波领域,为了对电磁波进行有效的观测,常常需要对曲面电磁波矢量进行可视化处理,但是,传统的场矢量可视化只考虑了单一方向矢量场的情形,并没有考虑正交场之间的耦合情况;无法动态的将电磁波的波动耦合动态地进行显示,电磁波波动显示速度无法调节,这就为电磁波的矢量动态可视化带来了诸多不便,特别是在涉及到复杂曲面的电磁波矢量时,传统显示方法无法动态显示复杂曲面上的电磁波的矢量波动情况。
发明内容
本发明的目的在于克服现有技术的不足,提供一种复杂曲面电磁波矢量动态可视化方法,能够将将电磁波的波动耦合动态地进行显示,有效实现了复杂曲面上的电磁波的矢量波动情况可视化。
本发明的目的是通过以下技术方案来实现的:复杂曲面电磁波矢量动态可视化方法,包括以下步骤:
S1.从nastran网格文件中读取复杂曲面模型的点列表[P]N和点连接列表形成三角面元列表[Tri]M;其中,N为点数,M为三角面元数;初始化一张维度为N×18的数据表格[Tbl]N,18;
S2.对三角面元列表中的每一个面元计算单位法向矢量ni,并计算与ni垂直的两正交单位向量bi与ti,并根据计算结果对数据表格[Tbl]N,18的部分数据进行填充;
S3.对步骤S2中得到的数据表格进行处理,计算数据表格[Tbl]N,18中的另一部分数据,对数据表格进行填充,得到完整的数据表格;
S4.根据完整的数据表格,对复杂曲面电磁波矢量进行动态可视化处理。
进一步地,所述步骤S2包括以下子步骤:
S201.遍历面元列表[Tri]M,得到第i个面元Tri[i]的三个点序号分别为i1,i2,i3,三个点坐标分别为Pi1(x1,y1,z1),Pi2(x2,y2,z2)和Pi3(x3,y3,z3);
其中,×表示矢量叉乘,|·|表示矢量的长度;
S203.计算与ni互为垂直的两正交单位向量:
S204.对[Tbl]N,18的部分数据进行填充:
S205.在i=1,2…,M时,循环执行步骤S201~步骤S204,直到遍历完面元列表中所有面元为止。
进一步地,所述步骤S3包括以下子步骤:
式中,|·|是复数幅度函数,而Angle是复数的辐角函数,单位为度;
S304.在iRow=1,2…,N时,重复执行S301~S303,直至完成数据表格[Tbl]N,18所有行的数据填充。
进一步地,所述步骤S4包括以下子步骤:
S401.设定电磁波初始显示相位Phs0=0,由用户根据动态刷新速度设定电磁波相位步进Pstep,由慢到快到非常快依次设定为k1、k2、k3、k4、k5、k6六档;初始化maxAc=0;
其中,α为箭头长度显示调整参数,设置为D1、D2、D3、D4、D5、D6六档,分别表示用户设置箭头相对模型的不同的显示比例;Δs为面元网格的平均边长;为由步骤S404得到;max(Ac)由步骤S402得到;
S407.在i=1,2…,N时,循环步骤S403~步骤S406,绘制Phs0相位下的所有箭头;
S408.等待100ms,将Phs0+Pstep赋值给Phs0;
S409.重复执行步骤S403~S408,直到用户停止电磁波矢量动态绘制为止,绘制过程结束。
式中,RAD=π/180,j=sqrt(-1)
A4、在i=1,2,…N时,循环执行步骤A1~A3,循环结束后输出maxAc。
本发明的有益效果是:1、电磁波在复杂曲面表面上的任意点的电磁场都是两正交极化场的合成,本申请相比传统的标量电磁场的可视化更能揭示电磁场在曲面表面上的场的分解合成现象,更接近电磁波耦合的本质;为复杂曲面上耦合电磁场的可视化提供了技术手段和方法。
2、本申请的可视化方法能够显示不同时刻下的复杂曲面表面上的电磁波动态分布情况,可以方便天线布局工程师从可视化中得到复杂曲面表面上电流的流动方向,方便天线孔径的位置和方向设计,减少了优化设计前的不断试错的机会和浪费的时间,提高了设计效率。
附图说明
图1为本发明的方法流程图;
图2为ni与垂直的两正交单位向量bi与ti的示意图;
图3为实施例中Phs0=30时的效果示意图;
图4为实施例中Phs0=31时的效果示意图;
图5为实施例中Phs0=32时的效果示意图;
图6为实施例中Phs0=33时的效果示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图1所示,复杂曲面电磁波矢量动态可视化方法,包括以下步骤:
S1.从nastran网格文件中读取复杂曲面模型的点列表[P]N和点连接列表形成三角面元列表[Tri]M;其中,N为点数,M为三角面元数;初始化一张维度为N×18的数据表格[Tbl]N,18;
S2.对三角面元列表中的每一个面元计算单位法向矢量ni,并计算与ni垂直的两正交单位向量bi与ti,并根据计算结果对数据表格[Tbl]N,18的部分数据进行填充;
S3.对步骤S2中得到的数据表格进行处理,计算数据表格[Tbl]N,18中的另一部分数据,对数据表格进行填充,得到完整的数据表格;
S4.根据完整的数据表格,对复杂曲面电磁波矢量进行动态可视化处理。
进一步地,所述步骤S2包括以下子步骤:
S201.遍历面元列表[Tri]M,得到第i个面元Tri[i]的三个点序号分别为i1,i2,i3,三个点坐标分别为Pi1(x1,y1,z1),Pi2(x2,y2,z2)和Pi3(x3,y3,z3);
其中,×表示矢量叉乘,|·|表示矢量的长度;
S203.如图2所示,计算与ni互为垂直的两正交单位向量:
S204.对[Tbl]N,18的部分数据进行填充:
S205.在i=1,2…,M时,循环执行步骤S201~步骤S204,直到遍历完面元列表中所有面元为止。
进一步地,所述步骤S3包括以下子步骤:
式中,|·|是复数幅度函数,而Angle是复数的辐角函数,单位为度;
S304.在iRow=1,2…,N时,重复执行S301~S303,直至完成数据表格[Tbl]N,18所有行的数据填充。
进一步地,所述步骤S4包括以下子步骤:
S401.设定电磁波初始显示相位Phs0=0,由用户根据动态刷新速度设定电磁波相位步进Pstep,由慢到快到非常快依次设定为k1、k2、k3、k4、k5、k6(例如0.01,0.1,1,3,5,10)六档;初始化maxAc=0;
其中,α为箭头长度显示调整参数,设置为D1、D2、D3、D4、D5、D6(例如0.1,0.2,0.5,1,1.3,2)六档,分别表示用户设置箭头相对模型的不同的显示比例;Δs为面元网格的平均边长;为由步骤S404得到;max(Ac)由步骤S402得到;
S407.在i=1,2…,N时,循环步骤S403~步骤S406,绘制Phs0相位下的所有箭头;
S408.等待100ms,将Phs0+Pstep赋值给Phs0;
S409.重复执行步骤S403~S408,直到用户停止电磁波矢量动态绘制为止,绘制过程结束。
式中,RAD=π/180,j=sqrt(-1)
A4、在i=1,2,…N时,循环执行步骤A1~A3,循环结束后输出maxAc。
在本申请的实施例中,利用均匀平面电磁波在复杂曲面上感应表面电磁波在各个显示相位上的箭头状态实例来阐述该方法的实际效果图3~6所示,其中图3为Phs0=30时的效果示意图,图4为Phs0=31时的效果示意图,图5为Phs0=32时的效果示意图,图6为Phs0=33时的效果示意图。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.复杂曲面电磁波矢量动态可视化方法,其特征在于:包括以下步骤:
S1.从nastran网格文件中读取复杂曲面模型的点列表[P]N和点连接列表形成三角面元列表[Tri]M;其中,N为点数,M为三角面元数;初始化一张维度为N×18的数据表格[Tbl]N,18;
S2.对三角面元列表中的每一个面元计算单位法向矢量ni,并计算与ni垂直的两正交单位向量bi与ti,并根据计算结果对数据表格[Tbl]N,18的部分数据进行填充;
S3.对步骤S2中得到的数据表格进行处理,计算数据表格[Tbl]N,18中的另一部分数据,对数据表格进行填充,得到完整的数据表格;
S4.根据完整的数据表格,对复杂曲面电磁波矢量进行动态可视化处理。
2.根据权利要求1所述的复杂曲面电磁波矢量动态可视化方法,其特征在于:所述步骤S2包括以下子步骤:
S201.遍历面元列表[Tri]M,得到第i个面元Tri[i]的三个点序号分别为i1,i2,i3,三个点坐标分别为Pi1(x1,y1,z1),Pi2(x2,y2,z2)和Pi3(x3,y3,z3);
其中,×表示矢量叉乘,|·|表示矢量的长度;
S203.计算与ni互为垂直的两正交单位向量:
S204.对[Tbl]N,18的部分数据进行填充:
S205.在i=1,2…,M时,循环执行步骤S201~步骤S204,直到遍历完面元列表中所有面元为止。
3.根据权利要求1所述的复杂曲面电磁波矢量动态可视化方法,其特征在于:所述步骤S3包括以下子步骤:
式中,|·|是复数幅度函数,而Angle是复数的辐角函数,单位为度;
S304.在iRow=1,2…,N时,重复执行S301~S303,直至完成数据表格[Tbl]N,18所有行的数据填充。
4.根据权利要求1所述的复杂曲面电磁波矢量动态可视化方法,其特征在于:所述步骤S4包括以下子步骤:
S401.设定电磁波初始显示相位Phs0=0,由用户根据动态刷新速度设定电磁波相位步进Pstep,由慢到快到非常快依次设定为k1、k2、k3、k4、k5、k6六档;初始化maxAc=0;
其中,α为箭头长度显示调整参数,设置为D1、D2、D3、D4、D5、D6六档,分别表示用户设置箭头相对模型的不同的显示比例;Δs为面元网格的平均边长;为由步骤S404得到;max(Ac)由步骤S402得到;
S407.在i=1,2…,N时,循环步骤S403~步骤S406,绘制Phs0相位下的所有箭头;
S408.等待100ms,将Phs0+Pstep赋值给Phs0;
S409.重复执行步骤S403~S408,直到用户停止电磁波矢量动态绘制为止,绘制过程结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111401043.7A CN114067079B (zh) | 2021-11-19 | 2021-11-19 | 复杂曲面电磁波矢量动态可视化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111401043.7A CN114067079B (zh) | 2021-11-19 | 2021-11-19 | 复杂曲面电磁波矢量动态可视化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114067079A true CN114067079A (zh) | 2022-02-18 |
CN114067079B CN114067079B (zh) | 2022-05-13 |
Family
ID=80276684
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111401043.7A Active CN114067079B (zh) | 2021-11-19 | 2021-11-19 | 复杂曲面电磁波矢量动态可视化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114067079B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115546015A (zh) * | 2022-10-27 | 2022-12-30 | 东莘电磁科技(成都)有限公司 | 一种网格曲面受平面电磁波时谐感应场特征图像生成方法 |
CN116719981A (zh) * | 2023-08-08 | 2023-09-08 | 东莘电磁科技(成都)有限公司 | 一种基于电磁领域的细线网格流场可视化方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8723870B1 (en) * | 2012-01-30 | 2014-05-13 | Google Inc. | Selection of object types with data transferability |
US20190347837A1 (en) * | 2018-05-14 | 2019-11-14 | Virtualitics, Inc. | Systems and Methods for High Dimensional 3D Data Visualization |
CN110706304A (zh) * | 2019-09-23 | 2020-01-17 | 中国科学院海洋研究所 | 一种极地水汽输送通量的可视化算法 |
CN113010829A (zh) * | 2021-03-31 | 2021-06-22 | 建信金融科技有限责任公司 | 一种数据分区可视化方法、装置、计算机设备及存储介质 |
CN113035283A (zh) * | 2021-01-12 | 2021-06-25 | 南京理工大学 | 基于分子动力学的多层富勒烯单向压缩的模拟方法 |
CN113115021A (zh) * | 2021-04-09 | 2021-07-13 | 天海欧康科技信息(厦门)有限公司 | 一种物流三维可视化场景中的相机位动态聚焦方法 |
US11069010B1 (en) * | 2016-05-20 | 2021-07-20 | Trsoft | Methods and systems for creating single pass optimized property search portal |
KR20210100347A (ko) * | 2020-02-06 | 2021-08-17 | 한국전자통신연구원 | 가시화 플랫폼에서 사전 샘플을 생성하는 방법 |
CN113282553A (zh) * | 2021-05-24 | 2021-08-20 | 北京航空航天大学 | 一种电磁环境远场数据的快速无损压缩存储方法 |
CN113642057A (zh) * | 2021-07-13 | 2021-11-12 | 长江勘测规划设计研究有限责任公司 | 一种自动生成Autocad格式流速矢量图的方法 |
-
2021
- 2021-11-19 CN CN202111401043.7A patent/CN114067079B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8723870B1 (en) * | 2012-01-30 | 2014-05-13 | Google Inc. | Selection of object types with data transferability |
US11069010B1 (en) * | 2016-05-20 | 2021-07-20 | Trsoft | Methods and systems for creating single pass optimized property search portal |
US20190347837A1 (en) * | 2018-05-14 | 2019-11-14 | Virtualitics, Inc. | Systems and Methods for High Dimensional 3D Data Visualization |
CN110706304A (zh) * | 2019-09-23 | 2020-01-17 | 中国科学院海洋研究所 | 一种极地水汽输送通量的可视化算法 |
KR20210100347A (ko) * | 2020-02-06 | 2021-08-17 | 한국전자통신연구원 | 가시화 플랫폼에서 사전 샘플을 생성하는 방법 |
CN113035283A (zh) * | 2021-01-12 | 2021-06-25 | 南京理工大学 | 基于分子动力学的多层富勒烯单向压缩的模拟方法 |
CN113010829A (zh) * | 2021-03-31 | 2021-06-22 | 建信金融科技有限责任公司 | 一种数据分区可视化方法、装置、计算机设备及存储介质 |
CN113115021A (zh) * | 2021-04-09 | 2021-07-13 | 天海欧康科技信息(厦门)有限公司 | 一种物流三维可视化场景中的相机位动态聚焦方法 |
CN113282553A (zh) * | 2021-05-24 | 2021-08-20 | 北京航空航天大学 | 一种电磁环境远场数据的快速无损压缩存储方法 |
CN113642057A (zh) * | 2021-07-13 | 2021-11-12 | 长江勘测规划设计研究有限责任公司 | 一种自动生成Autocad格式流速矢量图的方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115546015A (zh) * | 2022-10-27 | 2022-12-30 | 东莘电磁科技(成都)有限公司 | 一种网格曲面受平面电磁波时谐感应场特征图像生成方法 |
CN115546015B (zh) * | 2022-10-27 | 2023-09-22 | 东莘电磁科技(成都)有限公司 | 一种网格曲面受平面电磁波时谐感应场特征图像生成方法 |
CN116719981A (zh) * | 2023-08-08 | 2023-09-08 | 东莘电磁科技(成都)有限公司 | 一种基于电磁领域的细线网格流场可视化方法 |
CN116719981B (zh) * | 2023-08-08 | 2023-10-27 | 东莘电磁科技(成都)有限公司 | 一种基于电磁领域的细线网格流场可视化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114067079B (zh) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114067079B (zh) | 复杂曲面电磁波矢量动态可视化方法 | |
US8089480B2 (en) | Method for meshing a curved surface | |
CN107462154A (zh) | 机器人末端工具的位姿测量方法 | |
JPH01149176A (ja) | 図形表示システムにおける曲面情報生成方法 | |
TW201439968A (zh) | 使用誤差項之保守格柵化基元 | |
CN105631939B (zh) | 一种基于曲率滤波的三维点云畸变校正方法及其系统 | |
Wan et al. | Mitigating sign problem by automatic differentiation | |
CN114049465A (zh) | 适用于复杂曲面表面的电磁波动态三维展示方法 | |
CN109358852A (zh) | 一种使用shader实现车载GUI换肤的方法 | |
JP2023123337A (ja) | 量子システムのシミュレーション方法、コンピューティングデバイス、コンピューティング装置、及びコンピュータプログラム | |
CN103020351A (zh) | 一种飞机位姿的三维实时显示方法 | |
CN114943772A (zh) | 一种增强现实视觉显示器中的相机标定方法、系统及应用 | |
CN101625765A (zh) | 一种旋转椭圆弧的绘制方法 | |
CN109177138A (zh) | 一种玻璃与膜片的对位方法及装置 | |
CN115101140B (zh) | 确定分子的基态特征的方法、设备及存储介质 | |
CN114067089B (zh) | 基于支撑点拟合的复杂曲线电磁波矢量动态显示方法 | |
CN112254681B (zh) | 一种发散型多线激光投影测量仿真系统及其实现方法 | |
CN113065251A (zh) | 一种强耦合波导传播振幅获取方法及装置 | |
CN109740234B (zh) | 一种家居布局3d碰撞检测方法及系统 | |
CN117092149B (zh) | 一种溶析结晶在线监测系统 | |
CN101633134B (zh) | 数控机床系统中输入二维图形来实现三维图形加工的方法 | |
Zhang et al. | Multiquadrics collocation method for transient eddy current problems | |
WO2023245511A1 (zh) | 工艺参数调整方法、系统、生产系统及计算机设备 | |
Zhen-wei et al. | Effective factors analysis on camera calibration | |
CN115660110B (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 |