第四章子結(jié)構(gòu)

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ī)時(shí)，并且允許在比較有限的計(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ù)的部分生成超單元，然后將它拷貝到不同的位 置，這樣做可以節(jié)省大量的機(jī)時(shí)。子結(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)自由度和三個(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)分析有以下三個(gè)步驟：生成部分

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

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

5、幾點(diǎn)：文件名在子結(jié)構(gòu)分析中很有用處。有效地使用文件名，在三部分分析中可以省略很多文件處理操作用以下方法指定文件名:Comma nd: /FILENAMEGUI: Utility Men uFileCha nge Job name女口： /FILENAME,GEN將生成過(guò)程中所有文件名都定義為 GEN。缺省的文件名是FILE （或file） 或在進(jìn)入ANSYS后定義的任意文件名。單元類型一一ANSYS提供的絕大多數(shù)單元都可以用來(lái)生成超單元。唯一的 限制是單元必須是線性的。如果生成超單元時(shí)有雙線性單元的話，ANSYS將自動(dòng)作為線性單元處理。注意：在直接耦合中帶載荷向量的耦合單元是不能做子結(jié)構(gòu)分析的

6、?？梢杂猛N形狀的單元來(lái)替代。細(xì)節(jié)參看 ANSYS Coupled-Field An alysis Guide。材料特性一一定義所有必須的材料特性。例如，如果生成質(zhì)量矩陣，就必 須定義密度或其他形式的質(zhì)量；如果要生成熱傳導(dǎo)矩陣，就要定義比熱。同樣， 超單元是線性的，非線性材料將被忽略。模型生成一一在生成部分，主要生成模型的超單元部分。非超單元部分是 在以后的使用部分生成的。但是，在建模的開(kāi)始就需要對(duì)模型的兩個(gè)部分有所規(guī)劃，主要是確定超單元部分和非超單元部分如何連接。為了保證連接正確，應(yīng)該保證接觸部分結(jié)點(diǎn)號(hào)一致。(其他可以方便用戶的方法在本章“使用部分”一節(jié) 還有介紹。)要生成整體模型應(yīng)該這樣做

7、：將模型存儲(chǔ)在數(shù)據(jù)庫(kù)文件中，選擇子結(jié)構(gòu)部分進(jìn)行生成計(jì)算。在以后的使用部分，RESUME( Utility MenuFileResumefrom )數(shù)據(jù)庫(kù)文件，不選(unselect )子結(jié)構(gòu)，用超單元矩陣代替。第二步：施加邊界條件，生成超單元矩陣。生成部分的結(jié)果包含超單元矩陣。象其他分析一樣，用戶要定義分析類型 和分析設(shè)置，施加邊界條件，定義載荷步，開(kāi)始計(jì)算。如何完成這些工作見(jiàn)下面 的敘述：1 .進(jìn)入求解器：Comma nd: /SOLUGUI: Main Me nuSolution2 .定義分析類型和分析設(shè)置：分析類型一一選擇生成超單元使用下列方法：GUI: Main Men uSoluti

8、 on-A nalysis Type-New An alysis新的分析或重啟動(dòng)一一如果是開(kāi)始一個(gè)新的分析時(shí)，只要指定分析類型（如 上所述）即可。如果是重啟動(dòng)計(jì)算，必須在ANTYPE命令中設(shè)定STATUS=REST（Main MenuSolution-AnalysisType-Restart ）。如果要另外施加載荷時(shí)，可以用重啟動(dòng)。（重啟動(dòng)時(shí)，初始運(yùn)算后的Jobname.EMAT ,Jobname.ESAV 和Jobname.DB文件要存在。）超單元矩陣文件名指定超單元矩陣文件名（Sename ）。程序?qū)⒆詣?dòng)添加后綴.SUB，因此完整的文件名是Sename.SUB。缺省是使用工作文件名/FIL

9、ENAME定義超單元矩陣文件名，可以使用以下命令：Comma nd: SEOPTGUI: Mai n Men uSolutio nAn alysis Optio ns要生成的矩陣一一可以指定僅生成剛度矩陣（或傳導(dǎo)矩陣，電磁系數(shù)矩陣）； 生成剛度和質(zhì)量矩陣（或熱傳導(dǎo)矩陣等）；生成剛度，質(zhì)量和阻尼矩陣。質(zhì)量矩 陣用于結(jié)構(gòu)動(dòng)力學(xué)分析和在使用部分有慣性載荷的情況下。在熱分析中，只有瞬態(tài)熱分析才用到熱傳導(dǎo)矩陣。對(duì)于其他分析和阻尼矩陣也大同小異。用SEOPT命令或其GUI路徑來(lái)定義。輸出矩陣這個(gè)選項(xiàng)允許輸出超單元矩陣。可以指定輸出矩陣和載荷向量，也可以只輸出載荷向量。缺省值是不輸出任何矩陣。要輸出矩陣，用

