OS-T: 1365 NLSTAT Analysis of Solid Blocks in Contact

This tutorial demonstrates how to carry out nonlinear implicit small displacement analysis in OptiStruct, involving elasto-plastic materials, contact and continuing the nonlinear solution sequence from a preceding nonlinear loadcase.

Figure 1 illustrates the structural model used for this tutorial, which is two square solid blocks made of elasto-plastic steel material. The dimensions of the blocks and the material parameters are outlined below Figure 1.

In the first nonlinear subcase, pressure loading is be applied to the top solid block, the top corners of which are constrained in X and Y directions. The top solid is in contact with the bottom solid, the bottom corners of which are constrained in X, Y and Z directions. The second nonlinear subcase is to simulate the unloading and is a continuation of the nonlinear solution sequence from the previous loading subcase.

rd2100_nlstat_model
Figure 1. Model and Loading Description
Units
Length: mm; Time: s; Mass: Mgg; (Force: N; Stress: MPa)
Top block
72 mm x 72 mm
Bottom block
100 mm x 100 mm
Thickness of blocks
20. mm
Material
Steel, Elasto-plastic
Initial density ( ρ ): 7.90e-9 kg/mm3
Young's modulus (E): 210000 MPa
Poisson coefficient ( υ ): 0.3
Yield Stress ( σ 0): 850.0 MPa
Imposed pressure
1000.0 MPa, applied at the center of top block

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

Set Up the Model

Create Elasto-plastic Material

First, the stress versus plastic strain curve for the material needs to be defined.

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter stress-strain.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select TABLES1 from the drop-down menu.
  5. For TABLES1_NU, enter 2 (number of rows in the table), and press Enter.
  6. Click the Table icon table_pencil next to the Data field and enter the following values (x and y fields) in the pop-up window.

    os1365_tables1_num
    Figure 2.
  7. Click Close to close the dialog.

Now, the elasto-plastic material needs to be updated.

  1. In the Model Browser, click the material steel.
    The Entity Editor opens.
  2. Click on the checkbox next to MATS1 to define the elastic-plastic material for NLSTAT analysis.
  3. For TID, click Unspecified > Loadcol.
  4. In the Select Loadcol dialog, select the stress_strain load collector and click OK.
  5. Input the values, as shown below.
    TYPSTRN of 1 signifies specifying stress (Y) versus plastic strain (X).

    os1365_steel_13
    Figure 3.

Define Contact between the Blocks

The contact surfaces for the two blocks need to be defined.
  1. In the Model Browser, right-click and select Create > Set.
  2. For Name, enter top.
  3. For Card Image, select SET_ELEM from the drop-down menu.
  4. Leave the Set Type switch set to non-ordered type.
  5. For Entity IDs, click Unspecified > Property.
  6. In the Select Properties dialog, select the top solid block Solid1 and click OK.

    os1365_top_13
    Figure 4.
  7. Similarly create another set named bottom.
  8. Repeat steps 3 through 6 for bottom block select the bottom solid Solid2.

Next, the interface needs to be defined.

  1. In the Model Browser, right-click and select Create > Contact.
  2. For Name, enter SOLID_CONTACT.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select CONTACT from the drop-down menu.
  5. For MSID (main surface) and select Set from the extended selection menu.
  6. Click the yellow Set panel and select the bottom block bottom in the pop-up window and click OK.
  7. Similarly, for SSID (secondary surface), select the top set.
  8. For TYPE, select SLIDE from the drop-down menu.
  9. For MORIENT, select NORM from the drop-down menu.

    os1365_solid_contact
    Figure 5.

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. Enter the values as shown below:

    os1365_nlparm_13
    Figure 6.

Create the 1st Nonlinear Load Step

  1. In the Model Browser, right-click and select Create > Load Step.
  2. For Name, enter loading.
  3. Click Analysis type and select Nonlinear quasi-static from the drop-down menu.
  4. For SPC, click Unspecified > Loadcol.
  5. From the Select Loadcol dialog, select SPC from the list of load collectors and click OK.
  6. For LOAD, click Unspecified > Loadcol.
  7. From the Select Loadcol dialog, select pressure from the list of load collectors and click OK.
  8. For NLPARM, click Unspecified > Loadcol.
  9. From the Select Loadcol dialog, select nlparm from the list of load collectors and click OK.

    os1365_loading_13
    Figure 7.

Create the 2nd Nonlinear Load Step

  1. For Name, enter unload.
  2. The Analysis type should, again, be set to Nonlinear quasi-static from the drop-down menu.
  3. For SPC, select SPC from the list of load collectors.
  4. For NLPARM, select nlparm from the list of load collectors.

    Checkpoint: The unloading subcase (unload) does not contain the pressure load applied during the loading subcase (loading).


    os1365_unload_13
    Figure 8.

Define Output Control Parameters

  1. From the Analysis page, select control cards.
  2. Click on GLOBAL_OUTPUT_REQUEST.
  3. Below CONTF, DISPLACEMENT, STRAIN and STRESS, set Option to Yes.
  4. Under STRAIN, set TYPE(1) to PLASTIC.
  5. Click return twice to go to the main menu.

Submit the Job

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

    OS_1000_13_17
    Figure 9. Accessing the OptiStruct Panel
  2. Click save as.
  3. In the Save As dialog, specify location to write the OptiStruct model file and enter nlstat_complete 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 nlstat_complete.fem was written. The nlstat_complete.out file is a good place to look for error messages that could help debug the input deck if any errors are present.

View the Results

Using HyperView, plot the Displacement, the von Mises stress, plastic strains and contact pressure contours at the end of the 1st (loading) step.

rd2100_displacement_1
Figure 10. Contour of Displacements in Blocks Subject to Loading

rd2100_von_mises
Figure 11. Contour of von Mises Stress in Blocks Subject to Loading

rd2100_plastic_strain
Figure 12. Contour of Plastic Strains in the Blocks Subject to Loading

os1365_contact_pressure
Figure 13. Contour of Contact Pressure in the Block Interface After the 1st (Loading) Subcase

Next, change the subcase to the 2nd that is unloading subcase and plot the displacement contour to see the change in displacements in the blocks subject to unloading.


rd2100_displacement_2
Figure 14. Contour of Displacements in Blocks Subject to Unloading in 2nd Subcase