12-19-2017 10:10 AM
Hello,
Refer to the attached images. For same geometry, loading and boundary conditions; I get different values of max. nonlinear von Mises stress from nonlinear static and advanced nonlinear static anaylsis. Do you know why this is so?
Regards,
Su
12-19-2017 06:03 PM
Dear Su,
Without having the two FEMAP models in hand is difficult to know the details used in material model, analysis options, etc.. Please note the diferences are notable between NX NASTRAN basic nonlinear (SOL106) and advanced nonlinear (SOL601), aka ADINA. For instance, depending the resultant strain level the differences could be important.
Best regards,
Blas.
12-20-2017 04:53 PM
Dear Blas,
Thanks for your detailed response on SOL106 vs SOL601!
I have attached the Femap model. Could you please let me know if the difference in maximum stresses is due to some particular analysis settings in this case? The strain values are very close to each other.
Also, can I plot the stresses at Gauss points in Femap v11.2 using GUI?
Best regards,
Su
12-22-2017 02:37 PM
Dear Su,
I see some improvements are recommended to be performed in your FEMAP model, and also I found some errors as well, here you are:
HEX MESH
Here is extremely critical to mesh with HEXAEDRAL elements, specially in nonlinear analysis, forgot tetrahedral elements because TET are less accurate and the model size is in average around 10 times bigger for the same element size: 49437 nodes your TET10 model, 6872 nodes mine HEX8 only, size matter, specially in nonlinear analysis!!. I have recorded a video to explain how to mesh the geometry with solid HEX elements (I hope you love the Country music ...):
ERROR IN LOADINGS
I see you have applied a Moment load in surface, wrong!!, solid elements don't have rotational DOF, then is useless to apply a moment on a surface of a solid meshed with 3-D solid elements.
You need to apply a TORQUE load on surface. The command will ask you a vector defined by two points (base & tip) to apply the torque load. Once applied you can expand the load (use MODEL > LOAD > EXPAND command) to see the resulting force loads, this is the way how FEMAP will export the analysis file to the NX NASTRAN solver:
YOUR MATERIAL MODEL
I see you have defined a bi-linear material model with vonMises yield stress = 172 MPa and Work hardening slope H=1949 MPa, that is a slope of stress vs plastic strain.
My option here in general is clear: with STEEL because most of the times the Plasticity Modulus H in general is not know exactly, I always select PLASTIC model with vonMises Yield criterion (in your case yield stress= 172 MPa), which means that after reaching yield the line on the stress-strain chart is horizontal, ie, an elastic-perfectly plastic nonlinear material model.
YOUR NONLINEAR ANALYSIS PARAMETERS (SOL106)
I note in the BULK DATA Options that PARAM=LGDISP,1 is not activated, wrong!!. This means you are performing a linear geometric analysis (small displacements), the only nonlinearity is material. Please be aware of this error, this is critical to activate large displacements effect, then nonlinear geometric analysis will be performed, OK?.
Also, the most crucial data for successful nonlinear static solutions are contained in the NLPARM bulk data entry. NLPARM defines strategies for the incremental and iterative solution processes, surprise you have not setup any options!!.
Here you are my settings:
And here you are the nodal vonMises stress results of the nonlinear static (SOL106) analysis:
ADVANCED NONLINEAR (SOL601)
Again I note in the BULK DATA of your FEMAP model the parameter for LARGE DISPLACEMENTS EFFECT is not activated, then your nonlinear results are useless again. Here you are MY OPTIONS:
If I run the analysis here you are the results:
RESULTS COMPARISON
Well, not too much diferences, with SOL106 the material plastify locally (maximum displacement URES=0.0509 mm), with SOL601 the material seems to remain in elastic regime (maximum displacement URES=0.0509 mm), and the maximum resultant displacement is the same.
OUTPUT AT GAUSS POINTS
Regarding the question STRESS/STRAIN OUTPUT at GAUSS POINTS WITH SOLID ELEMENTS, it was explained here, apparently with low interest by the FEMAP & NX NASTRAN community:
Well, I hope the above explanation to be of help to study more in detail all the parameters of the nonlinear analysis of both Basic Nonlinear (SOL106) and ADVANCED NONLINEAR (SOL601) modules, OK?.
Best regards,
Blas.