ANSYS中的超单元解析

1、ANSYS中的超单元摘自htbbzzg的博客，网易从8.0版开始，ANSYS中增加了超单元功能，本文通过一些实际例子，探讨了ANSYS中超单元的具体使用。1. 使用超单元进行静力分析根据ANSYS帮助文件，使用超单元的过程可以划分为三个阶段(称为Pass)：(1) 生成超单元模型(GenerationPass)(2) 使用超单元数据(UsePass)(3) 扩展模型(ExpansionPass)下面以一个例子加以说明：块板，尺寸为20x40x2,材料为钢，一端固支，另一端承受法向载荷。首先生成原始模型se_all.db即按照整个结构进行分析，以便后面与超单元结果进行比较：首先生成两个矩形,尺寸

2、各为20x2。然后定义单元类型shell6；3定义实常数1为：2板(厚度)。材料性能：弹性模量E=201000；波松比u=0.3密度P=7.8e-9单位为mm-s-N-MPa。采用边长1划分单元；端设置位移约束all另一端所有(21个)节点各承受Z向力5计算模型如下图：图1几何模型和网格圏2约鬼及载商静力分析的计算结果如下:为了后面比较的方便，分别给出两个area上的结果:图5两段结构分别的位移云图.027®.2S01&?.111193u16790-9*91HB99-31959"?USM(JIVG|MT5-0SJMX，加占幻宵5HM-.iwiiaSMX-u60436

3、翫半段WG)U3UHRSY$>0DH>：250197$i«-tawi97后半段启半段图6两段结构分别的应力云图前半段3EQVUVC6.917159,0171OS.3171-2.117125-217159.017SEOVEAVG)DJK-.250197创卅-35u142莘IK-319r2Q5超单元部分，按照上述步骤操作如下：（1）生成超单元选择后半段作为超单元，前半段作为非超单元（主单元）。按照ANSYS使用超单元的要求，超单元与非超单元部分的界面节点必须一致（重合），且最好分别的节点编号也相同，否则需要分别对各节点对建立耦合方程，操作比较麻烦。实际上，利用ANSYS中提供

4、的mesh200单元，对超单元和非超单元的界面实体，按照同一顺序，先于所有其它实体划分单元，很容易满足界面节点编号相同的要求。对于多级超单元的情况，则还要结合其它操作（如偏移节点号等）以满足这一要求。对于本例，采用另一办法，即先建立整个模型，然后再划分超单元和非超单元。即：将上述模型分别存为se_1.db（超单元部分）和se_main.db（非超单元部分）两个文件，然后分别处理。对于se_l.db模型，按照超单元方式进行处理。由于模型及边界条件已建立，只需删除前半段上的划分，结果就是超单元所需的模型。然后直接进入创建超单元矩阵的操作，首先说明一下创建超单元矩阵的一般步骤：A进入求解模块：命令：

5、/SoluGUI：Mainmenu->SolutionB设置分析类型为子结构或部件模态综合“命令：ANTYPEGUI：Mainmenu->Solution->AnalysisType->NewAnalysis选择Substructuring/CMS子（结构或部件模态综合）C设置子结构选项命令：SEOPTGUI：Mainmenu->Solution->AnalysisType->AnalysisOptions设置内容有：方程求解方法波前法墩认，或稀疏矩阵法-推荐）；要求生成的矩阵，可以是：刚度矩阵；刚度和质量矩阵；刚度、质量和阻尼矩阵稀疏矩阵法不能生成阻

6、尼矩阵）；要求输出到output窗的矩阵矩阵和载荷矢量或仅输出载荷矢量，默认不输出）;扩展超单元的方法（后续操作时，对本超单元扩展时使用的方法），默认是回代法（backsubstitutionnethod）,需要存储主自由度的三角矩阵波前法为sebname.tri稀疏矩阵法为sename.LNxx）;完全求解方法（Fullresolvenethod不存储任何三角矩阵）。随着问题规模加大，三角矩阵可能很大；完全求解方法不需要三角矩阵，但花费时间较长，且不能重启动子结构分析（对于不同的载荷矢量）。质量矩阵公式：是否使用集中质量矩阵，推荐用一致质量矩阵。D选择主自由度命令：MGUI：Mainmenu

