[发明专利]一种模拟桥墩周期受力的数值方法

权利要求书

1.一种模拟桥墩周期受力的数值方法，其特征在于在浅水流动的动量方程中加入两种不同形式的力，分别是涡量在变化梯度方向的螺旋力和涡量在变化梯度方向的粘性耗散力其形式分别为

${\stackrel{\→}{B}}_1={\mathrm{\ϵ}}_1\▿\×\left(\frac{\mathrm{\φ}}{|\▿\mathrm{\φ}|}\stackrel{\→}{\mathrm{\ω}}\right);$

${\stackrel{\→}{B}}_2={\mathrm{\ϵ}}_2\frac{\mathrm{\φ}}{|\▿\mathrm{\φ}|}\left({\▿}^2\stackrel{\→}{V}\right),$

其中，代表涡量；φ是涡量的模；是涡量的模的梯度的模；是速度矢量；符号×是叉乘算子；是梯度算子；是拉普拉斯算子；ε₁是的放大系数；ε₂是的放大系数。

2.根据权利要求1所述的一种模拟桥墩周期受力的数值方法，其中所述的涡量在变化梯度方向的螺旋力的形式通过高斯定理，将计算网格单元内的力的积分形式转变为计算网格表面通量的积分形式，并且保留在动量方程等号的左边。

3.根据权利要求1所述的一种模拟桥墩周期受力的数值方法，其中所述的涡量在变化梯度方向的粘性耗散力的形式保留在动量方程右边的源项位置。