OS-T: 1560 3-Point Bending using RBODY

This tutorial demonstrates RBODY which is used in this nonlinear large displacement implicit analysis involving contacts using OptiStruct.



Figure 1. FE Model
The following steps are included:
  • Import the model into HyperMesh
  • Set up of contacts and RBODY
  • Set up nonlinear analysis
  • View the results in HyperView

Launch HyperMesh and Set the OptiStruct User Profile

  1. Launch HyperMesh.
    The User Profile dialog opens.
  2. Select OptiStruct and click OK.
    This loads the user profile. It includes the appropriate template, macro menu, and import reader, paring down the functionality of HyperMesh to what is relevant for generating models for OptiStruct.

Open the Model

  1. Click File > Open > Model.
  2. Select the 3Point_Bending_RBODY.hm file you saved to your working directory from the optistruct.zip file. Refer to Access the Model Files.
  3. Click Open.
    The 3Point_Bending_RBODY.hm database is loaded into the current HyperMesh session, replacing any existing data.

Set Up the Model

Create the Properties

  1. In the Model Browser, right-click and select Create > Property from the context menu.
    A default PSHELL property displays in the Entity Editor.
  2. For Name, enter Shell.
  3. For Material, click Unspecified > Material.
  4. In the Select Material dialog, select Steel and click OK.


    Figure 2. Select material steel for the Property shell
  5. For T (thickness of the plate), enter 1.


    Figure 3. Property Values for Shell
    A new property, Shell, has been created as a 2D PSHELL.
  6. In the Model Browser, click on the component Fixed Cylinder1.
    The component fields display in the Entity Editor.
  7. For Property, click Unspecified > Property.
  8. In the Select Property dialog, select Shell and click OK.
    The component Fixed Cylinder1 has been updated with a property of the same name and is currently the “Current Component” (see the box in the lower right for Fixed Cylinder1). This component uses the Shell property definition with a thickness value of 1.0. The material Shell is referenced by this component.


    Figure 4. Select Property Shell for the Component Fixed Cylinder1
  9. Simarly, assign the Shell property to the other components.

Create the Sets

  1. In the Model Browser, right-click and select Create > Set from the context menu.
  2. For Name, enter Fixed Cylinder1.
  3. Right-click on the component Fixed Cylinder1 and select Isolate Only.
  4. Click on the Entity IDs to select the elements contributing to the set.


    Figure 5. Select the Elements for the Set corresponding to Fixed_Cylinder1
    The Selection panel opens.
  5. Click elems > displayed.
    This selects all the elements corresponding to the component Fixed_Cylinder1.


    Figure 6. Element Selection Panel
  6. Click proceed.
  7. Similarly, create sets named Fixed Cylinder2 and Pressing Cylinder.

Create the Rigid Bodies

In this step, you will define the components Fixed Cylinder1, Fixed Cylinder2 and Pressing Cylinder as rigid bodies.

  1. In the Model Browser, right-click and select Create > Rigid Body from the context menu.
  2. For Name, enter Fixed Cylinder1.
  3. In the Entity Editor, click on Reference Node and select the node 181607.
    Tip: Switch the model to wireframe to locate the node.
  4. For Element Sets, select Small_Cylinder1.
  5. Similarly, repeat the above steps to create rigid bodies for Fixed Cylinder2 and Pressing Cylinder.
  6. For Fixed Cylinder2, select reference node 181608.
  7. For Pressing Cylinder, select reference node 181606.

Create Set Segments

The contact surfaces will be defined, which will be used later to define the contact groups.

  1. In the Model Browser, right-click and select Create > Set Segment from the context menu.
  2. For Name, enter Fixed Cylinder1.
  3. Click Entity IDs to select the elements corresponding to Fixed Cylinder1 that are in contact with the Bumper Beam.
    Note: Make sure to switch the selection panel from faces to elements.
  4. Similarly, create a contact surface for Fixed Cylinder2, Pressing Cylinder and Bumper Beam.

Create Contact Groups

Here the contact groups will be defined.

  1. In the Model Browser, right-click and select Create > Group from the context menu.
  2. For Name, enter Fixed Cylinder1-BumperBeam.
  3. For Card Image, select CONTACT.
  4. For Property type, select SLIDE.
  5. Expand PID and select Contact surfaces.


    Figure 7. Create a Contact Group
  6. Select the contact surface created earlier.
  7. For Secondary Entity IDs, select Fixed Cylinder1.
  8. For Main Entity IDs, select Bumper Beam.
  9. Similarly create the contact group between the Fixed Cylinder2-Bumper Beam and Pressing Cylinder-Bumper Beam.

Apply Loads and Boundary Conditions

In the following steps, you will constrain the nodes 181607 and 181608 in all degrees of freedom and 181606 in all, except x direction. A force of 500N is applied on the node 181606. Other load collectors required for Nonlinear analysis are also defined.

Create SPC Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector from the context menu.
    A default load collector displays in the Entity Editor.
  2. For Name, enter spcs.
  3. Click BCs > Create > Constraints to open the Constraints panel.
  4. Select the nodes 181607 and 181608.
  5. Constrain all dof and enter a value of 0.0 to all.


    Figure 8. Constrain all Degrees of Freedom of the Selected Nodes
  6. Click create.
    This applies the constraints to the selected nodes.
  7. Similarly select the node 181606 and deselect dof1 so that a load can be applied in the x-direction.
  8. Click create.

