The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, ...The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, are modeled by the finite elements, and the wave propagation properties of the far field extending to infinity are modeled by the infinite elements. One particular feature of the 2.5D approach is that it enables the computation of the three-dimensional response of the half-space, taking into account the load-moving effect, using only a two-dimensional profile. Although the 2.5D finite/infinite element approach shows a great advantage in studying the wave propagation caused by moving trains, attention should be given to the calculation aspects, such as the rules for mesh establishment, in order to avoid producing inaccurate or erroneous results. In this paper, some essential points for consideration in analysis are highlighted, along with techniques to enhance the speed of the calculations. All these observations should prove useful in making the 2.5D finite/infinite element approach an effective one.展开更多
A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full t...A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number (Ky) domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell's equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in展开更多
Single-grit grinding of a 2.5D woven composite was investigated by the finite-element method(FEM) using a unit-cell model. According to our hypotheses, the axis of the warp yarn was a sinusoidal curve and the cross se...Single-grit grinding of a 2.5D woven composite was investigated by the finite-element method(FEM) using a unit-cell model. According to our hypotheses, the axis of the warp yarn was a sinusoidal curve and the cross section of the weft yarn was shaped like a biconvex lens. AVUMAT subroutine was used to construct the constitutive model of the 2.5D woven composite. The grinding process of the composite was analyzed using an FEM simulation with the ABAQUS/Explicit software. A validation experiment was also carried out. The simulation results showed that a grinding crack was well simulated. Furthermore, the junctions between the warp yarn and weft yarn were found to be seriously damaged and cracks were observed to extend outward along the warp fiber during grinding, in good agreement with the experimental results. In addition,the strain of weft yarns was obviously greater than that of warp yarns when the grinding direction was perpendicular to the weft yarns and parallel to the axis of the warp yarns.These results demonstrate that the mesostructure strongly influences the grinding damage inflicted on woven composites.展开更多
To make better use of 2.5D C/SiC composites in industry, it is necessary to understand the mechanical properties. A finite element model'of 2.5D composites is established, by considering the fiber undulation and the ...To make better use of 2.5D C/SiC composites in industry, it is necessary to understand the mechanical properties. A finite element model'of 2.5D composites is established, by considering the fiber undulation and the porosity in 2.5D C/SiC composites. The fiber direction of warp is defined by cosine function to simulate the undulation of warp, and based on uniform strain assumption, analytical model of the elastic modulus and coefficient of thermal expansion (CTE) for 2.5D C/SiC composites were established by using dual- scale model. The result is found to correlate reasonably well with the predicted results and experimental results. The parametric study also demonstrates the effects of the fiber volume fraction, distance of warp yarn, and porosity in micro-scale on the mechanical properties and the coefficients of thermal expansion.展开更多
An efficient 2.5D finite element numerical modeling approach was developed to simulate wave motions generated in ground by high-speed train passages. Fourier transform with respect to the coordinate in the track direc...An efficient 2.5D finite element numerical modeling approach was developed to simulate wave motions generated in ground by high-speed train passages. Fourier transform with respect to the coordinate in the track direction was applied to re-ducing the three-dimensional dynamic problem to a plane strain problem which has been solved in a section perpendicular to the track direction. In this study, the track structure and supporting ballast layer were simplified as a composite Euler beam resting on the ground surface, while the ground with complicated geometry and physical properties was modeled by 2.5D quadrilateral elements. Wave dissipation into the far field was dealt with the transmitting boundary constructed with fre-quency-dependent dashpots. Three-dimensional responses of track structure and ground were obtained from the wavenumber expansion in the track direction. The simulated wave motions in ground were interpreted for train moving loads traveling at speeds below or above the critical velocity of a specific track-ground system. It is found that, in the soft ground area, the high-speed train operations can enter the transonic range, which can lead to resonances of the track structure and the sup-porting ground. The strong vibration will endanger the safe operations of high-speed train and accelerate the deterioration of railway structure.展开更多
The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to eff...The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to effectively invert these spectral parameters has become the focused area of the complex resistivity inversion. An optimized BP neural network (BPNN) approach based on Quantum Particle Swarm Optimization (QPSO) algorithm was presented, which was able to improve global search ability for complex resistivity multi-parameter nonlinear inversion. In the proposed method, the nonlinear weight adjustment strategy and mutation operator were used to enhance the optimization ability of QPSO algorithm. Implementation of proposed QPSO-BPNN was given, the network had 56 hidden neurons in two hidden layers (the first hidden layer has 46 neurons and the second hidden layer has 10 neurons) and it was trained on 48 datasets and tested on another 5 synthetic datasets. The training and test results show that BP neural network optimized by the QPSO algorithm performs better than the BP neural network without initial optimization on the inversion training and test models, and the mean square error distribution is better. At the same time, a double polarized anomalous bodies model was also used to verify the feasibility and effectiveness of the proposed method, the inversion results show that the QPSO-BP algorithm inversion clearly characterizes the anomalous boundaries and is closer to the values of the parameters.展开更多
This paper presents,for the first time,the consideration of three-dimensional(3D)oblique incident P and SV waves in calculating the 3D seismic response of a lined tunnel embedded in a half-space by the 2.5D finite/inf...This paper presents,for the first time,the consideration of three-dimensional(3D)oblique incident P and SV waves in calculating the 3D seismic response of a lined tunnel embedded in a half-space by the 2.5D finite/infinite element method(FIEM).Firstly,the applicability of the 2.5D FIEM for 3D seismic analysis is summarized.With the exact solutions obtained for the free field in the Appendix,the equivalent seismic forces are rationally computed for the near-field boundary,considering the horizontal and vertical excitations of the Chi-Chi Earthquake.By performing seismic analysis of the half space embedded with a tunnel using the 2.5D FIEM,the time-domain responses of the tunnel are obtained.The accuracy of the present solutions is verified against those of de Barros and Luco.Conclusions drawn from the parametric study include:(1)Stress concentration for the principal stress under oblique incident seismic waves occurs at the polar angles of 0(vault),90,180(inverted arch),and 270of the lining wall.(2)The vault and inverted arch are the weakest parts of the tunnel during earthquakes.(3)The accelerations of the tunnel during earthquakes can be regarded as of the rigid body type.(4)The responses of the tunnel lining caused by SV waves of an earthquake are much more critical than those by P waves.(5)For arbitrary seismic waves,the maximum longitudinal acceleration azmax is of the same order of magnitude as the maximum horizontal acceleration axmax.展开更多
基金Science Council Under Grant No.NSC 89-2211-E-002-020
文摘The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, are modeled by the finite elements, and the wave propagation properties of the far field extending to infinity are modeled by the infinite elements. One particular feature of the 2.5D approach is that it enables the computation of the three-dimensional response of the half-space, taking into account the load-moving effect, using only a two-dimensional profile. Although the 2.5D finite/infinite element approach shows a great advantage in studying the wave propagation caused by moving trains, attention should be given to the calculation aspects, such as the rules for mesh establishment, in order to avoid producing inaccurate or erroneous results. In this paper, some essential points for consideration in analysis are highlighted, along with techniques to enhance the speed of the calculations. All these observations should prove useful in making the 2.5D finite/infinite element approach an effective one.
文摘A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number (Ky) domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell's equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in
基金supported by the National Natural Science Foundation of China(No.51305296)
文摘Single-grit grinding of a 2.5D woven composite was investigated by the finite-element method(FEM) using a unit-cell model. According to our hypotheses, the axis of the warp yarn was a sinusoidal curve and the cross section of the weft yarn was shaped like a biconvex lens. AVUMAT subroutine was used to construct the constitutive model of the 2.5D woven composite. The grinding process of the composite was analyzed using an FEM simulation with the ABAQUS/Explicit software. A validation experiment was also carried out. The simulation results showed that a grinding crack was well simulated. Furthermore, the junctions between the warp yarn and weft yarn were found to be seriously damaged and cracks were observed to extend outward along the warp fiber during grinding, in good agreement with the experimental results. In addition,the strain of weft yarns was obviously greater than that of warp yarns when the grinding direction was perpendicular to the weft yarns and parallel to the axis of the warp yarns.These results demonstrate that the mesostructure strongly influences the grinding damage inflicted on woven composites.
基金Funded by the National Basic Research Program of China,National Natural Science Foundation of China(No.51075204)Aeronautical Science Foundation of China(No.2012ZB52026)+1 种基金Research Fund for the Doctoral Program of Higher Education of China(No.20070287039)NUAA Research Funding(No.NZ2012106)
文摘To make better use of 2.5D C/SiC composites in industry, it is necessary to understand the mechanical properties. A finite element model'of 2.5D composites is established, by considering the fiber undulation and the porosity in 2.5D C/SiC composites. The fiber direction of warp is defined by cosine function to simulate the undulation of warp, and based on uniform strain assumption, analytical model of the elastic modulus and coefficient of thermal expansion (CTE) for 2.5D C/SiC composites were established by using dual- scale model. The result is found to correlate reasonably well with the predicted results and experimental results. The parametric study also demonstrates the effects of the fiber volume fraction, distance of warp yarn, and porosity in micro-scale on the mechanical properties and the coefficients of thermal expansion.
基金Supported by the National Natural Science Foundation of China (Grant No. 10702063) the Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20070335086)
文摘An efficient 2.5D finite element numerical modeling approach was developed to simulate wave motions generated in ground by high-speed train passages. Fourier transform with respect to the coordinate in the track direction was applied to re-ducing the three-dimensional dynamic problem to a plane strain problem which has been solved in a section perpendicular to the track direction. In this study, the track structure and supporting ballast layer were simplified as a composite Euler beam resting on the ground surface, while the ground with complicated geometry and physical properties was modeled by 2.5D quadrilateral elements. Wave dissipation into the far field was dealt with the transmitting boundary constructed with fre-quency-dependent dashpots. Three-dimensional responses of track structure and ground were obtained from the wavenumber expansion in the track direction. The simulated wave motions in ground were interpreted for train moving loads traveling at speeds below or above the critical velocity of a specific track-ground system. It is found that, in the soft ground area, the high-speed train operations can enter the transonic range, which can lead to resonances of the track structure and the sup-porting ground. The strong vibration will endanger the safe operations of high-speed train and accelerate the deterioration of railway structure.
文摘The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to effectively invert these spectral parameters has become the focused area of the complex resistivity inversion. An optimized BP neural network (BPNN) approach based on Quantum Particle Swarm Optimization (QPSO) algorithm was presented, which was able to improve global search ability for complex resistivity multi-parameter nonlinear inversion. In the proposed method, the nonlinear weight adjustment strategy and mutation operator were used to enhance the optimization ability of QPSO algorithm. Implementation of proposed QPSO-BPNN was given, the network had 56 hidden neurons in two hidden layers (the first hidden layer has 46 neurons and the second hidden layer has 10 neurons) and it was trained on 48 datasets and tested on another 5 synthetic datasets. The training and test results show that BP neural network optimized by the QPSO algorithm performs better than the BP neural network without initial optimization on the inversion training and test models, and the mean square error distribution is better. At the same time, a double polarized anomalous bodies model was also used to verify the feasibility and effectiveness of the proposed method, the inversion results show that the QPSO-BP algorithm inversion clearly characterizes the anomalous boundaries and is closer to the values of the parameters.
文摘This paper presents,for the first time,the consideration of three-dimensional(3D)oblique incident P and SV waves in calculating the 3D seismic response of a lined tunnel embedded in a half-space by the 2.5D finite/infinite element method(FIEM).Firstly,the applicability of the 2.5D FIEM for 3D seismic analysis is summarized.With the exact solutions obtained for the free field in the Appendix,the equivalent seismic forces are rationally computed for the near-field boundary,considering the horizontal and vertical excitations of the Chi-Chi Earthquake.By performing seismic analysis of the half space embedded with a tunnel using the 2.5D FIEM,the time-domain responses of the tunnel are obtained.The accuracy of the present solutions is verified against those of de Barros and Luco.Conclusions drawn from the parametric study include:(1)Stress concentration for the principal stress under oblique incident seismic waves occurs at the polar angles of 0(vault),90,180(inverted arch),and 270of the lining wall.(2)The vault and inverted arch are the weakest parts of the tunnel during earthquakes.(3)The accelerations of the tunnel during earthquakes can be regarded as of the rigid body type.(4)The responses of the tunnel lining caused by SV waves of an earthquake are much more critical than those by P waves.(5)For arbitrary seismic waves,the maximum longitudinal acceleration azmax is of the same order of magnitude as the maximum horizontal acceleration axmax.