In terms of asymmetrical three-dimensional distribution(ID) of luminous intensity(LI) of light-emitting-diode(LED),a testing system was conducted in this study. Design and principle of the testing system were introduc...In terms of asymmetrical three-dimensional distribution(ID) of luminous intensity(LI) of light-emitting-diode(LED),a testing system was conducted in this study. Design and principle of the testing system were introduced. 31 photometers were placed on a concentric circle,and all of them were used to gather LI data of LED at the same time. The data acquisition card(DAC) was used to gather multichannel data and controlled motor. Experimental results indicated that the testing system had achieved the goal of testing three-dimensional distribution of LI. And each parameter could meet the requirements of industrial production and measurement.展开更多
This paper presents a formulation for three-dimensional elasto-dynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack...This paper presents a formulation for three-dimensional elasto-dynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack is determined by solving a Fredholm integral equation of the first kind. The results of this paper are very close to those given by the two-dimensional dual integral equation method.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
Large-scale topography, such as a seamount, substantially impacts low-frequency sound propagation in an ocean waveguide, limiting the application of low-frequency acoustic detecting techniques. A three-dimensional(3D)...Large-scale topography, such as a seamount, substantially impacts low-frequency sound propagation in an ocean waveguide, limiting the application of low-frequency acoustic detecting techniques. A three-dimensional(3D) coupledmode model is developed to calculate the acoustic field in an ocean waveguide with seamount topography and analyze the3D effect. In this model, a correction is introduced in the bottom boundary, theoretically making the acoustic field satisfy the energy conservation. Furthermore, a large azimuth angle calculation range is obtained by using the operator theory and higher-order Pade approximation. Additionally, the model has advantages related to the coupling mode and parabolic equation theory. The couplings corresponding to the effects of range-dependent environment are fully considered, and the numerical implementation is kept feasible. After verifying the accuracy and reliability of the model, low-frequency sound propagation characteristics in the seamount environment are analyzed. The results indicate lateral variability in bathymetry can lead to out-of-plane effects such as the horizontal refraction phenomenon, while the coupling effect tends to restore the abnormal sound field and produces acoustic field diffraction behind the seamount. This model effectively considers the effects of the horizontal refraction and coupling, which are proportional to the scale of the seamount.展开更多
In order to reduce the noise and vibration of the diesel engine,it is crucial to exactly identify the engine noise source character.Based on "two-microphone" method,the sound intensity measurement of a vehic...In order to reduce the noise and vibration of the diesel engine,it is crucial to exactly identify the engine noise source character.Based on "two-microphone" method,the sound intensity measurement of a vehicle four-stroke diesel engine was carried out in a hemi-anechoic chamber.Then the sound intensity contour maps were obtained from the measurement results and the main noise components of different frequencies on all the measurement surfaces were picked out to construct contour maps.By analysizing the relationship between the characteristics of contour maps and the space distribution of the engine compartment,the major sources of the exterior radiation noise of the diesel engine were identified.The results provided a creditable basis for improving the noise performance of the engine in the next phase.展开更多
The three-dimensional electron-electron correlation in an elliptically polarized laser field is investigated based on a semiclassical model. Asymmetry parameter α of the correlated electron momentum distribution is u...The three-dimensional electron-electron correlation in an elliptically polarized laser field is investigated based on a semiclassical model. Asymmetry parameter α of the correlated electron momentum distribution is used to quantitatively describe the electron-electron correlation. The dependence of α on elliptieity e is totally different in three directions. For the z direction （maJor polarization direction）, α first increases and reaches a maximum at ε = 0.275, then it decreases quickly. For the y direction in which the laser field is always absent, the ellipticity has a minor effect, and the asymmetry parameter fluctuates around α = -0.15. However, for the x direction （minor polarization direction）, α increases monotonously with ellipticity though starts from the same value as in the y direction when ε = 0. The behavior of α in the x direction actually indicates a transformation from the Coulomb interaction dominated correlation to the laser field dominated correlation. Therefore, our work provides an efficient way to control the three-dimensional electron electron correlation via an elliptically polarized intense laser field.展开更多
The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the rat...The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.展开更多
Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)h...Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)has covered whole south China.展开更多
For anti-jamming and anti-countermeasure techniques of the sonar receiver, the response characteristics of the automatic gain control （AGC） circuit and the survivability of the prime circuit under strong interferenc...For anti-jamming and anti-countermeasure techniques of the sonar receiver, the response characteristics of the automatic gain control （AGC） circuit and the survivability of the prime circuit under strong interference are analyzed by simulations and experiments. An AGC simulation model based on the voltage control amplifier VCA810 prototype is proposed. Then static and dynamic simulations are realized with single frequency signal and linear frequency modulated （LFM） signal commonly used in the active sonar. Based on intense sound pulse （ISP） interference experiments, the real-time response characteristics of each module of the receiver are studied to verify the correctness of the model as well as the simulation results. Simulation and experiment results show that, under 252 dB/20 μs ISP interference, the specific sonar receiver will produce sustained cut top oscillation above 30 ms, which may affect the receiver and block the regular sonar signal.展开更多
A three-dimensional（3D） parabolic equation（PE） model for sound propagation in a seismo-acoustic waveguide is developed in Cartesian coordinates, with x, y, and z representing the marching direction, the longitudin...A three-dimensional（3D） parabolic equation（PE） model for sound propagation in a seismo-acoustic waveguide is developed in Cartesian coordinates, with x, y, and z representing the marching direction, the longitudinal direction, and the depth direction, respectively. Two sets of 3D PEs for horizontally homogenous media are derived by rewriting the 3D elastic motion equations and simultaneously choosing proper dependent variables. The numerical scheme is for now restricted to the y-independent bathymetry. Accuracy of the numerical scheme is validated, and its azimuthal limitation is analyzed. In addition, effects of horizontal refraction in a wedge-shaped waveguide and another waveguide with a polyline bottom are illustrated. Great efforts should be made in future to provide this model with the ability to handle arbitrarily irregular fluid-elastic interfaces.展开更多
A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the ...A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.展开更多
The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental soluti...The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropic materials.展开更多
Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its face...Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its faces. Transform methods are used to reduce the boundary value problem to a single integral equation that can be solved by the Wiener-Hopf technique. The Cagniard-de Hoop method is employed to invert the transforms. An exact expression is derived for the mode I stress intensity factor as a function of time and position along the crack edge. Some features of the solution are discussed through numerical results.展开更多
Since the Xingtai （邢台） earthquake in 1966, China Earthquake Administration has carried out a survey campaign along more than thirty deep seismic sounding （DSS） profiles altogether about twenty thousand kilometer...Since the Xingtai （邢台） earthquake in 1966, China Earthquake Administration has carried out a survey campaign along more than thirty deep seismic sounding （DSS） profiles altogether about twenty thousand kilometers long in North China to study the velocity structure of the crust and the upper mantle in this region, and has obtained a great number of research findings. However, these researches have not provided a 3D velocity structure model of the crust of North China and cannot provide seismic evidence for the study of the deep tectonic characteristics of the crust of the whole region. Hence, based on the information from the published data of the DSS profiles, we have chosen 14 profiles to obtain a 3D velocity structure model of North China using the vectorization function of the GIS software （Arc/Info） and the Kriging data gridding method. With this velocity structure model, we have drawn the following conclusions： （1） The P-wave velocity of the uppermost crust of North China changes dramatically, exhibiting a complicated velocity structure in plane view. It can be divided into three velocity zones mainly trending towards north-west. In the research area, the lowest-velocity zones lie in the Haihe （海河） plain and Bohai （渤海） Bay. Although the geological structure of the sedimentary overburden in the study area is somewhat inherited by the upper crust, there are still several differences between them. （2） Generally, the P-wave velocity of the crust increases with depth in the study area, but there still exists local velocity reversion. In the east, low-velocity anomalies of the Haihe plain gradually disappear with increasing depth, and the Shanxi （山西） graben in the west is mainly characterized by relatively low velocity anomalies. Bounded by the Taihang （太行） Mountains, the eastern and western parts differ in structural trend of stratum above the crystalline basement. The structural trend of the Huanghuaihai （黄淮海） block in the east is mainly north-east, while that of the Shanxi block and the eastern edge of the Ordos block is mainly north-west. （3） According to the morphological features of Moho, the crust of the study area can be divided into six blocks. In the Shanxi block, Moho apppears like a nearly south-north trending depression belt with a large crustal thickness. In the southern edge of the Inner Mongolia block and the south of the Yanshan （燕山） block,the Moho exhibits a feature of fold belt, trending nearly towards east-west. In the eastern edge of the Ordos block, the structure of Moho is relatively complex, presenting a pattern of fold trending nearly towards north-west with alternating convexes and concaves. Beneath the Huanghuaihai block, the middle and northern parts of the North China rift zone, the Moho is the shallowest in the entire region, with alternating uplifts and depressions in its shape. For the anteclise zone in the west of Shandong （山东） Province, the Moho is discontinuous for the fault depression extending in the north-west direction along Zaozhuang （枣庄） -Qufu （曲阜）.展开更多
In this paper, a new semi-analytical and semi-engineering method of the closed form solution of stress intensity factors （SIFs） of cracks emanating from a surface semi-spherical cavity in a finite body is derived us...In this paper, a new semi-analytical and semi-engineering method of the closed form solution of stress intensity factors （SIFs） of cracks emanating from a surface semi-spherical cavity in a finite body is derived using the energy release rate theory. A mode of crack opening displacements of a normal slice is established, and the normal slice relevant functions are introduced. The proposed method is both effective and accurate for the problem of three-dimensional cracks emanating from a surface cavity. A series of useful results of SIFs are obtained.展开更多
Objective: The purpose of this study was to investigate the impact of intensity modulated radiotherapy (IMRT) on surface doses for brain, abdomen and pelvis deep located tumors treated with 6 MV photon and to evaluate...Objective: The purpose of this study was to investigate the impact of intensity modulated radiotherapy (IMRT) on surface doses for brain, abdomen and pelvis deep located tumors treated with 6 MV photon and to evaluate the skin dose calculation accuracy of the XIO 4.04 treatment planning system. Methods: More investigations for the influences of IMRT on skin doses would increase its applications for many treatment sites. Measuring skin doses in real treatment situations would reduce the uncertainty of skin dose prediction. In this work a pediatric human phantom was covered by a layer of 1 mm bolus at three treatment sites and thermoluminescent dosimeter (TLD) chips were inserted into the bolus at each treatment site before CT scan. Two different treatment plans [three-dimensional conformal radiation therapy (3DCRT) and IMRT] for each treatment sites were performed on XIO 4.04 treatment planning system using superposition algorism. Results: The results showed that the surface doses for 3DCRT were higher than the surface doses in IMRT by 1.6%, 2.5% and 3.2% for brain, abdomen and pelvis sites respectively. There was good agreement between measured and calculated surface doses, where the calculated surface dose was 15.5% for brain tumor calculated with 3DCRT whereas the measured surface dose was 12.1%. For abdomen site the calculated surface dose for IMRT treatment plan was 16.5% whereas the measured surface dose was 12.6%. Conclusion: The skin dose in IMRT for deep seated tumors is lower than that in 3DCRT which is another advantage for the IMRT. The TLD readings showed that the difference between the calculated and measured point dose is negligible. The superposition calculation algorism of the XIO 4.04 treatment planning system modeled the superficial dose well.展开更多
文摘In terms of asymmetrical three-dimensional distribution(ID) of luminous intensity(LI) of light-emitting-diode(LED),a testing system was conducted in this study. Design and principle of the testing system were introduced. 31 photometers were placed on a concentric circle,and all of them were used to gather LI data of LED at the same time. The data acquisition card(DAC) was used to gather multichannel data and controlled motor. Experimental results indicated that the testing system had achieved the goal of testing three-dimensional distribution of LI. And each parameter could meet the requirements of industrial production and measurement.
基金The project supported by the National Natural Science Foundation of China (K19672007)
文摘This paper presents a formulation for three-dimensional elasto-dynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack is determined by solving a Fredholm integral equation of the first kind. The results of this paper are very close to those given by the two-dimensional dual integral equation method.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Project supported by the National Natural Science Foundation of China(Grant No.11804360)the IACAS Frontier Exploration Project(Grant No.QYTS202103)the Key Laboratory Foundation of Acoustic Science and Technology(Grant No.2021-JCJQ-LB-066-08)。
文摘Large-scale topography, such as a seamount, substantially impacts low-frequency sound propagation in an ocean waveguide, limiting the application of low-frequency acoustic detecting techniques. A three-dimensional(3D) coupledmode model is developed to calculate the acoustic field in an ocean waveguide with seamount topography and analyze the3D effect. In this model, a correction is introduced in the bottom boundary, theoretically making the acoustic field satisfy the energy conservation. Furthermore, a large azimuth angle calculation range is obtained by using the operator theory and higher-order Pade approximation. Additionally, the model has advantages related to the coupling mode and parabolic equation theory. The couplings corresponding to the effects of range-dependent environment are fully considered, and the numerical implementation is kept feasible. After verifying the accuracy and reliability of the model, low-frequency sound propagation characteristics in the seamount environment are analyzed. The results indicate lateral variability in bathymetry can lead to out-of-plane effects such as the horizontal refraction phenomenon, while the coupling effect tends to restore the abnormal sound field and produces acoustic field diffraction behind the seamount. This model effectively considers the effects of the horizontal refraction and coupling, which are proportional to the scale of the seamount.
基金supported by programfor the Top Young Academic Leaders of Higher Learning Institutions of Shanxi(2009)Natural Science Foundation of Shanxi Province,China(No.2010011031-2)
文摘In order to reduce the noise and vibration of the diesel engine,it is crucial to exactly identify the engine noise source character.Based on "two-microphone" method,the sound intensity measurement of a vehicle four-stroke diesel engine was carried out in a hemi-anechoic chamber.Then the sound intensity contour maps were obtained from the measurement results and the main noise components of different frequencies on all the measurement surfaces were picked out to construct contour maps.By analysizing the relationship between the characteristics of contour maps and the space distribution of the engine compartment,the major sources of the exterior radiation noise of the diesel engine were identified.The results provided a creditable basis for improving the noise performance of the engine in the next phase.
基金Supported by the National Key Program for S&T Research and Development under Grant No 2016YFA0401100the National Basic Research Program of China under Grant No 2013CB922201the National Natural Science Foundation of China under Grant Nos 11504215,11374197,11334009 and 11425414
文摘The three-dimensional electron-electron correlation in an elliptically polarized laser field is investigated based on a semiclassical model. Asymmetry parameter α of the correlated electron momentum distribution is used to quantitatively describe the electron-electron correlation. The dependence of α on elliptieity e is totally different in three directions. For the z direction （maJor polarization direction）, α first increases and reaches a maximum at ε = 0.275, then it decreases quickly. For the y direction in which the laser field is always absent, the ellipticity has a minor effect, and the asymmetry parameter fluctuates around α = -0.15. However, for the x direction （minor polarization direction）, α increases monotonously with ellipticity though starts from the same value as in the y direction when ε = 0. The behavior of α in the x direction actually indicates a transformation from the Coulomb interaction dominated correlation to the laser field dominated correlation. Therefore, our work provides an efficient way to control the three-dimensional electron electron correlation via an elliptically polarized intense laser field.
基金the National Natural Science Foundation of China( No.K19672007)
文摘The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.
基金co-supported by the China Geological Survey Project(DD20190012 and DD20160082)
文摘Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)has covered whole south China.
基金supported by the National Natural Science Foundation of China (10974154)the National Innovation Project of China for Undergraduates (101069935)
文摘For anti-jamming and anti-countermeasure techniques of the sonar receiver, the response characteristics of the automatic gain control （AGC） circuit and the survivability of the prime circuit under strong interference are analyzed by simulations and experiments. An AGC simulation model based on the voltage control amplifier VCA810 prototype is proposed. Then static and dynamic simulations are realized with single frequency signal and linear frequency modulated （LFM） signal commonly used in the active sonar. Based on intense sound pulse （ISP） interference experiments, the real-time response characteristics of each module of the receiver are studied to verify the correctness of the model as well as the simulation results. Simulation and experiment results show that, under 252 dB/20 μs ISP interference, the specific sonar receiver will produce sustained cut top oscillation above 30 ms, which may affect the receiver and block the regular sonar signal.
基金Project supported by the National Nature Science Foundation of China(Grant Nos.11234002 and 11704337)the National Key Research Program of China(Grant No.2016YFC1400100)
文摘A three-dimensional（3D） parabolic equation（PE） model for sound propagation in a seismo-acoustic waveguide is developed in Cartesian coordinates, with x, y, and z representing the marching direction, the longitudinal direction, and the depth direction, respectively. Two sets of 3D PEs for horizontally homogenous media are derived by rewriting the 3D elastic motion equations and simultaneously choosing proper dependent variables. The numerical scheme is for now restricted to the y-independent bathymetry. Accuracy of the numerical scheme is validated, and its azimuthal limitation is analyzed. In addition, effects of horizontal refraction in a wedge-shaped waveguide and another waveguide with a polyline bottom are illustrated. Great efforts should be made in future to provide this model with the ability to handle arbitrarily irregular fluid-elastic interfaces.
文摘A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.
基金Project supported by the Program for New Century Excellent Talents in University of Henan Province (HANCET)
文摘The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropic materials.
基金The project supported by the Guangdong Provincial Natural Science Foundationthe Science Foundation of Shantou University
文摘Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its faces. Transform methods are used to reduce the boundary value problem to a single integral equation that can be solved by the Wiener-Hopf technique. The Cagniard-de Hoop method is employed to invert the transforms. An exact expression is derived for the mode I stress intensity factor as a function of time and position along the crack edge. Some features of the solution are discussed through numerical results.
基金This paper is supported by the National Natural Science Foundation of China (No.40434010)the Focused Subject Program of Beijing (No. XK104910589).
文摘Since the Xingtai （邢台） earthquake in 1966, China Earthquake Administration has carried out a survey campaign along more than thirty deep seismic sounding （DSS） profiles altogether about twenty thousand kilometers long in North China to study the velocity structure of the crust and the upper mantle in this region, and has obtained a great number of research findings. However, these researches have not provided a 3D velocity structure model of the crust of North China and cannot provide seismic evidence for the study of the deep tectonic characteristics of the crust of the whole region. Hence, based on the information from the published data of the DSS profiles, we have chosen 14 profiles to obtain a 3D velocity structure model of North China using the vectorization function of the GIS software （Arc/Info） and the Kriging data gridding method. With this velocity structure model, we have drawn the following conclusions： （1） The P-wave velocity of the uppermost crust of North China changes dramatically, exhibiting a complicated velocity structure in plane view. It can be divided into three velocity zones mainly trending towards north-west. In the research area, the lowest-velocity zones lie in the Haihe （海河） plain and Bohai （渤海） Bay. Although the geological structure of the sedimentary overburden in the study area is somewhat inherited by the upper crust, there are still several differences between them. （2） Generally, the P-wave velocity of the crust increases with depth in the study area, but there still exists local velocity reversion. In the east, low-velocity anomalies of the Haihe plain gradually disappear with increasing depth, and the Shanxi （山西） graben in the west is mainly characterized by relatively low velocity anomalies. Bounded by the Taihang （太行） Mountains, the eastern and western parts differ in structural trend of stratum above the crystalline basement. The structural trend of the Huanghuaihai （黄淮海） block in the east is mainly north-east, while that of the Shanxi block and the eastern edge of the Ordos block is mainly north-west. （3） According to the morphological features of Moho, the crust of the study area can be divided into six blocks. In the Shanxi block, Moho apppears like a nearly south-north trending depression belt with a large crustal thickness. In the southern edge of the Inner Mongolia block and the south of the Yanshan （燕山） block,the Moho exhibits a feature of fold belt, trending nearly towards east-west. In the eastern edge of the Ordos block, the structure of Moho is relatively complex, presenting a pattern of fold trending nearly towards north-west with alternating convexes and concaves. Beneath the Huanghuaihai block, the middle and northern parts of the North China rift zone, the Moho is the shallowest in the entire region, with alternating uplifts and depressions in its shape. For the anteclise zone in the west of Shandong （山东） Province, the Moho is discontinuous for the fault depression extending in the north-west direction along Zaozhuang （枣庄） -Qufu （曲阜）.
文摘In this paper, a new semi-analytical and semi-engineering method of the closed form solution of stress intensity factors （SIFs） of cracks emanating from a surface semi-spherical cavity in a finite body is derived using the energy release rate theory. A mode of crack opening displacements of a normal slice is established, and the normal slice relevant functions are introduced. The proposed method is both effective and accurate for the problem of three-dimensional cracks emanating from a surface cavity. A series of useful results of SIFs are obtained.
文摘Objective: The purpose of this study was to investigate the impact of intensity modulated radiotherapy (IMRT) on surface doses for brain, abdomen and pelvis deep located tumors treated with 6 MV photon and to evaluate the skin dose calculation accuracy of the XIO 4.04 treatment planning system. Methods: More investigations for the influences of IMRT on skin doses would increase its applications for many treatment sites. Measuring skin doses in real treatment situations would reduce the uncertainty of skin dose prediction. In this work a pediatric human phantom was covered by a layer of 1 mm bolus at three treatment sites and thermoluminescent dosimeter (TLD) chips were inserted into the bolus at each treatment site before CT scan. Two different treatment plans [three-dimensional conformal radiation therapy (3DCRT) and IMRT] for each treatment sites were performed on XIO 4.04 treatment planning system using superposition algorism. Results: The results showed that the surface doses for 3DCRT were higher than the surface doses in IMRT by 1.6%, 2.5% and 3.2% for brain, abdomen and pelvis sites respectively. There was good agreement between measured and calculated surface doses, where the calculated surface dose was 15.5% for brain tumor calculated with 3DCRT whereas the measured surface dose was 12.1%. For abdomen site the calculated surface dose for IMRT treatment plan was 16.5% whereas the measured surface dose was 12.6%. Conclusion: The skin dose in IMRT for deep seated tumors is lower than that in 3DCRT which is another advantage for the IMRT. The TLD readings showed that the difference between the calculated and measured point dose is negligible. The superposition calculation algorism of the XIO 4.04 treatment planning system modeled the superficial dose well.
