To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical c...To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.展开更多
Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with ...Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with theoretical method, the finite difference method has been verified to be feasible by a case study. It is found that under seismic loading, the dynamic response of anchorage system is synchronously fluctuated with the seismic vibration. The change of displacement amplitude of material points is slight, and comparatively speaking, the displacement amplitude of the outside point is a little larger than that of the inside point, which shows amplification effect of surface. While the axial force amplitude transforms considerably from the inside to the outside. It increases first and reaches the peak value in the intersection between the anchoring section and free section, then decreases slowly in the free section. When considering damping effect of anchorage system, the finite difference method can reflect the time attenuation characteristic better, and the calculating result would be safer and more reasonable than the dynamic steady-state theoretical method. What is more, the finite difference method can be applied to the dynamic response analysis of harmonic and seismic random vibration for all kinds of anchor, and hence has a broad application prospect.展开更多
The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of c...The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.展开更多
A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditiona...A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.展开更多
With the high-quality development of urban buildings,higher requirements are come up with for lateral bearing capacity of laterally loaded piles.Consequently,a more accurate analysis to predict the lateral response of...With the high-quality development of urban buildings,higher requirements are come up with for lateral bearing capacity of laterally loaded piles.Consequently,a more accurate analysis to predict the lateral response of the pile within an allowable displacement is an important issue.However,the current p-y curve methods cannot fully take into account the pile-soil interaction,which will lead to a large calculation difference.In this paper,a new analytical p-y curve is established and a finite difference method for determining the lateral response of pile is proposed,which can consider the separation effect of pile-soil interface and the coefficient of circumferential friction resistance.In particular,an analytical expression is developed to determine the compressive soil pressure by dividing the compressive soil pressure into two parts:initial compressive soil pressure and increment of compressive soil pressure.In addition,the relationship between compressive soil pressure and horizontal displacement of the pile is established based on the reasonable assumption.The correctness of the proposed method is verified through four examples.Based on the verified method,a parametric analysis is also conducted to investigate the influences of factors on lateral response of the pile,including internal friction angle,pile length and elastic modulus of pile.展开更多
Following the assumptions proposed by MESRI and ROKHSAR,the one-dimensional nonlinear consolidation problem of soil under constant loading is studied by introducing continuous drainage boundary.The numerical solution ...Following the assumptions proposed by MESRI and ROKHSAR,the one-dimensional nonlinear consolidation problem of soil under constant loading is studied by introducing continuous drainage boundary.The numerical solution is derived by using finite difference method and its correctness is assessed by comparing with existing analytical and numerical solutions.Based on the present solution,the effects of interface parameters,stress ratios(i.e.,final effective stress over initial effective stress,N_(σ))and the ratio c_(c)/c_(k)of compression index to permeability index on the consolidation behavior of soil are studied in detail.The results show that,the characteristics of one-dimensional nonlinear consolidation of soil are not only related to c_(c)/c_(k)and N_(σ),but also related to boundary conditions.In the engineering practice,the soil drainage rate of consolidation process can be designed by adjusting the values of interface parameters.展开更多
Based on the characteristic that the potential sliding surfaces of rock slope are commonly in the shape of either line or fold line,analysis thought of conventional pile foundation in the flat ground under complex loa...Based on the characteristic that the potential sliding surfaces of rock slope are commonly in the shape of either line or fold line,analysis thought of conventional pile foundation in the flat ground under complex load condition was applied and the upper-bound theorem of limit analysis was used to compute thrust of rock layers with all possible distribution shapes. The interaction of slope and pile was considered design load in terms of slope thrust,and the finite difference method was derived to calculate inner-force and displacement of bridge pile foundation in rock slope under complex load condition. The result of example shows that the distribution model of slope thrust has certain impact on displacement and inner-force of bridge pile foundation. The maximum displacement growth rate reaches 54% and the maximum moment and shear growth rates reach only 15% and 20%,respectively,but the trends of inner-force and displacement of bridge pile foundation are basically the same as those of the conventional pile foundation in the flat ground. When the piles bear the same level lateral thrust,the distribution shapes of slope thrust have different influence on inner-force of pile foundation,especially the rectangle distribution,and the triangle thrust has the smallest displacement and inner-force of pile foundation.展开更多
The phenomenon of heat accumulation and transportation in the composite materials is a very typical and critical issue during drilling process.In this study,a three-dimensional temperature field prediction model is pr...The phenomenon of heat accumulation and transportation in the composite materials is a very typical and critical issue during drilling process.In this study,a three-dimensional temperature field prediction model is proposed using finite difference method,based on the partly homogenization hypothesis of material,to predict temperature field in the process of drilling unidirectional carbon fiber/epoxy(C/E)composites.According to the drilling feed motion,drilling process is divided into four stages to study the temperature distributing characteristics.The results show that the temperature distribution predicted by numerical study has a good agreement with the experimental results.The temperature increases with increasing the drilling depth,and the burn phenomena is observed due to the heat accumulation,especially at the drill exit.Due to the fiber orientation,an elliptical shape of the temperature field along the direction is found for both numerical and experimental studies of C/E composites drilling process.展开更多
The pull-out capacities for soil nailing systems comprising of one single 29 mm diameter(type A) and four 16 mm diameter(type B) rebars with grouted cement were examined.A field test and numerical analysis for the typ...The pull-out capacities for soil nailing systems comprising of one single 29 mm diameter(type A) and four 16 mm diameter(type B) rebars with grouted cement were examined.A field test and numerical analysis for the type A and type B systems were carried out to investigate the pull-out capacities and the slope stability reinforcement efficiency in soil and rock slopes.The results of the pull-out tests show the mobilized shear force and load transfer characteristics with respect to soil depth.The load-displacement relationship was examined for both type A and type B systems.Slope stability analyses were carried out to study the relationships between soil and nail reinforcement and bending stiffness as well as combined axial tension and shear forces.Factors of safety were calculated in relation to the number of nails and their outside diameters.Both soil and rock slopes were included in this evaluation.展开更多
To further investigate the one-dimensional(1D)rheological consolidation mechanism of double-layered soil,the fractional derivative Merchant model(FDMM)and the non-Darcian flow model with the non-Newtonian index are re...To further investigate the one-dimensional(1D)rheological consolidation mechanism of double-layered soil,the fractional derivative Merchant model(FDMM)and the non-Darcian flow model with the non-Newtonian index are respectively introduced to describe the deformation of viscoelastic soil and the flow of pore water in the process of consolidation.Accordingly,an 1D rheological consolidation equation of double-layered soil is obtained,and its numerical analysis is performed by the implicit finite difference method.In order to verify its validity,the numerical solutions by the present method for some simplified cases are compared with the results in the related literature.Then,the influence of the revelent parameters on the rheological consolidation of double-layered soil are investigated.Numerical results indicate that the parameters of non-Darcian flow and FDMM of the first soil layer greatly influence the consolidation rate of double-layered soil.As the decrease of relative compressibility or the increase of relative permeability between the lower soil and the upper soil,the dissipation rate of excess pore water pressure and the settlement rate of the ground will be accelerated.Increasing the relative thickness of soil layer with high permeability or low compressibility will also accelerate the consolidation rate of double-layered soil.展开更多
The pile-soil system interaction computational model in liquefaction-induced lateral spreading ground was established by the finite difference numerical method.Considering an elastic-plastic subgrade reaction method,n...The pile-soil system interaction computational model in liquefaction-induced lateral spreading ground was established by the finite difference numerical method.Considering an elastic-plastic subgrade reaction method,numerical methods involving finite difference approach of pile in liquefaction-induced lateral spreading ground were derived and implemented into a finite difference program.Based on the monotonic loading tests on saturated sand after liquefaction,the liquefaction lateral deformation of the site where group piles are located was predicted.The effects of lateral ground deformation after liquefaction on a group of pile foundations were studied using the fmite difference program mentioned above,and the failure mechanism of group piles in liquefaction-induced lateral spreading ground was obtained.The applicability of the program was preliminarily verified.The results show that the bending moments at the interfaces between liquefied and non-liquefied soil layers are larger than those at the pile＇s top when the pile＇s top is embedded.The value of the additional static bending moment is larger than the peak dynamic bending moment during the earthquake,so in the pile foundation design,more than the superstructure＇s dynamics should be considered and the effect of lateral ground deformation on pile foundations cannot be neglected.展开更多
Aiming at accuracy control of the thermal crown of work rolls in cold rolling,new parameters such as regulation domain and control-efficiency factors were proposed and a numerical analysis model of the thermal crown o...Aiming at accuracy control of the thermal crown of work rolls in cold rolling,new parameters such as regulation domain and control-efficiency factors were proposed and a numerical analysis model of the thermal crown of work rolls was established using finite difference method to study roll's thermal deformation.Based on simulation results,the influences of control-efficiency factors on thermal crown are presented and the thermal crown of work rolls is analyzed after taking sub-cooling of sprinkling beam into consideration.It has been found that the control-efficiency factor of any position on the roll's surface is linear function of the temperature and the control ability of water temperature is stronger than other control parameters.In addition,the verification of the model has been carried out based on the producing technology data in some factories and the numerical simulation results coincide well with the experimental data.Therefore,this work has important value for on-line control of roll's crown in cold rolling.展开更多
Mid-infrared transmittance of submicron silver slit arrays was numerically studied with the finite difference time domain method. The slit width varies from 50 nm to 300 nm and a square feature may attach at either or...Mid-infrared transmittance of submicron silver slit arrays was numerically studied with the finite difference time domain method. The slit width varies from 50 nm to 300 nm and a square feature may attach at either or both slit sides. Although the side length of features is one or two orders of magnitude shorter than the wavelength, the attached nanoscale features can modify the transmittance significantly. The transmittance was also further investigated in detail by looking into the electromagnetic fields and Poynting vectors of selected slit geometries. The investigation results show that such change can be attributed to the cavity resonance effect inside the slit arrays. The work is of great importance to the wavelength-selective devices design in optical devices and thermal application fields.展开更多
In order to simulate the airflow in anhydrous case and the water-air flow in groundwater case, a numerical model of airflow in soil was developed. For the nonlinearity of the governing partial differential equation, t...In order to simulate the airflow in anhydrous case and the water-air flow in groundwater case, a numerical model of airflow in soil was developed. For the nonlinearity of the governing partial differential equation, the corresponding discretization and linearization methods were given. Due to the mass transfer between air-phase and water-phase, phase states of the model elements were constantly changing. Thus, parameters of the model were divided into primary ones and secondary ones, and the primary variables changing with phase states and the secondary variables can be obtained by their functional relationship with the primary variables. Additionally, the special definite condition of this numerical model was illustrated. Two examples were given to simulate the airflow in soil whether there was groundwater or not, and the effectiveness of the numerical model is verified by comparing the results of simulation with that of exoeriment.展开更多
This paper presents the one-dimensional(1D)viscoelastic consolidation system of saturated clayey soil under continuous drainage boundaries.The fractional-derivative Merchant(FDM)model has been introduced into the cons...This paper presents the one-dimensional(1D)viscoelastic consolidation system of saturated clayey soil under continuous drainage boundaries.The fractional-derivative Merchant(FDM)model has been introduced into the consolidation system to simulate the viscoelasticity.Swartzendruber’s flow law is also introduced to describe the non-Darcian flow characteristics simultaneously.The generalized numerical solution of the 1D consolidation under continuous boundaries is given by the finite difference scheme.Furthermore,to illustrate the effectiveness of the numerical method,two simplified cases are compared against the current analytical and numerical results.Finally,the effects of boundary parameters and model parameters on the viscoelastic consolidation were illustrated and discussed.The results indicated that the boundary parameters have a significant influence on consolidation.The larger the values of boundary parameters,the faster the whole dissipation of the excess pore-water pressure and soils’settlement rate.Fractional-order and viscosity parameter have little effect on consolidation,which are primarily significant in the middle and late consolidation phases.With the increase of the modulus ratio,the whole consolidation process becomes faster.Moreover,considering Swartzendruber’s flow delays the consolidation rate of the soil layer.展开更多
Thermal transport in porous media has stimulated substantial interest in engineering sciences due to increasing applications in filtration systems,porous bearings,porous layer insulation,biomechanics,geomechanics etc....Thermal transport in porous media has stimulated substantial interest in engineering sciences due to increasing applications in filtration systems,porous bearings,porous layer insulation,biomechanics,geomechanics etc.Motivated by such applications,in this article,a numerical study of entropy generation impacts on the heat and momentum transfer in time-dependent laminar incompressible boundary layer flow of a Casson viscoplastic fluid over a uniformly heated vertical cylinder embedded in a porous medium is presented.Darcy’s law is used to simulate bulk drag effects at low Reynolds number for an isotropic,homogenous porous medium.Heat line visualization is also included.The mathematical model is derived and normalized using appropriate transformation variables.The resulting non-linear time-dependent coupled governing equations with associated boundary conditions are solved via an implicit finite difference method which is efficient and unconditionally stable.The outcomes show that entropy generation and Bejan number are both elevated with increasing values of Darcy number,Casson fluid parameter,group parameter and Grashof number.To analyze the heat transfer process in a two-dimensional domain,plotting heat lines provides an excellent approach in addition to streamlines and isotherms.It is remarked that as the Darcy number increases,the deviations of heat lines from the hot wall are reduced.展开更多
A scheme for general purposed FDTD visual scientific computing software is introduced in this paper using object-oriented design (OOD) method. By abstracting the parameters of FDTD grids to an individual class and sep...A scheme for general purposed FDTD visual scientific computing software is introduced in this paper using object-oriented design (OOD) method. By abstracting the parameters of FDTD grids to an individual class and separating from the iteration procedure, the visual software can be adapted to more comprehensive computing problems. Real-time gray degree graphic and wave curve of the results can be achieved using DirectX technique. The special difference equation and data structure in dispersive medium are considered, and the peculiarity of parameters in perfectly matched layer are also discussed.展开更多
The method in which a source is set on the surface and electric potential is received in the borehole is called surface-borehole electric potential technique. Technique of surface-borehole electric potential was emplo...The method in which a source is set on the surface and electric potential is received in the borehole is called surface-borehole electric potential technique. Technique of surface-borehole electric potential was employed to study electric response of layered formation. The electric potential was obtained by solving Poisson equation with finite difference method. In the course of calculation, forward modeling wilth finite difference method was realized by adopting bandwidth non-zero storage technique and the incomplete Cholesky conjugate gradient method. The results show that method of surface-borehole can acquire anomalous electric potential corresponding tc geo-electric layers. In addition, application of appropriate mathematical operator can improve the resolution. Moreover, overburden low resistivity layers have severe influence on measuring results of surface-borehole electric potential. However, bottom low resistivity layers play a positive role in the measurement.展开更多
The radiative properties(absorptance, reflectance, and transmittance) of deep slits with five nanoscale slit profile variations at the transverse magnetic wave incidence were numerically investigated by employing the ...The radiative properties(absorptance, reflectance, and transmittance) of deep slits with five nanoscale slit profile variations at the transverse magnetic wave incidence were numerically investigated by employing the finite difference time domain method. For slits with attached features, their radiative properties can be much different due to the modified cavity geometry and dangled structures, even at wavelengths between 3 and 15 μm. The shifts of cavity resonance excitation result in higher transmittance through narrower slits at specific wavelengths and resonance modes are confirmed with the electromagnetic fields. Opposite roles possibly played by features in increasing or decreasing absorptance are determined by the feature position and demonstrated by Poynting vectors. Correlations among all properties of a representative slit array and the slit density are also comprehensively studied. When multiple slit types coexist in an array(complex slits), a wide-band transmittance or absorptance enhancement is feasible by merging spectral peaks contributed from each type of slits distinctively. Discrepancy among infrared properties of four selected slit combinations is explained while effects of slit density are also discussed.展开更多
The radiative properties of three different materials surfaces with one-dimensional microscale random roughness were obtained with the finite difference time domain method(FDTD) and near-to-far-field transformation.Th...The radiative properties of three different materials surfaces with one-dimensional microscale random roughness were obtained with the finite difference time domain method(FDTD) and near-to-far-field transformation.The surface height conforms to the Gaussian probability density function distribution.Various computational modeling issues that affect the accuracy of the predicted properties were discussed.The results show that,for perfect electric conductor(PEC) surfaces,as the surface roughness increases,the magnitude of the spike reduces and eventually the spike disappears,and also as the ratio of root mean square roughness to the surface correlation distance increases,the retroreflection becomes evident.The predicted values of FDTD solutions are in good agreement with the ray tracing and integral equation solutions.The overall trend of bidirectional reflection distribution function(BRDF) of PEC surfaces and silicon surfaces is the same,but the silicon's is much less than the former's.The BRDF difference from two polarization modes for the gold surfaces is little for smaller wavelength,but it is much larger for the longer wavelength and the FDTD simulation results agree well with the measured data.In terms of PEC surfaces,as the incident angle increases,the reflectivity becomes more specular.展开更多
基金Supported by Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)+3 种基金National Natural Science Foundation of China(12301556)Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)Basic Research Plan of Shanxi Province(202203021211129)。
文摘To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.
基金Projects(51308273,41372307,41272326) supported by the National Natural Science Foundation of ChinaProjects(2010(A)06-b) supported by Science and Technology Fund of Yunan Provincial Communication Department,China
文摘Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with theoretical method, the finite difference method has been verified to be feasible by a case study. It is found that under seismic loading, the dynamic response of anchorage system is synchronously fluctuated with the seismic vibration. The change of displacement amplitude of material points is slight, and comparatively speaking, the displacement amplitude of the outside point is a little larger than that of the inside point, which shows amplification effect of surface. While the axial force amplitude transforms considerably from the inside to the outside. It increases first and reaches the peak value in the intersection between the anchoring section and free section, then decreases slowly in the free section. When considering damping effect of anchorage system, the finite difference method can reflect the time attenuation characteristic better, and the calculating result would be safer and more reasonable than the dynamic steady-state theoretical method. What is more, the finite difference method can be applied to the dynamic response analysis of harmonic and seismic random vibration for all kinds of anchor, and hence has a broad application prospect.
文摘The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.
基金Project(Z110803)supported by the State Key Laboratory of Geomechanics and Geotechnical Engineering,ChinaProject(2008AA062303)supported by the National High Technology Research and Development Program of China
文摘A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.
基金Project(52068004)supported by the National Natural Science Foundation of ChinaProject(2018JJA160134)supported by the Natural Science Foundation of Guangxi Province,ChinaProject(AB19245018)supported by Key Research Projects of Guangxi Province,China。
文摘With the high-quality development of urban buildings,higher requirements are come up with for lateral bearing capacity of laterally loaded piles.Consequently,a more accurate analysis to predict the lateral response of the pile within an allowable displacement is an important issue.However,the current p-y curve methods cannot fully take into account the pile-soil interaction,which will lead to a large calculation difference.In this paper,a new analytical p-y curve is established and a finite difference method for determining the lateral response of pile is proposed,which can consider the separation effect of pile-soil interface and the coefficient of circumferential friction resistance.In particular,an analytical expression is developed to determine the compressive soil pressure by dividing the compressive soil pressure into two parts:initial compressive soil pressure and increment of compressive soil pressure.In addition,the relationship between compressive soil pressure and horizontal displacement of the pile is established based on the reasonable assumption.The correctness of the proposed method is verified through four examples.Based on the verified method,a parametric analysis is also conducted to investigate the influences of factors on lateral response of the pile,including internal friction angle,pile length and elastic modulus of pile.
基金Projects(51678547,41672296,51878634,51878185,41867034)supported by the National Natural Science Foundation of China。
文摘Following the assumptions proposed by MESRI and ROKHSAR,the one-dimensional nonlinear consolidation problem of soil under constant loading is studied by introducing continuous drainage boundary.The numerical solution is derived by using finite difference method and its correctness is assessed by comparing with existing analytical and numerical solutions.Based on the present solution,the effects of interface parameters,stress ratios(i.e.,final effective stress over initial effective stress,N_(σ))and the ratio c_(c)/c_(k)of compression index to permeability index on the consolidation behavior of soil are studied in detail.The results show that,the characteristics of one-dimensional nonlinear consolidation of soil are not only related to c_(c)/c_(k)and N_(σ),but also related to boundary conditions.In the engineering practice,the soil drainage rate of consolidation process can be designed by adjusting the values of interface parameters.
基金Project(50578060) supported by the National Natural Science Foundation of China
文摘Based on the characteristic that the potential sliding surfaces of rock slope are commonly in the shape of either line or fold line,analysis thought of conventional pile foundation in the flat ground under complex load condition was applied and the upper-bound theorem of limit analysis was used to compute thrust of rock layers with all possible distribution shapes. The interaction of slope and pile was considered design load in terms of slope thrust,and the finite difference method was derived to calculate inner-force and displacement of bridge pile foundation in rock slope under complex load condition. The result of example shows that the distribution model of slope thrust has certain impact on displacement and inner-force of bridge pile foundation. The maximum displacement growth rate reaches 54% and the maximum moment and shear growth rates reach only 15% and 20%,respectively,but the trends of inner-force and displacement of bridge pile foundation are basically the same as those of the conventional pile foundation in the flat ground. When the piles bear the same level lateral thrust,the distribution shapes of slope thrust have different influence on inner-force of pile foundation,especially the rectangle distribution,and the triangle thrust has the smallest displacement and inner-force of pile foundation.
基金Projects(51475073,51605076,51875079)supported by the National Natural Science Foundation of ChinaProject(2017YFB1301701)supported by the National Key Research and Development Program of China
文摘The phenomenon of heat accumulation and transportation in the composite materials is a very typical and critical issue during drilling process.In this study,a three-dimensional temperature field prediction model is proposed using finite difference method,based on the partly homogenization hypothesis of material,to predict temperature field in the process of drilling unidirectional carbon fiber/epoxy(C/E)composites.According to the drilling feed motion,drilling process is divided into four stages to study the temperature distributing characteristics.The results show that the temperature distribution predicted by numerical study has a good agreement with the experimental results.The temperature increases with increasing the drilling depth,and the burn phenomena is observed due to the heat accumulation,especially at the drill exit.Due to the fiber orientation,an elliptical shape of the temperature field along the direction is found for both numerical and experimental studies of C/E composites drilling process.
文摘The pull-out capacities for soil nailing systems comprising of one single 29 mm diameter(type A) and four 16 mm diameter(type B) rebars with grouted cement were examined.A field test and numerical analysis for the type A and type B systems were carried out to investigate the pull-out capacities and the slope stability reinforcement efficiency in soil and rock slopes.The results of the pull-out tests show the mobilized shear force and load transfer characteristics with respect to soil depth.The load-displacement relationship was examined for both type A and type B systems.Slope stability analyses were carried out to study the relationships between soil and nail reinforcement and bending stiffness as well as combined axial tension and shear forces.Factors of safety were calculated in relation to the number of nails and their outside diameters.Both soil and rock slopes were included in this evaluation.
基金Project(51578511)supported by the National Natural Science Foundation of China。
文摘To further investigate the one-dimensional(1D)rheological consolidation mechanism of double-layered soil,the fractional derivative Merchant model(FDMM)and the non-Darcian flow model with the non-Newtonian index are respectively introduced to describe the deformation of viscoelastic soil and the flow of pore water in the process of consolidation.Accordingly,an 1D rheological consolidation equation of double-layered soil is obtained,and its numerical analysis is performed by the implicit finite difference method.In order to verify its validity,the numerical solutions by the present method for some simplified cases are compared with the results in the related literature.Then,the influence of the revelent parameters on the rheological consolidation of double-layered soil are investigated.Numerical results indicate that the parameters of non-Darcian flow and FDMM of the first soil layer greatly influence the consolidation rate of double-layered soil.As the decrease of relative compressibility or the increase of relative permeability between the lower soil and the upper soil,the dissipation rate of excess pore water pressure and the settlement rate of the ground will be accelerated.Increasing the relative thickness of soil layer with high permeability or low compressibility will also accelerate the consolidation rate of double-layered soil.
基金Project(51109208)supported by the National Natural Science Foundation of ChinaProject(2013M531688)supported by the Postdoctoral Science Foundation of China+1 种基金Project(Z012009)supported by the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering(Institute of Rock and Soil Mechanics,Chinese Academy of Sciences)Project(CKSF2012054)supported by the Foundation of Changjiang River Scientific Research Institute,China
文摘The pile-soil system interaction computational model in liquefaction-induced lateral spreading ground was established by the finite difference numerical method.Considering an elastic-plastic subgrade reaction method,numerical methods involving finite difference approach of pile in liquefaction-induced lateral spreading ground were derived and implemented into a finite difference program.Based on the monotonic loading tests on saturated sand after liquefaction,the liquefaction lateral deformation of the site where group piles are located was predicted.The effects of lateral ground deformation after liquefaction on a group of pile foundations were studied using the fmite difference program mentioned above,and the failure mechanism of group piles in liquefaction-induced lateral spreading ground was obtained.The applicability of the program was preliminarily verified.The results show that the bending moments at the interfaces between liquefied and non-liquefied soil layers are larger than those at the pile＇s top when the pile＇s top is embedded.The value of the additional static bending moment is larger than the peak dynamic bending moment during the earthquake,so in the pile foundation design,more than the superstructure＇s dynamics should be considered and the effect of lateral ground deformation on pile foundations cannot be neglected.
基金Project(2007BAF02B12)supported by the National Science Technology Support Program of ChinaProjects(E2011203090,E2012203028)supported by the Natural Science Foundation of Hebei Province,China
文摘Aiming at accuracy control of the thermal crown of work rolls in cold rolling,new parameters such as regulation domain and control-efficiency factors were proposed and a numerical analysis model of the thermal crown of work rolls was established using finite difference method to study roll's thermal deformation.Based on simulation results,the influences of control-efficiency factors on thermal crown are presented and the thermal crown of work rolls is analyzed after taking sub-cooling of sprinkling beam into consideration.It has been found that the control-efficiency factor of any position on the roll's surface is linear function of the temperature and the control ability of water temperature is stronger than other control parameters.In addition,the verification of the model has been carried out based on the producing technology data in some factories and the numerical simulation results coincide well with the experimental data.Therefore,this work has important value for on-line control of roll's crown in cold rolling.
基金Project(N110402015) supported by the Fundamental Research Funds for the Central Universities of ChinaProject(2012M510075) supported by China Postdoctoral Science Foundation
文摘Mid-infrared transmittance of submicron silver slit arrays was numerically studied with the finite difference time domain method. The slit width varies from 50 nm to 300 nm and a square feature may attach at either or both slit sides. Although the side length of features is one or two orders of magnitude shorter than the wavelength, the attached nanoscale features can modify the transmittance significantly. The transmittance was also further investigated in detail by looking into the electromagnetic fields and Poynting vectors of selected slit geometries. The investigation results show that such change can be attributed to the cavity resonance effect inside the slit arrays. The work is of great importance to the wavelength-selective devices design in optical devices and thermal application fields.
基金Project(Y5080022) supported by the Natural Science Foundation of Zhejiang Province,ChinaProject(RC1202) supported by Scientific and Technological Program of Water Resources Department of Zhejiang Province in 2012,ChinaProject(Y201224384) supported by Scientific Research Program of Education Department of Zhejiang Province in 2012,China
文摘In order to simulate the airflow in anhydrous case and the water-air flow in groundwater case, a numerical model of airflow in soil was developed. For the nonlinearity of the governing partial differential equation, the corresponding discretization and linearization methods were given. Due to the mass transfer between air-phase and water-phase, phase states of the model elements were constantly changing. Thus, parameters of the model were divided into primary ones and secondary ones, and the primary variables changing with phase states and the secondary variables can be obtained by their functional relationship with the primary variables. Additionally, the special definite condition of this numerical model was illustrated. Two examples were given to simulate the airflow in soil whether there was groundwater or not, and the effectiveness of the numerical model is verified by comparing the results of simulation with that of exoeriment.
基金Projects(51879104,52078206)supported by the National Natural Science Foundation of China。
文摘This paper presents the one-dimensional(1D)viscoelastic consolidation system of saturated clayey soil under continuous drainage boundaries.The fractional-derivative Merchant(FDM)model has been introduced into the consolidation system to simulate the viscoelasticity.Swartzendruber’s flow law is also introduced to describe the non-Darcian flow characteristics simultaneously.The generalized numerical solution of the 1D consolidation under continuous boundaries is given by the finite difference scheme.Furthermore,to illustrate the effectiveness of the numerical method,two simplified cases are compared against the current analytical and numerical results.Finally,the effects of boundary parameters and model parameters on the viscoelastic consolidation were illustrated and discussed.The results indicated that the boundary parameters have a significant influence on consolidation.The larger the values of boundary parameters,the faster the whole dissipation of the excess pore-water pressure and soils’settlement rate.Fractional-order and viscosity parameter have little effect on consolidation,which are primarily significant in the middle and late consolidation phases.With the increase of the modulus ratio,the whole consolidation process becomes faster.Moreover,considering Swartzendruber’s flow delays the consolidation rate of the soil layer.
基金DST-INSPIRE (Code No. IF160028) for the grant of research fellowship
文摘Thermal transport in porous media has stimulated substantial interest in engineering sciences due to increasing applications in filtration systems,porous bearings,porous layer insulation,biomechanics,geomechanics etc.Motivated by such applications,in this article,a numerical study of entropy generation impacts on the heat and momentum transfer in time-dependent laminar incompressible boundary layer flow of a Casson viscoplastic fluid over a uniformly heated vertical cylinder embedded in a porous medium is presented.Darcy’s law is used to simulate bulk drag effects at low Reynolds number for an isotropic,homogenous porous medium.Heat line visualization is also included.The mathematical model is derived and normalized using appropriate transformation variables.The resulting non-linear time-dependent coupled governing equations with associated boundary conditions are solved via an implicit finite difference method which is efficient and unconditionally stable.The outcomes show that entropy generation and Bejan number are both elevated with increasing values of Darcy number,Casson fluid parameter,group parameter and Grashof number.To analyze the heat transfer process in a two-dimensional domain,plotting heat lines provides an excellent approach in addition to streamlines and isotherms.It is remarked that as the Darcy number increases,the deviations of heat lines from the hot wall are reduced.
基金This project was supported by the National Natural Science Foundation (No. 69831020).
文摘A scheme for general purposed FDTD visual scientific computing software is introduced in this paper using object-oriented design (OOD) method. By abstracting the parameters of FDTD grids to an individual class and separating from the iteration procedure, the visual software can be adapted to more comprehensive computing problems. Real-time gray degree graphic and wave curve of the results can be achieved using DirectX technique. The special difference equation and data structure in dispersive medium are considered, and the peculiarity of parameters in perfectly matched layer are also discussed.
基金Project supported by Global Center of Excellence in Novel Carbon Resource Sciences,Kyushu University,Japan
文摘The method in which a source is set on the surface and electric potential is received in the borehole is called surface-borehole electric potential technique. Technique of surface-borehole electric potential was employed to study electric response of layered formation. The electric potential was obtained by solving Poisson equation with finite difference method. In the course of calculation, forward modeling wilth finite difference method was realized by adopting bandwidth non-zero storage technique and the incomplete Cholesky conjugate gradient method. The results show that method of surface-borehole can acquire anomalous electric potential corresponding tc geo-electric layers. In addition, application of appropriate mathematical operator can improve the resolution. Moreover, overburden low resistivity layers have severe influence on measuring results of surface-borehole electric potential. However, bottom low resistivity layers play a positive role in the measurement.
基金Project(N130402006)supported by Fundamental Research Funds for the Central Universities,ChinaProject(51476024)supported by the National Natural Science Foundation of China
文摘The radiative properties(absorptance, reflectance, and transmittance) of deep slits with five nanoscale slit profile variations at the transverse magnetic wave incidence were numerically investigated by employing the finite difference time domain method. For slits with attached features, their radiative properties can be much different due to the modified cavity geometry and dangled structures, even at wavelengths between 3 and 15 μm. The shifts of cavity resonance excitation result in higher transmittance through narrower slits at specific wavelengths and resonance modes are confirmed with the electromagnetic fields. Opposite roles possibly played by features in increasing or decreasing absorptance are determined by the feature position and demonstrated by Poynting vectors. Correlations among all properties of a representative slit array and the slit density are also comprehensively studied. When multiple slit types coexist in an array(complex slits), a wide-band transmittance or absorptance enhancement is feasible by merging spectral peaks contributed from each type of slits distinctively. Discrepancy among infrared properties of four selected slit combinations is explained while effects of slit density are also discussed.
基金Project(2009AA05Z215) supported by the National High-Tech Research and Development Program of China
文摘The radiative properties of three different materials surfaces with one-dimensional microscale random roughness were obtained with the finite difference time domain method(FDTD) and near-to-far-field transformation.The surface height conforms to the Gaussian probability density function distribution.Various computational modeling issues that affect the accuracy of the predicted properties were discussed.The results show that,for perfect electric conductor(PEC) surfaces,as the surface roughness increases,the magnitude of the spike reduces and eventually the spike disappears,and also as the ratio of root mean square roughness to the surface correlation distance increases,the retroreflection becomes evident.The predicted values of FDTD solutions are in good agreement with the ray tracing and integral equation solutions.The overall trend of bidirectional reflection distribution function(BRDF) of PEC surfaces and silicon surfaces is the same,but the silicon's is much less than the former's.The BRDF difference from two polarization modes for the gold surfaces is little for smaller wavelength,but it is much larger for the longer wavelength and the FDTD simulation results agree well with the measured data.In terms of PEC surfaces,as the incident angle increases,the reflectivity becomes more specular.
