針對本科畢設(shè)中所涉及到的離心泵數(shù)值分析和性能計(jì)算，將用最簡單粗暴的方法，講解如何基于CFturbo、ICEM、FLUENT來開展離心泵水力設(shè)計(jì)和性能分析的計(jì)算機(jī)輔助（CAE）實(shí)現(xiàn)。離心泵的水力設(shè)計(jì)由CFturbo軟件實(shí)現(xiàn)；網(wǎng)格剖分由ICEM軟件實(shí)現(xiàn)；CFD數(shù)值計(jì)算由FLUENT軟件實(shí)現(xiàn)；并驗(yàn)證設(shè)計(jì)值是否達(dá)到。

這里是第三部分，FLUENT軟件實(shí)現(xiàn)泵內(nèi)流動(dòng)的數(shù)值模擬，并對模擬結(jié)果進(jìn)行分析，驗(yàn)證設(shè)計(jì)值是否達(dá)到參數(shù)要求。

1 讀入網(wǎng)格

先新建個(gè)文件夾作為FLUENT的工作目錄，并把之前畫好的三個(gè)網(wǎng)格文件fluent_in.msh、fluent_imp.msh、fluent_vol.msh全部拷進(jìn)去。

打開FLUENT軟件，在Fluent Launcher中設(shè)置Dimensions維數(shù)選3D三維，因?yàn)樵蹅兪侨S網(wǎng)格；

Options勾選Double Precision雙精度，現(xiàn)在也沒誰拿單精度Real算了，畢竟單精度計(jì)算誤差太大了；

Processiong Options處理器選項(xiàng)，可選Serial單核串行和Parallel多核并行，我這破機(jī)子，也只能Serial單核計(jì)算了；

在Working Directory工作目錄中指定到剛才新建的文件夾，即可確定FLUENT的工作目錄跟它相關(guān)了。

再次確認(rèn)設(shè)置參數(shù)后，OK確定打開FLUENT。

讀入網(wǎng)格，File->Read->Mesh，在打開的工作目錄中，可發(fā)現(xiàn)有三個(gè)剛才放進(jìn)去的網(wǎng)格文件，先打開fluent_imp.msh葉輪流道網(wǎng)格文件。

讀入完畢后，點(diǎn)左上方Display查看下網(wǎng)格是否正確。

在Mesh Display窗口，可選擇Surfaces中的面來逐個(gè)查看是否有錯(cuò)誤，注意別選int_body_imp，這是內(nèi)部區(qū)域的網(wǎng)格面，選上的話數(shù)量非常大，顯示起來極有可能把機(jī)子卡死，當(dāng)然你機(jī)子好的話，就當(dāng)我沒說。

再把進(jìn)口延長段的網(wǎng)格讀進(jìn)來，只是這次不能再用read mesh了，要用Setting Up Domain中的Zones中的Append下的Append Case File…，即再附加上另一個(gè)區(qū)域的方式來讀入新的網(wǎng)格。

選擇fluent_in.msh文件，OK讀入。

同樣的方式，把fluent_vol.msh也添加進(jìn)來，Display一下看看，確定各個(gè)面都是正確的。

2 縮放網(wǎng)格

ICEM輸出的網(wǎng)格長度單位是毫米mm，FLUENT讀入網(wǎng)格的長度單位是m，兩者差了1000倍，所以要進(jìn)行縮放的操作，不然的話，雷諾數(shù)就差了1000倍，結(jié)果完全不靠譜了。

Setting Up Domain -> Mesh -> Transform -> Scale…，打開網(wǎng)格縮放窗口。

在Scale Mesh窗口中，選擇Mesh Was Created In網(wǎng)格是用mm毫米創(chuàng)建的，然后點(diǎn)擊Scale來縮放，可看到左上方的區(qū)域范圍的值都變小了1000倍，縮放成功，Close關(guān)閉Scale Mesh窗口即可。

Setting Up Domain -> Mesh -> Info -> Size，查看網(wǎng)格信息，大概50萬節(jié)點(diǎn)不到，個(gè)人機(jī)還是可以跑跑的。

Mesh SizeLevel Cells Faces Nodes Partitions0 2655266 5375651 476076 13 cell zones, 21 face zones.

Setting Up Domain -> Mesh -> Check，檢查網(wǎng)格，沒有負(fù)體積，計(jì)算域范圍也縮放正確了，OK。