10、SEOPT 命令或其相應(yīng)的GUI路徑。7質(zhì)量矩陣形成一一只在想生成質(zhì)量矩陣時(shí)使用。用戶可以選擇缺省生成（取 決于所用單元類型）或集中質(zhì)量近似。對(duì)于絕大多數(shù)情況，推薦使用缺省生成的 方式。但是，在極薄構(gòu)件的分析中，如細(xì)長(zhǎng)桿或極薄殼體，集中質(zhì)量近似將得到 更好的結(jié)果。用下列方法指定集中質(zhì)量近似：Comma nd: LUMPMGUI: Mai n Men uSolutio nAn alysis Optio ns3 .用下列方法定義主自由度：Comma nd: MGUI: Main Men uSolutio n Master DOFsDefi ne在子結(jié)構(gòu)中，主自由度有四種作用：a.它們作為超單元與非

11、超單元的邊界。應(yīng)保證將超單元與非超單元接觸的結(jié)點(diǎn)自由度都定義為超單元（在 M命令中Lab=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 ( Ma in MeneSolutionAnalysisOptions )的使用部分需要主自由度，或者在使

12、用SETRAN命令(Mai nMenu PreprocessorCreateEleme nts-Supereleme nts-By CS Tranfer)時(shí)。在這些時(shí)候，所有主自由度的結(jié)點(diǎn)都要定義6個(gè)方向的自由度(UX，UY，UX，ROTX，ROTY，ROTZ)。4 .施加邊界條件。在生成部分可以施加所有的載荷類型，但有以下幾點(diǎn)情況需要注意：程序?qū)⑸梢粋€(gè)包括所有施加的載荷的等效載荷向量。每個(gè)載荷步一個(gè)載荷向量將寫入超單元矩陣文件中。載荷向量的最大允許值是31個(gè)。在生成部分可以使用非零的自由度約束并作為載荷向量的一部分。(在擴(kuò)展部分，如果被擴(kuò)展的載荷步中有非零的自由度約束時(shí)， 數(shù)據(jù)庫(kù)中必須有相

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

14、生成部分只能使用動(dòng)力選項(xiàng)（阻尼）。阻尼（動(dòng)力選項(xiàng)）一一只在生成阻尼矩陣時(shí)可用。指定質(zhì)量pha）阻尼使用下列方法：Comma nd: ALPHADGUI: Mai n Menu Solutio n Time/Freque ncDampi ng指定剛度（beta）阻尼使用下列方法：Comma nd: BETADGUI: Mai n Menu Solutio n Time/Freque ncDampi ng指定與材料有關(guān)的beta阻尼使用下列方法：Comma nd: MP , DAMPGUI: Mai n Menu PreprocessorMaterial Props-C on sta n-lsot

15、ropicMain Menu PreprocessorMaterial PropsPo lyno mial表4-1子結(jié)構(gòu)中可以施加的載荷Loddl Ham eLoad CategoryCom mands*Solid Model LoadsFinite Element LomIDisplacement Temperaiure Mag. Patentidil Etc.Go naira intoDK,DKLI8T, DKOELE,DL,DLLIST, DLDELE,DA, DALIST? DAOELE, DTRAND, D8VMM, DLIBTPDDELE,DSCALE, DCUMForce Hea

16、tFlaw RteIvkg, FluxEtc.FencesFK, FKLIST, FKDELEpFTRANF, FLIST, FDLE.F5CALEr FCUMPressureConveciion Maxell Surface Etc.Surface LoadsSFL, SFLLl3r SFLDELrSFA, 3FALI5, SFADEL,SFG口AD, SFTRAIMSF? 3FLISTr SFDELE? SFE, SFELISL SFEDEL, SFBEAW, SFGRAD, SFFUhl, SFSCALE, SFCUMTernperatkreHeat Generation RateCur

17、rent DensityEtc,Bo如BFK, BFKLIBLBFKDEL, BFTRANBFt BFLIST, BFDELE,BFE? BFELIS? BFEDELfBF5GALE, BFGUMGravi Linear and Angular Acc el eratio nInertia LoadsACEL, DOMfGA6 .存儲(chǔ)數(shù)據(jù)庫(kù)的備份文件。注一一這樣做的目的是需要在擴(kuò)展部分使用同樣的數(shù)據(jù)庫(kù)文件。用下列方法完成本步操作：Comma nd: SAVEGUI: Utility MenuFileSave as Jobname.db117 . 開(kāi)始計(jì)算：Comma nd: SOLVEGUI:

18、 Main MenuSolutionCurrent LS計(jì)算結(jié)果包括超單元矩陣文件，Sename.SUB , Sename是通過(guò)SEOPT指定的文件名或是工作文件名/FILENAME。矩陣文件包括根據(jù)施加的載荷計(jì)算出的載荷向量。（如果沒(méi)有施加載荷，載荷向量將為零。）8 .如有另外的載荷步，重復(fù)步驟 7來(lái)生成其他的載荷向量。載荷向量的號(hào)碼是遞增的，并添加到同一個(gè)超單元矩陣文件內(nèi)。 關(guān)于多載荷步的其他方法 見(jiàn) ANSYS Basic Analysis Procedures Guide9 . 退出 SOLUTION :Comma nd: FINISHGUI: Mai n MenuFinish使用部分

19、：使用超單元使用部分可以是任何ANSYS分析類型（FLOTRAN和顯式動(dòng)力分析除外）。和普通分析的區(qū)別就是一個(gè)或幾個(gè)單元是前面生成的超單元。每個(gè)單獨(dú)的分析指南中都有做不同分析的詳細(xì)介紹。在這一部分，我們主要介紹如何將超單元變成； 12模型的一部分。這個(gè)過(guò)程有以下幾個(gè)步驟:1 清除數(shù)據(jù)庫(kù)并指定一個(gè)新的工作文件名。2 .建立模型。3 .施加邊界條件并求解。第一步：清除數(shù)據(jù)庫(kù)并指定一個(gè)新的工作文件名。使用部分包含新的模型和新的邊界條件。因此，第一步是清除現(xiàn)存的數(shù)據(jù) 庫(kù)。這與退出并重新進(jìn)入 ANSYS的效果是一樣的。清除數(shù)據(jù)庫(kù)用下列方法：Comma nd: /CLEARGUI: Utility Me

20、nuFileClear&Start New缺省情況下，清除數(shù)據(jù)庫(kù)就會(huì)重新讀入START.ANS文件。（可以改變這個(gè)設(shè)置）注：如果通過(guò)命令行輸入來(lái)清除數(shù)據(jù)庫(kù)時(shí)，在 /CLEAR命令行中不能有其 他的命令。新定義的文件名要與生成部分使用的文件名不同。這樣，生成部分的文件 就不會(huì)被覆蓋。用下列方式之一定義新的工作文件名：Comma nd: /FILNAMEGUI: Utility Men uFileCha nge Job name第二步：建立模型本步是在PREP7中實(shí)現(xiàn)的。主要完成以下任務(wù)：1 定義MATRIX50 （超單元）為一種單元。用以下方法：Comma nd: ETGUI: Mai n Me

21、 nuPreprocessorEleme nt TypeAdd/Edit/Delete2 .定義其他非超單元的單元類型。非線性可以使用，取決于進(jìn)行分析的類型。3 .定義非超單元的單元實(shí)參和材料特性。非線性可以使用，取決于進(jìn)行分析的類型。4 .定義非超單元的幾何形狀。在定義與超單元接觸部分時(shí)應(yīng)非常注意其結(jié)點(diǎn)位置要精確重合。（見(jiàn)圖4-3 ）圖4-3 超單元與非超單元的接觸處結(jié)點(diǎn)應(yīng)與主自由度精確重合14O Nodes 1or contaci elemr Waster nodee cn sperele共有三種方法保證結(jié)點(diǎn)重合：使用與生成部分同樣的結(jié)點(diǎn)號(hào)碼。在生成部分的接觸部分結(jié)點(diǎn)和使用部分的接觸結(jié)點(diǎn)使

22、用相同的結(jié)點(diǎn)號(hào)碼增值（或平移）。（使用SETRAN，在第五步中說(shuō)明）將這兩部分結(jié)點(diǎn)所有的自由度固連起來(lái)。（用CP系列命令）在不能使用前兩種方法時(shí)，這種方法是有效的。定義自由度固連可以用下列方法：Comma nd: CPGUI: Mai n Men uPreprocessorCoupli ng/Ceq n Couple DOFs如果超單元與其他單元不相連時(shí)，在使用部分可以不定義任何單元。5 .用正確的單元類型定義超單元，讀入超單元矩陣。選擇正確的單元類型：Comma nd: TYPEGUI: Main Menu Preprocessor-Modeli ng-CreateEleme ntsElem

23、 Attributes用下列方法讀入超單元矩陣（可能要先用其他命令修改矩陣，見(jiàn)后）:Comma nd: SEGUI:Ma inMenu PreprocessorCreateEleme nts-Superleme nts-From .SUB Filea.如果模型中不包含非超單元，或是包含非超單元但接觸處的結(jié)點(diǎn)號(hào)與 主自由度結(jié)點(diǎn)號(hào)一致，那么可以用 SE命令直接讀入超單元：TYPE,！單元類型號(hào)SE, GEN!從文件GEN.SUB中讀入超單元SE命令的文件名域表示超單元矩陣的文件名。擴(kuò)展名 .SUB是給定的，因 此整個(gè)文件名就是Sename.SUB（GEN.SUB如上例）。超單元就給定下一個(gè)可用 的

24、單元號(hào)。b .如果模型中有非超單元存在，并且接觸處的結(jié)點(diǎn)號(hào)與主自由度的結(jié)點(diǎn) 號(hào)是有一個(gè)偏移量的話，就必須先用新的結(jié)點(diǎn)號(hào)生成一個(gè)新的超單元矩陣，然后 讀入新的矩陣。用下列方法生成新的超單元矩陣:Comma nd: SETRANGUI:MainMenu PreprocessorCreateEleme nts-Supereleme nts-By CS TransferMainMenu PreprocessorCreateEleme nts-Supereleme nts-ByGeom Offset用下列方法讀入新的矩陣：Comma nd: SEGUI:MainMenu PreprocessorCrea

25、teEleme nts-Supreeleme nts-From .SUB File例如，給出超單元矩陣文件 GEN.SUB和新的結(jié)點(diǎn)偏移量為2000，命令如下:SETTRAN,GEN,2000,GEN2,SUB!生成新的超單元GEN2.SUB，偏移量為2000TYPE,!單元類型號(hào)SE,GEN2!讀入新的超單元文件GEN2.SUBc.如果模型中包含非超單元而且接觸處結(jié)點(diǎn)與主自由度結(jié)點(diǎn)無(wú)任何關(guān)系17(一般是結(jié)點(diǎn)自由生成的情況)，那么先看下面的注意事項(xiàng):注意：生成部分的主自由度的結(jié)點(diǎn)號(hào)往往會(huì)覆蓋使用部分的模型結(jié)點(diǎn)號(hào)。這樣，讀入超單元將覆蓋使用部分的模型結(jié)點(diǎn)。為了避免這種覆蓋，在讀入超單元前使用結(jié)點(diǎn)

26、號(hào)偏移。在這種情況下，在使用SE命令前先存儲(chǔ)數(shù)據(jù)庫(kù)文件。因此，應(yīng)該先存儲(chǔ)數(shù)據(jù)庫(kù)文件SAVE，再用SETRAN命令用新的結(jié)點(diǎn)偏移 量生成新的超單元矩陣，最后用SE命令讀入新的矩陣。用CRINTF命令(Main Men uPreprocessorCoupli ng/Ce qn Coi ncide nt 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)偏移量為MAXNOD，文件名為 GEN2.SUBSE,GEN2NSEL,!

27、讀入新的超單元!選擇接觸處的所有結(jié)點(diǎn)CPINTF,ALL!固連所有結(jié)點(diǎn)的所有自由度NSEL,ALLd .如果要對(duì)超單元進(jìn)行移動(dòng)或拷貝到不同的位置，或進(jìn)行鏡向操作，必須 用 SETRAN 命 令 SESYMM 命 令 ( MainMenu PreprocessorCreateEleme ntsBy Reflection ）,并給出正確的結(jié)點(diǎn)偏移量，生成新的超單元矩陣文件，然后用 SE命令讀入新的矩陣。連接超單元與非超單元同樣用上面的方法一一用相同的 結(jié)點(diǎn)號(hào)，結(jié)點(diǎn)偏移量，或用 CPINTF命令。注一一如果用SETRAN命令將超單元傳遞到不同的坐標(biāo)系中， 缺省是超單 元的主自由度隨之旋轉(zhuǎn)。這一點(diǎn)在初

28、始的超單元主自由度發(fā)生旋轉(zhuǎn)時(shí)很有用， 比 如旋轉(zhuǎn)到柱坐標(biāo)系中。（在這種情況下，傳遞不影響超單元?jiǎng)偠染仃嚒＃┤绻跏?超單元沒(méi)有結(jié)點(diǎn)發(fā)生旋轉(zhuǎn)，那么傳遞的超單元也無(wú)須旋轉(zhuǎn)結(jié)點(diǎn)。用戶可以通過(guò)設(shè) 定SETRAN命令的NOROT域?yàn)?來(lái)防止結(jié)點(diǎn)旋轉(zhuǎn)。（在這種傳遞中，超單元的 剛度矩陣和載荷向量由程序自動(dòng)修改）。6 .用圖形顯示和列表功能驗(yàn)證超單元的位置是否正確。超單元用輪廓線表示，數(shù)據(jù)在生成部分寫到矩陣文件中。要生成圖形顯示，用下列命令：Comma nd: EPLOTGUI: Utility Me nuPlotEleme nts要生成列表，用下列命令：Comma nd: SELISTGUI: Utili

