ansys子模型介紹與應(yīng)用實例

1、精心整理第五章子模型何為子模型？子模型是得到模型部分區(qū)域中更加精確解的有限單元技術(shù)。在有限元分析中往往出現(xiàn)這種情 況，即 對于用戶關(guān)心的區(qū)域，如應(yīng)力集中區(qū)域，網(wǎng)格太疏不能得到滿意的結(jié)果，而對于這些區(qū)域 之外的部分， 網(wǎng)格密度已經(jīng)足夠了。見圖5-1 o圖5-1輪轂和輪輻的子模型a）粗糙模型，b）疊加的子模型要得到這些區(qū)域的較精確的解，可以采取兩種辦法：（a）用較細(xì)的網(wǎng)格重新劃分并分析整個模 型, 或（b）只在關(guān)心的區(qū)域細(xì)化網(wǎng)格并對其分析。顯而易見，方法a太耗費機(jī)時，方法b即為子模型技術(shù)。子模型方法又稱為切割邊界位移法或特定邊界位移法。切割邊界就是子模型從整個較粗糙的模型分 割開的邊界。整體模型切

2、割邊界的計算位移值即為子模型的邊界條件。子模型基于圣維南原理，即如果實際分布載荷被等效載荷代替以后，應(yīng)力和應(yīng)變只在載荷施加的位置附近有改變。這說明只有在載荷集中位置才有應(yīng)力集中效應(yīng)，如果子模型的位置遠(yuǎn)離應(yīng)力集中位置, 則子模型內(nèi)就可以得到較精確的結(jié)果。ANSYS?序并不限制子模型分析必須為結(jié)構(gòu)（應(yīng)力）分析。子模型也可以有效地應(yīng)用于其他分析中。 如在電磁分析中，可以用子模型計算感興趣區(qū)域的電磁力。除了能求得模型某部分的精確解以外，子模型技術(shù)還有幾個優(yōu)點： 它減少甚至取消了有限元實體模型中所需的復(fù)雜的傳遞區(qū)域。它使得用戶可以在感興趣的區(qū)域就不同的設(shè)計（如不同的圓角半徑）進(jìn)行分析。它幫助用戶證明網(wǎng)格

3、劃分是否足夠細(xì)。使用子模型的一些限制如下：只對體單元和殼單元有效。子模型的原理要求切割邊界應(yīng)遠(yuǎn)離應(yīng)力集中區(qū)域。用戶必須驗證是否滿足這個要求。如何作子模型分析子模型分析的過程包括以下步驟：1 .生成并分析較粗糙的模型。2 .生成子模型。3 .提供切割邊界插值。4 .分析子模型。5 .驗證切割邊界和應(yīng)力集中區(qū)域的距離應(yīng)足夠遠(yuǎn)。第一步：生成并分析較粗糙的模型第一個步驟是對整體建模并分析。（注-為了方便區(qū)分這個原始模型，我們將其稱為粗糙模型。這并不表示模型的網(wǎng)格劃分必須是粗糙的，而是說模型的網(wǎng)格劃分相對子模型的網(wǎng)格是較粗 糙的。）分析類型可以是靜態(tài)或瞬態(tài)的，其操作與各分析的步驟相同。下面列出了其他的一

4、些要記住的方面。文件名一一粗糙模型和子模型應(yīng)該使用不同的文件名。這樣就可以保證文件不被覆蓋。而且 在切割 邊界插值時可以方便地指出粗糙模型的文件。用下列方法指定文件名：Comma nd:/FILNAMEGUI:UtilityMe nu>File>Cha ngejob name單元類型一一子模型技術(shù)只能使用塊單元和殼單元。分析模型中可以有其他單元類型（如梁單元作 為加強(qiáng)筋），但切割邊界只能經(jīng)過塊和殼單元。一種特殊的子模型技術(shù)，稱為殼到體子模型技術(shù)，允許用戶用殼單元建立粗糙模型而用三維塊單元 建立子模型。本技術(shù)在后面還要討論。精心整理建模一一在很多情況下，粗糙模型不需要包含局部的細(xì)節(jié)如

