ansys高級(jí)分析子結(jié)構(gòu)adv4

1、第四章子結(jié)構(gòu)什么是子結(jié)構(gòu)？子結(jié)構(gòu)就是將一組單元用矩陣凝聚為一個(gè)單元的過(guò)程。這個(gè)單一的矩陣單 元稱為超單元。在ansys分析屮，超單元可以象其他單元類型一樣使用。唯一 的區(qū)別就是必須先進(jìn)行結(jié)構(gòu)生成分析以生成超單元。子結(jié)構(gòu)可以在 ansys/mutiphysics, ansys/mechanical 和 ansys/structural 中使用。使用子結(jié)構(gòu)主要是為了節(jié)省機(jī)吋，并且允許在比較有限的計(jì)算機(jī)設(shè)備資源 的基礎(chǔ)上求解超大規(guī)模的問(wèn)題。原因之一如a）非線性分析和帶宥大量重復(fù)幾何結(jié) 構(gòu)的分析。在非線性分析中，可以將模型線性部分作成子結(jié)構(gòu)，這樣這部分的單 元矩陣就不用在非線性迭代過(guò)程中重復(fù)計(jì)算。在有

2、重復(fù)兒何結(jié)構(gòu)的模型中（如有 四條腿的桌子），可以對(duì)于重復(fù)的部分生成超爭(zhēng)元，然后將它拷w到不同的位置， 這樣做可以節(jié)省大量的機(jī)吋。子結(jié)構(gòu)還用于模型宥大轉(zhuǎn)動(dòng)的情況卜。對(duì)于這些模型，ansys假定每個(gè)結(jié) 構(gòu)都是圍繞其質(zhì)心轉(zhuǎn)動(dòng)的。在三維情況下，子結(jié)構(gòu)冇三個(gè)轉(zhuǎn)動(dòng)h由度和三個(gè)平動(dòng) 自由度。在大轉(zhuǎn)動(dòng)模型中，用戶在使用部分之前無(wú)須對(duì)子結(jié)構(gòu)施加約束，因?yàn)槊?個(gè)子結(jié)構(gòu)都是作為一個(gè)單元進(jìn)行處理，是允許剛體位移的。另外一個(gè)原因b）個(gè)問(wèn)題就波前大小和需用磁盤空間來(lái)說(shuō)相對(duì)于一個(gè)計(jì)算 機(jī)系統(tǒng)太龐大了。這樣，用戶可以通過(guò)子結(jié)構(gòu)將問(wèn)題分塊進(jìn)行分析，每一塊對(duì)于 計(jì)算機(jī)系統(tǒng)來(lái)說(shuō)都是可以計(jì)算的。如何使用子結(jié)構(gòu)子結(jié)構(gòu)分析宥以卜'

3、;三個(gè)步驟：生成部分使用部分?jǐn)U展部分生成部分就是將普通的冇限元單元凝聚為一個(gè)超單元。凝聚是通過(guò)定義一 組主cl由度來(lái)實(shí)現(xiàn)的。主g由度用于定義超單元與模型中其他單元的邊界，提取 模型的動(dòng)力學(xué)特性。閣4-1是一個(gè)板狀構(gòu)件用接觸單元分析的示意。由于接觸單 元需要迭代計(jì)算，將板狀構(gòu)件形成子結(jié)構(gòu)將敁著地節(jié)省機(jī)時(shí)。本例中，主自由度 是板與接觸單元相連的自由度。圖4-1子結(jié)構(gòu)使用示例使用部分就是將超單元與模型整體相連進(jìn)行分析的部分。整個(gè)模型可以是 一個(gè)超單元，也可以象上例一樣是超單元與非超單元相連的。使用部分的計(jì)算只 是超單元的凝聚（自由度計(jì)算僅限于主自由度）和非超單元的全部計(jì)算。擴(kuò)展部分就是從凝聚計(jì)算結(jié)果

4、開始計(jì)算整個(gè)超笮元中所有的自由度。如果 在使用部分有多個(gè)超單元，那么每個(gè)超單元都需要有單獨(dú)的擴(kuò)展過(guò)程。圖4-2示出了整個(gè)子結(jié)構(gòu)分析的數(shù)裾流向和所用的文件。三個(gè)步驟的詳細(xì) 解釋見以后的敘述。圖4-2典型子結(jié)構(gòu)分析屮的數(shù)據(jù)流向生成部分：生成超單元本部分主要宥兩步：1. 建立模型。2. 施加邊界條件，生成超單元矩陣。第一步：建立模型在這一步中，指定文件名和分析名稱，用prep7定義單元類型，單元實(shí)參, 材料特性和模型幾何結(jié)構(gòu)。這些任務(wù)在ansys絕大多數(shù)分析屮都是通用的，在 ansys basic analysis procedures guide冇所敘述。在生成部分，需要記住以下幾 點(diǎn)：文件名在子

5、結(jié)構(gòu)分析中很有用處。有效地使用文件名，在三部分分析 中可以省略很多文件處理操作。用以下方法指定文件名：command: /filenamegui: utility menufilechange jobname如：/filename,gen將生成過(guò)程中所有文件名都定義為gen。缺省的文件名是file （或file） 或在進(jìn)入ansys后定義的任意文件名。單元類型ansys提供的絕大多數(shù)單元都可以用來(lái)生成超單元。唯一的 限制是單元必須是線性的。如采生成超單元時(shí)宥雙線性單元的話，ansys將cj 動(dòng)作為線性單元處理。注意：在直接耦合中帶載荷向量的耦合單元是不能做子結(jié)構(gòu)分析的?？梢?用同種形狀的單元來(lái)