29、ty Men uListOtherSuperelem Data7 .存儲(chǔ)整個(gè)模型文件：Comma nd: SAVEGUI: Utility MenuFileSave as Jobname.db退出PREP7:Comma nd: FINISHGUI: Main Menu Fi nish第三步：施加邊界條件并求解。本步在求解器中完成。取得使用部分結(jié)果的過(guò)程取決于所做的分析類型。以上提到，可以在絕大多數(shù)分析中使用超單元。用戶應(yīng)當(dāng)有相應(yīng)的通過(guò)生成部分生成的矩陣。例如，如果要做結(jié)構(gòu)動(dòng)力學(xué)分析，必須有質(zhì)量矩陣。具體過(guò)程如下:1 . 進(jìn)入 SOLUTION。Comma nd: /SOLUGUI: Main

30、Me nuSolution202 .定義分析類型和分析選項(xiàng)。對(duì)于大轉(zhuǎn)動(dòng)分析，打開(kāi)大位移效果NLGEOM , ON,并正確定義非線性分 析的子步數(shù)。3 . 在非超單元上施加邊界條件。主要包括自由度約束和對(duì)稱邊界，集 中載荷，平面載荷，體載荷，和慣性載荷。注意慣性載荷只有在通過(guò)生成部分生 成了質(zhì)量矩陣時(shí)才生效。注一一對(duì)大轉(zhuǎn)動(dòng)分析，在本步中要施加正確的約束條件。4 .用以下方法施加超單元載荷向量：Comma nd: SFEGUI:Mai nMen uSolutio n-Loads-ApplyLoadVectorForSupereleme nt在超單元矩陣文件中，每個(gè)載荷步對(duì)應(yīng)一個(gè)載荷向量，用參考號(hào)來(lái)