5、圓角等，見下圖。但是，有限元網(wǎng)格必 須細(xì)化到足以得到較合理的位移解。這一點很重要，因為子模型的結(jié)果是根據(jù)切割邊界的位移解插值得 到的。圖5-2粗糙模型可以不包括一些細(xì)節(jié)部分文件結(jié)果文件（Jobname. RST, Jobname. RM等）和數(shù)據(jù)庫文件（Jobname. DB包含幾何模型）在粗糙模型分析中是需要的。在生成子模型前應(yīng)存儲數(shù)據(jù)庫文件。用下列方法存儲數(shù)據(jù)庫： Comma nd:SAVEGUI:UtilityMenu>File>SaveasUtilityMenu>File>SaveasJobname, db第二步：生成子模型子模型是完全依靠粗糙模型的。因此在初始

6、分析后的第一步就是在初始狀態(tài)清除數(shù)據(jù)庫（另一種方法是退出并重新進(jìn)入ANSY） S o用下列方法清除數(shù)據(jù)庫：Command:/CLEARGUI:Ut i1i tyMenu>File>Clear&StartNew同時: 應(yīng)記住用另外的文件名以防止粗糙模型文件被覆蓋。用下列方法指定文件名：Command:/FILNAMEGUI:UtilityMenu>File>ChangeJobname然后進(jìn)入PREP曲建立子模型。應(yīng)該記住下列幾點：使用與粗糙模型中同樣的單元類型。同時應(yīng)指定相同的單元實參（如殼厚）和材料特性。（另一種子模型技術(shù)一一殼到體技術(shù)一一允許從粗糙模型的殼單元

7、轉(zhuǎn)換為體單元，見后。）子模型的位置（相對全局坐標(biāo)原點）應(yīng)與粗糙模型的相應(yīng)部分相同，見圖5-3o圖5-3疊加在粗糙模型上的子模型指定合適的結(jié)點旋轉(zhuǎn)位移。切割邊界結(jié)點的旋轉(zhuǎn)角在插值步驟一寫入結(jié)點文件時不應(yīng)改變（見第三步：生成切割邊界插值）。用下列方法指定結(jié)點旋轉(zhuǎn)：Command:NROTATGUI:MainMenu>Preprocessor>Create>Nodes>-RotateNodeCS-ToActiveCSMainMenu>Preprocessor>Move/Modify>-RotateNodeCS-ToACtiveCS注意結(jié)點旋轉(zhuǎn)角會因為施加結(jié)點

8、約束，傳遞線上約束或面上約束等操作而改變，同樣也會為更加明顯的操作如NROTAT和NM0DIF等改變。粗糙模型中結(jié)點旋轉(zhuǎn)角的出現(xiàn)或缺省并不影響子模型。子模型的載荷和邊界條件將在后面兩步中施加。第三步：生成切割邊界插值本步是子模型的關(guān)鍵步驟。用戶定義切割邊界的結(jié)點，ANSYSi序用粗糙模型結(jié)果插值方法計算這些點上的自由度數(shù)值（位移等）。對于子模型切割邊界上的所有結(jié)點，程序用粗糙模型網(wǎng)格中相 應(yīng)的單元確定自由度數(shù)值，然后這些數(shù)值用單元形狀功能插值到切割邊界上。在切割邊界插值中有下面幾步操作：1.指定子模型切割邊界的結(jié)點并將其寫入一個文件（缺省為Jobname. NODJE中?？梢栽赑REP7中選擇

9、切割邊界的結(jié)點，用下列命令將其寫入文件：Command:NWRITEGUI:MainMenu>Preprocessor>Create>Nodes>WriteNodeFile 下面是一個NWRIT命令的例子：NSEL,. !選擇切割邊界上的結(jié)點NWRITE 將其寫入 Jobname. NODE 圖5-4子模型切割邊界 在這里討論一下溫度插值的問題。在包含特性隨溫度變化的材料的分析中，或熱一應(yīng)力耦合分析 中，粗糙模型和子模型中的溫度分布是相同的。在這種情況下，必須將粗糙模型的溫度插值到子模型中 的所有結(jié)點上。要完成這步操作，要選擇子模型中所有結(jié)點并寫入另外一個文件中，精心整