7、->Solution->MasterDOFs->UserSelected->Define主自由度是超单元和其它单元之间的界面，必须确保将超单元与其它单元的界面上的所有节点自由度定义为主自由度。假如模型中除了超单元外没有其它单元，也必须定义主自由度。在将超单元用于动力分析时，主自由度的动力特性应能代表超单元的动力特性。对于在使用超单元时将被定义为约束或施加载荷的节点，应该定义为主自由度（如果它们也位于当前的超单元中）。对于大挠度问题或使用SETRAN命令的情况，主节点的6个自由度都应作为主自由度。本例中选择超单元外端面上所有节点的所有自由度作为主自由度。一般情况均推荐选择

8、相应主节点的Al自由度。E在超单元上施加载荷作为超单元的结构可以施加任意载荷，但要记住：程序会生成一个载荷矢量，其中包含了施加的所有载荷。对于每个载荷步，在超单元文件中写一个等效的载荷矢量，最多可以写31个载荷矢量。在超单元结构中可以使用非零位移约束，并作为载荷矢量的一部分。求解后对超单元进行结果扩展时，对于有非零位移的载荷步，其数据库必须匹配，否则需再次施加非零位移约束后再进行结果扩展。可以在使用超单元时再定义载荷和约束，此时应注意将相应节点选为主自由度节点。同样也可以在使用超单元时再定义惯性力（加速度或角速度），但在超单元部分要生成质量矩阵。如果在使用超单元时需要对其进行旋转，推荐这样做，

9、因为超单元的载荷矢量是随着超单元一起转动的。注意：如果生成质量矩阵，推荐在使用超单元时再施加约束，但要注意将约束节点选为主自由度。如果要施加加速度载荷，建议再生成超单元时施加，而不要施加到简缩的质量矩阵上。对于有大转动的问题，推荐在使用超单元时再施加约束。F设置载荷步选项在生成超单元的阶段，唯一有效的载荷步选项是设置阻尼。以上内容设置完毕，存储模型并求解。对于本例，首先选择新的求解类型：设置分析选项：rueluriitErjISwsprintedLatd¥eel*KkikkFrewLtwFront*seinerHSubstructuEiiif/CIS图8设置分靳选项ti$J1etbt

10、菇TtXffwifcftCiris)B&SLtf£Qun.tionSEKATRSEPKSESSTEKPKrEfEip临Ro髀贰莎*tiw4LlinCPJfJUxtlumpedRQ.35ppjrg?(ifn>3i±lMlrixift*be笔曲亡at芒An«iLysisTyp?irupercLcMcntm!tiefak«KaUfietst*bttwtr&ied詁"Slictf"s)?1【SEGE町JkutonAlie3-uji«iele««ht証九芒请忙*0©tenera曲口址

11、SDL耳即C4M41$£C?tSflbjtjfhxcfuyhnfRiutysiitlptlg-ASS4-IU.B4H汎t口b*4E£Ln4dto定义主自由度：然后将模型存盘并求解。求解后将生成超单元文件se_1.sub。(2) 使用超单元使用超单元的过程是：A建立新模型对于本例，将se_l.db文件存盘，点击File->Clear&StartNew菜单项，修改Jobname为se_main,Resume模型。此时读入的模型仍是个完整的模型，在前处理的Mesh下，清除后半段area上的网格，只保留前半段area上的网格以及已经施加的节点Z向载荷向载荷)。图10非

