320 likes | 680 Views
Comparative Analysis of Layered Elastic Computer Programs – LEAF, JULEA, and BISAR. Gordon F. Hayhoe FAA AAR-410. Need for Verifying the Accuracy of LEA Programs. LEA programs work well for simple structures. LEA programs can give very low accuracy solutions under certain conditions.
E N D
Comparative Analysis of Layered Elastic Computer Programs – LEAF, JULEA, and BISAR Gordon F. Hayhoe FAA AAR-410 1
Need for Verifying the Accuracy of LEA Programs • LEA programs work well for simple structures. • LEA programs can give very low accuracy solutions under certain conditions. • Very few closed form solutions which can be used to provide numerical checks. • Compare results from multiple programs and select the “correct” answers by consensus and reasonableness. 2
Layered Elastic Equations • is vertical stress, for example, where: • equation is inverse Hankel transform. • q = tire pressure • a = tire radius • r = evaluation point radius • z = depth from top of layer • a = Hankel domain variable a r z h i 4
Evaluation Points Close to the Top of the Structure • Load is constant pressure and discontinuous on the surface at the edge of the contact area. • Therefore, stresses are discontinuous on the surface. • Computed values close to the surface are unreliable because discontinuous functions cannot be properly represented by discrete transforms (using numerical integration). 6
Different Methods • BISAR changes to a special set of infinite series equations valid when the evaluation point is close to the surface. • JULEA aborts stress calculation when the evaluation point is close to the surface. JULEA extrapolates for vertical deflection (a special case). • LEAF shifts the origin vertically – improves the accuracy but does not eliminate the problem completely for stresses. 7
Different Methods • BISAR has two equation solvers, user selected. One is for bonded interfaces and the other is for unbonded interfaces. • LEAF has three equation solvers, automatically selected (on magnitude of residuals). • Part inverting – fastest, but inaccurate with ill-conditioned equations. • Partial pivoting (LU) – intermediate. • Full pivoting (Gauss-Jordan) – slow, but appears to give satisfactory accuracy under all critical conditions identified and tested. 8
Comparison with Boussinesq for a Uniform Half-Space zeval = distance from surface 9
Comparison with Boussinesq for a Uniform Half-Space zeval = distance from surface 10
Horizontal Stress Variation with Depth BISAR output precision is three significant figures in scientific format. 14
Widely Spaced Evaluation Points • LEDFAA 1.3 computes the vertical subgrade strain in flexible structures with the contributions from all wheels in the main gear of the aircraft. • This is 16 wheels for the B-747 and 20 wheels for the A380. • Need to be sure that the contributions of the wheels far away from the evaluation point are accurate and reliable. 17
Widely Spaced Evaluation Points • The contributions of the wheels far away are very small and a cut off distance could be used in many cases. • But an appropriate cut off radius would vary with structure type, thickness, subgrade strength, etc. • Therefore, use a consistent method. 19
Widely Spaced Evaluation Points • Flexible structure. • 5 in surface, 8 in stabilized base, 11 in subbase, 10 CBR subgrade. • 50,000 lb single-wheel load. • Seven evaluation points, 100 in intervals. 20
Conclusions • It is important to check the accuracy of LEA programs for extreme conditions. • The latest version of LEAF, used in LEDFAA and BAKFAA, provides accurate and reliable results for: • FWD backcalculation. • Unbonded concrete overlays. • Multiple-gear aircraft having widely spaced main gear wheels. 27