A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements f...A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.展开更多
In order to investigate the permanent deformation behavior of asphalt mixtures from discontinuity, the virtual rutting test of asphalt mixtures is developed by the discrete element method (DEM). A digital specimen g...In order to investigate the permanent deformation behavior of asphalt mixtures from discontinuity, the virtual rutting test of asphalt mixtures is developed by the discrete element method (DEM). A digital specimen generation procedure considering aggregate gradation and irregular shape is developed based on the probability theory and the Monte Carlo method. The virtual rutting test is then conducted based on the generated digital specimen. In addition, on the basis of the time-temperature superposition (TTS) principle, a calculation method is used to reduce the computation time of the virtual rutting test. The simulation results are compared with the laboratory measurements. The results show that the calculation method based on the TFS principle in the discrete element (DE) viscoelastic model can significantly reduce the computation time. The deformation law of asphalt mixtures in the virtual rutting test is similar to the laboratory measurements, and the deformation and the dynamic stability of the virtual rutting test are slightly greater than the laboratory measurements. The two-dimensional virtual rutting test can predict the permanent deformation performance of asphalt mixtures.展开更多
In order to study fracture behaviors of asphalt mixtures, virtual tests of the two-dimensional(2D) microstructure based on the discrete element method( DEM) are designed. The virtual structure of the 2D digital sp...In order to study fracture behaviors of asphalt mixtures, virtual tests of the two-dimensional(2D) microstructure based on the discrete element method( DEM) are designed. The virtual structure of the 2D digital specimen of asphalt mixture is generated based on a particle generation program, in which the gradation and the irregular shapes of aggregates are considered. With the 2D digital specimens, a DEM-based mixture model is established and center-point beam fracture simulation tests are conducted by the DEM. Meanwhile, a series of calibration tests are carried out in laboratory to evaluate the DEM model and validate the methods of virtual fracture tests. The test results indicate that the fracture intensity of asphalt mixtures predicted by the DEM matches very well with the intensity obtained in laboratory. It is concluded that the microstructural virtual tests can be used as a supplemental tool to evaluate fracture properties of asphalt mixtures.展开更多
The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete elem...The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.展开更多
Avirtual wall thicknessmethod is developed to simulate the temperature field of turbine bladeswith thermal barrier coatings(TBCs),to simplify the modeling process and improve the calculation efficiency.The results sho...Avirtual wall thicknessmethod is developed to simulate the temperature field of turbine bladeswith thermal barrier coatings(TBCs),to simplify the modeling process and improve the calculation efficiency.The results show that the virtualwall thickness method can improve themesh quality by 20%,reduce the number ofmeshes by 76.7%and save the calculation time by 35.5%,compared with the traditional real wall thickness method.The average calculation error of the two methods is between 0.21%and 0.93%.Furthermore,the temperature at the blade leading edge is the highest and the average temperature of the blade pressure surface is higher than that of the suction surface under a certain service condition.The blade surface temperature presents a high temperature at both ends and a low temperature in themiddle height when the temperature of incoming gas is uniformand constant.The thermal insulation effect of TBCs is the worst near the air film hole,and the best at the blade leading edge.According to the calculated temperature field of the substrate-coating system,the highest thermal insulation temperature of the TC layer is 172.01 K,and the thermal insulation proportions of TC,TGO and BC are 93.55%,1.54%and 4.91%,respectively.展开更多
A novel polygonal finite element method (PFEM) based on partition of unity is proposed, termed the virtual node method (VNM). To test the performance of the present method, numerical examples are given for solid m...A novel polygonal finite element method (PFEM) based on partition of unity is proposed, termed the virtual node method (VNM). To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.展开更多
In this paper, we discuss virtual element method (VEM) approximation of optimal control problem governed by Brinkman equations with control constraints. Based on the polynomial projections and variational discretizati...In this paper, we discuss virtual element method (VEM) approximation of optimal control problem governed by Brinkman equations with control constraints. Based on the polynomial projections and variational discretization of the control variable, we build up the virtual element discrete scheme of the optimal control problem and derive the discrete first order optimality system. A priori error estimates for the state, adjoint state and control variables in L<sup>2</sup> and H<sup>1</sup> norm are derived. The theoretical findings are illustrated by the numerical experiments.展开更多
This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional ...This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional boundary element method,and can be treated as a potential approach for solving the inverse problems of the heat conduction owing to the boundary-only discretization and semi-analytical algorithm.When the VBEM is applied to the inverse problems,the numerical instability may occur if a virtual boundary is not properly chosen.The method encounters a highly illconditioned matrix for the larger distance between the physical boundary and the virtual boundary,and otherwise is hard to avoid the singularity of the source point.Thus,it must adopt an appropriate regularization method to deal with the ill-posed systems of inverse problems.In this study,the VBEM and different regularization techniques are combined to model the inverse problem of three-dimensional heat conduction in orthotropic media.The proper regularization techniques not only make the virtual boundary to be allocated freer,but also solve the ill-conditioned equation of the inverse problem.Numerical examples demonstrate that the proposed method is efficient,accurate and numerically stable for solving the inverse problems of three-dimensional heat conduction in orthotropic media.展开更多
The second-order serendipity virtual element method is studied for the semilinear pseudo-parabolic equations on curved domains in this paper.Nonhomogeneous Dirichlet boundary conditions are taken into account,the exis...The second-order serendipity virtual element method is studied for the semilinear pseudo-parabolic equations on curved domains in this paper.Nonhomogeneous Dirichlet boundary conditions are taken into account,the existence and uniqueness are investigated for the weak solution of the nonhomogeneous initial-boundary value problem.The Nitschebased projection method is adopted to impose the boundary conditions in a weak way.The interpolation operator is used to deal with the nonlinear term.The Crank-Nicolson scheme is employed to discretize the temporal variable.There are two main features of the proposed scheme:(i)the internal degrees of freedom are avoided no matter what type of mesh is utilized,and(ii)the Jacobian is simple to calculate when Newton’s iteration method is applied to solve the fully discrete scheme.The error estimates are established for the discrete schemes and the theoretical results are illustrated through some numerical examples.展开更多
The virtual laminated element method (VLEM) can resolve structural shap e optimization problems with a new method. According to the characteristics of V LEM , only some characterized layer thickness values need be def...The virtual laminated element method (VLEM) can resolve structural shap e optimization problems with a new method. According to the characteristics of V LEM , only some characterized layer thickness values need be defined as design v ariables instead of boundary node coordinates or some other parameters determini ng the system boundary. One of the important features of this method is that it is not necessary to regenerate the FE(finite element) grid during the optimizati on process so as to avoid optimization failures resulting from some distortion grid elements. Th e thickness distribution in thin plate optimization problems in other studies be fore is of stepped shape. However, in this paper, a continuous thickness distrib ution can be obtained after optimization using VLEM, and is more reasonable. Fur thermore, an approximate reanalysis method named ″behavior model technique″ ca n be used to reduce the amount of structural reanalysis. Some typical examples are offered to prove the effectiveness and practicality of the proposed method.展开更多
The virtual element method(VEM)can be seen as an extension of the classical finite element method(FEM)based on Galerkin projection.It allows meshes with highly irregular shaped elements,including concave shapes.So far...The virtual element method(VEM)can be seen as an extension of the classical finite element method(FEM)based on Galerkin projection.It allows meshes with highly irregular shaped elements,including concave shapes.So far the virtual element method has been applied to various engineering problems such as elasto-plasticity,multiphysics,damage and fracture mechanics.This work focuses on the extension of the virtual element method to efficient modeling of nonlinear elasto-dynamics undergoing large deformations.Within this framework,we employ low-order ansatz functions in two and three dimensions for elements that can have arbitrary polygonal shape.The formulations considered in this contribution are based on minimization of potential function for both the static and the dynamic behavior.Generally the construction of a virtual element is based on a projection part and a stabilization part.While the stiffness matrix needs a suitable stabilization,the mass matrix can be calculated using only the projection part.For the implicit time integration scheme,Newmark-Method is used.To show the performance of the method,various two-and three-dimensional numerical examples in are presented.展开更多
应变-旋转(Strain-Rotation,S-R)和分解定理为分析几何非线性问题提供了合理可靠的理论基础,但用有限元求解时会遇到大变形发生后的网格畸变问题。近年提出的虚单元法(Virtual element method,VEM)适用于一般的多边形网格,因此,该文尝...应变-旋转(Strain-Rotation,S-R)和分解定理为分析几何非线性问题提供了合理可靠的理论基础,但用有限元求解时会遇到大变形发生后的网格畸变问题。近年提出的虚单元法(Virtual element method,VEM)适用于一般的多边形网格,因此,该文尝试使用一阶虚单元求解基于S-R和分解定理的二维几何非线性问题,以克服网格畸变的影响。基于重新定义的多项式位移空间基函数,推演获得一阶虚单元分析线弹性力学问题时允许位移空间向多项式位移空间的投影表达式;按照虚单元法双线性格式的计算规则,分析处理基于更新拖带坐标法和势能率原理的增量变分方程;进而建立离散系统方程及其矩阵表达形式,并编制MATLAB求解程序;采用常规多边形网格和畸变网格,应用该文算法分析均布荷载下的悬臂梁和均匀内压下的厚壁圆筒变形。结果与已有文献和ANSYS软件的对比表明:该文算法在两种网格中均可有效执行且具备足够数值精度。总体该文算法为基于S-R和分解定理的二维几何非线性问题求解提供了一种鲁棒方法。展开更多
基金supported by the National Natural Science Foundation of China(Nos.11271273 and 11271298)
文摘A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.
基金The National Natural Science Foundation of China(No.51108081)
文摘In order to investigate the permanent deformation behavior of asphalt mixtures from discontinuity, the virtual rutting test of asphalt mixtures is developed by the discrete element method (DEM). A digital specimen generation procedure considering aggregate gradation and irregular shape is developed based on the probability theory and the Monte Carlo method. The virtual rutting test is then conducted based on the generated digital specimen. In addition, on the basis of the time-temperature superposition (TTS) principle, a calculation method is used to reduce the computation time of the virtual rutting test. The simulation results are compared with the laboratory measurements. The results show that the calculation method based on the TFS principle in the discrete element (DE) viscoelastic model can significantly reduce the computation time. The deformation law of asphalt mixtures in the virtual rutting test is similar to the laboratory measurements, and the deformation and the dynamic stability of the virtual rutting test are slightly greater than the laboratory measurements. The two-dimensional virtual rutting test can predict the permanent deformation performance of asphalt mixtures.
基金The National High Technology Research and Development Program of China(863 Program)(No.2006AA11Z110)the Researchand Innovation Foundation for Graduate Students in Jiangsu Province(No.CX07B-156Z)the Excellent Doctoral Dissertation Foundation of Southeast University
文摘In order to study fracture behaviors of asphalt mixtures, virtual tests of the two-dimensional(2D) microstructure based on the discrete element method( DEM) are designed. The virtual structure of the 2D digital specimen of asphalt mixture is generated based on a particle generation program, in which the gradation and the irregular shapes of aggregates are considered. With the 2D digital specimens, a DEM-based mixture model is established and center-point beam fracture simulation tests are conducted by the DEM. Meanwhile, a series of calibration tests are carried out in laboratory to evaluate the DEM model and validate the methods of virtual fracture tests. The test results indicate that the fracture intensity of asphalt mixtures predicted by the DEM matches very well with the intensity obtained in laboratory. It is concluded that the microstructural virtual tests can be used as a supplemental tool to evaluate fracture properties of asphalt mixtures.
基金Project(51378006) supported by National Natural Science Foundation of ChinaProject(141076) supported by Huoyingdong Foundation of the Ministry of Education of China+1 种基金Project(2242015R30027) supported by Excellent Young Teacher Program of Southeast University,ChinaProject(BK20140109) supported by the Natural Science Foundation of Jiangsu Province,China
文摘The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.
基金supported by the National Science and Technology Major Project(J2019-IV-0003-0070)the National Natural Science Foundation of China(Grant No.12102320)+1 种基金the Advanced Aviation Power Innovation Workstation Project(HKCX2019-01-003)China Postdoc-toral Science Foundation(2021M692571).
文摘Avirtual wall thicknessmethod is developed to simulate the temperature field of turbine bladeswith thermal barrier coatings(TBCs),to simplify the modeling process and improve the calculation efficiency.The results show that the virtualwall thickness method can improve themesh quality by 20%,reduce the number ofmeshes by 76.7%and save the calculation time by 35.5%,compared with the traditional real wall thickness method.The average calculation error of the two methods is between 0.21%and 0.93%.Furthermore,the temperature at the blade leading edge is the highest and the average temperature of the blade pressure surface is higher than that of the suction surface under a certain service condition.The blade surface temperature presents a high temperature at both ends and a low temperature in themiddle height when the temperature of incoming gas is uniformand constant.The thermal insulation effect of TBCs is the worst near the air film hole,and the best at the blade leading edge.According to the calculated temperature field of the substrate-coating system,the highest thermal insulation temperature of the TC layer is 172.01 K,and the thermal insulation proportions of TC,TGO and BC are 93.55%,1.54%and 4.91%,respectively.
文摘A novel polygonal finite element method (PFEM) based on partition of unity is proposed, termed the virtual node method (VNM). To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.
文摘In this paper, we discuss virtual element method (VEM) approximation of optimal control problem governed by Brinkman equations with control constraints. Based on the polynomial projections and variational discretization of the control variable, we build up the virtual element discrete scheme of the optimal control problem and derive the discrete first order optimality system. A priori error estimates for the state, adjoint state and control variables in L<sup>2</sup> and H<sup>1</sup> norm are derived. The theoretical findings are illustrated by the numerical experiments.
基金This study was supported by“the Fundamental Research Funds for the Central Universities”(Grant No.2015B37814)the Postgraduate Research and Practice Innovation Program of Jiangsu Province(Grant No.KYLX15_0489)+1 种基金the National Natural Science Foundation of China(Grant No.51679081)“the Fundamental Research Funds for the Central Universities”(Grant No.2018B48514).
文摘This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional boundary element method,and can be treated as a potential approach for solving the inverse problems of the heat conduction owing to the boundary-only discretization and semi-analytical algorithm.When the VBEM is applied to the inverse problems,the numerical instability may occur if a virtual boundary is not properly chosen.The method encounters a highly illconditioned matrix for the larger distance between the physical boundary and the virtual boundary,and otherwise is hard to avoid the singularity of the source point.Thus,it must adopt an appropriate regularization method to deal with the ill-posed systems of inverse problems.In this study,the VBEM and different regularization techniques are combined to model the inverse problem of three-dimensional heat conduction in orthotropic media.The proper regularization techniques not only make the virtual boundary to be allocated freer,but also solve the ill-conditioned equation of the inverse problem.Numerical examples demonstrate that the proposed method is efficient,accurate and numerically stable for solving the inverse problems of three-dimensional heat conduction in orthotropic media.
基金supported by the National Natural Science Foundation of China(Grant No.12071100)by the Fundamental Research Funds for the Central Universities(Grant No.2022FRFK060019).
文摘The second-order serendipity virtual element method is studied for the semilinear pseudo-parabolic equations on curved domains in this paper.Nonhomogeneous Dirichlet boundary conditions are taken into account,the existence and uniqueness are investigated for the weak solution of the nonhomogeneous initial-boundary value problem.The Nitschebased projection method is adopted to impose the boundary conditions in a weak way.The interpolation operator is used to deal with the nonlinear term.The Crank-Nicolson scheme is employed to discretize the temporal variable.There are two main features of the proposed scheme:(i)the internal degrees of freedom are avoided no matter what type of mesh is utilized,and(ii)the Jacobian is simple to calculate when Newton’s iteration method is applied to solve the fully discrete scheme.The error estimates are established for the discrete schemes and the theoretical results are illustrated through some numerical examples.
文摘The virtual laminated element method (VLEM) can resolve structural shap e optimization problems with a new method. According to the characteristics of V LEM , only some characterized layer thickness values need be defined as design v ariables instead of boundary node coordinates or some other parameters determini ng the system boundary. One of the important features of this method is that it is not necessary to regenerate the FE(finite element) grid during the optimizati on process so as to avoid optimization failures resulting from some distortion grid elements. Th e thickness distribution in thin plate optimization problems in other studies be fore is of stepped shape. However, in this paper, a continuous thickness distrib ution can be obtained after optimization using VLEM, and is more reasonable. Fur thermore, an approximate reanalysis method named ″behavior model technique″ ca n be used to reduce the amount of structural reanalysis. Some typical examples are offered to prove the effectiveness and practicality of the proposed method.
基金The authors gratefully acknowledges support for this research by the“German Research Foundation”(DFG)in(i)the Collaborative Research Center CRC 1153 and(ii)the Priority Program SPP 2020.
文摘The virtual element method(VEM)can be seen as an extension of the classical finite element method(FEM)based on Galerkin projection.It allows meshes with highly irregular shaped elements,including concave shapes.So far the virtual element method has been applied to various engineering problems such as elasto-plasticity,multiphysics,damage and fracture mechanics.This work focuses on the extension of the virtual element method to efficient modeling of nonlinear elasto-dynamics undergoing large deformations.Within this framework,we employ low-order ansatz functions in two and three dimensions for elements that can have arbitrary polygonal shape.The formulations considered in this contribution are based on minimization of potential function for both the static and the dynamic behavior.Generally the construction of a virtual element is based on a projection part and a stabilization part.While the stiffness matrix needs a suitable stabilization,the mass matrix can be calculated using only the projection part.For the implicit time integration scheme,Newmark-Method is used.To show the performance of the method,various two-and three-dimensional numerical examples in are presented.
文摘应变-旋转(Strain-Rotation,S-R)和分解定理为分析几何非线性问题提供了合理可靠的理论基础,但用有限元求解时会遇到大变形发生后的网格畸变问题。近年提出的虚单元法(Virtual element method,VEM)适用于一般的多边形网格,因此,该文尝试使用一阶虚单元求解基于S-R和分解定理的二维几何非线性问题,以克服网格畸变的影响。基于重新定义的多项式位移空间基函数,推演获得一阶虚单元分析线弹性力学问题时允许位移空间向多项式位移空间的投影表达式;按照虚单元法双线性格式的计算规则,分析处理基于更新拖带坐标法和势能率原理的增量变分方程;进而建立离散系统方程及其矩阵表达形式,并编制MATLAB求解程序;采用常规多边形网格和畸变网格,应用该文算法分析均布荷载下的悬臂梁和均匀内压下的厚壁圆筒变形。结果与已有文献和ANSYS软件的对比表明:该文算法在两种网格中均可有效执行且具备足够数值精度。总体该文算法为基于S-R和分解定理的二维几何非线性问题求解提供了一种鲁棒方法。