6、替代。細(xì)節(jié)參看ansys coupled-field analysis guide。材料特性一一定義所冇必須的材料特性。例如，如果生成質(zhì)量矩陣，就必 須定義密度或其他形式的質(zhì)量：如采要生成熱傳導(dǎo)矩陣，就要定義比熱。同樣， 超單元是線性的，非線性材料將被忽略。模型生成一一在生成部分，主要生成模型的超單元部分。非超單元部分是 在以后的使用部分生成的。是，在建模的開始就需要對(duì)模型的兩個(gè)部分奮所規(guī) 劃，主要是確定超單元部分和非超單元部分如何連接。為了保證連接正確，應(yīng)該 保證接觸部分結(jié)點(diǎn)號(hào)一致。（其他可以方便用戶的方法在本章“使用部分” 一節(jié) 還有介紹。）要生成整體模型應(yīng)該這樣做：將模型存儲(chǔ)在數(shù)據(jù)厙文件

7、中，選擇子結(jié)構(gòu)部 分進(jìn)行生成計(jì)算。在以后的使用部分，resumecutility menu>file>resume from） 數(shù)據(jù)庫(kù)文件，不選（unselect）子結(jié)構(gòu)，用超單元矩陣代替。第二步：施加邊界條件，生成超單元矩陣。生成部分的結(jié)果含超單元矩陣。象其他分析一樣，用戶要定義分析類型 和分析設(shè)置，施加邊界條件，定義載荷步，開始計(jì)算。如何完成這些工作見下面 的敘述：1. 進(jìn)入求解器：command: /solugui: main menu>solution2. 定義分析類型和分析設(shè)置：分析類型選擇生成超單元使用下列方法：command: antypegui: main m

8、enusolution-analysis type-new analysis新的分析或重啟動(dòng)如果是開始一個(gè)新的分析時(shí)，只要指定分析類型（如 上所述）即可。如果是重扁動(dòng)計(jì)算，必須在antype命令屮設(shè)定status=rest（main menu>solution>-analysis type-restart）o 如來(lái)要另外施加載荷時(shí)，口j 以 用重啟動(dòng)。（重動(dòng)時(shí)，初始運(yùn)算后的jobname.emat, jobname.esav和 jobname.db文件要存在。）超單元矩陣文件名指定超單元矩陣文件名（sename）。程序?qū)⒆詣?dòng)添加 后綴.sub,因此完整的文件名是sename.sub

9、。缺省是使用工作文件名 /filename定義超單元矩陣文件名，可以使用以p命令：command: seoptgui: main menu>solution>analysis options要生成的矩陣可以指定僅生成剛度矩陣（或傳導(dǎo)矩陣，電磁系數(shù)矩陣）; 生成剛度和質(zhì)量矩陣（或熱傳導(dǎo)矩陣等）；生成剛度，質(zhì)量和阻尼矩陣。質(zhì)量矩 陣用于結(jié)構(gòu)動(dòng)力學(xué)分析和在使用部分宥慣性載荷的情況k。在熱分析中，只宥瞬 態(tài)熱分析j用到熱傳導(dǎo)矩陣。對(duì)于其他分析和阻尼矩陣也大同小異。用seopt 命令或其gui路徑來(lái)定義。輸出矩陣這個(gè)選項(xiàng)允許輸出超申元矩陣?？梢灾付ㄝ敵鼍仃嚭洼d荷向 量，也吋以只輸出載荷m量。

10、缺省值是不輸出任何矩陣。要輸出矩陣，用seopt 命令或其相應(yīng)的gui路徑。質(zhì)量矩陣形成只在想生成質(zhì)量矩陣時(shí)使用。用戶可以選擇缺杏生成（取 決于所用單元類型）或集中質(zhì)量近似。對(duì)于絕大多數(shù)情況，推薦使用缺省生成的 方式。但是，在極薄構(gòu)件的分析屮，如細(xì)長(zhǎng)桿或極薄殼體，集屮質(zhì)量近似將得到 更好的結(jié)果。用下列方法指定集中質(zhì)量近似：command: lumpmgui: main menu>solution>analysis options3. 用下列方法定義主自由度：command: mgui: main menu>solution>master dofs>define在子

11、結(jié)構(gòu)中，主cj由度有四種作用：a. 它們作為超單元與非超單元的邊界。應(yīng)保證將超單元與非超單元接 觸的結(jié)點(diǎn)自由度都定義為超單元（在m命令中l(wèi)ab=all）,如圖4-1所示。當(dāng)模 型屮只有超單元時(shí)同樣要定義主自由度。b. 如果在動(dòng)力學(xué)分析中使用超單元，那么主自由度規(guī)定了結(jié)果的動(dòng)力 學(xué)特性。在 ansys structural analysis guide 第三章的 “matrix reduction” 中吻 所說(shuō)明。c. 如果在使用部分要施加約束d或集中力f時(shí)，這些位置結(jié)點(diǎn)的自由 度也要定義為主自由度。d. 在大位移情況下nlgeom, on （ main mene>solution>

12、analysis options ）的使用部分需要主q由度，或者在使用setran命令（main menu>preprocessor>create>elements>-superelements-by cs tranfer）時(shí)。在這些 時(shí)候，所有主自由度的結(jié)點(diǎn)都要定義6個(gè)方向的自由度（ux, uy, ux, rotx, roty, rotz）o4. 施加邊界條件。在生成部分可以施加所冇的載荷類型，但有以下幾 點(diǎn)情況需耍注意：程序?qū)生成一個(gè)包括所有施加的載荷的等效載荷向量。每個(gè)載荷步一個(gè) 載荷向量將寫入超單元矩陣文件中。載荷向量的最大允許值是31個(gè)。在生成部分可以使用非