10、理使用WRITE, Filename, Ext o記住必須另外指定一個文件名，否則切割邊界結(jié)點文件將被覆蓋！第7步中說明了關(guān)于溫度插值的命令。2 .重新選擇所有結(jié)點并將數(shù)據(jù)庫存入Jobname.DB中，然后退出PREP7必須將數(shù)據(jù)庫寫入文 件，因為在后面子模型分析中要使用到。用下列命令重新選擇所有結(jié)點：Command:ALLSELGUI:UtilityMenu>Select>Everything用下列命令存儲數(shù)據(jù)庫：Command:SAVEGUI:Ut i1i tyMenu>Fi1e>SaveasJobname, db3 .要進(jìn)行切割邊界插值(和溫度插值)，數(shù)據(jù)庫中必須

11、包含粗糙模型的幾何特征。因此要用 下列一種方法讀入粗糙模型數(shù)據(jù)庫：Command:RESUMEGUI:UtilityMenu>File>Resumefrom女口，粗糙模型文件名為COARSE就輸入命令RESUME, COARSE<, DB4 .進(jìn)入 POST1 即通用處理器(/POST1 或 MainMenu>GeneralPostproc) o 插值只有在 POST1 中進(jìn)行。5 .指向粗糙模型結(jié)果文件(FILE 或 MainMenu>GeneralPostproc>Data&FileOpts )。6 . 讀入結(jié)果文件中相應(yīng)的數(shù)據(jù)(SET 或 Ma

12、inMenu>GeneralPostproc>-ReadResults-option )。7.開始切割邊界插值。用下列方法完成本步操作：Command:CBDOFGUI:MainMenu>GeneralPostproc>Submodeling>InterpolateDOF缺省狀態(tài)下，CBDO命令假定切割邊界結(jié)點在文件Jobname. NOD中。ANSYS §序?qū)⒂嬎闱懈钸吔绲腄OF數(shù)值并用D命令的形式寫入文件Jobname. CBDOPo用下列方法作溫度插值，但要保證文件包含所有子模型結(jié)點：Command:BFINTGUI:MainMenu>Gen

13、eralPostproc>Submodeling>InterpBodyForc溫度插值以BF命令的格式寫入文件Jobname. BFIN中。注一一如果數(shù)據(jù)包括實部和虛部的話，步驟 6和7就要作兩遍。先用SET命令讀入實部的數(shù)據(jù)并作插值CBDOF和/或BFINT,然后用SET命令將域設(shè)為1讀入虛部的數(shù)據(jù)并重新進(jìn)行插值，但這次將虛部插值寫入另一個文件。8.至此，所有的插值任務(wù)完成，退出POST1 FINISH并讀入子模型數(shù)據(jù)庫(RESUM或Ma i nMenu>F i1e >Re sumefrom)。第四步：分析子模型在本步中，用戶指定分析類型和分析選項，加入插值的DOF數(shù)

14、值(和溫度數(shù)值),施加其他的載荷和邊界條件，指定載荷步選項，并對子模型求解。第一步是進(jìn)入求解器(/SOLU或Ma inMen u>Solution )。然后定義分析類型(一般為靜態(tài))和分析選項。要施加切割邊界自由度約束，用下列命令讀入CBDO命令生成的由D命令組成的文件：Command:/INPUTGUI:UtilityMenu>File>ReadInputfrom要施加溫度插值，用下列命令讀入BFINT命令生成的由BF命令組成的文件：Command:/INPUTGUI:UtilityMenu>File>ReadInputfrom如果數(shù)據(jù)有實部和虛部，先讀入實部數(shù)

