scholarly journals Evaluation of Dynamic Soil-Pile-Structure Interactive Behavior in Dry Sand by 3D Numerical Simulation

2019 ◽  
Vol 9 (13) ◽  
pp. 2612 ◽  
Author(s):  
Sun Yong Kwon ◽  
Mintaek Yoo

A 3D numerical model based on finite-difference approximation was formulated to predict the dynamic soil-pile-structure interaction (SPSI) in dry sand. A non-linear elastic, Mohr–Coulomb plastic soil-constitutive model was adopted for the proposed methodology with a hysteretic damping model which can simulate nonlinear behavior of soil and an interface model which can predict separation and slippage between soil and pile according to the external load condition. Simplified continuum model was used to properly simulate the semi-infinite boundary and improve analysis efficiency. The proposed numerical model was validated by comparison with experimental results performed by Yoo (2013). Thereafter, a parametric study was also carried out to investigate the complex dynamic behavior of pile foundation under varying conditions. It was demonstrated that inertial force induced by superstructure is dominant for dynamic SPSI in dry sand whereas the kinematic force induced by soil deformation is relatively insignificant. Pile peak bending moment occurs at 30% of the pile length when pile length is no longer than 5 T and at about 30% of 5 T (1.6 T) when the pile length is longer than 5 T. The pile head fixity governed the peak bending moment profile of pile and affected the dynamic responses of the system in conjunction with other factors, such as pile rigidity.

1982 ◽  
Vol 28 (99) ◽  
pp. 239-265 ◽  
Author(s):  
Robert Bindschadler

AbstractA time-dependent numerical model of temperate glacier flow without sliding is developed and applied to the quiescent phase of surge-type Variegated Glacier, Alaska. The model is based on a one-dimensional continuity equation but the transverse channel shape is explicitly included allowing the complex geometries of real glaciers to be modelled. Velocities and volume fluxes are calculated from the glacier geometry. Transverse stress is taken into account by shape factors which are fitted to measurements of geometry and velocity and are chosen to be insensitive to changes in geometry. Longitudinal stress gradients are taken into account by use of a large-scale surface slope. A Crank-Nicholson finite-difference approximation is used and it is unconditionally stable when a small contribution from the local slope is added to the average slope.Model parameters are fitted to extensive data collected on Variegated Glacier in 1973 and 1974. Predictions of the model over a four year interval agree well with field measurements. Predictions of the current quiescent phase (1965–84) indicate depth increases in the upper glacier of more than 75 m with a twenty-fold increase in the volume flux. During this interval the base shear stress increases 40% in the upper glacier and decreases 20% in the lower glacier. During the mid to late quiescent phase, ice motion becomes more important than mass balance in the redistribution of mass over the central region of the glacier. If normal flow were to persist, the predicted steady-state profile would be an average of 100 m deeper and 41% more voluminous than in 1973.The predicted base shear-stress gradient is never negative enough to satisfy Robin and Weertman’s (1973) condition for blockage of subglacial water flow. The annual rate of water production by dissipation of mechanical straining at the bed remains two orders of magnitude below that produced by summer surface melt. The predicted fractional increase in base stress during the quiescent phase is a maximum in the region believed to be the trigger zone of the surges.


1982 ◽  
Vol 28 (99) ◽  
pp. 239-265 ◽  
Author(s):  
Robert Bindschadler

AbstractA time-dependent numerical model of temperate glacier flow without sliding is developed and applied to the quiescent phase of surge-type Variegated Glacier, Alaska. The model is based on a one-dimensional continuity equation but the transverse channel shape is explicitly included allowing the complex geometries of real glaciers to be modelled. Velocities and volume fluxes are calculated from the glacier geometry. Transverse stress is taken into account by shape factors which are fitted to measurements of geometry and velocity and are chosen to be insensitive to changes in geometry. Longitudinal stress gradients are taken into account by use of a large-scale surface slope. A Crank-Nicholson finite-difference approximation is used and it is unconditionally stable when a small contribution from the local slope is added to the average slope.Model parameters are fitted to extensive data collected on Variegated Glacier in 1973 and 1974. Predictions of the model over a four year interval agree well with field measurements. Predictions of the current quiescent phase (1965–84) indicate depth increases in the upper glacier of more than 75 m with a twenty-fold increase in the volume flux. During this interval the base shear stress increases 40% in the upper glacier and decreases 20% in the lower glacier. During the mid to late quiescent phase, ice motion becomes more important than mass balance in the redistribution of mass over the central region of the glacier. If normal flow were to persist, the predicted steady-state profile would be an average of 100 m deeper and 41% more voluminous than in 1973.The predicted base shear-stress gradient is never negative enough to satisfy Robin and Weertman’s (1973) condition for blockage of subglacial water flow. The annual rate of water production by dissipation of mechanical straining at the bed remains two orders of magnitude below that produced by summer surface melt. The predicted fractional increase in base stress during the quiescent phase is a maximum in the region believed to be the trigger zone of the surges.


Materials ◽  
2021 ◽  
Vol 14 (10) ◽  
pp. 2589
Author(s):  
Jung J. Kim

This study presents an explosion-resistant hybrid system containing a steel slab and a carbon fiber-reinforced polymer (CFRP) frame. CFRP, which is a high-strength material, acts as an impact reflection part. Steel slab, which is a high-ductility material, plays a role as an impact energy absorption part. Based on the elastoplastic behavior of steel, a numerical model is proposed to simulate the dynamic responses of the hybrid system under the air pressure from an explosion. Based on this, a case study is conducted to analyze and identify the optimal design of the proposed hybrid system, which is subjected to an impact load condition. The observations from the case study show the optimal thicknesses of 8.2 and 7 mm for a steel slab and a ϕ100 mm CFRP pipe for the hybrid system, respectively. In addition, the ability of the proposed hybrid system to resist an uncertain explosion is demonstrated in the case study based on the reliability methodology.


