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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

兩個(gè)數(shù)據(jù)集a,b. 要求比較第一列,保留a中重合部分到repeata, b中重合部分到repeatb 做這個(gè)工作主要是因?yàn)楹芏鄶?shù)據(jù)集要對(duì)比,比如根據(jù)股票的委托單號(hào),來(lái)比較同一委托的逐筆委托和逐筆成交,這時(shí)候就要用到這個(gè)方法 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;再打開(kāi)VS,對(duì)所有projects進(jìn)行build 其中有個(gè)報(bào)錯(cuò)error C3861: 'random': identifier notfound,換成window下的rand()函數(shù),目前對(duì)random的取值范圍和rand的取值范圍是否一致不太清楚 利用mex對(duì)DenseCRF進(jìn)行編譯時(shí),需要將相關(guān)的cpp和h文件進(jìn)行關(guān)聯(lián),...

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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

繼續(xù)閱讀>>

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