13、零的自由度約朿并作為載荷向量的一部分。（在 擴(kuò)展部分，如來(lái)被擴(kuò)展的載荷步中宥非零的g由度約束時(shí)，數(shù)裾庫(kù)中必須宥相應(yīng) 的自由度數(shù)值。如果沒有的話，自由度約束就耍在擴(kuò)展部分重新定義。自由度約束和集中載荷的施加可以推遲到使用部分，但在這些位置的主 自由度一定要先定義好。同樣，線或角加速度的施加也吋以推遲到使用部分，但只冇在生成質(zhì)量 矩陣時(shí)。如果在使用部分要旋轉(zhuǎn)超單元時(shí)，推薦使用這種作法，因?yàn)榇藭r(shí)載荷向 量的方向是“凍結(jié)”的，它隨著超單元旋轉(zhuǎn)。 maxwell力的標(biāo)志通常用于電磁分析中，用來(lái)標(biāo)記計(jì)算哪個(gè)單元面上的 電磁力分布。這個(gè)標(biāo)志在電磁子結(jié)構(gòu)分析屮無(wú)效，因此不要使用。注如果生成質(zhì)量矩陣，建議在使用

14、部分對(duì)主自由度（在生成部分定義） 施加約束。這樣就保證所宥的質(zhì)量都包含在子結(jié)構(gòu)中。5. 定義載荷步選項(xiàng)。子結(jié)構(gòu)生成部分只能使用動(dòng)力選項(xiàng)（阻尼）。阻尼（動(dòng)力選項(xiàng)）只在生成阻尼矩陣時(shí)可用。指定質(zhì)量（alpha）阻尼使用下列方法：command: alphadgui: main menu>solution>time/frequenc>damping指定剛度（beta）阻尼使用下列方法：command: betadgui: main menu>solution>time/frequenc>damping指定與材料冇關(guān)的beta阻尼使用下列方法：command: mp

15、, dampgui: main menu>preprocessor>material props>-constan-isotropicmain menu>preprocessor>materiai props>polynomiai表4-1子結(jié)構(gòu)中吋以施加的載荷load nameload categorycommands*solid model loadsfinite element loadsdisplacementtemperaturemag. potentialetc.conslraintsdk,dklist，dkdele，dl,dllist, dlde

16、le， da, dalist, dadele, dtrand, dsymm, dust,ddele,scale, dcumforceheat flow ratewag. fluxetc.forcesfk, fklist, fkdele, ftranf, fust, fdele,fscale, fcumpressureconvectionmaxwell surfaceetc.surface loadssfl, sfllis, sfldel, sfa, sfalis，sfadel, sfgrad，sftransf, sflist, sfdele,sfe, sfelis, sfedel,sf巳eaw

17、，sfgrad,sffun，sfscale,sfcuwtemperatureheat generation rate current densityetc.body loads巳f(wàn)k，bfklis,bfkdel,巳f(wàn)tranbf, bflist, bfdele,bfe, bfelis，bfedel, bfscale，bfcumgravity, linear and angular accelerationinertia loadsacel, domega6. 存儲(chǔ)數(shù)據(jù)庫(kù)的備份文件。注這樣做的目的是需耍在擴(kuò)展部分使用同樣的數(shù)據(jù)庫(kù)文件。用下列方法完成 本步操作：command: savegui:

18、utility menu>file>save as jobname.db7. 開始計(jì)算：command: solvegui: main menu>solution>current ls計(jì)算結(jié)果包括超爭(zhēng)元矩陣文件，sename.sub, sename是通過(guò)seopt1指定 的文件名或是工作文件名/filename。矩陣文件包括根據(jù)施加的載荷計(jì)算出的 載荷向量。（如果沒宥施加載荷，載荷向量將為零。）8. 如冇另外的載荷步，重復(fù)步驟7來(lái)生成其他的載荷向量。載荷向量 的號(hào)碼是遞增的，并添加到同一個(gè)超單元矩陣文件內(nèi)。關(guān)于多載荷步的其他方法 見 ansys basic analys

19、is procedures guide。9. 退出 solution：command: finishgui: main menufinish使用部分：使用超單元使用部分可以是任何ansys分析類型（flotran和顯式動(dòng)力分析除外）。 和普通分析的區(qū)別就是一個(gè)或幾個(gè)單元是前面生成的超單元。每個(gè)單獨(dú)的分析指 南屮都有做不同分析的詳細(xì)介紹。在這一部分，我們主要介紹如何將超單元變成 模型的一部分。這個(gè)過(guò)程有以下兒個(gè)步驟：1. 清除數(shù)據(jù)庫(kù)并指定一個(gè)新的工作文件名。2. 建立模型。3. 施加邊界條件丼求解。第一步：清除數(shù)據(jù)庫(kù)并指定一個(gè)新的工作文件名。使用部分包含新的模型和新的邊界條件。因此，第一步是清除

20、現(xiàn)存的數(shù)據(jù) 庫(kù)。這與退出并重新進(jìn)入ansys的效果是一樣的。清除數(shù)據(jù)庫(kù)用卜*列方法：command: /cleargui: utility menu>file>clear&start new缺省情況下，清除數(shù)據(jù)庫(kù)就會(huì)重新讀入start.ans文件。（可以改變這個(gè) 設(shè)置）注：如采通過(guò)命令行輸入來(lái)清除數(shù)據(jù)庫(kù)時(shí)，在/clear命令行中不能宥其 他的命令。新定義的文件名要與生成部分使用的文件名不同。這樣，生成部分的文件 就不會(huì)被覆蓋。用下列方式之一定義新的工作文件名：command: /filnamegui: utility menufilechange jobname第二步：建立

