ABAQUS軟件隨機(jī)振動(dòng)分析-final

1、ABAQUS軟件隨機(jī)振動(dòng)分析在工程中，結(jié)構(gòu)一般需要對(duì)它進(jìn)行隨機(jī)振動(dòng)分析。典型的例子是：通過(guò)機(jī)床的振動(dòng)響應(yīng)分析進(jìn)行機(jī)床的結(jié)構(gòu)設(shè)計(jì)，通過(guò)對(duì)結(jié)構(gòu)的地震響應(yīng)分析。在電子產(chǎn)品設(shè)計(jì)中，ABAQUS軟件不僅僅能對(duì)電子產(chǎn)品進(jìn)行沖擊、熱場(chǎng)、加工等過(guò)程進(jìn)行數(shù)值模擬，還可以對(duì)電子產(chǎn)品在隨機(jī)振動(dòng)下產(chǎn)品的響應(yīng)性能做出很好預(yù)測(cè)，以優(yōu)化產(chǎn)品設(shè)計(jì)。本例題就某電子產(chǎn)品在隨機(jī)激勵(lì)作用下的響應(yīng)結(jié)構(gòu)為例，采用如下圖所示的簡(jiǎn)化模型分析在特定隨機(jī)激勵(lì)（如圖2）中，分析該結(jié)構(gòu)的響應(yīng)。導(dǎo)入模型由于IGES文件給的是實(shí)體模型，我們?cè)谟?jì)算中產(chǎn)用shell模型，所以我們需要通過(guò)ABAQUS/CAE中對(duì)shell的編輯功能對(duì)模型進(jìn)行修改。導(dǎo)入IGE

2、S文件成Shell格式。在主菜單選擇File-Import-Part,進(jìn)入ImportPart對(duì)話框。選擇相應(yīng)的IGES文件，點(diǎn)擊按鈕。在彈出的CreatePartFromIGESFile對(duì)話框中，如下圖，對(duì)話框的Topology選擇SImportPartDirectory：_I3-16Name-RepairName二利用CAE_j.2.IGSPartname:randomFileName:FileFilter:RepairOptionsHConverttoanalyticalrepresentationHStitchEdgesIZZIConverttopreciserepresentatio

3、n編S修改模型11OWiresheii-RemoveFace,用鼠標(biāo)點(diǎn)擊選擇模型中的面，選上之后面在主菜單選擇Shape-會(huì)變紅色，點(diǎn)擊鼠標(biāo)中鍵，就可以去掉該面。重復(fù)操作，得到下圖模型。CarlZABAQUS/CAEVersion6-51一ModelDatabase:E:ftprootpartworkrandoin:二1FileModelViewportView已在匚亡ShapeFeatureToolsPlug-insHelpds申e入因口麗All自為曹ModSolidsShellExtrudeffll-昌ModelDatabase-黯Models-Mud亡.T+嗎Parts(2)+隆Mate

4、rials(1)+5feririR(4)WireCutBlendRevolveSweepLoftElanarFromSolidRemoveFace1:Model-1Part:P2-1etch模塊，如下圖所示。如上圖所示的表面，然后點(diǎn)2.在主菜單選擇Shape-Sh擊最左邊的一條邊，進(jìn)3需要增加3個(gè)外表面以及6個(gè)內(nèi)表面，利用ABAQUS的sketch模塊中的CreateLines功能，在sketch中畫(huà)出如下圖所示的線條。點(diǎn)擊圖形窗jFace，點(diǎn)擊Flipsg.y翅旳em4修補(bǔ)模型41切除表在主菜單選之后面會(huì)變紅色話框CreateParsketch模塊。在下圖的sk在主菜單選擇Shape同樣道理

5、，4.2延展表在主菜單選邊的短邊，點(diǎn)擊Type30Includetwist3pitch(Dist/Rev)Includedraft,angle(Degrees)按照上圖連接圖示的線條，點(diǎn)擊提示欄中的罰按鈕。進(jìn)入EditExtru下圖所示。EditExtrusionEndConditionUptoFacevDepthExtrudedirection:FlipOptionsNote:Twistanddraftcannotbespecifiedtogether.OK|Cane星|選擇Type選項(xiàng)為UptoFaCe,點(diǎn)擊五匚口按鈕，使得工作區(qū)的紅色箭頭指向如下圖所示。0K的示按鈕，選擇箭頭所指的面修補(bǔ)

