第一次見到老藍,小白都不太敢相信,對面那不修邊幅的糟老頭子會是自己要找的導(dǎo)師。嘴里叼著煙,牙都掉了好幾顆,穿著還算整齊,這是小白對老藍的第一印象,這印象并不太好,尤其是在小白發(fā)誓認真度過研究生三年時光之后。 那是一個秋風(fēng)蕭瑟的下午,沒有課的小白午休后還沒起床,忽然收到老藍的電話。“小白啊,下午沒課來我辦公...

繼續(xù)閱讀>>

雖然已經(jīng)是9月份了,但是西南地區(qū)的濕熱天氣依然將小白熱成了狗。獨自一個人來到這號稱天府之國的西南城市,小白卻一點都沒感受到這天府之國的友好。毒辣的陽光像利劍一樣照射到身上,讓人感覺如在蒸籠中一樣。從炎熱的華中火爐來到祖國大西南,卻沒感受到哪怕一絲的清涼。 今天是開學(xué)的日子,學(xué)校將研究生和本科生新生入學(xué)...

繼續(xù)閱讀>>

小白的CFD之旅 寫在前面 CFD是計算流體力學(xué)的英文簡稱,是計算機輔助工程(CAE)的主要分支,目前廣泛應(yīng)用與科學(xué)研究、工程設(shè)計中。這是一門綜合了數(shù)學(xué)、計算機及流體力學(xué)的綜合學(xué)科,涉及到眾多的專業(yè)理論,如果缺少相應(yīng)的專業(yè)基礎(chǔ),要想將CFD應(yīng)用于工程中則需要花費較多的時間。 CFD工程應(yīng)用主要涉及到以下幾個方面: 力學(xué)...

繼續(xù)閱讀>>

用過STAR-CD或STARCCM+的人一定對其中的多面體網(wǎng)格(Polyhedra)功能印象深刻。多面體網(wǎng)格相對于非六體網(wǎng)格的其它各種類型的網(wǎng)格來說,最大的優(yōu)勢在于其能大大減小網(wǎng)格的數(shù)量。至于其它方面的優(yōu)勢,可以查閱相關(guān)的資料。其實在FLUENT中也是可以應(yīng)用多面體網(wǎng)格的。FLUENT中可以通過將四面體、金字塔等非六面體網(wǎng)格轉(zhuǎn)化為...

繼續(xù)閱讀>>

流固耦合(Fluid-solidinteraction,FSI)計算,通常用于考慮流體與固體間存在強烈的相互作用時,對流體流場與固體應(yīng)力應(yīng)變的考察。FSI計算按數(shù)據(jù)傳遞方式可分兩類:單向耦合與雙向耦合。所謂單向耦合,主要是指數(shù)據(jù)只從流體計算傳遞壓力到固體,或者只從固體計算傳遞網(wǎng)格節(jié)點位移到流體。雙向耦合則在每一時刻都同時向?qū)Ψ桨l(fā)...

繼續(xù)閱讀>>

動態(tài)層方法(dynamiclayering)也許是fluent動網(wǎng)格中最容易理解與設(shè)置的方法了。動態(tài)層方法主要用于棱柱網(wǎng)格區(qū)域(如六面體網(wǎng)格、三棱柱網(wǎng)格等)或四邊形網(wǎng)格(2D幾何)中,用戶可以使用動態(tài)層方法增加或去除與運動邊界相鄰的單元網(wǎng)格層數(shù)。在ANSYSFLUENT中,允許用戶為每一個運動邊界指定一個理想高度作為網(wǎng)格層增加或去除...

繼續(xù)閱讀>>

FLUENT中的網(wǎng)格光順(smoothing)主要采用兩種方法:彈簧光順與擴散光順。 其中彈簧光順包含有:彈簧光順、邊界層光順、拉普拉斯光順。 擴散光順分為:基于壁面距離的擴散光順與基于單元體積的擴散光順。 一、彈簧/拉普拉斯/邊界層光順 1、彈簧光順需要設(shè)置的參數(shù) 主要有四個參數(shù)需要設(shè)置:SpringConstant Fa...