21、模型本步是在prep7屮實(shí)現(xiàn)的。主要完成以下任務(wù)：1. 定義matrix50 （超單元）為一種單元。用以下方法：command: etgui: main menu>preprocessor>element type>add/edit/delete2. 定義其他非超單元的單元類型。非線性可以使用，取決于進(jìn)行分析 的類型。3. 定義非超單元的單元實(shí)參和材料特性。非線性可以使用，取決于進(jìn) 行分析的類型。4. 定義非超單元的幾何形狀。在定義與超單元接觸部分時(shí)應(yīng)非常注意。 其結(jié)點(diǎn)位置要精確重合。（見圖4-3）圖4-3超單元與非超單元的接觸處結(jié)點(diǎn)應(yīng)與主自由度精確重合q nodes for

22、 contact elemer master nodes on superek共有三種方法保證結(jié)點(diǎn)重合： 使用與生成部分同樣的結(jié)點(diǎn)號(hào)碼。 在生成部分的接觸部分結(jié)點(diǎn)和使用部分的接觸結(jié)點(diǎn)使用相同的結(jié)點(diǎn) 號(hào)碼增值（或平移）。（使用setran,在第五步中說(shuō)明） 將這兩部分結(jié)點(diǎn)所有的自由度固連起來(lái)。（用cp系列命令）在不能使 用前w種方法時(shí)，這種方法是有效的。定義自由度固連可以用下列方法：command: cpgui: main menu>preprocessor>coupling/ceqn>couple dofs如果超單元與其他單元不相連時(shí)，在使用部分可以不定義任何單元。5. 用正

23、確的單元類型定義超單元，讀入超單元矩陣。選擇正確的爭(zhēng)元類型：command: typegui: main menu>preprocessor>-modeling-create>elements>elem attributes用下列方法讀入超單元矩陣（可能要先用其他命令修改矩陣，見后）：command: segui: main menu>preprocessor>create>elements>-superlements->from .subfilea. 如果模型中不包含非超單元，或是包含非超單元但接觸處的結(jié)點(diǎn)號(hào)與 主自由度結(jié)點(diǎn)號(hào)一致，那么可

24、以用se命令直接讀入超單元：type,.!單元類型號(hào)se, gen!從文件gen.sub屮讀入超爭(zhēng)元se命令的文件名域表示超單元矩陣的文件名。擴(kuò)展名.sub是給定的，因 此整個(gè)文件名就是sename.sub（gen.sub如上例）。超單元就給定k一個(gè)可用的 單元號(hào)。b. 如果模型中有非超單元存在，并且接觸處的結(jié)點(diǎn)號(hào)與主自由度的結(jié)點(diǎn) 號(hào)是有一個(gè)偏移量的話，就必須先用新的結(jié)點(diǎn)號(hào)生成一個(gè)新的超爭(zhēng)元矩陣，然后 讀入新的矩陣。用卜列方法生成新的超單元組陣：command: setrangui: main menu>preprocessor>create>elements>-sup

25、erelements-by cs transfermain menu>preprocessor>create>elements>-superelements-by geom offset 用下列方法讀入新的矩陣：command: segui: main menu>preprocessor>create>eiements>-supreeiements-from .sub file例如，給出超單元矩陣文件gen.sub和新的結(jié)點(diǎn)偏移量為2000,命令如下:settran，gen，2000，gen2,sub !生成新的超單元 gen2.sub,偏移 量

26、為2000type，.！單元類型號(hào)se，gen2!讀入新的超爭(zhēng)元文件gen2.subc. 如果模型中包含非超單元而且接觸處結(jié)點(diǎn)與主自由度結(jié)點(diǎn)無(wú)任何關(guān)系 （一般是結(jié)點(diǎn)cl由生成的情況），那么先看卞面的注意事項(xiàng)：注意：生成部分的主自由度的結(jié)點(diǎn)號(hào)往往會(huì)覆蓋使用部分的模型結(jié)點(diǎn)號(hào)。 這樣，讀入超單元將覆蓋使用部分的模型結(jié)點(diǎn)。為了避免這種覆蓋，在讀入超單 元前使用結(jié)點(diǎn)號(hào)偏移。在這種情況下，在使用se命令前先存儲(chǔ)數(shù)據(jù)庫(kù)文件。因此，應(yīng)該先存儲(chǔ)數(shù)據(jù)庫(kù)文件save,再用setran命令用新的結(jié)點(diǎn)偏移 量生成新的超單元矩陣，最用se命令讀入新的矩陣。用crintf命令（main menu>preprocess

27、or>coupling/ceqn>coincident nodes）將接觸處的結(jié)點(diǎn)連接。例 如，給出超單元矩陣文件gen.sub:*get，maxnod，node，num，max !maxnod=最大結(jié)點(diǎn)號(hào)setran，gen，maxnod，gen2，sub!生成新的超單元，結(jié)點(diǎn)偏移量為se,gen2 nsel， cpintf,allmaxnod,文件名為 gen2.sub !讀入新的超單元 !選擇接觸處的所脊結(jié)點(diǎn) !固連所冇結(jié)點(diǎn)的所冇自屮度nsel,alld. 如果要對(duì)超單元進(jìn)行移動(dòng)或拷w到不同的位置，或進(jìn)行鏡向操作，必 須用 setran 命令 sesymm 命令（main m

