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.展开更多
Dynamic responses of track structure and wave propagation in nearby ground vibration become significant when train operates on high speeds. A train-track-ground dynamic interaction analysis model based on the 2.5D fin...Dynamic responses of track structure and wave propagation in nearby ground vibration become significant when train operates on high speeds. A train-track-ground dynamic interaction analysis model based on the 2.5D finite element method is developed for the prediction of ground vibrations due to vertical track irregularities. The one-quarter car mode,1 is used to represent the train as lumped masses connected by springs. The embankment and the underlying ground are modeled by the 2.5D finite element approach to improve the computation efficiency. The Fourier transform is applied in the direction of train＇s movement to express the wave motion with a wave-number. The one-quarter car model is coupled into the global stiffness matrix describing the track-ground dynamic system with the displacement compatibility condition at the wheel-rail interface, including the irregularities on the track surface. Dynamic responses of the track and ground due to train＇s moving loads are obtained in the wave-number domain by solving the governing equation, using a conventional finite element procedure. The amplitude and wavelength are identified as two major parameters describing track irregularities. The irregularity amplitude has a direct impact on the vertical response for low-speed trains, both for short wavelength and long wavelength irregularities. Track irregularity with shorter wavelength can generate stronger track vibration both for low-speed and high-speed cases. For low-speed case, vibrations induced by track irregularities dominate far field responses. For high-speed case, the wavelength of track irregularities has very little effect on ground vibration at distances far from track center, and train＇s wheel axle weights becomes dominant.展开更多
The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, ...The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, are modeled by the finite elements, and the wave propagation properties of the far field extending to infinity are modeled by the infinite elements. One particular feature of the 2.5D approach is that it enables the computation of the three-dimensional response of the half-space, taking into account the load-moving effect, using only a two-dimensional profile. Although the 2.5D finite/infinite element approach shows a great advantage in studying the wave propagation caused by moving trains, attention should be given to the calculation aspects, such as the rules for mesh establishment, in order to avoid producing inaccurate or erroneous results. In this paper, some essential points for consideration in analysis are highlighted, along with techniques to enhance the speed of the calculations. All these observations should prove useful in making the 2.5D finite/infinite element approach an effective one.展开更多
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.展开更多
As a force-based finite element method (FEM), large increment method (LIM) has been developed in recent years. It has been shown that LIM provided prominent advantage of parallel computation with high efficiency and l...As a force-based finite element method (FEM), large increment method (LIM) has been developed in recent years. It has been shown that LIM provided prominent advantage of parallel computation with high efficiency and low time consumption for member structural system. To fully utilize its advantage in parallel computation, it is the time to extend LIM to 2D and 3D continua analysis. In this paper, a 2D finite element library with the capability of modeling arbitrary configurations is developed. Some illustrative numerical examples are solved by using the proposed library; the obtained results are compared with those obtained from both traditional displacement-based FEM and analytical solutions, which has clearly shown the advantages of LIM.展开更多
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.展开更多
A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full t...A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number （Ky） domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell＇s equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in展开更多
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.展开更多
With the application of 2.5D Woven Variable Thickness Composites(2.5DWVTC)in aviation and other fields,the issue of strength failure in this composite type has become a focal point.First,a three-step modeling approach...With the application of 2.5D Woven Variable Thickness Composites(2.5DWVTC)in aviation and other fields,the issue of strength failure in this composite type has become a focal point.First,a three-step modeling approach is proposed to rapidly construct full-scale meso-finite element models for Outer Reduction Yarn Woven Composites(ORYWC)and Inner Reduction Yarn Woven Composites(IRYWC).Then,six independent damage variables are identified:yarn fiber tension/compression,yarn matrix tension/compression,and resin matrix tension/compression.These variables are utilized to establish the constitutive equation of woven composites,considering the coupling effects of microscopic damage.Finally,combined with the Hashin failure criterion and von Mises failure criterion,the strength prediction model is implemented in ANSYS using APDL language to simulate the strength failure process of 2.5DWVTC.The results show that the predicted stiffness and strength values of various parts of ORYWC and IRYWC are in good agreement with the relevant test results.展开更多
Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is d...Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary（scattered） fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.展开更多
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.展开更多
To make better use of 2.5D C/SiC composites in industry, it is necessary to understand the mechanical properties. A finite element model＇of 2.5D composites is established, by considering the fiber undulation and the ...To make better use of 2.5D C/SiC composites in industry, it is necessary to understand the mechanical properties. A finite element model＇of 2.5D composites is established, by considering the fiber undulation and the porosity in 2.5D C/SiC composites. The fiber direction of warp is defined by cosine function to simulate the undulation of warp, and based on uniform strain assumption, analytical model of the elastic modulus and coefficient of thermal expansion （CTE） for 2.5D C/SiC composites were established by using dual- scale model. The result is found to correlate reasonably well with the predicted results and experimental results. The parametric study also demonstrates the effects of the fiber volume fraction, distance of warp yarn, and porosity in micro-scale on the mechanical properties and the coefficients of thermal expansion.展开更多
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.展开更多
The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to eff...The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to effectively invert these spectral parameters has become the focused area of the complex resistivity inversion. An optimized BP neural network (BPNN) approach based on Quantum Particle Swarm Optimization (QPSO) algorithm was presented, which was able to improve global search ability for complex resistivity multi-parameter nonlinear inversion. In the proposed method, the nonlinear weight adjustment strategy and mutation operator were used to enhance the optimization ability of QPSO algorithm. Implementation of proposed QPSO-BPNN was given, the network had 56 hidden neurons in two hidden layers (the first hidden layer has 46 neurons and the second hidden layer has 10 neurons) and it was trained on 48 datasets and tested on another 5 synthetic datasets. The training and test results show that BP neural network optimized by the QPSO algorithm performs better than the BP neural network without initial optimization on the inversion training and test models, and the mean square error distribution is better. At the same time, a double polarized anomalous bodies model was also used to verify the feasibility and effectiveness of the proposed method, the inversion results show that the QPSO-BP algorithm inversion clearly characterizes the anomalous boundaries and is closer to the values of the parameters.展开更多
基金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.
基金Project supported by the National Key Technology R&D Program of the Ministry of Science and Technology of China(No.2009BAG12A01-B12-3)the National Natural Science Foundation of China(No.51178418)the Technology Promotion Program from the Ministry of Railway of China(No.2008G005-D)
文摘Dynamic responses of track structure and wave propagation in nearby ground vibration become significant when train operates on high speeds. A train-track-ground dynamic interaction analysis model based on the 2.5D finite element method is developed for the prediction of ground vibrations due to vertical track irregularities. The one-quarter car mode,1 is used to represent the train as lumped masses connected by springs. The embankment and the underlying ground are modeled by the 2.5D finite element approach to improve the computation efficiency. The Fourier transform is applied in the direction of train＇s movement to express the wave motion with a wave-number. The one-quarter car model is coupled into the global stiffness matrix describing the track-ground dynamic system with the displacement compatibility condition at the wheel-rail interface, including the irregularities on the track surface. Dynamic responses of the track and ground due to train＇s moving loads are obtained in the wave-number domain by solving the governing equation, using a conventional finite element procedure. The amplitude and wavelength are identified as two major parameters describing track irregularities. The irregularity amplitude has a direct impact on the vertical response for low-speed trains, both for short wavelength and long wavelength irregularities. Track irregularity with shorter wavelength can generate stronger track vibration both for low-speed and high-speed cases. For low-speed case, vibrations induced by track irregularities dominate far field responses. For high-speed case, the wavelength of track irregularities has very little effect on ground vibration at distances far from track center, and train＇s wheel axle weights becomes dominant.
基金Science Council Under Grant No.NSC 89-2211-E-002-020
文摘The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, are modeled by the finite elements, and the wave propagation properties of the far field extending to infinity are modeled by the infinite elements. One particular feature of the 2.5D approach is that it enables the computation of the three-dimensional response of the half-space, taking into account the load-moving effect, using only a two-dimensional profile. Although the 2.5D finite/infinite element approach shows a great advantage in studying the wave propagation caused by moving trains, attention should be given to the calculation aspects, such as the rules for mesh establishment, in order to avoid producing inaccurate or erroneous results. In this paper, some essential points for consideration in analysis are highlighted, along with techniques to enhance the speed of the calculations. All these observations should prove useful in making the 2.5D finite/infinite element approach an effective one.
基金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.
基金the National Natural Science Foundation of China (No. 10872128)
文摘As a force-based finite element method (FEM), large increment method (LIM) has been developed in recent years. It has been shown that LIM provided prominent advantage of parallel computation with high efficiency and low time consumption for member structural system. To fully utilize its advantage in parallel computation, it is the time to extend LIM to 2D and 3D continua analysis. In this paper, a 2D finite element library with the capability of modeling arbitrary configurations is developed. Some illustrative numerical examples are solved by using the proposed library; the obtained results are compared with those obtained from both traditional displacement-based FEM and analytical solutions, which has clearly shown the advantages of LIM.
基金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.
文摘A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number （Ky） domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell＇s equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in
文摘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 National Science and Technology Major Project,China(No.2017-IV-0007-0044)National Natural Science Foundation of China(No.52175142),National Natural Science Foundation of China(No.52305170)Natural Science Foundation of Sichuan Province,China(No.2022NSFSC1885)。
文摘With the application of 2.5D Woven Variable Thickness Composites(2.5DWVTC)in aviation and other fields,the issue of strength failure in this composite type has become a focal point.First,a three-step modeling approach is proposed to rapidly construct full-scale meso-finite element models for Outer Reduction Yarn Woven Composites(ORYWC)and Inner Reduction Yarn Woven Composites(IRYWC).Then,six independent damage variables are identified:yarn fiber tension/compression,yarn matrix tension/compression,and resin matrix tension/compression.These variables are utilized to establish the constitutive equation of woven composites,considering the coupling effects of microscopic damage.Finally,combined with the Hashin failure criterion and von Mises failure criterion,the strength prediction model is implemented in ANSYS using APDL language to simulate the strength failure process of 2.5DWVTC.The results show that the predicted stiffness and strength values of various parts of ORYWC and IRYWC are in good agreement with the relevant test results.
基金supported by the Doctoral Fund Project of the Ministry of Education(No.20130061110060 class tutors)the National Natural Science Foundation of China(No.41504083)National Basic Research Program of China(973Program)(No.2013CB429805)
文摘Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary（scattered） fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.
基金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.
基金Funded by the National Basic Research Program of China,National Natural Science Foundation of China(No.51075204)Aeronautical Science Foundation of China(No.2012ZB52026)+1 种基金Research Fund for the Doctoral Program of Higher Education of China(No.20070287039)NUAA Research Funding(No.NZ2012106)
文摘To make better use of 2.5D C/SiC composites in industry, it is necessary to understand the mechanical properties. A finite element model＇of 2.5D composites is established, by considering the fiber undulation and the porosity in 2.5D C/SiC composites. The fiber direction of warp is defined by cosine function to simulate the undulation of warp, and based on uniform strain assumption, analytical model of the elastic modulus and coefficient of thermal expansion （CTE） for 2.5D C/SiC composites were established by using dual- scale model. The result is found to correlate reasonably well with the predicted results and experimental results. The parametric study also demonstrates the effects of the fiber volume fraction, distance of warp yarn, and porosity in micro-scale on the mechanical properties and the coefficients of thermal expansion.
基金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.
文摘The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to effectively invert these spectral parameters has become the focused area of the complex resistivity inversion. An optimized BP neural network (BPNN) approach based on Quantum Particle Swarm Optimization (QPSO) algorithm was presented, which was able to improve global search ability for complex resistivity multi-parameter nonlinear inversion. In the proposed method, the nonlinear weight adjustment strategy and mutation operator were used to enhance the optimization ability of QPSO algorithm. Implementation of proposed QPSO-BPNN was given, the network had 56 hidden neurons in two hidden layers (the first hidden layer has 46 neurons and the second hidden layer has 10 neurons) and it was trained on 48 datasets and tested on another 5 synthetic datasets. The training and test results show that BP neural network optimized by the QPSO algorithm performs better than the BP neural network without initial optimization on the inversion training and test models, and the mean square error distribution is better. At the same time, a double polarized anomalous bodies model was also used to verify the feasibility and effectiveness of the proposed method, the inversion results show that the QPSO-BP algorithm inversion clearly characterizes the anomalous boundaries and is closer to the values of the parameters.
