版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

bananaliu

木蟲 (正式寫手)

應(yīng)助: 28 (小學(xué)生)
金幣: 3046.4
紅花: 2
帖子: 524
在線: 1411.7小時(shí)
蟲號(hào): 1179664
注冊(cè): 2010-12-29
性別: GG
專業(yè): 巖土與基礎(chǔ)工程

[求助] Fluent與Abaqus聯(lián)合仿真計(jì)算的問題

小弟最近在編寫Fluent和Abaqus的接口程序，打算實(shí)現(xiàn)串行算法下的（Gauss-Seidel algorithm）的流固耦合（FSI）。主要是通過UDF和Scheme配套Abaqus的Co-simulation Engine (CSE) 實(shí)現(xiàn)。目前工作進(jìn)行到一半，碰到一些技術(shù)難題無法解決。
首先，我沒有Mpcci的權(quán)限，沒有權(quán)限無法發(fā)表論文，所以放棄。
碰到的問題就是，
1. 我不知道怎么實(shí)現(xiàn)Fluent計(jì)算過程中暫停，然后有條件（等待Abaqus CSE 傳來固體計(jì)算獲得的流域邊界，應(yīng)該是用define_on_demand & define_grid_motion實(shí)現(xiàn)我覺得）重啟計(jì)算。比如流體總共要計(jì)算600s，每計(jì)算100s暫停，輸出變量數(shù)據(jù)，等待網(wǎng)格更新信息，然后計(jì)算下100s。我想了一個(gè)權(quán)宜的辦法，把600s分成6個(gè)任務(wù)，每個(gè)任務(wù)計(jì)算100s，采用如下命令循環(huán)：

CODE:

(define do-simulation        (lambda (a b)

                (if (<= a b)

                        (begin 

                                ;;(rpsetvar 'Current_Time (* (/ Cosimu_time Num_exchange) (rpgetvar 'CurrNum_exchange)))

                                (display "\n1")

                                (%udf-on-demand "init::libudf")

                                ;;(ti-menu-load-string "/define/user-defined/execute-on-demand \"init::libudf\"")

                                (display "\n2")

                                (ti-menu-load-string (format #f "/solve/dual-time-iterate ~d ~d" Num_timestep Max_iter_per_timestep))

                          (set! a (+ a 1))

                          (rpsetvar 'CurrNum_exchange a)

                          (do-simulation a b)

                  )

                )

        )

)

就是兩個(gè)任務(wù)之間不再進(jìn)行初始化，可不知道這樣計(jì)算和一次性計(jì)算600s效果是否相同。而且這個(gè)重啟命令是在GUI情況下實(shí)現(xiàn)，類似于界面上點(diǎn)擊calculate后再cancel一樣，不知道能不能再dos命令狀態(tài)下實(shí)現(xiàn)（即命令 fluent xxx.cas ....）。
2. 如果在外部程序（比如A.cpp）調(diào)用批處理文件（.scm），scm文件中又調(diào)用UDF，如何實(shí)現(xiàn)A文件和UDF文件共享變量？畢竟UDF是加載到Fluent運(yùn)行環(huán)境中，而A是在外部直接運(yùn)行。
有經(jīng)驗(yàn)的同學(xué)能不能給點(diǎn)建議。多謝！

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

流固熱耦合

» 猜你喜歡

求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(qǐng)(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩301分求調(diào)劑已經(jīng)有7人回復(fù)

科研無止境，項(xiàng)目累死人。

1樓 2015-01-23 12:01:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

有關(guān)fluent并行計(jì)算的問題已經(jīng)有16人回復(fù)
做一個(gè)流體和固體的傳熱流固耦合CFX+ANSYS還是ABAQUS+FLUENT？已經(jīng)有3人回復(fù)
【求助】在用fluent計(jì)算內(nèi)波時(shí)，operating pressure以及operating density的設(shè)定問題已經(jīng)有6人回復(fù)
理論計(jì)算結(jié)果與仿真數(shù)值之間的關(guān)系已經(jīng)有8人回復(fù)
fluent非穩(wěn)態(tài)計(jì)算自動(dòng)保存的問題，你遇到?jīng)]？已經(jīng)有6人回復(fù)
abaqus膠接仿真中的相關(guān)屬性定義，求助各位朋友。。。已經(jīng)有5人回復(fù)
fluent 歐拉兩相流模擬穩(wěn)態(tài)計(jì)算松弛因子的問題已經(jīng)有4人回復(fù)
fluent計(jì)算多相流問題時(shí)，UDF能否提取單元內(nèi)分相速度，以及單元平均速度已經(jīng)有5人回復(fù)
【求助】請(qǐng)教一下ABAQUS的質(zhì)量擴(kuò)散模塊計(jì)算材料的水分吸濕問題已經(jīng)有6人回復(fù)
Fluent計(jì)算過程中出現(xiàn)的一個(gè)小問題（大家?guī)兔忉屢幌拢?/a> 已經(jīng)有11人回復(fù)

Fluent計(jì)算時(shí)的迭代問題已經(jīng)有7人回復(fù)
hypermesh的網(wǎng)格導(dǎo)入abaqus計(jì)算時(shí)提示存在2 unconnected regions 已經(jīng)有12人回復(fù)
Fluent計(jì)算時(shí)的迭代問題已經(jīng)有14人回復(fù)
求助}abaqus循環(huán)計(jì)算問題已經(jīng)有4人回復(fù)
在做fluent例子的時(shí)候遇到了問題，發(fā)現(xiàn)無法進(jìn)行迭代計(jì)算，求教大家已經(jīng)有9人回復(fù)
fluent的軸向?qū)嵊?jì)算問題已經(jīng)有4人回復(fù)
【求助】fluent的傳熱問題，計(jì)算結(jié)果有沒有問題，請(qǐng)大家?guī)兔纯础！疽呀鉀Q】已經(jīng)有10人回復(fù)
【求助】fluent中壁面怒謝爾特?cái)?shù)的計(jì)算問題已經(jīng)有15人回復(fù)
【求助】關(guān)于fluent多機(jī)并行計(jì)算連接問題已經(jīng)有6人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 bananaliu 的主題更新