Soil–pile interaction of pile embedded in deep-layered marine sediment under seismic excitation

This paper investigates the soil–pile interaction of a pile embedded in a deep multi-layered soil under seismic excitation considering both kinematic and inertial interaction effects. A comprehensive three-dimensional finite element model is developed and validated using existing results in the literature. The response of the pile in the deep multi-layered soil profile is investigated with respect to pile head response, deflection modes and maximum deflections along the pile. Results show that the pile exhibits complex deflection patterns and that the pile response is influenced by the properties of both the soil profile and the seismic excitation. It is also evident that kinematic interaction effects have a greater influence on the pile response than the inertial interaction effects.

A strain-based failure criterion for pillar stability analysis

Strain-based failure criteria have several advantages over stress-based failure criteria: they can account for elastic and inelastic strains, they utilise direct, observables effects instead of inferred effects (strain gauges vs. stress estimates), and model complete stress-strain curves including pre-peak, non-linear elasticity and post-peak strain weakening. In this study, a strain-based failure criterion derived from thermodynamic first principles utilising the concepts of continuum damage mechanics is presented. Furthermore, implementation of this failure criterion into a finite-element simulation is demonstrated and applied to the stability of underground mining coal pillars. In numerical studies, pillar strength is usually expressed in terms of critical stresses or stress-based failure criteria where scaling with pillar width and height is common. Previous publications have employed the finite-element method for pillar stability analysis using stress-based failure criterion such as Mohr-Coulomb and Hoek-Brown or stress-based scalar damage models. A novel constitutive material model, which takes into consideration anisotropy as well as elastic strain and damage as state variables has been developed and is presented in this paper. The damage threshold and its evolution are strain-controlled, and coupling of the state variables is achieved through the damage-induced degradation of the elasticity tensor. This material model is implemented into the finite-element software ABAQUS and can be applied to 3D problems. Initial results show that this new material model is capable of describing the non-linear behaviour of geomaterials commonly observed before peak strength is reached as well as post-peak strain softening. Furthermore, it is demonstrated that the model can account for directional dependency of failure behaviour (i.e. anisotropy) and has the potential to be expanded to environmental controls like temperature or moisture.

Revisiting borehole stresses in anisotropic elastic media : comparison of analytical versus numerical solutions

We present a rigorous validation of the analyticalAmadei solution for the stress concentration around arbitrarily orientated borehole in general anisotropic elastic media. First, we revisit the theoretical framework of the Amadei solution and present analytical insights that show that the solution does indeed contain all special cases of symmetry, contrary to previous understanding, provided that the reduced strain coefficients β11 and β55 are not equal. It is shown from theoretical considerations and published experimental data that the β11 and β55 are not equal for realistic rocks. Second, we develop a 3D finite-element elastic model within a hybrid analyticalnumerical workflow that circumvents the need to rebuild and remesh the model for every borehole and material orientation. Third, we show that the borehole stresses computed from the numerical model and the analytical solution match almost perfectly for different borehole orientations (vertical, deviated and horizontal) and for several cases involving isotropic and transverse isotropic symmetries. It is concluded that the analytical Amadei solution is valid with no restrictions on the borehole orientation or elastic anisotropy symmetry.

Response of segmented bored transit tunnels to surface blast

Increasing threat of terrorism highlights the importance of enhancing the resilience of underground tunnels to all hazards. This paper develops, applies and compares the Arbitrary Lagrangian Eulerian (ALE) and Smooth Particle Hydrodynamics (SPH) techniques to treat the response of buried tunnels to surface explosions. The results and outcomes of the two techniques were compared, along with results from existing test data. The comparison shows that the ALE technique is a better method for describing the tunnel response for above ground explosion with regards to modeling accuracy and computational efficiency. The ALE technique was then applied to treat the blast response of different types of segmented bored tunnels buried in dry sand. Results indicate that the most used modern ring type segmented tunnels were more flexible for in-plane response, however, they suffered permanent drifts between the rings. Hexagonal segmented tunnels responded with negligible drifts in the longitudinal direction, but the magnitudes of in-plane drifts were large and hence hazardous for the tunnel. Interlocking segmented tunnels suffered from permanent drifts in both the longitudinal and transverse directions. Multi-surface radial joints in both the hexagonal and interlocking segments affected the flexibility of the tunnel in the transverse direction. The findings offer significant new information in the behavior of segmented bored tunnels to guide their future implementation in civil engineering applications.

A thorough investigation of the interdependent behavior of prestressed system

