VM289

VM289
Pressurized Spherical Cavity in an Infinite Medium with Perfectly Matched Layers (PML)

Overview

Reference: Vurrujit, A., Soil Dynamics, Delft University of Technology, pp. 138-139
Analysis Type(s):
Static Analysis (ANTYPE = 0)
Element Type(s):
3D 20-Node Structural Solid (SOLID186)
Input Listing: vm289.dat

Test Case

A spherical cavity with a radius of 1.0 m is placed in an infinite medium. The spherical cavity is subjected to an internal pressure of 750 Pa. A perfectly matched layer (PML) is used to truncate the infinite boundary. Determine the displacement and stresses with respect to location (from X = 1 m to X = 4 m).

Figure 504: Finite Element/PML Mesh of a Pressurized Spherical Cavity in an Infinite Medium

Finite Element/PML Mesh of a Pressurized Spherical Cavity in an Infinite Medium

Material PropertiesGeometric PropertiesLoading
Young’s modulus, E = 107 Pa
Poisson’s ratio, ν = 0.33
Internal radius, a = 1 mInternal pressure = 750 Pa

Analysis Assumptions and Modeling Notes

The spherical cavity in the infinite medium and perfectly matched layers are modeled using SOLID186 elements. The PML absorbing condition is defined by setting Keyopt(15) = 1. Only one eighth of the sphere is modeled, using symmetric constraints on the three orthagonal planes. The displacements along the radial direction and the stresses along the radial and tangential directions are then obtained from /POST1 and compared against the analytical solution.

The analytic solution for radial displacement is:

The analytic solution for radial stress is:

The analytic solution for tangential stress is:

where is the internal applied pressure, is the internal radius, is the rigidity modulus, and is the radial location.

Results Comparison

 Location (m)TargetMechanical APDLRatio
Radial displacement (m)1.00000.49875 x 10-4 0.49852 x 10-4 1.000
1.94740.13151 x 10-4 0.13137 x 10-4 0.999
3.05260.53523 x 10-5 0.53338 x 10-5 0.997
3.84210.33787 x 10-5 0.33562 x 10-5 0.993
Radial stress (Pa)1.0000-750.000-736.6620.982
1.9474-101.554-101.4420.999
3.0526-26.366-26.5891.008
3.8421-13.224-13.4371.016
Tangential stress (Pa)1.0000375.000368.6290.983
1.947450.77750.4900.994
3.052613.18313.0600.990
3.84216.6126.4830.980