[发明专利]一种不同计算模型间的界面数据传递与插值方法

权利要求书

1.一种不同计算模型间的界面数据传递与插值方法，所述的数据传递是从模型A将界面数据传递到模型B，其特征在于，方法步骤如下：

步骤1：确定不同计算模型A与模型B之间的耦合几何界面；

所述的模型A与模型B既是不同物理域的计算模型，耦合界面即为物理域分界面；或是相同物理域不同结构尺度的模型，耦合界面即为局部模型的边界；

步骤2：提取模型A离散界面AS上的信息，并构成相应的外部数据库D^A；

所述的模型A离散界面AS上的信息包括单元、单元连接的节点编号、各节点的空间坐标与各变量的计算结果；所述的外部数据库包含三个，即模型A离散界面AS上的单元及其内部节点信息数据库，模型A离散界面AS上的节点及其空间坐标信息数据库，模型A离散界面AS上的节点在所有时刻下，所需传递的计算变量的数据库；所述的外部数据库，应采用相应软件便于调取的格式；

步骤3：依据问题类型，提取模型B离散界面BS上的信息，并构成相应的外部数据库D^B；

所述的模型B离散界面BS上的信息包括节点或积分点、各节点或积分点的空间坐标；所述的外部数据库为，模型B离散界面BS上的节点或积分点及其空间坐标信息数据库；所述的外部数据库，应采用相应软件便于调取的格式；

步骤4：针对模型B界面BS上的某一点P，在外部数据库D^A中寻找与该点空间位置满足相近规则的模型A界面AS上的网格集合E；

所述的相近规则是指模型A界面AS上的某一单元的中心与模型B界面BS上的某一点P在空间位置上差异在一定范围内；

步骤5：将模型B界面BS上一点P投影到网格集合E中的单元得到点P’，确定其局部坐标；

所述的局部坐标所在的坐标系以单元某一点为原点，与原点相连的某一条边作为X轴，以单元所在平面作为XY面而确定；

步骤6：利用局部坐标，在网格集合E中查找与点P匹配的单元，并构成匹配数据库；

在网格集合E中查找与点P匹配的单元，采用投影方法，判定投影后的点是否落在当前单元的内部而确定；所述的匹配数据库，应包含模型B界面BS上一点P所匹配的单元、及该单元所连接的节点、及点P在该单元内的面积坐标，并采用矩阵形式存储；

步骤7：调用匹配数据库及模型A计算结果数据库，将模型A某一时刻t的计算结果，传递到模型B界面BS某一点，得到插值后数据；

所述的传递到模型B界面BS某一点，是采用面积坐标结合形函数插值，计算得到模型B界面BS上一点处的数据；

步骤8：若所需得到的载荷f与模型B界面BS上的求解结果相关，则利用插值结果与模型B的实时计算结果，实现相应载荷f的插值；

所述的载荷，若不仅仅是由模型A界面AS上传递到模型B界面BS上的计算结果，也与模型B界面BS上的求解结果相关，则利用相应的关系，实时确定真实载荷；

步骤9：重复步骤7-8，在某一时刻t下实现由模型A界面AS向模型B界面BS上每个节点或积分点的数据传递；

步骤10：重复步骤9，实现所有时刻下，由模型A界面AS向模型B界面BS上每个节点或积分点的数据传递。

2.根据权利要求1所述的一种不同计算模型间的界面数据传递与插值方法，其特征在于，步骤1：依据模型A与模型B的单位制差异，两者之间的耦合几何界面在空间上具有等比例相似性，通过恰当的单位制转换，即缩放变换，在相同的坐标系SYS0下能够重合。

3.根据权利要求1所述的一种不同计算模型间的界面数据传递与插值方法，其特征在于，

步骤2：以AS采用三角形网格离散，所提取的数据包括离散化的网格与节点信息，不同时刻下节点的计算结果信息；所提取的网格信息采用矩阵形式存储，格式为:

其中，表示界面AS上第i个网格，表示单元所对应的节点编号，为界面AS上所有单元个数，代表由单元构成的数据库，所提取的节点信息采用矩阵形式存储，格式为：

其中，表示界面AS第i个节点，表示节点的几何坐标，为界面AS上所有节点个数，代表由节点构成的数据库，所提取的计算结果采用矩阵形式存储，格式为：

其中，r_i^A,j|_t表示节点处的第j个计算结果变量在t时刻的取值，为某个节点上所需传递计算结果变量的个数，代表所分析时刻的总数，代表由计算结果构成的数据库，由上述信息构成外部数据库：

步骤3：对模型B的界面进行离散得到BS，提取BS上的节点或积分点信息并变换到与AS相同的单位尺度，存储到外部数据库，所提取的节点信息采用矩阵形式存储，格式为：

其中，表示界面BS第i个节点，表示节点变换后的几何坐标，为界面BS上所有节点个数，代表由节点构成的数据库，所提取的积分点信息采用矩阵形式存储，格式为：

其中，表示界面BS第i个积分点，表示积分点变换后的几何坐标，为界面BS上所有节点个数，代表由积分点构成的数据库；

步骤4：针对模型B界面BS上的某一点或利用外部数据库寻找与该点空间位置满足相近规则的模型A界面AS上的网格集合E，所述的相近规则为：

或

其中L为的特征长度，k取1-2，M为单元所连接节点个数，三角形单元M＝3，四边形单元M＝4；

步骤5：将或投影到网格集合E中的单元得到点P’，确定其局部坐标，利用单元中的一个节点为原点O’，以单元作为XY平面，以为X轴构建局部直角坐标系SYS1；将的坐标或的坐标变换到SYS1中，获取其局部坐标，

或

其中，分别为SYS1的X、Y、Z在全局坐标SYS0下的方向向量，即：

步骤6：利用局部坐标，在网格集合E中查找与或匹配的单元并构成匹配数据库，在局部坐标系下，求解由P’与单元的节点构成的三角形面积，进而得到点P’在当前单元中的面积坐标N_i，分别为：

有：

则判定或与单元匹配，ε＝0.0001，所述的匹配数据库包含或及其对应匹配的单元编号单元所包含的节点，及含或在单元内的面积坐标：

其中上标^表示与或匹配；

步骤7：调用匹配数据库及模型A计算结果数据库，将模型A某一时刻t的计算结果，传递到模型B界面BS上或上，得到插值后数据，t时刻下变量j在模型B界面BS上或的插值数据记为其计算方法为：

步骤8：所需得到的载荷与模型B界面BS上的求解结果相关，CFD在计算热流中采用等温壁面，设壁面为常温，所计算得到的结果为冷壁热流，需要采用壁面温度实时计算结果进行修正，则利用插值结果与模型B的实时计算结果r_i^B,j|_t，实现相应载荷f的插值，即：

步骤9：重复步骤7-8，在某一时刻t下实现由模型A界面AS向模型B界面BS上每个节点或积分点的数据传递；

步骤10：重复步骤9，实现所有时刻下，由模型A界面AS向模型B界面BS上每个节点或积分点的数据传递。