31、區(qū)別。如，SFE，63，1，SELV,，0.75在63號(hào)單元上施加載荷向量，號(hào)為1，比例系數(shù)為0.75。因此，ELEM域 代表超單元的單元號(hào)，LKEY域代表載荷向量號(hào)（缺省為1）, Lab是SELV, VAL1 代表比例系數(shù)（缺省為0）。（參看SFE命令）21注載荷向量的方向是固定在超單元上的，因此如果超單元旋轉(zhuǎn)了，載荷向量就隨之旋轉(zhuǎn)。對(duì)于自由度方向也是如此（UX，UY，ROTY等），它們同樣 隨著超單元的旋轉(zhuǎn)而旋轉(zhuǎn)（除非 SETRAN命令中NOROT=1，此時(shí)結(jié)點(diǎn)坐標(biāo)系 不旋轉(zhuǎn)）。5 .定義與分析類型相適應(yīng)的載荷步選項(xiàng)。注用MATRIX50超單元時(shí)不要用PCG求解器。6. 開(kāi)始計(jì)算：Comm

32、a nd: SOLVEGUI: Mai n Men uSolutio n Curre nt LS本步計(jì)算包括非超單元的完整解和超單元的凝聚解一一主自由度解。非超單元的完整解記錄在結(jié)果文件中 （Jobname.RST,RTH或RMG ）,可以進(jìn)行普通 的后處理操作。凝聚解記錄在文件Jobname.DSUB中。可以通過(guò)如下方法查看這個(gè)文件：Comma nd: SEDLISTGUI: Main Men uGe neral PostprocList ResultsSuperelem DOFUtility Men uListResultsSuperelem DOF Solu如果想得到超單元中所有自由度的

