基于邊界打靶法的橫向流作用下懸臂梁靜態(tài)模態(tài)函數(shù)求解

王小文 周益嫻

摘? 要：該文以橫向流作用下的單根懸臂梁結(jié)構(gòu)為對(duì)象，研究其靜態(tài)偏折模態(tài)形狀函數(shù)，為后續(xù)反應(yīng)堆內(nèi)設(shè)備一階梁式振型下附加質(zhì)量的求解打下基礎(chǔ)。首先，采用大撓度非線性方程組進(jìn)行建模;其次，采用一種邊界打靶法來(lái)進(jìn)行上述非線性微分控制方程求解，得到懸臂梁的偏折曲線的數(shù)值解，并將數(shù)值計(jì)算結(jié)果和存在解析解的情況下的解析結(jié)果進(jìn)行對(duì)比。

關(guān)鍵詞：筒倉(cāng)卸載懸臂梁? ?靜態(tài)偏折模態(tài)形狀函數(shù)? ?邊界打靶法? ?梁式振型函數(shù)

中圖分類號(hào)：TU311? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1672-3791（2021）08（b）-0004-04

Solution of Static Deflection Mode Shape Function of Cantilever Beam under Transverse Flow Based on Boundary Shooting Method

WANG Xiaowen? ?ZHOU Yixian*

（Beijing Key Laboratory of Passive Nuclear Safety Technology， North China Electric Power University， Beijing， 102206 China）

Abstract： In this paper， it studies the cantilever beam structure under the action of transverse flow， and focus on the derivation of its static cantilever deflection mode shape function， in order to lay foundation for the calculation of hydrodynamic added mass for the beam mode vibration of nuclear reactor internal cantilever components. Firstly， a set of nonlinear differential equations for the analysis of the single cantilever beam has been proposed. Secondly， to solve the nonlinear differential equations， we used a boundary shooting framework. The results of the numerical simulation agree well with the analytical solution.

Key Words： Cantilever beam; Static deflection mode shape; Boundary shooting method; Beam mode function

反應(yīng)性事故發(fā)生后，堆內(nèi)構(gòu)件在流體中的振動(dòng)在反應(yīng)堆安全問(wèn)題中至關(guān)重要。堆內(nèi)構(gòu)件和高能流的接觸引發(fā)設(shè)備的振動(dòng)，進(jìn)而使得? ?浸沒(méi)于流動(dòng)流體中的設(shè)備，如設(shè)備支承筒、燃料棒或熱交換器管等發(fā)生疲勞，導(dǎo)致破損。隨著反應(yīng)堆規(guī)模的擴(kuò)大，堆內(nèi)構(gòu)件更易彎曲，堆內(nèi)冷卻劑流量增加，流致振動(dòng)問(wèn)題愈發(fā)嚴(yán)重。流致振動(dòng)問(wèn)題從1950年便開(kāi)始進(jìn)行研究，但由于該問(wèn)題涉及流固耦合和非線性科學(xué)，因此理論研究難度很大。目前，研究對(duì)象大多為不同堆型的燃料組件，其邊界條件一般為兩端固定，而對(duì)設(shè)備支承筒等懸臂結(jié)構(gòu)研究較少。且主要為實(shí)驗(yàn)和模擬研究，理論研究工作較少。實(shí)驗(yàn)和模擬方法主要針對(duì)特定的工程工況，所得結(jié)果精確，但耗時(shí)較長(zhǎng)。工程實(shí)踐中前期設(shè)計(jì)階段往往需要理論推導(dǎo)結(jié)果加以指導(dǎo)。

1? 理論分析

根據(jù)核電站內(nèi)構(gòu)件的工程特點(diǎn)，部分結(jié)構(gòu)可簡(jiǎn)化為橫向流作用下懸臂梁結(jié)構(gòu)振動(dòng)問(wèn)題[1-2]，由于高徑比較大，其一階梁式振型較易被激勵(lì)，有必要對(duì)一階梁式振型模態(tài)進(jìn)行單獨(dú)研究，且堆內(nèi)設(shè)備一階梁式振型下的附加質(zhì)量和頻率的求解需要知道其靜態(tài)偏折模態(tài)形狀函數(shù)，堆內(nèi)構(gòu)件一般可簡(jiǎn)化為具有一定剛度和彈性的梁[3-5]。因此，該文將以單根懸臂彈性梁結(jié)構(gòu)為研究對(duì)象，求解其偏折模態(tài)形狀函數(shù)。

考慮底端固定的懸臂梁，其他位置無(wú)任何支撐，并且有分布載荷作用于其上，懸臂梁的截面幾何形狀影響其剛度和彎曲機(jī)理[6-7]，我們假設(shè)其橫截面在垂直方向上不發(fā)生變化，最大彎矩和剪切應(yīng)力應(yīng)該位于懸臂梁的底部。假設(shè)只有一根懸臂梁完全暴露在水流中，而周圍沒(méi)有其他懸臂梁來(lái)改變垂直速度分布，懸臂梁受水流的拉力大小為：

式中，F(xiàn)D為流體對(duì)懸臂梁浸沒(méi)部分的拉力，CD為拉力系數(shù)，ρ為水的密度，u（z）為垂直方向的速度分布，A為物體浸入部分垂直于水流平面的投影面積。根據(jù)歐拉伯努利方程，該懸臂梁彎曲的微分控制方程如下：

（2）

（3）

（4）

式中，;;;。θ為偏折角度，s表示沿曲線的距離，M表示一個(gè)截面所受彎曲力矩，I表示慣性矩，E表示彈性模量，V表示作用在梁截面上的剪切應(yīng)力，w（z）表示水流對(duì)梁產(chǎn)生的分布載荷，x和z分別表示梁任意截面與其底部的水平距離和垂直距離。

