Abaqus/CAE(ver.6.
10)NonlinearBucklingTutorial
ProblemDescription
This is the NAFEMS1 proposed benchmark (Lees frame buckling) problem. The applied load is based on the normalized (EI/L2) value of F = 996.389 N. The analysis will investigate post-buckling nonlinear behavior of the frame at the applied load location. This tutorial will also describe x-y plotting capability in Abaqus/CAE, including combining variables to generate load-displacement plots. E = 71.74109 N/m2 = 0.0 L = 1.2 m
0.2 L
0.8 L
L
0.02
0.03
National Agency for Finite Element Methods and Standards, NAFEMS Non-Linear Benchmarks (Glasgow: NAFEMS, Oct., 1989., Rev. 1.) Test No. NL7.
2011 Hormoz Zareh & Donald Bell
Portland State University, Mechanical Engineering
AnalysisSteps
1. StartAbaqusandchoosetocreateanewmodeldatabase 2. InthemodeltreedoubleclickonthePartsnode(orrightclickonpartsandselectCreate)
3. IntheCreatePartdialogbox(shownabove)namethepartand a. Select2DPlanar b. SelectDeformable c. SelectWire d. Setapproximatesize=10 e. ClickContinue 4. Createthegeometryshownbelow(notdiscussedhere)
2011 Hormoz Zareh & Donald Bell
Portland State University, Mechanical Engineering
5. DoubleclickontheMaterialsnodeinthemodeltree
a. Namethenewmaterialandgiveitadescription b. ClickontheMechanicaltab Elasticity Elastic i. EnteraYoungsmodulusof71740000000,andPoissonsratioof0 c. ClickOK
2011 Hormoz Zareh & Donald Bell
Portland State University, Mechanical Engineering
6. DoubleclickontheProfilesnodeinthemodeltree a. NametheprofileandselectRectangular b. ClickContinue c. Enter0.03foraand0.02forb d. ClickOK
7. DoubleclickontheSectionsnodeinthemodeltree a. NamethesectionbeamandselectBeamforthecategoryandBeamforthetype b. ClickContinue c. Selecttheprofilecreatedabove(rect_beam)andthematerialcreatedabove(Material1)
2011 Hormoz Zareh & Donald Bell
Portland State University, Mechanical Engineering
d. ClickOK
8. ExpandthePartsnodeinthemodeltree,expandthenodeofthepartjustcreated,anddoubleclickon SectionAssignments a. SelecttheentiregeometryintheviewportandpressDoneinthepromptarea b. Selectthesectioncreatedabove(beam) c. ClickOK
9. InthetoolboxareaclickontheAssignBeamOrientationbutton a. Selectallthegeometry b. ClickDone
2011 Hormoz Zareh & Donald Bell
Portland State University, Mechanical Engineering
c. d. e. f.
Leavethedefaultvaluesof0.0,0.0,1.0 PresstheEnterkey Thebeamnormalsshouldbeorientedasshownbelow. ClickOKtoconfirm
10. ExpandtheAssemblynodeinthemodeltreeandthendoubleclickonInstances a. SelectDependentfortheinstancetype b. ClickOK
2011 Hormoz Zareh & Donald Bell
Portland State University, Mechanical Engineering
11. DoubleclickontheStepsnodeinthemodeltree a. Namethestep,settheproceduretoGeneral,andselectStatic,Riks,ClickContinue b. OntheBasictab i. Givethestepadescriptionand ii. Setgeometricnonlinearitieson(Nlgeom=ON) iii. UnderStoppingcriteriacheckMaximumloadproportionalityfactorandsetto30
c. OntheIncrementationtab, i. Settheinitialarclengthincrementsizeto0.1 ii. Setthemaximumarclengthincrementsizeto2 iii. Setthemaximumnumberofincrementsto200
2011 Hormoz Zareh & Donald Bell
Portland State University, Mechanical Engineering
d. ClickOK 12. DoubleclickontheBCsnodeinthemodeltree a. NametheboundaryconditionedPinnedandselectDisplacement/Rotationforthetype b. ClickContinue
c. SelectthetwofreeendsoftheframeandclickDone i. Note:toselectmultipleitems,holdtheshiftkey d. SelectU1andU2andsettozero,clickOK
2011 Hormoz Zareh & Donald Bell
Portland State University, Mechanical Engineering
13. DoubleclickontheLoadsnodeinthemodeltree a. NametheloadandselectConcentratedforceasthetype b. ClickContinue
c. Selectthepointalongthetopbeamnearthecorner,ClickDone d. SetCF1to0andCF2to996.389 e. ClickOK
2011 Hormoz Zareh & Donald Bell
Portland State University, Mechanical Engineering
14. InthemodeltreedoubleclickonMeshfortheframepart,andinthetoolboxareaclickontheAssign ElementTypeicon a. Selecttheentiregeometry b. SelectStandardforelementtype c. SelectLinearforgeometricorder d. SelectBeamforfamily e. Notethatthenameoftheelement(B21)anditsdescriptionaregivenbelowtheelementcontrols f. ClickOK
2011 Hormoz Zareh & Donald Bell
10
Portland State University, Mechanical Engineering
15. InthetoolboxareaclickontheSeedParticon
a. Enter0.08forApproximateglobalsize,clickOK
2011 Hormoz Zareh & Donald Bell
11
Portland State University, Mechanical Engineering
16. InthetoolboxareaclickontheMeshParticon,ClickYes
17. Createasetfortheuppercentervertex a. ExpandtheAssemblynodeinthemodeltree,andthendoubleclickonsets b. Namethesetandselectgeometryforthetype,ClickContinue c. Selectthevertexwheretheloadisapplied,ClickDone
18. ExpandtheHistoryOutputRequestsnodeinthemodeltree,andthenrightclickonHOutput1(HOutput1 wasautomaticallygeneratedwhencreatingthestep)andselectDelete
19. DoubleclickontheHistoryOutputRequestsnode a. NamethehistoryandselectContinue b. SetthedomaintoSetsandselectthesetcreatedabove
2011 Hormoz Zareh & Donald Bell
12
Portland State University, Mechanical Engineering
c. Leavethefrequencysettoeveryincrement(n=1) d. FortheoutputvariablesselecttheU2displacement
20. InthemodeltreedoubleclickontheJobnode a. Namethejobframe_buckle b. Givethejobadescription
2011 Hormoz Zareh & Donald Bell
13
Portland State University, Mechanical Engineering
21. InthemodeltreerightclickonthejobjustcreatedandselectSubmit a. WhileAbaqusissolvingtheproblemrightclickonthejobsubmitted,andselectMonitor
b. IntheMonitorwindowcheckthattherearenoerrorsorwarnings i. Ifthereareerrors,investigatethecause(s)beforeresolving ii. Iftherearewarnings,determineifthewarningsarerelevant,somewarningscanbesafely ignored iii. Inthefarrightcolumn,notehowAbaqusadjustedtheincrement
22. Inthemodeltreerightclickonthesubmittedandsuccessfullycompletedjob,andselectResults
2011 Hormoz Zareh & Donald Bell
14
Portland State University, Mechanical Engineering
23. InthemenubarclickonViewport ViewportAnnotationsOptions a. UnchecktheShowcompassoption b. ThelocationsofviewportitemscanbespecifiedonthecorrespondingtabintheViewport AnnotationsOptions
24. Displaythedeformedcontourofthe(Von)Misesstress a. Inthetoolboxareaclickonthefollowingicons i. PlotContoursonDeformedShape b. Notethatwhenincludingtheeffectsofgeometricnonlinearities,thedeformationscalefactor defaultstoavalueof1
2011 Hormoz Zareh & Donald Bell
15
Portland State University, Mechanical Engineering
25. Clickonthearrowsonthecontextbartochangethetimestepbeingdisplayed a. Clickonthethreesquarestobringuptheframeselectorsliderbar
26. Ontheresultstree,expandtheHistoryOutputnodeanddoubleclickonthedisplacementhistorycreated
a. NoticethatdisplacementitplottedagainstArcLength,notLoadorLoadProportionalityFactor. b. Toplotloadagainstdisplacement,wewillneedtoextractthevaluesforLoadanddisplacementfrom theFieldOutputs.
2011 Hormoz Zareh & Donald Bell
16
Portland State University, Mechanical Engineering
27. IntheToolboxareaclickontheCreateXYDataicon a. ChooseODBfieldoutputforSourceandclickContinue b. OntheVariablestab i. SelectUniqueNodalforPosition ii. ExpandCF:PointloadsandselectCF2 iii. ExpandU:SpatialdisplacementandselectU2
c. SelecttheElements/Nodestab i. SelectNodeSetsforMethod ii. SelectthesetcreatedearlierTop
2011 Hormoz Zareh & Donald Bell
17
Portland State University, Mechanical Engineering
d. ClickSave,thenOKonthenextwindow e. ClickDismissontheXYDatafromODBFieldOutputwindow 28. ExpandtheXYDatanodeontheresultstree. a. Thereshouldnowbetwosetsofdataunderthenodeasshown.
b. DoubleclicktheXYDatanode c. ForSourceselectOperateonXYdata d. FromtheOperatorslistselectcombine(X,X),Itshouldappearintheexpressionboxatthetopof thewindow. i. Thecombine(X,X)operatorcombinestwosetsofsavedXYdata ii. TheYvaluesofthefirstargumentbecometheXvaluesofthenewXYdata iii. TheYvaluesofthesecondargumentbecometheYvaluesofthenewXYdata iv. ThevaluesarecombinedwherevertheXvaluesofthetwoargumentsalign v. FormoredetailseeAbaqus/CAEUsersManualsection45.4.39,CombiningtwoXYdata objects
2011 Hormoz Zareh & Donald Bell 18 Portland State University, Mechanical Engineering
e. SelectU:U2P1:PART.FromtheXYDatasectionandclickAddtoExpression f. SelectCF:CF2PI:PART.FromtheXYDatasectionandclickAddtoExpression g. Sincetheloadanddisplacementbothincreaseinthenegativedirection,theyneedtobemultiplied by1tomakeloadanddisplacementincreaseinthepositivedirection. h. Thefinalexpressionshouldlooklike:
i. j.
ClickSaveAs,nameitloaddisplacement,clickOK ClosetheOperateonXYDatawindow
29. RightclickonloaddisplacementundertheXYDatanodeandselect Plot a. Thebucklingbehaviorcanbeseenintheplot.
2011 Hormoz Zareh & Donald Bell
19
Portland State University, Mechanical Engineering
30. ThisdatacanalsobecopiedintoExcelorotherprograms. a. RightclickonloaddisplacementundertheXYDatanodeandselectEdit b. Selectallthedataintheeditwindow,rightclickandchooseCopy c. OpenExcel,rightclickinanemptycellandchoosePaste
2011 Hormoz Zareh & Donald Bell
20
Portland State University, Mechanical Engineering