15、據(jù)文件，指定自由度約束數(shù)值和(/或)結(jié)點體載荷是否計算，然后讀入虛部數(shù)據(jù)文件。用下列方法指定計算自由度約束數(shù)值：精心整理Command:DCUM, ADDGUI:MainMenu>Preprocessor>Loads>-Loads-Settings>ContraintsMainMenu>Solution>-Loads-Settings>Constraints用下列方法指定計算結(jié)點體載荷數(shù)值：Command:BFCUM, ADDGUI:MainMenu>Preprocessor>Loads>-Loads-Settings>Noda

16、lBodyLdMainMenu>Solution>-Loads-Settings>NodalBodyLd注意在執(zhí)行DCUM 口 BFCUM命令時要先將其初始狀態(tài)設(shè)為初始值。重要的一點是要將粗糙模型上所有其他載荷和邊界條件復(fù)制到子模型上。比如對稱邊界條件，面力，慣性載荷（如重量），集中力等（見圖5- 5） o圖5-5子模型的載荷然后指定載荷步選項（如輸出控制）并開始計算：Command:SOLVEGUI:MainMenu>Solution>CurrentLSMainMenu>Solution>RunFLOTRAN在求解完成后，退出SOLUTION FIN

17、ISH子模型的數(shù)據(jù)流向（無溫度插值）見圖5-6o圖5-6子模型分析（無溫度插值）的數(shù)據(jù)流向第五步：驗證切割邊界和應(yīng)力集中位置的距離是否足夠最后一步是驗證子模型切割邊界是否遠(yuǎn)離應(yīng)力集中部分。可以通過比較切割邊界上的結(jié)果（應(yīng)力，磁通密度等）與粗糙模型相應(yīng)位置的結(jié)果是否一致來驗證。如果結(jié)果符合得很好，證明切割邊界的選取是正確的。如果不符合的話，就要重新定義離感興趣部分更遠(yuǎn)一些的切割邊界重新生成和計算子模型。一個比較結(jié)果的有效方法是使用云圖顯示和路徑顯示，見圖5 7和5-8o圖5-7比較結(jié)果時的云圖顯示圖5-8比較結(jié)果時的路徑顯示輸入示例下面列出了 一個子模型分析的輸入示例：!開始子模型分析/FILN

18、AME, coarse!工作文件名為 coarse/PREP7!進(jìn)入 PREP7.!生成粗糙模型FINISH/SOLU!進(jìn)入求解器ANTYPE,. . !分析類型和分析選項D, . !載荷和載荷步選項DSYMM,.ACEL,.SAVE粗糙模型數(shù)據(jù)庫文件coarse, db SOLVE求解粗糙模型!結(jié)果在文件coarse, rst （或rmg等） FINISH !生成子模型/CLEAR!清除數(shù)據(jù)庫(或退出ANSYS并重新進(jìn)入)精心整理/FILNAME, submod新工作文件名為submod/PREP7重新進(jìn)入PREP7.!生成子模型!進(jìn)行切割邊界插值NSEL,. !選擇切割邊界上的結(jié)點NWRI

19、TE將其寫入文件submod. nodeALLSEL讀入所有實體NWRITE, temp, node!將所有結(jié)點寫入文件temp, node(用于溫度插值)SAVE存儲子模型數(shù)據(jù)庫文件submod. db FINISHRESUME, coarse, db!讀入粗糙模型數(shù)據(jù)庫(coarse, db)/POST1!進(jìn)入 POST1FILE, coarse, rst!使用粗糙模型結(jié)果文件SET,.！讀入需要的結(jié)果數(shù)據(jù)CBDOF從submod. node中讀入切割邊界結(jié)點并將D命令寫入submod. cbdoBFINT, temps, node!從temps, node中讀入所有子模型結(jié)點并將BF命令