Domain Extents:x-coordinate: min (m) = -1.900081e-01, max (m) = 7.313037e-01y-coordinate: min (m) = -1.773689e-01, max (m) = 2.250622e-01z-coordinate: min (m) = -2.579141e-01, max (m) = 1.051617e-01Volume statistics:minimum volume (m3): 3.716923e-11maximum volume (m3): 3.629198e-08total volume (m3): 1.217895e-02Face area statistics:minimum face area (m2): 4.020561e-08maximum face area (m2): 3.415682e-05Checking mesh......................... Done.

Setting Up Domain -> Mesh -> Quality，檢查網(wǎng)格質(zhì)量，質(zhì)量確實(shí)不咋滴，都逼近下限了，咱也就大概算算，不深究了。

Mesh Quality:Minimum Orthogonal Quality = 1.05554e-02 (Orthogonal Quality ranges from 0 to 1, where values close to 0 correspond to low quality.)Maximum Ortho Skew = 9.87453e-01 (Ortho Skew ranges from 0 to 1, where values close to 1 correspond to low quality.)Maximum Aspect Ratio = 3.68410e+02

接下來，基本上照著左側(cè)的樹形菜單挨個(gè)整下來，就順序地完成了設(shè)置、計(jì)算和分析工作。

3 一般設(shè)置

左欄設(shè)置樹，Tree -> Setup -> General，一般設(shè)置，保持默認(rèn)的Pressure-Based壓力基求解器（適用于不可壓縮流動(dòng)問題，Density-Based密度基求解器則適用于可壓縮問題），Time中保持默認(rèn)的Steady定常問題（非定常計(jì)算會(huì)非常消耗資源，這里就是大概算算），Gravity重立項(xiàng)也不予考慮。

左欄設(shè)置樹，Tree -> Setup -> Models -> Viscous，打開粘性模型設(shè)置，把默認(rèn)的Laminar層流模型改成k-epsilon(2 eqn) $k??$兩方程湍流模型，該湍流模型的設(shè)置保持默認(rèn)即可，OK。

左欄設(shè)置樹，Tree -> Steup -> Materials，設(shè)置介質(zhì)，點(diǎn)開Fluid下面的air，再從Create/Edit Materials中打開Fluent Database…，從Fluent Database Materials中選擇water-liquid (h2o<l>)，Copy到當(dāng)前介質(zhì)中，Close關(guān)閉Fluent Database Materials窗口，檢查Create/Edit Materials窗口中是否改成了water水，注意水的密度是998.2 kg/m3，單擊Change/Create確認(rèn)更改，Close關(guān)閉Crete/Edit Materials窗口，看到Materials中Fluid下多了個(gè)water-liquid。

4 區(qū)域設(shè)置

該算例一共有三個(gè)區(qū)域，進(jìn)口延伸段、葉輪流道、蝸殼流道，咱們要把進(jìn)口延伸段和蝸殼流道定義成靜止區(qū)域，而把葉輪流道定義成旋轉(zhuǎn)區(qū)域，在區(qū)域設(shè)置中完成，同時(shí)還要指定區(qū)域中的流體介質(zhì)。

左欄設(shè)置樹，Tree -> Setup -> Cell Zone Conditions，設(shè)置區(qū)域條件，雙擊body_in設(shè)置進(jìn)口延伸區(qū)域，把Material Name從默認(rèn)的air空氣改成wate-liquid液態(tài)水，其余默認(rèn)不做任何勾選和修改，OK，表明該區(qū)域是靜止區(qū)域，且內(nèi)部介質(zhì)是液態(tài)水。

同樣設(shè)置body_vol蝸殼流道區(qū)域，其設(shè)置和body_in區(qū)域一樣。

