By using integral transform methods, the Green(s functions of horizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equ...By using integral transform methods, the Green(s functions of horizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above- mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous conditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green(s functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper.展开更多
In order to simplify the boundary conditions of pavement temperature field,the "Environment-Surface" system which considered the natural environment and pavement surface was established.Based on this system,...In order to simplify the boundary conditions of pavement temperature field,the "Environment-Surface" system which considered the natural environment and pavement surface was established.Based on this system,the partial differential equations of the one-dimensional heat conduction in the pavement were established on the basis of the heat transfer theory.Furthermore,the function forms of the initial and boundary conditions of the equations were created through the field experiments.The general solution of the pavement one-dimensional heat conduction partial differential equations was acquired by using Green's function,and the explicit expression of pavement temperature field under specific constraint conditions was derived.For the purpose of analysis,the pavement temperatures in different seasons were calculated using the explicit expression of pavement temperature field,and the calculation accuracy was analyzed through the comparison between measured and calculated values.Then,the relationship between fitting accuracy and calculation accuracy of pavement temperatures was analyzed.The analysis results show that: the usage of "Environment-Surface" system simplifies the calculation of pavement temperature field; the relative error between calculated and measured values is generally less than 7% and is seldom influenced by seasons; there is a positive correlation between the calculation accuracy and the fitting accuracy of pavement surface temperature; high fitting accuracy would result in less error of pavement temperature prediction.展开更多
The flow behavior in porous media with threshold pressure gradient(TPG) is more complex than Darcy flow and the equations of motion, and outer boundary and inner boundary with TPG are also different from Darcy flow fo...The flow behavior in porous media with threshold pressure gradient(TPG) is more complex than Darcy flow and the equations of motion, and outer boundary and inner boundary with TPG are also different from Darcy flow for unsteady flow of a producing well in a reservoir. An analytic method to solve this kind of problem is in a need of reestablishment. The classical method of Green's function and Newman product principle in a new way are used to solve the unsteady state flow problems of various shapes of well and reservoir while considering the TPG. Four Green's functions of point, line, band and circle while considering the TPG are achieved. Then, two well models of vertical well and horizontal well are built and simultaneously the function to calculate the moving boundary of each well model is provided. The results show that when considering TPG the pressure field is much different, which has a sudden pressure change, with a moving boundary in it. And the moving boundary of each well model increases with time but slows down rapidly, especially when the TGP is large.展开更多
By virtue of a complete set of two displacement potentials,an analytical derivation of the elastostatic Green’s functions of an exponentially graded transversely isotropic substrate–coating system is presented.Three...By virtue of a complete set of two displacement potentials,an analytical derivation of the elastostatic Green’s functions of an exponentially graded transversely isotropic substrate–coating system is presented.Three-dimensional point–load and patch–load Green’s functions for stresses and displacements are given in line-integral representations.The formulation includes a complete set of transformed stress–potential and displacement–potential relations,with utilizing Fourier series and Hankel transforms.As illustrations,the present Green’s functions are degenerated to the special cases such as an exponentially graded half-space and a homogeneous two-layered half-space Green’s functions.Because of complicated integrand functions,the integrals are evaluated numerically and for numerical computation of the integrals,a robust and effective methodology is laid out which gives the necessary account of the presence of singularities of integration.Comparisons of the existing numerical solutions for homogeneous two-layered isotropic and transversely isotropic half-spaces are made to confirm the accuracy of the present solutions.Some typical numerical examples are also given to show the general features of the exponentially graded two-layered half-space Green’s functions that the effect of degree of variation of material properties will be recognized.展开更多
This work presented the development and validation of an analytical method to predict the transient temperature field in the asphalt pavement.The governing equation for heat transfer was based on heat conduction radia...This work presented the development and validation of an analytical method to predict the transient temperature field in the asphalt pavement.The governing equation for heat transfer was based on heat conduction radiation and convection.An innovative time-dependent function was proposed to predict the pavement surface temperature with solar radiation and air temperature using dimensional analysis in order to simplify the complex heat exchange on the pavement surface.The parameters for the time-dependent pavement surface temperature function were obtained through the regression analysis of field measurement data.Assuming that the initial pavement temperature distribution was linear and the influence of the base course materials on the temperature of the upper asphalt layers was negligible,a close-form analytical solution of the temperature in asphalt layers was derived using Green's function.Finally,two numerical examples were presented to validate the model solutions with field temperature measurements.Analysis results show that the solution accuracy is in agreement with field data and the relative errors at a shallower depth are greater than those at a deeper one.Although the model is not sensitive to dramatic changes in climatic factors near the pavement surface,it is applicable for predicting pavement temperature field in cloudless days.展开更多
Based on Biot’s theory and considering the properties of a cavity,the boundary integral equations for the numerical simulation of wave scattering around a cavity with a circular cross-section embedded in saturated so...Based on Biot’s theory and considering the properties of a cavity,the boundary integral equations for the numerical simulation of wave scattering around a cavity with a circular cross-section embedded in saturated soil are obtained using integral transform methods.The Cauchy type singularity of the boundary integral equation is discussed.The effectiveness of the properties of soil mass and incident field on the dynamic stress concentration and pore pressure concentration around a cavity is analyzed.Our results are in good agreement with the existing solution.The numerical results of this work show that the dynamic stress concentration and pore pressure concentration are influenced by the degree of fluid–solid coupling as well as the pore compressibility and water permeability of saturated soil.With increased degree of fluid–solid coupling,the dynamic stress concentration improves from 1.87 to 3.42 and the scattering becomes more significant.With decreased index of soil mass compressibility,the dynamic stress concentration increases and its maximum reaches 3.67.The dynamic stress concentration increases from 1.64 to 3.49 and pore pressure concentration improves from 0.18 to 0.46 with decreased water permeability of saturated soil.展开更多
On July 20, 1995 an M_L= 4.1 earthquake occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, microseismicity sharp...On July 20, 1995 an M_L= 4.1 earthquake occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, microseismicity sharply increased in the basin. This earthquake sequence was recorded by Sino-European Cooperative Huailai Digital Seismograph Network (HDSN) and the hypocenters were precisely located. Using 8 aftershocks as empirical Green(s function (EGF) we calculated the source time function (STF) of the ML=4.1 event respectively, and by superposing these results, we obtained the average time function which has better signal to noise ratio. The final result shows that the ML=4.1 earthquake consisted of two events, one of which was stronger, the other was weaker, STFs from each station demonstrate an obvious (seismic Doppler effect(. Through a trial-and-error method, the following results were drawn: for the 1st event, the rupture length is 0.44 km, rupture velocity is 4.0 km/s, while for the 2nd event those parameters are 0.25 km, 3.0 km/s, the two events ruptured in the direction of 140° and 90° respectively (measured from the strike of the fault, anticlockwise positive), the two events were separated by 0.57 km in space, and 0.09 s in time. Using single EGF, we also retrieved the source time function of other 13 events with ML(2.1 in the Huailai earthquake sequence. These results were reached: each of the 13 events consisted of single event, with the duration of the source time functions ranging between 0.05～0.16 s. Using far-field ground displacement measurement method, the scalar seismic moment, stress drop and rupture radii of 25 events with 0.95(ML(4.1 were obtained. The trend that scalar seismic moment and stress drop increase with magnitude was obvious. In the magnitude range 0.9(ML(2.4, the obtained radii show no correlativity with magnitude, we tend to believe, in such small magnitude range, precise measurement for radii of such small events is beyond the data used in this study.展开更多
基金State Natural Science Foundation (59879012) and Doctoral Foundation from State Education Commission (98024832).
文摘By using integral transform methods, the Green(s functions of horizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above- mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous conditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green(s functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper.
基金Projects(2012zzts019,2012QNZT048)supported by the Fundamental Research Funds for the Central Universities of Central South University,ChinaProject(201306370121)supported by the State Scholarship Fund of China+3 种基金Project(JT20090898002)supported by Traffic Technology Fund of Hainan Province,ChinaProject(2012M521563)supported by the China Postdoctoral Science FoundationProject(51248006)supported by The National Natural Science Foundation of ChinaProject(511114)supported by the Natural Science Foundation of Hainan Province,China
文摘In order to simplify the boundary conditions of pavement temperature field,the "Environment-Surface" system which considered the natural environment and pavement surface was established.Based on this system,the partial differential equations of the one-dimensional heat conduction in the pavement were established on the basis of the heat transfer theory.Furthermore,the function forms of the initial and boundary conditions of the equations were created through the field experiments.The general solution of the pavement one-dimensional heat conduction partial differential equations was acquired by using Green's function,and the explicit expression of pavement temperature field under specific constraint conditions was derived.For the purpose of analysis,the pavement temperatures in different seasons were calculated using the explicit expression of pavement temperature field,and the calculation accuracy was analyzed through the comparison between measured and calculated values.Then,the relationship between fitting accuracy and calculation accuracy of pavement temperatures was analyzed.The analysis results show that: the usage of "Environment-Surface" system simplifies the calculation of pavement temperature field; the relative error between calculated and measured values is generally less than 7% and is seldom influenced by seasons; there is a positive correlation between the calculation accuracy and the fitting accuracy of pavement surface temperature; high fitting accuracy would result in less error of pavement temperature prediction.
基金Project(51304220) supported by the National Natural Science Foundation of ChinaProject(3144033) supported by the Beijing Natural Science Foundation,ChinaProject(20130007120014) supported by the Specialized Research Fund for the Doctoral Program of Higher Education of China
文摘The flow behavior in porous media with threshold pressure gradient(TPG) is more complex than Darcy flow and the equations of motion, and outer boundary and inner boundary with TPG are also different from Darcy flow for unsteady flow of a producing well in a reservoir. An analytic method to solve this kind of problem is in a need of reestablishment. The classical method of Green's function and Newman product principle in a new way are used to solve the unsteady state flow problems of various shapes of well and reservoir while considering the TPG. Four Green's functions of point, line, band and circle while considering the TPG are achieved. Then, two well models of vertical well and horizontal well are built and simultaneously the function to calculate the moving boundary of each well model is provided. The results show that when considering TPG the pressure field is much different, which has a sudden pressure change, with a moving boundary in it. And the moving boundary of each well model increases with time but slows down rapidly, especially when the TGP is large.
文摘By virtue of a complete set of two displacement potentials,an analytical derivation of the elastostatic Green’s functions of an exponentially graded transversely isotropic substrate–coating system is presented.Three-dimensional point–load and patch–load Green’s functions for stresses and displacements are given in line-integral representations.The formulation includes a complete set of transformed stress–potential and displacement–potential relations,with utilizing Fourier series and Hankel transforms.As illustrations,the present Green’s functions are degenerated to the special cases such as an exponentially graded half-space and a homogeneous two-layered half-space Green’s functions.Because of complicated integrand functions,the integrals are evaluated numerically and for numerical computation of the integrals,a robust and effective methodology is laid out which gives the necessary account of the presence of singularities of integration.Comparisons of the existing numerical solutions for homogeneous two-layered isotropic and transversely isotropic half-spaces are made to confirm the accuracy of the present solutions.Some typical numerical examples are also given to show the general features of the exponentially graded two-layered half-space Green’s functions that the effect of degree of variation of material properties will be recognized.
基金Project(2012zzts019)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(201306370121)supported by State Scholarship Fund of ChinaProject(51248006)supported by the National Natural Science Foundation,China
文摘This work presented the development and validation of an analytical method to predict the transient temperature field in the asphalt pavement.The governing equation for heat transfer was based on heat conduction radiation and convection.An innovative time-dependent function was proposed to predict the pavement surface temperature with solar radiation and air temperature using dimensional analysis in order to simplify the complex heat exchange on the pavement surface.The parameters for the time-dependent pavement surface temperature function were obtained through the regression analysis of field measurement data.Assuming that the initial pavement temperature distribution was linear and the influence of the base course materials on the temperature of the upper asphalt layers was negligible,a close-form analytical solution of the temperature in asphalt layers was derived using Green's function.Finally,two numerical examples were presented to validate the model solutions with field temperature measurements.Analysis results show that the solution accuracy is in agreement with field data and the relative errors at a shallower depth are greater than those at a deeper one.Although the model is not sensitive to dramatic changes in climatic factors near the pavement surface,it is applicable for predicting pavement temperature field in cloudless days.
基金Projects(50969007,51269021) supported by the National Natural Science Foundation of ChinaProjects(20114BAB206012,20133ACB20006) supported by the Natural Science Foundation of Jiangxi Province of China
文摘Based on Biot’s theory and considering the properties of a cavity,the boundary integral equations for the numerical simulation of wave scattering around a cavity with a circular cross-section embedded in saturated soil are obtained using integral transform methods.The Cauchy type singularity of the boundary integral equation is discussed.The effectiveness of the properties of soil mass and incident field on the dynamic stress concentration and pore pressure concentration around a cavity is analyzed.Our results are in good agreement with the existing solution.The numerical results of this work show that the dynamic stress concentration and pore pressure concentration are influenced by the degree of fluid–solid coupling as well as the pore compressibility and water permeability of saturated soil.With increased degree of fluid–solid coupling,the dynamic stress concentration improves from 1.87 to 3.42 and the scattering becomes more significant.With decreased index of soil mass compressibility,the dynamic stress concentration increases and its maximum reaches 3.67.The dynamic stress concentration increases from 1.64 to 3.49 and pore pressure concentration improves from 0.18 to 0.46 with decreased water permeability of saturated soil.
基金National Key Basic Research Project "Mechanism and Prediction of Continental Strong Earthquakes".
文摘On July 20, 1995 an M_L= 4.1 earthquake occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, microseismicity sharply increased in the basin. This earthquake sequence was recorded by Sino-European Cooperative Huailai Digital Seismograph Network (HDSN) and the hypocenters were precisely located. Using 8 aftershocks as empirical Green(s function (EGF) we calculated the source time function (STF) of the ML=4.1 event respectively, and by superposing these results, we obtained the average time function which has better signal to noise ratio. The final result shows that the ML=4.1 earthquake consisted of two events, one of which was stronger, the other was weaker, STFs from each station demonstrate an obvious (seismic Doppler effect(. Through a trial-and-error method, the following results were drawn: for the 1st event, the rupture length is 0.44 km, rupture velocity is 4.0 km/s, while for the 2nd event those parameters are 0.25 km, 3.0 km/s, the two events ruptured in the direction of 140° and 90° respectively (measured from the strike of the fault, anticlockwise positive), the two events were separated by 0.57 km in space, and 0.09 s in time. Using single EGF, we also retrieved the source time function of other 13 events with ML(2.1 in the Huailai earthquake sequence. These results were reached: each of the 13 events consisted of single event, with the duration of the source time functions ranging between 0.05～0.16 s. Using far-field ground displacement measurement method, the scalar seismic moment, stress drop and rupture radii of 25 events with 0.95(ML(4.1 were obtained. The trend that scalar seismic moment and stress drop increase with magnitude was obvious. In the magnitude range 0.9(ML(2.4, the obtained radii show no correlativity with magnitude, we tend to believe, in such small magnitude range, precise measurement for radii of such small events is beyond the data used in this study.