20、寫入文件submod. bfin (用于溫度插值)FINISH!結(jié)束插值過程RESUME!入子模型數(shù)據(jù)庫(submod. db)/SOLU!進(jìn)入求解器ANTYPE,. . !分析類型和選項/INPUT, submod, cbdo!切割邊界自由度/INPUT, submod, bfin!溫度插值DSYMM,. !其他載荷和載荷步選項ACEL,.SOLVE 子模型求解 FINISH /POST1!進(jìn)入 POST1.!驗證子模型數(shù)據(jù)FINISH殼到體子模型在殼到體子模型技術(shù)中，粗糙模型為殼模型而子模型為三維實體模型。圖5-9所示為三維實體子模型添加到粗糙殼模型上的例子。圖5-93-D實體模型疊加到殼

21、單元模型上殼到體子模型分析與體對體子模型分析大致一致。下面的幾點是要記住的：殼到體子模型分析是將DBD 0命令(MainMenu>GeneralPostproc>Submodeling>lnterpolateDOF )和 BFINT 命令(MainMenu>General>Submodeling>lnterpBodyForc )中 KSHS 域設(shè)為 1 實現(xiàn)的。本特性不適用于SHELL91 或 SHELL99 ( KEY0PTU)不等于 0)。精心整理子模型切割邊界應(yīng)為垂直于殼平面的端面（見圖5 10）。切割邊界上的結(jié)點寫入文件中NWRITE （MainMe

22、nu>Preprocessor>Create>Node>WriteNodeFile）。要確定切割邊界上結(jié)點的自由度數(shù)值CBD0F,程序首先將結(jié)點延伸到殼平面的最近的單元上， 該延伸結(jié)點的自由度數(shù)值就插值并賦值給相應(yīng)的結(jié)點。溫度插值是由計算最近殼單元的中面平均溫度得 到的。注一切割邊界上的結(jié)點位置必須在最近殼單元平均厚度的0. 75倍之間，見圖5-10o也就是 說， 子模型應(yīng)大致在粗糙模型的中間。圖5-10結(jié)點旋轉(zhuǎn):a） CBDOF命令之前,b） CBDOF命令之后在結(jié)構(gòu)分析中，切割邊界結(jié)點只計算平動自由度，但其數(shù)值是根據(jù)延伸結(jié)點的平動和轉(zhuǎn)動自由度得到的。而且，結(jié)點旋轉(zhuǎn)以

23、使結(jié)點的UY自由度始終垂直于殼平面，見圖510。UY自由度只有當(dāng)結(jié)點在殼平面平均厚度的10%之內(nèi)時才計算。這防止了子模型在反向的過度位移。由CBD0命令寫的.CBD0文件包括兩個部分：（1） 一組NM0DI命令（表示結(jié)點旋轉(zhuǎn)角度）和 DDELE命令（刪除UY約束），（2）一組D命令（施加自由度插值）。這兩個部分用/EOF命令和一 個:CBnn標(biāo)記分開（nn為結(jié)果序列迭代次數(shù)）。用戶必須將.CBDC文件讀入PREP沖，因為NMODI命令只能在PREP沖適用。要完成這步操 作，進(jìn)入 前處理器，然后用下列命令：Command:/INPUT GUI:UtilityMenu>File>Rea

24、d!nputfrom同時，要讀入.CBDO文件兩次，因為兩個命令部分被/EOF命令分開了。在第二次讀入文件時，用/INPUT命令的LINE域指定程序從:CBnn處開始讀入，見下：/PREP7LCBD0文件必須在PREP7中讀入/INPUT, cbdo!讀入 Jobname, cbdo 到/EOF 處 /INPUT, dbdo, :cbl!從：cbl處讀入同一文件 子模型分析實例（命令行格式）問題描述求解矩形平板中心開孔，承受橫向拉力時的應(yīng)力 集中情況。材料特性和模型的幾何形狀見下圖。/FILNAM, coarse/PREP7smrt, off/TITLE,STRESSCONCENTRATION

25、ATAHOLEINAPLATE/NOPRANTYPE, STATIC 靜力分析ET, 1, PLANE2MP, EX, 1, 30E6MP, NUXY, 1,0.3LProblem SketchCoarse Model - Area and KeypointsRefined Keypoint and Area Model Around Hole (Submodel)K, 1,6!關(guān)鍵點K, 2, 6, 6K, 3, 6K, 4,. 5K, 5,. 5K, 6L, 1,2L, 2, 3L, 3,4LESIZE, 3, 4,. 25!定義線3的分段數(shù)LARC, 4, 5, 6, 0. 5LESIZ