Due to the existing of many prestressed members in the structural system, the interdependent behavior of all prestressed members is the main concern in the analysis of the pretension process. A thorough investigation of this mutual effect is essential for an effective, reliable, and optimal analysis. Focus on this aspect, this paper presents an investigation of the interdependent behavior of all prestressed members in the whole structural system based on influence matrix (IFM). Four different types of IFM are introduced. Two different solving methods are brought forth to analyze the pretension process. The direct solving method solves for the accurate solution, whereas the iterative solving method repeatedly amends to achieve an approximate solution. A numerical example is then conducted. The result shows that various kinds of complicated batched and repeated tensioning schemes can be analyzed reliably, effectively, and completely based on IFM.

Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions

The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.

Influence of friction during forming processes-a study using a numerical simulation technique

Friction has an important influence in metal forming operations, as it contributes to the success or otherwise of the process. In the present investigation, the effect of friction on metal forming was studied by simulating compression tests on cylindrical Al-Mg alloy using the finite element method (FEM) technique. Three kinds of compression tests were considered wherein a constant coefficient of friction was employed at the upper die-work-piece interface. However, the coefficient of friction between the lower die-work-piece interfaces was varied in the tests. The simulation results showed that a difference in metal flow occurs near the interfaces owing to the differences in the coefficient of friction. It was concluded that the variations in the coefficient of friction between the dies and the work-piece directly affect the stress distribution and shape of the work-piece, having implications on the microstructure of the material being processed.

Association between aneurysm shoulder stress and abdominal aortic aneurysm expansion: A longitudinal follow-up study

Background Aneurysm expansion rate is an important indicator of the potential risk of abdominal aortic aneurysm (AAA) rupture. Stress within the AAA wall is also thought to be a trigger for its rupture. However, the association between aneurysm wall stresses and expansion of AAA is unclear. Methods and Results Forty-four patients with AAAs were included in this longitudinal follow-up study. They were assessed by serial abdominal ultrasonography and computed tomography scans if a critical size was reached or a rapid expansion occurred. Patient-specific 3-dimensional AAA geometries were reconstructed from the follow-up computed tomography images. Structural analysis was performed to calculate the wall stresses of the AAA models at both baseline and final visit. A nonlinear large-strain finite element method was used to compute the wall-stress distribution. The relationship between wall stresses and expansion rate was investigated. Slowly and rapidly expanding aneurysms had comparable baseline maximum diameters (median, 4.35 cm [interquartile range, 4.12 to 5.0 cm] versus 4.6 cm [interquartile range, 4.2 to 5.0 cm]; P=0.32). Rapidly expanding AAAs had significantly higher shoulder stresses than slowly expanding AAAs (median, 300 kPa [interquartile range, 280 to 320 kPa] versus 225 kPa [interquartile range, 211 to 249 kPa]; P=0.0001). A good correlation between shoulder stress at baseline and expansion rate was found (r=0.71; P=0.0001). Conclusion A higher shoulder stress was found to have an association with a rapidly expanding AAA. Therefore, it may be useful for estimating the expansion of AAAs and improve risk stratification of patients with AAAs.

Computed wall stress may predict the growth of abdominal aortic aneurysm

Growth rate of abdominal aortic aneurysm (AAA) is thought to be an important indicator of the potential risk of rupture. Wall stress is also thought to be a trigger for its rupture. However, stress change during the expansion of an AAA is unclear. Forty-four patients with AAAs were included in this longitudinal follow-up study. They were assessed by serial abdominal ultrasonography and computerized tomography (CT) scans if a critical size was reached or a rapid expansion occurred. Patient-specific 3-dimensional AAA geometries were reconstructed from the follow-up CT images. Structural analysis was performed to calculate the wall stresses of the AAA models at both baseline and final visit. A non-linear large-strain finite element method was used to compute the wall stress distribution. The average growth rate was 0.66cm/year (range 0-1.32 cm/year). A significantly positive correlation between shoulder tress at baseline and growth rate was found (r=0.342; p=0.02). A higher shoulder stress is associated with a rapidly expanding AAA. Therefore, it may be useful for estimating the growth expansion of AAAs and further risk stratification of patients with AAAs.

Impact of calcification and intraluminal thrombus on the computed wall stresses of abdominal aortic aneurysm

Background: Increased biomechanical stresses within the abdominal aortic aneurysm (AAA) wall contribute to its rupture. Calcification and intraluminal thrombus can be commonly found in AAAs, but the relationship between calcification/intraluminal thrombus and AAA wall stress is not completely described. Methods: Patient-specific three-dimensional AAA geometries were reconstructed from computed tomographic images of 20 patients. Structural analysis was performed to calculate the wall stresses of the 20 AAA models and their altered models when calcification or intraluminal thrombus was not considered. A nonlinear large-strain finite element method was used to compute the wall stress distribution. The relationships between wall stresses and volumes of calcification and intraluminal thrombus were sought. Results: Maximum stress was not correlated with the percentage of calcification, and was negatively correlated with the percentage of intraluminal thrombus (r = -0.56; P = .011). Exclusion of calcification from analysis led to a significant decrease in maximum stress by a median of 14% (range, 2%-27%; P < .01). When intraluminal thrombus was eliminated, maximum stress increased significantly by a median of 24% (range, 5%-43%; P < .01). Conclusion: The presence of calcification increases AAA peak wall stress, suggesting that calcification decrease the biomechanical stability of AAA. In contrast, intraluminal thrombus reduces the maximum stress in AAA. Calcification and intraluminal thrombus should both be considered in the evaluation of wall stress for risk assessment of AAA rupture.