28、enu>preprocessor>create>elements> by reflection）,并給出正確的結(jié)點(diǎn)偏移量，生成新的超單元矩陣文件，然后用 se命令讀入新的矩陣。連接超單元與非超單元同樣用上面的方法用相同的 結(jié)點(diǎn)號(hào)，結(jié)點(diǎn)偏移量，或用cpintf命令。注如果用setran命令將超單元傳遞到不同的坐標(biāo)系中，缺哲是超單 元的主自由度隨之旋轉(zhuǎn)。這一點(diǎn)在初始的超單元主自由度發(fā)生旋轉(zhuǎn)時(shí)很有用，比 如旋轉(zhuǎn)到柱華標(biāo)系屮。（在這種情況下，傳遞不影響超單元?jiǎng)偠染仃嚒＃┤绻跏?超單元沒冇結(jié)點(diǎn)發(fā)生旋轉(zhuǎn)，那么傳遞的超單元也無(wú)須旋轉(zhuǎn)結(jié)點(diǎn)。用戶可以通過(guò)設(shè) 定setran命令的noro

29、t域?yàn)?來(lái)防止結(jié)點(diǎn)旋轉(zhuǎn)。（在這種傳遞中，超單元的 剛度矩陣和載荷向量由程序自動(dòng)修改）。6. 用圖形敁示和列表功能驗(yàn)證超單元的位置是否正確。超單元用輪廓 線表示，數(shù)據(jù)在生成部分寫到矩陣文件屮。要生成圖形顯示，用下列命令：command: eplotgui: utility menuplotelements要生成列表，用下列命令：command: selistgui: utility menu>list>other>superelem data7. 存儲(chǔ)整個(gè)模型文件：command: savegui: utility menufilesave as jobname.db退出pre

30、p7：command: finishgui: main menu>finish第三步：施加邊界條件并求解。木步在求解器屮完成。取得使用部分結(jié)果的過(guò)程取決于所做的分析類型。 以上提到，可以在絕大多數(shù)分析中使用超單元。用戶應(yīng)當(dāng)冇相應(yīng)的通過(guò)生成部分 生成的矩陣。例如，如來(lái)耍做結(jié)構(gòu)動(dòng)力學(xué)分析，必須宥質(zhì)量矩陣。其體過(guò)程如卜、1. 進(jìn)入 solution。command: /solugui: main menu>solution2. 定義分析類型和分析選項(xiàng)。對(duì)于大轉(zhuǎn)動(dòng)分析，打開大位移效果nlgeom，on,并正確定義非線性分 析的子步數(shù)。3. 在非超單元上施加邊界條件。主要拈自由度約束和對(duì)稱邊

31、界，集 屮載荷，平面載荷，體載荷，和慣性載荷。注意慣性載荷從有在通過(guò)生成部分生成了質(zhì)量矩陣時(shí)才生效。注對(duì)大轉(zhuǎn)動(dòng)分析，在本步屮要施加正確的約束條件。4. 用以下方法施加超單元載荷向量：command: sfegui: main menu>solution>-loads-apply>load vector>for superelement在超單元矩陣文件中，每個(gè)載荷步對(duì)放一個(gè)載荷向量，用參考號(hào)來(lái)區(qū)別。 如，sfe, 63，1, selv，0.75在63號(hào)單元上施加載荷向量，號(hào)為1,比例系數(shù)為0.75。因此，elem域 代表超單元的單元號(hào)，lkey域代表載荷m量號(hào)（缺省為1）

32、, lab是selv, val1 代表比例系數(shù)（缺省為0）。（參看sfe命令）注一一載荷向量的方向是固定在超單元上的，因此如果超單元旋轉(zhuǎn)了，載 荷句量就隨之旋轉(zhuǎn)。對(duì)于自由度方句也是如此（ux, uy, roty等），它們同 樣隨著超申元的旋轉(zhuǎn)而旋轉(zhuǎn)（除非setran命令屮norot=1,此吋結(jié)點(diǎn)坐標(biāo) 系不旋轉(zhuǎn)）。5. 定義與分析類型相適應(yīng)的載荷步選項(xiàng)。注用matrix50超單元時(shí)不要用pcg求解器。6. 開始計(jì)算：command: solvegui: main menu>solution>current ls本步計(jì)算包括非超單元的完整解和超單元的凝聚解主cj由度解。非超 單元的完整

33、解記錄在結(jié)果文件屮（jobname.rst，rth或rmg）,可以進(jìn)行普通的 后處理操作。凝聚解兄呆在文件jobname.dsub屮?？梢酝ㄟ^(guò)如下方法查看這個(gè)文件：command: sedlistgui: main menugeneral postproc>list resultssuperelem dofutility menulistresultssuperelem dof solu如果想得到超單元中所有自由度的解，就要用到擴(kuò)展部分。下面還會(huì)說(shuō)明。7. 退出 solution。command: finishgui: main menu>finish擴(kuò)展部分：在超單元中擴(kuò)展求解結(jié)

34、果擴(kuò)展部分從使用部分的凝聚解計(jì)算出整個(gè)超早元的完整解。下面說(shuō)明擴(kuò)展 部分的過(guò)程。本部分要求生成部分的.emat，.esav，.sub，.tri，.db和.seld文件, 使用部分的.dsub文件存在。如來(lái)在使用部分中使用了結(jié)點(diǎn)偏移，在擴(kuò)展部分 中將自動(dòng)計(jì)入。1. 清除數(shù)據(jù)庠。和當(dāng)于退岀并重新進(jìn)入ansys。用下列方法：command: /cleargui: utility menu>file>clear&start new2. 將文件名切換到生成部分的文件名。這樣，程序就可以識(shí)別擴(kuò)展部 分所用的文件。用下列方法：command: /filenamegui: utility

35、menu>file>change jobname3. 讀入生成部分的數(shù)據(jù)庫(kù)文件。用不列方法：command: resumegui: utility menu>file>resume jobname.db4. 進(jìn)入solution。用下列方法：command: /solugui: main menu>solution5. 激活擴(kuò)展部分及其選項(xiàng)。擴(kuò)展部分開關(guān)選為on。command: expassgui: main menu>solution>expasionpass 被擴(kuò)展的超單元名指定sename。command: seexpgui: main men