26、E, 4, 6!定義線4的分段數(shù)L, 5,1LESIZE, 5, 4, 4!定義線5的分段數(shù)AL, 1, 2, 3, 4, 5ESIZE, 4!每條線分為4段AMESH,ALL/AUTO, 1/PLOPTS, INFO, 0/PLOPTS, WINS, 0/WINDOW, LTOPLSEL, S, LINE, 3, 5, 2 DL, ALL, SYMM LSEL, S, LINE, 1 NSLL, 1SF, ALL, PRES,-1OOO. !施加拉力LSEL, ALLNSEL, ALLCSYS, 1FINISH/SOLUSOLVEFINISHSAVE存儲文件VM142. DB/P0ST1SE

27、T, 1, 1NSORT, S, X, 3PRNSOL, S, COMP/WINDOW, 1, OFF/NOERASE/DSCALE, 2, 1/WINDOW, 2, RTOPPLNSOL, S, X* GET, CRSESTR, NODE, 18, S, X"STATUS* DIM, LABEL, CHAR, 1, 2* DIM, VALUE, ,1,3LABEL (1, 1)= MXSTR,LABEL (1,2)= CRSMODEL，*VFILL, VALUE(1, 1), DATA, 3018* VFILL, VALUE(1, 2), DATA, CRSESTR*VFILL,

28、VALUE(1, 3), DATA, ABS(CRSESTR/3018)SAVE, TABLE_1FINISH/CLEAR, NOSTART!CLEARTHEDATABASE/FILNAM, SUBMODE定義子模型文件名/PREP7 smrt, off /NOPR/TITLE, STRESSCONCENTRATIONATAHOLEINAPLATE ANTYPE, STATIC ET, 1, PLANE42MP, EX, 1, 30E6MP, NUXY, 1,0.3CSYS, 1K, 10,. 5,45K, 11,. 5, 90K, 12, 1. 5,45K, 13, 1.5,90A, 10,

29、 12, 13, 11ESIZE,8MSHK, 1MSHA, 0, 2DAMESH, 1/WINDOW, 2, OFF/NOERASE/PLOPTS, INFO, 0/PLOPTS, WINS, 0/WINDOW, 1, LTOP/USER!用戶定義的比例/DIST, 1, 3. 3/FOCUS, 1,3,3EPLOTLSEL, S, LINE, 1, 2NSLL, 1!選擇切割邊界的結(jié)點NWRITE將結(jié)點寫入文件SUBMODEL. NODELSEL, ALLNSEL, ALLFINISHSAVE存儲子模型數(shù)據(jù)庫文件SUBMODEL. DB/P0ST1RESUME, coarse, dbFIL

30、E, coarse, rstCBD0F,0,0!激活邊界條件插值FINISH/PREP7 smrt, offRESUME!文件SUBMODEL. D中讀入子模型/NOPR/INPUT, cbdo, :cbl!從文件 SUBMODEL. CBDO 讀入插值邊界條件 /GOPRLSEL, S, LINE, 3!施加其余的邊界條件DL, ALL, SYMMFINISH/SOLUSOLVEFINISH/P0ST1SET, 1, 1NSORT, S, X, 3PRNSOL, S, COMP/WINDOW, 1,OFF/AUTO, 3/WINDOW, 3, BOT/NOERASE/PLOPTS, MINM, 1/USER, 3/DIST, 3,. 2 /FOCUS, 3,. 2,. 5/CONTOUR, 3, AUTOPLNSOL,SX*GET, SUBSTR, NODE, 18, S, X *DIM, LABEL, CHAR, 1, 2 *DIM, VALUE, 1, 3LABEL (1, 1)=，MAXSTRS，LABEL(1, 2)=' SUBMOD' *VF