33、解，就要用到擴(kuò)展部分。下面還會(huì)說(shuō)明。7 . 退出 SOLUTION。Comma nd: FINISHGUI: Main Menu Fi nish擴(kuò)展部分：在超單元中擴(kuò)展求解結(jié)果擴(kuò)展部分從使用部分的凝聚解計(jì)算出整個(gè)超單元的完整解。下面說(shuō)明擴(kuò)展部分的過(guò)程。本部分要求生成部分的.EMAT,.ESAV,.SUB,.TRI,.DB和.SELD文件， 使用部分的.DSUB文件存在。如果在使用部分中使用了結(jié)點(diǎn)偏移，在擴(kuò)展部分 中將自動(dòng)計(jì)入。1 .清除數(shù)據(jù)庫(kù)。相當(dāng)于退出并重新進(jìn)入 ANSYS。用下列方法：Comma nd: /CLEARGUI: Utility Me nuFileClear&Start New

34、2 .將文件名切換到生成部分的文件名。這樣，程序就可以識(shí)別擴(kuò)展部分所用的文件。用下列方法：GUI: Utility Men uFileCha nge Job name3 .讀入生成部分的數(shù)據(jù)庫(kù)文件。用下列方法：Comma nd: RESUMEGUI: Utility MenuFileResume Jobname.db4 . 進(jìn)入SOLUTION。用下列方法：Comma nd: /SOLUGUI: Main Me nuSolution5 .激活擴(kuò)展部分及其選項(xiàng)。擴(kuò)展部分開(kāi)關(guān)一一選為ON。Comma nd: EXPASSGUI: Mai n Me nuSolutio nExpasio nPass被

