Download EN175 ABAQUS tutorial and more Summaries Engineering in PDF only on Docsity!
ENGN 1750: Advanced Mechanics of Solids
ABAQUS TUTORIAL
School of Engineering Brown University
This tutorial will take you all the steps required to set up and run a basic simulation using ABAQUS/CAE and visualize the results;
Background
1. The figure shows an FEA simulation of a rigid sphere rebounding of a soft rubber thin film. - The sphere is idealized as a near-incompressible neo-Hookean hyperelastic material (we will discuss what this means later in the class â for now, all you need to know is that this material model is often
used to model elastomers). The solid has shear modulus ÎŒ , bulk modulus K and mass density Ï.
We will focus on the limit where K >> Ό.
- The sphere has radius R, mass m , and is launched with initial velocity V 0 towards the rigid surface.
In this tutorial you will run a finite element simulation of this problem.
- Open ABAQUS/CAE (you will find a link on the Start menu). From the âStart sessionâ menu select âCreate Model Databaseâ with Standard/Explicit Model
- To set up and run a simulation in CAE, you can either work through the various âModulesâ in the dropdown menu at the top of the screen, or you can click on the menus in the workflow tree on the left hand side of the screen. We will do the former in this tutorial.
- The âPartâ module is used to define the geometry of the various interacting solids in the simulation. Here, we will create the elastomeric film and the rigid cylinder as two separate parts. Start by creating the sphere: a. Select Part>Create from the top menu bar, then use the menu that comes up to name the part âSphere;â select âAxisymmetricâ and âDiscrete Rigidâ for the Modeling Space and Type, and enter an Approximate Size of 0.75. (you can try âAnalytical Rigidâ as well if you like, but ABAQUS will sometimes throw an error when you try to make a semicircle with this option) b. Use the sketch tools from the left to draw a semicircle, centered at (0,0.05), with a radius of 0. c. To exit the sketch menu, click the circle sketch tool to unselect it (or you can click the sketch window with both mouse buttons and select âcancel procedureâ), then select âDoneâ from the menu below.
Rigid surfaces need a âreference pointâ that can be used to define the motion of the part, or forces acting on the part. To add a reference point to the cylinder Select âTools>Reference Pointâ in the Part module, and select the point at the center of the sphere.
Next, create the film: a. Select Part>Create; name the part âFilmâ, and select âAxisymmetricâ and âDeformable;â b. Draw a rectangle â the exact dimensions are not important, but place the top left corner of the rectangle at the origin. Exit the sketch menu.
If you need to make any changes to the parts, you can do so using the Model Database tree on the left- you can expand the various options and edit them as needed.
- Next, we will define material properties, and assign them to the parts. a. Switch the âModuleâ to âPropertyâ b. To create a hyperelastic material, select Material>Create c. Change the name of the material to âHyperelasticâ. To define properties, Select âGeneralâ on the popup window, and enter a density value of 1. d. Next select the Mechanical>Elastic>Hyperelastic. In the menu that comes up select âNeo Hookeâ for the strain energy potential, click the âCoefficientsâ radio button, and enter a value of 1.0 for C10. Leave the D field blank â if you do this, ABAQUS will assign an internal value that makes the material approximately incompressible. e. Next, you have to create a âsection.â This makes more sense for parts like beams or truss elements, where you need to define the geometry of the cross section, but you have to do it even for a 3D part. Select Section> Create, then select Homogeneous Solid section from the menu, set the options in the Edit Section menu as shown, then check that press OK from the menu shown f. Switch the âPartâ to Film, then select Assign>Section. Click on the rectangular film to select it (it should turn pink), then select âDoneâ from the bottom of the window. This has now assigned the hyperelastic material model to the part representing the film. If you want to make a part heterogeneous, you can partition it, and assign different properties to the different partition. g. Next, create the inertial properties of the rigid sphere. Switch the âPartâ to Sphere, then select Special> Inertial Properties.
Sphere-Ref-Point; select Geometry, and press Continue. Then, select the reference point on the sphere, and press Done. f. To create the output request use Output>History Output Requests> Create. Name the history output Sphere Motion and press continue. g. Then, in the Edit History Output Request menu, for the Domain select Set; and select Sphere-Ref-Point for the set. For Frequency, select âEvery n time incrementsâ and enter 1 for n. Check the âDisplacement/Velocity/Accelerationâ box, and then edit the menu to show U2,V2,A2 (the vertical displacement, velocity and acceleration). Finally, press OK.
- Next, we set up the properties of the contact between the sphere and the film. Select the âInteractionâ module to start a. First, define the properties of the contact. Select Interaction>Property>Create b. Name the property Hard-frictionless and select the Contact option, then press continue c. Select the Mechanical > Normal Behavior option and accept the defaults; d. Select the Mechanical > Tangential Behavior; check that the friction formulation says Frictionless, then press OK e. Next, specify the two surfaces that will contact: select Interaction>Create f. Name the interaction Sphere-Film-Interaction and apply it to the Impact step g. Select the Surface-to-surface contact option; press Continue h. Click the outline of the sphere in the window to select it, press Done, then select color of the arrow pointing towards the outer surface of the sphere (ABAQUS needs to know whether the contact will occur from inside the sphere or the outside) i. For the second surface, select the Surface option, then click the surface of the film, then press Done j. In the âEdit Interaction menu, check that âFinite Slidingâ is selected, the interaction property should be Hard-Frictionless. You can accept defaults for the other options. Press OK. You should see the interacting surfaces marked with small rectangles in the window Now that the interaction has been created, we can configure CAE to save results associated with the contact to the output database. To do this, return to the âstepâ module (i.e. part 7 above), then k. select Output>History Output Requests>Create l. Name the history request Contact-Force, and apply it to the Impact step m. Select âInteractionâ for the domain dropdown menu; select âEvery n time incrementsâ and set n=5 for the Frequency, and then activate the total forces due to contact pressure and the total contact area (see the figure)
- Next, we specify boundary conditions and initial conditions for the analysis. Switch to the âLoadâ module. a. To enforce the axial symmetry boundary condition for the film, select BC>Create; then in the Create Boundary Condition menu, name the boundary condition Film-Axial-Sym; select the initial step; the Mechanical category, and the Symmetry/Antisymmetry/Encastre type; then press Continue. b. Select the leftmost side of the rectangular film c. Select the XSYMM option, and press OK (note that the UR2 and UR3 degrees of freedom are not active for the film, even though they show up in the menu) d. To constrain the vertical velocity of the base of the film, select BC>Create; in the Create Boundary Condition menu name the boundary condition Film-Base; select the Initial step; select Displacement/Rotation, and press Continue e. Select the base of the film and press Done f. In the Edit Boundary Condition menu select U2 and press OK g. To enforce the axial symmetry boundary condition for the sphere, select BC>Create h. Name the boundary condition Sphere-Axial_Sym; select the Initial step, the Mechanical category and Symmetry/Antisymmetry/Encastre and press Continue i. Select the reference point (RP) on the sphere and press Done j. Select XSYMM on the Edit Boundary Condition menu and press OK k. To specify the initial downward velocity of the sphere, Predefined Field>Create l. Name the field Initial-Vel-Sphere; select the Initial step, the Mechanical category and the Velocity type. Press Continue m. Select the reference point (RP) on the sphere n. Enter -0.05 in the V2 box on the Edit Boundary Condition menu, and leave the V1 box blank. Press OK.
- ABAQUS also creates a large number of files during an analysis. You can find them in the ABAQUS working directory. The following files can be viewed using a text editor and may contain useful information (eg error messages if an analysis fails). Itâs worth checking each to see what they contain: a. Impact-Simulation.inp is an ASCII file that defines the finite element simulation. This is how ABAQUS/CAE (which is just a GUI) communicates with the ABAQUS finite element platform. We often edit this file and run ABAQUS with it directly when we use ABAQUS with user- element or user-material subroutines. See the last section of this tutorial for more details. b. Impact-Simulation.log summarizes the outcome of each phase of the analysis c. Impact-simulation.sta is updated during analysis to show its progress d. Impact-simulation.msg contains error messages e. Impact-simulation.dat contains both error messages and the output of an analysis. It is sometimes the only way to easily view results of an analysis that has user elements.
- Once the job was completed, you can view the results. To do so, start by selecting the Visualization module. a. You will need to load the âoutput databaseâ by using File>Open⊠then navigate to the ABAQUS working directory (if needed), and open the file called Impact-Simulation.odb b. You can now explore the results of the simulation. For example, to plot von Mises stress contours on the deformed mesh, select c. You can now use the buttons to scroll through the history of the simulation. d. You can view an animation of the impact by selecting Animate > Time History. You can change the behavior of the animation by selecting the Options > Animation menu. e. If you would like to plot contours of variables other than the Von Mises stress, go to Result>Field Output. The Field Output menu will allow you to select different variables. If you donât see the variable you want to plot that usually means ABAQUS did not save it to the odb. To get it, you will need to go back to the Step menu, and use the Output> menu to add the variable you are interested in. The analysis will need to be re-run. f. You can change the way the contour plots look by using the Options>Common menu, or Options>Contour. g. Some additional features of the viewport can be controlled using the Viewport>Viewport Annotation Options. h. You often want to change the background of the viewport to save plots for reports â you can do this using the View>Graphics Options menu i. ABAQUS can also plot graphs showing spatial or temporal variation of quantities of interest. For example, to plot the contact force acting between the sphere and the film as a function of time, you can use Result>History Output⊠and then select Total force due to contact pressure: CFN2 ASSEMBLY_M_SURF-1/ASSEMBLY_S_SURF-1 (you can select either of the CFN variables â one is the force on the sphere; the other is the force on the film. j. ABAQUS/CAE also allows you to post-process data. For example, you can save the contact force history to a dataset using Tools>XY Data>Create; then select Odb history from the XY Data menu. To change the data, you can use Tools>XY Data> Create, then select Operate on XY data. The menu that opens allows you to enter formulas that create new XY datasets by modifying existing ones. In practice itâs usually easier to post-process odb data by writing a python script that will read and modify data extracted from an odb, however. k. Similarly, you can plot the velocity or displacement of the center of the sphere, by using Result>History Output, and then select Spatial velocity: V2 Pl: SPHERE-1 Node nn in NSET SPHERE-REF-POINTâ; then press plot l. The visualization module has many other features; if you have time, the best way to learn them is just to explore all the various options and see what happens.
