OS-T: 1580 NLSTAT Analysis of Arthritic Finger

This tutorial demonstrates how to carry out nonlinear implicit analysis in OptiStruct involving hyperelastic material and contact.



Figure 1. FE Model
Figure 1 illustrates the structural model used for this tutorial: A hyperelastic Implant which is connected to bone on both sides, and a force of 9 Newtons is applied at one end of the model and the other end of the model is constrained with all DOFs.
You will be simulating an arthritic finger, which is modeled with hyperelastic material. The finger will be rotated approximately 90 degrees by applying a force of 9 Newtons. You will be running two iterations:
  • One with Tie contacts
  • One with Node2Surface contacts
The results, i.e. strains, displacements and stresses, will be reviewed 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 arthritis_finger.fem file you saved to your working directory from the optistruct.zip file. Refer to Access the Model Files.
  3. Click Open.
    The arthritis_finger.fem database is loaded into the current HyperMesh session, replacing any existing data.

Set Up the Model

Create Curves

Here you will create the curves for the hyperelastic material.

  1. In the Model Browser, right-click and select Create > Load Collector from the context menu.
  2. For Name, enter TABLES1100.
  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_NUM, enter 23, which is the number of rows in the table, and press Enter.
  6. Click on the Table icon table_pencil next to the Data field and enter the X (strain) and Y (stress) values:
    Table 1. Simple Tension Compression Data Values
    X (Strain) Y (Stress)
    1.1338 1.5506
    1.2675 2.4367
    1.3567 3.1013
    1.6242 4.2089
    1.8917 5.3165
    2.1592 5.981
    2.4268 6.8671
    3.051 8.8608
    3.586 10.6329
    4.0318 12.4051
    4.7898 16.1709
    5.3694 19.9367
    5.8153 23.481
    6.172 27.4684
    6.4395 31.0127
    6.707 34.557
    6.9299 38.3228
    7.0637 42.0886
    7.1975 45.6329
    7.3312 49.3987
    7.465 53.1646
    7.5541 56.9304
    7.6433 64.2405
  7. Click Close.
  8. Similarly, create the equi-biaxial tension data named TABLES1200 with the following X-Y data values:
    Table 2. Equi-biaxial Tension Data Values
    X (Strain) Y (Stress)
    1.02 0.9384
    1.06 1.59
    1.11 2.4087
    1.14 2.622
    1.2 3.324
    1.31 4.4278
    1.42 5.183
    1.68 6.6024
    1.94 7.7794
    2.49 9.7857
    3.03 12.6351
    3.43 14.6804
    3.75 17.4
    4.07 20.1058
    4.26 22.4502
    4.45 24.653
  9. Create the pure shear loading data named TABLES1400 with the following X-Y data values:
    Table 3. Pure Shear Loading Data Values
    X (Strain) Y (Stress)
    1.069 0.6
    1.1034 1.6
    1.1724 2.4
    1.2828 3.36
    1.4276 4.2
    1.8483 6
    2.3862 7.8
    3.0 9.6
    3.4897 11.12
    4.0345 12.96
    4.4483 14.88
    4.7793 16.58
    5.0621 18.2
  10. Create the volumetric data named TABLES1500 with the following X-Y data values:
    Table 4. Volumetric Data Values
    X (Strain) Y (Stress)
    0.9703 60
    0.9412 118.2
    0.9127 175.2
    0.8847 231.1

Create the Materials

This step requires creating the hyperelastic behavior of the implant material and the bone material.

  1. In the Model Browser, right-click and select Create > Material from the context menu.
  2. For Name, enter MATHE_2.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select MATHE from the drop-down menu.
  5. For MODEL, select ABOYCE from the drop-down menu.
  6. For TAB1, select load-curve (TABLES1100) created in Table 1.
  7. For TAB2, select load-curve (TABLES1200) created in Table 2.
  8. For TAB4, select load-curve (TABLES1400) created in Table 3.
  9. For TABD, select load-curve (TABLES1500) created in Table 4.
  10. In the Model Browser, right-click and select Create > Material.
  11. For Name, enter Bone.
  12. Click Color and select a color from the color palette.
  13. For Card Image, select MAT1 from the drop-down menu.
  14. For E, enter 14800.
  15. For NU, enter 0.3.