36、u>solution>expasionpass>expand superelem （完整的文件名假定為sename.sub）使用部分生成的凝聚解文件。用seexp命令（或其相應(yīng)的gui路徑）指 定該文件名。完整的文件名假定為usefil.dsub。位移的實(shí)部或虛部只冇在使用部分是諧波分析時(shí)使用。用seexp命令 （或其和應(yīng)的gui路徑）。被擴(kuò)展的解指定被擴(kuò)展的使用部分結(jié)果?？梢酝ㄟ^(guò)給出載荷步和子步， 也可以通過(guò)給出吋間或頻率來(lái)指定結(jié)果。用下列方法：command: expsolgui: main menusolutionexpansionpassby load step ma

37、in menu>solution>expansionpass>by time/freq注如果被擴(kuò)展的載荷步屮包含非零自由度約朿吋，數(shù)據(jù)庫(kù)中應(yīng)該有相 應(yīng)的自由度數(shù)值。如果沒冇，那么自由度約束必須在擴(kuò)展部分中重新聲明d。6. 指定載荷步選項(xiàng)。對(duì)于子結(jié)構(gòu)擴(kuò)展部分，只宥輸出控制選項(xiàng)是可用 的：輸出控制一一這呰選項(xiàng)控制打印輸出，數(shù)據(jù)庫(kù)和結(jié)果文件輸出和結(jié)果的插 值。如果在打印輸出文件（jobname.out）中包含某些內(nèi)容，使用下列方法： command: outprgui: main menu>solution>output ctrls>db/results file

38、 如果控制結(jié)果文件（jobname.rst）中的數(shù)據(jù)，使用下列方法：command: outresgui: main menu>solution>output ctrlsdb/results file如果想用將單元積分點(diǎn)上的數(shù)據(jù)拷貝到結(jié)點(diǎn)的方法代替插值方法（缺?。?時(shí)，用下列命令：command: eresxgui: main menu>solution>output ctrls>integration pt7. 開始擴(kuò)展部分：command: solvegui: main menu>solution>current ls8. 重復(fù)步驟5到7對(duì)其他使用

39、部分結(jié)果進(jìn)行擴(kuò)展。如果要擴(kuò)展不同超 單元的解，需要退出并重新進(jìn)入solution。9. 最jis,退出 solution。command: finishgui: main menufinish10. 用通用方法對(duì)超單元擴(kuò)展結(jié)果進(jìn)行后處理。 注擴(kuò)展部分不適用于功率譜密度分析。輸入示例下面給出一個(gè)子結(jié)構(gòu)分析的命令流文件示例。本例假定一個(gè)超單元與非超 單元連接。!生成部分!建模（超單元部分）/filname, gen/title, .!文件名為gen/prep7!進(jìn)入 i"rep7!生成超中元部分模型finish!加載并生成超單元矩陣 /sow!進(jìn)入 solutionantype, sub

40、st!子結(jié)構(gòu)分析seopt, gen,.!超單元名和其他子結(jié)構(gòu)分析選項(xiàng)9.!主0由度d,.!載荷。一個(gè)載荷向量將產(chǎn)生并寫入超單元矩 陣文件!載荷步選項(xiàng)save!存儲(chǔ)數(shù)據(jù)庫(kù)文件solve!開始求解生成gen.sub文件，包含超 單元矩陣載荷向量。!載荷。生成第二個(gè)載荷向量（d和m可以不變）solve!加入第二個(gè)載荷量。!重復(fù)加載和求解生成其他的載荷向量（最多達(dá)到31個(gè)）。finish!使用部分 !建模/clear!清除數(shù)據(jù)庫(kù)/filname, use!文件名為use/prep7!進(jìn)入prep7et, 1, matrix50! matrix50為超單元類型et, 2,.!非超單元的單元類型 !生成

41、非超單元部分模型type, 1!指向超單元類型setran, .!結(jié)點(diǎn)偏移量se,.!讀入由setran生成的超申元eplot nsel, . cpintf, all nsel, all finish !加載并求解 /soluantype, .sfe, .savesolvefinish!查看非超單元結(jié)果 !擴(kuò)展部分/clear/filname, genresume/soluexpass, on seexp, gen, usesolvefinish!檢杳超單元位置!選擇接觸處的結(jié)點(diǎn)!固連接觸處結(jié)點(diǎn)（在結(jié)點(diǎn)號(hào)不和等時(shí)使用）!進(jìn)入求解器 !分析類型和分析選項(xiàng)!非超單元上的載荷!施加超準(zhǔn)元載荷向ft!

42、載做步選項(xiàng)!存儲(chǔ)數(shù)據(jù)庫(kù)文件!開始求解，計(jì)算非超單元完整解和超單元凝聚 解。!清除數(shù)據(jù)庫(kù)!文件名切換到生成部分文件名!讀入生成部分?jǐn)?shù)據(jù)庫(kù)!進(jìn)入求解器!激活擴(kuò)展選項(xiàng)!要擴(kuò)展的超單元名!載荷步選項(xiàng)（主要是輸出控制）!丌始擴(kuò)展部分求解。超單元完整解寫入文件 gen.rst 或 rth, rmg 中。!查看超單元解請(qǐng)杳閱 antype, seopt, m, et, setran, se, cpintf, expass和seexp命令得到史加詳細(xì)的說(shuō)明。自頂而下子結(jié)構(gòu)分析上述的子結(jié)構(gòu)方法稱為自底而上的子結(jié)構(gòu)生成方法，每個(gè)超單元是通過(guò)獨(dú)立 的生成過(guò)程來(lái)生成，然盾通過(guò)使用過(guò)程組裝。這種方法適用于超大型的模型，