12、超单元部分的模型为了使用超单元，除了其它单元类型外，还需要定义单元类型Matrix50超单元）。图11定义超单元类型在对非超单元部分建模时，必须确保其界面节点与超单元模型的界面节点精确匹配且最好节点编号也完全一致，否则需要偏移节点编号或在超单元节点和非超单元节点之间建立耦合，或由老的超单元模型创建一个新的超单元模型(SETRAN命令），其界面节点与非超单元模型相匹配。对于本例，由于非超单元模型与超单元模型都是由同一个整体模型修改而来，故其界面节点的位置和编号都是完全相同的。B读入超单元矩阵首先设置单元属性为Superelement:命令：命令：TypeGUI:Mainmenu->Prep

13、rocessor->Modeling->Create->Elements->ElemAttributes->Type->Superelements图12定义单元属性为超单元然后读入超单元矩阵se_l.sub命令：命令：SEGUI:Mainmenu->Preprocessor->Modeling->Create->Elements->Superelements->From.SUBFile图13读入超早元矩阵读入se_l.sub文件，自动与非超单元部分组装到一起。可以通过Utilit菜单的List->Other->

14、SuperelemData列出超单元se_1的数据:EiUCfflTEH靑OFStflRELEHENT昶1FflUHCfT阳0硕RStMIftW:NUM删OF旳US值NUrtSES!(fHfilRUCES1NW«Ei!WEOCfPLATLINES-«0WICER胖DEPEESOf陂EOWPWHOOE-6WltES聊DECREEOFF煤闻1WOtCtS-曲NUMSEi!CFNODES<21闌降S!WLOffl饶tTtt第1wiEf腓nwtwntie-oBASEFILENWt>n1MUMSTlFFHESt0.1綢IE昭7TITLE-图14列表超单元数据现在应该施加超

15、单元载荷，然后将模型存盘并求解一般的静力求解）。求解后存盘，再进行后处理:图15前半段的位移去圏（与图5比较）ND£iAL3VLUTI0NSTEP"!TCHE-1USUHKTS*ODJK网MESHN-bZ4051BSHX*.flGZ436MIOSIS.3&S383"勺仞泊.filS13rUOOl.302951a412702-i.552695a*77565.8O243S图L6前半段的应力云图（与图6比较）（3）将结果扩展到超单元部分如果在生成超单元阶段产生的.EMAT、.ESAV、.SUB、.TRI、.DB、以及.SEID文件，和使用该超单元的模型求解时生成

16、的.DSUB文件都是可用的,则可以将结果扩展到超单元部分。清除当前模型，将Jobname改为超单元模型文件名(se_l),读入超单元模型，转换为Solution模块，激活ExpansionPass:命令：命令：EXPASSGUI:MainMenu->Solution->AnalysisType->ExpansionPass在EXPASS表单中，设置:ExpansionPass->ON然后读取超单元矩阵.SUB、以及非超单元模型求解时，为使用的超单元模型生成的数据文件.DSUB:命令命令:SEEXPGUI:MainMenu->Solution->LoadSte

17、pOpts->SingleExpand->ExpandSuperelements然后在SEEXP表单中分别输入超单元文件的.sub文件名和使用该超单元的主结构在求解时生成的.dsub文件名：然后进行求解->solve求解完毕，可以到后处理查看计算结果。注意此时并不生成.rs文件，但可以直接查看计算结果。对于超单元部分，相应的位移和VonMises应力云图分别见图18和图19。抽聊asolvtjon3TEP-1TIKE-JlVSVK(AVS)RSTS*DBKX-.250117STCC-.50197Du055599u166798.222397.02760&33M,I.36

18、9SS1959755019?E18超单元部分的位移分布与圏5比较）图19超单元部分的应力分布（与图6比较）KODALSDLUTiaV41丄召唱:SIS,m267.191STEPrlTIBEwi5EQ¥g如畑m«.isois?not«»S.142snix-31.20$29319G3IS.2(1乡137.1561G9.17111.14913,1£3如果存在完整的结构模型例如se_all.db）则可以显示整个结构的结果,具体操作方法是：在将Jobname改为se_1后，不是读取其本身，而是Resumese_all.dk）然后只选择超单元部分的单元和节

19、点，并按照前面相同的办法对结果进行扩展结果扩展完成后，显示整个结构的单元和节点，然后进行后处理。进入Postl后，首先执行:Data&FileOpts->ReadSingleResultFi选->择se_1.rst再点击ReadResult->First从se_1.rst中读取结果。然后执行：Data&FileOpts->ReadSingleResultFil选择>se_2.rst再点击ReadResult->First从se_2.rst中读取结果。读入这一结果不会删除从se_1.rst中读入的结果。然后可以显示或列表整个模型的结果:图20整

20、个模型的位穆云图itomlsolittidhSTEFFI普牺TIHT-lSEQVAVfi)PEDC旳2436牺耳*6.91?5HX-319.205ft,91776,3H15,711Ail.*15111.013160.41215-109204+5062勺9名。了3192。占图21整个模型的应力云图2. 使用超单元进行模态分析使用超单元进行模态分析的过程同样分为三个阶段：生成超单元模型、使用超单元数据和扩展超单元模型。使用静力分析的同一模型，约束条件与静力情况相同。首先对整个结构进行分析，相应的模态计算结果如下：NODALSOLUTIONSTEP"!趣“FP5fl-l«SlUS

21、UKN&PS¥3-0DK:56気DM31K54筑口西0"1.24,4*2S2,i393_379,3333135,76*朋皿SOLUTIONSTEM3U»-2FRE«-4$2SUSUK(AVG)P3YS-0DHXTSBu568SfK-758,$68HX宙3223169K67-335525t9.GMS19一阶模态的振型云图，fl=1051Hz口16.57133*7.141505.71267-3./厲？召勺.3S5ZS2r8562127585亨与77Sfi.SfiS图20二阶模态的振型云Eb=4525HzNOPALSULWIOSISTEMSUB-4US

22、U«(AVGI&HX-$42643HX£辭沦£4PC-：120r2SS6Z.1763.901301813斗22-号*9$43,2£4圏22四阶模态的振型云圏，U=8"83HzKmL50LUTI0MSTEP-XW耳RREQ»14734U3VHJlVR5T5-0-849r72fi0158,we56«o48475S.3136*197鼻£图23五阶模态的振型云E.f5=14734Hz图24六阶摸态的振型云图=18390Hz3皿LSOLLTIOirE.h“205.62£296J7EP-1sub-6HEq-l

23、S39Qusittimjnsrsopjft-辟字”豪SrtX筑4E41?.6139.7-1S9.SB使用超单元，按照上述步骤操作如下：（1）生成超单元操作方法与静力分析相同，先建立整个模型，然后选择后半段作为超单元，前半段作为非超单元，分别存为se_m1.db（超单元部分）和se_m2.db（非超单元部分）两个文件，然后分别处理。对于se_ml.db模型，按照超单元方式进行处理，现在要做的只需删除前半段area上的网格，结果就是建立超单元所需的模型。然后直接进入创建超单元矩阵的操作，具体步骤如下：A进入求解模块：命令：/SoluGUI：Mainmenu->SolutionB设置分析类型为

24、子结构或部件模态综合“命令：ANTYPEGUI：Mainmenu->Solution->AnalysisType->NewAnalysis选择Substructuring/CMS子（结构或部件模态综合）对于本例，首先选择新的求解类型：由于是设置超单元，故求解类型选择SUbstructuring/CMS。C设置子结构的分析选项命令：SEOPTGUI：Mainmenu->Solution->AnalysisType->AnalysisOptions设置内容有：方程求解方法波前法墩认，或稀疏矩阵法-推荐）；要求生成的矩阵：刚度和质量矩阵；要求输出到output窗的

25、矩阵默认不输出）；扩展超单元的方法（后续操作时，对本超单元扩展时使用的方法），采用默认的回代法（backsubstitutionmethod）质量矩阵公式：采用推荐的一致质量矩阵。SEfiiPT5-ubxtTuct'urariEAnkLyiaiOpt】on在HtBi«rob«vxiaE.n.-4-dWSESSI&«-Mih1+durjoeSuLtE?LWMUx«Lwp*dnaxiipproz?(ifMffa*1rjiLsgbtQritbqsi.-1Expinxtondkx白D选择主自由度命令：MGUI：Mainmenu->Solut

26、ion->MasterDOFs->UserSelected->Define选择超单元外端面上所有节点的所有自由度作为主自由度。图27定义主自由度E在超单元上施加载荷对于模态分析不需施加载荷F设置载荷步选项对于模态分析，无需载荷步设置。以上内容设置完毕，存储模型并求解，求解后将生成超单元文件se_m1.sub。(2) 使用超单元使用超单元的过程是：A建立新模型对于本例，将se_ml.db文件存盘，点击File->Clear&StartNe菜.单项，修改Jobname为se_m2，Resume模型。此时读入的模型仍是个完整的模型，在前处理的Mesh下，清除后半段上的

27、网格，得到非超单元部分的有限元模型。为了使用超单元，同样还需要定义单元类型Matrix50超（单元）。图29定义超单元类型在对非超单元部分建模时，必须确保其界面节点与超单元模型的界面节点精确匹配且最好节点编号也完全一致。对于本例，由于超单元模型与非超单元模型都是由同一个整体模型修改而来，故其界面节点的位置和编号都是完全相同的。B读入超单元矩阵首先设置单元属性为Superelement命令：TypeGUI：Mainmenu->Preprocessor->Modeling->Create->Elements->ElemAttributes->Type->S

28、uperelements然后读入超单元矩阵se_1.sub：命令：SEGUI：Mainmenu->Preprocessor->Modeling->Create->Elements->Superelements->From.SUBFile在Jobnameofmatrixfi中e输入se_m1,然后OK:读入se_ml.sub文件，自动与非超单元部分组装到一起。可以通过Utili菜单的List->Other->SuperelemData列出超单元se_m1的数据:Si强工TCflMand贋）Eili'ffiinTCtXrENT?cfSUfEL

29、ErtlH艸哑NTWTO=D商幽Sunnffii：NWBCTOFfcQUS=掏NUHBaGFlifiTUlCES=ZNUN®GFEDGEFL0711臓S0"卿GFDEC任出ffFREEUWPEE孤6NUUK®KWEE5FFRttMtlIKH®=«1NUHB&OFIVOESHNUrtEfSOFLOADMECW1跚删GF崗H利RHAUOUS.DBfiSEFILEIW1E=轴HHDU1$TIFFNE冷fl,199ME*fl7HILE图32超单元数据列表求解后存盘，再进行后处理:n闻昭souynarSTEJ-ISU34rsiOHfrsiUSUH

30、哪】P3¥3?0DIEC-5&9.gi43曲跆空4>0126.451675505,79931£BlSHZ.£1弓图33阶模态的振型云图，fi=121He图34二阶模态的振型云图，6=4528Hz駅soumoifSTE?-I.SUB临絆兰Q孙甜U34JHim)PjJYSsQDEC-759,57fi5!EC*7S9u1«7孕吕537驴R曲津b?2:23.19£50氣g糾675,170321,9*6S9CJ.?01S&uSiNAt34mriw5TEf-lSUS-3TOQ-SSSSU3U19LATC)P3¥SQDISC耳

31、目5!£C占姑心立Q0130.317S«a.$3$390.952521,2百客丄雪自imSuJlP百苫之雪,-?44S&.31142:&图35三阶模态的振型云图，E=出巧且卫图：邓四阶模态的振型云图，f4=8785HzSOLVnOM36Sr«3<63,£04SUP-1SUS-4科這Q翻毎usirai呻R3YS-QDIEC：-543.05SKX强九剜120.SOI.241.razSdu401l&l.2dz3C2.QG3图37五阶模态的振型云图，f5=15199Hz昭9QLUTIQX的吕，左目ST?-ISUS眄US-UH|AVG

32、RJYSsQDEC199.565399-.124的筑紳5T丸山*99.7*1毎空KT曹眼勺CH萌备4郃图38六阶模态的振型云图，&=：20S25Hz30UPTIQX151.613-Z3.225TEP-1515S略rao-i&issUSUHLATC)PSY3«ODHX-727.24S3KK巧纨讷43<.»1$4«.441住普粥乩驰4S就匕卫増肛表列出了整个结构和超单元方法求出的前20阶频率：No整个结构频率No超单元方法频率11051.211051.224525.324527.836551.936595.248782.648785.051473

