SimSolid Real World Validation Manual
SimSolid Real World Validation Manual
Stress Accuracy
Richard King 04/09/19
Abstract
This manual contains a series of SolidWorks models with static simulation studies that were run in Altair SimSolid™.
Some of them are complex, but an accurate solution was obtained in SolidWorks by the conventional procedure of
placing mesh controls in the regions of highest stress gradient and successively refining until the solution converged.
This provides a suite of accurate stress results to compare SimSolid’s results against. In all cases SimSolid’s stress results
compared to the SolidWorks comparison results within 10 percent.
Introduction
SimSolid is a very easy to use and fast product for structural analysis using novel alternative technology- a meshless
approach that based on similar theory to FEA. This approach can be shown theoretically to converge to correct
solutions, as discussed in the whitepaper on the technology [1], where it is shown that the mathematical foundation for
the method used in SimSolid is on as solid of a foundation as FEA. But in stress analysis the quantity being calculated is
the derivative of the directly calculated variables (displacements) so is more difficult to achieve accurate results. This is
true of FEA codes as well as SimSolid. In FEA codes, the solution is to have a fine enough mesh in regions of high stress
gradients. This can sometimes be done adaptively using automatic mesh refinement, but often requires one or more
user interventions, for example to create mesh controls, as described in appendix 1. On complex models this can be
tedious and time consuming.
SimSolid has an automated approach to address this issue. In a first step, the model is analyzed with automatically
chosen displacement functions to accurately calculate displacements, assure continuity at all interfaces and satisfaction
of boundary conditions. This captures overall load paths and reaction forces accurately [2]. In the second step, it adapts
to local features where there are high stress gradients by introducing additional functions it those areas. From a user’s
standpoint, this “two stage workflow”, as it is called in the training material, is very easy:
1. Assign materials, loads, and constraints directly to the CAD model. Hit the “Solve” button to run the model with
defaults settings. This run is typically very fast. When complete, examine the results to make sure the
deformation looks reasonable.
2. Turn on solution settings to adapt to features by clicking a single checkbox shown in the figure. If the model
contains thin shell or sheetmetal type parts, an additional button, “adapt to thin solids” can be checked. This is
all done for the entire model for a single part model. For an assembly, the user picks specific parts in which high
stress has been observed to use for adapting to features or thin solids.
I have used various finite element codes with adaptivity, including SolidWorks Simulate, Autodesk Fusion 360, and
Mechanica (of which I was co-inventor and one of the authors). In my opinion SolidWorks Simulate is a good quality
commercial code, with a robust automesher, so it was a good choice for comparing the accuracy of SimSolid. But all FEA
codes have the drawback of difficulties associated with meshing, so none of these is as easy to use in practice as the two
stage workflow for SimSolid described in the Altair SimSolid Fast Start Training Guide [2].
Model Location
All the models described in this document are in the folder “valmanModels”. There is a SimSolid workspace (.ssp) file for
each, as well as SolidWorks part and/or assembly files.
Summary
Solution Timing: Solution time was similar for all the smaller models. The last 3 models are relatively large and the timing
for SolidWorks vs SimSolid is described below. For these models the total solution time (including meshing time in
SolidWorks) is significantly less in SimSolid. For all models, considerably less user intervention was require in SimSolid to
achieve an accurate result.
Single Part Models
Model: pwh
Description: Plate with hole
Model: hanger-beam
Description: hanger-beam model from SimSolid Tutorial.
The default mesh from commercial FEA codes does a pretty good job of following the geometry of the model and will
often produce accurate displacement results. Stress results are more localized and may not be accurate unless manual
or adaptive mesh refinement is done in local regions of high stress gradient.
The following settings for an automesher are crucial as a prerequisite for achieving accurate stresses. They are turned on
by default in Solidworks Simulate, but not always in other codes (for example they are off by default in SolidEdge
Simulate), because they can make the automesher fail. But if they are not on there is little chance of accurate stresses.
• Project midside nodes on surfaces. If this is not on you may be modeling the wrong physical model, such as a
plate with a hexagonal hole instead of a circular hole.
• Curvature-based mesh refinement: more elements should be placed along highly curved edges of the model.
Here is an example of a decent default mesh near a curved cutout (this is actually the 14th validation model, the testing
machine assembly with a notched specimen:
This is a reasonable mesh. The geometry is followed well, there are more elements where the curvature is higher, and
the elements near the feature (the curved end of the notch) are not large compared the same size of the feature. They
are pretty large for a region of high stress gradient, however, which raises a red flag from the standpoint of accurate
stresses: elements used in FEA codes represent the displacements quadratically, so stresses, which are related to the
derivatives of displacement, are approximately linearly. But the stresses near a curved feature will vary sharply,
something like shown in the red line below, with distance from the edge of the notch:
Here is the mesh obtained using the manual mesh refinement procedure described below:
Now the high stress variation now occurs across several elements.
3. If the adaptive stress results agree well with the default, this is probably an accurate result. Of the 19 models in
this manual, that was true for only one, the plate with a hole. Otherwise further manual refinement is needed.
In the testing machine assembly, for example, the default solution and adaptive solution disagree by 12%, so I’d
recommend checking further with manual refinement.
4. Examine the mesh near the vicinity of maximum stress. Place a mesh control there, for example, the notch
surface:
Set the mesh density to high or “fine” in the control and remesh. Make sure the new mesh looks reasonable in
that vicinity:
If not, you can manually tighten the mesh density more. Once the mesh looks reasonable, run the model again.
In this case the first mesh from the mesh control still doesn’t look good with fairly large elements still adjacent
to the notch, so I reduce the element size on the control and try again:
The 5 step procedure I described was used for all the FEA results in this manual.
Now let’s look at this procedure in detail using the same example used in appendix 1, the testing machine assembly.
After solving the default analysis (“stage I”) we see the max stress in the model is in the notch of the specimen:
Now the procedure is to turn on adapt to features, but only for the part or parts where the stress is highest:
Note there is a pulldown for refinement level: standard, increased, or fine. The higher this level, the more accurate the
stresses. I set it to fine for all the models in this manual. Since we’re only doing this for a small subset of the model, it
does not add to the solution time too much. The other option on this form is “Max number of adaptive solutions” which
defaults to 3. I recommend leaving it at 3 and re-running:
If this agreed with the default solution within a few percent, we’d be confident in an accurate solution, but it is off
considerably. So now I recommend bumping up the Max number of adaptive solutions to 4 and trying again:
This did not change the result much so we can accept this result. Otherwise it would be necessary to adjust the number
of solutions again and repeat the process.
This procedure was done for all the models in this manual. Out of 19 models, one needed only 3 for the number of
solutions, thirteen needed 4, two needed 5, and two needed 6.
They can be caused by geometry, such as reentrant corners, or by loads, constraints, and special elements. Here are
some examples:
Singularities caused by loads can be avoided in SimSolid by applying the load to a “spot” or by using remote loads.
Reentrant corners (or similar situations where there is a kink where two surfaces meet) can be resolved with fillets.
If the maximum stress in your model occurs at a singularity it is best to rethink the true physical situation you are trying
to model and resolve it. For example, tube joints do not meet at an infinitely sharp corner, in reality there should be a
fillet there. If you are interested in the stress state elsewhere in your model, you may need to use probes to evaluate the
stress at that other location.
About The Author
Richard B King, Ph.D.
I received my Ph.D. in Applied Mechanics from Stanford University in 1980. My training included theoretical stress
analysis, numerical methods including the boundary element and the finite element methods, and variational methods
such as the Ritz and Galerkin methods (described in the SimSolid Technology Overview [1[, these are the predecessors of
both the finite element method and the external method used in SimSolid). I have decades of experience in stress
analysis.
In the early part of my career I used all of these methods in investigating the mechanical reliability of pipelines and naval
vessels at the National Bureau of Standards, and of disk drives at IBM research. At IBM it became clear that the state of
the art of finite element analysis was found to be clearly inadequate for design purposes, which led to the idea of using
the adaptive p-method of finite element analysis for this purpose. This led to my co-invented Mechanica and cofounding
Rasna Corporation.
Rasna was acquired by PTC in 1995, and I continued to work on enhancing the adaptive analysis capabilities in
Mechanica (later renamed Creo Simulate). This included advanced adaptive nonlinear features including snapthrough,
nonlinear thermal, large strain plasticity, large displacement and frictional contact.
After leaving PTC I did freelance consulting and teaching, and developed stressRefine, an accessory for p-adaptive
solution of existing FEA meshes, which is now available open-source. I also investigated the capabilities of commercial
FEA adaptive mesh refinement solutions include SolidWorks Simulation, Ansys, and Autodesk Fusion 360 Simulation.
More recently I have looked into SimSolid.
References
1. Altair SimSolid Technology Overview, Altair Corporation,
https://altairhyperworks.com/ResourceLibrary.aspx?category=Technical%20Papers&altair_products=SimSolid
2. Altair SimSolid Fast Start Training Guide,
https://altairhyperworks.com/ResourceLibrary.aspx?category=Technical%20Papers&altair_products=SimSolid
3. Timoshenko, S, and Goodier, J, Theory of Elasticity, McGraw Hill, 1970
4. Pilkey, W. D., Pilkey, D. F., & Peterson, R. E. (2008). Petersons stress concentration factors. Hoboken, NJ: John
Wiley.