In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model (CZM) is developed. Integrated with the CZM, a fatigue damage evolution model is es...In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model (CZM) is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.展开更多
Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the invest...Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the investigated powder Distaloy AE were determined as functions of relative density through typical mechanical property tests of powder. The model was implemented as a user subroutine USDFLD. Single sided compaction of a d20 ram^5 mm disk green compact of Distaloy AE was simulated, and the residual stress of the disk after ejection was predicted with FEA. The FEA results of the compaction process and the residual stress of the disk show good agreement with compaction experiments and X-ray diffraction measurements, which validates the model and its parameters. The results indicate that the compressive residual stresses exist mainly in a thin layer on the side surface, but the residual stresses are very small on the top and bottom surfaces.展开更多
IEC(国际电工委员会)/CISPR(国际无线电干扰特别委员会)2003年11月发布了重组后的CISPR16(无线电骚扰和抗扰度测量设备和方法规范》(Specification for radio distur-bance and immunity measuring apparatus and methods)系列标准。该...IEC(国际电工委员会)/CISPR(国际无线电干扰特别委员会)2003年11月发布了重组后的CISPR16(无线电骚扰和抗扰度测量设备和方法规范》(Specification for radio distur-bance and immunity measuring apparatus and methods)系列标准。该系列标准目前包括以下14个子标准,它们取代以前出版的CISPR16-1、CISPR16-2、CISPR16-3、CISPR16-4及其各修正案。 第一部分:无线电骚扰和抗扰度测量设备(Radio distur-bance and immunity measuring apparatus) CISPR 16-1-1:测量设备(Measuring apparatus)展开更多
Electromagnetic forming (EMF) is a high-speed forming method which can be quite effective in increasing the forming limits of metal sheet. However, the EMF process is complicated due to magnetic-structure coupling a...Electromagnetic forming (EMF) is a high-speed forming method which can be quite effective in increasing the forming limits of metal sheet. However, the EMF process is complicated due to magnetic-structure coupling analysis. Numerical simulation offers an opportunity to overcome the problem. Nevertheless, most present models for EMF process are limited to 2D axisymmetric model. So, a three-dimensional (3D) finite element model was established to analyze the electromagnetic sheet bulging. The contact between the sheet and the die and the effect of sheet deformation on the magnetic field analysis were both taken into consideration during the forming process. The simulation results of deflection at the sheet center and 20 mm away from the center were in agreement with the experimental ones. The plastic strain energy and plastic strain were analyzed.展开更多
The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of eros...The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of erosion and deposition processes during debris flow movement, no entrainment assumption is unreasonable. The numerical model considering these processes is used for simulating debris flow. Raster grid networks of a digital elevation model in GIS provide a uniform grid system to describe complex topography. As the raster grid can be used as the finite difference mesh, the numerical model is solved numerically using the Leap-frog finite difference method. Finally, the simulation results can be displayed by GIS easily and used to debris flow evaluation. To illustrate this approach, the proposed methodology is applied to the Yohutagawa debris flow that occurred on 2oth October 2010, in Amami- Oshima area, Japan. The simulation results that reproduced the movement, erosion and deposition are in good agreement with the field investigation. The effectiveness of the dam in this real-ease is also verified by this approach. Comparison with the results were simulated by other models, shows that the present coupled model is more rational and effective.展开更多
The use of a glass-fiber reinforced composite in marine structures is becoming more common, particularly due to the potential weight savings. The mechanical response of the joint between a glass-fiber reinforced polym...The use of a glass-fiber reinforced composite in marine structures is becoming more common, particularly due to the potential weight savings. The mechanical response of the joint between a glass-fiber reinforced polymer (GRP) superstructure and a steel hull formed is examined and subsequently modified to improve performance through a combined program of modeling and testing. A finite-element model is developed to predict the response of the joint. The model takes into account the contact at the interface between different materials, progressive damage, large deformation theory, and a non-linear stress-strain relationship. To predict the progressive failure, the analysis combines Hashin failure criteria and maximum stress failure criteria. The results show stress response has a great influence on the strength and bearing of the joint. The Balsawood-steel interface is proved to be critical to the mechanical behavior of the joint. Good agreement between experimental results and numerical predictions is observed.展开更多
The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at dif...The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at different soil strengths.With the increasing capacity of numerical computers and simulation software,finite element modeling of tire/terrain interaction seems a good approach for predicting the effect of change on the parameters.In this work,an elaborated 3D model fully complianning with the geometry of radial tire 115/60R13 was established,using commercial code Solidwork Simulation.The hyper-elastic and incompressible rubber as tire main material was analyzed by Moony-Rivlin model.The Drucker-Prager yield criterion was used to model the soil compaction.Results show that the model realistically predicts the laboratory tests outputs of the modeled tire on the soft soil.展开更多
In order to investigate Joule heating power,a three-dimensional finite element model(FEM) was developed to predict Joule heating power in the channels of double-loop inductor. The simulated results were compared with ...In order to investigate Joule heating power,a three-dimensional finite element model(FEM) was developed to predict Joule heating power in the channels of double-loop inductor. The simulated results were compared with experimental data from low load trials for a 400 kW inductor. The results,such as power factor and Joule heating power,show reasonable correlation with experimental data,and Joule heating rate reaches the maximum at the corners and the minimum at the centre of the cross-section area. With increasing relative permeability of iron core,length of coils,current frequency and resistivity of metal melt,the power factor and Joule heating power change. It is concluded that current frequency,the resistivity and length of the coil play a critical role in determining the power factor and Joule heating power,whereas relative permeability of the magnetic core shows no significant influence on them.展开更多
In the present analysis, several parameters used in a numerical simulation are investigated in an integrated study to obtain their influence on the process and results of this simulation. The parameters studied are el...In the present analysis, several parameters used in a numerical simulation are investigated in an integrated study to obtain their influence on the process and results of this simulation. The parameters studied are element formulation, friction coefficient, and material model. Numerical simulations using the non-linear finite element method are conducted to produce virtual experimental data for several collision scenarios. Pattern and size damages caused by collision in a real accident case are assumed as real experimental data, and these are used to validate the method. The element model study performed indicates that the Belytschko-Tsay element formulation should be recommended for use in virtual experiments. It is recommended that the real value of the friction coefficient for materials involved is applied in simulations. For the study of the material model, the application of materials with high yield strength is recommended for use in the side hull structure.展开更多
A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume...A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid (VOF) method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio BID was shown to have significant influence on the non-dimensional transverse wave force.展开更多
We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook (BGK) collision term in the LB equation, which...We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook (BGK) collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference (FLFD) scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including (i) The thermal Couette flow; (ii) One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.展开更多
Using the finite element method and Cole-Cole model for dual-frequency IP method to research numerical simulation, the authors introduced the fundamental principle of the dual-frequency IP method and the boundary valu...Using the finite element method and Cole-Cole model for dual-frequency IP method to research numerical simulation, the authors introduced the fundamental principle of the dual-frequency IP method and the boundary value problem and variational equations, then replaced the complex resistivity of the model with the Cole-Cole model's parameters under ignoring the EM effect. Through solving the last linear equations, electric potentials of all the model's points were obtained. With changing model's parameters, the authors got different curves of the Fs and phases. According to the results of the simulation, the algorithm is proved to be correct and adaptable.展开更多
Generally,tsunami waves become hazardous only when approaching the coast.Studying the runup and inundation of tsunami waves is important for understanding the tsunami evolution and for tsunami hazard assessment.Here,w...Generally,tsunami waves become hazardous only when approaching the coast.Studying the runup and inundation of tsunami waves is important for understanding the tsunami evolution and for tsunami hazard assessment.Here,we simulated the 1993 Hokkaido-Nansei-Oki(HNO) tsunami using a finite-difference model based on nonlinear shallow-water equations.We focused on the runup and inundation of tsunami wave propagation onto coastal area of Okushiri Island near Hokkaido,Japan,and investigate the relationship of different runup heights with the morphology and bathymetry of the seashore.In the simulation,a nested 4-layer grid system and moving boundary technique are adopted to study runup and inundation.The calculated tsunami heights and inundations in the region agreed well with field measurements.The local bathymetric and topographic characteristics had a first-order effect on the runup.Numerical experiments show that the focusing of certain local bathymetric features would amplify both wave height and current velocity remarkably.The results show that computation on dense grids is necessary to reproduce the observed runup heights,and inundation velocity is an important factor preventing tsunami devastation.In addition,we discussed the potential capability of sediment transport to illustrate the impact of tsunami waves on coastal geomorphology.展开更多
The relaxation property of both Eigen model and Crow-Kimura model with a single peak fitness landscape is studied from phase transition point of view. We first analyze the eigenvalue spectra of the replication mutatio...The relaxation property of both Eigen model and Crow-Kimura model with a single peak fitness landscape is studied from phase transition point of view. We first analyze the eigenvalue spectra of the replication mutation matrices. For sufficiently long sequences, the almost crossing point between the largest and seeond-largest eigenvalues locates the error threshold at which critical slowing down behavior appears. We calculate the critical exponent in the limit of infinite sequence lengths and compare it with the result from numerical curve fittings at sufficiently long sequences. We find that for both models the relaxation time diverges with exponent 1 at the error (mutation) threshold point. Results obtained from both methods agree quite well. From the unlimited correlation length feature, the first order phase transition is further confirmed. Finally with linear stability theory, we show that the two model systems are stable for all ranges of mutation rate. The Igigen model is asymptotically stable in terms of mutant classes, and the Crow-Kimura model is completely stable.展开更多
To overcome the deficiencies of conventional geosynthetic-reinforced and pile-supported (GRPS) embankment, a new improvement technique, fixed geosynthetic technique of GRPS embankment (FGT embankment), was developed a...To overcome the deficiencies of conventional geosynthetic-reinforced and pile-supported (GRPS) embankment, a new improvement technique, fixed geosynthetic technique of GRPS embankment (FGT embankment), was developed and introduced. Based on the discussion about the load transfer mechanism of FGT embankment, a simplified check method of the requirement of geosynthetic tensile strength and a mechanical model of the FGT embankment were proposed. Two conditions, the pile cap and pile beam conditions are considered in the mechanical model. The finite difference method is used to solve the mechanical model owing to the complexity of the differential equations and the soil strata. Then, the numerical procedure is programmed. Finally, a field test is conducted to verify the mechanical model and the calculated results are in good agreement with field measured data.展开更多
A mesoscopic model has been established to investigate the thermodynamic mechanisms and densification behavior of nickel-based superalloy during additive manufacturing/three-dimensional (3D) printing (AM/3DP) by n...A mesoscopic model has been established to investigate the thermodynamic mechanisms and densification behavior of nickel-based superalloy during additive manufacturing/three-dimensional (3D) printing (AM/3DP) by numerical simulation, using a finite volume method (FVM). The influence of the applied linear energy density (LED) on dimensions of the molten pool, thermodynamic mechanisms within the pool, bubbles migration and resultant densification behavior of AM/3DP-processed superalloy has been discussed. It reveals that the center of the molten pool slightly shifts with a lagging of 4 ktm towards the center of the moving laser beam. The Mar- angoni convection, which has various flow patterns, plays a crucial role in intensifying the convective heat and mass transfer, which is responsible for the bubbles migration and densification behavior of AM/3DP-processed parts. At an optimized LED of 221.5 J/m, the outward convection favors the numerous bubbles to escape from the molten pool easily and the resultant considerably high relative density of 98.9 % is achieved. However, as the applied LED further increases over 249.5 J/m, the convection pattern is apparently intensified with the formation of vortexes and the bubbles tend to be entrapped by the rotating flow within the molten pool, resulting in a large amount of residual porosity and a sharp reduction in densification of the superalloy. The change rules of the relative density and the corresponding distribution of porosity obtained by experiments are in accordance with the simulation results.展开更多
基金The Open Research Fund of Key Laboratory of Highway Engineering of Sichuan Province of Southw est Jiaotong University (No.LHTE002201102)
文摘In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model (CZM) is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.
基金Project(2009ZX04004-031-04) supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China
文摘Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the investigated powder Distaloy AE were determined as functions of relative density through typical mechanical property tests of powder. The model was implemented as a user subroutine USDFLD. Single sided compaction of a d20 ram^5 mm disk green compact of Distaloy AE was simulated, and the residual stress of the disk after ejection was predicted with FEA. The FEA results of the compaction process and the residual stress of the disk show good agreement with compaction experiments and X-ray diffraction measurements, which validates the model and its parameters. The results indicate that the compressive residual stresses exist mainly in a thin layer on the side surface, but the residual stresses are very small on the top and bottom surfaces.
文摘IEC(国际电工委员会)/CISPR(国际无线电干扰特别委员会)2003年11月发布了重组后的CISPR16(无线电骚扰和抗扰度测量设备和方法规范》(Specification for radio distur-bance and immunity measuring apparatus and methods)系列标准。该系列标准目前包括以下14个子标准,它们取代以前出版的CISPR16-1、CISPR16-2、CISPR16-3、CISPR16-4及其各修正案。 第一部分:无线电骚扰和抗扰度测量设备(Radio distur-bance and immunity measuring apparatus) CISPR 16-1-1:测量设备(Measuring apparatus)
基金Project (50875093) supported by the National Natural Science Foundation of China
文摘Electromagnetic forming (EMF) is a high-speed forming method which can be quite effective in increasing the forming limits of metal sheet. However, the EMF process is complicated due to magnetic-structure coupling analysis. Numerical simulation offers an opportunity to overcome the problem. Nevertheless, most present models for EMF process are limited to 2D axisymmetric model. So, a three-dimensional (3D) finite element model was established to analyze the electromagnetic sheet bulging. The contact between the sheet and the die and the effect of sheet deformation on the magnetic field analysis were both taken into consideration during the forming process. The simulation results of deflection at the sheet center and 20 mm away from the center were in agreement with the experimental ones. The plastic strain energy and plastic strain were analyzed.
基金finanicial support from the Global Environment Research Fund of Japan(S-8)from Grants-in-Aid for Scientific Research(Scientific Research(B),22310113,G.Chen)from Japan Society for the Promotion of Science
文摘The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of erosion and deposition processes during debris flow movement, no entrainment assumption is unreasonable. The numerical model considering these processes is used for simulating debris flow. Raster grid networks of a digital elevation model in GIS provide a uniform grid system to describe complex topography. As the raster grid can be used as the finite difference mesh, the numerical model is solved numerically using the Leap-frog finite difference method. Finally, the simulation results can be displayed by GIS easily and used to debris flow evaluation. To illustrate this approach, the proposed methodology is applied to the Yohutagawa debris flow that occurred on 2oth October 2010, in Amami- Oshima area, Japan. The simulation results that reproduced the movement, erosion and deposition are in good agreement with the field investigation. The effectiveness of the dam in this real-ease is also verified by this approach. Comparison with the results were simulated by other models, shows that the present coupled model is more rational and effective.
基金Foundation item: Supported by the National Natural Science Foundation of China (Grant No 61004008), the Central Universities under Grant HEUCFR1001 and LBH-10138 Higher Sliding Mode Control for Underactuated Surface Ship.
文摘The use of a glass-fiber reinforced composite in marine structures is becoming more common, particularly due to the potential weight savings. The mechanical response of the joint between a glass-fiber reinforced polymer (GRP) superstructure and a steel hull formed is examined and subsequently modified to improve performance through a combined program of modeling and testing. A finite-element model is developed to predict the response of the joint. The model takes into account the contact at the interface between different materials, progressive damage, large deformation theory, and a non-linear stress-strain relationship. To predict the progressive failure, the analysis combines Hashin failure criteria and maximum stress failure criteria. The results show stress response has a great influence on the strength and bearing of the joint. The Balsawood-steel interface is proved to be critical to the mechanical behavior of the joint. Good agreement between experimental results and numerical predictions is observed.
文摘The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at different soil strengths.With the increasing capacity of numerical computers and simulation software,finite element modeling of tire/terrain interaction seems a good approach for predicting the effect of change on the parameters.In this work,an elaborated 3D model fully complianning with the geometry of radial tire 115/60R13 was established,using commercial code Solidwork Simulation.The hyper-elastic and incompressible rubber as tire main material was analyzed by Moony-Rivlin model.The Drucker-Prager yield criterion was used to model the soil compaction.Results show that the model realistically predicts the laboratory tests outputs of the modeled tire on the soft soil.
基金Project(50876116) supported by the National Natural Science Foundation of ChinaProject(2007CK3077) supported by Innovative Program of Hunan Science and Technology AgencyProject(1343-77225) supported by the Graduate School of Central South University
文摘In order to investigate Joule heating power,a three-dimensional finite element model(FEM) was developed to predict Joule heating power in the channels of double-loop inductor. The simulated results were compared with experimental data from low load trials for a 400 kW inductor. The results,such as power factor and Joule heating power,show reasonable correlation with experimental data,and Joule heating rate reaches the maximum at the corners and the minimum at the centre of the cross-section area. With increasing relative permeability of iron core,length of coils,current frequency and resistivity of metal melt,the power factor and Joule heating power change. It is concluded that current frequency,the resistivity and length of the coil play a critical role in determining the power factor and Joule heating power,whereas relative permeability of the magnetic core shows no significant influence on them.
文摘In the present analysis, several parameters used in a numerical simulation are investigated in an integrated study to obtain their influence on the process and results of this simulation. The parameters studied are element formulation, friction coefficient, and material model. Numerical simulations using the non-linear finite element method are conducted to produce virtual experimental data for several collision scenarios. Pattern and size damages caused by collision in a real accident case are assumed as real experimental data, and these are used to validate the method. The element model study performed indicates that the Belytschko-Tsay element formulation should be recommended for use in virtual experiments. It is recommended that the real value of the friction coefficient for materials involved is applied in simulations. For the study of the material model, the application of materials with high yield strength is recommended for use in the side hull structure.
基金Supported by the Science and Technology Program on Transportation Construction in Western China Ministry of Communications under Grant No.2004-328-832-51
文摘A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid (VOF) method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio BID was shown to have significant influence on the non-dimensional transverse wave force.
基金Supported by the Science Foundations of LCP and CAEP under Grant Nos. 2009A0102005 and 2009B0101012National Natural Science Foundation of China under Grant Nos. 11075021, 11074300, and 11074303+3 种基金National Basic Research Program (973 Program) under Grant No. 2007CB815105Fundamental Research Funds for the Central University under Grant No. 2010YS03Technology Support Program of LangFang under Grant Nos. 2010011029/30/31Science Foundation of NCIAE under Grant No. 2008-ky-13
文摘We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook (BGK) collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference (FLFD) scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including (i) The thermal Couette flow; (ii) One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.
基金Project supported by the National Key Technology R &D Program(2006BAB01A07)
文摘Using the finite element method and Cole-Cole model for dual-frequency IP method to research numerical simulation, the authors introduced the fundamental principle of the dual-frequency IP method and the boundary value problem and variational equations, then replaced the complex resistivity of the model with the Cole-Cole model's parameters under ignoring the EM effect. Through solving the last linear equations, electric potentials of all the model's points were obtained. With changing model's parameters, the authors got different curves of the Fs and phases. According to the results of the simulation, the algorithm is proved to be correct and adaptable.
基金Supported by the Knowledge Innovation Program of the Chinese Academy of Sciences (No.KZCX1-YW-12-01)the National Natural Science Foundation of China (No.40576014)Science and Technology Planning Project of Guangdong Province,China (No.2005D33201006)
文摘Generally,tsunami waves become hazardous only when approaching the coast.Studying the runup and inundation of tsunami waves is important for understanding the tsunami evolution and for tsunami hazard assessment.Here,we simulated the 1993 Hokkaido-Nansei-Oki(HNO) tsunami using a finite-difference model based on nonlinear shallow-water equations.We focused on the runup and inundation of tsunami wave propagation onto coastal area of Okushiri Island near Hokkaido,Japan,and investigate the relationship of different runup heights with the morphology and bathymetry of the seashore.In the simulation,a nested 4-layer grid system and moving boundary technique are adopted to study runup and inundation.The calculated tsunami heights and inundations in the region agreed well with field measurements.The local bathymetric and topographic characteristics had a first-order effect on the runup.Numerical experiments show that the focusing of certain local bathymetric features would amplify both wave height and current velocity remarkably.The results show that computation on dense grids is necessary to reproduce the observed runup heights,and inundation velocity is an important factor preventing tsunami devastation.In addition,we discussed the potential capability of sediment transport to illustrate the impact of tsunami waves on coastal geomorphology.
基金Supported in part by the National natural Science Foundation of China under Grant No.10675170Major State Basic Research Developing Program under Gant No.2007CB815003
文摘The relaxation property of both Eigen model and Crow-Kimura model with a single peak fitness landscape is studied from phase transition point of view. We first analyze the eigenvalue spectra of the replication mutation matrices. For sufficiently long sequences, the almost crossing point between the largest and seeond-largest eigenvalues locates the error threshold at which critical slowing down behavior appears. We calculate the critical exponent in the limit of infinite sequence lengths and compare it with the result from numerical curve fittings at sufficiently long sequences. We find that for both models the relaxation time diverges with exponent 1 at the error (mutation) threshold point. Results obtained from both methods agree quite well. From the unlimited correlation length feature, the first order phase transition is further confirmed. Finally with linear stability theory, we show that the two model systems are stable for all ranges of mutation rate. The Igigen model is asymptotically stable in terms of mutant classes, and the Crow-Kimura model is completely stable.
基金Project(51278216) supported by the National Natural Science Foundation of ChinaProject(20091341) supported by the Scientific Research Foundation for Returned Overseas Chinese Scholars,Ministry of Education,ChinaProject(HF-08-01-2011-240) supported by the Graduates’ Innovation Fund of Huazhong University of Science and Technology,China
文摘To overcome the deficiencies of conventional geosynthetic-reinforced and pile-supported (GRPS) embankment, a new improvement technique, fixed geosynthetic technique of GRPS embankment (FGT embankment), was developed and introduced. Based on the discussion about the load transfer mechanism of FGT embankment, a simplified check method of the requirement of geosynthetic tensile strength and a mechanical model of the FGT embankment were proposed. Two conditions, the pile cap and pile beam conditions are considered in the mechanical model. The finite difference method is used to solve the mechanical model owing to the complexity of the differential equations and the soil strata. Then, the numerical procedure is programmed. Finally, a field test is conducted to verify the mechanical model and the calculated results are in good agreement with field measured data.
基金supported by the National Natural Science Foundation of China (51575267, 51322509)the Top-Notch Young Talents Program of China+9 种基金the Outstanding Youth Foundation of Jiangsu Province of China (BK20130035)the Program for New Century Excellent Talents in University (NCET-13-0854)the Science and Technology Support Program (the Industrial Part)Jiangsu Provincial Department of Science and Technology of China (BE2014009-2)the 333 high-level talents training project (BRA2015368)the Science and Technology Foundation for Selected Overseas Chinese Scholar, Ministry of Human Resources and Social Security of Chinathe Aeronautical Science Foundation of China (2015ZE52051)the Shanghai Aerospace Science and Technology Innovation Fund (SAST2015053)the Fundamental Research Funds for the Central Universities (NE2013103, NP2015206 and NZ2016108)the Priority Academic Program Development of Jiangsu Higher Education Institutions
文摘A mesoscopic model has been established to investigate the thermodynamic mechanisms and densification behavior of nickel-based superalloy during additive manufacturing/three-dimensional (3D) printing (AM/3DP) by numerical simulation, using a finite volume method (FVM). The influence of the applied linear energy density (LED) on dimensions of the molten pool, thermodynamic mechanisms within the pool, bubbles migration and resultant densification behavior of AM/3DP-processed superalloy has been discussed. It reveals that the center of the molten pool slightly shifts with a lagging of 4 ktm towards the center of the moving laser beam. The Mar- angoni convection, which has various flow patterns, plays a crucial role in intensifying the convective heat and mass transfer, which is responsible for the bubbles migration and densification behavior of AM/3DP-processed parts. At an optimized LED of 221.5 J/m, the outward convection favors the numerous bubbles to escape from the molten pool easily and the resultant considerably high relative density of 98.9 % is achieved. However, as the applied LED further increases over 249.5 J/m, the convection pattern is apparently intensified with the formation of vortexes and the bubbles tend to be entrapped by the rotating flow within the molten pool, resulting in a large amount of residual porosity and a sharp reduction in densification of the superalloy. The change rules of the relative density and the corresponding distribution of porosity obtained by experiments are in accordance with the simulation results.