OS-T: 1000 Linear Static Analysis of a Plate with a Hole

This tutorial demonstrates the creation of finite elements on a given CAD geometry of a plate with a hole. Further, application of boundary conditions and a finite element analysis of the problem are explained. Post-processing tools are used in HyperView to determine deformation and stress characteristics of the loaded plate.

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 plate_hole.hm file you saved to your working directory from the optistruct.zip file. Refer to Access the Model Files.
  3. Click Open.
    The plate_hole.hm database is loaded into the current HyperMesh session, replacing any existing data.

Set Up the Model

Create the Material

  1. In the Model Browser, right-click and select Create > Material from the context menu.
    A default material displays in the Entity Editor.
  2. For Name, enter steel.
  3. Set Card Image to MAT1.
  4. Enter the material values next to the corresponding fields.
    1. For E (Young's Modulus), enter 2.1E+05.
    2. For NU, (Poisson's Ratio), enter 0.3.
    3. For RHO (Mass Density), leave it undefined since only a static analysis is performed.


    Figure 1. Material Property Values for steel
A new material, steel, has been created. The material uses OptiStruct's linear isotropic material model, MAT1.

Create the Property

  1. In the Model Browser, right-click and select Create > Property from the context menu.
    A default property displays in the Entity Editor.
  2. For Name, enter plate_hole.
  3. Set Card Image to PSHELL.
  4. Enter the property values next to the corresponding fields.
    An empty Value field indicates that it is turned off. To edit these properties, click on the blank Value fields next to them and enter the required values.
    1. For Material, click Unspecified > Material. In the Select Material dialog, select steel and click OK.
    2. For T (thickness of the plate), enter 10.0.


    Figure 2. Property Values for plate_hole
A new property, plate_hole, has been created as a 2D PSHELL. Material information is also linked to this property.

Update the plate_hole Component

  1. In the Model Browser, click on the component plate_hole.
    The component fields are displayed in the Entity Editor below.
  2. For Property, click Unspecified > Property. In the Select Property dialog, select plate_hole and click OK.

    os_1000_13_06
    Figure 3.
The component plate_hole has been updated with a property of the same name, and is now the current component. This component uses the plate_hole property definition with a thickness value of 10.0. The material steel is referenced by this component.

Apply Loads and Boundary Conditions

In the following steps, the model is constrained so that two opposing edges of the four external edges cannot move. The other two edges remain unconstrained. A total load of 1000N is applied at the edge of the hole in the positive z-direction.

Create Load Collectors

  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 Color and select a color from the color palette.
  4. Set Card Image to None.
    A new load collector, spcs is created.


    Figure 4. Creating the spcs Load Collector
  5. Create another load collector.
    1. For Name, enter forces.
    2. For Card Image, select None.

Create Constraints

  1. In the Model Browser, Load Collectors folder, right-click on spcs and select Make Current to set spcs as the current load collector.

    os_1000_13_08
    Figure 5. Setting spcs as the Current Load Collector
  2. From the menu bar, click BCs > Create > Constraints to open the Constraints panel.

    os_1000_13_09
    Figure 6. Accessing the Constraints Panel
  3. Make sure nodes is selected from the entity selection switch.

    os_1000_13_10
    Figure 7. Menu after Clicking on the Entity Selection Switch
  4. Hold Shift while clicking-and-dragging your mouse to select the nodes on the two ends of the plate.

    OS_1000_13_23
    Figure 8. Nodes to Select for the Constraints
  5. Constrain dof1, dof2, dof3, dof4, dof5, and dof6 and set all of them to a value of 0.0.
    • DOFs with a check will be constrained while DOFs without a check will be free.
    • DOFs 1, 2, and 3 are x, y, and z translation degrees of freedom.
    • DOFs 4, 5, and 6 are x, y, and z rotational degrees of freedom.

    os1000_13_13
    Figure 9. Constraining all Degrees of Freedom of the Selected Nodes
  6. Click create.
    Constraints are applied to the selected nodes.
  7. Click return to go back to the main menu.

Create Forces on the Nodes around the Hole

  1. In the Model Browser, set your current load collector to forces.
  2. From the menu bar, click BCs > Create > Forces to open the Forces panel.
  3. Press Shift while left-clicking, then release your mouse button to access selection options. Select Circle Interior.

    os_1000_13_13
    Figure 10. Choosing a Circular (Inside of Circle) Selection Window
  4. Hold Shift while clicking-and-dragging your mouse to select the nodes around the hole.

    os1000_interior_circle_13
    Figure 11. Nodes Selected for the Application of Loads around the Hole
  5. Define settings in the Forces panel.
    1. Set the coordinate system toggle to global system.
    2. Set the vector definition switch to constant vector.
    3. In the magnitude= field, enter 21.277 (that is 1000 divided by the number of nodes 47).
    4. Set the direction definition switch, below magnitude =, to z-axis.

    os_1000_13_15
    Figure 12. Assign Direction and Magnitude to the Forces
  6. Click create.
    Point forces, with the given magnitude in the z-direction, are applied to the selected nodes about the hole.
  7. Click return to go back to the main menu.

Create Load Steps

  1. In the Model Browser, right-click and select Create > Load Step from the context menu.
    A default load step displays in the Entity Editor.
  2. For Name, enter lateral forces.
  3. Set Analysis type to Linear Static.
  4. Define SPC.
    1. For SPC, click Unspecified > Loadcol.
    2. In the Select Loadcol dialog, select spcs and click OK.
  5. Define LOAD.
    1. For LOAD, click Unspecified > Loadcol.
    2. In the Select Loadcol dialog, select forces and click OK.

An OptiStruct subcase has been created which references the constraints in the load collector spcs and the forces in the load collector forces.



Figure 13. Creating the lateral forces Loadstep

Submit the Job

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

    OS_1000_13_17
    Figure 14. Accessing the OptiStruct Panel
  2. Click save as.
  3. In the Save As dialog, specify location to write the OptiStruct model file and enter plate_hole for filename.
    For OptiStruct input decks, .fem is the recommended extension.
  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. Set the memory options toggle to memory default.
  8. Click OptiStruct to launch the OptiStruct job.
If the job is successful, new results files should be in the directory where the plate_hole.fem was written. The plate_hole.out file is a good place to look for error messages that could help debug the input deck if any errors are present.
The default files written to the directory are:
plate_hole.html
HTML report of the analysis, providing a summary of the problem formulation and the analysis results.
plate_hole.out
OptiStruct output file containing specific information on the file setup, the setup of your optimization problem, estimates for the amount of RAM and disk space required for the run, information for each of the optimization iterations, and compute time information. Review this file for warnings and errors.
plate_hole.h3d
HyperView binary results file.
plate_hole.res
HyperMesh binary results file.
plate_hole.stat
Summary, providing CPU information for each step during analysis process.

View the Results

Displacement and Stress results for linear static analyses are output from OptiStruct by default. The following steps describe how to view those results in HyperView.

HyperView is a complete post-processing and visualization environment for finite element analysis (FEA), multibody system simulation, video and engineering data.

View a Contour Plot of Stresses

  1. From the OptiStruct panel, click HyperView.
    HyperView is launched and the results are loaded. A message window appears to inform of the successful model and result files loading into HyperView.
  2. On the Results toolbar, click resultsContour-16 to open the Contour panel.
  3. Define settings in the Contour panel.
    1. Under Result type set the first first pull-down menu to Element Stresses (2D & 3D) (t) and set the second pull-down menu to vonMises.
    2. Set the Averaging method to None.

    OS_1000_13_20
    Figure 15. The Contour panel
  4. Click Apply.
    A contoured image representing von Mises stresses should be visible. Each element in the model is assigned a legend color, indicating the von Mises stress value for that element, resulting from the applied loads and boundary conditions.
  5. In the View Controls toolbar, click the XY Top Plane View icon to change the view the model.

    OS_1000_13_21
    Figure 16. The von Mises Stress Plot for the Given Subcase
Try to answer the following questions to test your understanding of the current problem.
  • What is the maximum von Mises stress value?
  • At what location does the model have its maximum stress?
  • Does this make sense based on the boundary conditions applied to the model?

View a Contour Plot of Displacements

  1. Under Result type set the first first pull-down menu to Displacement (v) and set the second pull-down menu to Mag.
  2. Click Apply.
The resulting contours represent the displacement field resulting from the applied loads and boundary conditions.
Try to answer the following questions to test your understanding of the current problem.
  • What is the maximum Displacement value?
  • At what location does the model have its maximum displacement?
  • Does this make sense based on the boundary conditions applied to the model?

View the Deformed Shape

  1. In the View Controls toolbar, click the Isometric View icon to display the isometric view of the model.
  2. Click the Deformed toolbar icon resultsDeformed-24.
  3. Define settings in the Deformed panel.
    1. Set Result type to Displacement(v).
    2. Set Scale to Scale factor.
    3. Set Type to Uniform.
    4. For Value, enter 500.
      This means that the displacement results of the analysis is multiplied by 500.
    5. Under Undeformed shape, set Show to Wireframe.
  4. Click Apply.
A deformed plot of the model with the displacement contour should be visible, overlaid on the original undeformed mesh in isometric view.

OS_1000_13_22
Figure 17. Isometric View of the Deformed Plot Overlaid on the Undeformed Mesh (Model Unit is Set to 500)