讓網(wǎng)格動(dòng)起來(lái)(9)—彈簧光順方法實(shí)例

2017-01-18  by:CAE仿真在線  來(lái)源:互聯(lián)網(wǎng)

上次講到,光順方法計(jì)算網(wǎng)格節(jié)點(diǎn)位移并不是依賴于網(wǎng)格質(zhì)量。因此,僅使用光順方法,只能處理一些小變形的情況 。倘若變形量過(guò)大,則會(huì)造成網(wǎng)格質(zhì)量降低,輕則影響收斂,重則出現(xiàn)負(fù)網(wǎng)格,導(dǎo)致計(jì)算出錯(cuò)。

1、問題描述

image

邊長(zhǎng)120mm的正方形空腔,內(nèi)有一個(gè)邊長(zhǎng)30mm的正方形固體。研究固體以角速度1rad/s旋轉(zhuǎn)時(shí)對(duì)內(nèi)部流場(chǎng)的擾動(dòng)。

1、網(wǎng)格方案

由于四邊形網(wǎng)格在剛體變形時(shí)會(huì)發(fā)生較大的扭曲,不太適合于使用彈簧光順方法,因此本例中使用三角形網(wǎng)格。使用全局尺寸2mm,生成網(wǎng)格如下圖所示。

image

2、指定運(yùn)動(dòng)

可以使用profile文件或DEFINE_CG_MOTION宏來(lái)定義運(yùn)動(dòng)。我們這里利用udf宏來(lái)定義運(yùn)動(dòng)。

UDF宏文件如下:

#include “udf.h”

DEFINE_CG_MOTION(RotationVel , dt , vel , omega , time , dtime)

{

omega[2] = 1;

}

omega[2]表示z方向的角速度 。本例中以恒定的角速度旋轉(zhuǎn)。

將內(nèi)部四條邊定義為rigid body。并把賦予運(yùn)動(dòng)速度。


3、動(dòng)網(wǎng)格參數(shù)

勾選Dynamic Mesh選項(xiàng),激活動(dòng)網(wǎng)格。本例只考慮網(wǎng)格光順,因此只勾選Smoothing。(說(shuō)明:通常動(dòng)網(wǎng)格應(yīng)用中不會(huì)只使用一種網(wǎng)格更新方法,此處只是考慮光順,所以才會(huì)這樣。正常情況下本例應(yīng)當(dāng)使用到網(wǎng)格重構(gòu))。

點(diǎn)擊setting…按鍵進(jìn)行smoothing設(shè)置面板 。如右下圖所示。

SNAGHTML1acfa4e
SNAGHTML1b06733

采用彈簧光順需要設(shè)置的參數(shù)一共有四個(gè):彈簧常數(shù)因子、邊界節(jié)點(diǎn)松弛、收斂精度及迭代次數(shù)。

通常需要設(shè)置的參數(shù)為彈簧常數(shù)因子及邊界節(jié)點(diǎn)松弛,其它兩個(gè)參數(shù)采用默認(rèn)即可。

(1)彈簧常數(shù)因子:該參數(shù)取值范圍[0,1],取值越大表示阻尼越大,也表示邊界運(yùn)動(dòng)對(duì)內(nèi)部網(wǎng)格的影響越小。

(2)邊界節(jié)點(diǎn)松弛:取值范圍[0,1],1表示無(wú)松弛,0表示禁止節(jié)點(diǎn)位移。

這里先將兩參數(shù)取中間值0.5。


3、其它設(shè)置

創(chuàng)建運(yùn)動(dòng)區(qū)域。設(shè)置中間四條邊為rigid body,設(shè)置運(yùn)動(dòng)udf(注意UDF必須為編譯型)。如下圖所示。

SNAGHTML1b74301

4、預(yù)覽網(wǎng)格運(yùn)動(dòng)

(1)觀察0.5s后的網(wǎng)格情況

SNAGHTML1bfdc24

image
image

可以看到,四角位置網(wǎng)格已經(jīng)扭曲得很厲害了。其它參數(shù)有興趣的同學(xué)可以自己嘗試。


5、計(jì)算結(jié)果演示






開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:讓網(wǎng)格動(dòng)起來(lái)(9)—彈簧光順方法實(shí)例 Fluent培訓(xùn) Fluent流體培訓(xùn) Fluent軟件培訓(xùn) fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學(xué)反應(yīng) fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

編輯
在線報(bào)名:
  • 客服在線請(qǐng)直接聯(lián)系我們的客服,您也可以通過(guò)下面的方式進(jìn)行在線報(bào)名,我們會(huì)及時(shí)給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國(guó)服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號(hào)金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項(xiàng)目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com