Home About us Contact | |||
Propagation Problems (propagation + problem)
Kinds of Propagation Problems Selected AbstractsMaterial Modelling of Porous Media for Wave Propagation ProblemsPROCEEDINGS IN APPLIED MATHEMATICS & MECHANICS, Issue 1 2003M. Schanz PD Dr.-Ing. Under the assumption of a linear geometry description and linear constitutive equations, the governing equations are derived for two poroelastic theories, Biot's theory and Theory of Porous Media (TPM), using solid displacements and pore pressure as unknowns. In both theories, this is only possible in the Laplace domain. Comparing the sets of differential equations of Biot's theory and of TPM, they show different constant coefficients but the same structure of coupled differential equations. Identifying these coefficients with the material data and correlating them leads to the known problem with Biot's ,apparent mass density'. Further, in trying to find a correlation between Biot's stress coefficient to parameters used in TPM yet unsolved inconsistencies are found. For studying the numerical effect of these differences, wave propagation results of a one-dimensional poroelastic column are analysed. Differences between both theories are resolved only for compressible constituents. [source] Extrinsic cohesive modelling of dynamic fracture and microbranching instability in brittle materialsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2007Zhengyu (Jenny) Zhang Abstract Dynamic crack microbranching processes in brittle materials are investigated by means of a computational fracture mechanics approach using the finite element method with special interface elements and a topological data structure representation. Experiments indicate presence of a limiting crack speed for dynamic crack in brittle materials as well as increasing fracture resistance with crack speed. These phenomena are numerically investigated by means of a cohesive zone model (CZM) to characterize the fracture process. A critical evaluation of intrinsic versus extrinsic CZMs is briefly presented, which highlights the necessity of adopting an extrinsic approach in the current analysis. A novel topology-based data structure is employed to enable fast and robust manipulation of evolving mesh information when extrinsic cohesive elements are inserted adaptively. Compared to intrinsic CZMs, which include an initial hardening segment in the traction,separation curve, extrinsic CZMs involve additional issues both in implementing the procedure and in interpreting simulation results. These include time discontinuity in stress history, fracture pattern dependence on time step control, and numerical energy balance. These issues are investigated in detail through a ,quasi-steady-state' crack propagation problem in polymethylmethacrylate. The simulation results compare reasonably well with experimental observations both globally and locally, and demonstrate certain advantageous features of the extrinsic CZM with respect to the intrinsic CZM. Copyright © 2007 John Wiley & Sons, Ltd. [source] Considerations of the discontinuous deformation analysis on wave propagation problemsINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 12 2009Jiong Gu Abstract In rock engineering, the damage criteria of the rock mass under dynamic loads are generally governed by the threshold values of wave amplitudes, such as the peak particle velocity and the peak particle acceleration. Therefore, the prediction of wave attenuation across fractured rock mass is important on assessing the stability and damage of rock mass under dynamic loads. This paper aims to investigate the applications of the discontinuous deformation analysis (DDA) for modeling wave propagation problems in rock mass. Parametric studies are carried out to obtain an insight into the influencing factors on the accuracy of wave propagations, in terms of the block size, the boundary condition and the incident wave frequency. The reflected and transmitted waves from the interface between two materials are also numerically simulated. To study the tensile failure induced by the reflected wave, the spalling phenomena are modeled under various loading frequencies. The numerical results show that the DDA is capable of modeling the wave propagation in jointed rock mass with a good accuracy. Copyright © 2009 John Wiley & Sons, Ltd. [source] Source signature and elastic waves in a half-space under a sustainable line-concentrated impulsive normal forceINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 4 2002Moche Ziv Abstract First, the response of an ideal elastic half-space to a line-concentrated impulsive normal load applied to its surface is obtained by a computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. Then, the geometry is determined of the obtained waves and the source signature,the latter is the imprint of the spatiotemporal configuration of the excitation source in the resultant response. Behind the dilatational precursor wave, there exists a pencil of three plane waves extending from the vertex at the impingement point of the precursor wave on the stress-free surface of the half-space to three points located on the other two boundaries of the solution domain. These four wave-arresting points (end points) of the three plane waves constitute the source signature. One wave is an inhibitor front in the behaviour of the normal stress components and the particle velocity, while in the behaviour of the shear stress component, it is a surface-axis wave. The second is a surface wave in the behaviour of the horizontal components of the dependent variables, while the third is an inhibitor wave in the behaviour of the shear stress component. An inhibitor wave is so named, since beyond it, the material motion is dying or becomes uniform. A surface-axis wave is so named, since upon its arrival, like a surface wave, the dependent variable in question features an extreme value, but unlike a surface wave, it exists in the entire depth of the solution domain. It is evident from this work that Saint-Venant's principle for wave propagation problems cannot be formulated; therefore, the above results are a consequence of the particular model proposed here for the line-concentrated normal load. Copyright © 2002 John Wiley & Sons, Ltd. [source] Algebraic preconditioning versus direct solvers for dense linear systems as arising in crack propagation problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 2 2005Erik Bängtsson Abstract Preconditioned iterative solution methods are compared with the direct Gaussian elimination method to solve dense linear systems Ax=b which originate from problems, discretized by boundary element method (BEM) techniques. Numerical experiments are presented and compared with the direct solution method available in a commercial BEM package, which show that the preconditioned iterative schemes are highly competitive with respect to both arithmetic operations required and memory demands. Copyright © 2004 John Wiley & Sons, Ltd. [source] An energy approach to space,time Galerkin BEM for wave propagation problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 9 2009A. Aimi Abstract In this paper we consider Dirichlet or Neumann wave propagation problems reformulated in terms of boundary integral equations with retarded potential. Starting from a natural energy identity, a space,time weak formulation for 1D integral problems is briefly introduced, and continuity and coerciveness properties of the related bilinear form are proved. Then, a theoretical analysis of an extension of the introduced formulation for 2D problems is proposed, pointing out the novelty with respect to existing literature results. At last, various numerical simulations will be presented and discussed, showing unconditional stability of the space,time Galerkin boundary element method applied to the energetic weak problem. Copyright © 2009 John Wiley & Sons, Ltd. [source] Modeling three-dimensional crack propagation,A comparison of crack path tracking strategiesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 9 2008P. Jäger Abstract The development of a new finite element technique for the simulation of discontinuous failure phenomena in three dimensions is the key objective of this study. In contrast to the widely used extended finite element technique, we apply a purely deformation-based strategy based on an independent interpolation of the deformation field on both sides of the discontinuity. This method has been applied successfully for two-dimensional crack propagation problems in the past. However, when it comes to three-dimensional failure phenomena, it faces the same difficulties as the extended finite element method. Unlike in two dimensions, the characterization for the three-dimensional failure surface is non-unique and the tracking of the discrete crack can be performed in several conceptually different ways. In this work, we review the four most common three-dimensional crack tracking strategies. We perform a systematic comparison in terms of standard algorithmic quality measures such as mesh independency, efficiency, robustness, stability and computational cost. Moreover, we discuss more specific issues such as crack path continuity and integratability in commercial finite element packages. The features of the suggested crack tracking algorithms will be elaborated by means of characteristic benchmark problems in failure analysis. Copyright © 2008 John Wiley & Sons, Ltd. [source] Non-reflecting artificial boundaries for transient scalar wave propagation in a two-dimensional infinite homogeneous layerINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 10 2003Chongbin Zhao Abstract This paper presents an exact non-reflecting boundary condition for dealing with transient scalar wave propagation problems in a two-dimensional infinite homogeneous layer. In order to model the complicated geometry and material properties in the near field, two vertical artificial boundaries are considered in the infinite layer so as to truncate the infinite domain into a finite domain. This treatment requires the appropriate boundary conditions, which are often referred to as the artificial boundary conditions, to be applied on the truncated boundaries. Since the infinite extension direction is different for these two truncated vertical boundaries, namely one extends toward x ,, and another extends toward x,- ,, the non-reflecting boundary condition needs to be derived on these two boundaries. Applying the variable separation method to the wave equation results in a reduction in spatial variables by one. The reduced wave equation, which is a time-dependent partial differential equation with only one spatial variable, can be further changed into a linear first-order ordinary differential equation by using both the operator splitting method and the modal radiation function concept simultaneously. As a result, the non-reflecting artificial boundary condition can be obtained by solving the ordinary differential equation whose stability is ensured. Some numerical examples have demonstrated that the non-reflecting boundary condition is of high accuracy in dealing with scalar wave propagation problems in infinite and semi-infinite media. Copyright © 2003 John Wiley & Sons, Ltd. [source] |