重頭戲來了，設(shè)置葉輪區(qū)域，雙擊body_imp，打開設(shè)置窗口，Material Name同樣把a(bǔ)ir改成water-liquid，勾選Frame Motion運(yùn)動(dòng)參考系，Rotation-Axis Origin旋轉(zhuǎn)軸心為坐標(biāo)原點(diǎn)默認(rèn)不動(dòng)，Rotation-Axis Direction旋轉(zhuǎn)軸方向默認(rèn)z軸也不更改，Rotational Velocity旋轉(zhuǎn)速度改為-2900rpm，負(fù)號表示它是繞z軸順時(shí)針旋轉(zhuǎn)的，（旋轉(zhuǎn)單位默認(rèn)是rad/s，可以在最上端的Setting Up Domain -> Mesh -> Units中將angular-velocity的單位Units從rad/s改為rpm），Translational Velocity平動(dòng)速度保持默認(rèn)的0就好，咱們這只有轉(zhuǎn)動(dòng)，沒有平動(dòng)。檢查無誤后OK確認(rèn)。

5 邊界設(shè)置

上一步設(shè)置的是區(qū)域條件，接下來設(shè)置邊界的條件。

左欄設(shè)置樹，Tree -> Setup -> Boundary Conditions，設(shè)置邊界條件。

5.1 旋轉(zhuǎn)壁面邊界

雙擊imp_blade_leading，設(shè)置葉片前緣面的邊界條件。

將Wall Motion壁面運(yùn)動(dòng)類型從Stationary Wall靜止壁面改為Moving Wall運(yùn)動(dòng)壁面；

Motion運(yùn)動(dòng)類型保持默認(rèn)的Relative to Adjacent Cell Zone相對于鄰近的單元區(qū)域，即給定的是其相對于其所屬的區(qū)域的運(yùn)動(dòng)（相對運(yùn)動(dòng)），它的Adjacent Cell Zone鄰近區(qū)域在最上面第2行已經(jīng)給你標(biāo)出來了是body_imp葉輪流道區(qū)域，運(yùn)動(dòng)類型把默認(rèn)的Translational平動(dòng)改成Rotational旋轉(zhuǎn)運(yùn)動(dòng)，轉(zhuǎn)速Speed中保持默認(rèn)的0，轉(zhuǎn)軸中心和轉(zhuǎn)軸依舊分別是原點(diǎn)和z軸。

Shear Conditions剪切條件默認(rèn)的No Slip無滑移壁面。

OK確定設(shè)置。

經(jīng)過這么一番折騰，就使得葉片的前緣面相對于葉輪流道區(qū)域做轉(zhuǎn)速是0的相對旋轉(zhuǎn)運(yùn)動(dòng)，換句話說，葉片此時(shí)和葉輪流道一樣做的是-2900rpm的絕對旋轉(zhuǎn)運(yùn)動(dòng)。

如果不做這個(gè)設(shè)置，那么葉輪前緣是固定不動(dòng)的壁面，這和實(shí)際情況是不符的。

同樣的，葉輪流道區(qū)域的葉片其他面、前后蓋板也要設(shè)置成這種相對轉(zhuǎn)速是0的運(yùn)動(dòng)壁面，如果一個(gè)個(gè)設(shè)置就太麻煩了，好在可以直接Copy的。

看到下面那個(gè)Copy…了么，點(diǎn)開來，在Copy Conditions窗口中，左邊選剛才設(shè)好的imp_blade_leading邊界，右邊選擇imp_blade_pressure、imp_blade_suction、imp_blade_tip、imp_blade_hub、imp_blade_shroud邊界，然后點(diǎn)最下面Copy完成復(fù)制邊界條件的操作。

如果不放心，再點(diǎn)開來看看它們設(shè)好了沒有吧。

5.2 無滑移壁面設(shè)置

還記得畫網(wǎng)格的時(shí)候講的葉輪流道前后蓋板延伸的環(huán)面，以及蝸殼進(jìn)口部分的側(cè)向延伸環(huán)面和徑向環(huán)面么？它們應(yīng)該如何設(shè)置呢？

實(shí)際上，它們也屬于壁面，即限制流體不能穿透，但是它們又不能給流體施加任何的切向力，所以，把它們設(shè)置成剪切力是零的滑移壁面就可以了，當(dāng)然，葉輪流道前后蓋板的延伸面可以定義成是旋轉(zhuǎn)面，也可以定義成是靜止面，沒有影響，因?yàn)樗鼈儔焊筒粫?huì)產(chǎn)生剪切應(yīng)力，也就不會(huì)對流體施加切向載荷的影響了。

只需要把Shear Condition剪切條件從No Slip無滑移改成Specified Shear指定剪切，同時(shí)保持Shear Stress剪切力為0，就可以了。