35、擴(kuò)展的超單元名一一指定 SENAME。Comma nd: SEEXPGUI: Mai n Menu Soluti on Expasi on PassExpa nd Superelem（完整的文件名假定為Sename.SUB）24使用部分生成的凝聚解文件。用SEEXP命令（或其相應(yīng)的GUI路徑）指定 該文件名。完整的文件名假定為 Usefil.DSUB。位移的實(shí)部或虛部一一只有在使用部分是諧波分析時(shí)使用。用SEEXP命令（或其相應(yīng)的GUI路徑）。被擴(kuò)展的解一一指定被擴(kuò)展的使用部分結(jié)果?？梢酝ㄟ^(guò)給出載荷步和子步，也可以通過(guò)給出時(shí)間或頻率來(lái)指定結(jié)果。用下列方法：Comma nd: EXPSOLGUI

36、: Main Menu Soluti on Expa nsion PassBy Load StepMain Menu Soluti on Expa nsion PassBy Time/Freq注一一如果被擴(kuò)展的載荷步中包含非零自由度約束時(shí)，數(shù)據(jù)庫(kù)中應(yīng)該有相 應(yīng)的自由度數(shù)值。如果沒(méi)有，那么自由度約束必須在擴(kuò)展部分中重新聲明D。6 .指定載荷步選項(xiàng)。對(duì)于子結(jié)構(gòu)擴(kuò)展部分，只有輸出控制選項(xiàng)是可用的：輸出控制一一這些選項(xiàng)控制打印輸出，數(shù)據(jù)庫(kù)和結(jié)果文件輸出和結(jié)果的插 值。如果在打印輸出文件（Jobname.OUT）中包含某些內(nèi)容，使用下列方法：GUI: Mai n Me nuSolutio n Outpu

37、t CtrlsDB/Results File如果控制結(jié)果文件(Jobname.RST)中的數(shù)據(jù)，使用下列方法：Comma nd: OUTRESGUI: Mai n Me nuSolutio n Output CtrlsDB/Results File如果想用將單元積分點(diǎn)上的數(shù)據(jù)拷貝到結(jié)點(diǎn)的方法代替插值方法(缺省) 時(shí)，用下列命令：Comma nd: ERESXGUI: Main Me nuSolutio n Output CtrlsI ntegratio n Pt7 .開(kāi)始擴(kuò)展部分：Comma nd: SOLVEGUI: Mai n Men uSolutio n Curre nt LS8 . 重

38、復(fù)步驟5到7對(duì)其他使用部分結(jié)果進(jìn)行擴(kuò)展。如果要擴(kuò)展不同超 單元的解，需要退出并重新進(jìn)入 SOLUTION。9 . 最后，退出SOLUTION。GUI: Mai n MenuFinish10 用通用方法對(duì)超單元擴(kuò)展結(jié)果進(jìn)行后處理。注一一擴(kuò)展部分不適用于功率譜密度分析。輸入示例F面給出一個(gè)子結(jié)構(gòu)分析的命令流文件示例。本例假定一個(gè)超單元與非超單元連接!生成部分!建模（超單元部分）/FILNAME，GEN/TITLE，/PREP7!文件名為GEN!進(jìn)入PREP7!生成超單元部分模型FINISH!加載并生成超單元矩陣/SOLU!進(jìn)入 SOLUTIONANTYPE, SUBST!子結(jié)構(gòu)分析SEOPT, G

39、EN,!超單元名和其他子結(jié)構(gòu)分析選項(xiàng)M，!主自由度D，!載荷。一個(gè)載荷向量將產(chǎn)生并寫入超單兀矩 陣文件!載荷步選項(xiàng)SAVE!存儲(chǔ)數(shù)據(jù)庫(kù)文件SOLVE!開(kāi)始求解一一生成 GEN.SUB文件，包含超單兀矩陣載荷向量。!載荷。生成第二個(gè)載荷向量（D和M可以不變）SOLVE!加入第二個(gè)載荷向量。!重復(fù)加載和求解生成其他的載荷向量（最多 達(dá)到31 個(gè)）。1 28 *讀萬(wàn)卷書行萬(wàn)里路FINISH!使用部分!建模/CLEAR!清除數(shù)據(jù)庫(kù)/FILNAME , USE!文件名為USE/PREP7!進(jìn)入PREP7ET, 1 , MATRIX50！ MATRIX50為超單元類型ET, 2，.!非超單兀的單兀類型!生