Create IMPLANT Property

The Implant property is created here.

  1. In the Model Browser, right-click and select Create > Property from the context menu.
  2. For Name, enter Implant.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select PLSOLID from the drop-down menu, and click Yes to confirm.
  5. For Material, select Unspecified > Material.
  6. In the Select Material dialog, select MATHE_2 from the list of materials and click OK.

Create BONE Property

Now you will create the Bone property.

  1. In the Model Browser, right-click and select Create > Property from the context menu.
  2. For Name, enter Bone.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select PSOLID from the drop-down menu and click Yes.
  5. For Material, click Unspecified > Material.
  6. In the Select Material dialog, select Bone from the list of materials and click OK.

Create PCONT Property

  1. In the Model Browser, right-click and select Create > Property from the context menu.
  2. For Name, enter Pcont.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select PCONT from the drop-down menu and click Yes.
  5. For STIFF_REAL_VAL, select STIFF = HARD.
  6. For MU1, enter 0.3.

Assign Properties

In this step, you will assign the properties to the respective components.

  1. In the Model Browser, click Component > Implant and enter the data, as shown below:


    Figure 2.
  2. Click Component > Bone1 and enter the data, as shown below:


    Figure 3.
  3. Click Component > Bone2 and enter the data, as shown below:


    Figure 4.

Apply Loads and Boundary Conditions

Create IMPLANT Set Segment

Here you will select the Implant faces and define the contact surface.

  1. Select the component Implant, right-click and select Isolate Only.
  2. In the Model Browser, right-click and select Create > Set Segment.
  3. For Name, enter Implant1.
  4. Click Color and select a color from the color palette.
  5. For Card Image, select SURF from the drop-down menu and click Yes.
  6. For Elements, select 0 Elements > Elements.
  7. In the Selection panel, select add solid faces from the first drop-down list.
  8. In the second drop-down list, select faces.


    Figure 5. Selection Panel
  9. Select all the faces of Implant using the Shift+Left mouse button.
  10. Click add and then return.

Create BONE Set Segment

Here you will select the Bone1 and Bone2 faces and define the contact surface.

  1. In the Entity Editor, select the components Bone1 and Bone2 together, right-click and select Isolate Only.
  2. In the Model Browser, right-click and select Create > Set Segment.
  3. For Name, enter Bone1.
  4. Click Color and select a color from the color palette.
  5. For Card Image, select SURF from the drop-down menu and click Yes.
  6. For Elements, select 0 Elements > Elements.
  7. In the Selection panel, select add solid faces from the first drop-down list.
  8. In the second drop-down list, select faces.


    Figure 6.
  9. Select the inside faces of Bone1 and Bone2 (selected faces are blue), as shown below:


    Figure 7. Bone1


    Figure 8. Bone2

Create Contacts

Here you will create a TIE contact and a N2S contact.

  1. In the Model Browser, right-click and select Create > Contact.
  2. For Name, enter Contact.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select TIE from the drop-down menu and click Yes.
  5. For SSID, select Unspecified > ContactSurf and select Implant1.
  6. For MSID, select Unspecified > ContactSurf and select Bone1.
  7. In the Model Browser, right-click and select Create > Contact.
  8. For Name, enter Contact.
  9. Click Color and select a color from the color palette.
  10. For Card Image, select CONTACT from the drop-down menu and click Yes.
  11. Toggle Property Option, and for PID, select Pcont.
  12. For SSID, select Unspecified > ContactSurf and select Implant1.
  13. For MSID, select Unspecified > ContactSurf and select Bone1.

Create NLPARM Load Collector

The nonlinear implicit parameters are defined.

  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.

Create NLADAPT Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter NLAdapt.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select NLADAPT from the drop-down menu.
  5. For NCUTS, enter 25.

