[发明专利]基于暴流强度的倾斜微下击暴流建模方法

权利要求书

1.一种基于暴流强度的倾斜微下击暴流建模方法，其特征在于：所述方法通过以下步 骤实现：

步骤一、建立倾斜涡环模型；

(一)、定义：

在距离水平地面±H的高度上布置两个对称的涡环，涡环强度为Γ，且水平地面以上为 正值，水平地面以下为负值，涡环半径为R，地面上方的涡环称为主涡环，与之对称的地面下 方的涡环称为镜像涡环；其中，主涡环和镜像涡环中心处称为涡丝；

(二)、建立垂直涡环模型：

中心轴与水平地面垂直时诱导出开口垂直向下的微下击暴流场，单一涡环在空间任一 任一点(x,y,z)处产生的诱导速度V_i＝[v_x,v_y,v_z]可由单一涡环的流线方程Ψ求得；

vx=(xr2)∂ψ∂zvy=(yr2)∂ψ∂zvz=(-1r)∂ψ∂r---(1)]]>

式中，为空间一点到涡环中心的距离；v_x,，v_y，v_z分别为涡环模型在任一点 (x,y,z)处诱导速度的三个分量；为偏微分运算；Ψ为单一涡环的流线方程，定义为：

ψ≈1.576Γπ[(r1+r2)λa21+31-λa2-(r1′+r2′)λs21+31-λs2]---(2)]]>

式(2)中，

r1=(z-H)2+(r-R)2---(3)]]>

r2=(z-H)2+(r+R)2---(4)]]>

r1′=(z+H)2+(r-R)2---(5)]]>

r2′=(z+H)2+(r+R)2---(6)]]>

λa=r2-r1r2+r1---(7)]]>

λs=r2′-r1′r2′+r1′---(8)]]>

其中，Γ为涡环强度，r₁为空间任一点(x,y,z)到主涡环的最小距离，r₂为空间任一点 (x,y,z)到主涡环的最大距离，r₁'为空间任一点(x,y,z)到镜像涡环的最小距离，r₂'为空间 任一点(x,y,z)到镜像涡环的最大距离；

(三)、将(二)中的垂直涡环模型按照倾斜角θ与朝向角进行旋转得到倾斜涡环模型：

(a)、设倾斜涡环模型中心轴过地面坐标系原点，倾斜涡环模型的主涡环的中心轴与地 面坐标系的z轴的夹角为θ，角θ范围为[0°,90°]，表示诱导出的微下击暴流场的倾斜程度， 称之为倾斜角，其中，0°表示主涡环的中心轴垂直于地面，90°表示主涡环的中心轴为水平；

x轴负半轴与中心轴在xoy平面的投影形成夹角为角范围为[0°,360°]，表示诱导出 的微下击暴流场的朝向，称之为朝向角；其中0°表示中心轴投影在y轴正半轴，90°表示投影 在x轴负半轴，180°表示投影在y轴负半轴，270°表示投影在x轴正半轴；

相应的，镜像涡环与主涡环关于xoy平面对称，则中心轴与z轴的夹角为-θ，x轴负半轴 与中心轴在xoy平面的投影形成夹角为这两个旋转角满足中心轴在空间中的任意点旋 转；

(b)、引入旋转矩阵和L(θ)实现主涡环模型转换，引入和L(-θ)实现镜像涡环 模型转换；

最终生成的倾斜涡环模型诱导风速的矢量形式为：

式中

L(θ)=cosθ0sinθ010-sinθ0cosθ---(10)]]>

v_sx,v_sy,v_sz为倾斜涡环模型在任一点(x,y,z)处诱导速度的三个分量，v_xa、v_ya、v_za为主 涡环诱导速度的三个分量，v_xsv_ysv_zs为镜像涡环诱导速度的三个分量；

(c)、旋转过程为：主涡环中心轴先绕y轴旋转θ角，再绕z轴旋转角，方向与右手螺旋相 反；相对的，镜像涡环旋转角为-θ、且通过令风速矩阵先乘倾斜角旋转矩阵，再乘朝向角 旋转矩阵的顺序，体现先作倾斜旋转再作朝向旋转的操作，通过倾斜角和朝向角完成仿真；

步骤二、根据微下击暴流场的强度对倾斜多涡环模型的参数{H,R,Γ}、倾斜角θ与朝向 角进行选择：

微下击暴流的强度是由其最大速度决定的，设待生成微下击暴流场的最大速度为V_max， 通过嵌套微分进化算法确定倾斜多涡环模型参数M＝{{Γ_i},{R_i},{H_i}}，i＝1,…,n，n为涡 环对总数，下式成立：

maxv(M,x,y,z)＝V_max(12)

v(M,x,y,z)为选定倾斜多涡环模型参数M后，微下击暴流模型在三维空间坐标为(x,y, z)处的综合诱导速度，用下式表示：

v=vx2+vy2+vz2=(Σi=1nvix)2+(Σi=1nviy)2+(Σi=1nviz)2---(13)]]>

其中，v_x、v_y和v_z为多涡环微下击暴流模型在任一点(x,y,z)处诱导速度的三个分量； v_ix、v_iy和v_iz为第i对涡环在任一点(x,y,z)处诱导速度的三个分量。