采用一個(gè)簡(jiǎn)單的雙向流固耦合算例對(duì)比lsdyna和workbench的計(jì)算結(jié)果，結(jié)果顯示dyna的計(jì)算結(jié)果不太理想，如有疑議，dyna計(jì)算的k文件在文中，請(qǐng)dyna高手賜教了。

模型尺寸如下：

厚度方向0.4m。

平板材料E= 2.5e06[Pa];密度2550 kg m^-3 ；泊松比0.35

流體是空氣，密度1kg m^-3 ；粘度0.2 kg/m-s

左側(cè)是速度入口邊界，沿x軸方向速度為4m/s；右側(cè)為壓力出口，壓強(qiáng)為0；前后及上側(cè)為自由滑移壁面；底面及平板為非自由滑移壁面。

幾何模型：

FSI_GEOM.zip

ansys workbench計(jì)算模型：

fsi_workbench.7z

lsdyna計(jì)算模型：

fsi_dyna.zip 算到5秒

fsi_dyna_to5_5.rar 計(jì)算到5.5秒

fsi_dyna_to6.zip 計(jì)算到6秒

采用雙精度求解器求解

由于lsdyna中的ICFD只支持四面體網(wǎng)格，所以網(wǎng)格與workbench中的六面體網(wǎng)格不同：

進(jìn)行流固耦合求解，經(jīng)常會(huì)遇到因?yàn)閱卧兌鴮?dǎo)致的計(jì)算終止

壓強(qiáng)與速度的時(shí)間歷程動(dòng)畫(huà)與workbench的差異較大，特別是dyna中平板的擺動(dòng)規(guī)律比較奇怪，感覺(jué)平板在那自激振動(dòng)了一樣。

在5秒時(shí)刻，對(duì)稱(chēng)面上的速度分布，ansys中的計(jì)算結(jié)果：

dyna第2個(gè)模型文件計(jì)算的結(jié)果：

dyna第3個(gè)模型文件計(jì)算的結(jié)果：

可以看出不同流體網(wǎng)格下dyna計(jì)算的結(jié)果有差異，速度值比ansys workbench中的更大，速度分布也非常紊亂。