6、好的模型如下圖所點(diǎn)擊EditEx確認(rèn)。用同樣的由于IGES文利，要將他們刪圈上畫(huà)一個(gè)圓，把原IGES文件中4.3修補(bǔ)底在主菜單選邊，點(diǎn)擊進(jìn)入Ske向負(fù)方1擊對(duì)話框中按鈕使得工作EditCutExtrusion點(diǎn)擊鼠色箭頭的方EndConditionType:Depth:ExtrudeOptionsNote:Twistanddraftcannotbespecifiedtogether.(Dist/Rev)Includetwistypitch:0(Degrees)Includedraftsangle:0按照同樣的方法把另外三個(gè)個(gè)孔和一個(gè)底座進(jìn)行修補(bǔ)。Cancel4.4分割模型。該模型中部分地方需要

7、加強(qiáng)，所用單元的殼的厚度不一樣，需要把它分割出來(lái)。在主菜單選擇Tools-Partition,進(jìn)入CreatePartition對(duì)話框，選擇type為Face,Method為Sketch，點(diǎn)擊對(duì)話框中的按鈕確認(rèn)。選擇如下圖所示區(qū)域，點(diǎn)擊提示欄中的按鈕，用鼠標(biāo)選擇區(qū)域最右的邊，點(diǎn)擊進(jìn)入Sketch模塊。Dont：利用Sketch按鈕完成。用同樣的方處理好的模區(qū)域，點(diǎn)擊鼠標(biāo)中鍵，1-1/7三添加材料屬性在Module右邊的列表框中選擇Property,進(jìn)入屬性設(shè)置模塊。因?yàn)镮GES文件中給的是米單位，所以屬性中所有的長(zhǎng)度量綱必須1.在主菜單選Elasticity選擇Mechanical-ateri

8、al-Create。進(jìn)入EditMateriaasfic。分別填入楊氏模量7e7和泊松比0.3。在EditMaterial對(duì)話框中，選擇General-Density。輸入密度為2.7e-6。點(diǎn)擊OK按鈕，退出EditMaterial對(duì)話框。在主菜單選擇Section-Create,取名為Section-l_6MM。進(jìn)入CreateSection對(duì)話框中，設(shè)置Category為Shell,點(diǎn)擊按鈕。進(jìn)入EditSection對(duì)話框,設(shè)置ShellThickness為1.6，點(diǎn)擊OK按鈕。重復(fù)4步驟，分別創(chuàng)建厚度為3mm,6mm,8mm的Section.照下圖選擇相應(yīng)的部分，按照?qǐng)D示把不同的截面

