Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis r...Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.展开更多
A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to ...A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to realize the optimal design of the butterfly-shaped linear ultrasonic motor. First, the operation principle of the motor was introduced. Second, the finite element parameterized model of the stator of the motor was built using ANSYS parametric design language and some structure parameters of the stator were selected as design variables. Third, the sample points were selected in design variable space using latin hypercube Design. Through modal analysis and harmonic response analysis of the stator based on these sample points, the target responses were obtained. These sample points and response values were combined together to build a response surface model. Finally, the simplex method was used to find the optimal solution. The experimental results showed that many aspects of the design requirements of the butterfly-shaped linear ultrasonic motor have been fulfilled. The prototype motor fabricated based on the optimal design result exhibited considerably high dynamic performance, such as no-load speed of 873 ram/s, maximal thrust of 27.5 N, maximal efficiency of 43%, and thrust-weight ratio of 45.8.展开更多
Texture evolution and inhomogeneous deformation of polycrystalline Cu during uniaxial compression are investigated at the grain scale by combining crystal plasticity finite element method(CPFEM) with particle swarm op...Texture evolution and inhomogeneous deformation of polycrystalline Cu during uniaxial compression are investigated at the grain scale by combining crystal plasticity finite element method(CPFEM) with particle swarm optimization(PSO) algorithm. The texture-based representative volume element(TBRVE) is used in the crystal plasticity finite element model, where a given number of crystallographic orientations are obtained by means of discretizing the orientation distribution function(ODF) based on electron backscattered diffraction(EBSD) experiment data. Three-dimensional grains with different morphologies are generated on the basis of Voronoi tessellation. The PSO algorithm plays a significant role in identifying the material parameters and saving computational time. The macroscopic stress–strain curve is predicted based on CPFEM, where the simulation results are in good agreement with the experimental ones. Therefore, CPFEM is a powerful candidate for capturing the texture evolution and clarifying the inhomogeneous plastic deformation of polycrystalline Cu. The simulation results indicate that the <110> fiber texture is generated finally with the progression of plastic deformation. The inhomogeneous distribution of rotation angles lays the foundation for the inhomogeneous deformation of polycrystalline Cu in terms of grain scale.展开更多
Due to the current trend towards lightweight design in automotive industry,hollow stepped gear shafts for automobile and its radial forging process are widely investigated.Utilizing coupled finite element thermo-mecha...Due to the current trend towards lightweight design in automotive industry,hollow stepped gear shafts for automobile and its radial forging process are widely investigated.Utilizing coupled finite element thermo-mechanical model,radial forging process of a hollow stepped gear shaft for automobile was simulated.The optimal combination of three process parameters including initial temperature,rotation rate and radial reduction was also selected using orthogonal design method.To examine the strain inhomogeneity of the forging workpiece,the strain inhomogeneity factor was introduced.The results reveal that the maximum effective strain and the minimum effective strain appeared in the outermost and innermost zones of different cross sections for the hollow stepped gear shaft,respectively.Optimal forging parameters are determined as a combination of initial temperature of 780°C,rotation rate of 21°/stroke and radial reduction of 3 mm.展开更多
In order to present the microstructures of dynamic recrystallization(DRX) in different deformation zones of hot extruded NiTi shape memory alloy(SMA) pipe coupling,a simulation approach combining finite element method...In order to present the microstructures of dynamic recrystallization(DRX) in different deformation zones of hot extruded NiTi shape memory alloy(SMA) pipe coupling,a simulation approach combining finite element method(FEM) with cellular automaton(CA) was developed and the relationship between the macroscopic field variables and the microscopic internal variables was established.The results show that there exists a great distinction among the microstructures in different zones of pipe coupling because deformation histories of these regions are diverse.Large plastic deformation may result in fine recrystallized grains,whereas the recrystallized grains may grow very substantially if there is a rigid translation during the deformation,even if the final plastic strain is very large.As a consequence,the deformation history has a significant influence on the evolution path of the DRX as well as the final microstructures of the DRX,including the morphology,the mean grain size and the recrystallization fraction.展开更多
The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES...The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.展开更多
A general and efficient parallel approach is proposed for the first time to parallelize the hybrid finiteelement-boundary-integral-multi-level fast multipole algorithm （FE-BI-MLFMA）. Among many algorithms of FE-BI-M...A general and efficient parallel approach is proposed for the first time to parallelize the hybrid finiteelement-boundary-integral-multi-level fast multipole algorithm （FE-BI-MLFMA）. Among many algorithms of FE-BI-MLFMA, the decomposition algorithm （DA） is chosen as a basis for the parallelization of FE-BI-MLFMA because of its distinct numerical characteristics suitable for parallelization. On the basis of the DA, the parallelization of FE-BI-MLFMA is carried out by employing the parallelized multi-frontal method for the matrix from the finiteelement method and the parallelized MLFMA for the matrix from the boundary integral method respectively. The programming and numerical experiments of the proposed parallel approach are carried out in the high perfor- mance computing platform CEMS-Liuhui. Numerical experiments demonstrate that FE-BI-MLFMA is efficiently parallelized and its computational capacity is greatly improved without losing accuracy, efficiency, and generality.展开更多
Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studi...Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studies performed with FEM requires sophisticated procedures and intensive computational effort,therefore simpler and accurate numerical approaches are always worthwhile to decrease armor development time.This study aims to apply a hybrid method using FEM simulation and artificial neural network(ANN) analysis to approximate ballistic limit thickness for armor steels.To achieve this objective,a predictive model based on the artificial neural networks is developed to determine ballistic resistance of high hardness armor steels against 7.62 mm armor piercing ammunition.In this methodology,the FEM simulations are used to create training cases for Multilayer Perceptron(MLP) three layer networks.In order to validate FE simulation methodology,ballistic shot tests on 20 mm thickness target were performed according to standard Stanag 4569.Afterwards,the successfully trained ANN(s) is used to predict the ballistic limit thickness of 500 HB high hardness steel armor.Results show that even with limited number of data,FEM-ANN approach can be used to predict ballistic penetration depth with adequate accuracy.展开更多
In engineering practice simplified methods are essential to the seismic design of embedded earth retaining walls,as fullydynamic numerical analyses are costly,time-consuming and require specific expertise.Recently dev...In engineering practice simplified methods are essential to the seismic design of embedded earth retaining walls,as fullydynamic numerical analyses are costly,time-consuming and require specific expertise.Recently developed pseudostatic methods provide earth stresses and internal forces,even in those cases in which the strength of the soil surrounding the structure is not entirely mobilised.Semiempirical correlations or Newmark sliding block method provide an estimate of earthquake-induced permanent displacements.However,the use of these methods is hindered by uncertainties in the evaluation of a few input parameters,affecting the reliability of the methods.This study uses 1 D site response analyses and 2 D fully-dynamic finite element analyses to show that simplified methods can provide a reasonable estimate of the maximum bending moment and permanent displacements for stiff cantilever walls embedded in uniform sand,providing that a few input parameters are evaluated through semiempirical correlations and a simple 1 D site response analysis.展开更多
A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element...A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.展开更多
Facing the problems lack of considering the non-uniform distribution of the static bias magnetic field and computing the panicle displacements in the simulation model of electromagnetic acoustic transducer （EMAT）, a...Facing the problems lack of considering the non-uniform distribution of the static bias magnetic field and computing the panicle displacements in the simulation model of electromagnetic acoustic transducer （EMAT）, a multi-field coupled model was established and the finite element method （FEM） was presented to calculate the entire transduction process. The multi-field coupled model included the static magnetic field, pulsed eddy current field and mechanical field. The FEM equations of the three fields were derived by Garlerkin FEM method. Thus, the entire transduction process of the EMAT was calculated through sequentially coupling the three fields. The transduction process of a Lamb wave EMAT was calculated according to the present model and method. The results show that, by the present method, it is valid to calculate the particle displacement under the given excitation signal and non-uniformly distributed static magnetic field. Calculation error will be brought about if the non-uniform distribution of the static bias magnetic field is neglected.展开更多
The aim of this work is to propose a 3D FE model of a theoretical assembling straight bevel gear pair to analyze the contact fatigue on the tooth surface and the bending fatigue in the tooth root. Based on the cumulat...The aim of this work is to propose a 3D FE model of a theoretical assembling straight bevel gear pair to analyze the contact fatigue on the tooth surface and the bending fatigue in the tooth root. Based on the cumulative fatigue criterion and the stress-life equation, the key meshing states of the gear pair were investigated for the contact fatigue and the bending fatigue. Then, the reliability of the proposed model was proved by comparing the calculation result with the simulation result. Further study was performed to analyze the variation of the contact fatigue stress and the bending fatigue stress under different loads. Furthermore, the roles of the driving pinion and the driven gear pair were evaluated in the fatigue life of the straight bevel gear pair and the main fatigue failure mode was determined for the significant gear. The results show that the fatigue failure of the driving pinion is the main fatigue failure for the straight bevel gear pair and the bending fatigue failure is the main fatigue failure for the driving pinion.展开更多
Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid...Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.展开更多
The extended finite element method(XFEM) is a numerical method for modeling discontinuities within the classical finite element framework. The computation mesh in XFEM is independent of the discontinuities, such that ...The extended finite element method(XFEM) is a numerical method for modeling discontinuities within the classical finite element framework. The computation mesh in XFEM is independent of the discontinuities, such that remeshing for moving discontinuities can be overcome. The extended finite element method is presented for hydro-mechanical modeling of impermeable discontinuities in rock. The governing equation of XFEM for hydraulic fracture modeling is derived by the virtual work principle of the fracture problem considering the water pressure on crack surface. The coupling relationship between water pressure gradient on crack surface and fracture opening width is obtained by semi-analytical and semi-numerical method. This method simplifies coupling analysis iteration and improves computational precision. Finally, the efficiency of the proposed method for modeling hydraulic fracture problems is verified by two examples and the advantages of the XFEM for hydraulic fracturing analysis are displayed.展开更多
Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Fi...Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Firstly,the single-edge notched beam test was used to analyze the temperature effect and calculate the material parameters.Then,the mechanical responses were studied using numerical analysis.It is concluded that 5℃ can be selected as the critical temperature that affects the fracture properties,and numerical simulations indicate that crack propagation is found to significantly affect the stress state of the epoxy asphalt mixture.The maximum principal stress at the crack surface exhibits different trends at various temperatures.Numerical solution of stress intensity factor can well meet the theoretical solution,especially when the temperature is lower than 5℃.展开更多
Deformation behavior of slab at the straightening stage during continuous casting was simulated by the explicit dynamic finite element method,and the stress distribution along the width direction of the slab and its c...Deformation behavior of slab at the straightening stage during continuous casting was simulated by the explicit dynamic finite element method,and the stress distribution along the width direction of the slab and its change regularity at slab center during continuous casting were obtained.The influence of distribution and change of stress on the propagation of longitudinal cracks on slab surface was discussed.The results show that the tensional stress appears on slab surface at the inner arc side and the compressive stress appears on slab surface at the outer arc side at stages 6-8 in straightening zone during continuous casting.Longitudinal cracks generally appear on slab top surface and do not appear on slab bottom surface,which are also observed in industry.展开更多
The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathe...The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathematical descriptions,namely,limit state functions of failure modes.Several problems are to be solved in the use of traditional methods for gravity dams.One is how to define the limit state function really reflecting the mechanical mechanism of the failure mode;another is how to understand the relationship among failure modes and enable the probability of the whole structure to be determined.Performing SFP analysis for a gravity dam system is a challenging task.This work proposes a novel nonlinear finite-element-based SFP analysis method for gravity dams.Firstly,reasonable nonlinear constitutive modes for dam concrete,concrete/rock interface and rock foundation are respectively introduced according to corresponding mechanical mechanisms.Meanwhile the response surface(RS) method is used to model limit state functions of main failure modes through the Monte Carlo(MC) simulation results of the dam-interface-foundation interaction finite element(FE) analysis.Secondly,a numerical SFP method is studied to compute the probabilities of several failure modes efficiently by simple matrix integration operations.Then,the nonlinear FE-based SFP analysis methodology for gravity dams considering correlated failure modes with the additional sensitivity analysis is proposed.Finally,a comprehensive computational platform for interfacing the proposed method with the open source FE code Code Aster is developed via a freely available MATLAB software tool(FERUM).This methodology is demonstrated by a case study of an existing gravity dam analysis,in which the dominant failure modes are identified,and the corresponding performance functions are established.Then,the dam failure probability of the structural system is obtained by the proposed method considering the correlation relationship of main failure modes on the basis of the mechanical mechanism analysis with the MC-FE simulations.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and st...The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.展开更多
文摘Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.
基金Projects(51275235, 50975135) supported by the National Natural Science Foundation of ChinaProject(U0934004) supported by the Natural Science Foundation of Guangdong Province, ChinaProject(2011CB707602) supported by the National Basic Research Program of China
文摘A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to realize the optimal design of the butterfly-shaped linear ultrasonic motor. First, the operation principle of the motor was introduced. Second, the finite element parameterized model of the stator of the motor was built using ANSYS parametric design language and some structure parameters of the stator were selected as design variables. Third, the sample points were selected in design variable space using latin hypercube Design. Through modal analysis and harmonic response analysis of the stator based on these sample points, the target responses were obtained. These sample points and response values were combined together to build a response surface model. Finally, the simplex method was used to find the optimal solution. The experimental results showed that many aspects of the design requirements of the butterfly-shaped linear ultrasonic motor have been fulfilled. The prototype motor fabricated based on the optimal design result exhibited considerably high dynamic performance, such as no-load speed of 873 ram/s, maximal thrust of 27.5 N, maximal efficiency of 43%, and thrust-weight ratio of 45.8.
基金Projects(51305091,51475101) supported by the National Natural Science Foundation of ChinaProject(20132304120025) supported by Specialized Research Fund for the Doctoral Program of Higher Education,China
文摘Texture evolution and inhomogeneous deformation of polycrystalline Cu during uniaxial compression are investigated at the grain scale by combining crystal plasticity finite element method(CPFEM) with particle swarm optimization(PSO) algorithm. The texture-based representative volume element(TBRVE) is used in the crystal plasticity finite element model, where a given number of crystallographic orientations are obtained by means of discretizing the orientation distribution function(ODF) based on electron backscattered diffraction(EBSD) experiment data. Three-dimensional grains with different morphologies are generated on the basis of Voronoi tessellation. The PSO algorithm plays a significant role in identifying the material parameters and saving computational time. The macroscopic stress–strain curve is predicted based on CPFEM, where the simulation results are in good agreement with the experimental ones. Therefore, CPFEM is a powerful candidate for capturing the texture evolution and clarifying the inhomogeneous plastic deformation of polycrystalline Cu. The simulation results indicate that the <110> fiber texture is generated finally with the progression of plastic deformation. The inhomogeneous distribution of rotation angles lays the foundation for the inhomogeneous deformation of polycrystalline Cu in terms of grain scale.
基金Projects(51774054,51974050)supported by the National Natural Science Foundation of China。
文摘Due to the current trend towards lightweight design in automotive industry,hollow stepped gear shafts for automobile and its radial forging process are widely investigated.Utilizing coupled finite element thermo-mechanical model,radial forging process of a hollow stepped gear shaft for automobile was simulated.The optimal combination of three process parameters including initial temperature,rotation rate and radial reduction was also selected using orthogonal design method.To examine the strain inhomogeneity of the forging workpiece,the strain inhomogeneity factor was introduced.The results reveal that the maximum effective strain and the minimum effective strain appeared in the outermost and innermost zones of different cross sections for the hollow stepped gear shaft,respectively.Optimal forging parameters are determined as a combination of initial temperature of 780°C,rotation rate of 21°/stroke and radial reduction of 3 mm.
基金Projects(51305091,51475101)supported by the National Natural Science Foundation of ChinaProject(20132304120025)supported by Specialized Research Fund for the Doctoral Program of Higher Education,China
文摘In order to present the microstructures of dynamic recrystallization(DRX) in different deformation zones of hot extruded NiTi shape memory alloy(SMA) pipe coupling,a simulation approach combining finite element method(FEM) with cellular automaton(CA) was developed and the relationship between the macroscopic field variables and the microscopic internal variables was established.The results show that there exists a great distinction among the microstructures in different zones of pipe coupling because deformation histories of these regions are diverse.Large plastic deformation may result in fine recrystallized grains,whereas the recrystallized grains may grow very substantially if there is a rigid translation during the deformation,even if the final plastic strain is very large.As a consequence,the deformation history has a significant influence on the evolution path of the DRX as well as the final microstructures of the DRX,including the morphology,the mean grain size and the recrystallization fraction.
基金funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under Grant number 107.02-2019.330。
文摘The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.
文摘A general and efficient parallel approach is proposed for the first time to parallelize the hybrid finiteelement-boundary-integral-multi-level fast multipole algorithm （FE-BI-MLFMA）. Among many algorithms of FE-BI-MLFMA, the decomposition algorithm （DA） is chosen as a basis for the parallelization of FE-BI-MLFMA because of its distinct numerical characteristics suitable for parallelization. On the basis of the DA, the parallelization of FE-BI-MLFMA is carried out by employing the parallelized multi-frontal method for the matrix from the finiteelement method and the parallelized MLFMA for the matrix from the boundary integral method respectively. The programming and numerical experiments of the proposed parallel approach are carried out in the high perfor- mance computing platform CEMS-Liuhui. Numerical experiments demonstrate that FE-BI-MLFMA is efficiently parallelized and its computational capacity is greatly improved without losing accuracy, efficiency, and generality.
基金Otokar Otomotiv ve Savunma Sanayi A.S. for the financial support
文摘Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studies performed with FEM requires sophisticated procedures and intensive computational effort,therefore simpler and accurate numerical approaches are always worthwhile to decrease armor development time.This study aims to apply a hybrid method using FEM simulation and artificial neural network(ANN) analysis to approximate ballistic limit thickness for armor steels.To achieve this objective,a predictive model based on the artificial neural networks is developed to determine ballistic resistance of high hardness armor steels against 7.62 mm armor piercing ammunition.In this methodology,the FEM simulations are used to create training cases for Multilayer Perceptron(MLP) three layer networks.In order to validate FE simulation methodology,ballistic shot tests on 20 mm thickness target were performed according to standard Stanag 4569.Afterwards,the successfully trained ANN(s) is used to predict the ballistic limit thickness of 500 HB high hardness steel armor.Results show that even with limited number of data,FEM-ANN approach can be used to predict ballistic penetration depth with adequate accuracy.
文摘In engineering practice simplified methods are essential to the seismic design of embedded earth retaining walls,as fullydynamic numerical analyses are costly,time-consuming and require specific expertise.Recently developed pseudostatic methods provide earth stresses and internal forces,even in those cases in which the strength of the soil surrounding the structure is not entirely mobilised.Semiempirical correlations or Newmark sliding block method provide an estimate of earthquake-induced permanent displacements.However,the use of these methods is hindered by uncertainties in the evaluation of a few input parameters,affecting the reliability of the methods.This study uses 1 D site response analyses and 2 D fully-dynamic finite element analyses to show that simplified methods can provide a reasonable estimate of the maximum bending moment and permanent displacements for stiff cantilever walls embedded in uniform sand,providing that a few input parameters are evaluated through semiempirical correlations and a simple 1 D site response analysis.
基金Project supported by the National Basic Research Program of China (973Project) (No.2002CB412709) and the National Natural Science Foundation of China (Nos.50278012,10272027,19832010)
文摘A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.
基金Project(10974115) supported by the National Natural Science Foundation of China
文摘Facing the problems lack of considering the non-uniform distribution of the static bias magnetic field and computing the panicle displacements in the simulation model of electromagnetic acoustic transducer （EMAT）, a multi-field coupled model was established and the finite element method （FEM） was presented to calculate the entire transduction process. The multi-field coupled model included the static magnetic field, pulsed eddy current field and mechanical field. The FEM equations of the three fields were derived by Garlerkin FEM method. Thus, the entire transduction process of the EMAT was calculated through sequentially coupling the three fields. The transduction process of a Lamb wave EMAT was calculated according to the present model and method. The results show that, by the present method, it is valid to calculate the particle displacement under the given excitation signal and non-uniformly distributed static magnetic field. Calculation error will be brought about if the non-uniform distribution of the static bias magnetic field is neglected.
基金Project(51105287) supported by the National Natural Science Foundation of ChinaProject(2012BAA08003) supported by the Key Research and Development Project of New Products and New Technologies of Hubei Province, ChinaProject(2011-P05) supported by the State Key Laboratory of Materials Processing and Die & Mould Technology, Huazhong University of Science and Technology,China
文摘The aim of this work is to propose a 3D FE model of a theoretical assembling straight bevel gear pair to analyze the contact fatigue on the tooth surface and the bending fatigue in the tooth root. Based on the cumulative fatigue criterion and the stress-life equation, the key meshing states of the gear pair were investigated for the contact fatigue and the bending fatigue. Then, the reliability of the proposed model was proved by comparing the calculation result with the simulation result. Further study was performed to analyze the variation of the contact fatigue stress and the bending fatigue stress under different loads. Furthermore, the roles of the driving pinion and the driven gear pair were evaluated in the fatigue life of the straight bevel gear pair and the main fatigue failure mode was determined for the significant gear. The results show that the fatigue failure of the driving pinion is the main fatigue failure for the straight bevel gear pair and the bending fatigue failure is the main fatigue failure for the driving pinion.
基金Projects(2006AA06Z105, 2007AA06Z134) supported by the National High-Tech Research and Development Program of ChinaProjects(2007, 2008) supported by China Scholarship Council (CSC)
文摘Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.
基金Project(2011CB013505)supported by the National Basic Research Program of ChinaProject(51279100)supported by the National Natural Science Foundation of China
文摘The extended finite element method(XFEM) is a numerical method for modeling discontinuities within the classical finite element framework. The computation mesh in XFEM is independent of the discontinuities, such that remeshing for moving discontinuities can be overcome. The extended finite element method is presented for hydro-mechanical modeling of impermeable discontinuities in rock. The governing equation of XFEM for hydraulic fracture modeling is derived by the virtual work principle of the fracture problem considering the water pressure on crack surface. The coupling relationship between water pressure gradient on crack surface and fracture opening width is obtained by semi-analytical and semi-numerical method. This method simplifies coupling analysis iteration and improves computational precision. Finally, the efficiency of the proposed method for modeling hydraulic fracture problems is verified by two examples and the advantages of the XFEM for hydraulic fracturing analysis are displayed.
基金Project(50578038)supported by the National Natural Science Foundation of China
文摘Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Firstly,the single-edge notched beam test was used to analyze the temperature effect and calculate the material parameters.Then,the mechanical responses were studied using numerical analysis.It is concluded that 5℃ can be selected as the critical temperature that affects the fracture properties,and numerical simulations indicate that crack propagation is found to significantly affect the stress state of the epoxy asphalt mixture.The maximum principal stress at the crack surface exhibits different trends at various temperatures.Numerical solution of stress intensity factor can well meet the theoretical solution,especially when the temperature is lower than 5℃.
基金Project(50634030) supported by the National Natural Science Foundation of ChinaProject(20090042120005) supported by the Doctorate Foundation of the Ministry of Education of ChinaProject(2006CB605208-1) supported by the State Basic Research Program of China
文摘Deformation behavior of slab at the straightening stage during continuous casting was simulated by the explicit dynamic finite element method,and the stress distribution along the width direction of the slab and its change regularity at slab center during continuous casting were obtained.The influence of distribution and change of stress on the propagation of longitudinal cracks on slab surface was discussed.The results show that the tensional stress appears on slab surface at the inner arc side and the compressive stress appears on slab surface at the outer arc side at stages 6-8 in straightening zone during continuous casting.Longitudinal cracks generally appear on slab top surface and do not appear on slab bottom surface,which are also observed in industry.
基金Projects(51409167,51139001,51179066)supported by the National Natural Science Foundation of ChinaProjects(201401022,201501036)supported by the Ministry of Water Resources Public Welfare Industry Research Special Fund,ChinaProjects(GG201532,GG201546)supported by the Scientific and Technological Research for Water Conservancy,Henan Province,China
文摘The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathematical descriptions,namely,limit state functions of failure modes.Several problems are to be solved in the use of traditional methods for gravity dams.One is how to define the limit state function really reflecting the mechanical mechanism of the failure mode;another is how to understand the relationship among failure modes and enable the probability of the whole structure to be determined.Performing SFP analysis for a gravity dam system is a challenging task.This work proposes a novel nonlinear finite-element-based SFP analysis method for gravity dams.Firstly,reasonable nonlinear constitutive modes for dam concrete,concrete/rock interface and rock foundation are respectively introduced according to corresponding mechanical mechanisms.Meanwhile the response surface(RS) method is used to model limit state functions of main failure modes through the Monte Carlo(MC) simulation results of the dam-interface-foundation interaction finite element(FE) analysis.Secondly,a numerical SFP method is studied to compute the probabilities of several failure modes efficiently by simple matrix integration operations.Then,the nonlinear FE-based SFP analysis methodology for gravity dams considering correlated failure modes with the additional sensitivity analysis is proposed.Finally,a comprehensive computational platform for interfacing the proposed method with the open source FE code Code Aster is developed via a freely available MATLAB software tool(FERUM).This methodology is demonstrated by a case study of an existing gravity dam analysis,in which the dominant failure modes are identified,and the corresponding performance functions are established.Then,the dam failure probability of the structural system is obtained by the proposed method considering the correlation relationship of main failure modes on the basis of the mechanical mechanism analysis with the MC-FE simulations.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Projects(40974077,41164004)supported by the National Natural Science Foundation of ChinaProject(2007AA06Z134)supported by the National High Technology Research and Development Program of China+2 种基金Projects(2011GXNSFA018003,0832263)supported by the Natural Science Foundation of Guangxi Province,ChinaProject supported by Program for Excellent Talents in Guangxi Higher Education Institution,ChinaProject supported by the Foundation of Guilin University of Technology,China
文摘The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.