imp_ext_hub和imp_ext_shroud葉輪前后蓋板延伸環(huán)面的設(shè)置如下。

vol_ext和vol_side蝸殼流道進(jìn)口周向延伸面和徑向延伸面的設(shè)置：

5.3 靜止壁面

進(jìn)口延伸段的外壁面，蝸殼的壁面都是靜止壁面，保持默認(rèn)的靜止壁面和無滑移條件設(shè)置就好。

in_walls和vol_wall的設(shè)置。

5.4 動(dòng)靜交界面設(shè)置和匹配

進(jìn)口延伸段的出口面，與，葉輪進(jìn)口面，前者靜止，后者運(yùn)動(dòng)，屬于內(nèi)部的動(dòng)靜交界面。

葉輪出口面，與，蝸殼進(jìn)口面，前者運(yùn)動(dòng)，后者靜止，屬于內(nèi)部的動(dòng)靜交界面。

先指定這四個(gè)面的類型為交界面。

先選中in_outlet進(jìn)口延伸段的出口面，然后從Type類型中下來選擇interface交界面并單擊確認(rèn)，它會(huì)讓你重命名，保持默認(rèn)名字就好，OK確認(rèn)，就把in_outlet面的類型從wall壁面修改成了interface交界面。

同樣的方法，把vol_inlet蝸殼進(jìn)口面、imp_inlet葉輪進(jìn)口面、imp_outlet葉輪出口面也設(shè)置成interface交界面類型。

雖然設(shè)好了它們是交界面，但是并沒有將它們匹配起來啊，好辦。

是不是設(shè)置樹，Tree中Setup中多了一個(gè)Mesh Interfaces網(wǎng)格界面的選項(xiàng)啊，沒錯(cuò)，就是它了，單擊Create/Edit…打開交界面匹配創(chuàng)建窗口。

在Create/Edit Mesh Interface窗口中，Mesh Interface中輸入名字“in_2_imp”表示進(jìn)口段到葉輪段的交界面，在Interface Zone Side 1中選擇第一個(gè)交界面in_outlet，在Interface Zone Side 2中選擇第二個(gè)交界面imp_inlet，單擊最下面Create創(chuàng)建，這樣便將進(jìn)口段出口和葉輪段進(jìn)口匹配為了一對交界面。

同樣創(chuàng)建另一個(gè)名為"imp_2_vol"葉輪到蝸殼的交界面，把imp_outlet和vol_inlet配對即可。

設(shè)置好后，Close關(guān)閉Create/Edit Mesh Interface窗口。

5.5 進(jìn)出口邊界條件設(shè)置

只剩下兩個(gè)邊界面，in_inlet進(jìn)口段的進(jìn)口面和vol_outlet蝸殼流道的出口面了，這兩個(gè)面也對應(yīng)著整個(gè)計(jì)算域的進(jìn)口邊界和出口邊界。

通常是把進(jìn)口指定總壓，出口指定流量。進(jìn)口總壓比較好理解，比如接個(gè)大氣管道什么的，那就是大氣壓了對吧，出口流量有點(diǎn)費(fèi)解，其實(shí)出口一般接了一堆亂七八糟的管道，那么會(huì)有阻力什么的，而做實(shí)驗(yàn)的時(shí)候呢，會(huì)用個(gè)節(jié)流閥什么的去調(diào)節(jié)出口流量，也就是說，出口流量實(shí)際上是由使用場所或者人為來調(diào)節(jié)的，所以是可以直接指定的（當(dāng)然了，如果指定一個(gè)非常高的不切實(shí)際的流量，那么算出來的結(jié)果也是不符合物理意義的），我們指定了這個(gè)流量，再去看能否達(dá)到所需要的揚(yáng)程就可以了，就跟自來水龍頭，可以調(diào)節(jié)流量是一樣的道理了。

將in_inlet的類型從wall改成pressure-inlet壓力進(jìn)口，設(shè)置其Guage Total Pressure總壓為0（這個(gè)是相對表壓），Supersonic/Initial Gauge Pressure初始靜壓為0（只是初始值而已，一旦計(jì)算過程中，這個(gè)面上有了速度，那么會(huì)把速度的動(dòng)壓減去來計(jì)算靜壓值，最終是保證動(dòng)壓和靜壓加起來的總壓值是0就OK了，所以算起來之后靜壓會(huì)成負(fù)值的，表明是進(jìn)口存在負(fù)壓差）。