Create NLMON Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter NLMon.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select NLMON from the drop-down menu.
  5. For ITEM, select DISP.
  6. For INT, select ITER.

Create NLOUT Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter NLOUT.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select NLOUT from the drop-down menu.
  5. For NINT, enter 10.
  6. Activate SVNONCNV and select YES for the VALUE.

Create CNTSTB Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter cntstb.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select CNTSTB from the drop-down menu.
  5. For S0, enter 0.01.
  6. For S1, enter 1e-05.

Create BCS Load Collector

Here the boundary condition is defined.

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter Bcs.
  3. Click Color and select a color from the color palette.
  4. Toggle Analysis panel and select constraints option.
  5. Select all the nodes of the Bone1 rear side face and activate all DOF and enter 0 for each.


    Figure 9.
  6. For Load Type, select SPC and click return.

Create Force

  1. In the Model Browser, right-click and select Create > Component.
  2. For Name, enter Rigid_Load.
  3. Click Color and select a color from the color palette.
  4. Select the 1D Panel and select rigids options.
  5. Toggle create.
  6. In the first drop-down, click multiple nodes and then click calculate node.
  7. In the second drop-down, select faces and click on rear face of Bone2 component.
  8. For elem types, enter RBE2 and click create.


    Figure 10.
  9. In the Model Browser, right-click and select Create > Load Collector.
  10. For Name, enter Force.
  11. Click Color and select a color from the color palette.
  12. Select the Analysis panel and click forces option.
  13. Toggle create and select center node of the Rigid_load element.
  14. For magnitude, enter 9 and select x-axis from the drop-down.
  15. For load types, enter FORCE.


    Figure 11.


    Figure 12.

Create Load Step

The nonlinear static load step is created.

  1. In the Model Browser, right-click and select Create > Load Step.
  2. For Name, enter load.
  3. Click Color and select a color from the color palette.
  4. For Analysis type, select Nonlinear static from the drop-down menu.
  5. For SPC, select Bcs from the list of load collectors.
  6. For LOAD, select Force.
  7. For NLPARM, select NLPARM.
  8. For NLADAPT, select NLAdapt.
  9. For NLOUT, select NLout.
  10. For CNTSTB, select cntstb.

Define Output Control Parameters

  1. From the Analysis page, select control cards.
  2. Click on GLOBAL_OUTPUT_REQUEST.
  3. Select ELFORCE, SPCF, STRAIN, and STRESS.
  4. In the control cards page, click GAPPRM, PARAM and FORMAT.
  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 13. Accessing the OptiStruct Panel
  2. Click save as.
  3. In the Save As dialog, specify location to write the OptiStruct model file and enter arthritis_finger.fem 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. Set the memory options toggle to memory default.
  8. Click OptiStruct to submit the job.
    If the job is successful, new results files should be in the directory where the arthritis_finger.out was written. The arthritis_finger.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:
    arthritis_finger.html
    HTML report of the analysis, providing a summary of the problem formulation and the analysis results.
    arthritis_finger.out
    ASCII based output file of the model check run before the simulation begins and gives nonlinear iteration history, as well as basic information on the results of the run.
    arthritis_finger.h3d
    HyperView compressed binary results file.
    arthritis_finger.stat
    Summary, providing CPU information for each step during analysis process.

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 14. Contour Panel
  6. Verify that the fields in the Contour panel match those in Figure 14 and click Apply.


    Figure 15. Contour of Displacement in Bone and Implant subject to Loading. TIE Contact Results


    Figure 16. Contour of Element Strains in Bone and Implant subject to Loading. TIE Contact Results


    Figure 17. Contour of Element Stresses in Bone and Implant subject to Loading . TIE Contact Results


    Figure 18. Contour of Displacement in Bone and Implant subject to Loading. N2S Contact Results


    Figure 19. Contour of Element Strains in Bone and Implant subject to Loading. N2S Contact Results


    Figure 20. Contour of Element Stresses in Bone and Implant subject to Loading . N2S Contact Results