43、分 解為小的超單元便于求解。如果對(duì)相對(duì)小一些的模型或具有統(tǒng)一的兒何外型控制和不同的結(jié)構(gòu)分析的 情況，可以用一種與上述方法有一點(diǎn)細(xì)微差別的技術(shù)，稱之為自頂而下的子結(jié)構(gòu) 分析。這種方法適于將非線性模型中線性部分（應(yīng)相對(duì)計(jì)算機(jī)系統(tǒng)人小比較適中） 作成子結(jié)構(gòu)。這種做法的優(yōu)點(diǎn)是可以在處理中將幾個(gè)子結(jié)構(gòu)的結(jié)采加以組合。 自頂而下子結(jié)構(gòu)分析的步驟如下：1. 首先逑立整體模型，包拈超單元和非超單元部分。將模型存儲(chǔ)在一 個(gè)命名的數(shù)據(jù)庫(kù)文件屮。這個(gè)整體的數(shù)據(jù)庫(kù)在后面擴(kuò)展部分屮還要使用。如果模型中包含非超單元部分，在使用部分中也耍用到這個(gè)文件。2. 將模型的一部分選出完成生成部分操作。因?yàn)檎w模型都已經(jīng)建立 了，用

44、戶所要做的就是選擇要生成超單元部分的單元，施加載份（目的是生成載 荷向量），然后 solve 命令生成超 ¥ 元（main menu>solution>-solve-current ls）。3. 完成使用部分。進(jìn)入prep7,讀入整體模型數(shù)據(jù)庫(kù)并選擇非超單元 部分。接卜*來(lái)定義超單元類型et，type,讀入相應(yīng)的超單元矩陣。絕大多數(shù) 情況下，用戶不必?fù)?dān)心接觸處的結(jié)點(diǎn)連接，因?yàn)樗鼈兪怯梢粋€(gè)模型生成的。進(jìn)入 solution,定義分析類型和分析選項(xiàng)。在非超單元上施加載荷，讀入載荷向量（如果有的話），指定載荷步選項(xiàng)，并開始使用部分的求解。4. 完成擴(kuò)展部分。首先讀入整體模型，需

45、包含所冇結(jié)點(diǎn)和單元，然后 對(duì)每個(gè)超單元做擴(kuò)展，每次應(yīng)定義不同的文件名，并每次進(jìn)入和退出 solution。用戶可以用通用的后處理過(guò)程查看每個(gè)超單元的結(jié)果。注使用 整個(gè)數(shù)據(jù)庫(kù)文件，可以讀入多個(gè)超單元結(jié)果：resume,full,db/post1file,gen 1 set,.file,gen2set，.！不清除上一個(gè)超單元的結(jié)果。以下給出一個(gè)自頂而下子結(jié)構(gòu)分析的輸入示例。本例假定模型冇一個(gè)超單元和其他非超單元部分。!自頂而下子結(jié)構(gòu)分析的輸入示例 !建立整體模型/filname, full/title, . /prep7savefinish !生成部分 /filname, gen /soluant

46、ype, subst seopt, gen,. esel, .nselm 9.d,.solve!文件名為full!進(jìn)入prep7!生成整體模沏，包括超單元部分和非超單元 部分!存儲(chǔ)整體模型，以備（使用部分和）擴(kuò)展部 分使用!文件名為gen !進(jìn)入 solution !子結(jié)構(gòu)分析 ! 了結(jié)構(gòu)分析選項(xiàng) !選擇超單元部分的單元和結(jié)點(diǎn)!主自由度!載荷。一個(gè)載荷a量將產(chǎn)生并寫入超單元矩 陣文件!載荷步選項(xiàng)!開始求解生成超單元矩陣文件gen.sub !載荷。生成第二個(gè)載荷向量（d和m可以不變）solve!加入第二個(gè)載荷向量。!重復(fù)加載和求解生成其他的載尙向量（最多 達(dá)到31個(gè)）。finish!使用部分/c

47、lear/filname, use resume, ffull, db esel，.nsel/prep7et,，matrtx50 type, .se, geneplotfinish!清除數(shù)據(jù)庫(kù)!文件名為use!讀入整體模型數(shù)據(jù)庫(kù) !選擇非超單元部分的單元和結(jié)點(diǎn)! matrix50為超單元類型 !指向超華元類型 !讀入超單元矩陣/solu!進(jìn)入求解器antype, .!分析類型和分析選項(xiàng)!非超單元上的載荷sfe, .!施加超爭(zhēng)元載荷向量-!載荷步選項(xiàng)solve!開始求解，計(jì)算非超單元完整解和超單元凝聚解。!擴(kuò)展部分/clear/filname, gen resume, full, /solufi

48、nish!清除數(shù)據(jù)庫(kù)!文件名切換到生成部分文件名 db !讀入整體模型數(shù)據(jù)庫(kù)antype, substr expass, on expsol, . seexp, gen, usesolvefinish!進(jìn)入求解器激活擴(kuò)展選項(xiàng)!指定要擴(kuò)展的解 !要擴(kuò)展的超單元名 !載荷步選項(xiàng)（主要是輸出控制）!開始擴(kuò)展部分求解。超單元完整解寫入文件 gen.rst 或 rthrmg 屮。!奔看超單元解請(qǐng)查閱 antype, seopt, m, et, setran, se, expass 和 seexp命令得到更加詳細(xì)的說(shuō)明。超單元嵌套在ansys屮一個(gè)強(qiáng)冇力的子結(jié)構(gòu)特性就是允許嵌套：允許一個(gè)子結(jié)構(gòu)屮 含另一個(gè)