33、4.515199.618390.620325.728307.728854.828457.831978.931867.932741.1033534.1034700.1136217.1141361.1238800.1249780.1346760.1364023.1454715.1466874.1560544.1577049.1670326.1685773.1770507.1792160.1875990.1898437.1976384.1999470.2085092.200.10180E-KJ6可以看到，超单元的频率结果比整个结构的计算结果要大，阶数越高，相差越大，即超单元的结果比整个结构要刚硬。这

34、是由于在使用超单元时，质量矩阵的处理误差比较大。(3) 将结果扩展到超单元部分如果在生成超单元阶段产生的.EMAT、.ESAV、.SUB、.TRI、.DB以及.SEID文件，和使用该超单元的模型求解时生成的.DSUB文件都是可用的，则可以将结果扩展到超单元部分。清除当前模型，将Jobname改为超单元模型文件名（se_ml），读入超单元模型，转换为Solution莫块，激活ExpansionPass:命令：EXPASSGUI:MainMenu->Solution->AnalysisType->ExpansionPass在出现的表单中:ExpansionPass->ON然

35、后读取超单元矩阵.SUB以及非超单元模型求解时，为使用的超单元模型生成的数据文件.DSUB:命令：SEEXPGUI:MainMenu->Solution->LoadStepOpts->SingleExpand->ExpandSuperelements然后在SEEXP表单中分别输入超单元文件的.sub文件名和使用该超单元的主结构在求解时生成的.dsub文件名：然后选择要扩展的模态：命令：EXPSOLGUI：MainMenu->Solution->LoadStepOpts->SingleExpand->ByLoadstep然后在EXPSOL表单中分别

