The element energy projection (EEP) method for computation of super- convergent resulting in a one-dimensional finite element method (FEM) is successfully used to self-adaptive FEM analysis of various linear probl...The element energy projection (EEP) method for computation of super- convergent resulting in a one-dimensional finite element method (FEM) is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton's method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation (ODE) of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.展开更多
Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted ...Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.展开更多
Based on the newly-developed element energy projection (EEP) method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite ele...Based on the newly-developed element energy projection (EEP) method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method (FEM) is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.展开更多
The newly proposed element energy projection(EEP) method has been applied to the computation of super_convergent nodal stresses of Timoshenko beam elements.General formulas based on element projection theorem were der...The newly proposed element energy projection(EEP) method has been applied to the computation of super_convergent nodal stresses of Timoshenko beam elements.General formulas based on element projection theorem were derived and illustrative numerical examples using two typical elements were given.Both the analysis and examples show that EEP method also works very well for the problems with vector function solutions.The EEP method gives super_convergent nodal stresses,which are well comparable to the nodal displacements in terms of both convergence rate and error magnitude.And in addition,it can overcome the “shear locking” difficulty for stresses even when the displacements are badly affected.This research paves the way for application of the EEP method to general one_dimensional systems of ordinary differential equations.展开更多
For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan's element energy projection (EEP) method has ...For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan's element energy projection (EEP) method has the accuracy O(h^min{2k,k+4}) The theoretical analysis coincides the reported numerical results.展开更多
对于结构动力分析中的离散系统运动方程,现有算法的计算精度和效率均依赖于时间步长的选取,这是时间域问题求解的难点.基于EEP(element energy projection)超收敛计算的自适应有限元法,以EEP超收敛解代替未知真解,估计常规有限元解的误...对于结构动力分析中的离散系统运动方程,现有算法的计算精度和效率均依赖于时间步长的选取,这是时间域问题求解的难点.基于EEP(element energy projection)超收敛计算的自适应有限元法,以EEP超收敛解代替未知真解,估计常规有限元解的误差,并自动细分网格,目前已对诸类以空间坐标为自变量的边值问题取得成功.对离散系统运动方程建立弱型Galerkin有限元解,引入基于EEP法的自适应求解策略,在时间域上自动划分网格,最终得到所求时域内任一时刻均满足给定误差限的动位移解,进而建立了一种时间域上的新型自适应求解算法.展开更多
将新近提出的C0有限元后处理中超收敛解答计算的单元能量投影(Element Energy Projection,简称EEP)法推广到一维C1类有限元。根据单元投影定理具体推导了一般梁单元的计算公式,并对两个有代表性的单元给出了数值算例。分析和算例表明,EE...将新近提出的C0有限元后处理中超收敛解答计算的单元能量投影(Element Energy Projection,简称EEP)法推广到一维C1类有限元。根据单元投影定理具体推导了一般梁单元的计算公式,并对两个有代表性的单元给出了数值算例。分析和算例表明,EEP法在一维C1类有限元中再次获得令人满意的效果,即对任一单元中的任一点,从位移一直到三阶导数(如梁的挠度、转角、弯矩、剪力),匀可获得与结点位移精度相当的超收敛结果,而且可精确满足自然边界条件。展开更多
利用单元能量投影(Element Energy Projection,简称EEP)法所计算的EEP超收敛解,在不改变有限元网格及其整体刚度矩阵的情况下,导出残差的等效结点荷载向量,只经回代过程即可得到具有更高阶精度的结点位移的误差估计,使结点位移精度得到...利用单元能量投影(Element Energy Projection,简称EEP)法所计算的EEP超收敛解,在不改变有限元网格及其整体刚度矩阵的情况下,导出残差的等效结点荷载向量,只经回代过程即可得到具有更高阶精度的结点位移的误差估计,使结点位移精度得到极大提高。该文以一般的二阶常微分方程边值和初值问题为例,给出算法和相应的数值算例。从中可以看出,本法十分简单而高效:对于m≥1次单元,采用EEP简约格式和凝聚格式修正后的结点位移,分别具有O(h^(2m+2))和O(h^(3m+mod(m,2)))的超常规的超收敛阶。该文给出了典型算例,并对该法的进一步拓展和应用作了讨论。展开更多
基于单元能量投影(element energy projection,EEP)法自适应分析在杆件静力问题以及离散系统运动方程组中所取得的成果,以直杆轴向受迫振动为例,研究并建立了一种在时间域和一维空间域同时实现自适应分析的方法.该方法在时间和空间两个...基于单元能量投影(element energy projection,EEP)法自适应分析在杆件静力问题以及离散系统运动方程组中所取得的成果,以直杆轴向受迫振动为例,研究并建立了一种在时间域和一维空间域同时实现自适应分析的方法.该方法在时间和空间两个维度都采用连续的Galerkin有限元法(finite element method,FEM)进行求解,根据半离散的思想,由空间有限元离散将模型问题的偏微分控制方程转化为离散系统运动方程组,对该方程组进行时域有限元自适应求解;然后再基于空间域超收敛计算的EEP解对空间域进行自适应,直至最终的时空网格下动位移解答的精度逐点均满足给定误差限要求.文中对其基本思想、关键技术和实施策略进行了阐述,并给出了包括地震波输入下的典型算例以展示该法有效可靠.展开更多
基金supported by the National Natural Science Foundation of China(Nos.51378293,51078199,50678093,and 50278046)the Program for Changjiang Scholars and the Innovative Research Team in University of China(No.IRT00736)
文摘The element energy projection (EEP) method for computation of super- convergent resulting in a one-dimensional finite element method (FEM) is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton's method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation (ODE) of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.
基金Project supported by the National Natural Science Foundation of China (No.50278046)
文摘Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.
基金the National Natural Science Foundation of China(No.50678093)Program for Changjiang Scholars and Innovative Research Team in University(No.IRT00736)
文摘Based on the newly-developed element energy projection (EEP) method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method (FEM) is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.
文摘The newly proposed element energy projection(EEP) method has been applied to the computation of super_convergent nodal stresses of Timoshenko beam elements.General formulas based on element projection theorem were derived and illustrative numerical examples using two typical elements were given.Both the analysis and examples show that EEP method also works very well for the problems with vector function solutions.The EEP method gives super_convergent nodal stresses,which are well comparable to the nodal displacements in terms of both convergence rate and error magnitude.And in addition,it can overcome the “shear locking” difficulty for stresses even when the displacements are badly affected.This research paves the way for application of the EEP method to general one_dimensional systems of ordinary differential equations.
基金Project supported by the National Natural Science Foundation of China (Nos. 10571046, 10371038)
文摘For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan's element energy projection (EEP) method has the accuracy O(h^min{2k,k+4}) The theoretical analysis coincides the reported numerical results.
文摘对于结构动力分析中的离散系统运动方程,现有算法的计算精度和效率均依赖于时间步长的选取,这是时间域问题求解的难点.基于EEP(element energy projection)超收敛计算的自适应有限元法,以EEP超收敛解代替未知真解,估计常规有限元解的误差,并自动细分网格,目前已对诸类以空间坐标为自变量的边值问题取得成功.对离散系统运动方程建立弱型Galerkin有限元解,引入基于EEP法的自适应求解策略,在时间域上自动划分网格,最终得到所求时域内任一时刻均满足给定误差限的动位移解,进而建立了一种时间域上的新型自适应求解算法.
文摘将新近提出的C0有限元后处理中超收敛解答计算的单元能量投影(Element Energy Projection,简称EEP)法推广到一维C1类有限元。根据单元投影定理具体推导了一般梁单元的计算公式,并对两个有代表性的单元给出了数值算例。分析和算例表明,EEP法在一维C1类有限元中再次获得令人满意的效果,即对任一单元中的任一点,从位移一直到三阶导数(如梁的挠度、转角、弯矩、剪力),匀可获得与结点位移精度相当的超收敛结果,而且可精确满足自然边界条件。
文摘利用单元能量投影(Element Energy Projection,简称EEP)法所计算的EEP超收敛解,在不改变有限元网格及其整体刚度矩阵的情况下,导出残差的等效结点荷载向量,只经回代过程即可得到具有更高阶精度的结点位移的误差估计,使结点位移精度得到极大提高。该文以一般的二阶常微分方程边值和初值问题为例,给出算法和相应的数值算例。从中可以看出,本法十分简单而高效:对于m≥1次单元,采用EEP简约格式和凝聚格式修正后的结点位移,分别具有O(h^(2m+2))和O(h^(3m+mod(m,2)))的超常规的超收敛阶。该文给出了典型算例,并对该法的进一步拓展和应用作了讨论。
文摘基于单元能量投影(element energy projection,EEP)法自适应分析在杆件静力问题以及离散系统运动方程组中所取得的成果,以直杆轴向受迫振动为例,研究并建立了一种在时间域和一维空间域同时实现自适应分析的方法.该方法在时间和空间两个维度都采用连续的Galerkin有限元法(finite element method,FEM)进行求解,根据半离散的思想,由空间有限元离散将模型问题的偏微分控制方程转化为离散系统运动方程组,对该方程组进行时域有限元自适应求解;然后再基于空间域超收敛计算的EEP解对空间域进行自适应,直至最终的时空网格下动位移解答的精度逐点均满足给定误差限要求.文中对其基本思想、关键技术和实施策略进行了阐述,并给出了包括地震波输入下的典型算例以展示该法有效可靠.