1.一種深水浮筒平臺的自激振蕩分析方法，建立的自激振蕩分析模型如下：

(m+ma)u··+cu·+ku=12ρD(U-u·)[(U-u·)C‾D(α)+|U-u·|CD′]]]>

(m+ma)v··+(c+ca)v·+kv=12ρD(U-u·)2C‾L(α)]]>

式中：m--浮筒平臺質量；

??????m_a--附加質量；

??????c--系泊系統阻尼；

??????c_a--附加阻尼；

??????k--系泊系統剛度；

??????u--平臺順流向運動位移；

--平臺順流向運動速度；

??????ü--平臺順流向運動加速度；

??????v--平臺橫流向運動位移；

--平臺橫流向運動速度；

--平臺橫流向運動加速度；

??????ρ--流體密度；

??????D--平臺直徑；

??????U--流速；

分別為自激阻力系數和自激升力系數，計算公式如下：

C‾D(α)=[CL′tanα+CD′]secα]]>

C‾L(α)=[CL′+CD′tanα]secα]]>

式中，C′_L--升力系數；

??????C′_D--拖曳力系數；

??????α--攻角；

采用迭代方法，計算規定時長內平臺自激振蕩的位移、速度和加速度。

2.如權利要求1所述的深水浮筒平臺的自激振蕩分析方法，其特征在于：所述的采用迭代方法進行計算的具體步驟如下：

1)給定計算時間、平臺順流向和橫流向速度和位移的初值：

t_j＝t₀＝0，u·j(i)=u·0(0)=0,]]>uj(i)=u0(0)=0,]]>v·j(i)=v·0(0)=0,]]>vj(i)=v0(0)=0]]>

式中：j--時間步數，計算開始時j＝0；

??????i--迭代次數，每個時間步開始時i＝0；

2)計算攻角

αj(i)=v·j(i)U]]>

式中：--第j時間步內第i次迭代的攻角；

--第j時間步內第i次迭代的平臺橫流向速度；

U--流體流速；

3)將步驟2)得到的攻角代入下面公式計算自激阻力系數和自激升力系數：

C‾D(αj(i))=[CL′tanαj(i)+CD′]secαj(i)]]>

C‾L(αj(i))=[CL′+CD′tanαj(i)]secαj(i)]]>

式中，C′_L--升力系數；

??????C′_D--拖曳力系數；

4)將步驟3)得到的自激阻力系數和自激升力系數代入下面公式：

(m+ma)u··j(i+1)+cu·j(i+1)+kuj(i+1)=12ρD(U-u·j(i))[(U-u·j(i))C‾D(αj(i))+|U-u·j(i)|CD′]]]>

(m+ma)v··j(i+1)+(c+ca)v·j(i+1)+kvj(i+1)=12ρD(U-u·j(i))2C‾L(αj(i))]]>

式中：m--浮筒平臺質量；

??????m_a--附加質量；

??????c--系泊系統阻尼；

??????c_a--附加阻尼；

??????k--系泊系統剛度；

--第j時間步內第i+1次迭代的平臺順流向運動位移；

--第j時間步內第i+1次迭代的平臺順流向運動速度；

--第j時間步內第i+1次迭代的平臺順流向運動加速度；

--第j時間步內第i+1次迭代的平臺橫流向運動位移；

--第j時間步內第i+1次迭代的平臺橫流向運動速度；

--第j時間步內第i+1次迭代的平臺橫流向運動加速度；

??????ρ--流體密度；

??????D--平臺直徑；

??????U--流速；

計算第j時間步內第i+1次迭代的平臺自激振蕩位移速度和加速度

5)如果max{|u··j(i+1)-u··j(i)|,]]>|u·j(i+1)-u·j(i)|}>ϵ]]>或max{|v··j(i+1)-v··j(i)|,]]>|v·j(i+1)-v·j(i)|}>ϵ,]]>ε為預先設定的計算精度，則繼續進行迭代計算，令：

u·j(i)=u·j(i+1),]]>u··j(i)=u··j(i+1),]]>v·j(i)=v·j(i+1),]]>v··j(i)=v··j(i+1)]]>

然后，重復步驟2)～5)的計算；

否則開始下一個時間步的計算，令：

t_j＝t_j+1＝t_j+Δt，uj(i)=uj+1(0)=uj(m),]]>u·j(i)=u·j+1(0)=u·j(m),]]>u··j(i)=u··j+1(0)=u··j(m),]]>

vj(i)=vj+1(0)=vj(m),]]>v·j(i)=v·j+1(0)=v·j(m),]]>v··j(i)=v··j+1(0)=v··j(m),]]>

Δt為時間增量，m為第j時間步的最大迭代次數；

重復步驟2)～5)的計算，直至計算時長滿足需要。