[发明专利]一种轴对称物理场的计算方法

说明书

技术领域

本发明涉及一种物理场的计算方法，更具体的，涉及一种计算单个物理源作用下场值遵循轴对称分布的物理场的方法。

背景技术

目前，关于物理场分布的计算多采用有限元仿真的方法，通过计算机根据物理场公式来计算每个有限元的场值，这种实现过程往往需要大量指数运算和幂运算，计算效率低，严重影响其在具体功能应用中的实时性。例如在肿瘤微波介入治疗中，通过实验测量及组合方程式计算出微波辐射天线在不同辐射功率、不同天线类型、不同功率作用时间及不同癌变组织的血流灌注条件下的实时三维凝固热场分布，再根据不同条件的排列组合，建立计算机模型，利用有限元分析求出三维凝固热场的分布数据库。专利号为02116887.3的中国发明专利公开了类似的计算方法。这种方法需要计算生物传热方程和微波辐射能力SAR的分布函数方程，计算费时，并且建立的三维凝固热场分布数据库规模巨大，不方便存储及维护。

发明内容

本发明的目的在于提供一种可以快速计算出单个物理源作用下场值遵循轴对称分布的物理场分布场值的方法。

为了实现上述目的，本发明采取如下的技术解决方案：

一种轴对称物理场的计算方法，所述物理场的场值轴对称分布，步骤如下：

步骤一、根据物理场公式计算出不同的轴向z和半径r对应的场值；

步骤二、将计算得到的场值存储至数据库，得到以轴向z和半径r为索引的物理场分布数据库；

步骤三、根据物理源的三维空间方位，计算出空间任意点P的轴向z和半径r：

其中，为原点O到P点在对称轴上的投影P’的矢量，为P点到其在对称轴上的投影P’的矢量；

步骤四、根据轴向z和半径r从数据库中读取对应的物理场值，便可计算P点处单物理源作用下的物理场分布值。

进一步的，所述步骤三中的其中，为原点O到P点的矢量，为对称轴的正方向。

进一步的，所述物理场为微波辐射热场，所述步骤一中的物理场公式为最终平面分布函数SAR＝ρce^-αr[X₀Z³+X₁Z²+X₂Z+X₃]，根据最终平面分布函数计算不同的轴向z和半径r对应的微波辐射热场值。

由以上技术方案可知，本发明方法在离线状态下先根据物理场公式按照固定步进计算出z方向上不同半径r对应的场值，并将计算结果存储在计算机数据库中；然后在实际应用计算过程中，通过简单的几何计算轴向z和半径r作为索引值用于数据库的索引，减少了计算量，运算速度快，大大提升了物理场分布的计算速度，能达到快速计算物理场分布的效果，解决了物理场分布场值计算效率低的问题。

附图说明

为了更清楚地说明本发明实施例，下面将对实施例或现有技术描述中所需要使用的附图做简单介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为单个物理源三维空间场值索引计算示意图；

图2为微波天线三维空间场值索引计算示意。

具体实施方式

为了让本发明的上述和其它目的、特征及优点能更明显，下文特举本发明实施例，并配合所附图示，做详细说明如下。

如图1所示，对于单个物理源作用下场值遵循轴对称分布的物理场来说，在柱坐标系下，某一点P处的场值只与轴向z和半径r有关，与角度θ无关。本发明方法的基本思路是：首先根据物理场公式按照固定步进计算出z方向上不同半径r对应的场值；将计算得到的场值存储至数据库，形成以轴向z和半径r为索引的物理场分布数据库；根据已知物理源的三维空间方位，计算出空间中任意一点P处物理场值的索引：轴向和半径以轴向z和半径r作为索引从数据库中读取对应的物理场值，便可快速计算出特定位置处单物理源作用下的物理场分布值。

以上是本发明的核心思想，下面将结合本发明实施例中的附图，对本发明实施例的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下获得的所有其他实施例，都属于本发明保护的范围。