Create Force Load Collector

This step will outline how to apply the force.

  1. In the Model Browser, right-click and select Create > Load Collector from the context menu.
  2. For Name, enter Forces.
  3. Click BCs > Create > Force to open the Force panel.
  4. Select node 181606.


    Figure 9. Select a Circular (Inside of Circle) Selection Window
  5. For coordinate system, toggle global system.
  6. Click the vector definition switch and select constant vector.
  7. For magnitude =, enter 500.
  8. Click the direction definition switch below magnitude =, and select x-axis from the pop-up menu.


    Figure 10. Assign Direction and Magnitude to the Forces
  9. Click create.
    This creates point forces, with the given magnitude in the x-direction, to be applied to the node 181606.
  10. Click return to go the main menu.

Create TABLED1 Curve

  1. In the Model Browser, right-click and select Create > Curve.
  2. For Name, enter TABLED1.
  3. For Card Image, select TABLED1 from the drop-down menu.
  4. In the Model Browser, right-click on the TABLED1 curve and select Edit enter the values:

    In the x(1) field, enter 0.0

    In the y(1) field, enter 1.0

    In the x(2) field, enter 1

    In the y(2) field, enter 1.0

  5. Click Close.
    The load collector TABLED1 that defines the time history of the loading has been created.

Create a TLOAD Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter TLOAD.
  3. For Card Image, select TLOAD1 from the drop-down list.
  4. For EXCITEID , click Unspecified > Loadcol.
  5. In the Select Loadcol dialog, select Force from the list of load collectors.
  6. For TID, select TABLED1.
  7. Click Close.

Create a DLOAD Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter DLOAD.
  3. For Card Image, select DLOAD from the drop-down list.
  4. For S (scale factor), enter 1.0.
  5. For L, select TLOAD created above.
  6. Click Close.

Create TSTEP Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter TSTEP.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select TSTEP from the drop-down menu.
  5. For TSTEP_NUM, enter 1 and press Enter.
  6. For N, enter the number of time steps as 100.
  7. For DT, enter the time increment of 0.1.
  8. For W4, enter 3e-5.
  9. Click Close.

Create NLPARM Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter NLPARM.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select NLPARM from the drop-down menu.
  5. For NINC, enter 10.
  6. For DT, enter 0.1.
  7. For CONV, select UP.
  8. For TTERM, enter 0.9.

Create NLADAPT Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter NLADAPT.
  3. For Card Image, select NLADAPT from the drop-down menu.
  4. For DTMAX, enter 0.01.

Create NLOUT Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter NLOUT.
  3. For Card Image, select NLOUT from the drop-down menu.
  4. For NINT, enter 200.
  5. Activate SVNONCNV.


    Figure 11. Create NLOUT Card

Create CNTSTB Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector from the context menu.
    A default load collector displays in the Entity Editor.
  2. For Name, enter CNTSTB.
  3. For Card Image, select CNTSTB.
  4. For S0, enter 0.01

Define Output Control Parameters

  1. From the Analysis page, select control cards.
  2. Click on GLOBAL_OUTPUT_REQUEST.
  3. Below DISPLACEMENT, ELFORCE, STRESS and STRAIN, set Option to Yes.
  4. Click return twice to go to the main menu.

Activate PARAM Control Cards

  1. From the Anaysis page, select Control Cards.
  2. For Control Cards, select PARAM.
  3. Activate HASHASSM, NLMON, and UNSYMSLV, enter YES.
  4. Activate LGDISP, enter 1.

Create Load Steps

  1. In the Model Browser, right-click and select Create > Load Step.
    A default load collector displays in the Entity Editor
  2. For Name, enter NLGEOMETRY.
  3. In the Select Loadcol dialog, select SPC from the list of load collectors and click OK.
    This selects the boundary conditions created above.
  4. In the Select Loadcol dialog, select DLOAD from the list of load collectors and click OK.
    This selects the boundary conditions created above.
  5. Similarly select the TSTEP, NLPARM, NLADAPT, and NLOUT, and assign respective load collectors.

Submit the Job

  1. From the Analysis page, click the OptiStruct panel.

    OS_1000_13_17
    Figure 12. Accessing the OptiStruct Panel
  2. Click save as.
  3. In the Save As dialog, specify location to write the OptiStruct model file and enter 3Point_Bending_RBODY.hm for filename.
  4. Click Save.
    The input file field displays the filename and location specified in the Save As dialog.
  5. Set the export options toggle to all.
  6. Set the run options toggle to analysis.
  7. Click OptiStruct to submit the job.

View the Results

  1. Once you receive the message Process completed successfully in the command window, click HyperView.
  2. Open the results and plot the displacement and the von Mises stress contour at 100% load.
  3. On the toolbar, click resultsContour-16 (Contour).
  4. Under Result type, from the first drop-down menu, select Element Stresses (2D & 3D)(t).
  5. Under Result type, from the second drop-down menu, select vonMises.


    Figure 13. Contour Panel
  6. Verify that the fields in the Contour panel match those in Figure 13 and click Apply.


    Figure 14. Displacement and Stress Result for the Analysis