Products: ABAQUS/Standard ABAQUS/Explicit
A cylindrical shell panel, firmly clamped on all four sides, is exposed to the detonation of a high explosive layer. The problem illustrates the use of initial velocity conditions to model sudden, impulsive loadings arising from the detonation. In the course of the analysis a strong plastic hinge forms along the edge of the detonation area. Both experimental and numerical results for this problem have been reported by Leech (1966) and Morino et al. (1971).
The panel is 319 mm (12.56 in) long and spans a 120° sector of a cylinder, with a midsurface radius of 74.6 mm (2.938 in) and a thickness of 3.18 mm (0.125 in). Only 60° of the panel is modeled because of the symmetry of the problem. Clamped boundary conditions are prescribed on three edges of the model, while the appropriate symmetry conditions are imposed along the remaining edge.
The shell is made from 6061-T6 aluminum alloy with a Young's modulus of 72.4 GPa (10.5 × 106 psi), a Poisson's ratio of 0.33, and a density of 2672 kg/m3 (2.5 × 104 lb sec2 in4). A von Mises elastic, perfectly plastic material model is used with a yield stress of 303 MPa (4.4 × 104 psi).
In the experiment the high explosive layer covers a 60° sector of the panel, extending 259 mm (10.21 in) from one end. Hence, there is no symmetry plane along the y-axis. All nodes in contact with the high explosive layer have been grouped in a node set named BLAST. The effect of the detonation is simulated by prescribing an initial inward radial velocity of 144 m/sec (5650 in/sec) to the nodes in this set.
All the relevant shell element types available in ABAQUS/Standard are used in the simulation for comparative purposes and as a gauge of the relative merits of each element type for this class of problem. An 8 × 16 mesh is used for first-order elements, and a 4 × 8 mesh is used for second-order elements.
The ABAQUS/Explicit analysis is performed using the finite-strain element, S4R, for three different mesh refinements (8 × 32, 16 × 32, and 32 × 64) and the small-strain elements, S4RS and S4RSW, for a 32 × 64 mesh. Geometrically equivalent analyses employing a shell offset with a value of 0.5 are performed using each of the quadrilateral shell elements in ABAQUS/Explicit for a 32 × 64 mesh refinement. In addition, an analysis is performed with a 16 × 32 mesh of S4R elements using ENHANCED hourglass control.
For the ABAQUS/Standard analysis we choose to set the time integration accuracy control parameter (HAFTOL) to a very large (essentially infinite) value. This implies that we are choosing automatic control for the time stepping, but we are not controlling the accuracy of the time integration. The time increments will be limited only by the ability of the Newton scheme to solve the nonlinear equilibrium equations. This is a common technique for obtaining low-cost solutions for highly dissipative, strongly nonlinear cases. It is effective because the nonlinearities limit the time increments, and the high level of dissipation quickly removes the high frequency content from the solution. In practice it is desirable to verify the results with a second, more expensive, analysis in which a realistic value of HAFTOL is used. Default controls are used in ABAQUS/Explicit.
In both the experimental results and the ABAQUS simulations, peak deflection occurs after about 400 s. Figure 1.3.31 shows deformed configuration plots for the S4R5 model and the S9R5 model after 400 s of response time. Figure 1.3.32 shows the deformed shapes at 400 s for the three meshes used in the ABAQUS/Explicit analysis.
The calculated values for the maximum deflection at a point midway along the centerline of the panel are reported for each of the analysis cases in Table 1.3.31. The experimental result for the maximum deflection reported by Morino et al. (1971) is also included for comparison. The mode of deformation in the problem is predominantly bending, and the second-order element models outperform the first-order element models for similar-cost analyses in ABAQUS/Standard. These meshes are quite coarse, and improved performance is observed in ABAQUS/Explicit upon mesh refinement. The results suggest that the 16 × 32 mesh of first-order elements provides a reasonably accurate solution for the maximum deflection. In addition, the results obtained using ENHANCED hourglass control closely match those obtained using the default hourglass control formulation.
S3R shell model.
S4 shell model.
S4R shell model.
S4R5 shell model.
S8R shell model.
S8R5 shell model.
S9R5 shell model.
STRI65 shell model.
S4R elements, fine mesh case.
S4R elements, 8 × 32 mesh.
S4R elements, 16 × 32 mesh.
S4R elements, 16 × 32 mesh, ENHANCED hourglass control.
S4RS elements.
S4RSW elements.
S4R analysis, shell offset.
S4RS analysis, shell offset.
S4RSW analysis, shell offset.
Additional high mesh refinement case included for the sole purpose of testing the performance of the code.
Leech, J. W., Finite-Difference Calculation Method for Large Elastic-Plastic Dynamically-Induced Deformations of General Thin Shells, Ph.D. Thesis, Dept. of Aeronautics and Astronautics, Massachussetts Institute of Technology, Cambridge, MA, 1966.
Morino, L., J. W. Leech, and E. A. Witmer, An Improved Numerical Calculation Technique for Large Elastic-Plastic Transient Deformations of Thin Shells: Part 2—Evaluation and Applications, Journal of Applied Mechanics, vol. 38, pp. 429436, 1971.
Table 1.3.31 Maximum deflection along center line of the panel at y=159.5 mm (6.28 in).
Code | Element Type | Mesh Size | POISSON | OFFSET | Maximum Deflection | |
---|---|---|---|---|---|---|
mm | in | |||||
ABAQUS/Standard | S4 | 8 × 16 | 0 | 0 | 28.70 | 1.13 |
ABAQUS/Standard | S4 | 8 × 16 | 0.5 | 0 | 27.43 | 1.08 |
ABAQUS/Standard | S4R | 8 × 16 | 0 | 0 | 28.96 | 1.14 |
ABAQUS/Standard | S4R | 8 × 16 | 0.5 | 0 | 27.69 | 1.09 |
ABAQUS/Standard | S4R5 | 8 × 16 | 0.5 | 0 | 30.48 | 1.20 |
ABAQUS/Standard | S3R | 8 × 16 | 0.5 | 0 | 31.24 | 1.23 |
ABAQUS/Standard | S8R | 4 × 8 | 0.5 | 0 | 31.24 | 1.23 |
ABAQUS/Standard | S8R5 | 4 × 8 | 0.5 | 0 | 31.24 | 1.23 |
ABAQUS/Standard | S9R5 | 4 × 8 | 0.5 | 0 | 31.49 | 1.24 |
ABAQUS/Standard | STRI65 | 4 × 8 | 0.5 | 0 | 31.49 | 1.24 |
ABAQUS/Explicit | S4R | 8 × 32 | — | 0 | 26.24 | 1.033 |
ABAQUS/Explicit | S4R | 16 × 32 | — | 0 | 30.81 | 1.213 |
ABAQUS/Explicit | S4R(enhanced hourglass) | 16 × 32 | — | 0 | 30.71 | 1.209 |
ABAQUS/Explicit | S4R | 32 × 64 | — | 0 | 30.86 | 1.215 |
ABAQUS/Explicit | S4R | 32 × 64 | — | 0.5 | 29.84 | 1.175 |
ABAQUS/Explicit | S4RS | 32 × 64 | — | 0 | 31.11 | 1.225 |
ABAQUS/Explicit | S4RS | 32 × 64 | — | 0.5 | 29.19 | 1.149 |
ABAQUS/Explicit | S4RSW | 32 × 64 | — | 0 | 31.21 | 1.229 |
ABAQUS/Explicit | S4RSW | 32 × 64 | — | 0.5 | 29.29 | 1.153 |
Experimental | — | — | — | — | 31.75 | 1.250 |