杆有限元模型的特征值反问题

1引言非均匀杆的轴向自由振动方程为其中L为杆的长度,A(x)为杆的横截面积,E(x)为弹性模量,ρ(x)是杆的密度．杆的主振动的一般形式为u(x,t)=u(x)sin(wt φ), (2)其中u(x)满足d/dx[E(x)A(x)du/dx] λρ(x)A(x)u=0,λ=ω2．(3)     

用有限元方法求解界面特征值问题

用有限元方法计算椭圆型界面特征值问题,实验数据显示近似特征值的变化规律:界面特征值问题中系数的间断性对协调和非协调Crouzeix-Raviart有限元特征值的收敛性并无影响,而且对协调有限元特征值外推以后得到高精度的解,相应的外推值还提供特征值下界;Crouzeix-Raviart元特征值提供特征值下界,这对一般有界区域如"镂空"型区域也成立.另外,还展示近似特征函数的图形.     

特征值问题有限元逼近位移超收敛的一个定理

本文建立了二阶椭圆特征值问题 m 次有限元逼近位移超收敛定理:‖Pu-u_h‖_(0,∞)≤ch~(m+2),m≥2‖Pu-u_h‖_(0,∞)≤ch~(m+1+(2/ξ),m≥3其中,1<ξ<2,u_h 是 m 次有限元特征函数,u 是某一精确特征函数,P 是 Ritz 投影算子。利用这个定理就可以把二阶线性椭圆方程有限元逼近位移超收敛的结果移植到特征值问题,     

数值积分对特征值有限元外推的影响

§1.引言在分片一致三角形剖分下,用线性有限元法解特征值问题求得近似特征值λ、λ~(h/2).[1]证明了对λ~h.λ~(h/2)作外推可提高收敛阶:     

有限元奇性问题的自适应处理

§1.引言 有限元奇性问题(包括凹角域问题以及方程本身的奇性等问题),由于奇点的存在,致使奇点附近收敛速度明显减慢。处理这一问题的有效方法是构作可靠实用的后验估计量来寻找收敛慢的坏单元,进行局部加密剖分,从而以最快的速度提高整体精度,但工作量最小。 本文的特点在于:1.对剖分条件限制较弱,只要求正规剖分(regular meshes)即     

对称三对角线特征值问题的一个并行修正拟Laguerre算法

在拟Laguerre算法的基础上，提出了用修正拟Lagureer算法来求求解对称三角线特征值问题，并给出了算法的并行实现。     

一维奇异问题有限元的插值校正

有限元的插值校正是一种高精度方法,当解较光滑时,理论上已经作了许多研究,见林群、杨一都的概括论文及陈传淼的研究新框架。对光滑解情形,由于校正量q_4很好地描述了有限元误差e=u_h-u的主部,因此它相当准确地提供了有限元的逐点误差,并可发     

反中心对称矩阵的广义特征值反问题

Given matrix X and diagonal matrix A , the anti-centrosymmetric solutions (A, B) and its optimal approximation of inverse generalized eigenvalue problem AX = BXA have been considered. The general form of such solutions is given and the expression of the optimal approximation solution to a given matrix is derived. The algorithm and one numerical example for solving optimal approximation solution are included.     

加权奇异拟线性椭圆方程特征值问题

在加权Sobolev空间中,利用Galerkin方法及推广的Brouwer定理,研究一类奇异拟线性椭圆方程高阶特征问题,得到了其非平凡弱解的存在性.     

非负矩阵最大特征值的平滑算法

１引　言 本文中Ａ＝（ａｉｊ）表示ｎ阶方阵,Ａ＞０表示Ａ为正矩阵,即ａｉｊ＞０（ｉ,ｊ＝１,２,…,ｎ）;Ａ≥０表示Ａ为非负矩阵,即ａｉｊ≥０（ｉ,ｊ＝１,２,…,ｎ）且至少有一个严格大于号成立,周知,当Ａ＞０时Ａ有一个正特征值λ满足λ＞｜λ｜,其中λ为Ａ的其它任一特征值;当Ａ≥０时Ａ有一个非负特征值λ满足λ≥｜λ｜,其中λ为Ａ的任一特征值．把这样的λ称为Ａ的最大特征值,为强调它属于Ａ,记作λ（Ａ）．同时,把与λ（Ａ）对应的Ａ的特征向量记作ｘ（Ａ）． 对Ａ≥０,记当Ｒｔ＞０（ｉ＝１,２,…,ｎ）时…     

THE ADAPTIVE IMMERSED INTERFACE FINITE ELEMENT METHOD FOR ELASTICITY INTERFACE PROBLEMS

In this paper, we propose adaptive finite element methods with error control for solving elasticity problems with discontinuous coefficients. The meshes in the methods do not need to fit the interfaces. We establish a residual-based a posteriori error estimate which is $λ$-independent multiplicative constants; the Lamé constant $λ$ steers the incompressibility. The error estimators are then implemented and tested with promising numerical results which will show the competitive behavior of the adaptive algorithm.     

A POSTERIORI ERROR ESTIMATES IN ADINI FINITE ELEMENT FOR EIGENVALUE PROBLEMS

In this paper, we discuss a posteriori error estimates of the eigenvalue $\lambda_h$ given by Adini nonconforming finite element. We give an asymptotically exact error estimator of the $\lambda_h$. We prove that the order of convergence of the $\lambda_h$ is just 2 and the $\lambda_h$ converge from below for sufficiently small $h$.     

振动杆有限差分模型的一类特征值反问题

1引言杆是重要的工程构件之一,具有分布质量的杆的纵向振动由下面的偏微分方程描述:     

二阶椭圆问题的一类广义有限元法

本文提出了求解二阶椭圆问题的一类广义有限元方法,分析了广义有限元方法的优越性,证明了二阶椭圆问题的广义有限元方法具有比标准的Galerkin有限元方法更高阶的收敛速度,根据插值算子的性质,进一步证明了有限元解的亏量迭代校正收敛到广义有限元解,并用数值例子说明广义有限元方法是有效的.     

非自共轭椭圆特征值问题有限元插值校正

本文研究非自共轭椭圆特征值问题有限元插值校正方案.基于插值校正和广义Rayleigh商加速技巧,用三角形线性元二次插值、双二次元双四次插值得到了较好的结果,并用三线性元的三二次捕值将捅值校正推广到三维.     

AN ITERATIVE HYBRIDIZED MIXED FINITE ELEMENT METHOD FOR ELLIPTIC INTERFACE PROBLEMS WITH STRONGLY DISCONTINUOUS COEFFICIENTS

An iterative algorithm is proposed and analyzed based on a hybridized mized finite element method for numerically solving two-phase generalized Stefan interface problems with strongly discontinuous solutions,conormal derivatives,and coefficients.This algorithm iteratively solves small problems for each single phase with good accuracy and exchange information at the interface to advance the iteration until convergence ,following the idea of Schwarz Alternating Methods,Error estimates are derived to show that this algorithm always converges provided that relaxation parameters are suitably chosen,Numeric exper-iments with matching and non-matching grids at the interface from different phases are performed to show the accuracy of the method for capturing discontinuities in the solutions and coefficients.In contrast to standard numerical methods,the accuracy of our method does not seem to deteriorate as the coefficient discontinuity increases.     

求解不适定问题的快速Landweber迭代法

本文从一般迭代法的级数形式出发，将一般迭代法的每一步分解为矩阵计算和求解两步，并对其中的矩阵计算部分进行了修改，在此基础上提出了快速迭代法，最后通过数值实验验证了我们的算法不仅提高了计算速度，同时也大大减少了计算量，是一种效率很高的算法。     

一类抛物型方程系数反问题的分裂算法

考虑利用终端时刻的温度u(x,T)=Z_T(x)反演热传导方程u_t-a~2u_(xx) q(x)u=0,x∈(0,1)中的未知系数q(x)的反问题.通过引进变换v(x,t)=(u_t(x,t)/u(x,t))将此非线性不适定问题的求解分解为两步.首先利用输入数据迭代求解一个非线性的正问题(该过程独立于未知系数),得到其迭代解v~(k)(x,t).其次利用q(x)与v(x,t)的关系式求出q(x)的近似解.对提出的反演方法,证明了采用的变换的可行性,得到了原反问题与由变换后的非线性正问题反演q(x)的等价性并且证明了迭代解的收敛性,给出了收敛速度.数值结果表明了该方法的有效性.     

A MODIFIED INVERSE ITERATION FOR A LARGE SPARSE SPD GENERALIZED EIGENPROBLEM

In this paper, an algorithm based on a shifted inverse power iteration for computing generalized eigenvalues with corresponding eigenvectors of a large scale sparse symmetric positive definite matrix pencil is presented. It converges globally with a cubic asymptotic convergence rate, preserves sparsity of the original matrices and is fully parallelizable. The algebraic multilevel itera-tion method (AMLI) is used to improve the efficiency when symmetric positive definite linear equa-tions need to be solved.