#112 MSC Nastran SOL 106<\/p>\n
MSC Nastran SOL 106<\/p>\n
The nonlinear effects in structures occur due to nonlinear material behavior and large deformations. Geometric nonlinearity becomes relevant when the structure is subjected to large displacement and rotation. Geometric nonlinearity effects are prominent in two aspects: geometric stiffening due to initial displacements and stresses, and follower forces due to a change in loads as a function of displacements. Material nonlinear effects may be classified into many categories: plasticity, nonlinear elasticity, creep, viscoelasticity. The solution 106 uses an iterative process based on the modified Newton\u2019s method. The solution operations includes gradual load or time increments Delta P, iterations with convergence tests for acceptable equilibrium error (R1, R2, …) and stiffness matrix updates.<\/p>\n
106 Solution Example:
\nSOL 106
\nCEND
\nTITLE = Test
\nDISP = ALL
\nSTRESS = ALL
\nSPC = 1
\nSUBCASE 1
\nSUBTITLE = Elastic – load at 85%
\nLABEL = Load in Elastic Field
\nLOAD = 50
\nNLPARM = 50
\nSUBCASE 2
\nSUBTITLE = Elastic – load at 100%
\nLABEL = Load in Plastic Field, the yield limit is passed
\nLOAD = 100
\nNLPARM = 100
\n$END OF CASE CONTROL DATA
\n$BULK DATA SECTION
\nBEGIN BULK
\nNLPARM 50 1 AUTO UPW NO
\nNLPARM 100 8 SEMI UPW NO
\n…
\n…
\n…
\nENDDATA<\/p>\n
Recommendations:
\n1) A simple model to understand the behavior of the structure.
\n2) The model has to be a compromise between accuracy and the efficiency of the analysis.
\n3) A finer mesh is advisable in the area of interest.
\n4) Avoid stiff elements like TRIA3 and TET4, the solution convergence could be more difficult.<\/p>\n
Nonlinear static analysis permits only one independent loading condition per run, while in a linear analysis subcases represent an independent loading condition and each subcase is distinct from all others. A non linear analysis is path dependent , the loads and position are incremented step by step, the changes are in the subcase are cumulative. The NLPARAM is the parameter which controls the nonlinear analysis, it defines strategies for the incremental and iterative solution processes. It is difficult to choose the optimal combination of all the options for a specific problem. The default option was intended to provide the best workable method for a general class of problems, users should start with the default option. When by default parameters the convergence is not reached, by some changes it is possible to get the solution.<\/p>\n
NLPARM Card
\n1 1 NLPARM
\n1 2 ID
\n1 3 NINC = number of increments
\n1 4 DT
\n1 5 KMETHOD = AUTO, SEMI, ITER
\n1 6 KSTEP
\n1 7 MAXITER = number of iterations allowed for each load increment.
\n1 8 CONV
\n1 9 INTOUT
\n2 1 nothing
\n2 2 EPSU
\n2 3 EPSP
\n2 4 EPSW
\n2 5 MAXDIV
\n2 6 MAXQN
\n2 7 MAXLS
\n2 8 FSTRESS
\n2 9 LSTOL
\n3 2 MAXBIS
\n3 6 MAXR
\n3 8 RTOLB
\n3 9 MINITER<\/p>\n
The MAXITER field is an integer representing the number of iterations allowed for each load increment. If the number of iterations reaches MAXITER without convergence, the load increment is bisected and the analysis is repeated. If the load increment cannot be bisected (i.e., MAXBIS is reached or MAXBIS = 0) and MAXDIV is positive, the best attainable solution is computed and the analysis is continued to the next load increment. If MAXDIV is negative, the analysis is terminated. MAXDIVis the limit on probable divergence conditions per iteration before the solution is assumed to diverge, it is based on a energy rate error by means an iteration is defined as diverged.<\/p>\n
The convergence test is performed at every iteration with the criteria specified in the CONV field. Any combination of U (for displacement), P (for load), and W (for work) may be specified. All the specified criteria must be satisfied to achieve convergence.<\/p>\n
INTOUT controls the output requests for displacements, element forces and stresses, etc. YES: Gives all outputs at every increment
\nNO: Gives output at last increment or load step. It is advisable to use YES to now the status step by step. It is useful to get results with all the steps used for the load increments.<\/p>\n
The convergence tolerances are specified in the fields EPSU, EPSP, and EPSW for U, P, and W criteria, respectively.<\/p>\n
Nonlinear iteration method:<\/p>\n
1) The AUTO method works very well as a starting point and is the default method in the nonlinear static solution. This method essentially
\nexamines the solution convergence rate and uses the quasi-Newton, line search and\/or bisection methods to perform the solution as efficiently as possible sometimes without a stiffness matrix update. Highly nonlinear behavior in some cases may not be handled effectively using the AUTO method.<\/p>\n
2) The SEMI method is also an efficient technique for nonlinear static iteration. It is similar to the AUTO method since it uses the quasi-Newton, line search and\/or bisection methods, but it differs from the AUTO method in one respect. The SEMI method forces a stiffness matrix update after the first iteration of a load increment. This update is performed irrespective of the convergence status of the solution. It is effective in many highly nonlinear problems where regular stiffness matrix updates help the solution to converge.<\/p>\n
3) The ITER method is known as the \u201cbrute force\u201d method in nonlinear static analysis because a stiffness matrix update is forced at every KSTEP-th iterations with a resultant increase in CPU time. Its best applications tend to be those involving highly nonlinear behavior for which using only multiple iterations may not be efficient.<\/p>\n
MAXBIS and Convergence Criteria:<\/p>\n
The incremental load value is bisected automatically if the convergence is not achieved in any incremental load step even after maximum number of iteration defined by the user. Bisection means cutting the time step size or load increment size by half. This bisection occurs mainly due to any large nonlinearity in the model or when the NINC or load increment for the particular load step is too large. Divergence processing uses a combination of stiffness updates and bisections to facilitate convergence. Note: MAXBIS is not written out directly from Patran. It has to be manually included if the default values need to be overwritten.<\/p>\n
Convergence criterion may be any combination of out-of-balance forces, change in displacements, and energy (work) error. Convergence criterion based on out-of balance force and change in displacement is used before yielding and convergence criterion based on out-of-balance force is used after yielding because displacement is very high after yielding. If a solution for particular load increment (or time step) cannot converge within MAXITER number of iterations, that increment is halved. Solution is attempted again, and if it does not converge again, another halving (bisection in NASTRAN speak) is done, up to a maximum number of bisections controlled by MAXBIS entry on NLPARM entry. Default is 5, max value is 10, but if you have to go higher, it most likely indicates either some irregularity with the model, or that NASTRAN is not the best tool for problem at hand.<\/p>\n
1) Increase NINC from default 10 to 15 or 20
\n2) Increase MAXBIS, up to 10
\n3) Change KMETHOD to ITER or SEMI and set KSTEP to 1 (caution: this will slow down the solution significantly, sometimes it is better to use the default value).<\/p>\n
Example:
\nNLPARM 1 20 AUTO 5 30 P YES +
\n+ 0.01<\/p>\n
To make that the solution is more fast, put this in the parameter section (it requests that the follower force stiffness will not be considered):
\nPARAM FOLLOWK NO<\/p>\n
Nonlinearities:<\/p>\n
Material Definition \u2013 MATS1 and TABLES1. To introduce the material nonlinearity 2 ways are generally adopted:
\n1) by H parameter or
\n2) MATS1 + TABLES. Elastoplastic material consists of elastic and a plastic portion in stress-strain curve. The elastic portion is defined by providing the Young’s Modulus, E and Poisson’s ratio, n (MAT1); and the plastic portion is defined using a stress-strain curve data points through tabular input (MATS1 + TABLES1 or ‘FIELDS’ in Patran).<\/p>\n
A couple of important points to consider while defining stress-strain curve is that the first point must be at the origin (X1 = 0, Y1 = 0)
\nand the second point (X2, Y2) must be at the initial yield point specified on the MATS1 entry. The slope of the line joining the origin to the yield stress must be equal to the value of E. Sometimes when the Stress-Strain curve is not available, is convenient to use the “Hardening slope” H parameter.<\/p>\n
Geometric Nonlinearity:
\nWe should introduce the following parameters:
\nPARAM LGDISP 1
\nPARAM FOLLOWK NO<\/p>\n
LGDISP = 1 all nonlinear element types will be assumed to have large displacements (updated coordinates and follower forces).
\nFOLLOWK = NO requests that all the follower force stiffnesses not be included.<\/p>\n
Postprocessing:
\nTwo output blocks in SOL106:<\/p>\n
1) Nonlinear stress\/strain (for non linear element only):
\nThis is printed as \u201cNONLI EAR STRESSES IN TETRAHEDRON SOLID ELEMENT S (TETRA)\u201d in the *.f06 and it is given by default case control
\ncommand NLSTRESS = ALL. The Patran stress output label NONLINEAR comes from this table. This is the \u201cTrue stress.\u201d (This table also contains strains, the strains labeled NONLINEAR come from this output.) The nonlinear stress data is output in the element coordinate system.<\/p>\n
2) Linear stress\/strain (for linear and nonlinear elements). This is printed as \u201cSTRESSES IN TETRAHEDRON SOLID ELEMENTS (CTETRA)\u201d in the *.f06
\nand it is given by default case control command STRESS = ALL. In Patran it is shown as \u201cSTRESS TENSOR\u201d. In short it is the \u201cequivalent linear stress\u201d. You can also ask for STRAIN=ALL, this is the STRAIN TENSOR.<\/p>\n
The stresses at center should match in both cases, however due to extrapolation, linear and nonlinear may not match. Many times the discrepancy is also caused by averaging stresses at node (in fringe plot).<\/p>\n
And sometimes if you don\u2019t have perfect elastic plastic the stresses could be higher than yield point as well. You can have a look at the
\n\u201cnonlinear stresses, stress tensor\u201d with the centroidal values. It will give you the best correlation with hand calculations.
\n[Nonlinear stresses are given for nodal as well as centroid]\n
Nastran nonlinear stress request (NLSTRESS) contains both \u2013 stress and strain. Whereas for linear elements (e.g.: CBAR), one has to request
\nstress and strain separately (STRESS=ALL\/STRAIN=ALL).<\/p>\n
\u201cSTRESS TENSOR\u201d is applicable for all the elements \u2013 (linear and nonlinear elements).
\n“NONLINEAR STRESSES, STRESS TENSOR” is applicable for nonlinear elements ONLY (e.g.: BEAM, TETRA (4&10)).
\n“NONLINEAR STRESSES, EQUIVALENT STRESSES” is applicable to Von-Mises stresses in Non-Linear elements.<\/p>\n
The results that you see marked as \u201cStress Tensor\u201d in Patran from your nonlinear run are linear stresses based on nonlinear displacements, and the results marked as \u201cNonlinear Stresses, Stress Tensor\u201d are the true nonlinear stresses. The ones marked \u201cNonlinear Stresses, Equivalent Stress\u201d are a von Mises stress calculation”.<\/p>\n
\u201cSTRAIN TENSOR\u201d is same as stress tensor, but Strain Values \u2013 (for both linear and nonlinear elements).
\nPLASTIC STRAIN is strain beyond the elastic limit.<\/p>\n
When post-processing strain results in PATRAN, the results menu under \u201ccreate Fringe plot\u201d allows for two strain values the user can plot:
\n1) Nonlinear strains, plastic strain.
\n2) Nonlinear strains, strain tensor.<\/p>\n
Choice (1) is consistent with the value from the .f06 file under the EFF. STRAIN\/ PLAS\/NLELAS column from the NLSTRESS case control request. Choice (2) is consistent with the value from the .F06 file under STRAIN output you got from STRAIN=ALL case control command.
\nThe results that you see marked as \u201cStress Tensor\u201d in Patran from your nonlinear run are linear stresses based on nonlinear displacements,
\nand the results marked as \u201cNonlinear Stresses, Stress Tensor\u201d are the true nonlinear stresses. The ones marked \u201cNonlinear Stresses,
\nEquivalent Stress\u201d are a Von Mises stress calculation.\u201d<\/p>\n
Convergence info in f06 file:
\nEUI: Normalized error in the displacement.
\nEPI: Normalized error in the load vector.
\nEWI: Normalized error in the energy.
\nLAMBDA: Rate of convergence. Solution is diverging if it is > 1.0.
\nDLMAG: Absolute norm of the load error vector.
\nFACTOR: Scale factor for line search method.
\nE-First: Initial error before line search begins.
\nE-FINAL: Final error after line search terminates.
\nN-QNV: Number of quasi-Newton correction vectors to be used in the current iteration.
\nN-LS: Number of line searches performed.
\nENIC: Expected number of iterations for convergence.
\nNDV: Number of occurrences of probable divergence during the iteration.
\nMDV: Number of occurrences of bisection conditions due to excessive increments in stress and rotations.<\/p>\n
Restriction and Limitations in SOL 106:
\n1) CROD, CBEAM, CGAP, CQUAD4, CQUAD8, CTRIA3, CTRIA6, CHEXA, and CTETRA elements may be defined with material or geometric nonlinear properties.
\nAll other elements will be treated as having small displacements and linear materials.
\n2) Constraints apply only to the nonrotated displacements at a grid point. In particular, multipoint constraints and rigid elements may cause problems if the connected grid points undergo large motions. However, also note that replacement of the constraints with overly stiff elements can result in convergence problems.
\n3) Large deformations of the elements may cause nonequilibrium loading effects. The elements are assumed to have constant length, area, and volume, except for hyperelastic elements.
\n4) Since the solution to a particular load involves a nonlinear search procedure, solution is not guaranteed. Care must be used in selecting the search procedures on the NLPARM data. Nearly all iteration control restrictions may be overridden by the user.
\n5) Follower force effects are calculated for loads which change direction with grid point motion.
\n6) Every subcase must define a new total load or enforced boundary displacement. It is necessary because error ratios are measured with
\nrespect to the load change.<\/p>\n
Common Errors with SOL 106:
\n1) Slope of elastic curve in MATS1\/TABLES1 & Young\u2019s Modulus, E in MAT1 should be same.
\n2) Avoid having large load steps. E.g.: If the total load to be applied is say 10000N, break it into 2 or preferably 3 subcases\/loadsteps.
\n3) In case that the analysis doesn\u2019t reach the convergence, sometimes could be due to an error in the loads applied. For example, in case of pressure applied, if it is wrong, too much high in some zone, the problem is unrealistically too nonlinear.
\n4) Since the units of convergence tolerance is independent of units, avoid tightening or loosening the tolerance in the first run.
\n5) Make sure enough system resources (memory and HDD space) are available at hand for large models.
\n6) \u201cUser Fatal Message 4676 (NMEPD) \u2013 ERROR EXCEEDS 20.00 PERCENT OF YIELD STRESS IN ELEMENT ID=XX\u201d may be encountered if there exists
\nlarge deformation of the nodes of an element or presence of severely distorted elements. In case of large deformation (caused due to extremely large loads), then it is basically a large strain problem (as opposed to small strain assumption in Nastran) which needs to be solved by another solution sequence like MSC Marc or SOL 600. In case of severely distorted element, it is advisable to inspect the element quality in the mesh of the model.
\n7) “USER FATAL MESSAGE 1221 (GALLOC) \u2013 THE PARTITION OF THE SCRATCH DBSET USED FOR DMAP-SCRATCH DATABLOCKS IS FULL”. This fatal message is
\ncommonly encountered due to low system resources. Note that this fatal will stop writing the results into xdb\/op2 completely and no output will be obtained. To avoid this fatal message, include \u201cNASTRAN SYSTEM(151)=1\u201d and increase the buff size \u201cNASTRAN BUFFSIZE=65537\u201d in the NASTRAN SYSTEM CELL SECTION of input file (above Executive Control Section).
\n8 ) When you use for the material the Hardening method, if you get this error: *** USER FATAL MESSAGE 5423 (SADD5) ATTEMPT TO ADD INCOMPATIBLE MATRICES, it means that the H is too small.<\/p>\n","protected":false},"excerpt":{"rendered":"
#112 MSC Nastran SOL 106 MSC Nastran SOL 106 The nonlinear effects in structures occur due to nonlinear material behavior and large deformations. Geometric nonlinearity becomes relevant when the structure is subjected to large displacement and rotation. Geometric nonlinearity effects are prominent in two aspects: geometric stiffening due to initial displacements and stresses, and follower […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"nf_dc_page":"","_et_pb_use_builder":"","_et_pb_old_content":"","_et_gb_content_width":"","_lmt_disableupdate":"yes","_lmt_disable":"","jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[22],"tags":[42,121],"class_list":["post-666","post","type-post","status-publish","format-standard","hentry","category-nastran","tag-nastran","tag-sol-106"],"modified_by":"gantovnik","jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"jetpack_shortlink":"https:\/\/wp.me\/p8bH0k-aK","jetpack_likes_enabled":true,"jetpack-related-posts":[{"id":936,"url":"https:\/\/gantovnik.com\/bio-tips\/2021\/06\/168-nastran-sol-106-vs-sol-400\/","url_meta":{"origin":666,"position":0},"title":"#168 Nastran SOL 106 vs SOL 400","author":"gantovnik","date":"2021-06-22","format":false,"excerpt":"#168 Nastran SOL 106 vs SOL 400 It is always challenging when it comes to nonlinear analysis, either it is geometric, material, or contact nonlinearity using SOL 106 in Nastran. For Advanced nonlinear analysis, SOL 400 is a new approach to solve the different types of nonlinear analysis. SOL 400\u2026","rel":"","context":"In "buckling"","block_context":{"text":"buckling","link":"https:\/\/gantovnik.com\/bio-tips\/category\/buckling\/"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":632,"url":"https:\/\/gantovnik.com\/bio-tips\/2020\/09\/question-1-param-cards-in-nastran\/","url_meta":{"origin":666,"position":1},"title":"#99: PARAM cards in Nastran","author":"gantovnik","date":"2020-09-25","format":false,"excerpt":"#99: PARAM cards in Nastran 1) K6ROT specifies the stiffness to be added to the normal rotation for CQUAD4 and CTRIA3 elements. This is an alternate method to suppress the grid point singularities and is intended primarily for geometric nonlinear analysis. A value between 1.0 and 100.0 is recommended to\u2026","rel":"","context":"In "HyperMesh"","block_context":{"text":"HyperMesh","link":"https:\/\/gantovnik.com\/bio-tips\/category\/hypermesh\/"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":648,"url":"https:\/\/gantovnik.com\/bio-tips\/2020\/09\/105-how-to-remove-offsets-defined-on-shell-and-beams\/","url_meta":{"origin":666,"position":2},"title":"#105: How to remove offsets defined on shell and beams?","author":"gantovnik","date":"2020-09-25","format":false,"excerpt":"#105: How to remove offsets defined on shell and beams? Method 1: Properties > Modify > GLOBAL Method 2: Utility > Property > Property Editor Property Words to Change: 4111 Plate Offset, Only Process If Exists: Check, Action: =, Value: 0 Offsets should not be used for the beam and\u2026","rel":"","context":"In "nastran"","block_context":{"text":"nastran","link":"https:\/\/gantovnik.com\/bio-tips\/category\/nastran\/"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":651,"url":"https:\/\/gantovnik.com\/bio-tips\/2020\/09\/106-nonlinear-buckling-sol106\/","url_meta":{"origin":666,"position":3},"title":"#106: Nonlinear Buckling in MSC Nastran (SOL106)","author":"gantovnik","date":"2020-09-25","format":false,"excerpt":"#106: Nonlinear Buckling (SOL106) 1) Linear buckling of Euler column. For clamped-free boundary conditions the critical load is: Pcrit = (pi**2)*E*I\/(4*(L**2)), where, E = 10.5E6, I = 8.333-5, L=10 means Pcrit = 21.59 2) Nonlinear buckling with PARAM,BUCKLE,2 In f06 result file search for following message (right after eigenvalue table).\u2026","rel":"","context":"In "buckling"","block_context":{"text":"buckling","link":"https:\/\/gantovnik.com\/bio-tips\/category\/buckling\/"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":646,"url":"https:\/\/gantovnik.com\/bio-tips\/2020\/09\/104-how-to-generate-new-bulk-data-using-updated-design-variables-from-sol200\/","url_meta":{"origin":666,"position":4},"title":"#104: How to generate new bulk data using updated design variables from SOL200","author":"gantovnik","date":"2020-09-25","format":false,"excerpt":"#104: How to generate new bulk data using updated design variables from SOL200 ECHO = PUNCH(NEWBULK) Step 1: To extract the updated (modified) element thicknesses, use ECHO=PUNCH (NEWBULK) in the case control section. This will writes the punch file with the updated thicknesses for each element which can be used\u2026","rel":"","context":"In "nastran"","block_context":{"text":"nastran","link":"https:\/\/gantovnik.com\/bio-tips\/category\/nastran\/"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":657,"url":"https:\/\/gantovnik.com\/bio-tips\/2020\/09\/108-nastran-memory-settings\/","url_meta":{"origin":666,"position":5},"title":"#108: Nastran Memory Settings","author":"gantovnik","date":"2020-09-28","format":false,"excerpt":"#108: Nastran Memory Settings Memory is allocated on the MEM parameter. The MEM parameter cannot exceed memorymax. memorymax defaults to 50% of RAM and should not exceed 80% of physical RAM. It may be redefined on the command line or the RC file. The memory used by the solver can\u2026","rel":"","context":"In "nastran"","block_context":{"text":"nastran","link":"https:\/\/gantovnik.com\/bio-tips\/category\/nastran\/"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]}],"_links":{"self":[{"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/posts\/666","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/comments?post=666"}],"version-history":[{"count":2,"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/posts\/666\/revisions"}],"predecessor-version":[{"id":2240,"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/posts\/666\/revisions\/2240"}],"wp:attachment":[{"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/media?parent=666"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/categories?post=666"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/gantovnik.com\/bio-tips\/wp-json\/wp\/v2\/tags?post=666"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}