3-D rigid visco-plastic finite element method (FEM) is used in the analysisof metal forming processes, including strip and plate rolling, shape rolling, slab edging, specialstrip rolling. The shifted incomplete Choles...3-D rigid visco-plastic finite element method (FEM) is used in the analysisof metal forming processes, including strip and plate rolling, shape rolling, slab edging, specialstrip rolling. The shifted incomplete Cholesky decomposition of the stiffness matrix with thesolution of the equations for velocity increment by the conjugate gradient method is combined. Thistechnique, termed the shifted ICCG method, is then employed to solve the slab edging problem. Theperformance of this algorithm in terms of the number of iterations, friction variation, shiftedparameter psi and the results of simulation for processing parameters are analysed. Numerical testsand application of this technique verify the efficiency and stability of the shifted ICCG method inthe analysis of slab edging.展开更多
In order to explore the mineral resources buried in sea mud,it is necessary to use seabed resistivity measuring equipment,which works closer to the sediments than ordinary ship-based geophysical measuring equipment. B...In order to explore the mineral resources buried in sea mud,it is necessary to use seabed resistivity measuring equipment,which works closer to the sediments than ordinary ship-based geophysical measuring equipment. Because of the harsh environment of seafloor,high pressure and highly conductive seawater,marine magnetotelluric method developed slowly. The sea floor environment is similar to the environment of logging, According to the design of dual lateral logging equipment,a new equipment for seafloor electrical resistivity measurement is designed. Four 3D FEM models that contain resistivity abnormal targets are built to test the ability of this equipment to locate different shape of shallow buried resistivity abnormal targets in sea mud. The authors propose the method to correct the response curve while the bottom surface of this equipment is suspended or not parallel to the seafloor. The resistivity of targets can be calculated accurately.展开更多
The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Lo...The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Longmaxi black shale samples using X-ray CT, we built detailed 3D digital images of cores with porosity properties and mineral contents. Next, we used finite-element （FE） methods to derive the elastic properties of the samples. The FE method can accurately model the shale mineralogy. Particular attention is paid to the derived elastic properties and their dependence on porosity and kerogen. The elastic moduli generally decrease with increasing porosity and kerogen, and there is a critical porosity （0.75） and kerogen content （ca. ≤3%） over which the elastic moduli decrease rapidly and slowly, respectively. The derived elastic moduli of gas- and oil-saturated digital cores differ little probably because of the low porosity （4.5%） of the Longmaxi black shale. Clearly, the numerical experiments demonstrated the feasibility of combining microstructure images of shale samples with elastic moduli calculations to predict shale properties.展开更多
Springback is one of important factors influencing the forming quality of numerical control （NC） bending of thin-walled tube. In this paper, a numerical-analytic method for springback angle prediction of the process...Springback is one of important factors influencing the forming quality of numerical control （NC） bending of thin-walled tube. In this paper, a numerical-analytic method for springback angle prediction of the process was put forward. The method is based on springback angle model derived using analytic method and simulation results from three-dimensional （3D） rigid-plastic finite element method （FEM）. The method is validated through comparison with experimental results. The features of the method are as follows： （1） The method is high in efficiency because it combines advantages of rigid-plastic FEM and analytic method. （2） The method is satisfactory in accuracy, since the field variables used in the model is resulting from 3D rigid-plastic FEM solution, and the effects both of axial force and strain neutral axis shift have been included. （3） Research on multi-factor effects can be carried out using the method due to its advantage inheriting from rigid-plastic FEM. The method described here is also of general significance to other bending processes.展开更多
The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress...The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress functions. The stress fields are initially assumed by means of the Lekhnitskii stress functions under the plane strain state. Applying the principle of complementary virtual work,the coupled ordinary differential equations are obtained in which the solutions can be obtained by solving a generalized eigenvalue problem. Then an iterative procedure is established to achieve convergent stress distributions. It should be noted that the stress function based extended Kantorovich method can satisfy both the traction-free and free edge stress boundary conditions during the iterative processes. The stress components near the free edges and in the interior regions are calculated and compared with those obtained results by finite element method（FEM）. The convergent stresses have good agreements with those results obtained by three dimensional（3D） FEM. For generality, various layup configurations are considered for the numerical analysis. The results show that the proposed polynomial stress function based extended Kantorovich method is accurate and efficient in predicting the local stresses in composite laminates and computationally much more efficient than the 3D FEM.展开更多
This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method （FEM） by recursive application of the one-dimensional （1D） element energy pro...This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method （FEM） by recursive application of the one-dimensional （1D） element energy projection （EEP） technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element （FE） model is reached. This conceptual dimension-by- dimension （D-by-D） discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional （2D） and three-dimensional （3D） problems can be obtained over the domain. Numerical examples of 3D Poisson＇s equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.展开更多
Backward erosion piping is an important failure mechanism for cohesive water retaining structures which are founded on a sandy aquifer. At present, the prediction models for safety assessment are often based on 2D ass...Backward erosion piping is an important failure mechanism for cohesive water retaining structures which are founded on a sandy aquifer. At present, the prediction models for safety assessment are often based on 2D assumptions. In this work, a 3D numerical approach of the groundwater flow leading to the erosion mechanism of backward erosion piping is presented and discussed. Comparison of the 2D and 3D numerical results explicitly demonstrates the inherent 3D nature of the piping phenomenon. In addition, the influence of the seepage length is investigated and discussed for both piping initiation and piping progression. The results clearly indicate the superiority of the presented 3D numerical model compared to the established 2D approach. Moreover, the 3D numerical results enable a better understanding of the complex physical mechanism involved in backward erosion piping and thus can lead to a significant improvement in the safety assessment of water retaining structures.展开更多
基金supported by Huo Yingdong Young Teachers Foundation,Ministry of State Education of ChinaNational Natural Science Foundation of China(No.59904003).
文摘3-D rigid visco-plastic finite element method (FEM) is used in the analysisof metal forming processes, including strip and plate rolling, shape rolling, slab edging, specialstrip rolling. The shifted incomplete Cholesky decomposition of the stiffness matrix with thesolution of the equations for velocity increment by the conjugate gradient method is combined. Thistechnique, termed the shifted ICCG method, is then employed to solve the slab edging problem. Theperformance of this algorithm in terms of the number of iterations, friction variation, shiftedparameter psi and the results of simulation for processing parameters are analysed. Numerical testsand application of this technique verify the efficiency and stability of the shifted ICCG method inthe analysis of slab edging.
文摘In order to explore the mineral resources buried in sea mud,it is necessary to use seabed resistivity measuring equipment,which works closer to the sediments than ordinary ship-based geophysical measuring equipment. Because of the harsh environment of seafloor,high pressure and highly conductive seawater,marine magnetotelluric method developed slowly. The sea floor environment is similar to the environment of logging, According to the design of dual lateral logging equipment,a new equipment for seafloor electrical resistivity measurement is designed. Four 3D FEM models that contain resistivity abnormal targets are built to test the ability of this equipment to locate different shape of shallow buried resistivity abnormal targets in sea mud. The authors propose the method to correct the response curve while the bottom surface of this equipment is suspended or not parallel to the seafloor. The resistivity of targets can be calculated accurately.
基金supported by the Chinese Academy of Sciences Strategic Leading Science and Technology projects(Grant No.XDB10010400)the China Postdoctoral Science Foundation(Grant No.2015M570142)
文摘The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Longmaxi black shale samples using X-ray CT, we built detailed 3D digital images of cores with porosity properties and mineral contents. Next, we used finite-element （FE） methods to derive the elastic properties of the samples. The FE method can accurately model the shale mineralogy. Particular attention is paid to the derived elastic properties and their dependence on porosity and kerogen. The elastic moduli generally decrease with increasing porosity and kerogen, and there is a critical porosity （0.75） and kerogen content （ca. ≤3%） over which the elastic moduli decrease rapidly and slowly, respectively. The derived elastic moduli of gas- and oil-saturated digital cores differ little probably because of the low porosity （4.5%） of the Longmaxi black shale. Clearly, the numerical experiments demonstrated the feasibility of combining microstructure images of shale samples with elastic moduli calculations to predict shale properties.
基金This work was supported by the National Natural Science Foundation of China for Distinguished Young Scholars （Grant No. 50225518）the Teaching and Research Award Program for 0utstanding Young Teachers in Higher Education Institution of M0E, PRCthe Aeronautical Science Foundation of China （Grant No. 04H53057）.
文摘Springback is one of important factors influencing the forming quality of numerical control （NC） bending of thin-walled tube. In this paper, a numerical-analytic method for springback angle prediction of the process was put forward. The method is based on springback angle model derived using analytic method and simulation results from three-dimensional （3D） rigid-plastic finite element method （FEM）. The method is validated through comparison with experimental results. The features of the method are as follows： （1） The method is high in efficiency because it combines advantages of rigid-plastic FEM and analytic method. （2） The method is satisfactory in accuracy, since the field variables used in the model is resulting from 3D rigid-plastic FEM solution, and the effects both of axial force and strain neutral axis shift have been included. （3） Research on multi-factor effects can be carried out using the method due to its advantage inheriting from rigid-plastic FEM. The method described here is also of general significance to other bending processes.
基金supported by the National Natural Science Foundation of China (Grants 11372145, 11372146, and 11272161)the State Key Laboratory of Mechanics and Control of Mechanical Structures (Nanjing University of Aeronautics and astronautics) (Grant MCMS-0516Y01)+1 种基金Zhejiang Provincial Top Key Discipline of Mechanics Open Foundation (Grant xklx1601)the K. C. Wong Magna Fund through Ningbo University
文摘The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress functions. The stress fields are initially assumed by means of the Lekhnitskii stress functions under the plane strain state. Applying the principle of complementary virtual work,the coupled ordinary differential equations are obtained in which the solutions can be obtained by solving a generalized eigenvalue problem. Then an iterative procedure is established to achieve convergent stress distributions. It should be noted that the stress function based extended Kantorovich method can satisfy both the traction-free and free edge stress boundary conditions during the iterative processes. The stress components near the free edges and in the interior regions are calculated and compared with those obtained results by finite element method（FEM）. The convergent stresses have good agreements with those results obtained by three dimensional（3D） FEM. For generality, various layup configurations are considered for the numerical analysis. The results show that the proposed polynomial stress function based extended Kantorovich method is accurate and efficient in predicting the local stresses in composite laminates and computationally much more efficient than the 3D FEM.
基金supported by the National Natural Science Foundation of China(Nos.51378293 and 51078199)
文摘This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method （FEM） by recursive application of the one-dimensional （1D） element energy projection （EEP） technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element （FE） model is reached. This conceptual dimension-by- dimension （D-by-D） discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional （2D） and three-dimensional （3D） problems can be obtained over the domain. Numerical examples of 3D Poisson＇s equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.
文摘Backward erosion piping is an important failure mechanism for cohesive water retaining structures which are founded on a sandy aquifer. At present, the prediction models for safety assessment are often based on 2D assumptions. In this work, a 3D numerical approach of the groundwater flow leading to the erosion mechanism of backward erosion piping is presented and discussed. Comparison of the 2D and 3D numerical results explicitly demonstrates the inherent 3D nature of the piping phenomenon. In addition, the influence of the seepage length is investigated and discussed for both piping initiation and piping progression. The results clearly indicate the superiority of the presented 3D numerical model compared to the established 2D approach. Moreover, the 3D numerical results enable a better understanding of the complex physical mechanism involved in backward erosion piping and thus can lead to a significant improvement in the safety assessment of water retaining structures.