49、子結(jié)構(gòu)。也就是說(shuō)，在生成超單元時(shí)，生成部分的其中一個(gè)單元可以是以前生成的超單元。例如，有一個(gè)超單元名為piston，可以在生成一個(gè)名為cylinder的超 單元，其中包含超單元piston。為了完成柱體及其內(nèi)部的活塞的分析，就要進(jìn) 行一個(gè)使用部分計(jì)算和兩個(gè)擴(kuò)展部分的計(jì)算。使用部分計(jì)算超單元cylinder 的主自由度凝聚解,第一個(gè)擴(kuò)展部分計(jì)算cylinder的完整解和超單元piston 的凝聚解，第二個(gè)擴(kuò)展部分計(jì)算piston的完整解。有預(yù)應(yīng)力的子結(jié)構(gòu)在建模屮正確的描述系統(tǒng)特性，很重要的一點(diǎn)就是考慮其應(yīng)力情況。應(yīng)力 情況會(huì)影響剛度矩陣項(xiàng)的數(shù)值。前一步結(jié)構(gòu)分析的應(yīng)力情況在超單元生成部分形 成剛度

50、矩陣吋會(huì)被計(jì)入。應(yīng)力剛化使得一般不能承受某種載荷的結(jié)構(gòu)加強(qiáng)。比如 說(shuō)，拉緊的繩索能夠承受法a力而松弛的繩索就不能。應(yīng)力剛化同樣吋以影響同 時(shí)宥模態(tài)和瞬態(tài)動(dòng)力問(wèn)題時(shí)系統(tǒng)的響應(yīng)頻率。冇兩種不同的方法能夠生成冇預(yù)應(yīng)力的子結(jié)構(gòu)。以下用方法a和方法b列 出：方法a:1. 建立模型，定義分析類型為靜力分析(antype, main menu>solution>newanalysis),施加剛化載荷。2. 指定計(jì)算預(yù)應(yīng)力效果(pstres 或 main menu>solution>analysis options )o3. 開始靜力計(jì)算(solve 或 main menu>so

51、lution>-solve-current ls)o4. 完成生成部分。用pstres命令或其等效的gui路徑計(jì)入靜力分析 的預(yù)應(yīng)力效果。(注在靜力分析和生成部分一定要打開預(yù)應(yīng)力效來(lái)。)5. 完成使用部分和擴(kuò)展部分。6. 查看結(jié)果。方法b:注本方法不需進(jìn)行整體模型的靜力計(jì)算。1. 建模井完成生成部分。耍通過(guò)將seopt(main menu>solution> analysis options)命令選項(xiàng)sesst設(shè)為1為應(yīng)力剛度矩陣保留空間。2. 加載并完成靜力的使用部分。3. 完成擴(kuò)展部分并用 pstres(main menu>solution>analysis

52、options)計(jì)入預(yù)應(yīng)力效果。4. 保待預(yù)應(yīng)力效果打開進(jìn)行其他的生成部分以生成新的超單元。5. 求解新的生成部分丼完成使用部分。6. 完成擴(kuò)展部分并查看結(jié)果。子結(jié)構(gòu)分析實(shí)例(命令行格式) 問(wèn)題描述圓盤承受沿垂直方向的對(duì)稱載蒔，材料特性和兒何特征如下，求其應(yīng)力應(yīng)變 情況。material propertiese = 30 x 106 psi v =0.3y.fproblem sketchgeometric propertiesd =2 int =0.2 inloadingf =20001bslailjl 1l 1 arepresentative finite element model (fo

53、r pune2, pune82 and shell93 only)/prep7smrt, off/device, vector, on/title, diametral compression of a disk antype, static!靜力分析et, 1, plane2, 3, 1et, 2, plane82, 3, 1 et, 3, plane82, 3, 2 r, 1,.2mp, ex, 1, 30e6 mp, nuxy, 1,0.3 csys, 1 k, 1, 1,90 k, 2, .5,90 k, 4,1k,5, 1, 50l, 1,5lesize, 1,"7,5 l

54、, 5, 4!選擇結(jié)點(diǎn)應(yīng)力輸出 !厚度!柱坐標(biāo)系!笛卡兒坐標(biāo)系lesize, 2, 4, 2 csys, 0 k, 3l, 3,4lesize, 3, 5 l, 2,3lesize, 4, 4, 2l, 2, 5lesize, 5, 5l, 1,2lesize, 6, 7, 5a, 1, 2, 5, 5a, 2, 3, 4, 5type, 2mshk, 1msha, 0, 2damesh,2eplottype, 1mshk, 0msha, 1, 2damesh, 1eplotnsel, s, loc, x, 0nsel, r, log, y, 0eslntype, 3emodif, allsa

55、ve, vm141,dbnsel, s, log, y, 0dsym, symm, ynsel, s, loc, x, 0!沿x軸對(duì)稱dsym, symm, xnsel, allesel, all!沿y軸對(duì)稱outpr, nsol, none!結(jié)點(diǎn)位移和支反力解輸出outpr, esol, all!單元解輸出fk, 1,fy,-1ooofinish/solusolvefinish/p0st1!施加載荷值的一半（對(duì)稱）nsel, s, loc, x, 0nsel, a, log, x, 0. 1nsel, r, loc, y, 0prnsol, s, comp!選擇感興趣的結(jié)點(diǎn)nsel, r, loc, x, 0nsel, a, log, x, 0. 2nsel, r, loc, y, 0*get, snod, node, num, min*get, fnod, node, num, maxnsel, allesel, allpath, stress 1，2,，48 ppath, 1, snod ppath, 2, fnod pdef, sy, s, y prange, 24 prpath, sy*get, s1,path, 0, min, synsrl, r, loc, x, 0nsel, a, loc, x, 0. 1nsel, r, log, y, 0*get, fn