Turbulence湍流模型的設(shè)置中，將設(shè)置方式改為Intensity and Hydraulic Diameter湍流強(qiáng)度和水力直徑，湍流強(qiáng)度給成5%（湍流度中等），水力直徑給成0.1259m，即進(jìn)口管道的直徑即可。

OK確定。

將vol_outlet的類型從wall改成velocity-inlet速度進(jìn)口（沒辦法，這里沒有流量出口，也沒有速度出口，就用速度進(jìn)口，把速度指定成負(fù)值來給流量算了，效果是一樣的）。

先得算下速度值是多少，出口圓截面的直徑為102.46mm，面積是$\pi d^2/4=0.00824515m^2$，泵的設(shè)計(jì)流量是$200m^3/h=0.0555556m^3/s$，出口速度應(yīng)該是$0.0555556/0.00824515=6.737967m/s$，別忘了加負(fù)號。

即把Velocity Magnitude速度幅值設(shè)為-6.737967m/s，初始背壓設(shè)為0Pa，這個(gè)初始背壓同樣只是初始值，會(huì)在計(jì)算過程中自動(dòng)做調(diào)整的。

湍流度同樣用強(qiáng)度和水力直徑給定，湍流強(qiáng)度指定為10%（出口一般湍流會(huì)更強(qiáng)烈一些的，10%依舊是中等湍流度并不是特別大），水力直徑為0.10246m為管道直徑。

OK確定。

至此，設(shè)置樹的Setup部分就設(shè)置完了。

6 求解設(shè)置和迭代求解

設(shè)置樹，Tree -> Solution -> Solution Methods，保持默認(rèn)求解方法即可。

設(shè)置樹，Tree -> Solution -> Solution Controls，保持默認(rèn)松弛因子即可。

設(shè)置樹，Tree -> Solution -> Monitors監(jiān)視器，Edit…打開設(shè)置。

把所有項(xiàng)的絕對殘差都從0.001改為0.00001（我其實(shí)也沒指望它能收斂到這么高的水平，只是改著玩玩罷了），OK確認(rèn)。

設(shè)置樹，Tree -> Solution -> Solution Initialization初始化。

選擇Standard Initialization，Compute from選擇in_inlet，用進(jìn)口值來做初始化，單擊Initialize完成初始化。

計(jì)算前，先保存下設(shè)置文件，File -> Write -> Case…，命名，保存。

設(shè)置樹，Tree -> Solution -> Run Calculation開始計(jì)算。

在Number of Iterations中輸入迭代步3000步，先跑3000步看看效果，如果收斂得不理想，再來幾千步看看。

單擊Calculate開始迭代計(jì)算。

拼機(jī)子的時(shí)刻到了……漫長的等待過程……

7 后處理

先看下殘差曲線，基本上殘差已經(jīng)水平，不會(huì)再往下降了，姑且認(rèn)為計(jì)算結(jié)束了吧，實(shí)際上，除非網(wǎng)格畫的非常精巧，設(shè)置得非常正確，問題非常簡單，才有可能讓殘差收斂到非常小的狀態(tài)……而咱們用這么渣渣的網(wǎng)格和模型，能收斂成這樣已經(jīng)是謝天謝地了。

計(jì)算完成后，可以做切面或者在已有面上查看各種壓力、速度云圖，可以查看流線……

但是咱們最關(guān)心的是泵的性能如何，即，如何獲取其揚(yáng)程、功率、效率等外特性。

7.1 揚(yáng)程計(jì)算

泵的揚(yáng)程 = (出口總壓 - 進(jìn)口總壓) / (介質(zhì)密度 * 重力加速度值) + 進(jìn)出口高度差，數(shù)學(xué)表達(dá)式為：
$$H=\frac{p_{out}?p_{in}}{\rho g}+\Delta z$$
左欄設(shè)置樹，Tree -> Results -> Reports -> Surface Integrals，打開Surface Integrals邊界面的相關(guān)量輸出設(shè)置窗口。

在Surface Integrals窗口中，選擇Report Type報(bào)告類型為Area-Weighted Average即面加權(quán)平均型，Field Variable場變量選擇Pressure壓力和Total Pressure總壓，Surfaces面則選擇in_inlet和vol_outlet，即整體計(jì)算域的進(jìn)口面和出口面，來計(jì)算這兩個(gè)面上的平均總壓值。單擊下方的Compute讓其計(jì)算。