公式（2）至公式（4） 控制懸臂梁彎折的角度，由于懸臂梁底端固定，因此有：

sinθ（z*=0）=0（5）

根據(jù)工程特點(diǎn)，考慮懸臂梁頂端不浸沒(méi)于水中的情況，因此露出水面位置懸臂梁所受彎矩和剪力為0，得到邊界條件：

M*（z*=1）=0（6）

V*（z*=1）=0（7）

2? 方法介紹

為了更好地利用邊界打靶法去求解上述的微分方程，我們引入了3個(gè)變量x1，x2，x3：

（8）

式（2）至式（4）改寫(xiě)為：

（9）

3個(gè)邊界條件為：

（10）

邊界打靶法是利用初值條件去代替已知邊界條件進(jìn)行求解，通過(guò)不斷改變初始值讓其逐漸逼近給定的邊界條件，因此調(diào)整邊界的公式選取十分重要，我們?cè)O(shè)定了兩個(gè)不斷調(diào)整初始條件的關(guān)系式公式（11）和公式（12）：

（11）

（12）

其中b2和b3分別代表代表頂端z*=1處x2和x3的邊界條件，α代表打靶系數(shù)，其值小于1，該值的選取取決于作用在懸臂梁的載荷大小。i表示的是迭代的步數(shù)。

3? 結(jié)果與分析

驗(yàn)證計(jì)算程序，我們采用Li等人[6]提供的數(shù)值進(jìn)行計(jì)算，并將結(jié)果與解析解進(jìn)行對(duì)比。懸臂梁和流場(chǎng)的參數(shù)為：抗彎剛度EI=2.5×10-7，懸臂梁尺寸D=2.4×10-5 m，水深h=0.05 m，水密度ρ=1×103 kg/m3，水流的平均速度，拉力系數(shù)CD=1.2。流場(chǎng)的分布滿足指數(shù)分布：

其中m=3.156，求解式（2）～（4）可以獲得沿懸臂梁的偏折角度、彎矩以及剪切力等分布情況。其解析解表達(dá)如下：

（14）

圖1給出了在指數(shù)分布流場(chǎng)中，懸臂梁的偏折曲線的數(shù)值解及解析解對(duì)比圖?？梢钥吹綌?shù)值解和解析解吻合較好，圖中最大相對(duì)誤差為3.02×10-6，可見(jiàn)計(jì)算精度很高，計(jì)算程序得到了驗(yàn)證。

接下來(lái)該文研究不同參數(shù)對(duì)懸臂梁彎折情況的影響。圖2（a）展示了不同抗彎剛度EI的懸臂梁的偏折情況，當(dāng)抗彎剛度大于5×10-7 N m2時(shí)，懸臂梁偏折很小或基本沒(méi)有偏折。而剛度較小時(shí)可產(chǎn)生較大的偏折。圖2（b）展示了不同來(lái)流平均速度u的懸臂梁的偏折情況。當(dāng)平均速度小于0.25 m/s時(shí)，懸臂梁偏折很小或基本沒(méi)有偏折，而速度較大時(shí)可產(chǎn)生較大的偏折。

4? 結(jié)語(yǔ)

該文以橫向流作用下的懸臂梁結(jié)構(gòu)為研究對(duì)象，采用理論分析和數(shù)值計(jì)算相結(jié)合方法，研究其偏折模態(tài)形狀函數(shù)。根據(jù)反應(yīng)堆內(nèi)組件的結(jié)構(gòu)特點(diǎn)，采用大撓度非線性方程組來(lái)進(jìn)行建模，論證了該模型處理懸臂梁偏折問(wèn)題的適用性。為求解該問(wèn)題，使用一種邊界打靶法來(lái)進(jìn)行上述非線性微分控制方程的求解，得到懸臂梁的偏折曲線的數(shù)值解，并將數(shù)值計(jì)算和解析解結(jié)果對(duì)比，二者吻合度極高，計(jì)算程序得到驗(yàn)證。后續(xù)的工作將深入研究迭代次數(shù)和打靶系數(shù)對(duì)結(jié)果的影響。

參考文獻(xiàn)

[1] FERRARI G，F(xiàn)RANCHINI G， BALASUBRAMANIAN P， et al.Nonlinear Vibrations of a Nuclear Fuel Rod Supported by Spacer Grids[J].Nuclear Engineering and Design，2020（361）：110503.

[2] GINEAU A，LONGATTE E，LUCOR D，et al.Macroscopic Model of Fluid Structure Interaction in Cylinder Arrangement using Theory of Mixture[J].Computers & Fluids，2020（202）：104499.

[3] DE SANTIS D，SHAMS A.Numerical Modeling of Flow Induced Vibration of Nuclear Fuel Rods[J].Nuclear Engineering and Design，2017（320）：44-56.

[4] DE SANTIS D，KOTTAPALLI S，SHAMS A.Numerical Simulations of Rod Assembly Vibration Induced by Turbulent Axial Flows[J].Nuclear Engineering and Design，2018（355）：94-105.

[5] DE SANTIS D，SHAMS A.An Advanced Numerical Framework for the Simulation of Flow Induced Vibration for Nuclear Applications[J].Annals of Nuclear Energy，2019（130）：218-231.

[6] LI C W， XIE J F.Numerical Modeling of Free Surface Flow over Submerged and Highly Flexible Vegetation[J].Advances in Water Resources，2011，34（4）：468-477.

[7] JING H F，CAI Y J，WANG W H，et al.Investigation of Open Channel Flow with Unsubmerged Rigid Vegetation by the Lattice Boltzmann Method[J].Journal of Hydrodynamics，2020（32）：771-783.