Assessment of carotid plaque vulnerability using structural and geometrical determinants

Background Because many acute cerebral ischemic events are caused by rupture of vulnerable carotid atheroma and subsequent thrombosis, the present study used both idealized and patient-specific carotid atheromatous plaque models to evaluate the effect of structural determinants on stress distributions within plaque. Methods and Results Using a finite element method, structural analysis was performed using models derived from in vivo high-resolution magnetic resonance imaging (MRI) of carotid atheroma in 40 non-consecutive patients (20 symptomatic, 20 asymptomatic). Plaque components were modeled as hyper-elastic materials. The effects of varying fibrous cap thickness, lipid core size and lumen curvature on plaque stress distributions were examined. Lumen curvature and fibrous cap thickness were found to be major determinants of plaque stress. The size of the lipid core did not alter plaque stress significantly when the fibrous cap was relatively thick. The correlation between plaque stress and lumen curvature was significant for both symptomatic (p = 0.01; correlation coefficient: 0.689) and asymptomatic patients (p = 0.01; correlation coefficient: 0.862). Lumen curvature in plaques of symptomatic patients was significantly larger than those of asymptomatic patients (1.50±1.0mm-1 vs 1.25±0.75 mm-1; p = 0.01). Conclusion Specific plaque morphology (large lumen curvature and thin fibrous cap) is closely related to plaque vulnerability. Structural analysis using high-resolution MRI of carotid atheroma may help in detecting vulnerable atheromatous plaque and aid the risk stratification of patients with carotid disease.

Does calcium deposition play a role in the stability of atheroma? Location may be the key

Background: Rupture of vulnerable atheromatous plaque in the carotid and coronary arteries often leads to stroke and heart attack respectively. The role of calcium deposition and its contribution to plaque stability is controversial. This study uses both an idealized and a patient-specific model to evaluate the effect of a calcium deposit on the stress distribution within an atheromatous plaque. Methods: Using a finite-element method, structural analysis was performed on an idealized plaque model and the location of a calcium deposit within it was varied. In addition to the idealized model, in vivo high-resolution MR imaging was performed on 3 patients with carotid atheroma and stress distributions were generated. The individual plaques were chosen as they had calcium at varying locations with respect to the lumen and the fibrous cap. Results: The predicted maximum stress was increased by 47.5% when the calcium deposit was located in the thin fibrous cap in the model when compared with that in a model without a deposit. The result of adding a calcium deposit either to the lipid core or remote from the lumen resulted in almost no increase in maximal stress. Conclusion: Calcification at the thin fibrous cap may result in high stress concentrations, ultimately increasing the risk of plaque rupture. Assessing the location of calcification may, in the future, aid in the risk stratification of patients with carotid stenosis.

Instability of laminated composite thin-walled open-section beams

Instability of thin-walled open-section laminated composite beams is studied using the finite element method. A two-noded, 8 df per node thin-walled open-section laminated composite beam finite element has been used. The displacements of the element reference axis are expressed in terms of one-dimensional first order Hermite interpolation polynomials, and line member assumptions are invoked in formulation of the elastic stiffness matrix and geometric stiffness matrix. The nonlinear expressions for the strains occurring in thin-walled open-section beams, when subjected to axial, flexural and torsional loads, are incorporated in a general instability analysis. Several problems for which continuum solutions (exact/approximate) are possible have been solved in order to evaluate the performance of finite element. Next its applicability is demonstrated by predicting the buckling loads for the following problems of laminated composites: (i) two layer (45°/−45°) composite Z section cantilever beam and (ii) three layer (0°/45°/0°) composite Z section cantilever beam.

Monoblock cylinder internally pressurised under asymmetric end constraints Some aspects of design

The finite element method is used to analyse stresses and displacements in a monoblock cylinder open at one end only. The cylinder is internally pressurised. The analysis shows that the minimum pressure required to cause yield in the cylinder decreases rapidly with increasing cylinder height until the height is about the same as the outer radius of the cylinder, beyond which the decrease is marginal. Introduction of a fillet at the internal corner enhances the design pressure substantially while a fillet at the outer corner affects this pressure only marginally.