在窗口中有如下輸出信息

Area-Weighted AverageTotal Pressure (pascal) -------------------------------- --------------------vol_outlet 877535.16in_inlet 0---------------- --------------------Net 350022.7

跟設(shè)置的和預(yù)想中的一樣，進(jìn)口的總壓值是$0Pa$，出口的總壓值是$877535.16Pa$，進(jìn)出口的高度差$\Delta z=0.7313m$，這個(gè)跟泵相關(guān)。考慮到水的密度為$\rho =998.2kg/m^3$，重力加速度為$g=9.8kg/m^2$，按揚(yáng)程計(jì)算公式可得
$$\begin{array}{rl}H& =\frac{p_{out}?p_{in}}{\rho g}+\Delta z\\ & =\frac{877535.16?0}{998.2\times 9.8}+0.7313\\ & =90.44m\end{array}$$

設(shè)計(jì)值是84m，算出來居然高達(dá)90.44m，我自己都不相信，顯然，這么粗糙的計(jì)算未必靠譜，好在兩者量級還是比較接近的。

7.2 功率和水力效率計(jì)算

泵的輸入功率，應(yīng)該等于葉片所有壁面的扭矩與轉(zhuǎn)動(dòng)角速度的乘積，即
$$P=T\omega$$
$\omega$根據(jù)轉(zhuǎn)速$2900rpm$很容易就能算出來，扭矩是要讓FLUENT也來report一下的。

左欄設(shè)置樹，Tree -> Results -> Reports -> Forces，打開Forces Reports力的相關(guān)量輸出設(shè)置窗口。

在Force Reports窗口中，Options選項(xiàng)設(shè)置為Moments力矩，Moment Center力矩中心設(shè)為原點(diǎn)$(0,0,0)$，Moment Axis力矩轉(zhuǎn)軸設(shè)為$z$軸，Wall Zones壁面區(qū)域中選擇葉輪的所有壁面，即imp_blade_leading、imp_blade_pressure、imp_blade_suction、imp_blad_tip、imp_hub、imp_shroud，單擊下方Print輸出。

窗口中有如下輸出信息。

整體的轉(zhuǎn)矩為$T=181.33732Nm$，算得其輸入功率為
$$P=T\omega =Tn\pi /30=181.33732\times 2900\times \pi /30=55069.839w=55.070kw$$

泵的有效功率為
$$P_e=\rho gQH=998.2\times 9.8\times (200/3600)\times 90.44=49150.924w=49.151kw$$

那么泵的水力效率為
$$\eta =P_e/P=49.151/55.070=89.25\%$$

當(dāng)然，也可以繼續(xù)設(shè)置不同的出口流速（出口流量$Q$），從而算得泵在不同流量下的揚(yáng)程、功率、效率，繪制出泵的性能曲線來。

8 后話

至此，基于商業(yè)軟件的離心泵的水力設(shè)計(jì)、網(wǎng)格劃分、數(shù)值計(jì)算、性能分析就講完了，這里只是簡單粗暴地走了一遍，并沒有對結(jié)果的可靠性做任何驗(yàn)證工作。

網(wǎng)格劃分部分用的是最簡單最粗糙的非結(jié)構(gòu)均勻網(wǎng)格，數(shù)值模擬部分用的也是最常規(guī)最普通的湍流模型，結(jié)果未必準(zhǔn)確，僅供參考。

網(wǎng)格若想劃分的漂亮，可能要用到邊界層+非結(jié)構(gòu)網(wǎng)格，或者分塊結(jié)構(gòu)化網(wǎng)格，或者嵌套網(wǎng)格等方法，需要耗費(fèi)較大的精力去折騰；而流場若想算得準(zhǔn)確，可能要用到大渦模擬湍流模型，非定常數(shù)值計(jì)算，這已經(jīng)遠(yuǎn)遠(yuǎn)超出了個(gè)人機(jī)的承受范圍……還有，實(shí)話實(shí)說，這么高端的操作我是不會(huì)的，哈哈。

總結(jié)

以上是生活随笔為你收集整理的离心泵CAE_3_FLUENT数值模拟的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。