9、賦給不同的部分。點(diǎn)擊OK按鈕，完成部件屬性的分配。在主菜單選擇Assign-Section,進(jìn)入EditSectionAssignment對(duì)話框。將鼠標(biāo)按IZZIAuto-offsetfromotherinstancesCancelNote：TochangeaDependentinstancesmeshjyoumustedititspartsmesh.在Module右邊的列表框中選擇Assepbly,進(jìn)入裝配模塊。1.InstanceType:二;期星F貝口.Q貝口工(云旦口.衛(wèi).藍(lán)：.?Independent(meshoninstance)在主菜單中選擇Instance-Create，進(jìn)入C

10、reateInstance對(duì)話框中，選擇random，并在InstanceType下面選擇Dependent(meshonPart)單選鈕。如右圖：點(diǎn)擊OK按鈕，關(guān)閉CreateInstance對(duì)話框，創(chuàng)建了PartT的Instance。在主菜單，選擇Tools-Partition，在彈出的CreatePartition對(duì)話框中，選擇Cell單選鈕。點(diǎn)擊皿按鈕，退出CreatePartition對(duì)話框四定義分析步在Module右邊的列表框中選擇Step，進(jìn)入分析步模塊，建立兩個(gè)分析步驟。由于在進(jìn)行隨機(jī)響應(yīng)步驟分析之間必須要進(jìn)行一次頻率分析，所以要建立兩個(gè)分析步驟(頻率分析和隨機(jī)響應(yīng)分析）。建立

11、頻率分析步驟：主菜單選擇Step-Create，進(jìn)入CreateStep對(duì)話框，Proceduretype選擇Linearperturbation項(xiàng)，子選項(xiàng)選擇Randomresponse選項(xiàng)。點(diǎn)擊按鈕。進(jìn)入EditStep對(duì)話框，在Numberofeigenvaluesrequestedeigen右邊的子選項(xiàng)中選擇Value,并在框中輸入20。建立隨機(jī)響應(yīng)分析步驟：Step-Create,進(jìn)入CreateStep對(duì)話框，Proceduretype選擇Static,Linearperturbation項(xiàng)，子選項(xiàng)選擇Frenquecy選項(xiàng)。點(diǎn)擊按鈕。進(jìn)入EditStep對(duì)話框，在對(duì)話框的Basi

12、c頁(yè)面中，Scale選擇Linear項(xiàng)，下面的LowerFrequency，UpperFrequency，NumberofPoints，以及Bias選項(xiàng)中分別寫(xiě)入：10,2000，20，3。在Damping頁(yè)面中，選上usingdirectdampingdata，在DirectModal選項(xiàng)中依次填寫(xiě)1，20，0.001。點(diǎn)擊0K按鈕，關(guān)閉EditStep對(duì)話框。五定義相互作用屬性在主菜單中，選擇Tools-ReferencePoint.,采用默認(rèn)值（0.，0.，0.）,點(diǎn)擊鼠標(biāo)中鍵確認(rèn)。在主菜單中，選擇Constraint-Create.,進(jìn)入CreateConstraint對(duì)話框中，Typ

13、e下面選擇Coupling，點(diǎn)擊按鈕。3.在圖形窗口中，點(diǎn)擊RP-1作為約束的控制點(diǎn)，點(diǎn)擊中鍵確認(rèn)。按鈕，選擇關(guān)心的Face）取消.扌窗口下方的3.5.點(diǎn)擊圖形點(diǎn)擊對(duì)話框中的1.受默認(rèn)設(shè)置2.Name:Type:Step:Region:EditBoundaryConditionBC-1Symmetry/Antisymmetry/EncastreInitial(Picked)EditR亡呂ion.OXSYMM(Ul=UR2=QYSYMM(U2=UR1=OZSYMM(U3=UR1=OXASYMMOYASYMM心ZASYMMOPINNED(U2=(Ul二(Ul二(Ul=U3=U3二U2二U2=UR3

14、UR3UR2UR1=0)=0)=0)=0；ABAQUS/StandardABAQUS/StandardABAQUS/StandardENCASTRE(UlUR2二0;UR3二0;U3=0)=U3=UR1=UR2=UR3=only)only)only)0)4.點(diǎn)擊Condition對(duì)話框。EditBoundaryCondition對(duì)話框中的卷聖占按鈕，關(guān)閉EditBoundary七劃分網(wǎng)格在Module右邊的列表框中選擇Mesh，進(jìn)入網(wǎng)格劃分模塊。在主菜單選擇Seed-Instance在圖形窗口中，按住鼠標(biāo)左鍵，拖動(dòng)鼠標(biāo)，選中所有三個(gè)部件，點(diǎn)擊圖形窗口下方的按鈕。對(duì)話框。4.Done按鈕，完成種

15、子的設(shè)定。3.在彈出的GlobalSeeds對(duì)話框中，Approximateglobalsize右邊的文本框中輸入3.5，如下圖：點(diǎn)擊圖形窗口下的5.選擇主菜單Mesh-Controls.，鼠標(biāo)選擇如下圖所指示的四塊肋板。點(diǎn)擊鼠標(biāo)中dominated選項(xiàng)如OKTechniqueAsis(J1Free匸:,StructuredshControls對(duì)話框，選擇ElementShape1為I按鈕確認(rèn)。亡shControlsElementShapeOQuadd-dominatedOTriAlgoritlunCjMedialazciAdvancing:ront6劃分的網(wǎng)格如下圖所示：八編輯關(guān)鍵HON由于

16、CAE不支持關(guān)鍵字*PS以需要手動(dòng)加入。主采單選擇Model-Editkeywords-ModelT,進(jìn)入Editkeywords對(duì)話框。在下圖方框所在處加入*PSD-DEFINITION,如下：*pSD-DEFINITION,NAME=pSD,TYpE=BASE,G=9.81e30.04,0.0,100.04,0.0,1490.065,0.0,1500.065,0.0,3490.094,0.0,3500.094,0.0,6990.130,0.0,7000.130,0.0,13990.072,0.0,14000.072,0.0,2000MEditkeywords,Iodel:Iodel-1Ma

17、terialjname=Material-lDensityTYP：EASE,G=9.Sle3N滄-阻Elastic76+07,0.3+PSD-DEFINITION,NAME=PSDJ0.0.0,100.00.0,1490.065,0.0,1500.065,0.0,349900u-9o-9o-on-_yo34od67112uoo門(mén)：uoI.ICICLCLclclJJJJJJ4:4o022y9337u-uiIiI-uI.ICICLCLcl7o-o-2在下圖方框*BASEMOTIO*CORRELATIO*襯BOUNDARYCONDITIONS所在處加入*BASEMOTION以及CORRELATION

18、,如下：NDOF=1LOADCASE=1Syrnmetry/Antisymmetry/EncastreNpsd=psd,type=uncorrelated按鈕確認(rèn)。九定義渝出AddAfterR亡moy亡|Diwc：3rdEditwOKDiscardAllEditsCancelrlMEditkeywords,Iodel:Iodel-1Outputjfield,variable=FRESELECT*EndStep*襯STEP:Step-2StepjnameStep-SjperturbationRandomResponse10.,2000.,20,3.,1ModalDamping1,20,0.001

19、*BASEMOTION,D0Y=l3LOADCASE=1CORRELATION,PSDFSD,TYFE=UNCORRELATED1,1.*OUTPUTREQUESTSOK襯FIELDOUTPUT:F-0utput-2JI1.響應(yīng)歷史輸出在該模型中，結(jié)構(gòu)的失效往往因?yàn)楸容^關(guān)鍵的部分失效，因而，在模擬該結(jié)構(gòu)在隨機(jī)激勵(lì)載荷作用下的響應(yīng)分布的時(shí)候，需要特別關(guān)系構(gòu)件中較薄弱的部分。在構(gòu)件中比較關(guān)鍵的部分，需要定義一個(gè)輸出set,在過(guò)程中輸出這個(gè)set在激勵(lì)過(guò)程中的響應(yīng)。在Module中選擇step模塊，主菜單中選擇ToolsSet-Create，進(jìn)入CreateSet對(duì)話框。認(rèn)。采用默認(rèn)的采用默認(rèn)的na

20、me,按下name，按下按鈕，選擇卜圖所示的選擇區(qū)域，點(diǎn)擊鼠標(biāo)中鍵確Frocedtire:RandomresponseL.ontinuE.按下H-0utput-2Procedure:RandomresponseDomain:U1JJ2,Hlu_,話框中的主菜單選擇框，選擇Step；fni_itpi_iTvariaoiesTranslationsRotationsTransIntiona1andr：:：tn11ona1ve1ocitiesTranslationa1審已1ocitiesFLotational審已locit1esTranslat1ona1andr：:：tat1ona1acc丘1已ra

21、tionsTranslationa1accelerationsName:按鈕，進(jìn)入EditHistoryOutputRequest對(duì)話框。aveoutputat:Thelastincrement(:Every1(SelectfromlistbelowIFresStrainsJUR1在該對(duì)話框中MDFJR3UR,VT,2定義場(chǎng)輸出。主菜單選擇Output-Fieldoutputrequests-Create進(jìn)入CreateField對(duì)話框,選擇Step為Step2。MCreateFieldName:|F-0utput-3Step:Step-2Procedure:RandomresponseCon

22、tinue.Continue.Cancel按下出量為如下圖所示選項(xiàng)。按鈕，進(jìn)入EditFieldOutputRequest對(duì)話框，默認(rèn)區(qū)域，選擇輸MEditFieldOutputRequestName:Step:Procedure:F-0utput-3Step-2RandomresponseDomain:WholemodelvincrementsAl按鈕確認(rèn)。點(diǎn)擊對(duì)話框中Outputforrebar十提交作業(yè)atshellybeamandlayeredsectionpoints:UsedefaultsOSpecify:二邊的列表框中選J，進(jìn)入作業(yè)模塊。用When在Module右邊OKCance

23、l1在主菜單選擇Job-Create.,在彈出的CreateJob對(duì)話框中，起名為random，其他按照默認(rèn)的設(shè)置，點(diǎn)擊按鈕。Saveoutputat/:ThelastincrementEvery1OutputVariablesSelectfromlistbelow.?PreselecteddefaultsAllSubmit-random3開(kāi)始運(yùn)算。在主菜單選擇Job-Manager.,在彈出的JobManager對(duì)話框中，點(diǎn)擊按鈕，打開(kāi)如下圖所示的監(jiān)控計(jì)算過(guò)程。5.模塊，進(jìn)Module中選擇Visualization模塊，進(jìn)行結(jié)果的后處理。1.頻率分析主菜單下Result-FieldOutp

24、ut.,進(jìn)入FieldOutput對(duì)話框。進(jìn)入Historyoutput對(duì)話框，選擇set2響應(yīng)分析主菜單進(jìn)入ModeModeModeModeModeModeModeModeModeModeModeMode=q-坦-g二-Q.一坦-g二-Q.一坦-g二-Q.一坦-15HU,y16H(_：i11716.7005,1小如AZUE4.Z皿1M156.6Z：迪.4如11ZE閃亍UMl.HonaidataLoiitlitue.Increment0:BaseStateMode1:Value=1.655B3E+07Freq=647.63Mode2:Value=3.4：3702E+07Freq=933.06M

25、ode3:Value=3.72883E+07Freq=971.87Mode4:Value=4.077B2E+07Freq=1016.3Mode5:Value=4.39393E+07Freq=1055.0Mode6:Value=5.56868E+07Freq=1187.7(cycles/time.)(.cycles/time.)(cycles/time.)(cycles/time.)(.cycles/time.)(cycles/time.)Obls-XYMdatea-CTeOtej疋，入CreateOXYFData對(duì)話框，如下圖。cycl!ODEhistoryoutput(J：*ODEfield

26、outputASCIIfile【：feyboardFath(.cycles/time.)(cycles/time.)(cycles/time.)(.cycles/time.)(cycles/time.)(cycles/time.)(.cycles/time.)(cycles/time.)(cycles/time.)(.cycles/time.)(cycles/time.)(cycles/time.)選擇Source如如圖所示，點(diǎn)擊一1中節(jié)點(diǎn)的U3分量作為響應(yīng)輸出，全部選上Set1中節(jié)點(diǎn)分量U3。如下圖所示。HistoryOutputyariablesSt亡口s/Frames*OutputVar

27、iablesSpatialdisplacement:U3PI:RAKDOM-1Node5S2inNSETSET-1Spatialdisplacement:U3FI:RANDOM-1Node583inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1Node5S4inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1Node5S5inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1Node586inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1Node58

28、7inNSETSET-1Spatialdisplacement:U3PI:RAKDOM-1Node5SBinNSETSET-1Spatialdisplacement:U3FI:RANDOM-1Node589inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1Node590inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1Node591inNSETSET-1Spatialdisplacement:U3PI:RAFDOM-1Node592inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1

29、Node593inNSETSET-1oatialdis口lau亡m亡rrt:U3PI：:RANDOM-1Node594inNSETSET-1Spatialdisplacement:U3FI:RANDOM-1Node595inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1Node596inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1Node597inNSETSET-1Spatialdisplacement:U3PI:RANDOM-1Node598inNSETSET-1Spatialdisplacement:U3PI:R

30、ANDOM-1Node599inNSETSET-1Namefilter:點(diǎn)擊plot按鈕。如下兩圖，圖4為節(jié)點(diǎn)集合set1在模型所設(shè)定激勵(lì)情況下響應(yīng)的頻率分布，圖5為S11應(yīng)力分量的均方根分布云圖。同時(shí)我們還可以進(jìn)行其他結(jié)果的顯示分析，在此不贅述。讀者可自己通過(guò)查相應(yīng)的手冊(cè)來(lái)完成。Ega圖3tSetEIS,RS1L0ailBGj(fraction=(AveCelt：751)乜402eH)4-1-l.Ze5aH)4*丄0S2eH)4呂.347eH)3L79aH)3.OLOcHlSabaquspythonpostprocess.py1在ABAQUSresponse文件夾postprocess.py1件，后處理計(jì)算H的均方根mises應(yīng)附錄：pusipluwss.py乂件fromodbAccessimport*fromabaqusConstantsimport*importsys導(dǎo)入系統(tǒng)模