36、输入要扩展的模态的载荷步和子步编号（可以在主结构中，后处理中由ResultsSummary查到），注意每次只能扩展一个模态:a!.dSingleSolutionbyLoadStep主结构求解眩，结果的载荷步和子步編号EKFSOLjE'zpandSingIe詁>1111ionbrLoadStepIdahtif7thesolutiontobe-espandedlstwmbsiwgLSTEP,SBSTEPE1e&lcCalculateele>resulis?OKCanc«lHlp图o选择要扩展的模态然后进行求解，求解完毕，可以到后处理查看计算结果。注意此时并不生

37、成.rs文件，但可以直接查看计算结果。以下是第五阶模态的形状，可以看到,与主结构结果相同，这是一个扭转模态：圏4超单元部分的第五阶模态W0DJ.1SCLUTiON470.225现37131自占勺lOS-SS1?15£.74Z23石LIZ&TEP-1SUB-1rsra-isi帕UM的创:I>MX-'POS.3373MX-leS.337如果存在完整的结构模型例如se_m_all.db）可以显示整个结构的模态形状，具体操作方法是：在将Jobname改为se_m1后，不是读取其本身，而是Resumese_m_all.db,然后只选择超单元部分的单元和节点，并按照前面相同

38、的办法对结果进行扩展。结果扩展完成后，显示整个结构的单元和节点，然后进行后处理。进入Postl后，首先执行:Data&FileOpts->ReadSingleResultFi选->择se_m1.rst;再点击ReadResult->First从se_m1.rst中读取结果。然后执行：Data&FileOpts->ReadSingleResultFi选择>se_m2.rst;再点击ReadResult->First从se_m2.rst中读取结果。读入这一结果不会删除从se_m1.rst中读入的结果。然后可以显示或列表整个模型的模态形状，注意每次