2012 ◽  
Vol 12 (1) ◽  
pp. 193-225 ◽  
Author(s):  
N. Anders Petersson ◽  
Björn Sjögreen

AbstractWe develop a stable finite difference approximation of the three-dimensional viscoelastic wave equation. The material model is a super-imposition of N standard linear solid mechanisms, which commonly is used in seismology to model a material with constant quality factor Q. The proposed scheme discretizes the governing equations in second order displacement formulation using 3N memory variables, making it significantly more memory efficient than the commonly used first order velocity-stress formulation. The new scheme is a generalization of our energy conserving finite difference scheme for the elastic wave equation in second order formulation [SIAM J. Numer. Anal., 45 (2007), pp. 1902-1936]. Our main result is a proof that the proposed discretization is energy stable, even in the case of variable material properties. The proof relies on the summation-by-parts property of the discretization. The new scheme is implemented with grid refinement with hanging nodes on the interface. Numerical experiments verify the accuracy and stability of the new scheme. Semi-analytical solutions for a half-space problem and the LOH.3 layer over half-space problem are used to demonstrate how the number of viscoelastic mechanisms and the grid resolution influence the accuracy. We find that three standard linear solid mechanisms usually are sufficient to make the modeling error smaller than the discretization error.


Author(s):  
Yawen Wang ◽  
Junyi Yang ◽  
Dong Guo ◽  
Teik C Lim

A generalized dynamic model of driveline system is formulated that includes the coupling effect and gyroscopic moments of the propeller shaft and hypoid gear rotor assembly. Firstly, the dynamic models with only gear-shaft coupling, with only gyroscopic effect, and with both gear-shaft coupling and gyroscopic effect are analyzed and compared. The results show that the combined effects of the gear-shaft interaction and gyroscopic behavior have considerable influence on the system dynamic responses surrounding gear bending resonances, especially for the bearing responses. However, the gear out-of-phase torsional modes still dominate the gear mesh frequency response. Secondly, the influence of pinion bending moment of inertia, propeller shaft stiffness and bearing stiffness on the system dynamic responses are examined. The system responses are then applied to perform further vibration and acoustic analysis for an axle housing structure. Computational results reveal that NVH (noise, vibration, and harshness) refinement can be achieved by tuning the pinion bearing rotational stiffness and pinion bending moment of inertia for the example considered. This study provides an understanding of the interaction between hypoid gear pair and propeller shaft, and can be employed to enhance driveline system design.


1966 ◽  
Vol 6 (03) ◽  
pp. 217-227 ◽  
Author(s):  
Hubert J. Morel-Seytoux

Abstract The influence of pattern geometry on assisted oil recovery for a particular displacement mechanism is the object of investigation in this paper. The displacement is assumed to be of unit mobility ratio and piston-like. Fluids are assumed incompressible and gravity and capillary effects are neglected. With these assumptions it is possible to calculate by analytical methods the quantities of interest to the reservoir engineer for a great variety of patterns. Specifically, this paper presentsvery briefly, the methods and mathematical derivations required to obtain the results of engineering concern, andtypical results in the form of graphs or formulae that can be used readily without prior study of the methods. Results of this work provide checks for solutions obtained from programmed numerical techniques. They also reveal the effect of pattern geometry and, even though the assumptions of piston-like displacement and of unit mobility ratio are restrictive, they can nevertheless be used for rather crude but quick, cheap estimates. These estimates can be refined to account for non-unit mobility ratio and two-phase flow by correlating analytical results in the case M=1 and the numerical results for non-Piston, non-unit mobility ratio displacements. In an earlier paper1 it was also shown that from the knowledge of closed form solutions for unit mobility ratio, quantities called "scale factors" could be readily calculated, increasing considerably the flexibility of the numerical techniques. Many new closed form solutions are given in this paper. INTRODUCTION BACKGROUND Pattern geometry is a major factor in making water-flood recovery predictions. For this reason many numerical schemes have been devised to predict oil recovery in either regular patterns or arbitrary configurations. The numerical solutions, based on the method of finite difference approximation, are subject to errors often difficult to evaluate. An estimate of the error is possible by comparison with exact solutions. To provide a variety of checks on numerical solutions, a thorough study of the unit mobility ratio displacement process was undertaken. To calculate quantities of interest to the reservoir engineer (oil recovery, sweep efficiency, etc.), it is necessary to first know the pressure distribution in the pattern. Then analytical procedures are used to calculate, in order of increasing difficulty: injectivity, breakthrough areal sweep efficiency, normalized oil recovery and water-oil ratio as a function of normalized PV injected. BACKGROUND Pattern geometry is a major factor in making water-flood recovery predictions. For this reason many numerical schemes have been devised to predict oil recovery in either regular patterns or arbitrary configurations. The numerical solutions, based on the method of finite difference approximation, are subject to errors often difficult to evaluate. An estimate of the error is possible by comparison with exact solutions. To provide a variety of checks on numerical solutions, a thorough study of the unit mobility ratio displacement process was undertaken. To calculate quantities of interest to the reservoir engineer (oil recovery, sweep efficiency, etc.), it is necessary to first know the pressure distribution in the pattern. Then analytical procedures are used to calculate, in order of increasing difficulty: injectivity, breakthrough areal sweep efficiency, normalized oil recovery and water-oil ratio as a function of normalized PV injected.


Sign in / Sign up

Export Citation Format

Share Document