40、成非超單元部分模型TYPE, 1!指向超單兀類型SETRAN，!結(jié)點(diǎn)偏移量SE,!讀入由SETRAN生成的超單元EPLOT!檢杳超單元位置NSEL,!選擇接觸處的結(jié)點(diǎn)CPINTF , ALL!固連接觸處結(jié)點(diǎn)（在結(jié)點(diǎn)號(hào)不相等時(shí)使用）129 *讀萬(wàn)卷書行萬(wàn)里路NSEL, ALLFINISH!加載并求解/SOLU!進(jìn)入求解器ANTYPE，!分析類型和分析選項(xiàng)D，!非超單元上的載荷SFE,!施加超單兀載何向量!載荷步選項(xiàng)SAVE!存儲(chǔ)數(shù)據(jù)庫(kù)文件SOLVE!開(kāi)始求解，計(jì)算非超單元完整解和超單元凝聚解。FINISH!查看非超單元結(jié)果30讀萬(wàn)卷書行萬(wàn)里路!擴(kuò)展部分/CLEAR!清除數(shù)據(jù)庫(kù)/FILNAME，G

41、EN!文件名切換到生成部分文件名RESUME!讀入生成部分?jǐn)?shù)據(jù)庫(kù)/SOLU!進(jìn)入求解器EXPASS，ON!激活擴(kuò)展選項(xiàng)SEEXP, GEN，USE!要擴(kuò)展的超單元名!載荷步選項(xiàng)（主要是輸出控制）SOLVE!開(kāi)始擴(kuò)展部分求解。超單元完整解寫入文件GEN.RST 或 RTH, RMG 中。FINISH!查看超單元解請(qǐng)查閱 ANTYPE , SEOPT, M , ET, SETRAN , SE, CPINTF , EXPASS和SEEXP命令得到更加詳細(xì)的說(shuō)明。31自頂而下子結(jié)構(gòu)分析上述的子結(jié)構(gòu)方法稱為自底而上的子結(jié)構(gòu)生成方法，每個(gè)超單元是通過(guò)獨(dú)立的生成過(guò)程來(lái)生成，然后通過(guò)使用過(guò)程組裝。這種方法適用

42、于超大型的模型，分 解為小的超單元便于求解。如果對(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)建

43、立了，用戶所要做的就是選擇要生成超單元部分的單元， 施加載荷（目的是生成載 荷向量），然后 SOLVE 命令生成超單元 （Main MenuSolution-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)，并開(kāi)始使用部分的求解。4 .完成擴(kuò)展部分。首先讀入整體模型，需包含所有結(jié)點(diǎn)和單

44、元，然后對(duì)每個(gè)超單元做擴(kuò)展，每次應(yīng)定義不同的文件名，并每次進(jìn)入和退出 SOLUTION。用戶可以用通用的后處理過(guò)程查看每個(gè)超單元的結(jié)果。注一一使 用整個(gè)數(shù)據(jù)庫(kù)文件，可以讀入多個(gè)超單元結(jié)果：RESUME,FULL,DB/POST1FILE,GEN1SET,FILE,GEN2SET,!不清除上一個(gè)超單元的結(jié)果。以下給出一個(gè)自頂而下子結(jié)構(gòu)分析的輸入示例。本例假定模型有一個(gè)超單元和其他非超單元部分。!自頂而下子結(jié)構(gòu)分析的輸入示例!建立整體模型/FILNAME , FULL!文件名為FULL/TITLE ,/PREP7!進(jìn)入PREP7!生成整體模型，包括超單元部分和非超單元部分SAVE!存儲(chǔ)整體模型，以

45、備（使用部分和）擴(kuò)展部分使用FINISH!生成部分/FILNAME，GEN!文件名為GEN/SOLU!進(jìn)入 SOLUTIONANTYPE， SUBST!子結(jié)構(gòu)分析SEOPT, GEN,!子結(jié)構(gòu)分析選項(xiàng)ESEL，!選擇超單元部分的單元和結(jié)點(diǎn)NSEL34SOLVESOLVE!主自由度!載荷。一個(gè)載荷向量將產(chǎn)生并寫入超單元矩陣文件!載荷步選項(xiàng)！開(kāi)始求解一一生成超單元矩陣文件 GEN.SUB!載荷。生成第二個(gè)載荷向量（D和M可以不變）!加入第二個(gè)載荷向量。!重復(fù)加載和求解生成其他的載荷向量（最多FINISH!使用部分/CLEAR /FILNAME RESUME,!清除數(shù)據(jù)庫(kù)，USE!文件名為USEFF

