Hwsolvers12.0.220 ReleaseNotes
Hwsolvers12.0.220 ReleaseNotes
Hwsolvers12.0.220 ReleaseNotes
MotionSolve
Improved workflow for co-simulation using Simulink CoderTM
The workflow for co-simulation between MotionSolve and Simulink via the Simulink CoderTM has been
revised for ease of use. In the current version, you can make use of a script to compile your Simulink
CoderTM model into a dynamic linked library.
The new workflow is illustrated below:
This script compiles and links the code generated by Simulink CoderTM automatically, making it ready to
use with MotionSolve for co-simulation. Please refer to the MV-7005 tutorial for more information.
The first change improves the robustness of models containing PTdSF elements. It was observed
that in the corrector iterations, the surface parameters U and/or V could go out of the range
defined for the surface and the simulation would stop. In the enhanced implementation, if the U or
V parameters begin to go out of range, MotionSolve prints out a warning and continues with the
simulation. The U or V parameters are held fixed at the start or end limits until they come back into
range. If the model is well defined, the U or V parameters will come back into range.
The second change increases the generality of PTdSF. PTdSFSUB the user subroutine definition
for the PTdSF model element has been enhanced to accept SYSFNC and SYSARY calls. This means
that the force in a PTdSF element may now be a function of other model states. This is particularly
useful for defining complex force models inside the PTdSFSUB. For instance, the coefficient of
friction may be calculated using a differential equation. This state dependent value may be
subsequently accessed in the PTdSFSUB.
RIGID elements in the flex h3d are now written out in the animation h3d for better visualization
Shell thickness information in the flex h3d is transferred to the animation h3d for better
visualization
Improved error messages if non-matching parenthesis are detected in any function expression
More understandable messages when redundant constraints are detected in the model
More understandable messages when a U or V parameter goes out of range in a PTCV, PTdCV,
PTSF, PTdSF etc. modeling element
More understandable messages when a negative stiffness or damping value is specified for a spring
damper or bushing modeling element
Modeling check for spring damper, bushing, beam and field elements
Previously, if you specified I and J markers belonging to the same body for the above model elements,
MotionSolve would fail during the analysis with little information about the cause. This made debugging
the model difficult. Within this release, a mechanism that checks this condition has been implemented. If
the I and J markers for the above elements belong to the same body, a relevant error message is printed
out and the simulation is stopped.
An equal but opposite force is applied on two bodies that these entities are defined between
The force applied depends on the displacement and velocity of one body relative to the other
Outputs specific to each entity are supported, for example tire lateral slip angle for the tire entity
The parameters (for example stiffness), used to compute the force are stored in a property file that
is independent of the MotionView MDL file and the MotionSolve input deck. The property files are:
o Text files and may be viewed and edited using any common text editor
o Compatible with Adams/CarTM property files from Adams/CarTM may be reused with the
Auto Entities in MotionView
o Read by MotionSolve just prior to analysis.
Function
Features
Elastically mount
suspension arms, links,
sub-frames, engines,
gearboxes, etc. to
vehicle chassis
Bump-stop
(autoBumpstop)
Uses
Suspension
kinematic and
compliance analysis
Vehicle Dynamics
Vehicle Dynamics,
Durability
Damper
(autoDamper)
Rebound-stop
Limit suspension
rebound travel
(autoReboundstop)
Spring
(autoSpring)
Limit forces in
suspension components
and isolate chassis by
absorbing energy from
vertical road inputs
Tire CTI
(autoTireCTI)
Tire CD
(autoTireCD)
Tire TNO
(autoTireTNO)
data
Linear Damping
Vehicle Dynamics
Ride Comfort
Durability
Vehicle Dynamics
Durability
Suspension
Kinematics and
Compliance
Vehicle Dynamics
Ride Comfort
Durability
Vehicle Dynamics,
Ride Comfort
Durability
MF-TYRE
SWIFT-TYRE (rigid-ring)
2D and 3D Road Models,
including CRG.
MF-Tyre:
Vehicle Dynamics
SWIFT-Tyre: Vehicle
Dynamics & Ride
Comfort
FTIRE:
Ride Comfort &
Durability
Pacejka 2000:
Vehicle Dynamics
ADAMS Compatibility
Additional support for model units
More options for choosing the model units are now available. Some of these include MegaNewton,
Poundal, US Ton, Yard, Nanosecond, Day etc.
Resolved Issues
Contact force visualization
While visualizing contact force vectors in HyperView, the force vectors are now correctly drawn to
originate from the point(s) of contact instead of the body CG.
Deactivating joints
Previously, if you deactivated a joint between simulations, any motions associated with that joint would
still be active which could lead to erroneous results. Now, in such a situation, any motions associated with
a deactivated joint are deactivated automatically and a warning message is printed out to alert the user of
the same.
Documentation
A number of errors in the documentation have been resolved for both model and command elements.
Additional Notes
HyperWorks Solver
12.0.212 Release Notes
Introduction
HyperWorks HW12.0.212-HWSolver contains enhancements and bug fixes to HyperXtrude Solver and
OptiStruct.
The details are documented below.
OptiStruct
The RANDOM ID is now written to the XYPUNCH header for PSDF output
Since multiple RANDOM data is allowed, the RANDOM data ID is now written to the XYPUNCH header for
PSDF output. This is so the output related to each RANDOM data can be identified. Note that if there is
only a single RANDOM datum, the ID will not be written to the header.
The DMIGMOD data can now be used to translate and rotate CMS Super Elements
In previous versions, the DMIGMOD data could only be used to translate CMS Super Elements. No
rotation was allowed. This restriction has been removed. The translation and rotation is accomplished by
specifying three grids on the CMS Super Element and three grids in the residual structure that will match
up with the moved CMS Super Element.
DMIGMOD: The DMIGMOD data now has specifications for three grids on the Super Elements
and the three grids that they will match up with in the residual structure when the Super
Element is moved into position.
MPCFORCE for RESVEC modes is now written to the .op2 and .pch files
MPCFORCE for normal modes when USET,U6 data is present is now written to the .op2 and .pch files.
addition, for modal frequency response and modal transient analysis, the MPCFORCE for RESVEC are
written to the .op2 and .pch files when the MPCFORCE(MODAL) command is used.
In
GPFORCE for RESVEC modes is now written to the .op2 and .pch files
GPFORCE for normal modes when USET,U6 data is present is now written to the .op2 and .pch files.
addition, for modal frequency response and modal transient analysis, the GPFORCE for RESVEC are
written to the .op2 and .pch files when the GPFORCE(MODAL) command is used. (212)
In
Modal STRAIN results are available for frequency response and transient
The STRAIN(MODAL) request will generate STRAIN output for the structural and RESVEC modes to the
.op2 and .pch files for frequency response and transient analysis.
PARAM,RBE3FREE can be used to disable ERROR 772 due to free spiders on RBE3 elements
Free spiders on RBE3 elements will singular in some cases. When AMSES or AMLS is used, this modeling
error will lead to invalid results as the free spider will be AUTOSPCed. These results will be invalid
because the RBE3 is now constrained due to the AUTOSPC. For Lanczos and static analysis, the
singularity will cause errors 155 or 153 to be issued. PARAM, RBE3FREE,YES will change ERROR 772 into
WARNING 824. If the free spider is singular, the run will fail with ERROR 155 or 153 for static analysis
and Lanczos eigenvalue solutions. For AMSES and AMLS solutions, the AUTOSPC output should be
carefully checked to make sure the RBE3 has not been constrained.
PARAM,RBE2FREE can be used to disable ERROR 725 due to rigid body rotation of RBE2
elements
RBE2 elements with rigid body rotations are singular. When AMSES or AMLS is used, the free rotational
DOF will be AUTOSPCed. For Lanczos and static analysis, the singularity will cause errors 155 or 153 to
be issued. PARAM, RBE2FREE,YES will change ERROR 725 into WARNING 825. In order to solve the
problem with Lanczos or static analysis, the user must use AMSES and note the RBE2 DOF that are
AUTOSPCed, and then manually SPC these DOF in the subsequent static analysis or Lanczos eigenvalue
analysis.
10
Resolved Issues
USET,U6 RESVEC are now calculated when Lanczos is used for normal modes analysis.
MFLUID results could be wrong if both sides of the panel were wetted surfaces. This is has
been resolved.
Unreferenced MFLUID data could lead to ERROR 2894. This issue has been resolved.
Negative SPC ID were displayed in the .out file. This issue has been resolved.
11
HyperWorks Solver
12.0.210 Release Notes
With the 12.0.210 release, a renaming of our structural solver products has taken place in order to
accommodate the growing trend of applying optimization technologies more intensively in the product
development process, instead of just single run analysis. Moving forward, our structural and thermal
analysis and optimization product for linear and non-linear events will be called OptiStruct; this includes
analysis disciplines such as linear and non-linear statics, noise and vibrations and thermal analysis. Our
product for highly-nonlinear, typically transient events such as crash, impact, forming or blast analysis will
be RADIOSS. Both codes are developed with a strong focus on the above mentioned disciplines, and will
continue to leverage technology mutually. OptiStruct will continue to leverage the widely used bulk data
input format, and RADIOSS will continue to use the well suited block format.
Moving forward, OptiStruct will also serve as the general framework for optimization involving multiple
physics including CFD and FSI with AcuSolve and multi-body dynamics with MotionSolve, and will also
offer API to integrate additional solvers and solutions within a general multi-model framework.
12
Introduction
Altair OptiStruct is the state-of-the-art finite element solver for the solution of a wide variety of
engineering problems such as linear and non-linear statics, frequency response and transient analysis,
complex eigenvalue analysis, and mechanical systems. This robust, multidisciplinary solver enables
designers to optimize performance related to durability, NVH, crash, safety, manufacturability, and fluidstructure interaction, in order to bring innovative products to market faster.
The 12.0.210 release of OptiStruct is yet another functionality-rich release. This release features the
support of frequency dependent material properties. Additionally, the ADJUST capability has been added
to the contact formulation to remove initial penetrations. Some key highlights of the 12.0.210 release are
listed below.
Bolt Pretension Forces are Now Written to the ASCII .pret File
Use the command PRETBOLT(OPTI)=YES to create the ASCII .pret file that contains the pretension force
for each bolt.
13
Gasket Tension Stiffness can be Specified as a Stabilization Coefficient or the Direct Tension
Modulus
The MGASK data EPLTYPE is used to specify if the EPL data is the tension stabilization coefficient (default),
or the direct tensile modulus.
The Thermal Field From the Last Time Step of a Transient Heat Transfer Analysis can be Used
for the Thermal Load
The thermal field from the last time step of a transient heat transfer analysis can be used as a static
thermal load for structural stress analysis. This is done by specifying the load set ID with the TSTRU data.
Solid Element Properties Expanded to Support Navier Stokes Viscosity and an Advanced
Tangent Matrix Estimation
The PSOLIDX data has been enhanced with the items LAMBDAV, MUV, and IHKT to support Navier Stokes
viscosity and an advanced tangent matrix estimation.
Composite Element Material MATX25 Enhanced to Support Failure Based on the Ratio of the
Number of Failed Layers to the Total Number of Layers
The MATX25 data has been enhanced with the item RATIO, which can be used to specify the ratio of the
number of failed layers to the total number of layers that signifies failure of the element.
14
Thermal Conduction Based on Contact Surface Opening Distance and Closed Pressure
In Heat Transfer Analysis, there can be conduction across open contact surfaces based on the distance
between the surfaces. Also, the conduction across surfaces in contact can be specified to be a function of
the contact pressure. TABLEDi data is used to specify the relationship between the open surface distance
and the thermal conduction coefficient as well as the contact pressure and thermal conduction coefficient.
This is specified on the PCONTHT data. The solution of this contact based heat transfer analysis is
performed using an iterative solution of heat transfer and non-linear structural analysis (to determine the
contact status).
New Parameters
PARAM,SNAPTHRU: PARAM,SNAPTHRU,YES is used to control the generation of a ForceDeflection curve from Snap Through analysis. The default is NO.
TSTRU: In addition to using the resulting thermal field from a steady state heat transfer
analysis for a thermal load on a structural analysis, the thermal field from the last time step of
a transient heat transfer analysis can be used for the thermal load.
PRETBOLT: The OPTI option has been added to generate an ASCII file containing the bolt
pretension forces.
15
PGAPHT: Used to specify the conduction coefficients for the open and closed gap as well as a
tabled relationship between the gap opening distance and the thermal conduction coefficient.
PCONTHT: Used to specify the conduction coefficients for the open and closed contact as well
as a tabled relationship between the contact opening distance and the thermal conduction
coefficient. Also, a tabled relationship between the contact pressure and the thermal
conduction coefficient can be specified.
CONTACT and TIE: The ADJUST option has been added to ignore initial penetrations.
MGASK: The MGASK data now accepts TABLES1 data to specify temperature dependent
material properties.
MGASK: The GPLTYPE option can be used to specify if the GPL data is the stress per unit
displacement (default), or force per unit area.
MGASK: The EPLTYPE option can be used to specify if the EPL data is the tension stabilization
coefficient (default), or the direct tensile modulus.
PSOLIDX: The PSOLIDX data has been enhanced with the items LAMBDAV, MUV, and IHKT to
support Navier Stokes viscosity and an advanced tangent matrix estimation.
MATX25: The MATX25 data has been enhanced with the item RATIO, which can be used to
specify the ratio of the number of failed layers to the total number of layers that signifies failure
of the element.
16
Sets of RIGID Elements can now be Specified with the Bulk Data SET Data
The set TYPE can now be specified as RIGID in the Bulk Data SET data. This is required when specifying
the set of RIGID elements attached to the connection points in the PFPATH data.
17
Run Time and Disk Space Usage Decreased When GPFORCE and MPCFORCE Output is Requested
for Frequency Response and Transient Analysis
The run time and disk space usage has been decreased when GPFORCE and MPCFORCE output is
requested for Frequency Response and Transient Analysis. The reduction is quite dramatic for large
problems.
The Output From Complex Eigenvalue Analysis can be Reduced by Only Outputting the Unstable
Modes
Instead of outputting all modes from a complex eigenvalue analysis, the command DISP(UNSTABLE)=ALL
can be use to specify only the output of the unstable modes.
The Creation Time of CDS Super Elements Has Been Dramatically Reduced
A new algorithm for CDSMETH Super Element generation dramatically reduces the calculation time.
PARAM,FASTFR: The parameter FASTFR is used to turn on or off the use of the new parallel
fast equation solver for large modal frequency response analysis. The default value is AUTO.
When set to AUTO, OptiStruct automatically determines if this solution algorithm will reduce the
run time.
18
PARAM,RFIOUT: When this PARAM is set to YES during a CMSMETH run using the CBN
method, a file is created with the extension .rfi. This file contains the modal super element to
be used with the RecurDyn Multibody Dynamics Software from FunctionBay.
DISP: The UNSTABLE option specifies that only the unstable mode shapes are output from
Complex Eigenvalue Analysis.
MATF1, MATF2, MATF3, MATF8, MATF9, and MATF10: Material properties can now be
made frequency dependent using the MATFx data that reference TABELDx that contains the
material property vs. frequency data.
CMSMETH: Using the DMIGDV continuation of the CMSMETH data, the contents of an include
file used for optimization of the Super Element modal frequencies and their corresponding
structural viscous damping values is specified.
SET: The SET data can now be used to specify a SET of Rigid Elements. This is used for onestep Transfer Path Analysis by specifying this SET on the PFPATH data.
CDSMETH: The OSET field can be used to specify the set of interior GRID for the super
element in which results will be recovered in the residual run.
19
Analysis Results can be Generated Using the Results from SIMPACK Analysis
After running a multibody dynamic analysis in SIMPACK, the resulting modal participation coefficients in
the .unv file can be used by OptiStruct to recover the dynamic displacements, velocities, accelerations,
stresses, and strains. The ASSIGN statement keyword SIMINP is used to specify the .unv file.
ASSIGN,SIMPINP: Used to specify the .unv file that contains the CMS flexbody modal
participation coefficients generated from a SIMPACK multibody dynamic analysis.
CMSMETH: A LOAD LID and SPC SID can be specified on the PRELOAD continuation to preload
the flexbody. The effect of the preload on the flexbody is taken into account when the flexbody
matrices are calculated.
CMSMETH: A LOAD LID can be specified on the LOADSET continuation to generate residual
vectors that can increase the accuracy of the flexbody results. The degrees of freedom of the
loads are used to create residual vectors.
General
Gap Element Forces can Now be Plotted in HyperView
Note that the model file must the .fem file and not the .h3d file, because the h3d file does not contain
information about the gap coordinate system.
20
Printing of Inertia Relief Rigid Body Force and Acceleration Controlled by PARAM,PRINFACC
By default, the inertia relief rigid body forces and accelerations are no longer printed to the .out file. To
turn on this printing, set PARAM,PRINFACC to 1.
The GEOM3, DPDT and GPL Data Blocks have Been Added to .op2 File for PARAM,POST,-2
When PARAM,POST,-2 is used, the GEOM3, GPDT and GPL data blocks are now included in the .op2 file.
These data blocks are required by some in-house post-processors.
21
Thin Bending Shell Elements can Now be Automatically Converted to Membrane Elements
Thin shell elements (thickness 0.001 or less) can cause numerical issues such as singularities when
bending material properties are specified using MID2 and MID3. These elements are generally used to
skin solid elements for surface stress recovery purposes. They should be converted to membrane
elements by removing MID2 and MID3. This can now be done automatically using
PARAM,SHL2MEM,value, where value is the upper bound thickness at which the bending elements will
automatically be converted to membrane elements.
Overly Stiff Bushings and Springs can be Automatically Set to a Reasonable Value
The parameters BUSHSTIF and ELASSTIF are used to specify upper bounds of stiffness for CBUSH and
CELAS elements. Any stiffness value above the value specified by the parameter will automatically be
reset to the value of the parameter.
CBAR and CBEAM can Optionally Reference PBEAM and PBAR Data, Respectively
In the past, CBAR data could reference PBEAM data and CBEAM could reference PBAR data. Whether or
not this is allowed is now controlled by the SYSSETTING command BARPROP. If BARPROP is set to
STRICT (the default), then, CBAR and CBEAM cannot reference PBEAM and PBAR data, respectively. If
BARPROP is set to MIXED, then CBAR and CBEAM can reference PBEAM and PBAR data, respectively.
Von Mises Stress and Strain results for Random Response Analysis based on the Segalman
Method
The von Mises stresses and strains can now be output for Random Response Analysis based on the
Segalman Method. The existing PSDF and RMS options on the STRESS I/O options entry now includes von
Mises stresses and strains for Random Response Analysis subcases.
22
PARAM,RECOVER can be used to Output the Full-structure Modes During Component Mode
Synthesis
The new parameter RECOVER can be used to output the full-structure modes, within a range of
frequencies specified by LB (Lower Bound) and UB (Upper Bound), during Component Mode Synthesis
(CMS) instead of the modes of the condensed system.
PARAM,PRINFACC: The new parameter PRINFACC controls the printing of the inertia relief
rigid body forces and accelerations to the .out file. By default, the inertia relief rigid body
forces and accelerations are no longer printed to the .out file.
PARAM,AUTOMSET: When the new parameter AUTOMSET is set to YES, any rigid elements
with double dependency issues are fixed by the conversion of dependent degrees-of-freedom to
independent degrees-of-freedom.
PARAM,SHL2MEM: The new parameter SHL2MEM controls conversion of thin shell bending
elements to pure membrane elements. All bending shells with thickness equal or less than the
value of this parameter are automatically converted to membrane elements by the removal of
the MID2 and MID3 data from the corresponding PSHELL data.
PARAM,BUSHSTIF: The new parameter BUSHSTIF controls the upper bound of stiffness for
CBUSH elements. Stiffness values greater than BUSHSTIF are reduced to the value of
BUSHSTIF.
PARAM,ELASSTIF: The new parameter ELASSTIF controls the upper bound of stiffness for
CELAS elements. Stiffness values greater than ELASSTIF are reduced to the value of ELASSTIF.
PARAM,RECOVER: The new parameter RECOVER can be used to output the full-structure
modes, within a specified frequency range, during Component Mode Synthesis (CMS) instead of
the modes of the condensed system.
PARAM,PRGPST: The parameter PRTPST can now be used to control the amount of AUTOSPC
information printed to the .out file.
23
Topology Optimization
DTPL Data MAXDIM Automatically Set if the Input Value is Too Low
If the input value of MAXDIM is too small, it will automatically be reset and an INFORMATION message
issued.
Optimization Output
Design Sensitivity Analysis for Frequency Response in H3D
The sensitivities for displacement, velocity, and acceleration from frequency response analysis are now
output to the .h3d file. These sensitivities can be plotted vs. frequency in HyperGraph and are used by the
Design Sensitivity Analysis utility.
24
Design Variable File for the Control of CMS Super Element Modal Information
For CMS Super Elements generated using the General Method, a design variable file can be generated with
the DMIGDV continuation data of the CMSMETH Bulk Data. This ASCII design variable file can be included
in an optimization run to design the modal (stiffness, mass, and damping) characteristics of the Super
Element.
25
PARAM,TRAKMTX: The new parameter TRAKMTX controls the printing of the mode tracking
matrix to the .out file.
PARAM,TRAKMETH: The new parameter TRACKMETH is used to specify the Mode Tracking
Method.
PARAM,MFILTER: The new parameter MFILTER is used to specify the threshold for the mode
tracking eigenvector correspondence.
Resolved Issues
The accuracy of the PENTA element has been improved through a new element formulation.
In certain situations, discrete design optimization was not activated when TMANUF is specified
on the PLY entry. This issue has been fixed and discrete design optimization is now activated
upon specification of TMANUF.
A possible internal programming error during modal transient analysis runs with static loads,
using the Lanczos solver, has now been fixed.
When a transient thermal analysis was run with very small time steps, negative temperatures
at some grid points were observed. This issue has now been resolved.
The compliance of SUBCOM subcases that linearly combine static subcases (with SPCD/nonzero SPC loading) was incorrect. This issue has now been resolved by triggering element-wise
compliance calculations.
An improved error message is now output when IDENT is specified on the ACMODL entry and
the fluid-structure interface grid points are controlled by shape design variables with non-zero
initial perturbations.
OptiStruct has now been enhanced with long integer support, allowing for the use of memory
sizes above 16 GB. This resolves an issue where memory overflow error occurs if AMPFFACT is
set to a high value.
An issue where OptiStruct automatically resets the user-defined negative value of GE on the
CELAS entry to -1 has now been fixed.
ESLM-MBD topology optimization runs werent converging due to incorrect ESL calculations
when penalties change over iterations. This issue has now been resolved.
von Mises ply stresses and strains are now output to the H3D format when CSTRESS/CSTRAIN
output is requested to H3D.
During ESL (NLGEOM) optimization runs, the incorrect time step was picked. This issue has now
been resolved.
26
In a B2B run, ply stresses/strains are now converted from elemental coordinate system to ply
coordinate system and output to the H3D file. These ply stresses/strains in the ply coordinate
system can be used to calculate the composite ply failure indices.
The function of INTOL on the ACMODL entry has been improved to fix incorrect peak amplitudes
at certain frequencies in an acoustic cavity.
The SURF entry definition has been updated to resolve an issue when nodes on the section
surface for PRETENS definition (3D Pretension) are the same as the nodes used to define the
contact between the bolt surface and the wall of the hole.
MODALSE output is now available regardless of whether any standard output is requested.
Residual vector output was not output if a DLOAD bulk data entry is referenced by a DLOAD
subcase information entry using AMSES or AMLS. This issue has now been resolved.
A possible internal programming error during frequency response analysis runs with panels
defined using ERPPNL and PFPANEL output request, has now been fixed.
Resolved an issue where the entire model was exported to the SIMPACK (.fbi) file regardless of
the specifications on the MODEL card.
An issue where OptiStruct incorrectly output an error when the M2GG was specified within a
subcase has now been resolved.
Previously, the solver GUI listed the last used solver options and input file regardless of the
selected solver type. This issue has now been fixed and the Solver GUI will start with last used
input file/options for the selected solver.
An error in the SUBSEQ entry while reading a single positive number at the start of a
continuation line has now been resolved.
The number of rigid body modes output from an AMSES run was lower than expected. This
issue has now been resolved.
An issue where the tensor of complex stress results from a frequency response analysis was
incorrectly visualized using the H3D file has been resolved by enforcing default systems in the
solver.
An error was incorrectly output when elements different from shells were input using the ELIST
entry. This issue has now been resolved by removing the error and allowing non-shell type
elements on the ELIST entry and these elements are skipped from the analysis.
An issue where SPCD unloading effects are not properly accounted for in geometric nonlinear
analysis has now been resolved.
A programming error occurred, for some models, when AMSES was run out-of-core. This issue
has now been resolved.
27
Incorrect sensitivities were output during free-size optimization runs for ERP responses (on
skinned solids), where the design space was modeled with TRIA6 elements. This issue has now
been resolved.
Incorrect sensitivities were output during topology optimization runs for ERP responses, where
the design space was modeled with PENTA15 elements. This issue has now been resolved.
28
Introduction
RADIOSS is a state-of-the-art finite element solver uniting implicit and explicit integration schemes for the
solution of a wide variety of engineering problems, from linear statics and linear dynamics to complex
nonlinear transient dynamics and mechanical systems. This robust, multidisciplinary solver enables
designers to maximize performance related to durability, NVH, crash, safety, manufacturability, and fluidstructure interaction, in order to bring innovative products to market faster.
AIRBAGS UP
29
AIRBAGS FVM
/DT/FVMBAG time step stability control for FVM to reduce time step drop
/FVMBAG/MODIF to redefine along the engine phase the merging criteria and parameters
originally defined in the starter which is particularly useful for airbags with internal surfaces like
dual chamber side airbags
Airbags
Several enhancements to improve results quality and solution robustness:
1.
2.
In addition to the above, the option ACOM (t) allows you to scale the area between communicating
chambers as function of time or relative pressure.
CONTACT TYPE24
Penalty formulation contact designed to work with large time-step. Main features are: large timestep (it
cant drop more than a factor 1.4x); zero gap; improved(*) robustness even with models with initial
penetrations and intersections.
Known limitations in 12.0.210 patch:
30
MPI parallelization is available but not fully optimized yet: load unbalancing might deteriorate
scalability on large number of cores
CONTACT TYPE7 Self Contact with Gap Larger than Mesh Size
Resolves self contact when mesh size is smaller than gap value; neighbor's nodes to a master
segment are automatically removed from the list of contact candidates
Icurve: now available for both concave and convex master surfaces
Incremental (stiffness) formulation allows larger time step and improves stability at the same
time
SPMD parallelization of the rad2rad has been achieved showing tremendous speed-up for high
number of cores compared to previous SMP rad2rad. It can also run in Hybrid SMP/SPMD mode.
Several post-processing improvements are also available, especially for TH files to make output
consistent between main and subdomain T01 files. Sections are now available for output.
XFEM (beta)
12.0.210 includes a brand new and more powerful X-FEM formulation; nevertheless actual implementation
has not been fully tested yet. Therefore it should be considered a beta option.
31
XFEM technology is available for multilayered 3 and single layer (PID1) shells but not mixed
together, yet.
The format was changed and simplified for a new XFEM formulation compatible with failure
criteria /FAIL: Teuler-Butcher /TBUTC; Johnson-Cook /JOHNS; Forming Limit Diagram FLD and
tabulated /TAB.
Known limitations:
PARITH/ON is not ensured for SPMD parallel mode
LEVSET is not plugged yet with the new XFEM formulation
32
Resolved Issues
Wrong direction for alpha values output in .sta files (LAW58 + PROP16).
/LAW59: nodes belonging to deleted elements still set the time step
TYPE2 + Spotflag=25: possible instability when distance between slave node and master
segment is large
/INTER/TYPE7: kinematic time step was computed using always Gapmin, ignoring variable gap
options
33
CONTACT TYPE21
Improved performance
Resolved Issues
LAW73: engine stop when extrapolation from tables was giving yield < 0; now Radioss always
consider yield >= 0
Porous material law to take into account the effect of porosity on the equation of states
MATERIAL LAW51
/ANIM/BRIC/VFRAC - allows to track material evolution; its targeted to replace global density
contour which was not always adequate.
/ALE/DONEA
34
/ALE/DISP
/ALE/SPRING
/ALE/ZERO
/ALE/STANDARD
Resolved Issues
LAW51: multi-material outlet doesn't work when detonation occurs in air
RADIOSS - General
SENSORS Improvements
Activation/deactivation of CONTACTS
TRANSFORM MATRIX
Compatible with transformations actually available in CAD systems and or multi-body models.
35
The list of nodes with highest velocity is output in case of mass or energy error. Scope is an
easier identification of the nodes (parts) responsible for the divergence.
/DAMP/INTER
Sensitivity analysis
Integrated Beams
In /PROP/TYPE18 it is possible to refer to predefined sections (rectangular or circular)
AMS
Performance improvement for /DT/INTER/AMS
Compatibility between /DT/AMS and /DT/NODA/CST allowing the application of AMS only to a group of
parts and a classical time step control to the remaining parts of a model for optimized global
performances.
Resolved Issues
Domain decomposition fails if NUMSPH*KVOISPH > 2exp31 for large SPH models
Running in SPMD a first job with /DT/NODA/CST and then switching to AMS caused engine
failure
36
37
All Solvers
Performance Enhancements
HyperXtrude Solvers can now handle large models on both Windows and LINUX 64-bit platforms. In the
previous release, size was limited to 450,000 nodes on Windows machines. Now the solver can handle up
to 4 million nodes on both the platforms. In addition to this, there is a significant speed improvement.
Metal Extrusion
Improvements to Solver Accuracy
Accuracy of the nose cone prediction of profiles, especially those with large variation in thicknesses, is
improved.
Mandrel Offset
The recommended procedure is to flush the mandrel offset before creating the mesh. This is accurate and
also will result in good quality mesh in the welding chamber region. In this release, if desired, user can
also mesh the model with the mandrel offset and solver will account for this step in the computation and
the computed solution will be equally accurate. The advantage of including mandrel offset is savings in
38
time spent in CAD correction; however, there may be a disadvantage in creating good quality mesh in this
region.
Resolved Issues
Error in the units printed for the surface area of boundary faces in the output balance tables in now
corrected.
Billet Forging
Grain Size Computations
Solver now supports computation of grain size using static, dynamic, and metadynamic recrystallization
models.
Material Models
Flow stress of the workpiece can now be computed using the following constitutive models.
Johnson-Cook
Norton-Hoff
Hensel-Spittel
Flow stress data in 3D table format (as function of strain rate, strain, temperature)
In addition, users have the choice of using the existing Sellars-Tegart and metal forming power law
models.
39
Resolved Issues
A bug in resin front calculation is now resolved.
40