繼續(xù)閱讀>>

Fluent的動網(wǎng)格設(shè)置,主要是通過dynamicmesh面板進行設(shè)置。動網(wǎng)格面板如下圖所示。 主要包括兩部分內(nèi)容:(1)MeshMethods (2) Dynamic MeshZones。另外還包括網(wǎng)格域運動預(yù)覽及網(wǎng)格運動預(yù)覽。整個動網(wǎng)格面板通過勾選Dynamic Mesh選項而激活。 1、MeshMethods(網(wǎng)格更新方法) 該部分主要用于網(wǎng)格或...

繼續(xù)閱讀>>

本次實例采用的場景來自于流體中高速飛行的物體。如子彈、火箭、導(dǎo)彈等。這里只是為了說明profile在動網(wǎng)格運動定義中的應(yīng)用,因此為了計算方便不考慮高速問題。問題描述如下圖所示: 如圖所示,1為運動剛體,2為計算域。由于不考慮也沒辦法考慮剛體的變形,因此在構(gòu)建面域的時候,將1中的部分通過布爾運算去除。計算...

繼續(xù)閱讀>>

在固體有限元計算中,網(wǎng)格運動實非什么稀奇事兒。而且在絕多數(shù)固體計算的基本物理量是網(wǎng)格的節(jié)點位移,所以,固體計算中,網(wǎng)格節(jié)點運動是對的,沒有運動反而不正常了。也可以這么說:正因為計算域內(nèi)部節(jié)點間的相對運動,才導(dǎo)致了內(nèi)應(yīng)力的產(chǎn)生。 流體計算與固體完全不同。其根源在于它們使用的網(wǎng)格類型不同。當(dāng)前固體有限...

繼續(xù)閱讀>>

學(xué)習(xí)一款軟件,其最基本的要求是了解操作界面。必須清楚的知道軟件操作界面上每一個按鈕、每一個文本框所代表的含義,才可能得心應(yīng)手的解決我們自己的問題。 FLUNET是一款CFD求解器。其讀入用戶輸入的模型、邊界及求解控制信息,求解內(nèi)嵌的控制方程,最后輸出計算結(jié)果數(shù)據(jù)。它對于用戶來說就是一個黑盒子,除了軟件開發(fā)...

繼續(xù)閱讀>>

對于三角形或四面體風(fēng)格區(qū)域,通常都會使用到彈簧光順方法。然而當(dāng)邊界位移遠大于局部網(wǎng)格尺寸時,網(wǎng)格質(zhì)量會下降甚至?xí)霈F(xiàn)退化的網(wǎng)格單元。這可能會使網(wǎng)格實效(出現(xiàn)負體積網(wǎng)格),并因此導(dǎo)致計算收斂問題。 為什么彈簧光順方法會導(dǎo)致此類問題?主要原因在于:彈簧光順方法并不將網(wǎng)格質(zhì)量作為網(wǎng)格更新的判斷因素,而只...

繼續(xù)閱讀>>

上次講到,光順方法計算網(wǎng)格節(jié)點位移并不是依賴于網(wǎng)格質(zhì)量。因此,僅使用光順方法,只能處理一些小變形的情況。倘若變形量過大,則會造成網(wǎng)格質(zhì)量降低,輕則影響收斂,重則出現(xiàn)負網(wǎng)格,導(dǎo)致計算出錯。 1、問題描述 邊長120mm的正方形空腔,內(nèi)有一個邊長30mm的正方形固體。研究固體以角速度1rad/s旋轉(zhuǎn)時對內(nèi)部流...

繼續(xù)閱讀>>

與光順方法類似,動態(tài)層方法同樣不是以網(wǎng)格質(zhì)量作為網(wǎng)格更新判據(jù)的。動態(tài)層方法原理很簡單,且非常適合于結(jié)構(gòu)網(wǎng)格在某一方向上具有主導(dǎo)運動的情況下。本例即是這么一種情況。 1、問題描述 如下圖所示,計算域空間尺寸160mmx50mm,一個長方形固體塊在大的空間沿x軸勻速運動,速度0.05m/s,運動時間2s,可計算出運動位移0...

