An extended displacement discontinuity method (EDDM) is proposed to analyze the stress wave propagation in jointed viscoelastic rock mass (VRM).The discontinuities in a rock mass are divided into two groups.The primar...An extended displacement discontinuity method (EDDM) is proposed to analyze the stress wave propagation in jointed viscoelastic rock mass (VRM).The discontinuities in a rock mass are divided into two groups.The primary group with an average geometrical size larger than or in the same order of magnitude of wavelength of a concerned stress wave is defined as 'macro-joints',while the secondary group with a high density and relatively small geometrical size compared to the wavelength is known as 'micro-defects'.The rock mass with micro-defects is modeled as an equivalent viscoelastic medium while the macro-joints in the rock mass are modeled explicitly as physical discontinuities.Viscoelastic properties of a micro-defected sedimentary rock are obtained by longitudinally impacting a cored long sedimentary rod with a pendulum.Wave propagation coefficient and dynamic viscoelastic modulus are measured.The EDDM is then successfully employed to analyze the wave propagation across macro-joint in VRM.The effect of the rock viscosity on the stress wave propagation is evaluated by comparing the results of VRM from the presented EDDM with those of an elastic rock mass (ERM) from the conventional displacement discontinuity method (CDDM).The CDDM is a special case of the EDDM under the condition that the rock viscosity is ignored.Comparison of the reflected and transmitted waves shows that the essential rock viscosity has a significant effect on stress wave attenuation.When a short propagation distance of a stress wave is considered,the results obtained from the CDDM approximate to the EDDM solutions,however,when the propagation distance is sufficiently long relative to the wavelength,the effect of rock viscosity on the stress wave propagation cannot be ignored.展开更多
An unstructured finite volume time domain method (UFVTDM) is proposed to simulate stress wave propagation. The original variables of displacement and stress are solved based on the dynamic equilibrium equations. An Eu...An unstructured finite volume time domain method (UFVTDM) is proposed to simulate stress wave propagation. The original variables of displacement and stress are solved based on the dynamic equilibrium equations. An Euler explicit and unstructured finite volume method is used for time and spacial terms respectively. The displacements are stored on the cell vertex and a vertex based finite volume is formed with the integral surface and the stress is assumed as uniform in the cell. This is some similar with the stager grid method in computational fluid dynamics. Several cases are used to show the capability of the algorithm.展开更多
An experimental-numerical method for measuring dynamic crack propagatingvelocities under stress wave loading is established in this paper. The experiments of thethree-point bend specimen are done on the improved Hopki...An experimental-numerical method for measuring dynamic crack propagatingvelocities under stress wave loading is established in this paper. The experiments of thethree-point bend specimen are done on the improved Hopkinson bar. Deflection of loading point,dynamic load and instantaneous crack length are measured, then crack propagating velocities arecalculated. Experiments on 40Cr steel show that the results given by this method have a goodagreement with that obtained by the resistance fracture gage method. Therefore this method isfeasible for measuring crack propagating velocities under high loading rate and will have wideapplication.展开更多
Using the complex variable method and conformal mapping,scat- tering of flexural waves and dynamic stress concentrations in Mindlin's thick plates with a cutout have been studied.The general solution of the stress...Using the complex variable method and conformal mapping,scat- tering of flexural waves and dynamic stress concentrations in Mindlin's thick plates with a cutout have been studied.The general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of cutouts is obtained. Applying the orthogonal function expansion technique,the dynamic stress problem can be reduced into the solution of a set of infinite algebraic equations.As examples, numerical results for the dynamic stress concentration factor in Mindlin's plates with a circular,elliptic cutout are graphically presented in sequence.展开更多
By using the complex variable method and conformal mapping, the diffraction of flexural waves and dynamic stress concentrations in thick plates with a cavity have been studied. A general solution of the stress problem...By using the complex variable method and conformal mapping, the diffraction of flexural waves and dynamic stress concentrations in thick plates with a cavity have been studied. A general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of an arbitrary cavity is obtained. By employing the orthogonal function expansion technique, the dynamic stress problem can be reduced to the solution of an infinite algebraic equation series. As an example, the numerical results for the dynamic stress concentration factor in thick plates with a circular, elliptic cavity are graphically presented. The numerical results are discussed.展开更多
Boundary Element Method (BEM) is employed to run theoretical analsis and numerical calculation of dif-fraction of elastic wave and dynamic stress concentration in an infinite then plate with a cireular hole. Based on ...Boundary Element Method (BEM) is employed to run theoretical analsis and numerical calculation of dif-fraction of elastic wave and dynamic stress concentration in an infinite then plate with a cireular hole. Based on the work equivalent law of dynamics,boundary integral equation is established for flexural waves of thin plate. Calculation formulas of influence coefficients are derived using Mathematica software and numerical results are obtained for dynam-ic stress conoentration factors in a then plate with a circular hole.展开更多
In this paper, based on complex variables and conformal mapping methods, using the refined dynamic equation of plates, elastic wave scattering and dynamic stress concentrations in plates with two cutouts were studied....In this paper, based on complex variables and conformal mapping methods, using the refined dynamic equation of plates, elastic wave scattering and dynamic stress concentrations in plates with two cutouts were studied. Applying the orthogonal function expansion method, the problem to be solved can be reduced into the solution of a set of infinite algebraic equations. According to free boundary conditions, numerical results of dynamic moment concentration factors in thick plates with two circular cutouts analyze that: there will be more complex interaction changes between two-cutout situation than single cutout situation. In the case of low frequency or high frequency and thin plate, the hole-spacing in the absence of coupling interactions was larger or smaller. The numerical results and method can be used to analyze the dynamics and strength of plate-like structures.展开更多
A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)t...A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.展开更多
The indirect boundary element method （IBEM） is developed to solve the scattering of plane SH-waves by a lined tunnel in elastic wedge space. According to the theory of single-layer potential, the scattered-wave fiel...The indirect boundary element method （IBEM） is developed to solve the scattering of plane SH-waves by a lined tunnel in elastic wedge space. According to the theory of single-layer potential, the scattered-wave field can be constructed by applying virtual uniform loads on the surface of lined tunnel and the nearby wedge surface. The densities of virtual loads can be solved by establishing equations through the continuity conditions on the interface and zero-traction conditions on free surfaces. The total wave field is obtained by the superposition of free field and scattered-wave field in elastic wedge space. Numerical results indicate that the IBEM can solve the diffraction of elastic wave in elastic wedge space accurately and effi- ciently. The wave motion feature strongly depends on the wedge angle, the angle of incidence, incident frequency, the location of lined tunnel, and material parameters. The waves interference and amplification effect around the tunnel in wedge space is more significant, causing the dynamic stress concentration factor on rigid tunnel and the displacement amplitude of flexible tunnel up to 50.0 and 17.0, respectively, more than double that of the case of half-space. Hence, considerable attention should be paid to seismic resistant or anti-explosion design of the tunnel built on a slope or hillside.展开更多
A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks.Numerical results are presented for the dynamic stress in- tensit...A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks.Numerical results are presented for the dynamic stress in- tensity factors.The effects of the wave type,wave frequency,wave incidence angle,and crack spacing on the dynamic stress intensity factors are analyzed in detail.展开更多
A virtual Taylor impact of cellular materials is analyzed with a wave propagation technique, i.e. the Lagrangian analysis method, of which the main advantage is that no pre-assumed constitutive relationship is require...A virtual Taylor impact of cellular materials is analyzed with a wave propagation technique, i.e. the Lagrangian analysis method, of which the main advantage is that no pre-assumed constitutive relationship is required. Time histories of particle velocity, local strain, and stress profiles are calculated to present the local stress-strain history curves, from which the dynamic stress-strain states are obtained. The present results reveal that the dynamic-rigid-plastic hardening （D-R-PH） material model introduced in a previous study of our group is in good agreement with the dynamic stress-strain states under high loading rates obtained by the Lagrangian analysis method. It directly reflects the effectiveness and feasibility of the D-R-PH material model for the cellular materials under high loading rates.展开更多
文摘An extended displacement discontinuity method (EDDM) is proposed to analyze the stress wave propagation in jointed viscoelastic rock mass (VRM).The discontinuities in a rock mass are divided into two groups.The primary group with an average geometrical size larger than or in the same order of magnitude of wavelength of a concerned stress wave is defined as 'macro-joints',while the secondary group with a high density and relatively small geometrical size compared to the wavelength is known as 'micro-defects'.The rock mass with micro-defects is modeled as an equivalent viscoelastic medium while the macro-joints in the rock mass are modeled explicitly as physical discontinuities.Viscoelastic properties of a micro-defected sedimentary rock are obtained by longitudinally impacting a cored long sedimentary rod with a pendulum.Wave propagation coefficient and dynamic viscoelastic modulus are measured.The EDDM is then successfully employed to analyze the wave propagation across macro-joint in VRM.The effect of the rock viscosity on the stress wave propagation is evaluated by comparing the results of VRM from the presented EDDM with those of an elastic rock mass (ERM) from the conventional displacement discontinuity method (CDDM).The CDDM is a special case of the EDDM under the condition that the rock viscosity is ignored.Comparison of the reflected and transmitted waves shows that the essential rock viscosity has a significant effect on stress wave attenuation.When a short propagation distance of a stress wave is considered,the results obtained from the CDDM approximate to the EDDM solutions,however,when the propagation distance is sufficiently long relative to the wavelength,the effect of rock viscosity on the stress wave propagation cannot be ignored.
基金supported by the Fundamental Research Funds for the Central Universities (HEUCF100307)
文摘An unstructured finite volume time domain method (UFVTDM) is proposed to simulate stress wave propagation. The original variables of displacement and stress are solved based on the dynamic equilibrium equations. An Euler explicit and unstructured finite volume method is used for time and spacial terms respectively. The displacements are stored on the cell vertex and a vertex based finite volume is formed with the integral surface and the stress is assumed as uniform in the cell. This is some similar with the stager grid method in computational fluid dynamics. Several cases are used to show the capability of the algorithm.
文摘An experimental-numerical method for measuring dynamic crack propagatingvelocities under stress wave loading is established in this paper. The experiments of thethree-point bend specimen are done on the improved Hopkinson bar. Deflection of loading point,dynamic load and instantaneous crack length are measured, then crack propagating velocities arecalculated. Experiments on 40Cr steel show that the results given by this method have a goodagreement with that obtained by the resistance fracture gage method. Therefore this method isfeasible for measuring crack propagating velocities under high loading rate and will have wideapplication.
基金The project supported by the National Natural Science Foundation of China
文摘Using the complex variable method and conformal mapping,scat- tering of flexural waves and dynamic stress concentrations in Mindlin's thick plates with a cutout have been studied.The general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of cutouts is obtained. Applying the orthogonal function expansion technique,the dynamic stress problem can be reduced into the solution of a set of infinite algebraic equations.As examples, numerical results for the dynamic stress concentration factor in Mindlin's plates with a circular,elliptic cutout are graphically presented in sequence.
文摘By using the complex variable method and conformal mapping, the diffraction of flexural waves and dynamic stress concentrations in thick plates with a cavity have been studied. A general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of an arbitrary cavity is obtained. By employing the orthogonal function expansion technique, the dynamic stress problem can be reduced to the solution of an infinite algebraic equation series. As an example, the numerical results for the dynamic stress concentration factor in thick plates with a circular, elliptic cavity are graphically presented. The numerical results are discussed.
文摘Boundary Element Method (BEM) is employed to run theoretical analsis and numerical calculation of dif-fraction of elastic wave and dynamic stress concentration in an infinite then plate with a cireular hole. Based on the work equivalent law of dynamics,boundary integral equation is established for flexural waves of thin plate. Calculation formulas of influence coefficients are derived using Mathematica software and numerical results are obtained for dynam-ic stress conoentration factors in a then plate with a circular hole.
文摘In this paper, based on complex variables and conformal mapping methods, using the refined dynamic equation of plates, elastic wave scattering and dynamic stress concentrations in plates with two cutouts were studied. Applying the orthogonal function expansion method, the problem to be solved can be reduced into the solution of a set of infinite algebraic equations. According to free boundary conditions, numerical results of dynamic moment concentration factors in thick plates with two circular cutouts analyze that: there will be more complex interaction changes between two-cutout situation than single cutout situation. In the case of low frequency or high frequency and thin plate, the hole-spacing in the absence of coupling interactions was larger or smaller. The numerical results and method can be used to analyze the dynamics and strength of plate-like structures.
基金Project supported by the National Natural Science Foundation of China(Nos.59978005 and 10232024)the National Distinguished Youth Fund of China(No.10025212).
文摘A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.
基金National Natural Science Foundation of China under Grants (51278327)the Tianjin Research Program of Application Foundation and Advanced Technology (14JCYBJC21900)
文摘The indirect boundary element method （IBEM） is developed to solve the scattering of plane SH-waves by a lined tunnel in elastic wedge space. According to the theory of single-layer potential, the scattered-wave field can be constructed by applying virtual uniform loads on the surface of lined tunnel and the nearby wedge surface. The densities of virtual loads can be solved by establishing equations through the continuity conditions on the interface and zero-traction conditions on free surfaces. The total wave field is obtained by the superposition of free field and scattered-wave field in elastic wedge space. Numerical results indicate that the IBEM can solve the diffraction of elastic wave in elastic wedge space accurately and effi- ciently. The wave motion feature strongly depends on the wedge angle, the angle of incidence, incident frequency, the location of lined tunnel, and material parameters. The waves interference and amplification effect around the tunnel in wedge space is more significant, causing the dynamic stress concentration factor on rigid tunnel and the displacement amplitude of flexible tunnel up to 50.0 and 17.0, respectively, more than double that of the case of half-space. Hence, considerable attention should be paid to seismic resistant or anti-explosion design of the tunnel built on a slope or hillside.
基金The project supported bythe Committee of Science and Technology of Shanghai and Tongji University
文摘A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks.Numerical results are presented for the dynamic stress in- tensity factors.The effects of the wave type,wave frequency,wave incidence angle,and crack spacing on the dynamic stress intensity factors are analyzed in detail.
基金supported by the National Natural Science Foundation of China(11372308 and 11372307)the Fundamental Research Funds for the Central Universities(WK2480000001)
文摘A virtual Taylor impact of cellular materials is analyzed with a wave propagation technique, i.e. the Lagrangian analysis method, of which the main advantage is that no pre-assumed constitutive relationship is required. Time histories of particle velocity, local strain, and stress profiles are calculated to present the local stress-strain history curves, from which the dynamic stress-strain states are obtained. The present results reveal that the dynamic-rigid-plastic hardening （D-R-PH） material model introduced in a previous study of our group is in good agreement with the dynamic stress-strain states under high loading rates obtained by the Lagrangian analysis method. It directly reflects the effectiveness and feasibility of the D-R-PH material model for the cellular materials under high loading rates.
