...
Answer: Numerical instability warnings indicate that digits of accuracy were lost during computation. CSI software performs computation with 16 digits of accuracy, some of which may be lost when analysis involves systems of significantly different stiffness, in which lower-order values are truncated to suit operations with higher-order values. This is normal for any computational software. Different levels of information are provided depending on the solver used:
- Standard solver provides details on the number of digits lost and the location of the joint and degree-of-freedom (DOF) associated with the computation
- Advanced solver indicates that a significant number of digits were lost but does not provide details on the number of digits or the location
- Multi-threaded solver provides no information on the accuracy of the computations
When this warning occurs, ETABS reports the location of the occurrence such that the user may edit the model if necessary.Reviewing animation of deformation for gravity self weight case , and first mode using the Standard solver, the number of digits lost for each equation solved is reported as follows:
- Less than 6 lost digits is not reported because the solution is sufficiently accurate.
- Between 6 and 11 lost digits is reported for informational purposes. The results of the analysis may still be acceptable but the user should carefully check the results, especially when the number of digits lost is larger.
- Over 11 digits lost is reported and also generates an instability warning. In this case the structure is either unstable or so poorly conditioned as to produce inaccurate results, and the user should check and revise the model.
The Multi-threaded solver should only be used after the stability of the model has been determined by previous analyses. If the Advanced solver indicates instabilities, re-run with the Standard solver for more details, or consider the methods below to check the model. If the Standard solver or other checks indicate sufficient stability and accuracy, then warnings from the Advanced solver can be ignored.
All three solvers report the number of negative stiffness eigenvalues found during linear equation solving ("number of eigenvalues below the shift"). This should always be zero except when performing an eigen modal analysis with a positive shift, which is not common. If the number of negative stiffness eigenvalues is not zero, a serious problem exists. This could be due to instability, very poor conditioning, or to buckling under P-delta.
A model should always be checked for stability and ill-conditioning without P-delta effects first. If the lost digits or negative stiffness eigenvalues occur only when performing P-delta analysis, please refer to the following information: P-Delta FAQ and troubleshooting.
Regardless of the solver used, it is always recommended to carefully review the deformed shape for the gravity load case and the first few eigen-modes of vibration to make sure they are as expected. This can help determine problematic areas. Some situations where numerical instability may occur include:
...
- When multiple releases are assigned to a joint, causing joint DOF to become orphaned (without stiffness). In this instance, numerical instability warnings indicate a modeling problem which should be corrected.
Orphan
- There exist orphan nodes or elements which that are not correctly attached to the rest of the structure, laterally un-braced members, or laterally unbraced other mechanisms. Reviewing deformation the deformed-shape animation under gravity cases and modal cases usually helps locate these orphan elements, which will be deflecting by large amounts and completely out-of-phase with rest of the structure
...
P-Delta FAQ and troubleshooting
Depending on the number of digits lost, actions which result include the following:
- Less than 6 does not generate a response in that values and solution are sufficiently accurate.
- Between 6 and 11 lost digits of accuracy may or may not create numerical problems, therefore warning message is reported which reads: The results of the analysis may still be acceptable but the user should carefully check the results, especially the global sums of loads.
- Over 11 digits lost causes an Ill-condition warning at the end of the analysis and results may not be sufficiently accurate. ETABS indicates the location of the numerical problem, in this case the user needs to check and revise the model.
ETABS 2013 and above, offers different options for equation solvers, when getting an ill-condition message while running Advanced or Multi-threaded solver, the user must re-run analysis using standard solver (Analyze>Advanced Sap fire), to get complete information about numerical issues, including location coordinates within the model(default coordinates are inches or mm), then check for instabilities as indicated above at these coordinates.
...
- The auto-merge tolerance may be too small, resulting in different connectivity than expected. The default is 0.1 inch or 1mm, and rarely needs to be changed
...
- . It can be reviewed using the command Options>Tolerances.
Note that Eigen modal analysis is recommended for finding instabilities and numerical problems and to better understand your structure, even though Ritz analysis is preferred for performing subsequent dynamic analyses (response-spectrum and modal time-history).
Recommendations for Large Models
In the case of very large models, running the standard solver can take an impractical amount of time , in . In this case we recommend:
- Make a copy of the model, and add unit mass to ALL all 6 DOF at every joint. Rotational inertial is necessary to pick up torsional modes for frame elements.
- Make sure the mass source is defined as from self-weigh, weight and additional masses, includes vertical mass, and is not lumped at stories
- Run an Eigen modal analysis. Any solver will do for this, so use the fastestUse the Advanced or Multi-threaded solver. A dozen modes is enough to detect problems
- Loss of accuracy more than 11 digits usually represents a rigidRigid-body mode or mechanism, and the eigen modes usually pick these up as having large periods
It is also good to look for the report on the number of eigen-values below the shift, as this may pick up problems due to P-delta that may not be near zero-frequency (large period)
The standard-solver is a good way to pick up problems like 6 to 10 digits lost due to large stiffness differences that cause inaccuracy but as indicated above this are not necessarily instabilitiesmodes or mechanisms will usually generate eigen-modes with very large periods (positive or negative). If you cannot see the deformed shape (such as for torsional column modes), generate a table of joint deflections and sort to find the largest values.
- Once the shape of the mechanism or rigid-body mode is determined, correct the connectivity or other contributing problems and check the model again before continuing with further analysis and design.