繼續(xù)閱讀>>

模型選擇 對于某些問題,某個輻射模型可能比其它模型更適用。在確定使用何種輻射模型時,需要考慮的因素如下:?? 光學(xué)深度:光學(xué)深度αL 是確定選擇輻射模型較好的指標(biāo)。其中,L 為計算域大致的長度標(biāo)尺。例如,對于燃燒室內(nèi)的流動, L為燃燒室的直徑。如果光學(xué)深度》1,那么,最好的選擇是使用P-1 或Rosseland 輻射模型。P-1模...

繼續(xù)閱讀>>

這里要說的PROFILE文件,只針對利用于網(wǎng)格運動定義之用的瞬態(tài)profile文件。其它類型的profile文件,可參閱FLUENT用戶文檔,里頭有詳細的描述。 在ANSYSFLUENT中,有兩種方式用于指定瞬態(tài)網(wǎng)格區(qū)域及邊界條件: (1)與標(biāo)準(zhǔn)profile格式相同的瞬態(tài)profile文件 (2)表格格式的瞬態(tài)profile文件。 對于這兩種方式,網(wǎng)格...

繼續(xù)閱讀>>

雖然瞬態(tài)PROFILE文件可以在一定程度上定義網(wǎng)格運動,然而其存在著一些缺陷。最主要的一些缺陷存在于以下一些方面: (1)PROFILE無法精確的定義連續(xù)的運動。其使用離散的點值進行插值。如果想獲得較為精確的運動定義,勢必要定義很多點。 (2)一些情況下無法使用profile。比如穩(wěn)態(tài)動網(wǎng)格。 在FLUENT中定義網(wǎng)格運動,更多的是采...

繼續(xù)閱讀>>

FLUENTUDF中關(guān)于動網(wǎng)格的DEFINE宏均要求為編譯型。FLUENT的UDF運行環(huán)境有兩種:解釋型與編譯型。雖然兩種方式都是采用編譯,但是編譯方式及編譯后獲得的結(jié)果是不同的。編譯型UDF需要通過Complied進行編譯,編譯后得到相應(yīng)的dll文件。而解釋型UDF則是通過Interprete進行加載,不會生成dll文件。據(jù)我的直覺,編譯型UDF可能...

繼續(xù)閱讀>>

FLUENT是一款對使用者要求較高的CFD軟件,不管是從前處理的網(wǎng)格生成、計算模型選取及參數(shù)設(shè)置、求解參數(shù)設(shè)置,還是后處理對計算結(jié)果的科學(xué)解讀,無不要求使用者對其所要解決的物理現(xiàn)象有深刻的認識。 從本質(zhì)上來講,FLUENT僅僅只是一款CFD軟件,其無法對我們的現(xiàn)實世界做出正確的判斷,關(guān)于物理現(xiàn)象的一切信息均由使用者輸...

繼續(xù)閱讀>>

uigetfile命令的應(yīng)用 此函數(shù)的用法為 [FileName,PathName,FilterIndex] =uigetfile(FilterSpec,DialogTitle,DefaultName) FileName:返回的文件名 PathName:返回的文件的路徑名 FilterIndex:選擇的文件類型 FilterSpec:文件類型設(shè)置 DialogTitle:打開對話框的標(biāo)題 DefaultName:默認指向的文件名...

繼續(xù)閱讀>>

兩個數(shù)據(jù)集a,b. 要求比較第一列,保留a中重合部分到repeata, b中重合部分到repeatb 做這個工作主要是因為很多數(shù)據(jù)集要對比,比如根據(jù)股票的委托單號,來比較同一委托的逐筆委托和逐筆成交,這時候就要用到這個方法 a=[1,11,111;2,22,222;3,33,333]; b=[2,121,222;3,323,323;4,434,443;5,565,656]; 法一: index=ismem...

繼續(xù)閱讀>>

