OS-T: 5100 Stress Response based on Neuber Correction Method

The tutorial demonstrates Neuber shape optimization to reduce the stress concentration at the fillet regions by changing its shape.

The following exercises for shape optimization based on the Neuber correction method are:
  • Set up the objective, constraints and design variables in HyperMesh
  • Submit the job in OptiStruct
  • View results in HyperView
Figure 1 illustrates the structural model used for this tutorial. This tutorial has a half symmetric valve, which is modeled using first order tetra elements and has a force applied at a specific location. A shape ID is created using HyperMorph and the Shape design variables is created. A static stress response is created with Neuber correction and a mass response defined by an equation to minimize the change in the mass. A design constraint is defined for static-stress response.


Figure 1. Model and Loading Description
The optimization problem for this tutorial is stated as:
Objective
Minimize the change in mass.
Constraints
Minimize the static stress to a value less than 330 MPa using Neuber correction method
Design Variables
Shape design variables.

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

Set Up the Model

Create Property

Here you will assign the property, as the Valve is a solid part.

  1. In the Model Browser, right-click and select Create > Property.
  2. For Name, enter PSOLID.
  3. For Card Image, select PSOLID.
  4. For Material, select Aluminum.
  5. From the Model Browser, Component folder, right-click on Valve and assign the property PSOLID.

Create FORCE Load Collector

First you will isolate the component and then apply the force.

  1. In the Model Browser, right-click the component axis and select isolate only.
  2. In the Model Browser, right-click Create > Load Collector.
  3. For Name, enter Force.
  4. From the menu bar, click BCs > Create > Forces.
  5. Toggle to nodes and select the middle node of AXIS element (RBE2 element).
  6. Assign the values, as shown in the image.


    Figure 2. Apply Force
  7. Click Create.

Create SPC Load Collector

  1. In the Model Browser, right-click Create > Load Collector.
  2. For Name, enter SPC.
  3. From the menu bar, click BCs > Create > Constraints.
  4. For the first drop-down, toggle to nodes and for the second drop-down menu, select faces.
  5. Select only dof2, deselect all other dof's and click Create.
  6. Select the side face of the Valve and select only dof1, dof2 and dof3.
  7. Click Create.
  8. Similarly, select the other side faces of the Valve and select only dof1 and dof2.


    Figure 3. Apply Constraints

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 Static.
  3. For Type, select linear static from the drop-down menu.
  4. For SPC, select SPC.
  5. For load definition, select force.

Create Design Equations

  1. From the Analysis page, click optimziation > dequations panel.
  2. In the dequation = field, enter Delta.
  3. In the edit equation field, enter delta(x)=x-9.892e-02.


    Figure 4. Design Equation
  4. Click Create.
    Note: Delta(x) is changed in mass; where, x represents the actual mass with respect to change in shape and mass of the valve is 9.892e-02 kg.
  5. Select Tools > mass Calc > comps, select valve and click calculate.


    Figure 5. Mass of the Valve

Create Shape Optimization Design Variables

  1. From the Analysis page, click optimization > shape panel.
  2. Toggle desvar subpanel.
  3. For desvar =, enter Variable as the name of the design variable.
  4. For shape =, enter Fillet_shape.
  5. For lower bound, enter 0 and for upper bound, enter 1.
  6. Click create and return.


    Figure 6. Create Design Variables
    Note: Potential variation in the shape of the Valve is increased in the fillet radius.

Create Optimization Responses

Create Mass Response Type

  1. From the Analysis page, click optimization > responses panel.
  2. For response =, enter Mass.
  3. In the drop-down, select response type as mass and click return.


    Figure 7. Mass Response

Create Static Stress Response Type

  1. From the Analysis page, click optimization > responses panel.
  2. For response, enter Neuber.
  3. From the drop-down, select response type as static stress.
  4. Click elems > by collector > Valve.
  5. Select neuber correction, as shown in the image.


    Figure 8. Static Stress Response
  6. Click create.

Create Function Response Type

  1. In the Model Browser, select Optimization View.
  2. Right-click and select Window > Create > Responses.
  3. For Name, enter Delta_Response.
  4. For Function, select Response type.
  5. For Function ID, select Delta, which is Dequation created in previous step.
  6. Expand Response List and select Optiresponse as the Mass.


    Figure 9. Function Response Type
    Note: To create change in mass response, Function response is used here.

Create Design Constraints

  1. Click the dconstraints panel.
  2. In the constraint= field, enter Constraints.
  3. Click response type and select neuber, which is a static stress response type.
  4. Select loadstep 1 - Static.
  5. Select upper bound and enter a value of 330.

Define the Objective Function

  1. Click the objective panel.
  2. Verify that min is selected.
  3. Click response type and select Delta_Response.
  4. Click Create.

Define the SHAPE Card

Only displacement and stress results are available in the .h3d file by default. In order to look at displacement/stress results on top of a shape change that was applied to the model in HyperView, a SHAPE card needs to be defined.

  1. From the Analysis page, click the control cards panel.
  2. In the Card Image dialog, click SHAPE.
  3. Set FORMAT to H3D.
  4. Set TYPE to ALL.
  5. Set OPTION to ALL.
  6. Click return twice to go back to the main menu.

Run the Optimization

  1. From the Analysis page, click OptiStruct.
  2. Click save as.
  3. In the Save As dialog, specify location to write the OptiStruct model file and enter neuber_based_opimization 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 optimization.
  7. Set the memory options toggle to memory default.
  8. Click OptiStruct to run the optimization.
  9. Click Close.

View the Results

The following steps demonstrate how to review the contour plot with the optimized shape in HyperView.

It is helpful to view the deformed shape, stresses in the model using Neuber correction method .

  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. Go to the Results tab.
  3. In the Results tab, select Iteration 23, the last iteration.


    Figure 10. Iterations Selection
  4. On the Results toolbar, click resultsContour-16 to open the Contour panel.
  5. Set Result type to Shape change (v) and click Apply to contour the elements.


    Figure 11. Shape Change Results
  6. Set Result type to Neuber stress.
  7. Click Apply.


    Figure 12. Element Stress Results
  8. Open a new HyperView window and select File > Open > Session.
  9. Select <filename>_hist.mwv to view the change in mass and the optimized mass after changing the shape.


    Figure 13. Session Plot Results
  10. View the plot results for objective function and total mass.


    Figure 14. Objective Function and Total Mass Plot