46、ULL，DB !讀入整體模型數(shù)據(jù)庫(kù)35達(dá)到31 個(gè)）。ESEL,!選擇非超單元部分的單元和結(jié)點(diǎn)NSEL/PREP7ET,，MATRIX50! MATRIX50為超單元類型TYPE,!指向超單兀類型SE, GEN!讀入超單元矩陣EPLOTFINISH/SOLU!進(jìn)入求解器ANTYPE，!分析類型和分析選項(xiàng)D，!非超單元上的載荷SFE,!施加超單兀載何向量36讀萬(wàn)卷書行萬(wàn)里路SOLVEFINISH!擴(kuò)展部分/CLEAR/FILNAME , GENRESUME, FULL, DB/SOLUANTYPE, SUBSTREXPASS, ON!EXPSOL,SEEXP, GEN, USESOLVE!開(kāi)始求

47、解，計(jì)算非超單元完整解和超單元凝聚解!載荷步選項(xiàng)!清除數(shù)據(jù)庫(kù)!文件名切換到生成部分文件名!讀入整體模型數(shù)據(jù)庫(kù)!進(jìn)入求解器激活擴(kuò)展選項(xiàng)!指定要擴(kuò)展的解!要擴(kuò)展的超單元名!載荷步選項(xiàng)（主要是輸出控制）!開(kāi)始擴(kuò)展部分求解。超單元完整解寫入文件GEN.RST 或 RTH,RMG 中。A1占、讀萬(wàn)卷書行萬(wàn)里路FINISH!查看超單元解請(qǐng)查閱 ANTYPE , SEOPT, M , ET, SETRAN , SE, EXPASS和 SEEXP命 令得到更加詳細(xì)的說(shuō)明。超單元嵌套在ANSYS中一個(gè)強(qiáng)有力的子結(jié)構(gòu)特性就是允許嵌套：允許一個(gè)子結(jié)構(gòu)中 包含另一個(gè)子結(jié)構(gòu)。也就是說(shuō)，在生成超單元時(shí)，生成部分的其中一個(gè)

48、單元可以 是以前生成的超單元。例如，有一個(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)力；38情況會(huì)影響剛度矩陣項(xiàng)的數(shù)值。前一步結(jié)構(gòu)分析的應(yīng)力情況在超單元生成部分形 成剛度矩陣時(shí)會(huì)被計(jì)入。應(yīng)力剛化使得一般不能承受某

49、種載荷的結(jié)構(gòu)加強(qiáng)。比如 說(shuō)，拉緊的繩索能夠承受法向力而松弛的繩索就不能。應(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 , MainMenuSolutionNew Analysis),施加剛化載荷。2 .指定計(jì)算預(yù)應(yīng)力效果(PSTRES或 Main MenuSolutionAnalysisOptions )。3 .開(kāi)始靜力計(jì)算(SOLVE 或 Ma in Men uSolutio n-Solve-Curre ntLS)。4 .完成生成部分。用PSTRES

50、命令或其等效的GUI路徑計(jì)入靜力分析的預(yù)應(yīng)力效果。(注一一在靜力分析和生成部分一定要打開(kāi)預(yù)應(yīng)力效果。)5 .完成使用部分和擴(kuò)展部分。6.查看結(jié)果。方法B:注本方法不需進(jìn)行整體模型的靜力計(jì)算。1 建模并完成生成部分。要通過(guò)將 SEOPT(Mai n Men uSolutio nAnalysis Options)命令選項(xiàng)SESST設(shè)為1為應(yīng)力剛度矩陣保留空間。2 .加載并完成靜力的使用部分。3 . 完成擴(kuò)展部分并用 PSTRES(Ma in Me nuSolutio nA nalysisOptions)計(jì)入預(yù)應(yīng)力效果。4 .保持預(yù)應(yīng)力效果打開(kāi)進(jìn)行其他的生成部分以生成新的超單元。5 .求解新的生成部

51、分并完成使用部分。6 .完成擴(kuò)展部分并查看結(jié)果。子結(jié)構(gòu)分析實(shí)例(命令行格式)問(wèn)題描述圓盤承受沿垂直方向的對(duì)稱載荷，材料特性和幾何特征如下，求其應(yīng)力應(yīng)變40情況Material Properties E - SO x 10 psi v =0.3Geometric Prop rties d = Jin t = 0.2 inLoadingF = 20001bRepresentalive Finite Elerrierit K/lodel x (for PLATJEZ, PLANE82 and SIHELL93 only)flr/PREP7smrt,off/DEVICE,VECTOR,ON/TITLE, DIAMETRAL COMPRESSION OF A DISKANTYPE,STATIC!靜力分析ET,1,PLANE2,3,1ET,2,PLANE82,3,1ET,3,PLANE82,3,2!選擇結(jié)點(diǎn)應(yīng)力輸出R,1,.2!厚度MP,EX,1,30E6MP,NUXY,1,0.3CSYS,1!