利用CMake直接加載相應(yīng)文件夾,利用已有的CMake文件generate生成solution;再打開VS,對所有projects進行build 其中有個報錯error C3861: 'random': identifier notfound,換成window下的rand()函數(shù),目前對random的取值范圍和rand的取值范圍是否一致不太清楚 利用mex對DenseCRF進行編譯時,需要將相關(guān)的cpp和h文件進行關(guān)聯(lián),...

繼續(xù)閱讀>>

dlmwrite 功能簡介 將矩陣寫入ASCII分隔的文件。 使用方法 1.dlmwrite(filename, M) 使用默認的分隔符(')將矩陣M寫入ASCII格式的文件中。在目標(biāo)文件filname中,數(shù)據(jù)是從第一行的第一列開始寫的。輸入的filename是使用單引號括起來的字符串。matlab2.dlmwrite(filename, M, 'D')matabc ...

繼續(xù)閱讀>>

1、查看Matlab是否配置完全 在Command Window窗口中輸入mbuild–setup查看編譯器,如果沒有,將編輯器設(shè)置為vs 2、安裝MCRinstaller.exe (這步很重要,否則無法編譯生成dll) 位置在安裝目錄下的toolbox\compiler\deploy\win64文件夾中。 3、注冊dll 對mwcomutil.dll進行注冊。這個Matlab的dll文件位于安裝目錄下的\bin\w...

繼續(xù)閱讀>>

函數(shù)rand用來生成(0,1)之間均勻分布的偽隨機數(shù),調(diào)用格式: ? Y=rand:生成一個隨機數(shù) ? Y=rand(n):生成n*n的隨機矩陣 ? Y=rand(m,n):生成m*n的隨機矩陣 ? Y=rand(size(A)):生成與矩陣A大小相同的隨機矩陣 ? Y=rand(m,n,p,…):生成m*n*p*…的隨機數(shù)多維數(shù)組 對多維的具體理解 Y=rand(m,...

繼續(xù)閱讀>>

可以用MATLAB和1stOpt進行曲線擬合,我總結(jié)了5種方法: MATLAB 1、polyfit/polyconf命令 2、regress:多元線性回歸 3、畫圖窗口的基本擬合工具 fig格式的圖片可以進行一些簡單的曲線擬合,位置在fig窗口菜單欄-工具-基本擬合 4、cftool(推薦):多種類型的線性、非線性曲線擬合gui工具。 可...

繼續(xù)閱讀>>

耐磨鋼板是通過高溫加熱和壓力變形實現(xiàn)基層材(碳鋼或低合金鋼)與覆層材(耐磨鋼板)冶金結(jié)合,使其成為一種具有特殊性能的整體疊層材料?；谀臀g性能考慮,寬厚耐磨鋼板在工程機械、礦山、軍事等高端產(chǎn)業(yè)領(lǐng)域有著廣闊的應(yīng)用前景。相對于爆炸復(fù)合法和爆炸-軋制復(fù)合法而言,采用層間真空軋制生產(chǎn)耐磨鋼板的方法更能適應(yīng)大規(guī)...

繼續(xù)閱讀>>

近來搞了個人臉檢測程序。用到matlab打包技術(shù),網(wǎng)上百度了一下,挺有幫助的。 具體過程是: 1.在matlab的命令窗口,輸入mbuild -setup.之后就有提示,選擇C++編譯器吧。點擊了一下下面的提示。mex -setup C++什么的。 2.據(jù)說也可以用mcc -m function1.mfunction2.m...來編譯,但是我是失敗了。所以我在命令行輸入deploytool ...

繼續(xù)閱讀>>

m=size(A,1)返回矩陣的行數(shù); n=size(A,2)返回矩陣的列數(shù); sum(A)矩陣的每一列求和;

繼續(xù)閱讀>>

reply = 'Y'; while reply == 'Y' score=input('Please input your score: '); if isempty(score) score = 60; end if (score < 60) && (score >= 0) disp('E'); elseif (score >= 60) && (score < 70) disp('D'); elseif (score >= 70) && (score < 80) disp('C'); elseif (...

繼續(xù)閱讀>>