OptiStruct 2017 Tutorials and Examples
OptiStruct 2017 Tutorials and Examples
Special Notice: Pre-release versions of Altair software are provided ‘as is’, without warranty of any
kind. Usage is strictly limited to non-production purposes.
The Altair web site is a valuable online companion to Altair software. Visit www.altairhyperworks.com
for tips and tricks, training course schedules, training/tutorial videos, and other useful information.
Altair training courses provide a hands-on introduction to our products, focusing on overall
functionality. Courses are conducted at our main and regional offices or at your facility. If you are
interested in training at your facility, please contact your account manager for more details. If you do
not know who your account manager is, please send an e-mail to training@altair.com and your account
manager will contact you.
When contacting Altair support, please specify the product and version number you are using along
with a detailed description of the problem. Many times, it is very beneficial for the support engineer to
know what type of workstation, operating system, RAM, and graphics board you have, so please have
that information ready. If you send an e-mail, please specify the workstation type, operating system,
RAM, and graphics board information in the e-mail.
To contact an Altair support representative, reference the following table or the information available
on the HyperWorks website:
http://www.altairhyperworks.com/ClientCenterHWSupportProduct.aspx
For questions or comments about this help system, send an email to hwsupport@altair.com.
In addition, the following countries have resellers for Altair Engineering: Colombia, Czech Republic,
Ecuador, Israel, Russia, Netherlands, Turkey, Poland, Singapore, Vietnam, Indonesia
Official offices with resellers: Canada, China, France, Germany, India, Malaysia, Italy, Japan, Korea,
Spain, Taiwan, United Kingdom, USA
See www.altair.com for complete contact information.
Table of Contents
Verification Problems...................................................................................................................................... 3
Scordelis-Lo Roof................................................................................................................................... 34
Example Problems for Topology Optimization Using Minimum Member Size Control .......................... 982
Example Problems for Topography Optimization Using Pattern Grouping .......................................... 1005
The following verification problems using OptiStruct can be found in this section:
Raasch Challenge - curved strip hook
NAFEMS Test Problem (LE3) - Radial Point Load on a Hemisphere
NAFEMS Frequency Benchmarks
MacNeal-Harder Test Problems
Summary
Raasch challenge is a curved strip hook problem with a tip in-plane shear load, posed in 1990 by
Ingo Raasch of BMW in Germany. The problem poses a significant challenge to shell elements
because of the inherent coupling between three modes of deformation: bending, extension, and
twist. OptiStruct is benchmarked against the Raasch challenge to assure its shell elements
performance on linear static analysis.
Benchmark Model
Three types of elements are used for linear static analysis of the Raasch's hook. They are tria-
shell, quad-shell, and hexa-solid elements. For each type of element, five different mesh
densities (1x9, 3x17, 5x34, 10x68, and 20x136) are used to investigate the convergence of the
solution. Two elements are used in the thickness direction for solid elements, thus, the mesh
densities become 1x9x2, 3x17x2, 5x34x2, 10x68x2, and 20x136x2. The target is the z-
directional displacement at point A, the upper end corner of the hook. The target value is 4.9366
inches.
Linear Static Analysis Results
All results are normalized with the target value.
Element
Type
1x9 3 x 17 5 x 34 10 x 68 20 x 136
Reference
Knight, Jr. N. F., Raasch Challenge for Shell Elements, AIAA Journal, Vol. 35, No. 2, February
1997.
Summary
This is the NAFEMS test problem LE3. The model is a hemispherical shell subjected to
concentrated radial loads at its free edges. It examines the performance of the three-dimensional
shell to model local bending behavior under conditions where the deformations are primarily due
to bending.
Benchmark Model
4-node, first order CQUAD4 elements are benchmarked in LE3. The hemisphere is 10m in radius
and 0.04m in radial thickness. Two pairs of identical loads, 4000N, are applied at the free edge of
the hemisphere, and are at right angles to each other. One pair of the loads is directed inwards
(toward the center) of the hemisphere, while the second pair is directed outward from the center,
producing deformation of compression in one direction and elongation in another. Since both the
geometry and loads are symmetrical, only a quarter of the hemisphere is modeled. Symmetric
boundary constraints are applied on edges AE and CE. The z-translation at point E is fixed, and
all displacements on edge AC are free. The test also requires the mesh of the hemisphere to have
equally spaced nodes on edges AC, CE, EA, BG, DG, and FG. The target is x-translation at point
A, with a target value of 0.185m. The material properties for the hemisphere are: E = 68.25 GPa
and = 0.3.
Summary
Test No. FV2 is a pin-ended double cross, in-plane vibration problem. OptiStruct is used to
investigate the coupling between flexural and extensional behavior, as well as the repeated and
closed eigenvalues of the double cross in normal modes analysis.
Benchmark Model
The 2-node simple beam elements are used to model the double cross. Each arm of the cross
consists of four elements. The x- and y-translation displacements are fixed at the end of all arms.
The material properties for the cross beams are:
E = 200x109 N/m2 and = 8000 kg/m3
File Location
<install_directory>/demos/hwsolvers/optistruct/fv2.fem
Summary
Test No. FV4 is a cantilever beam with off-center point masses. The problem is set up to test the
behavior of coupling between torsional and flexural, off-center inertial axis, and close
eigenvalues. OptiStruct is used in this normal modes analysis.
Benchmark Model
The 2-node simple beam elements are used to model the problem. The cantilever beam consists
of five elements. Two rigid elements are used to connect mass 1, 10000 kg, and mass 2, 1000
kg, at the end of the cantilever beam. All degrees of freedom at end A are constrained. The
material properties for the cantilever beam are:
File Location
<install_directory>/demos/hwsolvers/optistruct/fv4.fem
Reference
NAFEMS Ltd, The Standard NAFEMS BENCHMARKS TNSB Rev. 3, NAFEMS Ltd, Scottish Enterprise
Technology Park, Whitworth Building, East Kilbride, Glasgow, United Kingdom, 1990.
Summary
Test No. FV5 is a simply-supported beam problem. OptiStruct tests the repeated eigenvalues,
shear deformation, and rotary inertia. The possibility of missing extensional modes is also
investigated.
Benchmark Model
The 2-node beam elements are used to model the simply-supported beam which consists of five
elements. The displacements in x, y, and z direction, as well as the rotation in x direction are
fixed at the end A. In addition, the displacements in y and z direction are constrained at end B.
The material properties for the beam are:
E = 200x109 N/m2, = 0.3 and = 8000 kg/m3
File Location
<install_directory>/demos/hwsolvers/optistruct/fv5.fem
Reference
NAFEMS Ltd, The Standard NAFEMS BENCHMARKS TNSB Rev. 3, NAFEMS Ltd, Scottish Enterprise
Technology Park, Whitworth Building, East Kilbride, Glasgow, United Kingdom, 1990.
Summary
Test No. FV12 is a thin square plate model, which contains three rigid body modes. OptiStruct
investigates the repeated eigenvalues and kinematically incomplete suppressions of the thin plate
in normal modes analysis.
Benchmark Model
The 4-node quad elements with 8x8 mesh configuration are used to model the plate. The x-
translation, y-translation displacements and z-rotational displacements are fixed for all nodes.
The material properties for the thin plate are:
E = 200x109 N/m2, = 0.3 and = 8000 kg/m3
File Location
<install_directory>/demos/hwsolvers/optistruct/fv12.fem
Reference
NAFEMS Ltd, The Standard NAFEMS BENCHMARKS TNSB Rev. 3, NAFEMS Ltd, Scottish Enterprise
Technology Park, Whitworth Building, East Kilbride, Glasgow, United Kingdom, 1990.
Summary
Test No. FV15 is a thin rhombic plate problem. The behavior of distorted thin elements in normal
modes analysis is examined using OptiStruct.
Benchmark Model
The 4-node quad elements with 12x12 mesh configuration are used to model the thin rhombic
plate. The x-translation, y-translation displacements, and z-rotational displacements are fixed for
all of the nodes. In addition, the z-translation, x-rotational and y-rotational are fixed for the
nodes along the four edges. The material properties for the rhombic plate are:
E = 200x109 N/m2, = 0.3 and = 8000 kg/m3
File Location
<install_directory>/demos/hwsolvers/optistruct/fv15.fem
Reference
NAFEMS Ltd, The Standard NAFEMS BENCHMARKS TNSB Rev. 3, NAFEMS Ltd, Scottish Enterprise
Technology Park, Whitworth Building, East Kilbride, Glasgow, United Kingdom, 1990.
Summary
Test No. FV22 is a thick rhombic plate problem. The behavior of distorted, thick elements in
normal modes analysis is examined using OptiStruct.
Benchmark Model
The 4-node quad elements with 10x10 mesh configuration are used to model the thick rhombic
plate. The x-translation, y-translation displacements, and z-rotational displacements are fixed for
all of the nodes. In addition, the z-translation, x-rotational and y-rotational are fixed for the
nodes along the four edges. The material properties for the thick plate are:
File Location
<install_directory>/demos/hwsolvers/optistruct/fv22.fem
Reference
NAFEMS Ltd, The Standard NAFEMS BENCHMARKS TNSB Rev. 3, NAFEMS Ltd, Scottish Enterprise
Technology Park, Whitworth Building, East Kilbride, Glasgow, United Kingdom, 1990.
Summary
Test No. FV32 is a tapered membrane problem with irregular mesh. The geometry and mesh
symmetry are maintained along the x-direction. The shear behavior of membrane elements in the
normal modes analysis is examined using OptiStruct.
Benchmark Model
The 4-node quad elements with 16x8 mesh configuration are used in test FV32. The
displacements in z direction are fixed on all nodes and the x- and y-translation displacements are
fixed for the nodes along the y-axis. The material properties for the membrane are:
E = 200x109 N/m2, = 0.3 and = 8000 kg/m3
File Location
<install_directory>/demos/hwsolvers/optistruct/fv32.fem
Reference
NAFEMS Ltd, The Standard NAFEMS BENCHMARKS TNSB Rev. 3, NAFEMS Ltd, Scottish Enterprise
Technology Park, Whitworth Building, East Kilbride, Glasgow, United Kingdom, 1990.
Summary
Test No. FV52 is a well-established solid square plate, which contains three rigid modes with the
given boundary condition. The test examines the performance of OptiStruct normal modes
analysis on 3D solid elements.
Benchmark Model
The HEXA 8-node and HEXA 20-node solid elements are used in test FV52. An 8x8x3 mesh
configuration is used for HEXA 8-node elements, and a 4x4x1 mesh configuration is used with
HEXA 20-node elements. The z-directional displacement is constrained at Z = -5m plane along
the four edges of the plate. The material properties for the solid plate are:
E = 200x109 N/m2, = 0.3 and = 8000 kg/m3
f* 193.59 Hz f* 206.19 Hz
HEXA20 1.0017 HEXA20 1.0022
HEXA8 1.0164 HEXA8 1.0164
File Location
<install_directory>/demos/hwsolvers/optistruct/fv52_HOE.fem
<install_directory>/demos/hwsolvers/optistruct/fv52_LOE.fem
Length = 6.0; Height = 0.2; Depth = 0.1; E = 1.0E7: Poisson's ratio = 0.3.
Summary
This is a straight cantilever beam solved with solid and shell elements. Three models
(rectangular, parallelogram, trapezoidal) are made with each element's type to investigate the
effect of distorted elements with a high aspect ratio.
Benchmark Model
Six types of elements are used for this problem. They are tria-shell, quad-shell, and hexa-solid
elements, each with 1st and 2nd order. Four loading cases are used for each model; extension,
in-plane bending, transverse bending, and twist. For the extension and bending load cases, unit
loads are applied in a consistent fashion over all of the nodes at the tip of the beam. For the twist
load cases, a unit moment is applied at the tip.
Theoretical solutions for the deflections at the tip, computed by beam theory, are as follows.
extension UX 0.00003
(a) Rectangular
(b) Parallelogram
(c) Trapezoidal
File Location
<install_directory>/demos/hwsolvers/optistruct/CBhex.fem
<install_directory>/demos/hwsolvers/optistruct/CBpent1.fem
<install_directory>/demos/hwsolvers/optistruct/CBpen2.fem
<install_directory>/demos/hwsolvers/optistruct/CBquad.fem
<install_directory>/demos/hwsolvers/optistruct/cbquad2.fem
<install_directory>/demos/hwsolvers/optistruct/cbquad4.fem
Reference
MacNeal, R.H., and Harder, R.L., A Proposed Standard Set of Problems to Test Finite Element
Accuracy, Finite Elements in Analysis and Design, 1 (1985) 3-20.
Summary
This is a curved cantilever beam solved with solid and shell elements. A model is made with each
element's type to investigate the effect of distorted elements with a high aspect ratio.
Benchmark Model
Six types of elements are used for this problem. They are tria-shell, quad-shell, and hexa-solid
elements, each with 1st and 2nd order. Two loading cases are used for each model; in-plane
bending, transverse bending. For both load cases, unit loads are applied in a consistent fashion
over all of the nodes at the tip of the beam.
Theoretical solutions for the deflections at the tip, computed by beam theory, are as follows.
File Location
<install_directory>/demos/hwsolvers/optistruct/CBhex.fem
<install_directory>/demos/hwsolvers/optistruct/CBpent1.fem
<install_directory>/demos/hwsolvers/optistruct/CBpen2.fem
<install_directory>/demos/hwsolvers/optistruct/CBquad.fem
<install_directory>/demos/hwsolvers/optistruct/cbquad2.fem
<install_directory>/demos/hwsolvers/optistruct/cbquad4.fem
<install_directory>/demos/hwsolvers/optistruct/CBtri1.fem
<install_directory>/demos/hwsolvers/optistruct/CBtri2.fem
Reference
MacNeal, R.H., and Harder, R.L., A Proposed Standard Set of Problems to Test Finite Element
Accuracy, Finite Elements in Analysis and Design, 1 (1985) 3-20.
Length = 12.0; Width = 1.1; Depth = 0.32; E = 29.0E6; Poisson’s ratio = 0.22; Loading = unit forces at tip.
Summary
This is a twisted cantilever beam solved with solid and shell elements. A model is made with each
element's type to investigate the effect of distorted elements with a high aspect ratio.
Benchmark Model
Six types of elements are used for this problem. They are tria-shell, quad-shell, and hexa-solid
elements, each with 1st and 2nd order. Two loading cases are used for each model; in-plane
bending, transverse bending. For both load cases, unit loads are applied in a consistent fashion
over all of the nodes at the tip of the beam.
Theoretical solutions for the deflections at the tip, computed by beam theory, are as follows.
File Location
<install_directory>/demos/hwsolvers/optistruct/CBhex.fem
<install_directory>/demos/hwsolvers/optistruct/CBpent1.fem
<install_directory>/demos/hwsolvers/optistruct/CBpen2.fem
<install_directory>/demos/hwsolvers/optistruct/CBquad.fem
<install_directory>/demos/hwsolvers/optistruct/cbquad2.fem
<install_directory>/demos/hwsolvers/optistruct/cbquad4.fem
<install_directory>/demos/hwsolvers/optistruct/CBtri1.fem
<install_directory>/demos/hwsolvers/optistruct/CBtri2.fem
Reference
MacNeal, R.H., and Harder, R.L., A Proposed Standard Set of Problems to Test Finite Element
Accuracy, Finite Elements in Analysis and Design, 1 (1985) 3-20.
R = 25.0; Length = 50.0, t = 0.25; E = 4.32E8; Poisson's ratio = 0.0; Ux = Uz = 0 on curved edge.
Summary
The Scordelis-Lo Roof is a classical benchmark problem for shell elements. Analytical and
experimental investigations were initially performed by Scordelis and Lo.
Benchmark Model
The roof structure is supported on both ends, and loaded by self weight of 90 pounds per square
foot, with only one quadrant modeled. Six types of elements are used for this problem. They are
tria-shell, quad-shell, and hexa-solid elements, each with 1st and 2nd order. Each element type
was benchmarked with different mesh density to check the convergence. As the original study
was related to concrete structure, Poisson's ratio is set to zero. The structure is curved and its
solution will exhibit membrane and bending behavior. Target solution for the vertical deflection at
the midpoint of the free edge is 0.3024.
2 4 6 8 10
QUAD4 1.487 1.082 1.035 1.020 1.013
QUAD8 1.436 1.179 1.139 1.116 1.102
File Location
<install_directory>/demos/hwsolvers/optistruct/CBhex.fem
<install_directory>/demos/hwsolvers/optistruct/CBpent1.fem
<install_directory>/demos/hwsolvers/optistruct/CBpen2.fem
<install_directory>/demos/hwsolvers/optistruct/CBquad.fem
<install_directory>/demos/hwsolvers/optistruct/cbquad2.fem
<install_directory>/demos/hwsolvers/optistruct/cbquad4.fem
<install_directory>/demos/hwsolvers/optistruct/CBtri1.fem
<install_directory>/demos/hwsolvers/optistruct/CBtri2.fem
Reference
MacNeal, R.H., and Harder, R.L., A Proposed Standard Set of Problems to Test Finite Element
Accuracy, Finite Elements in Analysis and Design, 1 (1985) 3-20.
Exercise
plate.out ASCII output file containing specific information on the file set up, the
set up of your optimization problem, estimate for the amount of RAM
and disk space required for the run, information for each optimization
iteration, and compute time information. Review this file for warnings
and errors.
plate.stat Summary of analysis process, providing CPU information for each step
during analysis process.
While still in HyperMesh, launch HyperView after the job has finished from the OptiStruct panel
by clicking HyperView. HyperView opens and automatically loads the H3D file from the
OptiStruct job for post-processing.
To execute a check run to validate your input deck and determine how much RAM and
disk space is necessary for the run, at the command prompt, enter:
$HWSDIR/<solver_name> plate.fem -check
Information regarding memory requirements is written to the file plate.out.
See the Running OptiStruct section of the OptiStruct User's Guide for more detailed information.
Exercise
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
2. Select the plate_hole.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
3. Click Open. The plate_hole.hm database is loaded into the current HyperMesh session,
replacing any existing data. The database only contains geometric data.
When building models, it is encouraged to create the material and property collectors before
creating the component collectors. This is the most efficient way of setting up the file since
components need to reference materials and properties.
If a material property does not display a value next to it, it is turned off. To edit these
material properties, click in the blank Value fields next to them and enter the required
values.
A new material, steel, has been created. The material uses OptiStruct's linear isotropic
material model, MAT1. For E (Young's Modulus), enter 2E+05 and for NU (Poisson's Ratio),
enter 0.3. It is not necessary to define a density value since only a static analysis is
performed. Density values may be required, however, for other solution sequences.
At any time, the card image for this collector can be modified by clicking on steel in the
Model browser or by using the Card Editor .
If a property field does not display a value next to it, it is turned off. To edit these properties,
click on the blank Value fields next to them and enter the required values.
A new property, plate_hole, has been created as a 2D PSHELL. Material information is also
linked to this property.
6. In the Model browser, click on the component plate_hole. The component fields are
displayed in the Entity Editor below.
7. For Property, click Unspecified > Property.
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.
5. Create another load collector named forces using steps 5.1 to 5.3.
2. Click BCs > Create > Constraints to open the Constraints panel, as shown in Figure 8.
4. Use Shift + Left Click + Hold to select the nodes, as shown in Figure 10.
5. Constrain dof1, dof2, dof3, dof4, dof5, and dof6 and set all of them to a value of 0.0, as
shown in Figure 11.
Dofs with a check are constrained, while dofs without a check are 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.
4. Use SHIFT + Left Click + Hold to select the nodes shown in Figure 13 using the circular
selection window.
5. Set the coordinate system toggle to global system, as shown in Figure 14.
6. Click the vector definition switch and select constant vector.
7. Click magnitude = and enter 21.277 (that is 1000 divided by the number of nodes 47).
8. Click the direction definition switch below magnitude =, and select z-axis from the pop-
up menu.
2. Click save as following the input file: field. A Save file browser window opens.
3. Select the directory where you would like to write the OptiStruct model file and enter the
name for the model, plate_hole.fem, in the File name: field. For OptiStruct input decks
.fem is the recommended extension.
4. Click Save.
The name and location of the plate_hole.fem file displays in the input file: field.
5. Set the export options: toggle to all.
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.stat Summary, providing CPU information for each step during analysis
process.
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), multi-body system simulation, video and engineering data.
Figure 21: The vonMises Stress plot for the given subcase
At this stage, you can try to answer the following questions to test the 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?
Figure 22: Isometric view of the deformed plot overlaid on the undeformed mesh (model units is set to 500)
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
2. Select the coffee_lid.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
3. Click Open.
The coffee_lid.hm database is loaded into the current HyperMesh session, replacing any
existing data. The database only contains geometric data.
When building models, it is encouraged to create the material and property collectors before
creating the component collectors. This is the most efficient way of setting up the file since
components need to reference properties and materials.
6. Repeat steps 4.1 to 4.5 to update the property PSHELL1 and assign the material plastic to
it.
The property collectors and component collectors, PSHELL and PSHELL1, now reference the
material plastic. The component collectors that reference the corresponding properties are
automatically updated with the specified material. If you access the Entity Editor and edit
either of these property or component collectors, notice that the Material fields are now all set
to plastic(1).
Thermal loading has already been applied to the model. In the following steps, constraints will be
applied to the model.
lid_complete.out OptiStruct output file containing specific information on the file set
up, the set up of your optimization problem, estimates for the
amount of RAM and disk space required for the run, and compute
time information. Review this file for warnings and errors.
Displacement and Stress results are output from OptiStruct for linear static analyses by default.
The following steps describe how to view those results in HyperView.
In this analysis, the region around the hinges may be a concern. There are relatively high stress
values that must be resolved. For instance, if testing shows that the coffee pot lid wears out
around the hinge area over time, these thermal stresses could possibly cause that fatigue.
Exercise
3. Similarly, input the material values next to the corresponding fields as shown in Figure 1.
Click ENTER after inputting the values. A new material, steel, has been created.
The material uses OptiStruct linear isotropic material model, MAT1. For E (Young's Modulus),
enter 2E+05; for NU (Poisson's Ratio), enter 0.3 and for RHO (material density), enter 7.85E-
09. A material density is required for the normal modes solution sequence.
If a material property does not display a value next to it, it is turned off. To edit these
material properties, click on the blank Value fields next to them and enter the required
values.
At any time, the card image for this collector can be modified by clicking on steel in the
Model browser or by using the Card Editor .
4. In the Model browser, click on the property design to display the property entry in the Entity
Editor. Change the thickness value (T) from 0.3 to 0.25 as shown in Figure 2.
5. Click on the field next to the Material field and then click on the Material selector.
6. In the Select Material dialog, select the newly created steel material and click OK. This
changes the material from gn to steel on the property design.
7. Similarly, on the property nondesign, update the thickness value (T) from 0.3 to 0.25 and
change the material from gn to steel.
Figure 2: Updating the thickness value for design and nondesign property entries
The model is to be constrained using SPCs at the bolt locations, as shown in the following figure.
The constraints are organized into the load collector 'constraints'.
To perform a normal modes analysis, a real eigenvalue extraction (EIGRL) card needs to be
referenced in the subcase. The real eigenvalue extraction card is defined in HyperMesh as a load
collector with an EIGRL card image. This load collector should not contain any other loads.
4. Click Color and select a new color from the color palette.
5. For V2, enter the value 200.000.
6. For ND, enter the value 6.
2. Use the browser to select set the directory in which you want to save the file.
3. For the File name:, enter sshield_complete.hm.
4. Click Save.
Eigenvector results are output by default, from OptiStruct for a normal modes analysis. This
section describes how to view the results in HyperView.
Step 10: Load the Model and Result Files into the Animation Window
In this step, you will load a file into the HyperView animation window.
1. When the message 'ANALYSIS COMPLETED' is received in the HyperWorks Solver View
window, click Results. HyperView is launched and the results are loaded.
2. Click Close to exit the Message Log menu if it appears.
1. Click the animation selector switch in the lower toolbar and select Set Modal
Animation Mode .
6. Click Apply.
7. Under Undeformed shape:, set Show: to Wireframe. A deformed plot of the model
overlaid on the original undeformed mesh is displayed in the graphics area.
8. In the Results browser pull-down menu, you can change the view between various subcases
using the Load Case and Simulation Selection drop-down menus, as shown below:
Figure 1: Structural model with static loads and support constraints applied
Exercise
Step 1: Launch HyperMesh, set the OptiStruct User Profile and retrieve
the model file
1. Launch HyperMesh. The User Profiles dialog appears.
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.
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
4. Select the ie_carm.hm file you saved to your working directory from the optistruct.zip file.
Refer to Accessing the Model Files.
5. Click Open. The ie_carm.hm database is loaded into the current HyperMesh session,
replacing any existing data.
Step 2: Create load collectors used to conduct the inertia relief analysis
In this step you create two collectors: one for static loads and the other for constraints.
1. In the Model browser, right-click and select Create > Load Collector. A default load
collector template displays in the Entity Editor.
2. For Name, enter static_loads.
3. Click Color and select another color from the color palette. A new load collector named
static_loads is created.
4. Create another load collector named SPCs using steps 2.1 to 2.3.
Note: Use Next and Prev to browse through the different control
card pages.
2. Click TITLE and enter a title for this inertia relief analysis and click return.
3. Click PARAM and turn on INREL.
4. Under INREL_V1, toggle the selection to be -1. This requests that an inertia relief analysis
be performed.
5. Click return twice to go to the main menu.
ie_carm.html Web-based file which gives a summary of the input deck and
simulation conducted.
ie_carm.out ASCII based output file of the model check run before the
simulation begins and gives some basic information on the
results of the run.
ie_carm.stat Detailed breakdown on the CPU time used for each significant
stage in the analysis.
OptiStruct provides contour information for all of the loadsteps that were run. The following steps
describe the process for viewing those results in HyperView.
Step 10: View a von Mises stress contour of the static loadcase
1. Click the Contour icon on the toolbar.
2. Select Element Stresses (2D & 3D) as the Result type:.
3. The stress type should be set to vonMises.
4. Click Apply and notice the graphical display of stresses.
5. Once you are finished viewing, select File from the menu and select Exit to exit HyperView.
Overview
In this tutorial the steps required to perform a buckling analysis using OptiStruct are covered.
The figure below illustrates the structural model used for this tutorial.
This tutorial uses the following exercises to set up the structural model for a buckling analysis:
Create boundary conditions for buckling analysis
Post-process results
Step 1: Launch HyperMesh, set the OptiStruct User Profile and retrieve
the file
1. Launch HyperMesh. The User Profiles dialog appears.
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.
3. Click File > Open > Model. An Open Model browser window opens.
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
4. Select the buckling.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
5. Click Open.
The structural model has already been set up to contain the necessary elements, parts,
property, and material data.
buckling.out ASCII based output file of the model check run before the simulation
begins and gives some basic information on the results of the run.
buckling.stat Detailed breakdown on the CPU time used for each significant stage in
the analysis.
OptiStruct gives you contour information for all of the loadsteps that were run. This section
describes the process for viewing those results in HyperView.
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
4. Select the dissimilar.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
5. Click Open. The dissimilar.hm database is loaded into the current HyperMesh session,
replacing any existing data.
11. In the Model browser, right-click and select Create > Property. A default PSHELL property
template displays in the Entity Editor.
23. Click dest component = and select membrane_coarse from the list of components.
24. Click move. The elements are now part of the membrane_coarse component.
25. Click return.
26. On the Tool page, select faces.
27. Click delete faces.
28. Click return to exit the panel.
29. In the Model browser, right-click and select Create > Component.
30. For Name, enter membrane_fine.
31. For Property, click Unspecified > Property.
32. In the Select Property dialog, select membrane_fine and click OK.
33. On the Tool page, select faces.
34. Click comps, select the solid_fine component, and click select > find faces. Membrane
elements are created on the faces of solid_fine component and they appear on the graphics
window.
35. Click return.
36. Right-click ^faces component and click Isolate. This displays only ^faces component in the
graphics window.
37. From the Tool page, select the organize panel.
38. Select only the elements that lie on the matching face as shown below. (Use a method similar
to the one mentioned in Step 2.22 of this section).
39. Click dest component= and select membrane_fine from the list of components.
40. Click move. The elements are now part of the membrane_fine component.
41. Click return.
42. From the Tool page, select the faces panel.
43. Click delete faces and click return to return to the main menu.
44. In the Model browser, right-click and select Components > Hide.
45. Click the icon highlighted in red below to keep only the membrane elements in display.
11. Click Load Model and select the file control.h3d you saved to your working directory
from the optistruct.zip file as both the model and results file.
12. Click Apply.
Torsion frame
Step 1: Launch HyperMesh, set the OptiStruct User Profile and retrieve
the file
1. Launch HyperMesh.
2. Select OptiStruct from the User Profiles dialog and click OK.
3. Click File > Open. An Open Model browser window opens.
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
4. Select the frame.hm file you saved to your working directory from the optistruct.zip file.
Refer to Accessing the Model Files.
5. Click Open. The frame.hm database is loaded into the current HyperMesh session, replacing
any existing data.
Components in this model that have names starting with the word "Flange" represent junctions in
which different components are connected together. While reviewing, closely watch the flange
area formed by the Skin and Rib components (highlighted in the following figure). Review the
ply lay-up of the Skin_inner, Rib, Flange1_Rib_Skin, and Flange2_Rib_Skin components
(laminate layout is shown in the bottom portion of the following figure). Note that few plies are
common for the Skin_inner, Flange1_Rib_Skin, Flange2_Rib_Skin, and Skin_outer
components, but appear in different stacking sequence in each component. For example, the 4th
ply in Skin_inner is the 3rd ply in Flange2_Rib_Skin and the 2nd ply in Skin_outer
components.
1. From the 2D page, click HyperLaminate to enter the Graphic User Interface (GUI).
This opens the HyperLaminate (GUI) in which the ply lay-up information can be defined,
reviewed and edited. Material properties and design variables can also be created and edited
here.
2. Expand the Laminates portion of the tree structure on the left-hand side of the screen.
3. Select the Skin_inner PCOMP.
Details of the laminate appear in the GUI.
4. Verify that the lay-up definition for Skin_inner matches the first 5 entries of the table below,
which is the lay-up information of Flange1_Rib_Skin component.
5. Select the Rib PCOMP and verify that the 3rd and 4th lay-up definition for Rib matches the
6th and 7th entries in the following table.
6. Select the Flange1_Rib_Skin PCOMP to view the ply lay-up definitions. Verify that the lay-
up definition for Flange1_Rib_Skin matches the following table.
Observe that the first 5 layers are the same as Skin_inner lay-ups and that the last two lay-
ups are the same as the 3rd and 4th lay-up of Rib, as shown in the last figure. You can verify
how other flanges are modeled.
Laminate properties of Flange1_Rib_Skin:
Ply
Material Thickness T Orientation SOUT
Number
Ply
Material Thickness T Orientation SOUT
Number
7. You can also review the other components. Once the review is completed, select Exit from
the File menu. Exit the HyperLaminate GUI and return to HyperMesh.
8. Click Apply. This contours the maximum principle stress for the 3rd ply of all the components
in the model.
9. Click the Isometric View icon in the Standard Views toolbar to see the model, as shown
in the following figure.
The stress value does not vary gradually in the top face region, but suddenly decreases to a
lower value across the Flange2_Rib_Skin component. Looking at the table of laminate
properties of Flange1_Rib_Skin again, observe that the 3rd ply property of the
Flange2_Rib_Skin component is of a matrix material and the third plies in the components
adjacent to it (Flange1_Rib_Skin and Skin_outer) are of a carbon fiber material. The sudden
changes in the stress values occur because we are looking at stress on two different
materials. This example shows that, for the results to be meaningful during post-processing of
the PCOMP results, you have to correlate the ply results to their corresponding ply property.
This highlights that, during the post-processing of PCOMP components, plotting results based
on just the ply number is not sufficient. You have to keep track of ply properties (material,
thickness, orientation, failure index, etc.) on your own during post-processing with this
method. In cases that use large and complex models, it becomes tedious to track the
individual ply properties during post-processing.
This drawback to using PCOMP can be avoided with the use of the PCOMPG card for property
definition. Using the PCOMPG card, you can assign a global ply number for each ply and post-
process the results based on global ply number. The following steps explain the procedure to
redefine the model with PCOMPG property.
2. From the 2D page, select the HyperLaminate panel. This opens the HyperLaminate GUI in
which the ply lay-up information can be defined, reviewed and edited.
Now you create new PCOMPG components with global ply numbers defined as shown in the
above figure. As discussed earlier, the 4th ply in Skin_inner is the 3rd ply in Flange2_Rib_Skin
and the 2nd ply in Skin_outer components. Therefore, all of these plies will be defined with
the same global ply ID 4. Similarly, all other plies are to be defined, as shown in the above
figure.
6. Expand the laminates portion of the tree structure on the left-hand side of the screen.
7. Right-click PCOMPG. A menu appears.
8. Click New. This creates new component, which is named NewLaminate1 by default, and the
tree structure is expanded.
9. Rename the component to Skin_inner_GPLY by right-clicking and select Rename in the text
field and overwrite the default component name.
10. In the Add/Update plies: section under the field GPLYID, enter 1.
11. Select the pull-down menu below Material and select carbon_fiber.
12. Below the Thickness T1 field, enter 1.2.
13. Below the Orientation field, enter 45.
14. Select the pull-down menu below SOUT and select YES.
15. Click Add New Ply to add the ply information.
16. Repeat this procedure to add 4 more plies with the properties shown in the table:
GPLY
Material Thickness T Orientation SOUT
ID
Referring to the figure showing laminate information with global ply number above, you will
create the Flange1_Rib_Skin_GPLY component.
19. Right-click Skin_inner_GPLY and select Duplicate from the menu to create an identical
component.
20. Rename the component as Flange1_Rib_Skin_GPLY by right-clicking and select rename in
the text field and overwrite the component name.
21. Add 2 more plies with the properties shown in the following table using the Add New Ply
feature.
GPLY
Material Thickness T Orientation SOUT
ID
The new component Flange1_Rib_Skin_GPLY was created. Its first 5 plies are the same as
Skin_inner_GPLY and its last 2 plies are the 3rd and 4th plies of the Rib component.
To reduce the number of steps in this tutorial, the ply lay-up information of other components
is already defined with PCOMPG property and appropriate laminate information in the
updated_PCOMPG_properties.fem file you saved to your working directory from the
optistruct.zip file. This file is imported into HyperMesh to update (overwrite) the properties
instead of manually updating them.
The updated_PCOMPG_properties.fem file is saved in OptiStruct input file format. Open this
in any text editor to review how the components are defined with PCOMPG properties. A
section of the file is shown below.
22. Click File > Exit. This will exit the HyperLaminate GUI and return to HyperMesh.
23. Click File > Import > Solver Deck.
24. Toggle and expand the Import options and check the box next to FE overwrite.
This option overwrites the old PCOMP properties with PCOMPG properties defined in the
updated_PCOMPG_properties.fem file.
25. Click on the folder icon next to File: and select the updated_PCOMPG_properties.fem file
and click Import.
26. Click Close.
Post-processing the results based on global ply number eliminates the need to track the ply
number and corresponding ply properties on the components. The results are displayed based on
the global ply number, irrespective of the ply order, so you can chose any one global ply number
and view results across the whole component. If a particular ply is not present in any given
region, no result is displayed.
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
2. Select the axi-symmetry_full_geometry.hm file you saved to your working directory from
the optistruct.zip file. Refer to Accessing the Model Files.
3. Click Open.
The axi-symmetry_full_geometry.hm database is loaded into the current HyperMesh
session, replacing any existing data. The database only contains geometric data.
You will find that the structural model has already been set up with the necessary elements,
boundary conditions, property, and material data so that it is ready to solve. Pressure load is
applied on the top face of the geometry and constraints are defined at the bottom face. Note that
the model is symmetrical about the z-axis and that loads and boundary conditions are
symmetrical about the same axis as well. These represent the conditions necessary for modeling
axi-symmetry problems. First, obtain the result for the full model and then you model a small
part of the model with boundary conditions suitable to enforce the axi-symmetric behavior.
Finally, you compare the results of the axi-symmetric model with the full model results.
Solving the job and post-processing the results is explained in the next sectioin.
6. Click the Load model icon on the toolbar and load the axi-symmetry_model.h3d. This
loads the complete path of the selected .h3d file in the field. Also note that the same file path
is loaded next to the field Load results.
7. Click Apply.
8. Click XZ Left Plane View to display the Left view.
9. Click the Contour icon on the toolbar and contour the displacements.
10. Compare the displacement results of the axi-symmetry model with the result from the full
model. The results should match, as shown in the below picture. Similarly, stress and other
results will also match.
Step 1: Load the OptiStruct User Profile and import the model
1. Launch HyperMesh.
2. Select OptiStruct in the User Profile dialog 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.
3. Click the Import Solver Deck panel toolbar icon . An Import panel opens. The File type:
is OptiStruct.
4. Click the open file icon in the File: field. A Select OptiStruct file browser window opens.
5. Select the pipe.fem file you saved to your working directory from the optistruct.zip file.
Refer to Accessing the Model Files.
6. Click Open. The location of pipe.fem displays in the file: field.
7. Click Import. The pipe.fem database is loaded into the current HyperMesh session.
8. Click Close.
5. In the Model browser, right-click and select Create > Property. A default PSHELL property
displays in the Entity Editor.
6. For Name, enter solid.
7. For Material, click Unspecified > Material.
8. In the Select Material dialog, select steel and click OK.
9. For Card Image, select PSOLID from the drop-down menu and click Yes to confirm.
The property of the solid steel pipe has been created as 3D PSOLID. Material information is
linked to this property.
A structural constraint spc_struct is applied on the RBE2 element to fix the pipe on the ground.
Two empty load collectors, spc_heat and heat_flux have been pre-created. In this section, the
thermal boundary conditions and heat flux are applied on the model and saved in spc_heat and
heat_flux, respectively.
7. Click the Table icon to the right of Data: Comments and enter the following text in the
first row of the pop-out table.
TEMP=1
Note the TEMP ID used above could be different from your model. Make sure the ID of the
heat transfer subcase is selected for TEMP.
8. Click Close.
Gradient temperatures and flux contour results for the steady state heat conduction analysis and
the stress and displacement results for the structural analysis are computed from OptiStruct.
HyperView is used to post-process the results.
5. In the Contour panel, select the first pull-down menu below Result type: and select Element
Fluxes (V).
6. Click Apply. A contoured image representing thermal fluxes should be visible.
7. Select the first pull-down menu below Result type: and select Grid Temperatures (s).
8. Click Apply. Both flux and temperature results are shown below.
Step 1: Load the OptiStruct User Profile and import the model
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog 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.
3. Click the Import Solver Deck panel toolbar icon . An Import window opens. The File
type: is OptiStruct.
4. In the File: field, click the open file icon . A Select OptiStruct file browser window opens.
5. Select the thermal.fem file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
6. Click Open. The location of thermal.fem displays in the File: field.
7. Click Import > Close. The thermal.fem database is loaded into the current HyperMesh
session.
A new material, steel, is created with both structural and thermal properties.
5. In the Model browser, right-click and select Create > Property. A default PSHELL property
displays in the Entity Editor.
6. For Name, enter solid.
7. For Card Image, select PSOLID and click Yes to confirm.
8. For Material, click Unspecified > Material.
9. In the Select Material dialog, select steel and click OK.
The property of the solid steel pipe has been created as 3D PSOLID. Material information is
linked to this property.
In this exercise the thermal boundary conditions are applied on the model and saved in a
predefined load collector spc_temp. A predefined node 4679 specifies the ambient temperature. A
predefined node set node_temp contains the nodes on the inside surface of the pipe.
9. Click on the switch button beside elems and select faces from the list.
10. Click the highlighted solid elems and select by sets from the selection menu.
13. Select 4 nodes on the surface face of a solid element, as shown in the following figure.
Figure 4: Selected surface nodes on the solid element outside the pipe.
9. Activate the OPTION fields for both outputs and select ALL.
The FORMAT and OUTPUT fields for THERMAL output may open up a new window. Click on the
first field in the window to select the corresponding values.
FLUX and THERMAL output can also be requested in the Control cards panel on the Analysis
page.
Gradient temperatures and flux contour results for the steady state heat conduction analysis and
the stress and displacement results for the structural analysis are computed from OptiStruct.
HyperView will be used to post process the results.
Figure 1: Extended surface heat transfer fin for convective and conductive transient heat transfer.
The extended surface heat transfer fin shown in Figure 1 is meshed with CHEXA elements in
HyperMesh and a transient heat transfer analysis is performed in HyperMesh using the Altair
OptiStruct solver. A typical heat flux load of 100 KW/m 2 is applied to the face connected to the
outer surface of the system. An ambient temperature of 25°C is assumed and all material
properties are assumed to remain constant with temperature and time. Free (Natural) convection
is assumed over the entire surface of the material, wherein heat transfer between the surface of
the fin and the surrounding air occurs due to a complex mechanism of density differences as a
result of temperature gradients.
Tutorial Prerequisites
1. The latest version of Altair HyperMesh, Altair HyperView and Altair OptiStruct software
installations. Transient heat transfer analysis is available only in Altair HyperMesh version-
12.0.110, Altair HyperView version-12.0.110 and Altair OptiStruct version-12.0.202 and later.
2. The heat_transfer_fin.fem solver deck is available from the optistruct.zip file. Refer to
Accessing the Model Files.
Figure 2: Heat exchanger fin model for Transient Heat Transfer Analysis.
(1)
Where,
[C] is the heat capacity matrix
[K] is the conductivity matrix
[H] is the boundary convection matrix due to free convection
Checkpoint
Steady-state heat transfer analysis, generally, is sufficient for a wide variety of applications.
However, in situations where the system properties vary significantly over time the transient
nature of heat transfer must be considered. Some examples are the relatively slow heating up of
airplane gas turbine compressor disks compared to the turbine casing leading to aerodynamic
issues during takeoff or the analysis of the time taken for the onset of frostbite in fingers or toes.
The following exercises are included:
Creating the thermal material and the solid property for the given component
Assigning the material and property to the component
Creating flux and convective loads and boundary conditions for the model
Submitting the job to Altair OptiStruct
Post-processing the results using Altair HyperView
Step 1: Loading the OptiStruct User Profile and importing the model into
HyperMesh
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. Click OK. This loads the OptiStruct 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.
6. Click the open file icon in the File: field. A Select OptiStruct file browser window opens.
7. Select the heat_transfer_fin.fem file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
8. Click Open. The location of heat_transfer_fin.fem displays in the File: field.
9. Click Import. The heat_transfer_fin.fem database is loaded into the current HyperMesh
session.
Since you are conducting a purely heat transfer analysis, structural isotropic properties (for
example, MAT1 card) are not required. Also, it is assumed that the thermal material
properties (MAT4) are temperature independent.
A new material, steel, is created with thermal properties necessary for a transient heat
transfer analysis.
Now, create the solid property for this model referencing the PSOLID entry and connect the
material, steel, to this property; the property can then be assigned to the existing
component.
5. In the Model browser, right-click and select Create > Property. A default PSHELL property
displays in the Entity Editor.
6. For Name, enter solid.
7. For Card Image, select PSOLID and click Yes to confirm.
8. For Material, click Unspecified > Material.
9. In the Select Material dialog, select steel and click OK.
5. Click and enter the number of time steps (N) = 50 and set each time increment (DT) to
10.
This encompasses a total time period of 500 seconds in which to capture the behavior of the
system.
6. Click Close.
Step 5: Creating initial conditions for the transient heat transfer analysis
Since the temperature profile of the system varies over time, the initial grid point temperature
profile must be set to specify the starting point for the analysis. You assume that the
temperature of the entire system is equal to 25°C at T=0 seconds, the TEMPD bulk data entry
sets the initial temperatures.
1. In the Model browser, right-click and select Create > Load Collector.
2. For Name, enter Initial Conditions.
3. For Card Image, select TEMPD.
4. For T1, enter a value of 25.
Ambient temperature thermal boundary conditions is applied on the model by creating specific
load collectors for each. The ambient temperature is controlled using an SPCD entry, as this will
allow an ambient temperature variation over time to help mimic such physical requirements (if
any).
6. Next, create the amplitude (constant part) of the time variant ambient temperature using an
SPCD data entry. Click BCs > Create > Constraints.
10. Click create/edit and enter 25.0 in the D field on the SPCD data entry. This creates an SPCD
referencing the ambient node specifying a temperature of 25°C.
11. Click return twice to go back to the Analysis page.
12. Next, create another load collector to define the time variant nature of the ambient
temperature. This is done by specifying a TABLED1 entry referenced by the previously created
TLOAD1 entry.
13. Create a new load collector (refer to Step 6.1) named Ambient SPCD Table, and set Card
Image, as TABLED1.
14. For TABLED1_NUM, enter 2 and press ENTER.
15. In the table, enter x(1) = 0.0, y(1) = 1.0, x(2) = 500.0, and y(2) = 1.0.
16. Click Close.
Checkpoint
The SPCD and its corresponding table are linked to the previously created TLOAD1 entry.
17. In the Model browser, click the Ambient SPCD TLOAD1 load collector. It is displayed in the
Entity Editor.
18. For EXCITEID, select the Ambient SPCD load collector from the menu.
19. For TYPE, select DISP, then click TID and select the Ambient SPCD Table load collector
menu.
All entities referenced by SPCD entries should also be constrained by SPC data entries. The
value of the corresponding SPC referencing an ambient point controlled via an SPCD by
TLOAD1/2 entries should be equal to zero (0.0).
20. Create a new load collector (refer to Step 6.1) named Ambient SPC and for Card Image,
select None.
21. Make sure that the newly created load collector Ambient SPC is current and click BCs >
Create > Constraints and click nodes in the Constraints panel. Select by id from the
extended menu, enter 5672 and press ENTER. The ambient node is highlighted in the
graphics window above the structure.
22. Enter 5.0 in the size= field and uncheck the boxes beside all the degrees of freedom (dof1
through dof6) and enter 0.0 in all the fields next to the dof#.
23. For load types =, select SPC.
24. Click create/edit and enter 0.0 in the D field on the SPC data entry. This creates an SPC
referencing the same ambient node that is controlled by the SPCD data entry.
25. Click return twice to go back to the Analysis page.
Ambient temperature thermal boundary conditions have been assigned to the model and heat
flux load from the outer surface of the engine (to which the fin is attached) is applied on the
model. A time-varying heat flux load of 0 to 0.1 W/mm2 from 0 to 500 seconds is used for the
analysis of this fin. This load is applied on the model by creating specific load collectors for the
corresponding TLOAD1, QBDY1 and TABLED1 entries similar to the procedure used for the
ambient temperature SPCD definition.
4. An interface is now created between the heat flux source and the solid elements on the
surface of the fin. This is done by clicking BCs > Create > Interfaces and specifying
conduction_interface in the Name field of the Create group dialog.
5. For Card Image, select CONDUCTION from the drop-down menu and click Yes to confirm.
6. For Slave Entity IDs, click on the yellow Elements panel. A panel appears under the
graphics area.
7. Click on the switch button beside elems and select faces from the list.
8. Click the highlighted solid elems and select by sets from the selection menu.
9. Check the box next to Element_set_Flux and click select. The predefined element set is
now highlighted in white on the model.
Tip: The break angle helps find adjacent solid faces for the same element
set, however, since this surface element set generation requires only
one face, the value of the break angle is not germane in this situation.
10. Click nodes and select the nodes in the Figure 11.
Figure 11: Selecting the nodes on the highlighted surface for conduction surface element creation.
Figure 12: Newly generated surface elements are displayed in blue as shown above.
14. Next, create the amplitude (constant part) of the time variant heat flux using a QBDY1 data
entry. Do this by clicking on BCs > Create > Flux.
15. Click elems, select by group and select conduction_interface. The newly created surface
elements are highlighted in white on the model.
21. Click next to Data. In the pop-out window, enter x(1) = 0.0, y(1) = 0.0, x(2) = 500.0
and y(2) = 1.0.
Checkpoint
The QBDY1 flux load and its corresponding table are linked to the previously created TLOAD1
entry.
22. In the Model browser, click the Heat Flux TLOAD1 load collector. The entry is displayed in
the Entity Editor.
Free convection is assigned in a similar manner to the procedure used for the creation of the
conduction interface. Free convection is, however, automatically assigned to all heat transfer
subcases and the PCONV and CONV entries should refer to the material, steel, and the ambient
temperature. The ambient temperature calculates the amount of heat transferred through free
convection.
7. Click nodes next to face nodes:, this highlights the same element set in black. Click any four
nodes on all highlighted faces of the model as shown in Figures 17 and 18 and specify a break
angle of 89°.
Tip: The break angle helps find adjacent solid faces for the same element set, all
adjacent faces with the angle between surface normals less than the
specified break angle are selected for surface element creation.
Figure 17: Selecting the nodes on four of the seven highlighted surfaces for convection surface element creation.
The newly created CHBDYE surface elements are displayed in yellow, as shown in Figure 19
below.
Figure 19: Newly generated CHBBDYE surface elements are displayed in yellow on the model.
Step 9: Combining the two TLOAD1 entries into one DLOAD entry.
Two different TLOAD1 entries have been defined and since they are to be referenced in the same
subcase they should be combined using a DLOAD data entry.
1. Create a new load collector (refer to Step 6.1) named Combined Flux and Convection and
select DLOAD as the Card Image.
2. For S, enter 1.0.
3. As only a simple linear addition of the two TLOAD1 entries are required, for DLOAD_NUM,
enter 2 and press ENTER.
4. Click next to Data below the DLOAD_NUM field. In the DLOAD_NUM pop-up window,
enter S(1) = 1.0 and S(2) = 1.0.
5. For L(1), select Ambient SPCD TLOAD1 from the menu and for L(2), select Heat Flux
TLOAD1 from the menu.
Checkpoint
The DLOAD entry is created as a linear combination of two TLOAD1 entries – Heat Flux
TLOAD1 and Ambient SPCD TLOAD1.
6. Click Close.
Grid temperatures, element temperature gradients and flux contour results are computed for the
transient heat transfer analysis and HyperView are used to post-process the results.
Step 12: Viewing the results of the transient heat transfer analysis
1. When the message Process completed successfully is received in the command window, 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. Click Close to close the message window, if one appears.
3. Click the Contour toolbar icon .
4. Select the first pull-down menu below Result type: and select Grid Temperatures(s).
5. Click Apply, select Time = 5.0000000E+02 from the Results browser. A contour plot of
grid temperatures at the final time step is created as shown in Figure 22.
Checkpoint
In Figure 22, this is the grid point temperature plot after 500 seconds. The system is input a
linearly increasing heat flux from 0 to 0.1 W/mm 2 from 0 to 500 seconds respectively. Therefore,
a physical correlation can be the effect of starting an IC engine to full capacity wherein the flux
transmitted to the outer surface linearly increases with time. Note that the flux patterns in
actuality may be different and may fluctuate based on the duration of the power cycles. The
maximum temperature of 81.3°C predictably occurs at the elements closest to the heat flux
loading site and the minimum temperature of 29.5°C occurs at elements farthest from the heat
source.
6. Click Apply, select Time = 2.0000000E+01 from the Results browser. A contour plot of grid
temperatures is created, as shown in Figure 23.
Figure 23: Grid temperature contour plot after 20 seconds – WITH FREE CONVECTION.
Figure 24: Grid Temperature Contour Plot after 500 seconds - WITHOUT FREE CONVECTION
The maximum temperature at the outer surface of the heat source system is 125.3°C which is an
increase of 44°C in 500 seconds. Therefore using an extended surface fin is a very effective way
to reduce the temperature of a system.
Printed Circuit Boards (PCB’s) are used in electronic components to both mechanically support
and provide electrical connections between components. Construction involves etching a thin
copper layer that has been deposited onto a non-conductive, glass-fiber/epoxy composite
substrate. Electrical components are then mounted to the board and connected to the copper
traces with electrical solder.
The concentrated, intense heating that occurs during the soldering process creates stresses in
the substrate material. In this exercise, you will simulate this process and determine if the
stresses and strains resulting from this process are acceptable or not.
The model makes use of solid hexahedral (CHEXA8) elements with a thin skin of shell elements
(CQUAD4) on the outside faces.
The consistent unit system used in this simulation are: kg, mm, GPa, kN and °C
Problem Statement
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
2. Select the circuit_board.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
3. Click Open. The circuit_board.hm database is loaded into the current HyperMesh session,
replacing any existing data. The database only contains geometric data.
5. Enter the following values for the thermal expansion rates and reference temperature:
A1 A2 A3 TREF
7. Repeat steps 1 to 6 to create another property with name Solids, with Card Image set as
PSOLID and Material as PCB_solids.
8. In the Model browser, click the pcb_solids component. The component entry is displayed in
the Entity Editor below.
9. For Property, click Unspecified > Property.
10. In the Select Property dialog, select Solids and click OK to complete the property selection.
11. Repeat steps 8 to 10 for both solder_pads and shell_faces selecting shell for the property
name.
5. Click BCs > Create > Constraints to open the Constraints panel.
6. Click nodes > by sets.
7. Select the constrain_nodes entity set and click select.
8. Leave all 6 degrees of freedom selected and click create.
9. Click return to go back to the main menu.
Step 9: Add OUTPUT, STRAIN and SCREEN control cards to the analysis
deck
1. Click Setup > Create > Control Cards to open the Control Cards panel.
2. Click next to advance until OUTPUT is available, click OUTPUT to add card requesting output
results format.
3. For the number_of_outputs field on the lower part of the panel, enter 2.
4. Set one of the KEYWORD to OP2 to request the OP2 format results file, and set the second
output as H3D format. The frequency (FREQ) of the output can be set as ALL.
5. Click return to go back to the Control Cards panel.
6. Click next to advance to the second page of control cards, then once more to go to the third
page.
Problem Statement
The channel_brkt_modal.hm file is used for this tutorial. It contains the bracket and channel
assembly pictured below.
To complete the setup of the model for a modal analysis with OptiStruct, you need to define a
normal modes SUBCASE, containing METHOD and SPC statements.
Exercise
Step 1: Select User Profile and retrieve the HyperMesh model file
1. Launch HyperMesh. The User Profiles dialogs appears.
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.
3. Click File > Open. An Open Model browser window opens.
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
4. Select the channel_brkt_modal.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
4. Click BCs > Create > Constraints to open the Constraints panel.
5. Switch the entity selector to lines.
6. Select the six lines on the perimeter of the channel’s bottom surface.
To view the selected lines clearly, switch to the Transparent Elements mode, as shown
below:
Step 7: Define the load step (OptiStruct loadsteps) for the modal load
case
Use the Load Step Entity Editor in this step. Define the loadstep to contain the load collectors
constraints and modal.
1. In the Model browser, right-click and select Create > Load Step.
2. For Name, enter normal_modes.
3. For Analysis type, select Normal modes.
4. For METHOD(STRUCT), select modal.
5. For SPC, select the load collector constraints.
The following file found in the optistruct.zip file is needed to perform this tutorial. Refer to
Accessing the Model Files.
direct_response_flat_plate_input.fem Original ASCII OptiStruct input deck
Exercise
7. Uncheck all the dof's; except dof3, and click the = to the right of dof3 and enter a value of
20.
8. Click load types= and verify that DAREA is selected from the extended entity selection
menu.
9. Click create, and then click return.
6. Click the Table icon below TABLED1_NUM and enter x(1) = 0.0, y(1) = 1.0, x(2) =
1000.0 and y(2) = 1.0 in the pop-out window.
7. Click Close. This provides a frequency range of 0.0 to 1000.0 with a constant 1.0 over this
range.
6. Click and enter F1= 20.0, DF = 20.0, NDF = 49, in the pop-out window.
7. Click Close. This provides a set of frequencies beginning with 20.0, incremented by 20.0 and
49 frequencies increments.
8. Select the Next Page icon again to display page 3 containing Subcase 1 (subcase1) -
Displacement of grid 19 (see figure).
Exercise
2. Click the Open file icon in the File: field. A Select OptiStruct File file browser window
opens.
3. Select the modal_response_flat_plate_input.fem file you saved to your working directory
from the optistruct.zip file and click Open. Refer to Accessing the Model Files.
The location of the file displays in the File: field.
4. Click Import to import the model file.
5. Click Close to close the Import tab.
5. In the Model browser, right-click and select Create > Load Collector.
6. For Name, enter unit-load.
7. Click Color and select a color from the color palette. A new load collector, unit-load is
created.
8. Click the Display Numbers icon to open the Numbers panel.
9. Click nodes > displayed.
10. Check the box next to display.
11. Select the green on button. All of the node numbers on the flat plate should now be
displayed.
4. Uncheck all dofs except dof3, and click the = to the right of dof3 and enter a value of 1.
6. Click below TABLED1_NUM and enter x(1) = 0.0, y(1) = 1.0, x(2) = 1000.0 and y(2) =
1.0 in the pop-out window.
7. Click Close. This gives a frequency range of 0.0 to 1000.0 with a constant 1.0 over this
range.
6. Click and enter F1= 20.0, DF = 20.0, NDF = 49, in the pop-out window.
7. Click Close. This gives a set of frequencies beginning with 20.0, incremented by 20.0 and 49
frequencies increments.
7. Click the Axis toolbar icon . Select Logarithmic option and use the parameters shown
below (see figure) to make logarithmic plots of the results.
There are two sets of results on this page. The top graph shows Phase Angle verses
Frequency (log). The bottom graph shows Magnitude verses Frequency (log) (see figure) for
Displacement of grid 15.
8. Directly underneath the blue graph border, click the Next Page icon . Page 2 displays,
which shows Subcase 1 (subcase1) - Displacement of grid 17.
The bracket is constrained at the bottom of the two legs. Transient dynamic loads are to be
applied at the grid points of the top, flat surface of the bracket around the hole in the negative z
direction. The time history of the loading is shown in the next figure. The direct transient analysis
is run for a total time of 4 seconds with the time being divided into 800 increments (that is time
step is 0.005). Structural damping has been considered for the model. A concentrated mass
element is defined at the center of the spider and z displacements are monitored at the
concentrated mass at the center of this hole.
This tutorial uses the following exercises to set up a direct transient dynamic analysis:
Create the time dependent dynamic load or the variation of load vs time
Create the time step for transient analysis
Create the grid point forces on the top flat surface of the bracket
Create the transient response dynamic excitation
Create the subcase to include all the necessary loads as defined above
Specify structural damping and output requests
Run direct transient dynamic analysis
Post-process results using Altair HyperGraph
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
2. Select the bracket_transient.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
6. Click the table icon below TABLED1_NUM and enter the values in the pop-out window, as
shown in the figure below.
Step 4: Create a TSTEP (Transient time step to define the time step
intervals at which solution is generated and output)
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 on the pop-up menu.
5. For TSTEP_NUM, enter 1 and press ENTER.
6. For N, enter the number of time steps as 800.
7. For DT, enter the time increment of 0.005. The total time applied to the load is: 800 x 0.005
= 4 seconds.
NO has a default value of 1.0. This is the time step at which output is requested.
8. Click Close.
4. Click Yes to close the message window. This file automatically creates plots for the
displacement results contained in the file.
Since the loading is applied only in the z-direction, you are interested in the z-displacement
time history of node 395.
5. Click on the Curve Attributes toolbar icon and turn off the curves X Trans and Y Trans.
This can be done by selecting the individual curves (X Trans and Y Trans) and then by clicking
the line attributes Off, as shown below:
Z-displacement time history of the concentrated mass at center of spider for direct transient dynamic analysis
As can be observed from the above image, the displacements of node 395 are in the negative
z-direction as the loading is in the –z direction too. The displacements eventually damp out
due to the structural damping present in the model.
The bracket is constrained at the bottom of the two legs. Transient dynamic loads are to be
applied at the grid points of the top, flat surface of the bracket around the hole in the negative z-
direction. The time history of the loading is shown in the next figure. The modal transient
analysis is run for a total time of 4 seconds with the time being divided into 800 increments (that
is time step is 0.005). Modal damping has been defined as 2% critical damping for all the modes.
Modes up to 1000 Hz have been considered. A concentrated mass element is defined at the
center of the spider and z-displacements are monitored at the concentrated mass at the center of
this hole.
This tutorial uses the following exercises to set up a modal transient dynamic analysis:
Define the time dependent dynamic load or the variation of load vs time
Define the time step for transient analysis
Define the grid point forces on the top flat surface of the bracket
Define modal damping table
Define load collector to extract normal modes up to 1000 Hz using the Lanczos method
Define the transient response dynamic excitation
Define the subcase to include all the necessary loads as defined above
Specify output requests
Run modal transient dynamic analysis
Post-process results using Altair HyperGraph
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
2. Select the bracket_transient.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
3. Click Open. The bracket_transient.hm database is loaded into the current HyperMesh
session, replacing any existing data.
6. Click the Table icon below TABLED1_NUM and enter the values in the pop-out window, as
shown in the figure below.
Step 4: Create a TSTEP card (the transient time step to define the time
step intervals at which solution is generated and output)
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 on the pop-up menu.
5. For TSTEP_NUM, enter the number 1 and press ENTER.
6. In the dialog, under N, enter the number of time steps as 800.
7. To specify the time increment, enter 0.005 under DT. The total time applied to the load is:
800 x 0.005 = 4 seconds. NO has a default value of 1.0. This is the time step at which output
is requested.
8. Click Close to go back to the Entity Editor.
Step 5: Create DAREA card to define forces on the top surface of the
bracket
1. In the Model browser, right-click and select Create > Load Collector.
2. For Name, enter darea.
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 >> by sets from the pop-up menu. Two sets are displayed.
7. Select force and click select. The nodes that belong to the set force are selected.
8. Unselect (right-click) all degrees of freedom (dof); except dof3 indicating that dof3 is the
only active degree of freedom.
9. For dof3, enter a value of -1500.
10. Set load types = to DAREA.
11. Click create. This creates a force of 1500 units applied to the selected nodes in the negative
z-direction.
12. Click return.
6. Click below TABDMP1_NUM and enter the values in the pop-out window, as shown in the
figure below.
7. Populate the frequency and damping values for frequencies 0 and 1000 Hz and damping to be
0.02, as shown below. This provides a table of damping values for the frequency range of
interest.
Step 9: Create the load step to perform the modal transient dynamic
analysis
1. In the Model browser, right-click and select Create > Load Step.
2. For Name, enter transient.
3. For Analysis type, select Transient (modal) from the drop-down menu.
4. For SPC, select spc from the Select Loadcol pop-out window.
5. For DLOAD, select tload1.
6. For TSTEP(TIME), select tstep.
7. For METHOD (STRUCT), select the load collector eigrl created previously.
8. For SDAMPING (STRUCT, select the load collector tabdmp1 created previously. A subcase
is created that specifies the loads, boundary conditions, and damping for modal transient
dynamic analysis.
Z-displacement time history of the concentrated mass at center of spider for direct transient dynamic analysis
As can be observed from the above image, the displacements of node 395 are in the negative z-
direction as the loading is in the –z direction too. The displacements eventually damp out due to
the structural damping present in the model.
There are four shell components in the model: the mounting flange, the web, the top and bottom
flanges, and the lug. Gap elements have already been defined in the model and they connect the
web to the lug. Coupling forces are applied to the lug and pressure loading has been defined on
the top and bottom flanges of the rib joint. The mounting flange is constrained in all degrees of
freedom at the four mounting hole locations and the lug is constrained for the z-displacements
and rotations to prevent rigid body motion.
The following exercises are included:
Step 1: Load the User Profile and Retrieve the Model File
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog. 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.
3. Click File > Open. An Open Model browser window opens.
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
4. Select the rib.hm file you saved to your working directory from the optistruct.zip file.
Refer to Accessing the Model Files.
5. Click Open. The rib.hm database is loaded into the current HyperMesh session, replacing any
existing data. The database only contains geometric data.
The location of rib.hm displays in the file: field.
7. Click Mesh > Create > 1D Elements > Gaps to open the Gaps panel.
8. Select the update subpanel.
9. Click elems >> by collector.
10. Select gap by checking the box beside it.
11. Click the green select button.
12. Click property= and click on gap_prop.
8. Click the XY Top Plane View icon to display a top view of the Web.
9. Click Apply. This should show the contour of stresses on the Web component under the
coupled loading.
Analysis Review
Even though the deformation patterns are similar for both linear and nonlinear analyses, the
stress patterns differ. Though the horizontal loads are in opposing directions in the lug, the stress
distribution in the web for the linear run are the same around both the lug holes which is not
correct. This happens as all the gaps are in a closed condition for the linear analysis. Nonlinear
gap analysis gives more accurate representation. The gap status, open or closed, depending on
loading condition can also be observed from the .out file (shown below):
ITERATION 0
NONLINEAR ITERATION SUMMARY Subcase 1
Exercise
Step 1: Launch HyperMesh, set the OptiStruct User Profile and retrieve
the structural model
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. 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.
4. Click File > Open. An Open Model browser window opens.
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
5. Select the direct_psd.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
6. Click Open. The direct_psd.hm database is loaded into the current HyperMesh session,
replacing any existing data. The database only contains model data.
6. Click Close.
Step 3: Add the RANDOM subcase information entry and output request
The RANDOM subcase information entry needs to be added to the frequency analysis model and
the output commands for RMS and PSD results will be added as well.
1. From the Analysis page, click control cards.
2. Go to GLOBAL_CASE_CONTROL panel.
3. Check the box in front of RANDOM and double-click the highlighted ID button and select
randps.
4. Return to Control Cards and click GLOBAL_OUTPUT_REQUEST. Check the box for
STRESS to activate the card edit panel.
5. Select OUTPUT2 as the FORMAT, PSDF under RANDOM, and YES under OPTION. RMS
and PSDF stress are output to a .op2 file.
6. Click return to go back to the Control Cards panel.
7. Select CASE_UNSUPPORTED_CARDS and add the following cards:
XYPLOT,DISP,PSDF/ 19(T3)
OptiStruct will output the PSDF for the translational displacement in z direction at node 19.
8. Click OK, and then click return.
Also, the following files will be output and which are specific to the random response analysis.
direct_psd.peak ASCII result file, containing RMS and peak values of PSD.
direct_psd_rand.mvw HyperView script file. This file will automatically create the
plot of PSD over the frequency for the results contained in
.rand file.
The half car model is excited at the bottom of the car, as shown by a red constraint symbol
(triangle) in the above figure. The excitation provided is with the application of a unit load along
the direction of the height of the car (Z-axis).
Process Overview
Step 1: Launch HyperMesh, set the OptiStruct User Profile and retrieve
the file
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. 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.
4. Click File > Open. An Open file browser opens.
Note: If HyperMesh Desktop was launched, use: File > Open >
Model.
5. Select the Half_Car.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
6. Click Open.
Step 2: Create isotropic and fluid materials and properties and assign
them to their structural and fluid elements
1. In the Model browser, right-click and select Create > Material.
2. For Name, enter MAT1_shells.
3. For Card Image, select MAT1 from the drop-down menu.
4. Fill in the fields for E, Nu and Rho respectively as 2.1e04, 0.33 and 8.0e-10.
6. Click the Table icon next to Data field and enter x(1) = 0.0, y(1) = 1.0, x(2) = 200.0
and y(2) = 1.0 in the pop-out window.
This provides a frequency range of 0.0 to 200.0 with a constant 1.0 over this range and it
would appear, as shown below on the GUI.
7. Click next to the Data field and enter, F1= 0.0, DF= 1.0,and NDF= 200.
This provides a set of frequencies beginning with 0.0, incremented by 1.0 and 200
frequencies increments and the card appears as shown below on the GUI.
Half_car.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 all
optimization iterations, and compute time information. Review this file
for warnings and errors.
Half_car.stat Summary of analysis process, providing CPU information for each step
during analysis process.
There are two sets of results on this page. The top graph shows Phase Angle verses
Frequency (log). The bottom graph shows Magnitude verses Frequency (log) (see figure
below) for Pressure at grid 18881.
Figure 1: Low Cycle and High Cycle regions on the S-N curve
Since S-N theory deals with uniaxial stress, the stress components need to be resolved into one
combined value for each calculation point, at each time step, and then used as equivalent
nominal stress applied on the S-N curve.
In OptiStruct, various stress combination types are available with the default being “Absolute
maximum principle stress”. In general “Absolute maximum principle stress” is recommended for
brittle materials, while “Signed von Mises stress” is recommended for ductile material. The sign
on the signed parameters is taken from the sign of the Maximum Absolute Principal value.
A flowchart of the fatigue setup in HyperMesh can be described as shown in the image below.
The three aspects to the fatigue definition are the fatigue material properties, the fatigue
parameters and the loading sequence and event definitions.
The fatigue material properties: (S-N curve)
FATPARM: Defines the parameters for the fatigue analysis. These include stress combination
method, mean stress correction method (Figure 4), Rainflow parameters, Stress
Units.
FATSEQ: Defines the loading sequence for the fatigue analysis. This card can refer to another
FATSEQ card or a FATEVNT card.
FATEVNT: Defines loading events for the fatigue analysis.
FATLOAD: Defines fatigue loading parameters.
The following files found in the optistruct.zip file are needed to perform this tutorial. Refer to
Accessing the Model Files.
ctrlarm.fem, load1.csv and load2.csv
In this tutorial, a control arm loaded by brake force and vertical force is used, as shown in Figure
6. Two load time histories acquired for 2545 seconds with 1 HZ, shown in Figure 7(a) and 7(b),
are adopted. The SN curve of the material used in the control arm is shown in Figure 8. Because
a crack always initiates from the surface, a skin meshed with shell elements is designed to cover
the solid elements, which can improve the accuracy of calculation as well.
(a)
Figure 8: SN curve
5. Click on the Table icon next to the Data field and select FATLOAD1 for FLOAD(1) and
FATLOAD2 for FLOAD(2) in the pop-out window.
5. Click next to the Data field and select shell for PID(1), and pfat for PFATID(1) in the
pop-out window.
6. Click Close.
Step 11: Submitting the OptiStruct analysis and reviewing the results
1. From the Analysis page, go to the OptiStruct panel.
2. Click save as following the input file: field. A Save As browser window opens.
3. In File name field, enter the name ctrlarm_hm.fem.
4. Click Save. The field input file: is set to the location of ctrlarm_hm.fem.
5. Submit the analysis by clicking on OptiStruct.
6. When the analysis process completes, click HyperView to launch the results.
7. In the Results tab, select Subcase 3 (Fatigue) from the subcase field.
8. Go to the Contour panel in HyperView.
9. Set Result type to Damage and click Apply to contour the elements.
Figure 1: Low Cycle and High Cycle regions on the S-N curve
Since S-N theory deals with uniaxial stress, the stress components need to be resolved into one
combined value for each calculation point, at each time step, and then used as equivalent
nominal stress applied on the S-N curve (Figure 2).
In OptiStruct various stress combination types are available, with the default being “Absolute
maximum principle stress”. In general “Absolute maximum principle stress” is recommended for
brittle materials, while “Signed von Mises stress” is recommended for ductile material. The sign
on the signed parameters is taken from the sign of the Maximum Absolute Principal value.
In this tutorial, you will be able to evaluate fatigue life with the S-N method through process
manager step by step. In brief, the following contents are covered:
Launching Fatigue Process Manager
Import a model
Creating fatigue subcase
Defining fatigue analysis parameters
Defining fatigue elements and S-N properties
Defining load-time history and loading sequence
Submitting the job
Viewing results summary and launch HyperView for post-processing
The following files found in the optistruct.zip file are needed to perform this tutorial. Refer to
Accessing the Model Files.
A control arm loaded by brake force and vertical force is used, as shown in Figure 3. Two load
time histories acquired for 2545 seconds with 1 HZ, shown in Figures 4(a) and 4(b), are applied.
The material of the control arm is aluminum, whose S-N curve is shown in Figure 5. Because a
crack always initiates from the surface, a skin meshed with shell elements is designed to cover
the solid elements, which can improve the accuracy of calculation as well.
(a)
Figure 4: Load time history (a) for vertical force (b) for braking force
Gate 0.0
4. Click Apply.
8. Click the Show SN curve definition icon . An SN method description window introducing
how to generate the SN material parameter opens.
9. Click Close.
10. For Material type, select Aluminum Alloys and click Estimate. All the data for SN curve
definition are automatically estimated.
11. Click Plot SN Curve at the bottom of the window to show the SN curve.
ctrlarm_fpmtut.h3d Hyper 3D binary results file, with both static analysis results and
fatigue analysis results.
ctrlarm_fpmtut.out OptiStruct output file containing specific information on the file set
up, the set up of your fatigue problem, compute time information,
etc. Review this file for warnings and errors.
Figure 1: Low Cycle and High Cycle regions on the S-N curve
Since E-N theory deals with uniaxial strain, the strain components need to be resolved into one
combined value for each calculation point, at each time step, and then used as equivalent
nominal strain applied on the E-N curve (Figure 2).
In OptiStruct various strain combination types are available with the default being “Absolute
maximum principle strain”. In general “Absolute maximum principle stain” is recommended for
brittle materials, while “Signed von Mises strain” is recommended for ductile material. The sign
on the signed parameters is taken from the sign of the Maximum Absolute Principal value.
A flowchart of the fatigue setup in HyperMesh can be described as shown in the image below.
The three aspects to the fatigue definition are the fatigue material properties, the fatigue
parameters and the loading sequence and event definitions.
FATDEF: Defines the elements and associated fatigue properties that will be used for the
fatigue analysis.
PFAT: Defines the finish, treatment, layer and the fatigue strength reduction factors for the
elements.
MATFAT: Defines the material properties for the fatigue analysis. These properties should be
obtained from the material’s E-N curve (Figure 2). The E-N curve, typically, is
obtained from completely reversed bending on mirror polished specimen.
FATPARM: Defines the parameters for the fatigue analysis. These include stress combination
method, mean stress correction method (Figure 4), Rainflow parameters, and Stress
Units.
The fatigue sequence and event definition
The following files found in the optistruct.zip file are needed to perform this tutorial. Refer to
Accessing the Model Files.
ctrlarm.fem, load1.csv and load2.csv
Exercise
In this tutorial, a control arm loaded by brake force and vertical force is used, as shown in Figure
6. Two load time histories acquired for 2545 seconds with 1 HZ, shown in Figure 7(a) and 7(b),
are adopted. The material of the control arm is aluminum, whose E-N curve is shown in Figure 8.
Because a crack always initiates from the surface, a skin meshed with shell elements is designed
to cover the solid elements, which can improve the accuracy of calculation as well.
(b)
Figure 7: Load time history (a) for vertical force (b) for braking force
5. Click on the Table icon next to the Data and select FATLOAD1 for FLOAD(1) and
FATLOAD2 for FLOAD(2) in the pop-out window.
Step 11: Submitting the OptiStruct analysis and reviewing the results
1. From the Analysis page, go to the OptiStruct panel.
2. Click save as following the input file: field. A Save file browser window opens.
3. Enter the name ctrlarm_hm.fem in File name field.
4. Click Save.
5. Submit the analysis by clicking on OptiStruct.
6. When the analysis process completes, click HyperView to launch the results.
7. Go to the Results tab.
8. Change the Load Case to Subcase 3 – fatigue.
9. Go to the Contour panel in HyperView.
10. Set Result type to Damage and click on Apply to contour the elements.
Figure 1: Low Cycle and High Cycle regions on the S-N curve
Since E-N theory deals with uniaxial strain, the strain components need to be resolved into one
combined value for each calculation point, at each time step, and then used as equivalent
nominal strain applied on the E-N curve (Figure 2).
In OptiStruct, various strain combination types are available with the default being “Absolute
maximum principle strain”. In general “Absolute maximum principle stain” is recommended for
brittle materials, while “Signed von Mises strain” is recommended for ductile material. The sign
on the signed parameters is taken from the sign of the Maximum Absolute Principal value.
In this tutorial, you will be able to evaluate fatigue life with the E-N method, which is a new
feature in HW10.0, through Fatigue Process Manager (FPM) step by step. In brief, the following
contents are covered:
Launching Fatigue Process Manager
Importing model
Creating fatigue subcase
Defining fatigue analysis parameters
Defining fatigue elements and E-N properties
Defining load-time history and loading sequence
Submitting the job
Viewing results summary and launch HyperView for post-processing
The following files found in the optistruct.zip file are needed to perform this tutorial. Refer to
Accessing the Model Files.
A control arm loaded by brake force and vertical force is used, as shown in Figure 3. Two load
time histories acquired for 2545 seconds with 1 HZ, shown in Figures 4(a) and 4(b), are applied.
The material of the control arm is aluminum, whose E-N curve is shown in Figure 5. Because a
crack always initiates from the surface, a skin meshed with shell elements is designed to cover
the solid elements, which can improve the accuracy of calculation as well.
(a)
1. Make sure the task Fatigue Subcase is selected in the Fatigue Analysis tree.
2. In the Create new fatigue subcase field, enter fatsub_fpmtut.
3. Click Create.
4. For the Select existing fatigue subcase:, select the newly created fatigue subcase
fatsub_fpmtut.
fatsub_fpmtut is selected as the active fatigue subcase. Definitions in the following
processes (analysis parameters, fatigue elements and properties, loading sequences, etc.) will
be for this subcase.
5. Click Apply. This saves the current definitions and guides you to the next task Analysis
Parameters of the Fatigue Analysis tree.
Gate 0.0
8. Click the Show EN curve definition icon . An EN method description window introducing
how to generate the EN material parameter opens.
9. Click Close.
Figure 1 below illustrates the structural model used for this tutorial: A 1mm thick cylindrical
gasket is sandwiched between two co-axial steel cylindrical tubes. The outer cylinder is subjected
to a pressure of 300MPa on the outer surface as shown. Using symmetry boundary conditions,
only a quarter of the geometry has been modeled. The gasket is connected to the inner and outer
cylinders using contact.
Exercise
6. Click the Table icon next to the Data field and enter the following values (X (closure) and
Y (pressure) fields) in the pop-out window.
7. Click Close.
For details on pressure-closure definitions of gaskets, refer to the HyperWorks 2017 online
help.
X Y
0.0 0.0
0.005 200.0
0.05 450.0
0.135 700.0
0.22 820.0
0.287 830.0
X Y
0.08 0.0
0.12 140.0
9. Next, create the second unloading curve named unload-curve2 with the following X-Y data:
X Y
0.17 0.0
0.2 250.0
0.22 820.0
10. Finally, create the third unloading curve named unload-curve3 with the following X-Y data:
X Y
0.23 0.0
0.265 360.0
0.287 830.0
15. Click next to the Data field and select the following:
For TABLU(1), select unload-curve1
For TABLU(2), select unload-curve2
9. Next, assign this property to the gasket component. Click on the component GASKET in the
Model browser.
10. For Property, select gasket_prop property.
gasket_complete.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.
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
can be obtained in the table below.
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 un-loading and is a continuation of the nonlinear solution sequence
from the previous loading subcase.
Top block 72 mm x 72 mm
Thickness of 20. mm
blocks
Note: If HyperMesh Desktop was launched, use: File > Open > Model.
5. Select the nlstat.hm file you saved to your working directory from the optistruct.zip file.
Refer to Accessing the Model Files.
6. Click Open.
6. Click the Table icon next to the Data field and enter the following values (x and y fields) in
the pop-up window.
nlstat_complete.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.
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.
Step 1: Launch HyperMesh, set the OptiStruct User Profile and retrieve
the structural model
1. Download both brake.fem and DMIG.pch files and save to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
2. Launch HyperMesh.
3. Select OptiStruct in the User Profiles dialog 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.
4. Click the Import Solver Deck panel toolbar icon . An Import tab appears. The File type:
is OptiStruct.
5. Click the open file icon in the File: field. A Select OptiStruct file browser window appears.
6. Select brake.fem from your working directory and click Open. The model file brake.fem is
loaded in the File field.
7. Click Import to import the file into the current HyperMesh session.
Step 3: Retrieve the friction data from DMIG data and define analysis
parameters
1. Go to the Analysis page, then click control cards.
2. Click INCLUDE_BULK.
3. Input the name of the include file, DMIG.pch.
4. Click return to go back to control cards.
5. Click K2PP.
6. In K2PP panel, set number_of_k2pps = 1.
7. In the field of K2PP=, enter KF.
KF is the name of the DMIG data entry.
8. Click return and back to control cards.
Retrieve the friction coefficients from the DMIG.pch file.
9. Click PARAM.
10. Check the small box in front of parameter G.
As you can see, the 5th mode was divergent while all of the other modes were stable. The
friction coefficient parameter can be reduced by setting the PARAM,FRIC factor from a value
of 0.05 to 0.01, and all roots become stable. It illustrates that there is a stability threshold
between the friction factor 0.05 and 0.01. It can be determined by resetting the scale factor
of PARAM, FRIC and rerunning the model till the damping value of this mode approaches zero.
2. Load the brake_complex.h3d file into HyperView to review complex eigenvectors.
Step 1: Launch HyperMesh, set the OptiStruct User Profile and retrieve
the structural model
1. Download brsq.fem file from the optistruct.zip file and save to your working directory.
Refer to Accessing the Model Files.
2. Launch HyperMesh.
3. Select OptiStruct in the User Profile dialog 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.
4. Click the Import Solver Deck panel toolbar icon . An Import tab appears. The File type: is
OptiStruct.
5. Click the open file icon in the File: field. A Select OptiStruct file browser window appears.
6. Select brsq.fem from your working directory and click Open. The model file brsq.fem is
loaded in the File field.
7. Click Import to import the file into the current HyperMesh session.
Step 1: Launch HyperMesh, set the OptiStruct User Profile and retrieve
the structural model
1. Launch HyperMesh.
2. Select OptiStruct in the User Profile dialog 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.
3. Click the Import Solver Deck panel toolbar icon . An Import tab appears. The File type:
is OptiStruct.
4. Click the open file icon in the File: field. A Select OptiStruct file browser window appears.
5. Select the rotor.fem file you saved to your working directory from the optistruct.zip file.
6. Click Open. The model file rotor.fem is loaded in the File field.
7. Click Import to import the file into the current HyperMesh session.
Step 3: Define Grids for the Rotor Line Model using the ROTORG set
entry and Reference Rotor speeds using the RSPINR entry
1. Right-click in the Model browser and select Create > SET.
2. Click Name and enter ROTORG_SET.
2. Load the rotor_sync.out file in a text editor. The Frequencies which you get from the
Synchronous Rotor dynamic analysis give you the critical frequencies. The complex modes
contain the imaginary part, which represents the cyclic frequency, and the real part which
represents the damping of the mode. If the real part is negative, then the mode is said to be
stable. If the real part is positive, then the mode is unstable. The eigenvalues of the complex
modes are shown below:
Step 1: Launch HyperMesh Desktop, set the OptiStruct User Profile and
retrieve the structural model
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. 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.
Note: If HyperMesh Desktop was launched, use File > Open >
Model.
8. Check all dofs (that is dof1 to dof6) with the value 0.000, confirm load types is set to SPC,
and click create. The constraints are created as shown in the figure below.
Step 4: Define the input response spectrum using the TABLED1 card
1. Go to the Utility tab. If the Utility menu is not displayed, select View > Browsers >
HyperMesh > Utility.
2. At the bottom of the Utility menu, click the FEA panel.
3. Under Tools, click TABLE Create.
4. Select Import Table under Options: and TABLED1 under Tables.
5. Click Next.
6. Under Options, select Create New Table.
7. For Name, enter tabled1_card.
8. Click Browse.
9. For Files of type: change to CSV (*.csv), select the file sourceFileTABELD1.csv (which
contains the ‘x’ and ‘y’ values to define the input response spectrum, with frequency plotted
on the x-axis and acceleration on the y-axis) located in your working directory from the
optistruct.zip file.
10. Click Open. If the Import TABLED1 GUI is minimized, click on it on the taskbar.
6. Click the Table icon next to the Data field. In the pop-out window, select tabled1_card
for TID(1) and enter 0.02 for DAMP(1). The damping value is in the units of fraction of
critical damping.
8. Click next to Data. In the pop-out window, select dti_card for the DTISPEC field, and for
SCALE, enter the value 9800.0.
9. Since the direction of excitation for the structure is the Global Z direction, enter 0.0 for X(0),
0.0 for X(1), and 1.0 for X(2) respectively.
10. Click Close to exit the window.
6. Click next to the Data field and enter the values 0.0, 0.02, 50.0 and 0.02 for f(1),
g(1), f(2) and g(2), respectively in the window.
7. Click Close to exit the window.
Note: Save the file in a folder different from the folders under
HyperWorks installation folder.
4. Click Save.
Note: The name and location of the file displays in the input file:
field.
Step 1: Load the OptiStruct User Profile and import the model
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. 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.
5. Click the Open file icon in the File: field. A Select OptiStruct file browser window opens.
6. Select the cover.fem file you saved to your working directory from the optistruct.zip file.
Refer to Accessing the Model Files.
7. Click Open. The location of cover.fem displays in the File: field.
8. Click Import. The cover.fem database is loaded into the current HyperMesh session.
9. Click Close.
The grid contribution of each grid in addition to the ERP results for the panel is output to an
.h3d file. ERP panel results are output to punch format.
5. Click return twice to return to the Analysis page.
Step 1: Load the OptiStruct User Profile and import the model
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. 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.
4. Click the Import Solver Deck icon . An Import window opens.
The Import type: is FE Model and the File type: is OptiStruct.
5. Click the Open file icon in the File: field. A Select OptiStruct file browser window opens.
6. Select the Rings.fem file you saved to your working directory from the optistruct.zip file.
Refer to Accessing the Model Files.
7. Click Open. The location of Rings.fem displays in the File: field.
8. Click Import. The Rings.fem database is loaded into the current HyperMesh session.
9. Click Close.
As described at the beginning of the tutorial, the heat transfer boundary condition SPC_heat and
heat flux input flux are predefined in the model. An OptiStruct steady state heat conduction
loadstep, referring to SPC_heat and flux, has been defined as well. The heat transfer results are
requested in loadsteps panel. Refer to tutorial OS-1080 for the details on how to define heat
transfer boundary condition, heat flux, and the output request.
Note that without PGAPHT, the heat cannot be transferred through the gap elements. In this
case, the outer ring remains zero degree and the inner ring takes all heat, as shown in the
temperature results in Figure 4. Run the completed model and compare the results with PGAPHT
to the results below.
Figure 1: Illustration of the model showing the cylinder head, engine block and head bolts.
Checkpoint
The model consists of eight predefined components along with their corresponding property and
material allocations. A contact surface (PT_Surf) has been defined, which is used for 3D
pretensioning of an existing pretension surface. The pretension sections for 1D pretensioning
have also been created on two of the four bolts and the sectioned bolts are reconnected using 1D
beam elements (via rigids). A predefined visualization aid is also available under View, which
allows the user to easily look at the pretensioned sections of the four bolts. Contact surfaces and
Contact Interfaces (TYPE=FREEZE) between the various parts have also been created so you can
focus on the Pretensioning aspect of the tutorial.
Tutorial Prerequisites
1. The latest version of Altair HyperMesh, Altair HyperView and Altair OptiStruct software
installations. Pretensioned Bolt analysis is available only in Altair HyperMesh version-
12.0.110, Altair HyperView version-12.0.110 and Altair OptiStruct version-12.0.202 and later.
2. The pretension.hm solver deck is located in the optistruct.zip file. Refer to Accessing the
Model Files.
Brief Introduction to Pretensioned Bolt Analysis
Step 1: Selecting the OptiStruct User Profile and loading the model into
HyperMesh
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. Click OK. This loads the OptiStruct 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.
User Profiles can also be accessed from the Preferences menu on the toolbar.
4. Click the Open Model panel toolbar icon and browse to the location that contains the
model.
5. Select the Pretension.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
6. Click Open. The model is now loaded into HyperMesh.
Checkpoint
This tutorial helps the you apply 1D and 3D bolt pretensioning to the four head bolts (two of
each) and then apply a pressure load to the constrained system. The applied pressure load
models the pressure on the inside walls of an IC engine due to combustion. Pressure within the
engine compartment varies with time (transient); however, you capture the response of the
system at a specific instant frozen in time. A constant single-valued pressure load of 1 Pascal is
applied to the inner walls of the cylinder head and the engine block.
Gasket behavior is nonlinear and it may undergo cycles of loading and unloading which lead to
changes in its properties at each step. In this tutorial, which focuses on 1D and 3D pretensioning,
the loading and unloading paths for the gasket material are pre-populated in the MGASK data
entry via the TABLES# entries referenced by corresponding load collectors. As a quasi-static
analysis is running, the initial applied pressure load is compared with corresponding values within
the loading/unloading path tables and the initial material properties of the gasket are
determined. The nonlinear properties of the gasket via the MGASK data entry are a function of
pressure and the closure distance (Refer to MGASK Bulk Data Entry for more information).
FREEZE contact has been predefined for all parts in contact.
8. Click on Gask_Load, in the Load Collector folder and then click the table icon next to the
Data: field.
9. Make sure that the values on the TABLES1 bulk data entry defining the gasket loading paths
are input, as shown below.
11. The gasket normal direction is now reviewed by clicking on normals in the Tool panel.
12. To select the gasket component, use the Show/Hide tool (Figure 12) to hide the cylinder
head thereby exposing the gasket to view.
13. Click on the Show/Hide icon, shown in Figure 12 and right-click on the cylinder head to hide
it from view. The gasket should now be visible.
14. In a similar fashion, hide (right-click) the engine block from view to be able to better visualize
the gasket normals.
15. Click the Show/Hide icon again to deselect it and select the gasket directly from the
graphics area and click display normals.
16. The gasket normals can be seen in the graphics window, as shown in Figure 14. Notice that all
the normals point in the negative Z direction.
Checkpoint
This concludes the review section of the tutorial. You will now focus on generating contact
interfaces, contact surfaces and applying pretensioning to the head bolts.
Tip: Care must be taken not to use CTRL+left mouse click while
zooming in and positioning the elements in the graphics area
for selection. Using CTRL+left mouse click can lead to the
model being rotated about an axis and thus disengaging from
the Y-Z plane of View1. It is recommended to only use
CTRL+right mouse click (dragging action) while working in
View1.
5. Select both fields under the Load Type column in the Pretension Manager window (Click on the
first field and then while holding down the CTRL key, click on the second field). Click on the
downward facing arrow next to the second field and select Force from the drop-down menu.
6. In a similar fashion, enter 4500.0 for both bolts in the Load Magnitude column.
7. Click Apply. A pretensioning force of 4500.0 N is applied to both 1D bolts, as shown in Figure
19.
8. Click on Add 3D Bolts and select Select Existing Surface from the drop-down menu.
9. Click on the Wireframe elements skin only icon to view the predefined contact surface
PT_Surf on the third bolt.
11. Select Force under the Load Type column and enter 4500.0 N for the Load Magnitude
column and click Apply.
A pretensioning force of 4500.0 N is applied normal to the PT_Surf surface, as shown in
Figure 21.
12. Click on Add 3D Bolts and select Create New Surface from the drop-down menu.
13. Utilize the click and drag technique (while holding down the shift key) described previously to
select the top of the fourth bolt, as shown in Figure 22.
17. Click OK in the Pretension Manager to view all four bolts with their respective pretensioning
forces, as shown in Figure 25.
5. Select Nonlinear quasi-static from the drop-down menu next to Loadstep type: in the
Loadstep Type tab.
6. Switch to the Load References tab and click on NLPARM in the list of subcase entries.
7. Click on Nlparm in the Available nonlinear parameters: section and then click on the right
facing arrow to add it to the selected nonlinear parameter: section.
8. Similarly, click on SPC in the subcase entry list and add the Available SPC constraint to the
Selected SPC constraints: section.
9. Follow the instructions in Steps 6 or 7 to add PRETENS_1 to the list from the PRETENSION
subcase entry section.
10. Click OK after all three subcase entries are added to the Pretension loadstep.
11. Right-click on Loadsteps in the Load Step browser and select New loadstep.
12. In the Loadstep name: field, enter Pressure and click Create.
13. Select Nonlinear Quasi-static from the drop-down menu next to Loadstep type: in the
Loadstep Type tab.
14. Switch to the Load References tab and click on NLPARM in the list of subcase entries.
15. Click on Nlparm in the Available nonlinear parameters: section and then click on the right
facing arrow to add it to the selected nonlinear parameter: section.
16. Similarly, click on SPC in the subcase entry list and add the Available SPC constraint to the
Selected SPC constraints: section.
Displacements, Element stresses, Contact pressures, Gasket Sealing Status, Gasket thickness
direction closure and Gasket thickness direction pressures are calculated and can be plotted using
the Contour panel in HyperView.
Figure 29: Displacement contour for the Pressure subcase after pretensioning.
Checkpoint
In Figure 29, the displacement plot after running the pressure subcases can be seen. The
maximum displacement is around 0.089 mm and it occurs in the region near the pretensioned
bolt heads.
6. Select Gasket Thickness-direction Pressure in the Contour panel and click Apply.
A contour plot of gasket pressure in the thickness direction is created, as shown in Figure 30.
The other components are hidden to be able to better view the pressure variation on the
gasket.
Checkpoint
The maximum pressure on the Gasket in the thickness direction is equal to 0.21 MPa.
Summary
The model consists of two cubes in contact and enforced displacement on the top compressing
the structure.
The following exercises are included:
Launch HyperMesh and importing the model
Define the surface for each part, which will be in contact with the other
Define the interface, which links the two surfaces together
Request contact related output
Perform the N2S and S2S analyses
Review and compare the results
Check the contacts that have been analyzed
Step 1: Selecting the OptiStruct User Profile and loading the model into
HyperMesh
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. In the next panel, select the settings to request contact related output, as shown below.
9. Click Next to locate the GAPPRM control card and click HMGAPST.
10. Set the VALUE to YES.
Outputs the open/closed status of the CGAPG elements.
Displacements, Element Stresses, Contact Force, Contact Deformation, Contact Status and
Contact Traction are calculated and can be plotted using the Contour panel in HyperView. Only
compare the Contact Traction between the N2S and the S2S run.
4. Select the page window layout icon to split the page into two windows.
5. Load the other model in the new window by clicking and selecting contact_N2S.h3d.
6. Click the Contour toolbar icon in one of the two windows.
7. For Result type, select Contact Traction/Normal(s).
8. Click Apply.
10. Right-click in the window that shows the contour and select Apply Style > Current Page >
All selected to view the same results for both models.
A contour plot of normal contact traction shows for both runs. The traction for the S2S run is
much more uniform than for N2S by comparing the maximum and minimum values.
Figure 12: Applying the setup in one window to the rest of the page.
7. Repeat Steps 8.3 through 8.6 for the component top and the set ^MasterGrids_Contact_#1.
Step 10: Reviewing the contact Status in HyperView for N2S and S2S
1. The contact status for both N2S and S2S can also be reviewed in HyperView, if the model
contained the contact result output request CONTF. To view this, repeat Step 7 while choosing
Contact Status/Normal(s) as the contour plot.
Summary
The model consists of two circular parts where the inner one is heated and the outer one cooled
down, leading to contact between the two. The effect of using several important contact settings
such as contact stabilization, clearance, adjust and smoothing on both the results and the
convergence behavior is considered.
Launching HyperMesh and importing the model
Create and assign the contact property PCONT
Run the model
Stabilization and clearance
Adjusting and running smoothing
Load and compare the results in HyperView
Prerequisites
1. The latest version of HyperMesh, HyperView and OptiStruct software installations. Some of
the presented features are available only in HyperMesh, HyperView, and OptiStruct version
14.0 and later.
2. The solver deck wheels_contact.fem is located in the optistruct.zip file. Refer to
Accessing the Model Files.
Step 1: Selecting the OptiStruct User Profile and loading the model into
HyperMesh
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. Click OK.
This loads the OptiStruct 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.
User Profiles can also be accessed from the Preferences menu on the toolbar.
4. Click the Import Solver Deck panel toolbar icon and browse to the location that contains
the model.
5. Select the wheels_contact.fem file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
6. Click Import.
The model is now loaded into HyperMesh.
Figure 3: Changing the PCONT contact property through the Entity browser.
6. In the Model browser, select the interface cont_interf to assign the property to the interface.
7. In the Entity browser, select Property Id as the property and change the PID to
cont_prop.
Displacements, Element Stresses, Contact Force, Contact Deformation, Contact Status and
Contact Traction are calculated and can be plotted using the Contour panel in HyperView. Only
compare the Contact Traction between the N2S and the S2S run.
10. In the Model browser, unselect the outer part of the structure.
Only the results on the contact surface will be visible
11. Right-click in the window that shows the contour and select Apply Style > Current Page >
All selected to view the same results for both models.
12. A contour plot of normal contact traction shows for both runs.
The traction for the runs with clearance and adjust are more uniform than they are for the
model with stabilization only. The surface smoothing leads to a more uniform contour. In
addition, the peaks are much lower for these three models. The reason why is, the traction is
much higher for adjust than it is for the clearance, and the adjust run is that for adjust, the
gap is closed initially, which leaves less room for stress free thermal expansion as for the
other runs.
Figure 10: Normal contact traction contour for the four different runs.
Model Description
The figure below illustrates the structural model used for this tutorial: a long rectangular plate
which is supported at one end and distributed load is applied on the top surface. The dimensions
of the plate and material parameters can be obtained in the table below.
Length 1000 mm
Width 200 mm
Thickness 4.0 mm
11. Select the degrees of freedom dof1, dof2 and dof3. Deselect all others.
12. For load types, select SPC.
13. Click create to create the boundary constraints.
14. Click return to go to the main menu.
13. Click the Table icon next to Data and enter the values shown below in the pop-out
window.
plate.out ASCII based output file of the model check run before the
simulation begins and gives some basic information on the results
of the run.
Figure 1: Illustration of the models: Follower Loads: Concentrated Forces – Beam (left), Pressures – Rubber Disk Model
(Right).
Summary
Beam Model
The model is a curved steel beam constructed with CHEXA elements. A Force of 100 N is applied
to the top cross-section of the beam. The bottom of the beam is constrained by single point
constraints (SPC).
The following exercises are included:
Launch HyperMesh and importing the model
Run the model without Follower Load activation
Activate Follower Load (Concentrated force)
Learn how to adjust nonlinear adaptive criteria (NLADAPT)
Request intermediate nonlinear results using NLOUT.
Review and compare results with and without follower loads. Also visualize the correct
animation of the deforming model as a result of NLOUT activation.
Step 1: Selecting the OptiStruct User Profile and loading the model into
HyperMesh
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. Click OK.
This loads the OptiStruct 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.
The User Profiles dialog can also be accessed from the Preferences menu on the toolbar.
4. Click the Import Solver Deck panel toolbar icon and browse to the location that contains
the model.
5. Select the beam_fllwer.fem file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
6. Click Import.
The model is now loaded into HyperMesh.
7. Select the NLADAPT load collector in the Select Loadcol dialog and click OK.
7. Select the NLOUT load collector in the Select Loadcol dialog and click OK.
Displacements and Element Stresses are calculated by default and can be plotted using the
Contour panel in HyperView.
Step 8: Selecting the OptiStruct User Profile and loading the model into
HyperMesh
1. Similar to Step 1, select the disk_fllwer.fem file you saved to your working directory from
the optistruct.zip file.
2. Click Import.
The model is now loaded into HyperMesh.
Step 9: Creating Follower Load Bulk Data Entries for the three options
In the beam model in the first section of this exercise, follower loads for concentrated forces do
not depend on the area of the element face to whose grids they are applied. Therefore, the
results will remain the same for any activation option chosen on either the FLLWER bulk data
entry or the PARAM,FLLWER entry. Additionally, since you only had one subcase, select the
parameter PARAM,FLLWER to activate follower loading for this model.
In this current disk model, the loads are pressure loads which may depend on the area of the
element faces to whom they are applied. Additionally, you have multiple subcases to showcase
the effect of different FLLWER options. In the FLLWER bulk data entry in the OptiStruct help, you
will see:
5. Similarly, create load collectors FLLWER_2 and FLLWER_3 and select OPT=2 and OPT=3
respectively for each entry.
Step 10: Referencing the FLLWER bulk entries in the Subcase section
The created FLLWER Bulk Data entries should now be selected in the Subcase section.
1. Click on the fllwer_1 Load Step to open it in the Entity Editor.
2. Check the box next to FLLWER in the Subcase Options section of the Entity Editor.
3. Click on the field next to ID and then on the yellow Loadcol button.
4. Choose the FLLWER_1 load collector in the Select Loadcol dialog and click OK.
8. Similarly, select the NLOUT load collector for fllwer_2, fllwer_3, and NO_fllwer subcases.
Step 12: Running the Disk model with Follower Loads activation
The model now consists of loads which have been identified as follower loads. Additionally, you
have learned how to request results at intermediate increments.
1. Follow the instructions provided in Step 2 to save and launch OptiStruct.
2. For File name, enter disk_fllwer_ON.fem.
If the job is successful, the new results files will appear in the directory from which
disk_fllwer_ON.fem was selected. Look in the disk_fllwer_ON.out file for error messages
that could help debug the input deck, if any errors are present.
Step 13: Comparing the Displacement Results between models with and
without Follower Load activation
1. Launch Hyperview.
2. Select the page window layout icon to split the page into four windows.
3. In the first window, load the disk_fllwer_ON.h3d file by clicking and navigating to your
working directory.
4. Select Subcase 1 (fllwer_1) and Load Factor = 1.000000E+00 in the Results browser.
5. Similarly load the disk_fllwer_ON.h3d file again in the second window and select Subcase 2
(fllwer_2). Repeat this process until the four windows contain results from four different
subcases.
6. Activate/Select the first window by clicking in it and then click the Contour toolbar icon .
7. For Result type, select Displacement (v).
8. Click Apply.
9. Activate/Select the second window by clicking in it and load the displacement results similarly.
However, since you have requested results for intermediate iterations via NLOUT, you will see
results for all intermediate iterations. This is the reason you navigate the final load factor to
compare results.
10. The Displacement contour results are now displayed as in the figure below. You can clearly
see that the activation of follower forces has modified the displacement profile significantly.
Additionally, you can see that since subcase 3 (OPT=3) updates the load direction but not the
area, the force (Pressure*Area) distribution at the grid points is low for Subcase 3 and the
Displacement results reflect this when compared to OPT=1 and 2.
Tutorial prerequisites
1. The latest version of the Altair HyperMesh, Altair HyperView and Altair OptiStruct software
installations. Finite sliding of Rack and Pinion Gear model is available only in Altair HyperMesh
version-14.0.120, Altair HyperView version-14.0.120 and Altair OptiStruct version-14.0.230
and later.
2. The following file found in the optistruct.zip file is needed to perform this tutorial. Refer to
Accessing the Model Files.
finite_sliding.hm
Step1: Selecting the OptiStruct User Profile and loading the model into
HyperMesh
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog.
3. Click OK.
This loads the OptiStruct 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.
User Profiles can also be accessed from the Preferences menu on the toolbar.
4. Click the Open Model panel toolbar icon and browse to the location that contains the
model.
5. Select the finite_sliding.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
6. Click Open.
The model illustrated in Figure 1 is now loaded into HyperMesh.
Checkpoint
This tutorial helps you to define finite sliding contact between the circular gear and rack. The
gear is held fixed at the center in all dof while the rack is given displacement in x dof but
constrained in all other dof. All constraints and enforced displacements have already been
defined in model. Contact surfaces to define the slave and master surfaces are also pre-defined
in the model. Contact stabilization has been defined for the contact to help stabilize any rigid
body motion before contact gets established. A very tiny end-of-subcase stabilization also has
been specified to overcome any temporary instabilities that may sometimes occur at end-of-
analysis.
Tip
You can review, in a similar manner, the remaining pre-defined data entries, like properties and
load collectors. The procedure for load collector review is not as straight forward as shown above
in some cases; however, this has been thoroughly illustrated in various other tutorials for your
benefit.
Checkpoint
This concludes the review of the model section of the tutorial. You will now focus on reviewing
contact surfaces and generating finite sliding contact interface.
1. Click Analysis > contactsurfs panel to review the already created contact surfaces.
2. Click solid faces > name and select rack.
3. Click review.
11. For slave choose the rack contact surface and click update.
12. Similarly, for master choose the gear contact surface and click update.
13. To review the interface, click review, shown in Figure 8.
14. Click on card image and then click edit to edit the contact interface.
15. Set TYPE as SLIDE.
16. Set DISCRET as S2S.
17. Set TRACK as FINITE.
Surface-to-surface a finite sliding contact without friction have been defined, shown in Figure
9.
18. Click return to exit the interfaces panel. The finite sliding contact definition is now complete.
4. Click Save.
Note the name and location of the rack_pinion.fem file displays in the input file: field.
5. Set the memory options: toggle to memory default.
6. Set the run options: toggle to analysis.
7. Set the export options: toggle to all.
8. Click OptiStruct.
This launches the job. If the job is successful, you should see new results files in the directory
from which rack_pinion.fem was selected. The rack_pinion.out file is a good place to look for
error messages that could help debug the input deck if any errors are present.
Checkpoint
You can animate the results using the Set Transient Animation Mode in HyperView.
1. Select other Result type in the Contour panel and click Apply.
Model Description
The figure below illustrates the fluid structural model used for this tutorial: the dimensions of the
beam are shown in Figures 1 and 2.
The AcuSolve fluid model (slab_dcfsi.inp) and the OptiStruct structural beam model
(Slab.fem) are located in the fsi_models.zip file. Refer to Accessing the Model Files.
Exercise
5. Click the icon and input the values, as shown below. See NLPARM Bulk Data Entry for more
information.
Step 11: Prepare and Submit the AcuSolve job with AcuSolve command
prompt
1. Open the AcuSolve input file (slab_dcfsi.inp) in a text editor and change the socket_host
parameter in the EXTERNAL_CODE block to your machines hostname and save the file.
2. Open the AcuSolve Cmd Prompt application and enter the command: acuRun-pb
slab_dcfsi -np 8 (shown below).
If the job is successful, you will see new results files in the directory where HyperMesh was
invoked. The Slab.out file is where you will find error messages that will help you debug your
input deck, if any errors are present.
The default files that will be written to your directory are:
slab.html HTML report of the analysis, giving a summary of the problem formulation
and the analysis results.
slab.out ASCII based output file of the model check run before the simulation
begins and gives some basic information on the results of the run.
slab.stat Summary of analysis process, providing CPU information for each step
during the process.
Model Description
The figure below shows the Fluid structural model used for this tutorial:
Exercise
2. Open the AcuSolve Cmd Prompt application and enter the command: acuRun-pb
FSI_AS_MANIFOLD -np 8 (shown below).
If the job is successful, you will see new results files in the directory where HyperMesh was
invoked. The FSI_OS_MANIFOLD.out file is where you will find error messages that will help
you debug your input deck, if any errors are present.
FSI_OS_MANIFOLD.out ASCII based output file of the model check run before the
simulation begins and gives some basic information on the
results of the run.
In this tutorial, two revolute joints, one ball joint, and one universal joint are created to constrain
the degrees of freedom (shown in the following figure), such that the remaining degree of
freedom will be just 1.
DOF = 3*6 – (5+5+4+3) = 1
Removes Removes Removes total
Type of Joint
translational dof rotational dof number of dof
Revolute 3 2 5
Universal 3 1 4
Ball (Spherical) 3 0 3
5. Select node ID 12910 as a node for first orientation. The vector 12319 to 12910 defines the
axis of rotation of the revolute joint.
6. Click create.
Next, the revolute joint at the lower left corner of body1 will be created.
17. Click the joint type: selection switch and select ball.
18. Select node ID 11104 as first terminal which belongs to body1.
19. Select node ID 7578 as second terminal which belongs to body2.
20. Click create.
This section describes how to view the results in HyperView which will be launched from within
the OptiStruct panel of HyperMesh.
HyperView is a complete post-processing and visualization environment for finite element
analysis (FEA), multi-body system simulation, video and engineering data.
1. Click on the Open Model icon . An Open Model browser window opens.
2. Select the slider_crank.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
3. Click Open. The model has five components and a few free nodes that will be used to create
bodies and joints for the MBD model.
4. Click elems and select two RBE2 elements that are inside a hole on the connecting rod.
Alternatively you could also use ‘elems by id’ and input IDs 18795 and 18796 for the two RBE2
elements.
Revolute 3 2 5
Fixed 3 3 6
Translational 2 3 5
5. Click create.
Next, the revolute joint between support and crank will be created.
6. Click the joint type: selection switch > revolute.
7. Select node ID 25472 as a first terminal (see the following figure).
8. Select node ID 15124 as a second terminal.
9. Select coordinates as first orientation and enter x= 0.0, y= 0.0, z= 1.0. The z-axis will be
the axis of rotation of revolute joint.
A revolute joint between the crank and connecting rod will be created next.
11. Click the joint type: selection switch > revolute.
12. Select node ID 25229 as a first terminal (see the following figure).
13. Select node ID 25231 as a second terminal.
14. Select the coordinates as the first orientation and enter x= 0.0, y= 0.0, z= 1.0. The z-axis
will be the axis of rotation of revolute joint.
20. Click create. A translational joint between the sliding block and ground will be created next.
21. Click the joint type: selection switch > translational.
22. Select node ID 14519 as a first terminal (see the figure below).
23. Select node ID 25228 as a second terminal.
24. Select coordinates as first orientation and enter x= 1.0, y= 0.0, z= 0.0.
X is the direction of translation.
25. Click create > return.
The gravity force that applies to the model and MBSIM bulk data card, which is to specify the
parameter for multi body simulation, will be created in this step.
This section describes how to view the results in HyperView which will be launched from within
the OptiStruct panel of HyperMesh.
HyperView is a complete post-processing and visualization environment for finite element
analysis (FEA), multi-body system simulation, video and engineering data.
7. The Animation Controls are in the panel next to the playback controls.
In this tutorial, you learn how to create a JOINT, a PFBODY, an MBMNTC and a multi-body
dynamics subcase.
The beam model consists of 10 different flexible bodies (PFBODY) and each body is to be
connected through a fixed joint and rigid element (RBE2).
There are two RBE2’s defined at the boundary of each body (one for each body at this boundary).
The fixed joint will be created using coincident nodes which are independent nodes of each of the
RBE2s.
These exercises show how to:
Load a model into HyperMesh
Create a component for JOINT
Create a fixed joint
Create a PFBODY
Create load collectors for MBMNTC and MBSIM
Create an MBMNTC
Create an MBSIM
Create a multi-body dynamics subcase
Submitting the job
Post-process the results in HyperView
Load the model
Step 1: Retrieve the structural model and define the OptiStruct template
1. Launch HyperMesh Desktop.
2. Set the User Profiles to OptiStruct and click OK.
3. Click File > Open > Model.
4. Select the cantilever_beam_MBD.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
5. Click Open.
1. Click Geometry > Create > Nodes > XYZ to open the Create Nodes panel.
2. Move the cursor to the upper left corner of body_1 and click three times (see the following
figure). You will get the nodal coordinate (x=, y=, z=) of that node.
3. Click create.
4. Repeat the last two steps to create another coincident node at lower left corner of body_1.
9. Under second terminal:, click node and select the other coincident nodes at the same
location.
Fixed joint
7. Repeat this process to create a fixed joint for the boundary of each body.
8. Click return.
The moment applied at the end of the beam, the gravity force that applies to the model and
MBSIM bulk data card, which is to specify the parameter for multi body simulation, will be
created in this step.
15. Select the 2 nodes at right tip of a beam, as shown in the following figure .
MBMNTC
16. Click the switch beside magnitude= and change to curve, components.
17. In the z comp= field, enter 1.0.
7. The Animation Controls are in the panel next to the playback controls.
8. With the animation running, use the bottom slider bar to adjust the speed of the animation.
9. Click the Start/Pause Animation icon again to stop the animation.
10. For Result type:, select Element Stresses (2D & 3D) (t).
11. Click Apply.
Step 2: Create load collectors that will be used to conduct the flexible
body reduction
In this step, two collectors will be created; one for the ASET that defines the connecting degrees
of freedom of the flexible body and the other for the method and parameters for the component
mode synthesis.
1. In the Model browser, right-click and select Create > Load Collector.
2. For Name, enter ASET.
3. Click Color and select any color.
4. Set Card Image to None. This load collector is used to define connecting degrees of freedom
of the flexible body to the multi-body system.
5. In the Model browser, right-click and select Create > Load Collector.
6. For Name, enter CMS.
7. For Card Image, select CMSMETH from the list. This creates the load collector used to
define the component mode synthesis method and parameters. The card image can now be
edited.
8. The Craig-Bampton (CB) method is selected by default.
Either UB_FREQ (upper bound on frequency) or NMODES (number of modes) must be
defined.
9. Under NMODES, enter 10.
10. Create the fourth constraint by selecting dof1, dof2, and dof3 and choosing the top node in
the rigid on the boss to the right.
11. Click create > return.
Step 6: Define the analysis type, unit system, create a title, and create
the output control for the component mode synthesis
1. Click on DTI_UNITS in the list of Control Cards to define the units system for the flex body
output. The units should be defined consistent with the material properties of the material
defined for this model. This way, you will not need to take care of the units of Multi-body
Dynamics Analysis.
For this model, the material used is Steel. Since the values of
the material properties provided are consistent with
Megagram, Newton, Millimeter, Second, the sequence: MGG
N MM S is selected for this control card.
2. Click return.
3. From the Control Cards panel, click GLOBAL_OUTPUT_REQUEST > STRESS.
4. Check the box preceding STRESS.
5. Select the options shown below.
6. Click return.
7. Click next and then click OUTPUT and set the options shown in figure below.
susp_sla.h3d H3D file that contains the flexible body data for use in Altair
MotionSolve.
susp_sla.out ASCII-based output file of the model check run performed before
the simulation begins which gives some basic information on the
results of the run.
susp_sla.stat Detailed breakdown on the CPU time used for each significant
stage in the analysis.
Exercise
Step 1: Retrieve the structural model and define the OptiStruct template
1. Launch HyperMesh Desktop.
2. Select OptiStruct in the User Profiles dialog and click OK.
3. Click File > Open > Model. An Open file browser window opens.
4. Select the for_contact_tutorial.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
5. Click Open. The model has five components and a few free nodes that will be used to create
bodies and joints for the MBD model.
9. Click return.
Here, all the necessary joints required for this model will be defined. Five joints for the model are
required, as shown in the table below.
Type of Joint Body 1 Body 2
5. Click create.
Next the fixed joint between the slider and the circle is created.
6. Select node ID 19399 as the first terminal.
7. Select node ID 19400 as the second terminal.
19. Click the joint type: selector switch and select translational.
20. Select the node ID 2 as the first terminal.
21. Select the node ID 4 as the second terminal.
22. Select the vector option and for the first orientation, select x-axis.
23. Click create.
24. Click return to exit the panel.
Note: These element sets are defined from the Analysis page, entity sets
by choosing a set of elements. The set of elements on the face of
the pentagon body is named master and the sets elements on the
face of the circle body is named slave.
10. Click the card image radio button and click edit.
11. Select POISSON for CNFTYPE and enter the values, as shown in the image below:
The motion which drives the mechanism, the gravity force that applies to the model and MBSIM
bulk data card, which is to specify the parameter for multi body simulation, will be created in this
step.
Note: There are a few more files written to the directory with
the name for_contact_tutorial_mbd.
8. With the animation running, use the bottom slider bar to adjust the speed of the animation.
9. Click the Start/Pause Animation icon again to stop the animation.
In this tutorial, a Curved Pentagon Positive Return Cam system with the help of a CVCV
constraint is modeled.
Exercise
Step 1: Retrieve the structural model and define the OptiStruct template
1. Launch HyperMesh Desktop.
2. Select OptiStruct in the User Profiles dialog and click OK.
3. Click File > Open > Model. An Open Model browser window opens.
9. Click return.
Here, you will define all the necessary joints including the CVCV joint. Five joints for the model
are needed. The first joint is the fixed joint between the stand and ground body. The second joint
is a revolute joint between the stand and driver, the third joint is the translational joint that
connects the slider to the stand, the fourth joint is the revolute joint between the driver and the
pentagon, and the fifth joint is the fixed joint between the slider and the circle.
5. Click create.
Next, the fixed joint between the slider and the circle is created.
6. Select node ID 19399 as the first terminal.
7. Select node ID 19400 as the second terminal.
8. Click create.
To create the revolute joints between (stand-driver and driver-pentagon) in the model, refer
to the image below and follow the steps.
18. Click the joint type: selector switch, and select translational.
Note: These curves are defined from the Analysis page, entity sets
by choosing a set of nodes. The curve on the pentagon body
is named master and the curve on the circle body is named
slave.
The motion which drives the mechanism, the gravity force that applies to the model and MBSIM
bulk data card, which is to specify the parameter for multi body simulation, will be created in this
step.
Note: You can find a new load collector (auto1) added to the
model after completing this step. The motion is assigned to
this load collector and will be used to reference in the
OptiStruct subcase.
8. With the animation running, use the bottom slider bar to adjust the speed of the animation.
9. Click the Start/Pause Animation icon again to stop the animation.
In this tutorial, you learn how to create JOINTS (Fixed, PTDCV), a PFBODY, a PRBODY, an
MBDCRV, and a multi-body dynamics subcase.
Model in HyperMesh
Exercise
Step 1: Define the OptiStruct User Profile and retrieve the finite element
model
1. Launch HyperMesh Desktop and select OptiStruct in the User Profiles dialog.
2. Click File > Open > Model. An Open Model browser window opens.
3. Select the flex_cable.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
4. Click Open. The flex_cable.hm database is loaded into the current HyperMesh session,
replacing any existing data.
5. Click create.
Next, the fixed joint between the other end of the Cable and Ground is created.
5. Click return.
The motion which drives the mechanism, the gravity force that applies to the model and MBSIM
bulk data card, which is to specify the parameter for multi body simulation, will be created in this
step.
flex_cable.html HTML report of the analysis, giving a summary of the problem formulation
and the results from the final iteration.
flex_cable.out OptiStruct output file containing specific information on the file set up,
estimates for the amount of RAM and disk space required for the run, and
compute time information. Review this file for warnings and errors.
flex_cable.log Log file containing the information on the joints and markers, simulation
etc., which are specific to MBD analysis.
flex_cable.stat Summary of analysis process, providing CPU information for each step
during analysis process.
This step describes how to view the results in HyperView which will be launched from within the
OptiStruct panel of HyperMesh.
HyperView is a complete post-processing and visualization environment for finite element
analysis (FEA), multi-body system simulation, video and engineering data.
View the animation
1. While in the OptiStruct panel of the Analysis page, click HyperView.
Note: That the path and file name for flex_cable.h3d appears in
the fields to the right of Load model and Load results. This is
fine because the .h3d format contains both model and results
data.
8. With the animation running, use the bottom slider bar to adjust the speed of the animation.
9. Click the Start/Pause Animation icon again to stop the animation.
15. In the Entity Editor, for Property, click Unspecified > Property.
4. Click return.
The finite element model, consisting of shell elements, element properties, material properties,
and loads and boundary conditions has been defined. Now a topology optimization will be
performed with the goal of minimizing the amount of material to be used. Typically, removing the
material in an existing volume with the same loads and boundary conditions makes the model
less stiff and more prone to deformation. Therefore, you need to track the displacements (which
represent the stiffness of the structure) and constrain the optimization process such that the
least material necessary is used and overall stiffness is also achieved.
The forces in the structure are applied on the outer nodes of the opening of the clip, making
those two nodes critical locations in the mesh where the maximum displacement is likely to
occur. In this tutorial, you will apply a displacement constraint on the nodes so that they would
not displace more than 0.07 in the y-axis.
cclip_complete.stat Contains information about the CPU time used for the
complete run and also the break-up of the CPU time for
reading the input deck, assembly, analysis, convergence,
etc.
OptiStruct provides element density information for all iterations, and also gives displacement
and von Mises stress results (linear static analysis) for the starting and last iterations. This
section describes how to view those results in HyperView.
9. Move the slider below Current value: to change the density threshold.
You will see the iso value in the graphics window update interactively when you scroll to a
new value. Use this tool to get a better look at the material layout and the load paths from
OptiStruct.
16. Change the Result type to Element Stresses (2D & 3D) (t).
17. For Averaging method:, select Simple and click Apply.
18. Click with the right button on the first window; select Apply Style, then Current Page, and
select Contour.
Exercise
This optimization problem applies the technique of minimum member size control on the
elements of the model to achieve a discrete solution.
The following exercise is an OptiStruct topology optimization and covers the following processes:
Import the .fem file into HyperMesh
Setup the minimum member size control problem using HyperMesh
Step 1: Launch the HyperMesh Desktop, Load the User Profile, and
Retrieve a File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK.
3. Click the Import Solver Deck icon . The Import tab appears.
4. Set the File type: to OptiStruct.
5. Select the cclip_complete.fem file you saved to your working directory from the
optistruct.zip file and click Open. Refer to Accessing the Model Files.
The location of the file displays in the File: field.
6. Click Import to import the model file.
7. Compare this image to the one you achieved in the previous optimization without the
application of minimum member size control, OS-2000.
The iso value plot displayed is similar to the one previously. Notice the smaller members in
the original iso surface plot are replaced by a more discrete rib pattern. This design is easier
to manufacture.
Finite element mesh containing designable (blue) and non-designable (yellow) material.
A finite element model representing the designable and non-designable material (shown in
figure) is imported into HyperMesh. Appropriate properties, boundary conditions, loads, and
optimization parameters are defined and the OptiStruct software determines the optimal material
distribution. The results (the material layout) are viewed as contours of a normalized density
value ranging from 0.0 to 1.0 in the design space. Isosurfaces are also used to view the density
results. Areas that need reinforcement will tend towards a density of 1.0.
The optimization problem for this tutorial is stated as:
Step 1: Launch the HyperMesh Desktop, Set the User Profile and Retrieve
the File
1. Launch HyperMesh Desktop.
2. Choose the OptiStruct in the User Profiles dialog 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.
3. From the File menu on the toolbar, select Open > Model.
4. Select the carm.hm file you saved to your working directory from the optistruct.zip file.
Refer to Accessing the Model Files.
5. Click Open.
11. In the Select Material dialog, select Steel as the material and click OK to complete the
material selection.
12. In the Model browser, right-click and select Create > Property.
13. For Name, enter nondesign_prop.
14. For Card Image, select PSOLID as the property type.
15. For Material, click Unspecified > Material.
16. In the Select Material dialog, select Steel as the material and click OK to complete the
material selection.
17. Open the Component folder to update the component collectors with the newly created
properties.
8. Select the node at the other end of the bushing (see the following figure) by clicking on it in
the graphics window.
9. Constrain dof2 and dof3; make sure only dofs 2 and 3 are checked.
10. Click create.
A constraint is created. A constraint symbol (triangle) appears in the graphics window at the
selected node. The number 23 is written beside the constraint symbol, indicating that dof2
and dof3 are constrained.
11. Click nodes and select by id from the extended entity selection window.
12. Enter the value 3239 and press ENTER.
13. This selects node ID 3239 (see the next figure), which corresponds to the shock absorber
mounting location.
14. Constrain only dof3.
15. Click create. A constraint is created. A constraint symbol (triangle) appears in the graphics
window at the selected node. The number 3 is written beside the constraint symbol, indicating
that dof3 is constrained.
Three separate forces in load collectors: brake, corner, and pothole with the component "design" turned off using the
display panel.
5. Click props , select design_prop from the list of props, and click
select.
6. Set type: as PSOLID.
7. Click create. A topology design space definition, design_prop, has been created. All
elements organized in this design property collector are now included in the design space.
8. Click return.
carm_complete.oss OSSmooth file with a default density threshold of 0.3. You may
edit the parameters in the file to obtain the desired results.
carm_complete.sh Shape file for the final iteration. It contains the material density,
void size parameters and void orientation angle for each
element in the analysis. The .sh file may be used to restart a
run and, if necessary, run OSSmooth files for topology
optimization.
Element density results are output to the carm_complete_des.h3d file from OptiStruct for all
iterations. In addition, Displacement and Stress results are output for each subcase for the first
and last iterations by default into carm_complete_s#.h3d files, where # specifies the sub case
ID. This section describes how to view those results in HyperView.
Step 15: View an Iso Value Plot on Top of the Element Densities Contour
This plot provides the information about the element density. Iso Value retains all of the
elements at and above a certain density threshold. Pick the density threshold providing the
structure that suits your needs.
1. From the Results menu, click Plot > Iso.
2. From the panel area (lower portion of the GUI), select Element Densities as the Result
type.
3. Click Apply.
4. Set the Current Value: to 0.15.
5. Move the slider below Current value: to change the density threshold.
You will see the iso value in the graphics window update interactively when you scroll to a
new value. Use this tool to get a better look at the material layout and the load paths from
OptiStruct.
Finite element mesh containing designable (red) and non-designable (blue) material.
Step 1: Launch the HyperMesh Desktop, Load the User Profile, and
Retrieve a File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK.
3. From the File menu on the toolbar, select Import > Solver Deck. An Import tab is added to
your tab menu.
sshield_analysis_frames.html HTML file used to post-process the .h3d with HyperView Player
using a browser. It is linked with the _menu.html file.
sshield_analysis_menu.html HTML file to post-process the .h3d with HyperView Player using a
browser.
Eigenvector results are output from OptiStruct for a normal modes analysis by default. This
section describes how to view the results in HyperView.
Type: Uniform
Value: 10
8. Click Apply.
9. Click the Animation Controls toolbar icon.
10. Move the Max Frame Rate: slider between 60 and 1 to increase or decrease the animation
speed.
You can also change the default values for Angular Increment: to refine your animation.
11. Click the Play/Pause Animation icon to start the animation. An animation of the mode
shape should be seen for the first frequency.
12. Click the Play/Pause Animation icon again to stop the animation.
13. On the Page Control toolbar, click the Page Delete icon to delete the HyperView page as
seen in the figure below.
sshield_optimiz The OptiStruct history file containing iteration number, objective function
ation.his_data values, and percent of constraint violation for each iteration.
sshield_optimiz HyperMesh command file used to organize elements into components based
ation.HM.comp.c on their density result values. This file is only used with OptiStruct topology
mf optimization runs.
sshield_optimiz HyperMesh command file used to organize elements into entity sets based on
ation.HM.ent.cm their density result values. This file is only used with OptiStruct topology
f optimization runs.
sshield_optimiz HTML report of the optimization, giving a summary of the problem formulation
ation.html and the results from the final iteration.
sshield_optimiz OSSmooth file with a default density threshold of 0.3. You may edit the
ation.oss parameters in the file to obtain the desired results.
sshield_optimiz OptiStruct output file containing specific information on the file setup, the
ation.out setup of your optimization problem, estimates for the amount of RAM and disk
space required for the run, information for all optimization iterations, and
compute time information. Review this file for warnings and errors.
sshield_optimiz Shape file for the final iteration. It contains the material density, void size
ation.sh parameters and void orientation angle for each element in the analysis. The
.sh file may be used to restart a run and, if necessary, to run OSSmooth files
for topology optimization.
sshield_optimiz Summary of analysis process, providing CPU information for each step during
ation.stat analysis process.
sshield_optimiz HTML file used to post-process the .h3d with HyperView Player using a
ation_frames.ht
browser. It is linked with the _menu.html file.
ml
sshield_optimiz HTML file used to post-process the .h3d with HyperView Player using a
ation_menu.html
browser.
sshield_optimiz HyperView binary results file that contains: Density results from topology
ation_des.H3D optimizations, Shape results from topography or shape optimizations and
Thickness results from size and topology optimizations.
sshield_optimiz HyperView binary results file that contains: Displacement results from linear
ation_s1.H3D static analysis, Element strain energy results from normal mode analysis
and Stress results from linear static analysis, etc.
The name and location of the sshield_optimization.fem file displays in the input file: field.
With topology optimization of shell elements, Element Density and Element Thickness results are
output from OptiStruct for all iterations. In addition, Eigenvector results are output for the first
and last iterations by default. This section describes how to view those results in HyperView.
Each element of the model is assigned a legend color, indicating the thickness of each
element for the selected iteration.
Have most of your elements converged to a thickness close to 1 or 0?
Contour plot of element densities at iteration 6 with all components except the designable component turned off. (top
view)
7. On the Page Control toolbar, click the Page Delete icon to delete the HyperView page.
Based on the topology results obtained above, a number of ribs were added to the model. The
new design sshield_newdesign.fem, which includes these ribs can be found in the
optistruct.zip file.
5. Click Close.
Type: Uniform
Value: 10
7. Click Apply.
8. Click Pause/Play Animation to start the animation . An animation of the mode shape
should be seen for the first frequency.
9. Click Pause/Play Animation again to stop the animation.
What is the percentage increase in frequency for your first mode (sshield_analysis.fem vs.
sshield_newdesign)?
You have seen that the frequency of the structure for the first mode has increased from 43.63 Hz
to 84.88 Hz.
How much mass has been added to the part (check the mass of your ribs in the mass calc panel
in the Tool page)?
What is the percentage increase in mass?
Step 1: Launch the HyperMesh Desktop and Load the User Profile
1. Launch HyperMesh Desktop. The User Profiles dialog appears.
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.
controlarm_opt.sh Shape file for the final iteration containing the material
density, void size parameters, and void orientation angle
for each element in the analysis. The .sh file may be used
to restart a run and, if necessary, run OSSmooth files for
topology optimization.
OptiStruct provides the Element density information for all of the iterations. OptiStruct will also
show Displacement and von Mises stress results of a linear static analysis for iteration 0 and
iteration 38. This section describes how to view those results in HyperView. First, the HyperMesh
binary results file needs to be loaded.
8. You can also move the slider bar to change current value of Element densities.
The stress results are also available for the respective iterations.
7. Similarly, results for Load Case 2 can be viewed on page 4 using Next Page.
Final configuration.
The following file can be found in the optistruct.zip file and copied to your working directory.
Refer to Accessing the Model Files.
Step 1: Launch the HyperMesh Desktop, Set the User Profile, and Import
the File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK.
3. From the File menu on the toolbar, select Import > Solver Deck. An Import tab is added to
your tab menu.
4. Click the Select file icon to locate and open the no_repeat.fem file you saved to your
working directory from the optistruct.zip file. Refer to Accessing the Model Files.
5. Click Import. The no_repeat.fem OptiStruct input file is loaded into the current HyperMesh
session.
6. Click Close.
Step 1: Launch HyperMesh, Set the User Profile, and Import a File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog 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.
The User Profiles GUI can also be accessed from the Preferences menu on the toolbar.
3. Click the Import Solver Deck icon and select the carm_draw_symm.fem file you saved to
your working directory from the optistruct.zip file. Refer to Accessing the Model Files.
4. Click Import.
This forces the diameter or thickness of any structural member to be higher than 16 mm; if
this is not user-set, OptiStruct automatically selects a minimum member size based on the
average mesh size (if a manufacturing constraint is selected).
5. Click update to confirm the minimum member size set up.
6. Select the draw subpanel to set up the draw direction.
7. Set the selector under draw type: to single.
8. Select the anchor node and the first node as shown below. Together, these two nodes
define a vector in the positive Z direction. This defines that the die draw direction is along the
positive Z direction.
Step 5: View an Iso Value Plot on Top of the Element Densities Contour
This plot provides the information about the element density. Iso Value retains all of the
elements at and above a certain density threshold. Pick the density threshold providing the
structure that suits your needs.
1. From Graphics menu, click Iso, and select Element Densities as the Result type.
2. Set the Current Value: to 0.2.
3. Click Apply.
An iso value plot is displayed in the graphics window. The parts of the model with densities
greater than the specified value of 0.2 are shown in the figure below.
Review
Design variables: The density for each element in the design space.
1. Click the Open Model icon and load the cantilever_full.hm file you saved to your
working directory from the optistruct.zip file. Refer to Accessing the Model Files.
Complete model.
Generate a Superelement
Elements to be deleted.
The matrices are written to the .pch file with the same format as the DMIG bulk data entry.
They are defined by a single header entry and one or more column entries. By default, the
name of the stiffness matrix is KAAX, the mass is MAAX, and the load is PAX. Since mass
matrix is not used in this tutorial, it is not written to .pch file.
The I/O Option entry, DMIGNAME, provides you with control over the name of the matrices.
Step 10: Delete the Superelement Part Reduced out Using DMIG
Since the matrices for the superelement part will be replaced by DMIG, the bulk data entries for
the nodes and elements, as well as all loads and boundary conditions that are in the
superelement, should be deleted.
1. Press F2 to open the Delete panel.
2. Click the entity selection switch and select elems from the pop-up menu.
3. Click elems and select by window from the pop-up menu.
4. Move the mouse over to the graphics area and draw a window around the elements, as shown
in the following figure.
Element density and element thickness results are output from OptiStruct for all iterations. In
addition, displacement and stress results are output for the first and last iteration by default. This
section describes how to view those results in HyperView.
The resulting contours represent the element densities field resulting from the applied loads
and boundary conditions.
6. The current value can be changed to see more material. Enter a value of 0.3 in the Current
value: field.
You can also move the slider bar to change current value of Element densities.
The isosurface post-processing feature in HyperView is an excellent tool for viewing the
density results from OptiStruct.
Click and move the slider bar (currently pointing to a value representing 0.3) for your density
to change the isosurface.
The isosurface is available in the graphics window interactively update when you change it to
a new value. Use this tool to get a better look at the material layout and the load paths from
OptiStruct.
FEA model
The structural model is loaded into HyperMesh Desktop. The constraints, loads, subcases and
material properties are already defined in the model. The topology design variables and the
optimization problem setup will be defined using HyperMesh, and OptiStruct is used to determine
the optimal material layout. The results can then be reviewed in HyperView.
The optimization problem is stated as:
3. Click the Select file icon to locate and open the hook.fem file you saved to your working
directory from the optistruct.zip file. Refer to Accessing the Model Files.
4. Click Open.
5. Click Import > Close.
3. Click to fit the model to the screen. This displays only the components in the graphics
area.
5. Click create. A response, mass, is defined for the total volume of the model.
6. Click return to go back to the Optimization panel.
The isosurface post-processing feature is an excellent tool to use for viewing the density results from OptiStruct
(Current value=0.4528).
You will see the isosurface in the graphics window interactively update when you change it to
a new value. Use this tool to get a better look at the material layout and the load paths from
OptiStruct.
Notice that there are some local regions where the stresses are still high; this is because
topology stress constraints should be interpreted as global stress control or global stress
target.
The functionality has some ways to filter out the artificial or local stresses caused by point
loading or boundary conditions, but those artificial stresses will not be completely removed
unless the geometry is changed by shape optimization.
Note: There might still be high local stress regions which can be
improved more effectively with local shape and size
optimization.
The DTPL (Design Variable for Topology Optimization) card is used for this optimization.
The finite element mesh of the curved beam is shown in the following.
Finite element mesh of the curved beam showing loads and boundary conditions
Step 1: Launch HyperMesh, Set the User Profile, and Retrieve a File
1. Launch HyperMesh Desktop.
2. Click the Import Solver Deck icon . An Import tab is added to the tab menu.
3. Select the rail_complete.fem file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
4. Click Import > Close.
5. Click props , select new_solid from the list of props, and click
select.
6. Set type: as PSOLID.
7. Click Create. A topology design space definition, design_solid, has been created. All
elements organized in this design property collector are now included in the design space.
8. Click return twice to go back to the main menu.
A line of nodes starting from 71559 and ending with node 70001 should be highlighted,
indicating the extrusion path.
It is not required to select as many nodes to define the curve. This is an exercise to show that
the nodes by path option can also be used.
13. Click return to go back to the optimization panel.
Constraint set up
3. Go to the Iso icon and set the Result type: to Element Densities.
4. Click Apply.
5. Set Current Value: to 0.3.
6. Click Apply.
Isosurface plot of a curved beam rail layout of the topology optimization with extrusion constraints
As expected, the result with manufacturing extrusion constraints permits a constant cross
section for the entire length of the model.
As expected, the result with manufacturing extrusion constraints shows constant cross section
through the length of the model.
Plate model
Step 1: Launch the HyperMesh Desktop, Set the User Profile and Import
a File
1. Launch HyperMesh Desktop.
2. A User Profiles GUI will appear, select OptiStruct.
3. 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.
The User Profiles GUI can also be accessed from the Preferences menu on the toolbar.
4. Click the Import Solver Deck icon .
5. Click the File icon and select the frf_response_input.fem file you saved to your working
directory from the optistruct.zip file. Refer to Accessing the Model Files.
6. Click Import. The frf_response_input.fem OptiStruct input file is loaded into the current
HyperMesh session.
7. Click Close to exist the Import panel.
Step 11: Create a Set of Outputs and Include Damping for Frequency
Response Analysis
1. On the Analysis page, click the control cards panel and enter the
GLOBAL_OUTPUT_REQUEST panel.
2. Check DISPLACEMENTS. A new set of options appears in the work area screen.
3. Click the field box for FORM(1) and select PHASE from the pop-up menu.
4. Click the field box for OPTION(1) and select SID from the pop-up menu.
A new field appears in yellow.
5. Double-click the yellow SID(1) box and select SETA from the pop-up selection on the bottom
left corner.
A value of 1 now appears below the SID field box. This sets the output for only the nodes in
set 1.
5. Click the Define Curves button and delete the X Trans and Y Trans curves.
The excitation is applied on Z direction then, the main effect will be detected on this direction.
8. Click the Axes button and change to Vertical to edit the vertical axis.
9. Click the Scale and Tics (Magnitude) tab and change it to Logarithmic.
10. Change the Min: to 5 and the Max: to 200000.
11. Click the Scale and Tics (Phase) tab and change the Tics per axis: to 7.
12. Click Horizontal X or change it using the drop-down list and change the Min: to 5 and the
Max: to 1000.
15. Click the maximum button to see the maximum Y-magnitude ~ 15055 in the table. The peak
displacement of the baseline model.
Element Density and Element Thickness results are output from OptiStruct for all iterations. This
section describes how to view an iso-plot of the new design in HyperView.
3. Click the Build Plots icon to prepare to add curves on top of the existing analysis
information.
4. For the Data File:, load the optimization file with the final iteration analysis in it:
frf_response_optimization_s2.h3d.
5. For Subcase:, select the last iteration. Leave the X Type: as Frequency. Select Y Type:
Displacement (Grids), Y Request: N3, and Y Component: X,Y, and Z.
6. Click Apply to overlay the new information onto the original plot.
Excavator model
Design variables: Element density of elements in the lower arm (flexible body)
component.
Step 1: Launch HyperMesh, Set the User Profile, and Retrieve the File
1. Launch HyperMesh Desktop.
2. A User Profiles GUI will appear, 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.
5. Close the HyperView window using the Delete Page icon and return back to HyperMesh.
7. Click update.
A message appears in the lower left corner to indicate that the body has been update to a
new type.
8. Click return.
5. Click props and choose lowerarm from the list of props; click
select.
6. Choose type: PSOLID.
7. Click create.
8. Activate the parameters option.
9. Click the toggle next to minmemb off to activate mindim= and set the value to 0.05.
10. Click update.
Minimum member size control of 0.05 meter has been applied on this topology optimization
problem.
11. Click return to go back to the optimization panel.
Step 12: View element density results for the last iteration
1. When the message OPTIMIZATION HAS CONVERGED is received in the command window,
close the DOS window.
2. Click HyperView from within the OptiStruct panel. HyperView is launched and the results are
loaded.
A message window appears to inform about the successful loading of the model and result
files into HyperView.
12. You can also move the slider bar to change current value of Element densities to look at
material distribution for different density threshold values.
Finite element model of the design space with loads and constraints.
A finite element model (shown in the above figure) is loaded into HyperMesh. The constraints,
load, material properties, and subcase (loadstep) of the model are already defined. Topography
design variables and optimization parameters are defined and the OptiStruct software determines
the optimal reinforcement patterns. The results are viewed as animations of the contours of
shape changes of the design space. Finally, the use of the grouping patterns is shown; based on
the shape changes suggested by OptiStruct, a possible pattern is chosen for ease of
manufacturing.
The optimization problem for this tutorial is stated as:
Step 1: Launch the HyperMesh Desktop, Load the User Profile and
Retrieve the File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK. This loads the OptiStruct 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.
User Profiles can also be accessed from the Preferences menu on the toolbar.
3. From the File menu, click Open
4. Select the torsion_plate.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
5. Click Open. The torsion_plate.hm database is loaded into the current HyperMesh session,
replacing any existing data.
torsion_plate_s1_h3d HyperView binary results file for displacement and stress results
for subcase 1.
torsion_plate.sh Shape file for the final iteration. It contains the material density,
void size parameters, and void orientation angle for each
element in the analysis. The .sh file may be used to restart a run
and, if necessary, run OSSmooth files for topology optimization.
Shape contour information is output from OptiStruct for all iterations. In addition, Displacement
and Stress results are output for the first and last iteration by default. This section describes how
to view those results using HyperView.
3. Click the Contour icon on the toolbar to enter the contour panel.
5. Click to choose the last iteration from the Simulation list. A deformed plate appears.
6. Click Apply.
Is the max= field showing 4.0e + 00?
In this case, it is. If it is not, your optimization has not progressed far enough. Decrease the
OBJTOL parameter (set in the opti control panel on the optimization panel). This value,
4.0e+00, comes from the draw height defined earlier.
Contour plot showing the reinforcement pattern at the last iteration (converged solution)
1. Click the forward arrow icon to go to the next page. This page has the subcase information
from the torsion_plate_s1.h3d file.
2. Verify that the animate mode menu is set to Linear Static. For a better visual of what it
happening with this model, turn on mesh lines and contour the results.
3. Click the Deformed icon on the toolbar to enter the panel.
4. Set Result type: to Displacement(v).
5. Select the first iteration (Iteration 0) in the Results browser, as shown below.
Contour plot showing the reinforcement pattern with pattern grouping constraint at the last iteration
L-bracket layout.
Step 1: Launch the HyperMesh Desktop, Set the User Profile, and
Retrieve the File
1. Launch HyperMesh Desktop.
2. Select OptiStruct in the User Profiles dialog 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.
3. From the File menu on the toolbar, select Open.
4. Select the Lbkttopog.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
5. Click Open.
lbkttopog.hgdata HyperGraph file containing data for the objective function, constraint
functions, design variables, and response functions for each iteration.
lbkttopog.hist An OptiStruct output file for xy plotting containing the iteration history
of the objective function, maximum constraint violation, design
variables, DRESP1 type responses, and DRESP2 type responses.
lbkttopog.oss OSSmooth file with a default density threshold of 0.3. You may edit the
parameters in the file to obtain the desired results.
lbkttopog.out The OptiStruct output file containing specific information on the file set
up, the set up of the optimization problem, estimate for the amount of
RAM and disk space required for the run, information for each
optimization iteration, and compute time information. Review this file
for warnings and errors that are flagged from processing the
lbkttopog.fem file.
Lbkttopog_s1_h3d HyperView binary results file for displacement and stress results for
subcase 1.
lbkttopog.sh Shape file for the final iteration. It contains the material density, void
size parameters and void orientation angle for each element in the
analysis. The .sh file may be used to restart a run and, if necessary,
run OSSmooth files for topology optimization.
lbkttopog.stat Summary of analysis process, providing CPU information for each step
during analysis process.
3. Select the directory where you would like to write the OptiStruct model file and enter the
name for the model, Lbkttopog.fem, in the File name: field. The .fem extension is used for
OptiStruct input decks.
4. Click Save.
The name and location of the Lbkttopog.fem file displays in the input file: field.
5. Set the export options: toggle to all.
6. Click the run options: switch, and select optimization.
7. Set the memory options: toggle to memory default.
8. Click OptiStruct. This launches the OptiStruct job.
If the job was successful, new results files can be seen in the directory where the OptiStruct
model file was written. The lbkttopog.out file is a good place to look for error messages that
will help to debug the input deck if any errors are present.
Shape contour information is output from OptiStruct for all iterations. In addition, Eigenvector
results are output for the first and last iteration by default. This section describes how to view
those results in HyperView.
Topography results
L-bracket layout
Step 1: Launch HyperMesh, Set the User Profile, and Import a File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK.
User Profiles can also be accessed from the Preferences menu on the toolbar.
3. Click the Import Solver Deck . An Import tab is added to your tab menu.
4. Set the Import type: to FE Model.
5. Choose the File type: OptiStruct.
6. For File:, click the Select file icon to open the OptiStruct file browser.
7. Select the Lbkttopog_bead.fem file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
8. Click Import. The Lbkttopog_bead.fem OptiStruct input file is loaded into the current
HyperMesh session.
9. Click Close to exit the Import tab.
Shape contour information is output from OptiStruct for all iterations. In addition, eigenvector
results are output for the first and last iteration by default. This section describes how to view
those results in HyperView.
4. Ensure that the animation mode is set to Transient and click the Play icon to start the
animation.
An animation of the shape changes over the course of the optimization is displayed.
5. To slow down the animation, move the Animation Controls slider underneath the Current
Notice the topography optimization yields an almost 100% increase in the frequency of the
first mode by reviewing the Mode 1-F value in the Simulation list. To animate the model, click
the modal icon, next to the director's chair.
2. For file:, select the OptiStruct base input file from which to extract the final geometry.
3. For output:, select the IGES output format of the final geometry.
The default output format is STL. Other format options are: Mview, Nastran, IGES, and
H3D.
If you select IGES as the output format, select the output unit type. The default is mm
(millimeters).
4. Select load geom to load the new geometry into the current HyperMesh session.
5. Check the box next to autobead and enter a value of 0.3 for the bead threshold:.
6. Leave the rest of the options at their default settings.
7. Click OSSmooth.
8. Click Yes to overwrite.
The new geometry will be automatically loaded into the existing HyperMesh file, turn off the
display of all the elements to view the new concept geometry.
9. OSSmooth can automatically create geometry based on the new mesh. Click FE > Surf to
generate new geometry from the optimization results.
10. Click Save and Exit to continue.
11. Using the Mask tab, click Isolate for Geometry and on Hide for Load Collectors.
Model review
Step 1: Load the OptiStruct User Profile and Import the Model
1. Launch HyperMesh.
2. Choose OptiStruct in the User Profiles dialog and click OK. This loads the OptiStruct user
profile. It includes the appropriate template, macro menu, and import reader, paring down
the functionality of HyperMesh to what is relevant for generating OptiStruct models.
User Profiles can also be accessed from the Preferences menu on the toolbar.
3. Click the Import panel toolbar icon . The File type: is OptiStruct.
4. Click the open file icon . A Select OptiStruct file browser window opens.
5. Select the panel.fem file you saved to your working directory from the optistruct.zip file.
Refer to Accessing the Model Files.
6. Click Open. The location of panel.fem displays in the file: field.
7. Click Import. The panel.fem database is loaded into the current HyperMesh session.
HyperView is used to view the bead design generated from the topography optimization.
“XYPUNCH, ACCE, PSDF/67(T1RM)” was used to output the PSD accelerations to punch files. The
PSD plot from punch output can be viewed with HyperGraph. The RMS and peak PSD values are
output to the .peak file and can be viewed with text editor.
Problem Statement
Perform combined topology and topography optimization on a disk drive slider suspension to
maximize the stiffness and weighted mode. The lower bound constraint on the seventh mode is
12 cycles/ms.
Step 1: Set the User Profile and Import a Finite Element Model
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog 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.
3. From the File menu on the toolbar, select Import > Solver Deck. An Import tab is added to
your tab menu.
4. Choose the proper File type: OptiStruct.
5. Click the Select Files icon and open the combined.fem file you saved to your working
directory from the optistruct.zip file. Refer to Accessing the Model Files.
6. Click Open.
7. Click Import.
8. Click Close to close the Import tab menu.
1 1.0
2 2.0
3 1.0
4 1.0
5 1.0
6 1.0
Problem Definition
The finite element model of the underbelly fairing was generated in HyperMesh, as shown in the
following image. Material properties for carbon-fiber were considered and represented using an
orthotropic material (MAT8) for two dimensional elements. The fairing was modeled with four ply
orientations (0°, 90°, 45° and -45°) of uniform thickness. The SMEAR option is applied in the
PCOMP card to eliminate stack biasing.
Two load cases were defined to represent the operating conditions – an internal uniform pressure
loading of 0.02MPa and an external gravity loading of 6.75g. The fairing was considered to be
riveted along its edges to the surrounding structure. Two equipment masses, weighing 2Kg and
3Kg each, were mounted to the fairing through RBE3 elements. The fairing has been designed
considering two major performance criteria: the first natural frequency is at least 20Hz, and the
maximum strain is less than 1000 micro-strain.
Model overview
Step 1: Load the OptiStruct user profile and import the model into
HyperMesh
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog 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 OptiStruct.
User Profiles can also be accessed from the Preferences menu on the toolbar.
3. Click the Import Solver Deck icon .
4. The File type: is OptiStruct.
5. Click the open file icon in the File: field. A Select OptiStruct file browser window opens.
6. Select the fairing.fem file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
7. Click Open. The location of fairing.fem displays in the file: field.
8. Click Import. The fairing.fem database is loaded into the current HyperMesh session.
6. Verify that the fields in the create subpanel match the fields shown above and click create.
The design variable fairing is created for the free-size optimization. The manufacturing
constraints on ply percentage and ply balance will be defined next.
7. Go to composites subpanel. Make sure fairing is selected as the design variable.
8. Click edit and enter the DSIZE panel to define the manufacturing constraints on ply
percentage and ply balance.
9. Check the box in front of PLYPCT.
10. Set Ply Percentage Options to BYANG.
11. In the DSIZE_NUMBER_OF_PLYPCT = field, enter the value of 1. A PLYPCT continuation
line is added to the DSIZE data entry.
12. Check the box in front of BALANCE.
13. In the field of DSIZE_NUMBER_OF_BALANCE=, enter the value of 1. A BALANCE
continuation line is added to the DSIZE data entry.
14. Check the box in front of PLYMAN. A PLYMAN continuation line is added to the DSIZE data
entry.
15. Define the PLYPCT, BALANCE and PLYMAN constraints, as shown below.
7. Click create.
8. Enter volfrac in the response= field.
9. Set the response type: to volumefrac.
10. Make sure the total is selected and the toggle is set to no regionid.
5. Click create.
6. Click return.
4. Click create.
5. Click return twice to go back to the Analysis page.
Figure 7: Requesting CSTRAIN and CSTRESS results output to the .h3d file
4. Click return.
5. Click next and click the OUTPUT card.
6. In the OUTPUT panel, enter 1 as the number_of_outputs.
Figure 8: Requesting the free-size to size (FSTOSZ) optimization output file for Phase 2.
3. Select the last iteration in the Load Case and Simulation Selection drop-down in the
Results browser.
4. Click Apply and then XY Top Plane View to view the results in the X-Y plane.
The element thickness results from the free-size optimization are shown in the following
image. The regions indicated in red or in colors tending towards red (from the legend) can be
interpreted as thicker regions, while those in blue or tending towards blue are thinner regions.
The contour plot indicated above is the total thickness distribution that includes contributions
from each ply orientation, i.e. a thickness contribution from the 0s, +/-45s and the 90s. It
2. Select the last iteration in Load Case and Simulation Selection drop-down (similar to
Figure 12).
3. Click Apply.
The thickness distribution of 0 degree super ply is generated and shown in Figure 15. It
represents the ply shapes and patch locations of the 0 degree ply bundles.
4. Repeat step 1 through 3 to create the ply thickness contours for super-ply 2 (45°), 3 (-45°),
and 4 (90°) by selecting Layers 2, 3 and 4, respectively in the Contour panel.
5. Figures 16 and 17 represent the ply shapes and patch locations of +/-45 and 90° ply bundles.
Due to the balance constraint applied, the thickness distribution of the +45° and the -45°
super ply are the same.
Figure 16: Ply Thickness Contour plot of the -45/+45 degree super-plies
Figure 18: Element set 5 representing ply bundle 1 of the +45 degree super ply
10. Repeat step 10.9 to review the element sets 6 though 8 (Figures 19 through 21).
Figure 19: Element set 6 representing ply bundle 2 of the +45 degree super ply
Figure 21: Element set 8 representing ply bundle 4 of the +45 degree super ply
The shapes of the plies as indicated through the element set can be used as-is in design
Phase 2, or modified easily by updating the element sets in HyperMesh to improve the
manufacturability. In this case, the element sets are used as-is.
In the second design phase, a size optimization is performed to fine tune the thicknesses of the
optimized ply bundles from Phase 1. To ensure that the optimization design meets the design
requirements, additional performance criteria on natural frequencies and composite strains are
incorporated into the problem formulation. A normal modes analysis load case is added to
calculate the natural frequencies of the fairing under assembled conditions. The optimization
setup is also modified to factor in these additional performance targets, among others.
The following is the modified optimization setup:
Design Ply thicknesses, which have been defined in the size input deck from Phase 1
variables:
Manufacturing constraints are preserved and transferred to the DCOMP card. A minimum
manufacturable ply thickness of 0.1, defined in Phase 1, is transferred to the PLY card. It allows
for the optimal ply bundle thicknesses to be a multiple of the minimum ply thickness value, and
helps in calculating the total number of plies required per fiber orientation.
Step 1: Load the OptiStruct user profile and import the model
Follow Step 1 in the free-size phase (Phase 1) to load the fairing_freesize.*.fem file into
HyperMesh.
5. Click create.
6. For the response = field, enter freq1.
7. Set frequency as the response type.
8. In the Mode Number field, enter 1 to create the frequency response for the first mode and
make sure that the toggle is set to no regionid.
9. Click create.
10. Repeat step 6.6 through 6.9 to create frequency responses for mode 2, 3, 4, and 5.
11. For response =, enter cstrain.
12. Select composite strain as the response type.
13. Toggle the selection option to plies.
14. Click the highlighted plies and select all plies.
15. Select maj. Principal as the strain type.
The cstrain response should be set, as shown below.
6. Click create.
7. Repeat steps 7.2 through 7.6 to create constraints freq2, freq3, freq4, and freq5
respectively with the same lower bound of 0.02.
8. For constraint=, enter cstrain.
9. Click response= and select cstrain.
10. Activate the upper bound= and enter 0.001.
11. Click loadsteps and select loadcases gravity and pressure.
4. Click create.
5. Click return twice to go back to the Analysis page.
Figure 32: HyperWorks Solver View window to show the progress of the optimization
If the job was completed successfully, new results files can be seen in the same directory
where fairing_size.fem was saved.
Figure 33: Element thickness contour plot (final iteration) after phase-2 size optimization
3. Set the Result type to Orientation Thicknesses (s) to create the thickness contour for
each ply orientation.
This algorithm is aimed at providing a global view of what the optimal stacking sequence could
be. An input deck for the ply stacking sequence optimization, fairing_size_shuffling.*.fem,
was generated from a previous design stage. Each ply bundle is divided into multiple PLYs whose
thickness is equal to the manufacturable thickness (0.1 in this case), and the STACK card is
updated accordingly. In this design phase, composite plies are shuffled to determine the optimal
stacking sequence.
It is important that design performances are preserved. Hence, the optimization problem is
retained as previously formulated in the size optimization phase. Two manufacturing constraints
are applied:
The maximum successive number of plies of a particular orientation does not exceed 4
plies.
The + 45s and – 45s are reversed paired.
Step 1: Load the OptiStruct user profile and import the model
Follow Step 1 from the free-size phase to load the fairing_size_shuffling.*.fem file in
HyperMesh.
6. Click update.
7. Click return.
If the job was completed successfully, new results files can be seen in the same directory
where fairing_shuffling.fem was saved.
The following result files are generated:
The above image shows the history of the shuffling optimization. The columns represent the
global trend of the ply stacking sequence at a particular iteration, with the last column being
the final solution.
Problem Statement
In the model, the non-design regions surround the load-bearing points and without the design
space, the non-design area is discontinuous. This allows the optimization to determine the most
efficient structure to transition the load and stress throughout the structure.
The Model browser shows that there is already an optimization set up.
Lattice optimization differs from topology as a concept-level optimization in that elements with
intermediate densities will be mathematically interpreted as solid, lattice, or void, depending on
your settings in the setup phase. You will add continuation cards and parameters in a text editor
to change this from standard topology to lattice optimization.
5. Click the open file icon in the File: field. A Select OptiStruct file browser window
opens.
6. Select the controlarm.fem file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
7. Click Open. The location of controlarm.fem displays in the File: field.
8. Click Import.
Since DOPTPRM, LATLB, CHECK was applied, OptiStruct checked the drop-in compliance, due to
the removal of elements with density below LB. The drop is relatively small and this LB can be
considered a good choice.
Step 2: Review the model to verify the changes to the original model
1. The beginning of the deck shows that the new file has retained your original objective, but
added a new optimization parameter indicating that this size optimization is Phase II of a
lattice structure optimization.
3. Each CBAR element has its own PBAR property definition, listed later in the deck, which
indicates each element can be individually optimized.
5. At the end of the file are the design variable property relationships (DVPREL1), which relate
properties to their associated design variable for the size optimization in the second phase.
These are followed by the responses, the constraints, and the objective.
Note: The volume fraction response in the original model has been replaced by total volume in
the size optimization phase, which will now be minimized and a new response has been
added to constrain the von Mises stress in the lattice region to 200. This was created from
the value in the last field of the LATTICE continuation card in the first phase.
Note: The difference in lattice penalty (lower than 1.8, set by the optimization parameters)
causes the compliance of the final Phase I model to differ from the initial Phase II
model. This compliance difference is also affected by solid elements retained in the
Phase II model, which recover their full density/stiffness. For this reason, post-
processing a lattice optimization requires that you analyze changes to the model
compliance between the Phase I final optimization compliance and the Phase II initial
compliance calculation and again at the end of Phase II.
5. Load the controlarm_lattice_optimized file in the HyperWorks Solver Run Manager.
6. Click Run for verification purposes.
Since the optimization removes CBAR elements of small radius after the last optimization, the
compliance for the last optimized run should be confirmed against the
controlarm_lattice_optimized.fem file, an analysis of the optimized structure provided by
OptiStruct.
7. Compare the final compliance between the controlarm_lattice.out and
controlarm_lattice_optimized.out files.
Detail of the front of the optimized lattice/solid model showing the variation in CBAR radius
6. Using the Contour panel , select CBAR/CBEAM Stresses (ROD) and set Subtype to
NORMAL S1N(A).
7. Click Apply to contour the beam element visualization.
Problem Definition
The composite design optimization methodology presented within this tutorial was developed to
solve very complex composite design optimization problems. The methodology breaks down the
complex composite design optimization problem, which is not solvable by itself, into several
simpler composite design optimization problems, which are solvable by themselves. The
cumulative solution to each of the simpler composite design optimization problems provides a
solution to the complex design optimization problem. This process of breaking down complex
problems into several simpler problems is consistent with the engineering method.
The following set of steps completes the analysis setup of the initial model and provides a
baseline analysis for comparison with the final optimized structure.
Step 1: Load the OptiStruct user profile and open the model into
HyperMesh
1. Launch HyperMesh.
2. Select OptiStruct in the User Profiles dialog 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 OptiStruct.
User Profiles can also be accessed from the Preferences menu on the toolbar.
Step 3: Create a new element set containing the initial shape of each ply
1. Right-click in the Model browser and select Create > Set to create a new set.
2. In the Entity Editor, for Name, enter ply_shape.
3. Set the value for Card Image to SET_ELEM.
4. Click Entity IDs and click the yellow Elements entity selector.
5. Use the yellow elems entity selector to select all of the elements in the model.
6. Click proceed to continue.
By Prop .
This represents each of the plies in the model according to the color of its ply as shown in the
Model browser. If all of the plies in the model are the same color, change the ply colors in the
Model browser so that each is different to help differentiate the plies in the graphics area.
6. Verify that the fields in the create subpanel match the fields shown above and click create.
The design variable free-size is created for the free-size optimization. The manufacturing
constraints on ply percentage and ply balance will be defined next.
7. Go to the composites subpanel. Make sure free-size is selected as the design variable.
8. Click edit and enter the DSIZE panel to define the manufacturing constraints on ply
percentage, ply balance, and ply drop-off.
9. Check the box in front of PLYPCT.
10. Set Ply Percentage Options to BYANG.
11. In the DSIZE_NUMBER_OF_PLYPCT= field, enter the value of 2.
Two PLYPCT continuation lines are added to the DSIZE data entry.
Figure 4: DSIZE data entry fields for the PLYDRP card using the PLYSLP method
5. Click create.
6. Click return twice to go back to the Analysis page.
With this keyword, OptiStruct automatically generates a sizing model after free-size
optimization.
5. Click return twice to go back to the Analysis page.
6. Save the file as oht_opti_ph1.hm.
oht_opti_ph1_s1(2).h3d Hyper 3D binary results file, with the static analysis results
of both subcases.
2. On the page with the results for oht_opti_ph1_des.h3d (open file name is shown in the
upper right corner by default), go to the Contour panel and select the plot options.
3. In the Results browser, select the last iteration in the Load Case and Simulation
Selection drop-down.
4. Click Apply and then XY Top Plane View to view the results in the X-Y plane.
The element thickness results from the free-size optimization are shown in the image below.
The regions indicated in red or in colors tending towards red (from the legend) can be
interpreted as thicker regions, while those in blue or tending towards blue are thinner regions.
2. Select the last iteration in the Load Case and Simulation Selection drop-down.
3. Click Apply.
The thickness distribution of 0 degree super ply is generated. It represents the ply shapes and
patch locations of the 0 degree ply bundles.
4. Repeat step 1 through 3 to create the ply thickness contours for super-ply 2 (45°), ply 3 (-
45°), and ply 4 (90°) by selecting Layers 2, 3 and 4, respectively in the Contour panel.
5. The figures below represent the ply shapes and patch locations of +/-45 and 90° ply bundles.
Due to the balance constraint applied, the thickness distribution of the +45° and the -45°
super ply are the same.
Figure 14: Ply Thickness Contour plot of the -45/+45 degree super-plies
The shapes of the plies as indicated through the element set can be used as-is in design
Phase 2, or modified easily by updating the element sets in HyperMesh or using ply
smoothing to improve the manufacturability. Ply smoothing operations are shown in the next
section.
In the second design phase, a size optimization is performed to fine tune the thicknesses of the
optimized ply bundles from Phase 1. To ensure that the optimization design meets the design
requirements, additional performance criteria on may be incorporated into the problem
formulation. These new criteria will be fiber strain, matrix strain, and mass.
The following is the modified optimization setup:
Design Ply thicknesses, which have been defined in the size input deck from Phase
variables: 1
Manufacturing constraints previously applied are preserved and transferred to the DCOMP card.
Step 3: Edit the manufacturable thickness and initial value of each ply
1. In the Model browser, select all plies as shown below.
Figure 19: A view of the Model browser showing all plies selected
Figure 22: Element thickness contour plot (final iteration) after phase-2 size optimization
This brings us to the third and final phase of the design process in which you try to identify a
proposal for the optimal stacking sequence of the plies.
This algorithm is aimed at providing a global view of what the optimal stacking sequence could
be. An input deck for the ply stacking sequence optimization, oht_opti_ph2_shuffling.*.fem,
was generated from a previous design stage. Each ply bundle is divided into multiple PLYs whose
thickness is equal to the manufacturable thickness (0.01 in this case), and the STACK card is
updated accordingly. In this design phase, composite plies are shuffled to determine the optimal
stacking sequence.
It is important that design performances are preserved. Hence, the optimization problem is
retained as previously formulated in the size optimization phase. Two manufacturing constraints
are applied:
The maximum successive number of plies of a particular orientation does not exceed 4
plies.
The outermost four layers of the layup must be -45, 0, 45, 90.
Step 1: Load the OptiStruct user profile and import the model into a new
session
Follow Step 1 from the free-size phase to load the oht_opti_ph2_shuffling.*.fem file into a
new HyperMesh session.
Step 3: Update the design variable to input the constraints for the
shuffling optimization
1. In the optimization panel, click on composite shuffle to enter the composite shuffling
panel.
2. Click review and select free-size to review the design variable.
3. In the parameters subpanel, click edit to enter the DSHUFFLE card.
4. Check MAXSUCC and set MSUCC to 4.
5. Check COVER and set NUMBER_OF_VANG to 4.
6. In the DSHUFFLE card editor, set VREP to 1, VANG(1) to -45, VANG(2) to 0, VANG(3) to
45, and VANG(4) to 90.
This sets the outermost layer of the shuffling sequence to maintain one repetition of a set of -
45, 0, 45, 90 degree plies. This ply set is applied to both faces of this optimization, since the
laminate is symmetric.
7. Click return to exit the editor and click update to update the design variable.
8. Click return to return to the optimization page.
The above image shows the history of the shuffling optimization. The columns represent the
global trend of the ply stacking sequence at a particular iteration, with the last column being
the final solution.
The weight of the part has not been changed during the shuffling design phase, rather the
plies were reordered to obtain the maximum performance.
This light weight design therefore meets all of the performance requirements, is feasible and
manufacturable.
The following set of steps provides a final optimized structure for comparison with the initial
model.
Step 1: Load the OptiStruct user profile and open the model into
HyperMesh
1. Launch a new session of HyperMesh.
2. Select OptiStruct in the User Profiles dialog 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 OptiStruct.
User Profiles can also be accessed from the Preferences menu on the toolbar.
Step 2: Import the final optimization properties over the existing model
1. In the Import tab, browse for the oht_opti_ph3.prop file.
2. Expand the Import options section by using the drop-down toggle .
3. Check the box for FE overwrite and click Import to update the properties of the model with
the optimized parameters.
By Prop .
This represents each of the plies in the model according to the color of its ply as shown in the
Model browser. If all of the plies in the model are the same color, change the ply colors in the
Model browser so that each is different to help differentiate the plies in the graphics area.
Figure 24: Visualization of model thickness in HyperMesh with half of the elements in the model masked
Figure 25: A plot of maximum Normal Y strain across all plies, showing that no element exceeds 0.007 me
The structural model, shown above, is loaded into HyperMesh. The constraints, loads, material
properties, and subcases (loadsteps) are already defined in the model. Size design variables and
optimization parameters are defined, and OptiStruct determines the optimal gauges for the
components. The results are then reviewed in HyperView.
The optimization problem for this tutorial is stated as:
Constraints: A given maximum nodal displacement at the loading grid point for
two loading conditions.
Step 1: Launch the HyperMesh Desktop, set the User Profile and Retrieve
the File
1. Launch HyperMesh Desktop.
2. Choose the OptiStruct in the User Profiles dialog 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.
User Profiles GUI can also be accessed from the Preferences menu on the toolbar.
joint_sizeOPT.prop OptiStruct property output file containing all updated property data
from the last iteration for size optimization.
Displacement and stress results are output by OptiStruct (by default) for linear static analyses.
This section describes how to view those results in HyperView. Size optimization results from
OptiStruct are given in the .h3d files and joint_sizeOPT.out.
joint_sizeOPT_s1.h3d Contains displacement and stress results for the linear static
analysis for iteration 0 and iteration 4 of subcase with ID 1
(subcase Force_X).
joint_sizeOPT_s2.h3d Contains displacement and stress results for the linear static
analysis for iteration 0 and iteration 4 of subcase with ID 2
(subcase Force_Z).
The results contained in the HyperView binary results file will be examined first. Then the gauge
history in the joint_sizeOPT.out file will be reviewed.
8. Click Apply.
A contoured image representing shell thickness should be visible. Each element in the model
is assigned a legend color, indicating the thickness value for that element for the current
iteration.
The x-displacement value for 3143 (center of rigid spider, where loading is applied) is shown
in the graphic area. The x-displacement is larger than the upper bound constraint, which was
defined earlier, of 0.9.
14. Click the Next Page icon again to move to the fourth page.
The fourth page shows results loaded from the joint_sizeOPT_s2.h3d file. The name of the
page is displayed as Subcase 2 – Force_Z to indicate that the results correspond to subcase
2.
15. Click the Contour toolbar icon .
16. Select the first drop-down menu below Result type: and select Displacement [v].
17. Select the second drop-down menu and select Z.
18. Click Apply. The resulting contours represent the z component displacement field resulting
from the applied loads and boundary conditions.
19. Repeat steps 8 through 14 to measure and display the z-displacement value for node 3143.
From the UNIX or MSDOS shell, open the joint_sizeOPT.out file in a text editor. Review all five
iterations, noting the volume, constraint information, and gauge at each iteration.
Has the volume been minimized for the given constraints?
Have the displacement constraints been met?
What are the resulting gauges for the rail and tube?
The structural model, as shown, is loaded into HyperMesh. The constraints, loads, material
properties, and subcases (loadsteps) are already defined in the model. Size design variables and
optimization parameters are defined and the OptiStruct software determines the optimal gauges.
The results are then reviewed in HyperMesh.
The optimization problem is stated as:
Constraints: Maximum von Mises Stress of the brackets < 100 Mpa.
Step 1: Launch the HyperMesh Desktop, Set the User Profile and Retrieve
the Database File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK. This loads the OptiStruct user
profile. It includes the OptiStruct template, macro menu, and import reader, paring down the
functionality of HyperMesh to what is relevant for use with OptiStruct.
User Profiles can also be accessed from the Preferences menu on the toolbar.
bracket_size.hgdata HyperGraph file containing data for the objective function, percent
constraint violations and constraint for each iteration.
bracket_size. prop OptiStruct property output file containing all updated property data
from the last iteration for size optimization.
bracket_size.sh Shape file for the final iteration. It contains the material density,
void size parameters and void orientation angle for each element
in the analysis. The .sh file may be used to restart a run and, if
necessary, run OSSmooth files for topology optimization.
Size optimization results from OptiStruct are given in two places. The bracket_size.out file
contains gauge and volume information for all iterations. The bracket_size.h3d file contains the
element thickness for all five iterations and Displacement and Stress results for the linear static
analysis for iteration 0 and iteration 3. This section describes how to review all results. The
results contained in the HyperMesh binary results file will be examined first. The gauge history in
the bracket_size.out file will then be reviewed.
This will now contour your final iteration of that loadcase. Review the stress to see that it is
under the proper constraints.
The .out file contains a summary of the optimization process. From the information in the .out
file, you can see how the objective, constraints, and design variables are changing from one
iteration, to the next.
Has the volume been minimized for the given constraints?
Have the stress constraints been met?
What are the resulting gauges for the two parts?
Did the design variable linking work?
Step 1: Launch the HyperMesh Desktop and Set the User Profile
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog 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.
Step 2: Create Load Collectors for the Loads and Boundary Conditions
1. In the Model browser, right-click and click Create > Load Collector.
2. For Name, enter crank.
3. Click Color and select any color.
4. For Card Image, select None.
This creates a new load collector, crank.
5. Similarly, create another load collector named spcs.
Maximum constraint violation (% [disp > 1.8 mm]) for each iteration.
The next pages show the Design Variables (DVs) which are grouped together making it
possible to compare the behavior of the different plies. This plot can be created by opening
the bicycle_frameOPT.hgdata file.
The structural model, as shown in the figure, is loaded into HyperMesh. The constraints, loads,
material properties, and subcases (loadsteps) are already defined in the model. Size design
variables and optimization parameters are defined, and OptiStruct determines the optimal
gauges. The results are then reviewed in HyperView.
The optimization problem is stated as:
Constraints: Maximum von Mises stress of the brackets < 120 MPa.
Step 1: Launch the HyperMesh Desktop, Set the User Profile, and
Retrieve the Database File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK.
User Profiles can also be accessed from the Preferences menu on the toolbar.
3. Click the Next Page toolbar icon to move to the third page.
The third page has the results loaded from the discrete_bracket_size_s1.h3d file. The
name of the page is displayed as Subcase 1 – STEP to indicate that the results correspond to
subcase 1.
4. Click the Contour toolbar icon .
5. Select the first drop-down menu below Result type: and select Element Stresses [2D &
3D] (t).
6. From the second drop-down menu, select vonMises.
7. Select None in the field below Averaging method:.
8. 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. If you did not change the Iteration
step you should be contouring the stress of the initial step. To contour the final step, set the
last iteration of that loadcase using the icon on the animation toolbar, or use the Results
browser, as follows.
9. Use the bottom drop-down menus in the Results browser to select the last Iteration # in the
Simulation list.
Review
The .out file contains a summary of the optimization process. From the information in the .out
file, you can see how the objective, constraints, and design variables are changing from one
iteration to the next.
Has the volume been minimized for the given constraints?
Have the stress constraints been met?
What are the resulting gauges for the two parts?
Hints
Go to the des.h3d page, clear the contour if one was applied, set to the last simulation step and
apply the Element Thickness contour.
Append discrete_bracket_size.mvw to review objective, constraints, and other information.
Problem Statement
Design variables: Beam width, beam thickness, beam depth, and shell thickness
Step 1: Launch the HyperMesh Desktop, Set the User Profile, and Import
an FE Model File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles 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.
4. From the File menu on the toolbar, select Import > Solver Deck. An Import tab is added to
your tab menu.
7. Click the Select Files icon and select the shredder.fem file you saved to your working
directory from the optistruct.zip file. Refer to Accessing the Model Files.
8. Click Open.
9. Click Import.
12. Click the contour icon and click Apply to plot the Eigen Mode contour.
13. Click the Note icon and remove the first two lines, leaving the rest as below, and click
Apply.
14. Click the Page Window Layout icon to choose the four window layout .
19. Click the play icon in the animation controls to start the animation, and click again to stop
the animation.
The 3rd and 4th mode (~ 3.9 and 4.8 Hz) has a transversal shape that can reduce the
performance of the shredder when it gets excited. The objective, then, is to get the minimum
mass to greater than 7Hz.
20. From the File menu, click Save as > Report Template.
21. Change Save as type: to Report definition (*tpl).
22. For File name:, enter report.tpl.
Beamwide 50 30 90
Beamthck 10 5 15
1
Beamthck 20 15 30
2
5. From the Model browser, select the last Iteration # and click OK, as shown in the image
below.
It is assumed that the tail is cantilevered about its inboard section. Three loading scenarios are
considered; one where the tail experiences pressure loads of 0.25 psi on the bottom skin, a
second where the tail experiences a tip load of 400 lbs, and a third where the tail experiences
both the pressure load and tip load simultaneously. The applied loading is represented in the
following figure.
The materials available for this part are described in the following table. The optimum design
should be as light as possible without failing or buckling under the given loading conditions.
G12 800ksi (800000 psi) 3ksi (3000 psi) Rho 0.1 lb/in3
G1,Z 800ksi (800000 psi) 4ksi (4000 psi) Yield 50ksi (50000 psi)
Step 1: Launch HyperMesh Desktop, Set the User Profile, and Retrieve
the File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog 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.
3. From the File on the main menu, select Open > Model.
4. Select the tail_baseline.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
5. Click Open.
A property collector named Ribs has been created. It has a PSHELL definition with a thickness
of 1.0. It also references the Aluminum 2024-T3 material definition and the component name
Ribs.
The definition of the Inboard_section_top laminate is now complete. The following figure
shows the laminate as displayed in the right-hand side Review panel.
Inboard_section Laminate
12. Right-click on Tail and select Isolate Only. Only the elements forming the ribs which are in
the tail collector should now be displayed.
13. From Tool page, select the organize panel.
14. Click elems and select displayed.
15. Click Destination = and select Ribs from the list of component collectors.
16. Click Move.
17. In the Model browser, right-click on Component and select Show.
18. Click return to return to the main menu.
19. Press F2 on the keyboard.
20. Set the entity selection to comps.
21. Click preview empty and delete entity to clear any empty components (the tail component
in this case).
22. Click Return.
Three loading scenarios are to be considered in this exercise: one where the tail experiences
pressure loads on the bottom skin, a second where the tail experiences a tip load, and a third
where the tail experiences both the pressure load and tip load simultaneously.
In previous steps, a load collector containing the pressure loads and another containing the tip
load were created, but a load collector containing both together is still needed. Next is to create a
load collector which is a combination of the load collectors pressure and tip_load.
Step 10: Request Stress, Strain, and Failure Results for Composite
Laminates
Stress, strain, and failure results are not output by default for composite laminates, but need to
be requested.
1. In the Model browser, right-click on the Outboard_section_top property and select Card
Edit.
The PCOMP card image for the Outboard_section_top laminate appears in the lower portion of
the display area. For more details on the PCOMP card image, refer to the OptiStruct online
documentation.
2. If HILL does not appear beneath [FT], click [FT] notice HILL appears beneath.
This activates failure theory calculation. If HILL is selected, a list of other failure theories
appear - use the Hill failure criteria for this exercise.
3. Click [SB] in the card image window and enter 3,500 in the field beneath it.
This is the interlaminate shear strength of the laminate, which is the bonding material shear
strength. 3.5ksi is an assumed value, as no material data was provided.
4. Click the button beneath SOUT(1) and select YES from the pop-up menu. This requests
stress and strain results to be output for ply1.
5. Set all other plies, i.e. SOUT(2), SOUT(3), SOUT(4) and SOUT(5) to YES also.
6. Click return to keep the changes you made to the card image.
7. Repeat steps 3 through 9 for the other composite laminates.
ANALYSIS RESULTS :
------------------
ITERATION 0
After calculating the failure indices for individual plies, OptiStruct calculates the potential
failure index for the composite shell element. This is based on the premise that failure of a
single layer qualifies as failure of the composite. Thus, a failure index for composite elements
is calculated as a maximum of all computed ply and bonding failure indices (note that only
plies with requested stress output are taken into account here).
15. Change Entity with Layers: to Max to have the maximum index for the laminate.
Repeat this process to have the maximum failure index for all loadsteps.
MAX FAILURE INDEX = 3.73 e-3 (Combo Loadstep)
Now return to HyperMesh to set up the optimization problem. The first step in this process is to
define the design variables. The design variables for this exercise are the rib thicknesses and the
laminates used in the composite skins.
HyperMesh Desktop allows you to use one HyperMesh page and multiple pages from the
HyperView, HyperGraph, MotionView, and MediaView clients without having to switch
applications. To delete the HyperView page and return to the HyperMesh client, click the Delete
Page icon . To keep the page open but return to the HyperMesh client page, click the Previous
Page or Next Page icons until the HyperMesh client returns.
9. Click create.
10. Click return twice to go to the main page.
Twelve total composite design variables now exist, one for the thickness of the glass fabric for
each composite laminate component, and the other for the thickness of the core for each
composite laminate component. As the laminates are symmetric, the glass fabric will
reference the same design variables on either side of the core.
Step 18: Create the Mass, Static Stress and Composite Failure Responses
1. From the Analysis page, select the optimization panel.
2. Select responses.
3. Click response= and enter mass.
4. Click the response type: switch and select mass from the pop-up menu.
5. Click create. The optimization response mass, which is the total mass of the structure, is
created.
6. Click response= and enter vm_strs.
7. Click the response type: switch and select static stress from the pop-up menu.
8. Click props and select the Ribs collector.
9. Click Select. A new selector switch appears next to comps.
10. Make sure that the switch is set to von Mises.
11. Click the switch below von Mises and set it to both surfaces.
12. Click create. The optimization response vm_strs, which is the von Mises stress for the
metallic ribs, is created.
13. Click response= and enter hl_ist.
14. Click the response type: switch and select composite failure from the pop-up menu.
15. Click props and select the Inboard_section_top collector.
16. Click the switch next to props and select hill.
17. Click the switch below hill and select all plies.
18. Click create. The optimization response hl_ist is created.
This is the hill failure criteria for all plies of the composite skins of the Inboard_section_top
component collector.
19. Repeat steps 12 through 18 to create optimization responses for the hill failure criteria for the
plies of the other composite laminate skins. The responses could be similarly named: hl_osb,
hl_ost, hl_msb, hl_mst, and hl_isb.
2. Change the current client to HyperView using the client selector drop-down .
3. Click the down arrow to the right of the Select application toolbar icon and select
HyperView .
4. Click Yes to continue.
5. Create a two pane view by clicking on the Page Window Layout icon and selecting the two
view icon from the pop-up menu.
4. Click the file folder icon at the end of the File: field and select the tail_opt.fem file from
the directory where you ran the optimization.
5. Click Import. This loads the *.fem that the optimization was run with into HyperMesh.
Conclusion
The objective of this tutorial was to achieve the lightest design by varying the laminate properties
and rib thicknesses. Experimenting with other materials and other laminate configurations could
lead to a lighter design.
There are four shell components in the model: the mounting flange, the web, the top and bottom
flanges, and the lug. The web is connected to the lug by gap elements. Appropriate properties,
loads, boundary conditions, and nonlinear subcases have already been defined in the model. The
design region is the web and the rest of the components are non-design. Since a large portion of
aerospace components are shell structures which are manufactured by machining or milling
operations, free-sizing optimization is very suitable for those components. To understand the
limitations of topology optimization for such applications, a nonlinear gap topology optimization
will also be done on the wing rib model.
Design variables for free Thickness of each shell element in the design space.
sizing optimization:
Design variables for Element density of each element in the design domain.
topology optimization:
Exercise
Step 1: Launch HyperMesh Desktop, Set the User Profile, and Retrieve
the File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK.
3. From the File menu on the toolbar, select Open Model. An Open Model browser window
opens.
4. Select the rib_complete.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
5. Click Open.
rib_freesize_hist.mvw This file is a HypeView session file and may be opened from
the File menu in HyperView or HyperGraph. The file
automatically creates individual plots for each of the results
(objectives, constraints) contained in the .hist file. Each
plot occupies its own page within HyperView (HyperGraph).
rib_freesize_frame.html The file contains two frames. The top frame opens one of
the .h3d files using the HyperView Player browser plug-in.
The .h3d file opened depends on the results selected for
display in the bottom frame. The bottom frame opens the
_menu.html file, which facilitates the selection of results to
be displayed.
rib_freesize.oss The file contains default settings for running OSSmooth after
a successful optimization.
rib_freesize.sh Shape file for the final iteration. The .sh file may be used to
restart a run.
rib_freesize.fsthick The element definitions for those elements that were part of
a free size design space. The optimized thickness of these
elements is provided as nodal thickness values (Ti).
rib_freesize.mvw This file is a HypeView session file and may be opened from
the File menu in HyperView. The file automatically creates
individual load the optimization results (dens.h3d) and the
loadstep results (s#.h3d).
Element thickness distributions are output from OptiStruct for all iterations. In addition,
Displacement and Stress results are output for each subcase for the first and last iteration by
default. This section describes how to view those results in HyperView.
1. From the OptiStruct panel, click HyperView.
This should open a new window with the HyperView client and load the rib_freesize.h3d,
reading the model and the results.
2. Click close to close the message window.
3. Click the Entity Attributes icon on the toolbar and undisplay all of the components,
except Web. This is accomplished by activating the Auto apply mode: (Display OFF) and
then clicking on the component that you want turned off in the GUI.
4. Click the Mesh:, shaded mesh option .
5. Click the Web component to get a shaded mesh.
6. Go to the Contour panel and set the Result type: to Element Thicknesses.
7. In the loadcase selection area above the Results browser, select the last iteration listed in the
Simulation list and click OK.
Thickness contour from free-sizing nonlinear gap optimization, on the Web of plate thickness 0.1mm
rib_topology_frame.html The file contains two frames. The top frame opens one of
the .h3d files using the HyperView Player browser plug-in.
The .h3d file opened depends on the results selected for
display in the bottom frame. The bottom frame opens the
_menu.html file, which facilitates the selection of results to
be displayed.
rib_topology_menu.html This file facilitates the selection of the appropriate .h3d file
for the HyperView Player browser plug-in in the top frame
of the _frames.html file, based on chosen results.
rib_topology.sh Shape file for the final iteration. It contains the material
density, void size parameters, and void orientation angle
for each element in the analysis. The .sh file may be used
to restart a run.
Element density results are output from OptiStruct for all iterations. In addition, displacement
and stress results are output for each subcase for the first and last iteration by default. This
section describes how to view those results in HyperView.
1. From the OptiStruct panel, click HyperView. This opens new pages with the HyperView client
and loads the session file, rib_topology.mvw that is linked with .h3d files where the model
and results are defined.
2. Click close to close the message window.
3. Click the Entity Attributes icon on the toolbar and undisplay all of the components, except
the Web component. This is accomplished by activating the Auto apply mode: (to Display
Off) and clicking on the components that you want turned off in the GUI.
4. Click the Mesh: panel shaded mesh option.
5. Click the Web component to get a shaded mesh.
6. Go to the Contour panel and set the Result type: to Element Densities.
7. Click in the bottom right portion of the GUI to activate the Load Case and Simulation
Selection dialog.
8. Select the last iteration listed in the Simulation list and click OK.
9. Click Top in the view controls to get a top view of the Web.
Contour of element density on the Web component from topology nonlinear gap optimization
The results from topology optimization show very discrete results as expected.
Results from the topology optimization show a truss type design with extensive cavities and
voids, while the results from free-sizing optimization tend to come up with shear panels. While
solid/void density distribution is the only choice for solid elements; for shell structures,
intermediate densities can be interpreted as different thicknesses and penalizing then could result
in potentially inefficient shell structures. Moreover, since a large portion of aerospace structures
are shell structures, a shear panel type design is often desirable for manufacturing purposes
especially for machine milled shell structures. Free-sizing optimization can prove to be very
beneficial in those situations.
Rotating bar
Exercise
Step 1: Launch HyperMesh Desktop, Set the User Profile and Retrieve the
File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK.
3. From the File menu on the toolbar, select Open > Model.
4. Select the rotating_bar_design.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
5. Click Open.
9. Make sure all dofs (dof1 to dof6) are checked and their values are set to 0.00.
10. Click create.
11. Click return.
12. Click create. A response, Stress for all the bar elements, is defined.
13. Click return to go to the optimization panel.
One of the advantages with composite materials is that the structural performance can be
controlled precisely by choosing the appropriate ply thickness, ply orientation, stacking sequence,
ply materials, and so on. The ability to vary many different parameters provides greater
flexibility, but at the same time it is tougher to optimize the part as the number of design
variables increases many fold. OptiStruct has the ability to directly or indirectly optimize the ply
thickness, ply orientation and stacking sequence for composite structures.
Free-size optimization handles the thickness of each ply in each element as a design variable and
optimizes the structure by determining the optimal thickness distribution for each ply in the
laminate. Refer to the Free-Size Optimization and Composite Topology and Free-size Optimization
pages of the OptiStruct User's Guide help for more details.
For several reasons, every composite manufacturer has their own manufacturability standards for
the laminated composites. These additional manufacturing constraints are to be included with
free-size optimization to achieve an acceptable manufacturing solution. OptiStruct supports
different manufacturability constraints that can be defined with free-size optimization. This
tutorial helps explain the procedure used to define the manufacturing constraints in the free-size
optimization of composite structures.
The optimization problem for this tutorial is stated as:
Exercise
Step 1: Launch HyperMesh Desktop, Load the User Profile and Retrieve
the File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog 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.
User Profiles GUI can also be accessed from the Preferences menu on the toolbar.
3. From the File menu on the toolbar, select Open.
4. Select the Composite_Wing.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
5. Click Open to load the model into the HyperMesh session.
1. From the 2D page, select the HyperLaminate panel. This opens the HyperLaminate GUI in
which the ply lay-up information can be defined, reviewed and edited. Size design variables
can also be set up in this panel for performing size optimization.
2. Expand the Laminates portion of the tree structure on the left-hand side of the screen.
3. Select the TopSkin component for review by clicking on the component name. This loads the
TopSkin component properties in the Laminate definition section and Review section. The
Laminate definition section shows the ply material, ply thickness, ply orientation, etc. which is
shown graphically under the Review section.
The same layup information is also defined for the BottomSkin component. In this tutorial,
you will define free-size optimization on the TopSkin and BottomSkin components.
4. From the File menu, select Exit.
Exiting the HyperLaminate GUI and returning to HyperMesh.
10. In the first field next to COMP PLYPCT and below PANGLE, enter 0 to define that ply
percentage constraints are defined for the ply with ply angle of 0 degree. Enter 0.1 for
PPMIN and 0.6 for PPMAX. The PTMAN field can be left blank.
The above step defines that for each element, the thickness of the ply with ply angle 0,
should be no less than 10% or more than 60% of the total thickness.
11. Similarly define the same constraints of 0.1 for PPMIN and 0.6 for PPMAX for the plies with
ply angle 45, -45 and 90 as shown in the image above.
12. Similarly, activate BALANCE by clicking on the box, and define 1 in the field for
DSIZE_NUMBER_OF_BALANCE=.
The BALANCE constraint ensures that two plies will always be of equal thickness.
BALA BALANCE flag indicating that a balancing constraint is applied.
NCE
BANG First ply orientation to which the BALANCE constraint is applied.
LE1
BANG Second ply orientation to which the BALANCE constraint is applied.
LE2
13. In the field below BANGLE1, enter 45 and in the field below BANGLE2 enter -45.
The above step defines that the plies with ply angle of 45 and -45 will always have the same
thickness.
14. Click return to return from the panel.
15. Click update to update the above defined manufacturing constraints to the free-size design
variable.
16. Click return to return from the free size panel.
OptiStruct provides element thickness, ply thickness information for all iterations, and also writes
out displacement and von Mises stress results for the linear static analysis. This section describes
how to view the results in HyperView.
3. Click the Isolate Shown icon and then click BottomSkin and TopSkin in the component
tree.
This leaves only these two components displayed.
9. From the Results menu, activate the menu over Plot and click Contour. Choose Ply
Thickness as the Result type and leave Thickness selected.
10. For Entity with Layer:, select 1, and click Apply.
This shows the first ply thickness contour. You can repeat these steps to plot the thickness for
Ply 2, Ply 3, and Ply 4 or Max, etc.
11. Verify if all the manufacturing constraints (ply percentage, balance and minimum laminate
thickness) are satisfied. Additionally, open the Wing_FreeSize_with_PLYPCT.out file in a text
editor and verify that the displacement constraints are satisfied in the last iteration.
The structural model shown above is loaded to HyperMesh. The materials, shell properties, loads
and boundary conditions are already defined in this model. The thicknesses of the three
components are identified as design variables. The von Mises stress of element numbers 58 and
59 (elements located on the circumference of the hole) are defined as responses, and a total
volume response is defined as well.
The von Mises stress of elements 58 and 59 are passed as inputs to the HyperMath script which
in turn, returns two values: the sum of the two von Mises stresses, and the average value of the
two elemental von Mises stresses.
Constraints: Constraints on the sum of the von Mises stresses and the average von Mises
stress
Exercise
Step 1: Launch HyperMath and HyperMesh Desktop, Set User Profile, and
Retrieve the File
1. Launch HyperMath.
2. Under the File menu, select open and select the dresp3_simple_h.hml file you saved to
your working directory from the optistruct.zip file. Refer to Accessing the Model Files.
Examine the HyperMath script to calculate external responses.
4. Click return.
5. Click prev.
Note: A ‘$’ symbol indicates a comment and the following data will
not be read by the solver.
This defines two external responses: the sum of the von Mises stresses of elements 58 and 59
(SUMH) and the average von Mises stress for elements 58 and 59 (AVGH).
The DRESP3 responses have different IDs from the DRESP1 responses and point to the library
called HLIB defined in Step 2. Also, the function MYSUM is the same function name in the
dresp3_simple_h.hml script.
This completes linking of the DRESP3 cards with the HyperMath Script.
The structural model, shown in Figure 1, is loaded into HyperMesh and is used to generate and
run a shape optimization of the cantilever beam. Shape perturbation vectors are generated using
HyperMorph, which is accessed, through the HyperMesh interface. The OptiStruct software
determines the optimal shape. The results are then reviewed in HyperView.
The optimization problem for this tutorial is stated as:
Constraints: Given maximum nodal displacement at the end of the beam < 3.0
mm.
Step 1: Load the OptiStruct User Profile and Retrieve the Geometry
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog 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.
User Profiles can also be accessed from the Preferences menu on the toolbar.
3. From the File menu, select Open > Model.
4. Select the beamshape.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
5. Click Open. The beamshape.hm is loaded into the current HyperMesh session, replacing any
existing data. The database only contains geometric data.
The location of beamshape.hm displays in the file: field.
Figure 4. New handle added at the center of the top edge of the beam.
40. Select the set biasing subpanel using the radio buttons on the left-hand side of the panel.
41. Click handles.
49. Select the save shape subpanel using the radio buttons on the left-hand side of the panel.
50. Click shape = and enter shape2.
51. Select a color or leave the default.
52. Set the toggle under shape= to as node perturbations.
53. Click save.
54. Click Yes in the dialog box that prompts if you want to Save perturbations for nodes at
global and morph volume handles.
This shape is now saved, later it can be associate to a design variable.
55. Click undo all. The model returns to its original shape.
56. Click return twice to return to the OptiStruct panel.
This section describes how to view results in HyperView, which will be launched from within the
OptiStruct panel of HyperMesh.
HyperView is a comprehensive post-processing and visualization environment for finite element
analysis (FEA), multi-body system simulation, video and engineering data.
4. Select last iteration from Simulation using the drop-down selection, as seen in the figure
below.
Step 10: View a Contour Plot of the Displacement on Top of the Shape
Optimized Model
1. Click the Next Page arrow to move to page 2 using the page navigation arrow icon.
2. Select the last iteration from Simulation using the drop-down selection, as seen in the image
below.
3. Go to the Deformed panel and set Result type: to Shape Change (v).
4. Click Apply. The optimized shape of the beam is shown on the screen now.
5. Go to the Contour panel and set Result type: to Displacement (v).
6. Select Mag using the toggle.
Exercise
Step 1: Launch the HyperMesh Desktop, Set the User Profile and Retrieve
the File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog 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.
3. From the File menu on the toolbar, select Open > Model.
4. Select the Lbeamshape.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
5. Click Open.
23. Repeat steps 9 through 22 for the local handles 3, 4 and 5. Referring to the figure above,
translate handles 3 and 4 by x=-10 and handle 5 by y=-10. Save the shapes after morphing
each handle as shape2, shape3 and shape4, respectively.
24. Click return twice to go to the optimization panel.
7. Select the response node by clicking on the node, shown in the above figure.
Shape contour information is output from OptiStruct for all iterations. In addition, displacement
and stress results are output for the first and last iterations by default. This section describes
how to view those results in HyperView.
4. After reviewing the animation, click the Stop icon to stop the animation and move Current
time: back to 0.
6. Click Apply.
A plot of the displacements on your final shape should be displayed. The maximum
displacements for the last Iteration #, is still below 2.0.
Constraints: No Constraints
5. Click create.
6. Click return to go to the main menu.
This section describes how to view the results in HyperView which will be launched from within
the OptiStruct panel of HyperMesh.
HyperView is a complete post-processing and visualization environment for finite element
analysis (FEA), multi-body system simulation, video and engineering data.
4. Select the last iteration (Iteration 8) using the drop-down selection, seen in the image below.
2. Select the last iteration (iteration 8) using the drop-down selection, as seen in the image
below.
3. Go to the Contour panel and select Element Stresses [2D & 3D] as the Result type:.
4. For stress type, select von Mises.
In the previous run, no constraints were applied on the movement of the DSHAPE grids.
Therefore, grids are free to move and the part thickness increases, as shown in the figure.
In practice, however, there will be some sort of constraints imposed upon the movement of grids
due to manufacturability. For this tutorial model, thickness must be unchanged to avoid any
interference with other parts.
The constraints on the curved part will be created using a local rectangular coordinate system
(the other constraints on the flat part do not need a local coordinate system). Therefore, a local
rectangular coordinate system (z-axis will point to normal to DSHAPE surface) needs to be
created first.
1. Use the page navigation controls to move back to Page 1 and the HyperMesh client.
2. Back in HyperMesh, click return and go to the 1D page.
3. Click systems.
4. Select the create by axis direction subpanel.
5. Click nodes and select node ID 20999 (See the following figure).
6. Click origin and select the same node (ID 20999) as nodes.
7. Click x-axis and select node ID 15989.
8. Click xy-plane and select node ID 19462.
9. Click create.
10. Click return.
11. From the Analysis page, click optimization.
12. Select the free shape panel.
13. Select the gridcon subpanel.
The constraints on the flat part will be created first without any coordinate system.
14. Click desvar= and select shape.
15. Select constraint type as planar.
16. Select nodes, shown in the following figure.
Follow the previously described steps on how to post-process the results (optimization results
without constraints) using HyperView, and compare the final shape change and stress results.
Problem Statement
This problem will perform a size and shape optimization on a structural rail to prevent buckling in
the rail structure.
Step 1: Launch the HyperMesh Desktop, Set the User Profile, and
Retrieve the File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog and click OK.
3. From the File menu on the toolbar, select Open > Model.
4. Select the os_buckle_original.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
5. Click Open. The os_buckle_original.hm database is loaded into the current HyperMesh
session, replacing any existing data.
OptiStruct will search for the three lowest eigenvalues below 100.
4. Click Analysis type and select Linear buckling from the drop-down menu.
5. For STATSUB, select the pre-existing LINEAR subcase.
6. For METHOD(STRUCT), select buckling from the list of load collectors.
4. Click return.
9. With the animation running, use the slider bar to adjust the speed of the animation.
10. Click the Animation Controls icon next to the seek slider and activate the Bounce option to
review a back and forth animation of the impact.
11. Stop the animation and use the Current time: slider bar to manually control the animation.
For complete details on the various options available in this panel, refer to the online help.
Rail joint
Problem Statement
Constraint: Maximum von Mises stress of the joint < 200 MPa
Step 1: Launch HyperMesh Desktop, Set the User Profile and Retrieve the
File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog 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.
User Profiles GUI can also be accessed from the Preferences menu on the toolbar.
3. From the File menu on the toolbar, select Open > Model.
4. Select the rail_joint_original.hm file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
7. Take note of the Maximum von Mises Stress of the joint and use the Delete Page icon to
delete the page with the HyperView client.
This should return you to Page 1 and the HyperMesh client.
8. Back in HyperMesh, click return to exit the panel.
Rail domains
Step 6: Split the Circular Edge Domains Around the Opening of the Rail
The following steps show the procedure to split each of the two circular domains (as seen in the
previous figure) into four curved edge domains.
19. Under edges only, click reset . This will clean up previous selection from buffer.
20. Click domains and select the red edge curves, as shown the following figure.
21. Click morph. A new curvature is applied to the selected eight edge domains. See the
following figure below.
22. Select the radio button save shape.
23. Click shape =, enter the name sh2.
29. Select the radio button apply shapes. In HyperMorph, a new shape can be created as a
linear combination of existing shapes.
30. Click shapes and select both sh1 and sh2.
31. Click Select.
32. Verify that the multiplier is 1.0.
33. Click apply.
34. Select the radio button save shapes.
35. Click shape = and enter sh3.
36. Make sure that the toggle is set to node perturbations. The new shape sh3 includes
influences from both sh1 and sh2 shapes, as shown in the next figure.
37. Click save.
38. In the Model browser, right-click Shape and select Hide.
Do NOT click undo all at this moment, because one more shape based on this third shape
change will be created.
An additional shape variable is created using the shape created in the previous step.
39. In the Model browser, right-click on the component PSHELL and click Hide. These
components are turned off for ease of visualization.
40. Toggle the radio button to alter dimensions.
41. Under edges only, click reset . This will clean up previous selection from buffer.
42. Switch the top selector from curve ratio to distance =. This feature allows you to shorten
the distance between selected domains.
43. Switch the end a: selector from two handles to nodes and handles.
44. Click node a and pick node, as shown in the next figure.
45. Click node b and pick node, as shown in the next figure.
Once nodes a and b are selected, the distance between node a and node b is measured
automatically and appears in distance = field.
The distance between node a and node b is about 43.
46. Click handles under node a and select the 8 handles shown by the downward pointing
arrows in the previous figure.
To select, click the handles on the screen until they are highlighted.
47. Click handles under node b and similarly as in the previous step, select the 8 handles near
the opposite face of the tube.
48. Toggle the bottom selector and select hold middle. These components are turned on for ease
of visualization
49. Click distance = and enter 20.
50. Click morph. A rectangular shape appears to the joint, as shown in the next figure.
51. Toggle the button to save shape.
52. Click shape = and enter sh4.
53. Make sure that the toggle is set to node perturbations.
54. Click save.
55. Click undo all to restore the mesh to the baseline configuration.
56. In the Model browser, right-click Shape and select Hide.
57. Click return three times to go the main menu.
Step 10: Define the Shape Design Variables and Review by Animation
1. From the Analysis page, select the optimization panel.
2. Click the shape panel.
3. Make sure the radio button is set to desvar and create.
4. Toggle the switch to select multiple desvars.
5. Click shape = and select sh1, sh2, sh3 and sh4.
6. Click select.
7. Click initial value = and enter 0.0.
8. Click lower bound = and enter -1.0.
9. Click upper bound = and enter 1.0.
10. Click create.
This creates four design variables with the same initial value, lower bound, and upper bound.
HyperMesh automatically links the design variables to each shape, respectively and assigns
names to each design variable the same as its associated shapes.
11. Click animate.
12. Click simulation = SHAPE – sh1 (1).
13. Make sure that data type = is set to Perturbation vector.
14. Click modal to animate the first shape variable.
This section describes how to view the results in HyperView which will be launched from within
the OptiStruct panel of HyperMesh.
von Mises stress for the last iteration (Max < 200 MPa)
Is your design objective of minimizing the volume obtained? If not, can you explain why?
Are your design constraints satisfied?
Which shape has the most influence in this problem setup?
What is the percentage decrease in compliance?
Can size optimization be introduced to the joint?
4 bar linkage.
Step 1: Launch HyperMesh Desktop, Set the User Profile, and Retrieve
the File
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles 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.
3. From the File menu on the toolbar, select Open > Model.
4. Select the 4bar_design.hm file you saved to your working directory from the optistruct.zip
file. Refer to Accessing the Model Files.
The 4bar_design.hm database is loaded into the current HyperMesh session.
13. In the Measure Groups list, check the box beside Dynamic MinMax Result.
14. Click Max at the list below Resolved in.
15. Click Create Curve, modify Place curve on: to Existing Plot. This will expose a list of plot
windows on this report.
16. Check the Live link box, select Window 2 and click Apply. This will plot the Maximum von
Mises stress (t). Now this file will be saved as a template to be applied on the optimized
result.
17. From the File menu, select Save As > Session and change the Save as type: to Report
definition (*.tpl).
18. Enter Stress_report.tpl click Save.
19. Use the page navigation icons to return to Page 1 and the HyperMesh client.
10. From the graphics area, click the center of the lower spider of the component Left. You should
see two node numbers at one location; choose node number 288.
11. From the graphics area, click the center of the left spider of the component Middle and choose
node number 143.
12. From the graphics area, click the center of the lower spider of the component Right and
choose node number 441.
13. Make sure the Load Type= is set to SPC.
14. Click create.
Open the file 4bar_design_opt.dsvar to see how OptiStruct changed the DVs during the
optimization process.
This will show that all DVs for the right and mid arms went to the limit of 1.2, showing that
minimizing the mass of these two arms are key to reducing the Stress.
Exercise
Step 1: Launch HyperMesh Desktop, Set the User Profile, and Retrieve
the Model
1. Launch HyperMesh Desktop.
2. Choose OptiStruct in the User Profiles dialog 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.
User Profiles can also be accessed from the Preferences menu on the toolbar.
3. From the File menu, select Open > Model.
6. Click create.
7. Click the parameters subpanel and select the direction as grow.
8. Set mvfactor= to 0.5.
9. Set nsmooth to 10 and click update.
This section describes how to view the results in HyperView which will be launched from within
the OptiStruct panel of HyperMesh.
HyperView is a complete post-processing and visualization environment for finite element
analysis (FEA), multi-body system simulation, video and engineering data.
Step 12: View a Contour Plot of the Stress on Top of the Shape
Optimized Model
1. Go to the Contour panel and select Element Stresses (2D & 3D) (t) as the Result
type:.
2. Select von Mises as the stress type.
3. Click Elements, click By set, and select the set stress.
4. Click Add and close.
5. Click Apply.
Figure 2(a). Load time histories for vertical force Figure 2(b). Load time histories for braking force
The three files found in the optistruct.zip file are required to perform this tutorial. Refer to
Accessing the Model Files.
tarm_fatigue.fem
tarm_loadY.csv
tarm_loadX.csv
Exercise
2. Click the Open Model file icon . An Open file browser window opens.
3. Select the tarm_fatigue.fem file you saved to your working directory from the
optistruct.zip file. Refer to Accessing the Model Files.
4. Click Import.
This loads the control arm model. It includes a whole definition of two static subcases,
elements sets, and material static properties, and so on.
5. Click Apply.
This guides you to the next task of the fatigue process tree, Fatigue Subcase.
fatsub is selected as the active fatigue subcase. Definitions in the following processes
(analysis parameters, fatigue elements and properties, loading sequences, etc.) will be for
this subcase.
5. Click Apply.
This saves the current definitions and guides you to the next task, Analysis Parameters, of the
fatigue process tree.
Certainty of 0.5
survival:
4. Click Apply.
This saves the current definitions and guides you to the next task of the fatigue process tree,
Elements and Materials. For details, consult the HyperWorks online help.
8. Click the Show SN curve definition icon . An SN Method description window introducing
how to generate the SN material parameter pops up.
9. Click Close.
This saves the current definitions and guides you to the next task in the fatigue process tree,
Load-Time History.
5. Click the Open load-time file icon . An Open file browser window opens.
6. Browse for tarm_loadY.csv.
7. Click Open.
8. Click Import.
9. Click Plot L-T to show the load-time history.
10. Click Save to write the new load-time history into HyperMesh database.
11. Close the Load Time History window.
12. Follow steps 2 - 11 to create another load-time history, hist_x, by importing the file
tarm_loadX.csv.
13. Click Apply.
This saves the current definitions and guides you to the next task, Loading Sequences, of the
fatigue process tree.
This completes the fatigue problem set-up. Next, you will set-up the optimization problem.
This section describes how to set-up the free-shape design variables, objective, and fatigue life
constraint.
6. Click create.
7. Click return.
2. Click the Save .fem file as icon . A Save file browser window opens.
3. Set the directory in which to save the file, and in File name, enter tarm_fatigue_opti.fem.
4. Click Save to close the window.
5. Click Save to save the OptiStruct model file.
6. From Run Options:, select optimization.
7. Click Submit.
tarm_fatigue_opti_s4.h3d Hyper 3D binary results file, with both static analysis results and
fatigue free-shape optimization results.
This launches HyperView to load the tarm_fatigue_opti_s4.h3d results file for life results on
top of shape results.
3. Click on last iteration (Iteration 11) from the Simulation list of the Load Case and Simulation
Selection Manager.
A regular shape optimization has been defined in the model. The formulation of this optimization
is listed below.
Objective: Minimize the maximum (minmax) displacement at the node where the
excitation load was applied.
Step 1: Load the OptiStruct user profile and import the model
1. Launch HyperMesh Desktop.
2. Select OptiStruct in the User Profiles dialog 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 OptiStruct.
User Profiles can also be accessed from the Preferences menu on the toolbar.
3. From the File menu, select Import > Solver Deck.
4. The File type: is OptiStruct.
5. Click the open file icon in the File: field. A Select OptiStruct file browser window opens.
2. To review the shape design variables, click shape on the optimization panel.
3. Click animate and enter the Deformed Shape panel. One of the shapes should be displayed
in simulation=.
4. Click linear, you will see the animation of that shape.
5. The other shapes can be reviewed by clicking the next or prev icon.
6. Return to the optimization panel.
Model review
Step 1: Load the OptiStruct user profile and import the model
1. Launch HyperMesh Desktop.
2. Select OptiStruct in the User Profiles dialog 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 OptiStruct.
User Profiles can also be accessed from the Preferences menu on the toolbar.
3. From the File menu, select Import > Solver Deck.
4. The File type: is OptiStruct.
5. Click the open file icon in the File: field. A Select OptiStruct file browser window opens.
Downloading If you are using the examples via the Altair website, you will need to
Model Files download the model files before beginning. Access them by clicking:
http://www.altairhyperworks.com/hwhelp/Altair/hw2017/index.aspx
Please note that a User ID and password is required to access this area.
Follow the instructions provided to obtain the login information.
Summary
Altair OptiStruct was used to run a free-free normal modes analysis on a gas turbine compressor
wheel solid FE model.
Model
A solid hexa model of a gas turbine compressor wheel with fins was generated with Altair
HyperMesh. The model contained 44,369 grids and 31,022 solid hexa elements. A free-free
modal analysis was conducted using OptiStruct to find the first twelve natural frequencies.
Mode # Frequency
7 2.46
8 2.46
9 3.33
10 3.62
11 3.62
12 3.83
File Location
<install_directory>/demos/hwsolvers/optistruct/compressor.fem
Summary
Two overlapping plates (aluminum) are connected by a rivet (titanium) forming a lap joint. The
aluminum and titanium materials are both defined by piece-wise linear elasto-plastic law. The
plates and the rivet are meshed with solid elements. The free end of the bottom plate is
constrained and the free end of the top plate is pulled (by applying imposed displacement) to
shear the joint. An all inclusive contact is defined such that all the components in the model are
master and all nodes of the model are slave.
This example is considered a static problem and the nonlinear implicit solver is used.
Title
Lap joint
Number
40.1
Brief Description
A lap joint is fixed at one end and pulled at the other to shear the joint.
Keywords
Nonlinear large displacement analysis (NLGEOM)
Contact definition (CONTACT)
Plasticity and Piece-wise linear elasto-plastic material (MATX36 and TABLES1)
OptiStruct Options
Parameters for Geometric Nonlinear Implicit Static Analysis Control (NLPARMX)
Boundary conditions (SPC)
Imposed displacement (SPCD and NLOAD1)
Solid element (PSOLIDX)
Contact property for NLGEOM analysis (PCONTX)
Overview
The material used for the aluminum plates have the following properties:
Density: 1.2e-9 Mg/mm3
Young’s modulus: 71700 MPa
Poisson’s ratio: 0.33
Yield stress: 350 MPa
The stress vs plastic strain plot for aluminum is shown in Figure 2.
The material used for the titanium rivet has the following properties:
Density: 7.8e-9 Mg/mm3
Young’s modulus: 112000 MPa
Poisson’s ratio: 0.34
Yield stress: 907 MPa
The stress vs plastic strain plot for titanium is shown in Figure 3.
The imposed displacements are defined in FEM file using NLOAD1 card:
………………
NLOAD1 7 3 DISP 8
CONTACT 6 7 2 7 OPENGAP
A small physical gap of around 0.02mm has been introduced between the top and bottom plates
and also between the plates and the rivet. The minimum gap specified (0.022) for the contact is
slightly higher than the physical gap for contact to take effect. A static Coulomb friction of 0.05 is
defined for the interface.
PCONT 7 AUTO
PCONTX 7 0.05 0.022 0
+ CONST
The plasticity and contact causes major nonlinearities; therefore, a static nonlinear analysis is
performed using the arc-length displacement strategy. The time step is determined by a
displacement norm control.
The nonlinear implicit parameters used are:
Tolerance: 0.01
Maximum convergence 20
iteration number:
The input nonlinear implicit options set in the FEM file are defined by NLPARMX:
NLPARM 9 100 P
+ 0.01
+ 6 1.02 20 0.8
Refer to the OptiStruct manual for more details about implicit options.
The nonlinear large deformation analysis has to be defined through a subcase. An NLPARM
statement, as well as ANALYSIS=NLGEOM has to be present in the subcase. The termination time
of 1.0s is defined through the TTERM entry.
SUBCASE 1
ANALYSIS NLGEOM
SPC = 10
NLPARM = 9
NLOAD = 7
TTERM = 1.000
Animations
The displacement, stresses (mises) and plastic strain results after the shearing of the joint are
shown in the following figures.
Four cantilever beams are analyzed. The objective is to evaluate the response of the beams
under (i) linear static analysis (small displacements) and (ii) geometric nonlinear analysis (large
displacements) with and without the application of follower forces.
Considering the example is a static problem, the nonlinear implicit solver is used.
Title
Follower Force
Number
41.1
Brief Description
Cantilever beams.
Keywords
Nonlinear large displacement analysis (NLGEOM)
Termination time (TTERM)
OptiStruct Options
Parameters for Geometric Nonlinear Implicit Static Analysis Control (NLPARMX)
Boundary conditions (SPC)
Applied forces (FORCE, TABLED1, and NLOAD1)
Default shell element parameters (XSHLPRM)
Fixed coordinate system (CORD2R)
Moving coordinate system (CORD3R)
Input File
Follower force: <install_directory>/demos/hwsolvers/optistruct/follow_force.fem
Modeling Methodology
The mesh is a regular shell mesh with average element size of 5mm.
The beams have been modeled using first order reduced integrated shell elements as specified in
the default definition of shell element properties (XSHLPRM):
PSHELL 1 11.0 1 1 0.0
XSHLPRM 24 2 VAR NEWT 5
The loads and boundary conditions applied in the model are shown in Figure 2.
The cantilever beams are constrained at one end and all dof and forces are applied at the other
end (15N are applied at the the two outer nodes as 30N is applied at the center node).
SPC 6 1 1234560.0
………………
NLOAD1 7 4 8
Tolerance: 0.01
Maximum convergence 20
iteration number:
The input nonlinear implicit options set in the FEM file are defined by NLPARMX:
NLPARM 1 2
Refer to the OptiStruct manual for more details about implicit options.
The nonlinear large deformation analysis has to be defined through a subcase. An NLPARM
statement, as well as ANALYSIS=NLGEOM has to be present in the subcase. The termination time
of 1.0s is defined through the TTERM entry. The first subcase is linear static and the second
subcase is geometric nonlinear.
$
SUBCASE 1
SPC = 6
LOAD = 2
$
$HMNAME LOADSTEP 2"nlgeom" 15
$
SUBCASE 2
ANALYSIS NLGEOM
SPC = 6
NLPARM = 1
NLOAD = 7
TTERM = 1.0
Animations
The displacement contour on the beams is shown in Figure 3. As expected, beam 1 for which
static analysis has been performed, shows the largest deformation. Beam 2 (for which loading
has been defined in a fixed coordinate system) and beam 4 (for which loading has been defined
in the default global coordinate system) show the exact same deformations. Beam 3 for which
NLGEOM analysis has been performed with follower forces, shows higher deformations than
beams 2 and 4, and the end where load is applied bulges out into a spherical shape.
Figure 4 shows the differences in deformation characteristics with and without the application of
follower forces for geometrically linear and geometrically nonlinear analyses.
Whether follower force should be applied or not depends on the application. For situations where
the applied force rotates with the rotation of the load application point, follower forces should be
defined for correct representation of the physical situation. In all other situations where the
direction of the force remains constant, follower forces do not need to be considered.
Design Space.
A design for the socket wrench is generated using the OptiStruct solution. The solution shows
major load paths running from the tip, which is loaded with a torque and surface pressure, to the
rigidly held areas farthest on the handle. Translators provided with OptiStruct can be used to
generate IGES surfaces which can be imported into any CAD system for developing a complete
design concept. A finite element model of the concept can also be developed for analysis and, if
necessary, further optimization using size and shape optimization methods can be performed.
Execution
This example is analyzed using the one-file setup with the file, michell.fem. The OptiStruct
batch job is submitted using the command shell script, % optistruct michell.
Results
The optimization converges in 29 iterations. The results are requested in HyperMesh binary
format and written to the file, michell.res. The shape of the solution at the final iteration is
visualized by creating a contour plot of the density results at the 29th iteration in the HyperMesh
contour panel.
Execution
This example is analyzed in the one-file setup with the file, bridge.fem. The OptiStruct batch job
is submitted using the command shell script, % optistruct bridge.
Results
The optimization converges in 24 iterations. The solution is well defined with discrete truss
members connecting the load carrying arch to the load applied points. The results are requested
Design topology for suspension bridge with all loads weighted equally.
PSOLID 1 1 1
The 1 in the 9th field denotes that the component is designated as design material. The example
is run for 20 iterations.
By running the file, beam.HM.comp.cmf, as a command file in HyperMesh, the elements are
grouped into sets according to their final material density values. The set labeled "0.0 - 0.1"
contains all of the elements in which densities range from 0% to 10%. The set labeled "0.1 - 0.2"
contains all of the sets in which densities range from 10% to 20%. The elements in the sets that
have material densities less then 30% are masked so that the solution is easier to visualize. The
material densities of the remaining elements are plotted as shown.
Execution
This example is analyzed in the one-file setup with the file, beam.fem. The OptiStruct batch job is
submitted using the command shell script, % optistruct beam.
Results
The optimization runs for 20 iterations. The results are requested in HyperMesh binary format
and written to the file, beam.res. The shape of the solution at the final iteration is visualized by
creating an assign plot of the density results at the 20th iteration in the HyperMesh contour
panel. By removing components labeled "0.1 - 0.2", "0.2 - 0.3", "0.3 - 0.4", and "0.4 - 0.5" from
the display, a concept of the optimized beam can be visualized.
ITERATION 0
Density plot of air conditioner bracket for combined eigenvalue and compliance objective.
Two thick ribs extend from the engine bolts to the lower air conditioner attachment and one thin
rib extends from the middle of the upper main rib to the upper air conditioner attachment. There
is webbing between the two main ribs and a wide, half-height rib running through the upper half
of the design space.
Many elements did not converge to either 100% dense or 0% dense. These elements can be
forced toward 0% and 100% by increasing the discreteness parameter. This solution is
acceptable because webbings and half-height ribs are manufacturable within a casting process.
Optimized models for linear static and frequency are analyzed separately.
The two main ribs from the combined solution are present in both solutions although they vary in
thickness between the cases. The third rib, which is full height in the combined load case,
appears only in the static solution and is thinner in the combined solution. The webbing between
the two main ribs appears as several discrete crossribs in both separate solutions. The upper rib
is half-height in the combined solution and appears in both the static and eigenvalue solutions.
This rib is in a different position in each analysis. The thicker half-height rib in the combined
solution is a compromise between the two.
Michell-truss Example
For a Michell-truss type solution, a point load is applied at the right hand tip, and nodes
along the circle are fixed. The upper half of the system is modeled because of the
symmetry (Dimension of the FE model = 100*30).
CHECKER,2.
MINMEMB,6.0,1.
MINMEMB,6.0,2.
In this example, an MBB-Beam with a point load in the upper middle point and two roller
point supports at the two lower corners is used. The right half of the beam is modeled
(dimension of design domain = 900*300).
DOPTPRM CHECKER 2
In this example, an arch with three single point loads as three individual load cases on the
lower edge is used. Two roller supports are placed at the two lower corners. The dimension
of the domain equals 500*250.
CHECKER, 2.
MINMEMB,40.0,1.
In this example, a 3D bridge model with a fixed support at both ends is used. Three load
cases: a uniform load p=1 on the entire surface of the bridge and a uniform load p=0.3 on
the lanes in each direction separately. The dimension of the design domain:
length/width/depth = 500/100/100. Target volume fraction = 0.2.
CHECKER, 1.
MINMEMB, 20.0,1.
The reason for running CHECKER=1 in two stages with DISCRETE=1, and then 2, is that
after the convergence with DISCRETE = 1, the solution is far from DISCRETE. The two-
stage solution provides a fair comparison to results with minimum member control since
the latter contains an automatic increase of the DISCRETE value by an increment of 1.
Topography optimization divides the design region into smaller areas, each with its own shape
variable. OptiStruct performs this process using the parameters defined by the user. In this
example, the design space consists of the entire plate minus the areas near the loads and
constraints. The smaller areas can each move upward. The inner portions of these smaller areas
are shown fully deflected in Figure 1.2.
The potential reinforcement pattern can be any combination of these variables deflected at any
height between zero and the user-defined maximum height. OptiStruct creates a reinforcement
pattern of any shape by manipulating the 174 discrete shape variables. The pattern could
resemble an X, an oval, a series of straight beads, or any number of the millions of potential
designs. By setting various parameters, the user can ensure that any design OptiStruct creates is
manufacturable.
Once OptiStruct generates the shape variables, it begins the optimization of the plate. The
objective in this example is simple, the stiffness of the plate under the given torsion load is to be
maximized. OptiStruct performs a series of analysis runs to evaluate the stiffness of the plate,
determines what variable value changes will improve the stiffness of the plate, and applies those
changes to the model. After much iteration, OptiStruct reaches the maximum design at a point
where the stiffness can no longer be improved. The solution is shown in Figure 1.3 where the
colors display the heights of the bead reinforcement patterns. A finite element model based on
that solution is shown in Figure 1.4. The symmetry of the solution is due to the use of an
OptiStruct design symmetry plane feature.
The face of the plate is covered with X-shaped cross-beads that work well in torsion. None of the
beads run completely across the plate in a straight line which would reduce their effectiveness.
Finite element analysis of the plate revealed a well distributed stress pattern and low deflection
at the load point.
The traditional method is to design a reinforcement pattern in the form of raised "beads" across
the surface of the plate, test the stiffness of the plate, and increase the stiffness of the plate until
it meets the design requirements. The efficiency of the plate in terms of cost and weight is
strongly dependent upon how good the reinforcement pattern is, so it is critical to generate a
good one. Two examples of conventional bead reinforcement patterns for the plate in twisting
model are shown in Figure 1.5. These patterns are ones which would commonly be found in
commercial products.
The plate generated by OptiStruct using topography optimization is far stiffer than both of the
conventional plates shown in Figure 1.5. Peak deflection for the topography plate is 0.83mm. For
the conventional plates, the peak deflections are 1.27mm for the one on the left and 6.47mm for
the one on the right. The plate developed with OptiStruct is 35% stiffer than a good conventional
design and far better than a poor one. The poor design, while following the conventional wisdom
of using an X-shaped reinforcement pattern, uses beads that run completely across the plate in a
straight line which are susceptible to kinking when loaded. Such design mistakes are caught and
corrected by OptiStruct during the optimization process allowing it to yield a superior design.
For the input file sample, see
<install_directory>/demos/hwsolvers/optistruct/twistplate.fem.
The design domain includes the entire box with the exception of the filling hole on the top (shown
in red). All of the elements in the design domain are placed in the same component and
reference the same material property. The normal vectors for all of the elements in the design
domain are pointing outward. The topology variables are set up with the following DTPG card:
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 1 PSHELL 1
DESOBJ(MIN) 1
DESSUB 101
The objective is to minimize compliance for the pressure load case, which is the same as
minimizing the strain energy of the entire model. The displacement of the center point of each of
the five loaded surfaces was constrained to be less than a given value.
OptiStruct topography solution for pressure vessel, front and rear views.
All of the optimization constraints are met for the model. The red areas represent the bead
reinforcements that OptiStruct created to increase the stiffness of the model. Circular or oval
reinforcement beads are generated for the large side panels and the bottom panel of the box.
Circular and oval beads are very effective in stiffening the panels against a distributed or central
load. This is due to the fact that bending in the central areas of the panels is occurring in two
directions, both vertically and horizontally. Straight beads provide stiffness for bending in one
direction, but would be vulnerable to kinking in the other. Round beads provide stiffness without
kinking. Bulbous beads are created at the eight corners of the model, anchoring the sides of the
box together and allowing each side of the box to gain support from the adjacent sides.
By changing the draw height on the DTPG card from 4.5 to –4.5, the direction of the beads can
be changed from outward to inward.
The bead pattern for the inward bead model is almost a complete reflection of the one that
OptiStruct created for the outward bead model. The areas that are not pushed out in the outward
bead model are pushed in for the inward bead model, generating the same basic structure.
For the input file sample, see
<install_directory>/demos/hwsolvers/optistruct/pressurebox.fem.
The elements shown in green are included in the design space. The nodes in the top layer can
move upward and the nodes on the bottom layer can move downward. Symmetry is used to
force the bead pattern on the top to match the one on the bottom. The DTPG card is generated
as follows:
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 2 DVGRID 1
The plane of symmetry is normal to the vertical axis and is positioned running through the center
plane of the model.
The optimization objective function is simply defined as minimizing the sum of the weighted
compliance of all three load cases.
The solution for the model is shown in the figures below. The first figure shows the symmetry of
the solution about the vertical axis.
The solution shows the importance of adding vertical bending stiffness in the area around the
shock absorber attachment point. OptiStruct creates a large bead running from the upper frame
attachment point, past the shock attachment, and up to the spindle attachment (second figure).
This bead supports most of the bending load. In addition, there is vertical bending which runs in
the perpendicular direction. OptiStruct creates a bead running from the shock attachment point
to the lower frame attachment (second figure). This second bead is not as pronounced because
there is less bending in that direction compared to the primary direction.
For the input file sample, see
<install_directory>/demos/hwsolvers/optistruct/stampedarm.fem.
The red areas are excluded from the design domain. The blue area is open for OptiStruct to add a
bead reinforcement pattern. The bead is drawn upward with respect to the plate orientation. The
DTPG card used is as follows. Four different runs were made using different values for the draw
height. The first run was made with a draw height of 20mm, the second with 40mm, the third
with 60mm, and the fourth with 80mm.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 1 PSHELL 5
The optimization is set up to maximize the frequencies of the first six modes (minimizing the sum
of the weighted inverse eigenvalues) and to ensure that the first three modes were above certain
design constraints. This is accomplished by placing the following cards in the subcase definition:
DESOBJ(MIN) 1
DESSUB 101
Setting constraints on the first three modes results in separations between the frequency values
of the modes and prevents OptiStruct from falling into local minimums when optimizing the
modes. This approach ensures that a minimum performance criterion is satisfied. Note that for
the 40mm, 60mm, and 80mm draw height runs, the constrained frequencies are higher than
those shown above.
The solutions generated for the plate runs are shown in the following figures:
The reinforcement patterns have a similar shape, but runs with a higher maximum draw height
use more levels of draws throughout the plate. All of the solutions made good engineering sense,
connecting the weak areas of the plate with beads running primarily across the short span of the
plate. These beads were fluidly connected together across the long span of the plate allowing the
beads to reinforce each other.
For the input file sample, see
<install_directory>/demos/hwsolvers/optistruct/eigenplate.fem.
To use topography optimization to optimize the surface of a solid model, the areas to be
optimized must first be covered with shell elements (see gray elements in figure below).
Solid control arm model partially covered in shell elements (gray areas).
The side surfaces of the structural members are optimized in this example.
OptiStruct solution for surface optimized control arm, areas to be thinned are red.
OptiStruct undeformed model for surface optimized control arm, top view.
Loads and constraints for a plate to be optimized for both shape and size.
The design shape variables are generated externally to OptiStruct. Size variation is defined by
creating a single design variable that controls the thickness of the plate. Shape variation is
defined by creating a single design variable that controls the height of the plate relative to the
constrained nodes.
Models showing the size and shape variations at maximum thickness and height.
The two design variables are included in the finite element deck submitted to OptiStruct. Each
variable has a single DESVAR card and a DVGRID for each node in the model affected by that
variable. DTPG cards added to the deck reference the ID number of the DESVAR cards as follows:
DTPG 3 DVGRID 1
DTPG 4 DVGRID 2
3.0 60.0 NO
DVGRID* 1 4 0 1.0
DVGRID* 1 5 0 1.0
Entries in fields 15 through 18 are not required because the draw vectors are already defined for
each node on the DVGRID cards. A buffer zone is requested for the shape variable that ensures a
smooth transition between the constraints and the rest of the structure. Lengthwise symmetry is
enforced (the left side of the plate mirrors the right side).
The compliance of the part was minimized with the mass of the part constrained to be below a
given value.
The solution is shown in the figures below.
OptiStruct solution for shape and size optimization of plate in bending, top view.
OptiStruct creates a thick rib around the perimeter of the plate, providing a solid foundation for
the attachment of the remaining design features. The highest degree of bending occurs at the
center of the part, with the largest bending component running in the longitudinal direction. The
plate is stiffened in this direction when OptiStruct creates a central rib and a W-shaped
reinforcement pattern by raising the thick areas of the plate and lowering the thin areas of the
plate. This shape employs the entire structure, not just the parts with ribs.
Figure 1.1: Loads and constraints for the spot welded tube.
The initial cross-section of the tube is arbitrarily chosen to be roughly circular. Two similar
topography variables are defined to allow the shape of the tube to change. The two variables
assign the blue and green pieces to be in the design domain while the flanges remain in their
original shape. The topography variables are in the bulk data section using the DTPG card as
shown below.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 1 PSHELL 1
The draw height of 5.0, combined with the upper and lower bounds of 2.0 and –2.0, allow a bead
height of 10.0 model units in both directions (inward and outward) and a bead growth direction
normal to the surface of the tube. A linear type pattern grouping is applied for this problem. This
means that the beads formed during the topology optimization will be constant along a line
parallel to the direction defined -- the central axis (X-axis) of the tube, in this case. By setting
this pattern grouping option, the cross-section of the tube will be allowed to change, but will
remain constant through the length of the tube.
The solution is not completely smooth, but the basic shape of the tube is clear. The lower half of
the tube has been lowered to increase the bending stiffness of the section while the upper half of
the tube runs directly from one flange to the other to support the shear force generated by the
twist load. The cross-sections of the model before and after optimization are shown in Figure 1.3.
The objective for the model is to increase the sum of the frequencies of the first four normal
modes of the disc. This is achieved by using the WFREQ response type. The WFREQ response is
the sum of the inverse eigenvalues of the chosen frequencies. This is done in order to assign
higher weight to the earlier modes than the latter ones.
The first manufacturing method to be considered is turning the disc on a lathe. This restricts the
bead reinforcements to being circular. A circular pattern grouping type is define. The DTPG card
for this card is shown below.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 1 PSHELL 1
The second manufacturing method to be considered is stamping with radial beads only. The
topography set up for this approach is almost identical to that of the one for the circular approach
with the only change being the radial pattern grouping type. The DTPG card for this configuration
is shown below.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 1 PSHELL 1
The solution for the radial reinforcement pattern is shown in Figure 2.3. OptiStruct generated a
series of eight evenly spaced wide radial beads. Note that the beads are not drawn to their full
height. This demonstrates that full height beads are not necessarily optimal. This reinforcement
pattern also doubled the frequencies of the first three modes and increased the fourth and fifth
frequencies by more than 70%, but was not quite as efficient as the circular pattern.
The third manufacturing method to be considered was stamping with no constraints on the bead
shapes. No pattern grouping was defined for this approach. The DTPG card for this set up is
shown below.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 1 PSHELL 1
The solution for this setting is shown in Figure 2.4. OptiStruct generated a reinforcement pattern
that appears to be a combination of the circular and radial reinforcement patterns. The inside of
the disc has a series of radial beads while the outside has a roughly circular reinforcement. This
reinforcement pattern had the best increases in the modal frequencies for all of the patterns
compared.
Since the default reinforcement pattern is almost symmetric, a cyclically symmetric pattern can
be used to clean up the solution. Choosing four wedges for cyclical symmetry appears to be very
close to the pattern created with the default settings. The cyclical symmetry pattern grouping
type is applied. The DTPG card for the cyclical symmetry method is shown below.
DTPG 1 PSHELL 1
For cyclical symmetry, both the first vector (fields 12, 13, and 14) and second vector are defined
in the same plane as the disc. The anchor node is located at the center of the disc. This cyclical
pattern is set up to have four symmetric wedges with each wedge being symmetric about itself.
The solution for the cyclical pattern grouping option is shown in Figure 2.6. OptiStruct generated
a reinforcement pattern very similar to the default pattern grouping pattern but with fourfold
symmetry. The modal frequencies were slightly less than those for the default pattern.
The frequency results (in Hz) for the first six normal modes for the baseline model and four
reinforcement patterns are shown in the following table:
The entire box is to be used as the design domain with the exception of the filling hole on the top
shown in red. All of the elements in the design domain are placed in the same component and
reference the same material property. The normal vectors for all of the elements in the design
domain are pointing outward. The topology variables are set up with the following DTPG card:
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 1 PSHELL 1
Three orthogonal planes of symmetry are defined (see Figure 3.2). The anchor node is placed at
the center of the box. The first and second vectors are defined parallel to the X and Y axes. The
first vector is defined pointing away from the filler cap. This ensures that the automatically
generated variables will cover the entire surface of the box. If the first vector was pointing in the
other direction, the symmetry method would reflect the lack of variables in the area of the filler
cap across all three planes of symmetry.
The objective is to minimize compliance for the pressure load case, which is the same as
minimizing the strain energy of the entire model. The displacement of the center point of each of
the five loaded surfaces was constrained to be less than a given value.
OptiStruct generated the following solution for the pressure box shown in Figure 3.3.
The OptiStruct solution met all of the optimization constraints and yielded a good design. The
areas shown in red (In Figure 3.3) are the bead reinforcements that OptiStruct created to
increase the stiffness of the model. The solution is unconventional, but makes a good deal of
engineering sense. For the large side panels and the top and bottom panels of the box,
OptiStruct has generated large, rounded, centrally located reinforcement beads. These types of
beads are very effective in stiffening the panels against a distributed or central load. This is due
to the fact that bending in the central areas of the panels is occurring in two directions, both
vertically and horizontally. The rounded beads create stiffness in both directions and are weak in
neither. At the eight corners of the model, OptiStruct created beads that anchor the sides of the
box together allowing each side of the box to gain support from the neighboring sides.
A finite element model was created from this reinforcement pattern (see Figure 3.4).
Figure 3.4: Finite element model of the OptiStruct solution for the pressure box.
The OptiStruct model is superior in stiffness to both of the two conventional models. The
maximum deflection of the OptiStruct model was 30% less than the lightly reinforced
conventional model (the one on the left), and 46% less than the heavily reinforced model (the
one on the right). The lightly reinforced model was stiffer than the heavily reinforced model,
which goes against the assumption that more reinforcements result in increased stiffness. With
bead type reinforcements that assumption is not always true, which demonstrates the
effectiveness of topography optimization. OptiStruct delivers an optimized first design,
eliminating the need to do a series of re-designs where the second, third, fourth, etc., model
does not always result in an improvement.
Manufacturing constraints can be accounted for in the pressure box model using other pattern
grouping options. The set up is done easily through the HyperMesh interface. In order to
manufacture the pressure box using a two piece die mold, bead reinforcements that run laterally
would need to be eliminated or else they would cause a die lock condition. Topography
optimization can be used to generate reinforcement beads on the sides of the box that run
vertically only. This is done by using separate topography variables for the side walls, front and
back panels, and the top and bottom panels. Topography variables with planar symmetry with
one plane symmetry grouping type are defined for the side panels and the front/back panels,
respectively. As in the earlier case, a three plane symmetry topography variable is assigned for
the top/bottom panels. The two cards are shown below and differ only in the direction of the first
vector.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 1 PSHELL 1
DTPG 1 PSHELL 1
The planes for the planar pattern grouping run vertically and perpendicular to the Y-axis, which
causes OptiStruct to generate vertical beads. Because the planes run through both sides of the
box, there will be symmetry between opposing sides. Also, symmetry on either side of the anchor
node in the direction of the first vector (Y-axis) is forced with the one plane symmetry option
(pattern grouping option 13). See Figure 3.6.
Figure 3.6: Planar pattern grouping planes for sides of pressure box.
With planar symmetry enforced for the sides of the pressure box and three-plane symmetry
enforced on the top and bottom of the box, OptiStruct generated the solution shown in Figure
3.7. Even without the presence of lateral beads on the sides of the box, the OptiStruct solution
shown below had a maximum deflection 6% less than the lightly reinforced conventional model.
Figure 4.1: Loads and constraints for the stamped hat section.
It is preferable to have the size of the reinforcements able to run deeper than the height of a
single element. To ensure that this will not cause a problem with the element mesh, three shape
variables are created using HyperMesh and are added to the deck. The shape variables for the
face and top side of the hat are shown in Figure 4.2.
Note in Figure 4.2 that the first three rows of elements adjacent to the elements being fully
deflected are a part of the user-defined shape variable for that side. Also, the draw depth is equal
to one and a half times the average element size.
It is desired to create this hat section using a stamping process which means that reinforcing
features on the sides of the hat must be constant (from top to bottom), or else a die lock
condition will occur. Pattern grouping can be used to create variables that ensure
manufacturability. For the three variables created for the hat section optimization, planar pattern
grouping was selected with the planes running perpendicular to the length of the section (X-axis).
The DTPG card and associated DESVAR card for one of the variables are shown below.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 4 DVGRID 1
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Additionally, a plane of symmetry was used to force both halves of the section to be the same.
OptiStruct generates variables which allow for great flexibility in the reinforcement possibilities,
but which prevent a die lock condition as shown in Figure 4.4. Note that the area where the load
is applied is left out of the face variable.
The objective is to minimize the compliance for the applied load. OptiStruct generated the shape
shown in Figure 4.5.
The solution generated by OptiStruct is manufacturable using a stamping process. Also, the
solution is very well behaved and needs little refinement to turn it into a production-ready
design.
The optimized hat section increases the stiffness of the part by more than eightfold from the
initial condition with no beads. The eight ‘square’ beads for the hat section, especially the four at
the ends of the beam, are the key to bolstering the beam against shear collapse. Those beads
also serve to prevent the flanges from folding under the bending load. OptiStruct has generated a
strong design that supports both torsion and bending with restricted reinforcement possibilities.
The shape and placement of the reinforcements are optimized resulting in a very efficient
solution.
For the input file sample, see <install_directory>/demos/hwsolvers/optistruct/hat.fem.
The first step is to generate shape variables for the part. HyperMesh is an excellent tool for doing
this. For this model, four basic shape variables were created which control the sizes of the two
legs and rear beam of the control arm. These are shown below in Figures 5.2 and 5.3. Care was
taken to ensure that the final design would have no internal cavities which would prevent it from
being manufacturable.
Figure 5.3: Shape variables controlling the thickness and height of the control arm rear beam.
Figure 5.5 shows the variables created from the planar pattern grouping option. Note that Figure
5.5 only shows the deflections of the centers of the variables. Because of the way that
topography optimization works, if all three variables were fully deflected, the legs of the control
arm would be uniformly at the minimum height.
Planar variable pattern grouping was also used for the height and thickness of the rear beam. For
both variables cards, TYP = 13 was used to ensure a beam shape that was symmetric about the
center. The DTPG card for variable #4 (beam height) and the associated DESVAR card are shown
below.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 8 DVGRID 4
+ 15.0 60.0 NO
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Since shape variables are being assigned to topography, the height and the bounds in the DTPG
card do not need to be defined, as they are controlled by the shape variable. Even if these values
are defined, OptiStruct will ignore them.
The objective for the optimization was to minimize the mean compliance of the control arm for all
three load cases combined. The mass of the control arm was constrained to be less than 10%
below its initial mass. OptiStruct generated the following solution (see Figure 5.6). The red shows
areas where material was added, the blue shows areas where material was removed.
OptiStruct optimizes the shape of the legs by shifting mass from the centers to the top and
bottom, creating C-shaped sections. These sections add vertical bending stiffness while leaving
the shear and axial stiffness intact. Additionally, the height of the legs was increased along their
length to further increase the bending stiffness of the sections. OptiStruct greatly reduced the
size of the rear beam in both height and thickness, indicating that it was oversized or perhaps
even unnecessary. Overall, the maximum deflection for all three load cases was reduced by 5%
while the mass was reduced by 10%. Figure 5.7 shows a refined finite element model of the
OptiStruct solution.
Figure 6.1: Loads and constraints for the solid block model.
Two shape variables are generated using HyperMesh to control the height and width of the block.
These are shown in Figures 6.2 and 6.3.
It is preferable to manufacture the resulting part using a casting process. This can be
accomplished by using a linear pattern grouping in the casting draw direction and a planar
pattern grouping perpendicular to the draw. This will ensure that there are no cavities that would
create a die lock situation.
Thinking ahead, it is predictable that the cross-section of the solution will be roughly an I-shaped
section with the web running vertically. This prediction establishes the draw direction as being
horizontal, which corresponds to variable #1 (block width), thus variable #1 will be split using
linear pattern grouping and variable #2 will be split using planar pattern grouping. The DTPG
cards and associated DESVAR cards are shown below:
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 3 DVGRID 1
+ 2.0 60.0 NO
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
DTPG 4 DVGRID 2
+ 20.0 60.0 NO
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
The linear variable dispersion pattern for variable #1 allows OptiStruct to control the thickness of
the block at numerous points across its side giving the solution a great deal of flexibility. The
planar variable dispersion pattern for variable #2 allows OptiStruct to control the height of the
cross-sections along the length of the block. The objective was to minimize the mean compliance
of the block under the given load. The mass was constrained to be below one fourth of the initial
mass of the block. OptiStruct generated the solution, shown in Figure 6.4.
The cross-section of the block is roughly I shaped, concentrating the material at the top and
bottom of the end and center areas where the bending moment is the greatest. The design is flat
and tall in areas where shear is dominant. The solution is manufacturable by use of a casting
process since there are no cavities or die lock conditions. The efficiency of the solution can be
seen by looking at the stress contours. The nearly uniform stress distribution, shown in Figure
6.5, indicates that almost every part of the structure is being used to its fullest potential.
Figure 6.6: OptiStruct solution for solid block with 2.5 times smaller cross-section.
The basic shape of the block is the same in the reduced dimension model, but has more
pronounced features. The I shaped cross-sections in the center and at the ends have wider
flanges, and the shear carrying areas in between are thinner. This makes sense considering the
smaller dimensions increase the need for bending stiffness more than the need for shear
stiffness.
For the input file sample, see <install_directory>/demos/hwsolvers/optistruct/block.fem
and <install_directory>/demos/hwsolvers/optistruct/blocklong.fem.
Structural Model
In this example, the thickness for each element is varied. Therefore, five variables are defined in
the analysis deck. The initial thickness of all elements is 2mm. All thickness variables have a
lower bound of 1mm and an upper bound of 3mm.
The optimization problem, objective and constraint functions, is defined in the same manner as
for the other types of structural optimization. OptiStruct then goes through an iteration process
to solve the optimization problem. OptiStruct converged after 3 iterations. The table below shows
the optimization results of the cantilever beam, which met the constraints and minimizes the
mass (decrease of 6.6%).
Structural Model
In this example, the cross-sectional area for each element is varied. Therefore, 10 variables are
defined in the analysis deck. The initial cross-sectional area of all the elements is 20, with a lower
bound of 0.2 and an upper bound of 200.
The optimization problem, objective and constraint functions, is defined in the same manner as
for the other types of structural optimization. OptiStruct then goes through an iteration process
to solve the optimization problem. OptiStruct converged after thirteen iterations. The table below
shows the optimization results of the ten bar truss structure, which met the constraints and
minimizes the mass (decrease of 38.7%).
Bar1 Bar2 Bar3 Bar4 Bar5 Bar6 Bar7 Bar8 Bar9 Bar10 Volume
Initia 20 20 20 20 20 20 20 20 20 20 83929
l
Final 30.4 38.4 24.5 15.4 0.2 0.2 8.8 21.5 20.1 0.4 51448
DEQATN,111,A(D)=PI(1)*D**2/4
$MOMENT OF INERTIA
DEQATN,122,I(D)=PI(1)*D**4/64
$TORSIONAL CONSTANT
DEQATN,133,J(D)=PI(1)*D**4/32
Using these equations, DVPREL2 statements are used to assign each design variable to the
respective PBAR property. The statements that assign the diameter of the first bar element to the
cross-sectional area of that element look like:
DESVAR,1,Diam1,10,1,20,0.5
DVPREL2,11,PBAR,1,4,,,111
+,DESVAR,1
The optimization problem to be solved is the minimization of the tip displacement with a volume
constraint of 4000. Convergence was achieved after four iterations.
The result is given in the following table.
Constraint: The surface of the cube should be between 2.0 and 3.0
Design Variables: A, B, C
The volume and surface are defined as equations using DRESP2 and DEQATN:
$
$ VOLUME
$
DEQATN 1 VOL(W,L,H)=W*L*H
$
$ SURFACE
$
DEQATN 2 AREA(W,L,H)=2.0*(W*H+L*H+W*L)
$
DRESP2 1 VOLUME 1
DESVAR 1 2 3
DRESP2 2 SURFACE 2
DESVAR 1 2 3
$
DESVAR 1 W 1.1 0.1 10.0
DESVAR 2 L 0.9 0.1 10.0
DESVAR 3 H 2.0 0.1 10.0
$
Then, in the optimization problem, the objective and constraint are global responses (for
example, DESOBJ and DESGLB are used outside of a SUBCASE).
To trick OptiStruct into solving this problem, a dummy finite element model must be provided.
Here, a single shell element with some load is used.
As expected, the solution yields a cube with even sides of about 0.707, a surface of 3.0, and a
volume of 3.53.
The input file can be found in <install_directory>/demos/hwsolvers/optistruct/box.fem.
Structural Model.
The shape of the beam is defined using a linear combination of the two basis shapes below. The
linear factors are the design variables in the optimization problem.
Basis Shape 1.
The basis shapes are defined using shape basis vectors. These can be generated using AutoDV,
which is part of HyperMesh. The output file of AutoDV contains the definition of the DESVAR and
DVGRID cards. This file can then be included in the OptiStruct input file, bulk data section, via
the INCLUDE statement.
The optimization problem, objective and constraint functions, is defined in the same manner as
the other types of structural optimization. OptiStruct then goes through an iteration process to
solve the optimization problem. The figures below show the optimization result of the cantilever
beam.
Final Shape.
The design domain is subdivided into two design elements; the web and the flange. Six design
variables are defined using the design elements and vectors (Fig. 2). For shape optimization, the
shape of the beam is defined using the nodal positions of the original shape and a linear
combination of the six shape perturbations associated with the design variables.
The linear factors are the design variables in the optimization problem. The shape
of the beam appears as:
.
Figure 3 shows the shape of the beam perturbed by the first design variable, which is a linear
perturbation. Figure 4 shows the quadratic perturbation caused by design variable 4.
The perturbation vectors need to be provided in the format of the DVGRID cards
using AutoDV (part of HyperMesh). These cards can be generated automatically. The output of
AutoDV also includes the design variable definition DESVAR. The output file Beam_shape.dat can
be incorporated into the bulk data section of the OptiStruct input deck via an include statement.
The definition of the optimization problem is included in the case control section of the input
deck. Figure 5 shows the section of the OptiStruct input file that includes the definition of the
optimization problem and the inclusion of the AutoDV output.
All optimization constraints are met for the model. The final shape is shown in Figure 5.
Cantilever beam. input data
Rotating bar
SPC1 10 123456 1
:
:
ENDDATA
The input file is just like one for an ordinary min-max problem. The maximum normal stress of
the flexible body of the subcase 1 is to be minimized by using DOBJREF and MINMAX.
Notable points include:
Because the stress in multi-body dynamics systems is a time variant quantity, the
minimization of stress in multi-body dynamics analysis subcases should be a min-max
problem.
The SPC1 card fixes only 6 dof of node 1 of the flexible body in order to remove 6 rigid
body motions. If you fix more than 6 dof of the flexible body, the additional fixed dof
become constraints of the flexible body.
Stress output request is placed above the first subcase. If you place output request inside
subcase, your output request will be ignored.
The input file can be found in
<install_directory>/demos/hwsolvers/optistruct/rotating_bar_design.fem.
Slider crank
Rotating shell
After analyzing the initial model, the time history of stress using HyperView can be seen.
:
:
Now you can find out the behavior of the structure as in the analysis result image.
2. According to the results post-processed by HyperView, the maximum stresses are developed
at around 0.3 seconds. Increase the number of time steps around 0.3 seconds. That is, divide
the time period of 0.28 seconds – 0.34 seconds into 200 steps. Increasing the number of time
steps in this period provides the optimizer with more information.
The element that has the maximum stress and the corresponding time can be changed as the
design changes. Thus, Step 2 does not always work. If the time when maximum stress is
developed and corresponding element are expected to change dramatically as the design
changes, it is best to consider the changed peak time and corresponding element as much as
possible.
3. Replace the previous single MBSIM card with multiple MBSIM cards as the following.
MBSIM 1 TRANS END 0.28 NSTEPS 50
+ VSTIFF
MBSIM 2 TRANS END 0.34 NSTEPS 200
+ VSTIFF
MBSIM 3 TRANS END 0.70 NSTEPS 50
+ VSTIFF
MBSIM 4 TRANS END 1.0 NSTEPS 100
+ VSTIFF
MBSEQ 10 1 2 3
:
:
MBSIM 1 TRANS END 0.28 NSTEPS 50
+ VSTIFF
MBSIM 2 TRANS END 0.34 NSTEPS 200
+ VSTIFF
MBSIM 3 TRANS END 0.70 NSTEPS 50
+ VSTIFF
MBSEQ 10 1 2 3
Using the above steps, the design process convergence can be enhanced.
The input file can be found in
<install_directory>/demos/hwsolvers/optistruct/rotating_shell_design.fem.
In order for the system level responses to be made available to the optimization process,
MBREQE or MBREQM entries must be defined so that the MBD solver can generate the
appropriate output.
$------|--MID--|--GID--|
MARKER 55 9929617
$------|--RSID-|--RID--|--ITEM-|--MID--|
MBREQM 99 999 VEL 55
Changing the length of a body is a common method of achieving desired system level responses.
This change in shape can be defined using DVGRIDs. In this case, the DVGRIDs need to be
defined carefully, as it is necessary to ensure that the joints remain in a valid configuration while
changing the length of the bodies.
In the figure above, Nodes A and B must remain coincident after applying the shape
perturbation.
The correct configuration after applying shape perturbation vector should be:
Correct way of defining shape change that involves change of joint locations