In this paper, the high-frequency response of isotropic-laminated cylindrical shells is investigated using a layer-wise theory. The cylindrical shell is discretized in an arbitrary number of layers in the radial direction, and a three-dimensional stress state is assumed in each layer. Approximate numerical results obtained by the layer-wise theory are compared with the exact wave-dispersion analytical results. The very good agreement between approximate and exact results indicates that the layer-wise theory can accurately describe of the dynamic response of cylindrical shells in the high-frequency (short-wavelength) range.

This article presents the first effort to develop a two-dimensional model using the extended finite element method (XFEM) for the simulation of discrete fracture networks, in which the mesh does not conform to the natural fracture network. The model incorporates contact, cohesion, and friction between blocks of rock. Shear dilation is an important mechanism impacting the overall nonlinear response of naturally fractured rock masses and is also included in the model; physics previously not simulated within an XFEM context. Here, shear dilation is modeled by means of a linear dilation model, capped by a dilation limiting displacement. Highly nonlinear problems involving multiple joint sets are investigated within a quasi-static context. An explicit scheme is used in conjunction with the dynamic relaxation technique to obtain equilibrium solutions in the face of the nonlinear constitutive models from contact, cohesion, friction, and dilation. The numerical implementation is verified and its convergence is illustrated using a shear test and a biaxial test. The model is then applied to the practical problem of the stability of a slope of fractured rock. KEYWORDSextended finite element method, discrete fracture network, shear dilation Int

No abstract

This article presents a numerical methodology for the simulation of mineral dissolution which couples brine flow, dissolved mineral transport, and cavity evolution. The proposed model considers both (1) the varying density brine flow within the cavity governed by the compressible Navier‐Stokes equations and (2) the evolution of the cavity boundary using a sharp interface model with a physically‐derived dissolution rate equation. The proposed nonlinear multi‐physics model can capture complex flow patterns such as the generation of a vortex in the cavity. The impact of those complex flow patterns on the cavity development can be studied because of the coupling of brine flow and dissolution front movement. The model employs a new strategy to explicitly track the dissolution front, which results in low computational cost for long‐term dissolution simulations. The proposed model is verified through a convergence analysis, showing both spatial and temporal convergence. Numerical simulations of mineral dissolution in horizontal cavities are conducted to investigate the flow velocity, mass fraction of dissolved mineral, cavity shape evolution, and dissolution rate over time. Additionally, a discussion on the effect of Peclet number on mineral dissolution in the cavity is undertaken.

No abstract