39、只能处理一个模态：图20整个模型的第5阶模态形状S9B耀&3TEP-1SQ?FREO-JSIS?USWf(AVCJRSYS-JDHX$9目卫昶.$9礼盟苹E20整个模型的第6阶模态形状mbktsolutioit141百1323直怎8Du809242.41S04.02«S£S£鼻百727.24JTEESUB<1FKEd>3带UStJttt*V&5S5T5-0m*w.胡韦SFK-737,2463. 使用超单元进行瞬态响应分析使用超单元进行瞬态响应分析的过程同样分为三个阶段：生成超单元模型、使用超单元数据和扩展超单元模型。使用静力分析的同一模

40、型，约束条件与静力情况相同。载荷为另一端各节点上的冲击力100*sin（360*10作用时间为0.05”的半正弦波）。相应的激励力时程曲线为：t(s)00.0050.010.0150.020.025F(N)030.958.7880.995.11100t(s)0.030.0350.040.0450.050.2F(N)95.1180.958.7830.900采用瑞利阻尼相当粘性阻尼E=0.03，频率取10Hz）：alpha=1.885beta=0.00048首先对整个结构进行分析，相应的响应计算结果如下：SOCIALSOLUTIDH7.1OI1吕.97©STEFFI5DB-Z6JT

41、63;q-.QZ6UZAVCRSTSbDDID-：-LS.S7SSHX-IS.9783,S51丄.1T7S呂.32S10£52£S7612.4Z7图42最大位移响应，t=0.026”BTO&ALsonrriossmp-iSUBbSIraEfl-.QSlAZ(AV)DM>：x.aiuilSFK-fcJLS7E+Q7<3B8126776*4总呻牝各,11.0E40?.132E-MJ71T4E+0"?-1S4E+07誌鶴監MT图43绝（对值）最大加速度响应，t=0.051”图44最大Mises应力响应，t二0.026”图45节点532位移响应时程曲线

42、节点12Mises应力响应时程曲线图46节点532加速度响应时程曲线VALU图47节点72Mises应力响应时程曲线使用超单元，按照上述步骤操作如下：（1）生成超单元操作方法与静力分析相同，先建立整个模型，然后选择后半段作为超单元，前半段作为非超单元，分别存为se_r1.db（超单元部分）和se_r2.db（非超单元部分）两个文件，然后分别处理。对于se_rl.db模型，按照超单元方式进行处理，现在要做的只需删除前半段area上的网格，结果就是建立超单元所需的模型。然后直接进入创建超单元矩阵的操作，具体步骤如下：A进入求解模块：命令：/SoluGUI：Mainmenu->Solution

43、B设置分析类型为子结构或部件模态综合“命令：ANTYPEGUI：Mainmenu->Solution->AnalysisType->NewAnalysis选择Substructuring/CMS子（结构或部件模态综合）对于本例，首先选择新的求解类型：由于是设置超单元，故求解类型选择SUbstructuring/CMS。图48选择求解类型C设置子结构的分析选项命令：SEOPTGUI：Mainmenu->Solution->AnalysisType->AnalysisOptions设置内容有：方程求解方法波前法墩认，或稀疏矩阵法-推荐）；要求生成的矩阵：刚度矩阵

44、，质量矩阵和阻尼矩阵；要求输出到output窗的矩阵默认不输出）；扩展超单元的方法（后续操作时，对本超单元扩展时使用的方法），采用默认的回代法（backsubstitutionmethod）质量矩阵公式：采用推荐的一致质量矩阵。图49SEOPT命令表单，设置分析选项D选择主自由度命令：MGUI：Mainmenu->Solution->MasterDOFs->UserSelected->Define选择超单元外端面上所有节点的所有自由度作为主自由度。=宣墅竖.置图50定义主自由度以上内容设置完毕，存储模型并求解，求解后将生成超单元文件se_r1.sub(2)使用超单元使用

45、超单元的过程是：A建立新模型对于本例，将se_rl.db文件存盘，点击File->Clear&StartNew菜单项，修改Jobname为se_r2，Resume模型。此时读入的模型仍是个完整的模型，在前处理的Mesh下，清除后半段上的网格，得到非超单元部分的有限元模型。图51非超单元部分的模型为了使用超单元，同样还需要定义单元类型Matrix50超（单元）。在对非超单元部分建模时，必须确保其界面节点与超单元模型的界面节点精确匹配且最好节点编号也完全一致。对于本例，由于超单元模型与非超单元模型都是由同一个整体模型修改而来，故其界面节点的位置和编号都是完全相同的。B读入超单元矩阵首

46、先设置单元属性为Superelement：命令：TypeGUI：Mainmenu->Preprocessor->Modeling->Create->Elements->ElemAttributes->Type->Superelements图52定义单元属性为超单元然后读入超单元矩阵se_1.sub：命令：SEGUI：Mainmenu->Preprocessor->Modeling->Create->Elements->Superelements->From.SUBFile在Jobnameofmatrixfi中输入se

47、_m1，然后OK:AddLtdlteadinSnpMfleMrtPfob.EtHixPilbGSBB”4fUtt卵比卞皿日亡工«Ad-StiMillsfsIvinfuertLcflh<-ntFroiihtzixFbl«(.JobtLimtofA*.trixfiITomTclfF*ngtZ图53读入超单元矩阵se_r1.sub读入se_rl.sub文件，自动与非超单元部分组装到一起。可以通过Utili菜单的List->Other->SuperelemDaJO出超单元se_r1的数据:图54超单元数据列表C施加瞬态激励在模型端部的节点上施加瞬态激励，与整个结构

48、的情况相同。图55施加瞬态激励力求解后存盘，再进行后处理，相应的响应计算结果如下:EODJLLSDLUTiaiT?,2fia022Ba505LDBgfie13-471FREq«.GZ6TIE(A冏R5ira«oUK1乩虻4not-4,761SBC-15»S4图56最大位移响应，t=0.02”6ITOILSOLUTIONSTEP"!SUB»51FFZQ日051JlZ(AVC)DKX-.Z411373HK-.i«£4!D70苗JL目幹4的臨陀-127E+O7159112477536TOSSfiflK111E+D7.143E+07图57绝对值）最大加速度响应，t=0.051节点幻2位移晌应时程曲线TIME图58节点532位移响应时程曲线15QQ12501000hTIME节点532加速度响应时程曲线(x10*3)VALU250*250-5