式（5）是基于質(zhì)量平衡的方程，利用式（5）可以求解流體流過(guò)各向異性多孔介質(zhì)的的問(wèn)題，壓力是求解的目的參數(shù)。解方程的方法很多，本研究中采用有限元/控制體積方法（FEM/CV）。
    2  充模過(guò)程模擬分析優(yōu)化
    2．1 風(fēng)機(jī)葉片RTM工藝模擬
    根據(jù)上述數(shù)學(xué)模型自行開(kāi)發(fā)了數(shù)值計(jì)算軟件，此軟件與PATRAN相結(jié)合構(gòu)成了基于有限元/控制體積方法（FEM/CV）的RTM工藝模擬仿真系統(tǒng)。利用此模擬仿真系統(tǒng)對(duì)風(fēng)機(jī)葉片的充模過(guò)程進(jìn)行了分析和工藝優(yōu)化。
    風(fēng)機(jī)葉片的木制模型如圖1所示，利用有限元分析軟件PATRAN（對(duì)風(fēng)機(jī)葉片（尺寸：1300mm×173mm×43mm）進(jìn)行建模（三維）。采用四結(jié)點(diǎn)四邊形單元對(duì)整體零件進(jìn)行有限元部分，得到的有限元模型如圖2所示，模型中包含結(jié)點(diǎn)3225個(gè)，單元3215個(gè)。

    采取線注射方式，注口位置如圖3所示，將復(fù)合材料盒體的有限元模型數(shù)據(jù)輸入三維計(jì)算軟件。在計(jì)算中設(shè)定滲透率為1.0×10-9m2,纖維體積含量35%:

      樹(shù)脂粘度0.3Pa.s，注射壓力為0.15MPa，在PATRAN中處理計(jì)算所得數(shù)據(jù)，并顯示不同時(shí)刻的流動(dòng)前峰位置和模腔內(nèi)的壓力分布，如圖4所示。分析流動(dòng)前峰位置和壓力分布圖時(shí)，以圖左側(cè)的色柱為參考，色柱左側(cè)標(biāo)有單位（時(shí)間和壓力），數(shù)值由下至上增大，色柱上色塊的顏色由下至上從白色開(kāi)始，黑度依次變化。流動(dòng)前峰圖書館時(shí)（圖4（a）），充模過(guò)程中最先填充的地方即注口處顯白色，然后構(gòu)件上條紡的顏色沿樹(shù)脂流動(dòng)方向以色柱上色塊的黑度由下至上變化，從而可以考察充模過(guò)程中樹(shù)脂流動(dòng)前鋒的位置。壓力分布圖中（圖4（b）），最先填充的地主壓力最大，因而此處的顏色是其色柱上頂端色塊的顏色，然后條紋的顏色沿樹(shù)脂流動(dòng)動(dòng)方向以色柱上色塊的黑度由上至下變化，從而考察模腔內(nèi)的壓力分布。最終為確定注射口和溢料口的位置提供依據(jù)。實(shí)際工藝過(guò)程：注射壓力0.15MPa，纖維體積含量35%，樹(shù)脂粘度0.3Pa.s。實(shí)際充模的時(shí)間為在180～300s，與模擬充模時(shí)間176.1s比較接近。可見(jiàn)模擬充模時(shí)間對(duì)實(shí)際工藝有一定的指導(dǎo)意義，而更重要的是通過(guò)分析模擬結(jié)果可以優(yōu)化工藝。因整個(gè)充模過(guò)程在封閉不透明的模具中進(jìn)行，可得的實(shí)驗(yàn)數(shù)據(jù)只有充模時(shí)間。