OS-T: 1570 Nonlinear Transient Analysis of Crank Slider

In this tutorial, an existing finite element model of a slider crank is used to demonstrate how to perform Nonlinear Transient Analysis using OptiStruct.



Figure 1. Model and loading description

Figure 1 illustrates the structural model used for this tutorial: A slider crank is a four-link mechanism with three revolute joints and one sliding joint. The mechanism is completely constrained at one end and transient load are applied to the crank. The Nonlinear transient analysis is run for a total time of 1 second with the time being divided into 200 increments (that is step time equal to 0.005). A concentrated mass is defined at the center of the slider, both the ends of the crank and at the hinge.

You will be simulating the Slider Crank model, which is modeled using beam/solid elements. The crank will be rotated by 180 degrees about the z-axis by applying appropriate constraints on other links.

Two iterations will be run:
  • One with beam elements
  • One with solid elements
The results 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.

Import the Model

  1. Click File > Import > Solver Deck.
    An Import tab is added to your tab menu.
  2. Select the crank_slider.fem file you saved to your working directory from the optistruct.zip file. Refer to Access the Model Files.
  3. Click Import, then click Close to close the Import tab.

Set Up the Model

Create the Material

Here you will define the beam material.

  1. In the Model Browser, right-click and select Create > Material from the context menu.
  2. For Name, enter material1.
    A new material, material1 has beeen created.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select MAT1.
  5. For E, enter 200000000.
  6. For N, enter 0.3.
  7. For RHO, enter 7200.

Create the Property

Here you will create the beam property.

  1. In the Model Browser, right-click and select Create > Property from the context menu.
  2. For Name, enter property1.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select PBEAM and click Yes to confirm.
  5. For Material, click Unspecified > Material.
  6. In the Select Material dialog, select material1 and click OK.
  7. For Beam Section, click Unspecified > Beamsection.
  8. In the Select Beamsection dialog, select rectangle section 2 from the list and click OK.
  9. In the Model Browser, click the component BEAM.
  10. Under property window for BEAM, for Property, select property1 that you created.
    The Material option is auto-filled.

Apply Loads and Boundary Conditions

Create TABLED1 Curve

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


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

Create TSTEP Load Collector

Here you will define the hyper elastic behavior of the implant.

  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 200.
  7. For DT, enter the time increment of 0.005.
  8. Click Close.

Create SPC Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter spc.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select NONE.
  5. Click BCs > Create > Constraints to open the Constraints panel.
  6. Click nodes, then select the hinge node (Figure 3).
  7. Check all Degrees of Freedom (DoFs) and enter 0 for each DoF.
  8. For load types=, select SPC.
  9. Click create.


    Figure 3. Constraints on the hinge node
  10. Select the center node of the slider (Figure 4) and check all DoFs; except DoF1.
    This indicates that DoF1 is the only active Degree of Freedom.


    Figure 4. Constraints applied to the center node of the slider
  11. Click create and return.
    The 2 nodes are constrained in one load collector.

Create SPC1 Load Collector

One more SPC load collector is created here.

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter spc1.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select NONE.
  5. Click BCs > Create > Constraints to open the Constraints panel.
  6. Select the first node of crank (Figure 5).
  7. Check all Degrees of Freedom (DoFs); except DoF4.
    This indicates that DoF4 is the only active Degree of Freedom.


    Figure 5. Constraints applied to the first node of the crank
  8. For load types=, select SPC.
  9. Click create and return.

Create SPCD Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter spcd.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select NONE.
  5. Click BCs > Create > Constraints to open the Constraints panel.
  6. Select the first node of crank.
  7. Check all Degrees of Freedom (DoFs); except dof4 and enter a value of 0 against each checked DoF.
  8. For DoF4, enter the value 3.14.
    This indicates that DoF4 is the only active Degree of Freedom.
  9. For load types=, select SPCD.
  10. Click create and return.

Create SPCADD Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter spcadd.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select SPCADD.
  5. For SPCADD_Num_Set, enter 2.
  6. Click on the Table icon table_pencil next to the Data field and enter the following values in the pop-out window as:


    Figure 6.

Create GRAV Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter grav.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select GRAV from the drop-down menu.
  5. For G, enter 9.8.
  6. For N3, enter -1.0.

Create a TLOAD Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter tload.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select TLOAD1 from the drop-down menu.
  5. For EXCITEID , click Unspecified > Loadcol.
  6. In the Select Loadcol dialog, select spcd from the list of load collectors.
  7. For Type, select VELO.
  8. For TID, select table.

Create tload_grav Load Collector

Here you will create one more TLOAD load collector.

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter tload_grav.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select TLOAD from the drop-down menu.
  5. For EXCITED, click Unspecified > Loadcol.
  6. In the Select Loadcol dialog, select grav from the list of load collectors.
  7. For TYPE, select LOAD.
  8. For TID, select table.

Create a DLOAD Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter dload.
  3. Click Color and select a color from the color palette.
  4. For Card Image, select DLOAD from the drop-down list.
  5. For S (scale factor), enter 1.0.
  6. For DLOAD_NUM, enter 2.
  7. For data, select tload and tload_grav and enter 1 as the scale factor for both loads.
  8. Click Close.

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.
  5. For NINC, enter 200.
  6. For DT, enter 0.005.

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 DTMAX, enter 0.005.
  6. For DTMIN, enter 1e-05.
  7. For NCUTS, enter 5.

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 200.

Create Load Steps

The nonlinear transient load step is created here.

  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 loadstep1.
  3. Click Color and select a color from the color palette.
  4. For Type, select Nonlinear transient from the drop-down menu.
  5. For SPC, select spcadd.
  6. For TSTEP, select tstep.
  7. For NLPARM(LGDISP), select nlparm.
  8. For DLOAD, select dload from the list of load collectors.
  9. For NLADAPT, select NLAdapt from the list of load collectors.
  10. For NLOUT, select NLout.
  11. Under SUBCASE OPTION, toggle ANALYSIS, and select TYPE DTRAN.
  12. Activate NLMON and select nlmon.

Define Output Control Parameters

  1. From the Analysis page, select control cards.
  2. Click on GLOBAL_OUTPUT_REQUEST.
  3. Below DISPLACEMENT, set Option to Yes and select H3D for Output format.
  4. For Format, select PLOT.
  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 7. Accessing the OptiStruct Panel
  2. Click save as.
  3. In the Save As dialog, specify location to write the OptiStruct model file and enter crank_slider.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 crank_slider.fem was written. The crank_slider.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:
    slider_crank.html
    HTML report of the analysis, providing a summary of the problem formulation and the analysis results.
    slider_crank.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.
    slider_crank.h3d
    HyperView compressed binary results file.
    slider_crank.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 8. Contour Panel
  6. Verify that the fields in the Contour panel match those in Figure 8 and click Apply.


    Figure 9. Contour of Displacement in Crank Slider model. Time = 0.98 sec