This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is e...This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.展开更多
Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass ...Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.展开更多
The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good ac...The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.展开更多
To study the influence of roll casting process parameters on temperature and thermal-stress fields for the AZ31 magnesium alloy sheets,three-dimensional geometric and 3D finite element models for roll casting were est...To study the influence of roll casting process parameters on temperature and thermal-stress fields for the AZ31 magnesium alloy sheets,three-dimensional geometric and 3D finite element models for roll casting were established based on the symmetry of roll casting by ANSYS software.Meshing method and smart-sizing algorithm were used to divide finite element mesh in ANSYS software.A series of researches on the temperature and stress distributions during solidification process with different process parameters were done by 3D finite element method.The temperatures of both the liquid-solid two-phase zone and liquid phase zone were elevated with increasing pouring temperature.With the heat transfer coefficient increasing,the two-phase region for liquid-solid becomes smaller.With the pouring temperature increasing and the increase of casting speed,the length of two-phase zone rises.The optimized of process parameters（casting speed 2 m/min,pouring temperature 640 ℃ and heat transfer coefficient 15 kW/（m2·℃） with the water pouring at roller exit was used to produce magnesium alloy AZ31 sheet,and equiaxed grains with the average grain size of 50 μm were achieved after roll casting.The simulation results give better understanding of the temperature variation in phase transformation zone and the formation mechanism of hot cracks in plates during roll casting and help to design the optimized process parameters of roll casting for Mg alloy.展开更多
Applying stiffness migration method,a 3D finite element mechanical model is established to simulate the excavation and advance processes.By using 3D nonlinear finite element method,the tunnel boring machine(TBM) excav...Applying stiffness migration method,a 3D finite element mechanical model is established to simulate the excavation and advance processes.By using 3D nonlinear finite element method,the tunnel boring machine(TBM) excavation process is dynamically simulated to analyze the stress and strain field status of surrounding rock and segment.The maximum tensile stress of segment ring caused by tunnel construction mainly lies in arch bottom and presents zonal distribution.The stress increases slightly and limitedly in the course of excavation.The maximum and minimum displacements of segment,manifesting as zonal distribution,distribute in arch bottom and vault respectively.The displacements slightly increase with the advance of TBM and gradually tend to stability.展开更多
Flexible roll forming is a promising manufacturing method for the production of variable cross section products. Considering the large plastic strain in this forming process which is much larger than that of uniform d...Flexible roll forming is a promising manufacturing method for the production of variable cross section products. Considering the large plastic strain in this forming process which is much larger than that of uniform deformation phase of uniaxial tensile test, the widely adopted method of simulating the forming processes with non-supplemented material data from uniaxial tensile test will certainly lead to large error. To reduce this error, the material data is supplemented based on three constitutive models. Then a finite element model of a six passes flexible roll forming process is established based on the supplemented material data and the original material data from the uniaxial tensile test. The flexible roll forming experiment of a B pillar reinforcing plate is carried out to verify the proposed method. Final cross section shapes of the experimental and the simulated results are compared. It is shown that the simulation calculated with supplemented material data based on Swift model agrees well with the experimental results, while the simulation based on original material data could not predict the actual deformation accurately. The results indicate that this material supplement method is reliable and indispensible, and the simulation model can well reflect the real metal forming process. Detailed analysis of the distribution and history of plastic strain at different positions are performed. A new material data supplement method is proposed to tackle the problem which is ignored in other roll forming simulations, and thus the forming process simulation accuracy can be greatly improved.展开更多
The paper focuses on the combination of the Finite Element simulation and optimization to improve process or product quality. Three different examples to illustrate the developed genetic approach are given. In all th...The paper focuses on the combination of the Finite Element simulation and optimization to improve process or product quality. Three different examples to illustrate the developed genetic approach are given. In all three examples is-DYNA3D is used to simulate the process and a general aptimiza- tion sensitivity based strategy is utilized to improve the design. The included examples are: 1) stretch bending of tubes, 2) bulging of tubes, and finally 3) hydromechanical deep drawing. these examples clearly illustrate the potential of systematic optimization in the area of metal processing.展开更多
The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying ge...The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.展开更多
In order to study the failure patterns and strength of 3D braided composites from the microscopic view, the damage propagation under tensile loading steps in three kinds of unit cells is simulated. The homogenization ...In order to study the failure patterns and strength of 3D braided composites from the microscopic view, the damage propagation under tensile loading steps in three kinds of unit cells is simulated. The homogenization formula of micro-stress and the solving approach of finite element method are given firstly. A criterion is presented to determine the damage and its pattern of each element, and then the stiffness degradation method based on Murakami's geometric damage theory is used to simulate the status of damage under tensile loading steps for three kinds of unit cells. It can be seen that the damage percentage and damage pattern of damaged unit cell are totally different for different kind of unit cells. More damaged elements are observed for face cell and corner cell than that for body cell. It is also observed that the damage firstly occurs at the area of face cell, which agrees well with experimental results. It is verified that considering the effects of face and corner cells are important for the damage and strength analysis of 3D braided composites.展开更多
In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wa...In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.展开更多
In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform ...In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.展开更多
On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate an...On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.展开更多
Laser welding (LW) becomes one of the most economical high quality joining processes. LW offers the advantage of very controlled heat input resulting in low distortion and the ability to weld heat sensitive components...Laser welding (LW) becomes one of the most economical high quality joining processes. LW offers the advantage of very controlled heat input resulting in low distortion and the ability to weld heat sensitive components. To exploit efficiently the benefits presented by LW, it is necessary to develop an integrated approach to identify and control the welding process variables in order to produce the desired weld characteristics without being forced to use the traditional and fastidious trial and error procedures. The paper presents a study of weld bead geometry characteristics prediction for laser overlap welding of low carbon galvanized steel using 3D numerical modelling and experimental validation. The temperature dependent material properties, metallurgical transformations and enthalpy method constitute the foundation of the proposed modelling approach. An adaptive 3D heat source is adopted to simulate both keyhole and conduction mode of the LW process. The simulations are performed using 3D finite element model on commercial software. The model is used to estimate the weld bead geometry characteristics for various LW parameters, such as laser power, welding speed and laser beam diameter. The calibration and validation of the 3D numerical model are based on experimental data achieved using a 3 kW Nd:Yag laser system, a structured experimental design and confirmed statistical analysis tools. The results reveal that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects on the weld quality. The results show great concordance between predicted and measured values for weld bead geometry characteristics, such as depth of penetration, bead width at the top surface and bead width at the interface between sheets, with an average accuracy greater than 95%.展开更多
The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for th...The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for the solution of 3 D frictional contact problems, the key part is the determination of the tangential slip states at the contact points, and a great amount of computing work is needed for a high accuracy result. A new method based on a combination of programming and iteration methods, which are respectively known as two main kinds of methods for contact analysis, was put forward to deal with 3 D elastic-plastic contact problems. Numerical results demonstrate the efficiency of the algorithm illustrated here.展开更多
Current developments in magnetohydrodynamic(MHD)convection and nanofluid engineering technology have have greatly enhanced heat transfer performance in process systems,particularly through the use of carbon nanotube(C...Current developments in magnetohydrodynamic(MHD)convection and nanofluid engineering technology have have greatly enhanced heat transfer performance in process systems,particularly through the use of carbon nanotube(CNT)–based fluids that offer exceptional thermal conductivity.Despite extensive research on MHD natural convection in enclosures,the combined effects of complex obstacle geometries,magnetic fields,and CNT nanofluids in three-dimensional configurations remain insufficiently explored.This research investigates MHD natural convection of carbon nanotube(CNT)-water nanofluid within a three-dimensional cavity.The study considers an inclined cross-shaped hot obstacle,a configuration not extensively explored in previous works.The work aims to elucidate the combined effects of CNT nanofluid concentration,magnetic field strength,and obstacle inclination on fluid flow patterns and heat transfer characteristics.Numerical simulations are performed using the finite element method(FEM)based on the Galerkin Weighted Residual approach.The analysis systematically considers variations in Rayleigh number(Ra),Hartmann number(Ha),nanoparticle volume fraction(Φ),and obstacle inclination angle(θ).Results show that increasing Ra from 103 to 106 enhances convective heat transfer by up to 228%,while raising the CNT volume fraction to 4.5%improves heat transfer by about 64%.In contrast,strengthening the magnetic field from Ha=0 to Ha=100 suppresses fluid motion and reduces heat transfer by nearly 67%,whereas varying the obstacle inclination from 0○to 45○leads to a 4.6%decrease in efficiency.The addition of nanoparticles slightly increases viscosity,reducing flow intensity by 8.3%when Ha=0.Furthermore,a novel multiparametric correlation is proposed,accurately predicting the average Nusselt number as a function of Ra,Ha,ϕ,andθ,with an R2 of 0.98.These findings provide new insights into the role of geometry,magnetic effects,and nanofluids in heat transfer enhancement,offering practical guidance for the design and optimization of advanced thermal systems.展开更多
Backward erosion piping is an important failure mechanism for cohesive water retaining structures which are founded on a sandy aquifer. At present, the prediction models for safety assessment are often based on 2D ass...Backward erosion piping is an important failure mechanism for cohesive water retaining structures which are founded on a sandy aquifer. At present, the prediction models for safety assessment are often based on 2D assumptions. In this work, a 3D numerical approach of the groundwater flow leading to the erosion mechanism of backward erosion piping is presented and discussed. Comparison of the 2D and 3D numerical results explicitly demonstrates the inherent 3D nature of the piping phenomenon. In addition, the influence of the seepage length is investigated and discussed for both piping initiation and piping progression. The results clearly indicate the superiority of the presented 3D numerical model compared to the established 2D approach. Moreover, the 3D numerical results enable a better understanding of the complex physical mechanism involved in backward erosion piping and thus can lead to a significant improvement in the safety assessment of water retaining structures.展开更多
The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy ...The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.51890912,51979025 and 52011530189).
文摘This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.
基金supported by the National Natural Science Foundation of China(Grant No.11304160)the Natural Science Foundation of Jiangsu Provincial Higher Education Institutions,China(Grant No.13KJB140008)the Foundation of Nanjing University of Posts and Telecommunications,China(Grant No.NY213018)
文摘Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.
文摘The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.
基金Project(CSTC 2010BB4301) supported by Natural Science Foundation Project of Chongqing,ChinaProject supported by the Open Fund for Key Laboratory of Manufacture and Test Techniques for Automobile Parts of Ministry of Education Chongqing University of Technology,2003,China
文摘To study the influence of roll casting process parameters on temperature and thermal-stress fields for the AZ31 magnesium alloy sheets,three-dimensional geometric and 3D finite element models for roll casting were established based on the symmetry of roll casting by ANSYS software.Meshing method and smart-sizing algorithm were used to divide finite element mesh in ANSYS software.A series of researches on the temperature and stress distributions during solidification process with different process parameters were done by 3D finite element method.The temperatures of both the liquid-solid two-phase zone and liquid phase zone were elevated with increasing pouring temperature.With the heat transfer coefficient increasing,the two-phase region for liquid-solid becomes smaller.With the pouring temperature increasing and the increase of casting speed,the length of two-phase zone rises.The optimized of process parameters（casting speed 2 m/min,pouring temperature 640 ℃ and heat transfer coefficient 15 kW/（m2·℃） with the water pouring at roller exit was used to produce magnesium alloy AZ31 sheet,and equiaxed grains with the average grain size of 50 μm were achieved after roll casting.The simulation results give better understanding of the temperature variation in phase transformation zone and the formation mechanism of hot cracks in plates during roll casting and help to design the optimized process parameters of roll casting for Mg alloy.
基金Supported by National Natural Science Foundation of China(No.90815019)National Key Basic Research Program of China("973" Program,No.2007CB714101)Key Project in the National Science and Technology Pillar Program during the Eleventh Five-Year Plan Period(No.2006BAB04A13)
文摘Applying stiffness migration method,a 3D finite element mechanical model is established to simulate the excavation and advance processes.By using 3D nonlinear finite element method,the tunnel boring machine(TBM) excavation process is dynamically simulated to analyze the stress and strain field status of surrounding rock and segment.The maximum tensile stress of segment ring caused by tunnel construction mainly lies in arch bottom and presents zonal distribution.The stress increases slightly and limitedly in the course of excavation.The maximum and minimum displacements of segment,manifesting as zonal distribution,distribute in arch bottom and vault respectively.The displacements slightly increase with the advance of TBM and gradually tend to stability.
基金Supported by National Natural Science Foundation of China(Grant Nos.51205004,51475003)Beijing Municipal Natural Science Foundation of China(Grant No.3152010)Beijing Municipal Education Committee Science and Technology Program,China(Grant No.KM201510009004)
文摘Flexible roll forming is a promising manufacturing method for the production of variable cross section products. Considering the large plastic strain in this forming process which is much larger than that of uniform deformation phase of uniaxial tensile test, the widely adopted method of simulating the forming processes with non-supplemented material data from uniaxial tensile test will certainly lead to large error. To reduce this error, the material data is supplemented based on three constitutive models. Then a finite element model of a six passes flexible roll forming process is established based on the supplemented material data and the original material data from the uniaxial tensile test. The flexible roll forming experiment of a B pillar reinforcing plate is carried out to verify the proposed method. Final cross section shapes of the experimental and the simulated results are compared. It is shown that the simulation calculated with supplemented material data based on Swift model agrees well with the experimental results, while the simulation based on original material data could not predict the actual deformation accurately. The results indicate that this material supplement method is reliable and indispensible, and the simulation model can well reflect the real metal forming process. Detailed analysis of the distribution and history of plastic strain at different positions are performed. A new material data supplement method is proposed to tackle the problem which is ignored in other roll forming simulations, and thus the forming process simulation accuracy can be greatly improved.
文摘The paper focuses on the combination of the Finite Element simulation and optimization to improve process or product quality. Three different examples to illustrate the developed genetic approach are given. In all three examples is-DYNA3D is used to simulate the process and a general aptimiza- tion sensitivity based strategy is utilized to improve the design. The included examples are: 1) stretch bending of tubes, 2) bulging of tubes, and finally 3) hydromechanical deep drawing. these examples clearly illustrate the potential of systematic optimization in the area of metal processing.
基金supported by a grant from the Research Grant Council of Hong Kong Special Administrative Region(Project No.11207724).
文摘The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.
基金National Natural Science Foundation of China (10772115)
文摘In order to study the failure patterns and strength of 3D braided composites from the microscopic view, the damage propagation under tensile loading steps in three kinds of unit cells is simulated. The homogenization formula of micro-stress and the solving approach of finite element method are given firstly. A criterion is presented to determine the damage and its pattern of each element, and then the stiffness degradation method based on Murakami's geometric damage theory is used to simulate the status of damage under tensile loading steps for three kinds of unit cells. It can be seen that the damage percentage and damage pattern of damaged unit cell are totally different for different kind of unit cells. More damaged elements are observed for face cell and corner cell than that for body cell. It is also observed that the damage firstly occurs at the area of face cell, which agrees well with experimental results. It is verified that considering the effects of face and corner cells are important for the damage and strength analysis of 3D braided composites.
文摘In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.
基金China Postdoctoral Science Foundation Under Grant No.20100480321National Basic Research Program of China Under Grant No. 2007CB714200
文摘In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.
基金supported by the Key Projects of the National Science Foundation of China(Nos.52178190,52078250,11832013)
文摘On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.
文摘Laser welding (LW) becomes one of the most economical high quality joining processes. LW offers the advantage of very controlled heat input resulting in low distortion and the ability to weld heat sensitive components. To exploit efficiently the benefits presented by LW, it is necessary to develop an integrated approach to identify and control the welding process variables in order to produce the desired weld characteristics without being forced to use the traditional and fastidious trial and error procedures. The paper presents a study of weld bead geometry characteristics prediction for laser overlap welding of low carbon galvanized steel using 3D numerical modelling and experimental validation. The temperature dependent material properties, metallurgical transformations and enthalpy method constitute the foundation of the proposed modelling approach. An adaptive 3D heat source is adopted to simulate both keyhole and conduction mode of the LW process. The simulations are performed using 3D finite element model on commercial software. The model is used to estimate the weld bead geometry characteristics for various LW parameters, such as laser power, welding speed and laser beam diameter. The calibration and validation of the 3D numerical model are based on experimental data achieved using a 3 kW Nd:Yag laser system, a structured experimental design and confirmed statistical analysis tools. The results reveal that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects on the weld quality. The results show great concordance between predicted and measured values for weld bead geometry characteristics, such as depth of penetration, bead width at the top surface and bead width at the interface between sheets, with an average accuracy greater than 95%.
基金theNationalKeyBasicResearchSpecialFoundation (G1 9990 3 2 80 5 ) the FoundationforUniversityKeyTeacherbytheMinistryofEducationo
文摘The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for the solution of 3 D frictional contact problems, the key part is the determination of the tangential slip states at the contact points, and a great amount of computing work is needed for a high accuracy result. A new method based on a combination of programming and iteration methods, which are respectively known as two main kinds of methods for contact analysis, was put forward to deal with 3 D elastic-plastic contact problems. Numerical results demonstrate the efficiency of the algorithm illustrated here.
基金Deputyship for Research&Innovation,Ministry of Education in Saudi Arabia for funding this research work through the project number RI-44-0451.
文摘Current developments in magnetohydrodynamic(MHD)convection and nanofluid engineering technology have have greatly enhanced heat transfer performance in process systems,particularly through the use of carbon nanotube(CNT)–based fluids that offer exceptional thermal conductivity.Despite extensive research on MHD natural convection in enclosures,the combined effects of complex obstacle geometries,magnetic fields,and CNT nanofluids in three-dimensional configurations remain insufficiently explored.This research investigates MHD natural convection of carbon nanotube(CNT)-water nanofluid within a three-dimensional cavity.The study considers an inclined cross-shaped hot obstacle,a configuration not extensively explored in previous works.The work aims to elucidate the combined effects of CNT nanofluid concentration,magnetic field strength,and obstacle inclination on fluid flow patterns and heat transfer characteristics.Numerical simulations are performed using the finite element method(FEM)based on the Galerkin Weighted Residual approach.The analysis systematically considers variations in Rayleigh number(Ra),Hartmann number(Ha),nanoparticle volume fraction(Φ),and obstacle inclination angle(θ).Results show that increasing Ra from 103 to 106 enhances convective heat transfer by up to 228%,while raising the CNT volume fraction to 4.5%improves heat transfer by about 64%.In contrast,strengthening the magnetic field from Ha=0 to Ha=100 suppresses fluid motion and reduces heat transfer by nearly 67%,whereas varying the obstacle inclination from 0○to 45○leads to a 4.6%decrease in efficiency.The addition of nanoparticles slightly increases viscosity,reducing flow intensity by 8.3%when Ha=0.Furthermore,a novel multiparametric correlation is proposed,accurately predicting the average Nusselt number as a function of Ra,Ha,ϕ,andθ,with an R2 of 0.98.These findings provide new insights into the role of geometry,magnetic effects,and nanofluids in heat transfer enhancement,offering practical guidance for the design and optimization of advanced thermal systems.
文摘Backward erosion piping is an important failure mechanism for cohesive water retaining structures which are founded on a sandy aquifer. At present, the prediction models for safety assessment are often based on 2D assumptions. In this work, a 3D numerical approach of the groundwater flow leading to the erosion mechanism of backward erosion piping is presented and discussed. Comparison of the 2D and 3D numerical results explicitly demonstrates the inherent 3D nature of the piping phenomenon. In addition, the influence of the seepage length is investigated and discussed for both piping initiation and piping progression. The results clearly indicate the superiority of the presented 3D numerical model compared to the established 2D approach. Moreover, the 3D numerical results enable a better understanding of the complex physical mechanism involved in backward erosion piping and thus can lead to a significant improvement in the safety assessment of water retaining structures.
基金Projects(41074085,41374118)supported by the National Natural Science Foundation of ChinaProject(20120162110015)supported by Doctoral Fund of Ministry of Education of ChinaProject(NCET-12-0551)supported by Program for New Century Excellent Talents in University,China
文摘The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.
