[发明专利]一种基于弹性变形的区域网格动态重构方法

权利要求书

1.一种基于弹性变形的区域网格动态重构方法，其特征在于，包括以下步骤：

步骤1、计算变形前内边界到外边界之间的距离；

步骤2、计算变形前空间网格点到物面的距离；

步骤3、计算空间网格点的分布情况；

步骤4、物面网格点更新；

步骤5、基于弹性形变完成区域空间网格点动态重构；

所述步骤1中变形前内边界到外边界之间的距离计算方法为：内边界到外边界之间的距离即为物面上的网格点沿法向到外边界点形成的线段长度,设为S(i,j_max,k)：

其中，i表示流向、j表示法向、k表示周向,(i,j,k)共同确定一个网格点，x(i,j_max,k)则表示该网格点在x方向的坐标值，y(i,j,k)表示该网格点在y方向的坐标值,z(i,j,k)表示该网格点在z方向的坐标值；

所述步骤2中计算变形前空间网格点到物面的距离计算方法为：空间网格点到物面的距离即为物面上的网格点沿法向到任意空间网格点之间的距离，设为S(i,j,k)：

所述步骤3中空间网格点的分布情况计算方法为：dS(i,j,k)＝S(i,j,k)/S(i,j_max,k)；

所述步骤4中物面网格点通过外部输入舵偏的变形角度进行的更新；

所述步骤5中的动态重构具体方法为：在物面网格点更新之后，外边界保持不动，对内外边界构成的区域内网格点进行位置调整，基于弹性拉伸保持网格点在网格线中的相对位置不变，调整后的位置为：

x′(i,j,k)＝dS(i,j,k)*(x(i,j_max,k)-x′(i,1,k))+x′(i,1,k)

y′(i,j,k)＝dS(i,j,k)*(y(i,j_max,k)-y′(i,1,k))+y′(i,1,k)

z′(i,j,k)＝dS(i,j,k)*(z(i,j_max,k)-z′(i,1,k))+z′(i,1,k)。