首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 203 毫秒
1.
张然 《计算数学》2020,42(1):1-17
本文考虑弱有限元(简称WG)方法在线弹性问题中的应用.WG方法是传统有限元方法的推广,用于偏微分方程的数值求解.和传统有限元一样,它的基本思想源于变分原理.WG方法的特点是使用在剖分单元内部和剖分单元边界上分别有定义的分片多项式函数(即弱函数)作为近似函数来逼近真解,并针对弱函数定义相应的弱微分算子代入数值格式进行计算.除此之外,WG方法允许在数值格式中引进稳定子以实现近似函数的弱连续性.WG方法具有允许使用任意多边形或多面体剖分,数值格式与逼近函数构造简单,易于满足相应的稳定性条件等优点.本文考虑WG方法在求解线弹性问题中的应用.围绕线弹性问题数值求解中常见的三个问题,即:数值格式的强制性,闭锁性,应力张量的对称性介绍WG方法在线弹性问题求解中的应用.  相似文献   

2.
基于平方和松弛和有理向量恢复,提出了一种符号数值混合计算方法来构造多项式Lyapunov函数以判定非线性混成系统的稳定性,首先,为Lyapunov函数预定一个给定次数的多项式模板,则Lyapunov函数构造问题可转化为相应的带参数的多项式优化问题,然后运用平方和松弛方法求得一个近似的数值多项式Lyapunov函数,再应用高斯-牛顿精化和有理向量恢复将数值多项式转化为验证的有理多项式Lyapunov函数.  相似文献   

3.
陈明卿  谢小平 《计算数学》2021,43(3):279-300
本文针对带有随机杨氏模量和荷载的平面线弹性问题,提出了一类随机弱Galerkin有限元方法.先利用Karhunen-Loève展开把随机项参数化,将方程转化为一个确定性问题;再采用弱Galerkin有限元法和$k$-/$p$-型方法分别离散空间区域和随机场.在弱Galerkin离散中,用分片$s(s\geqslant 1$)和$s+1$次多项式逼近单元内部的应力和位移,用分片$s$次多项式逼近位移在单元边界上的迹.证明了该方法关于空间网格尺度最优且与Lamé常数$\lambda$一致无关的误差估计.最后通过数值算例验证了理论结果.  相似文献   

4.
该文在弱双代数$H$上给出了扭曲积$(H^\sigma,\cdot_\sigma)$成为弱双代数的充分必要条件.设$[B, H, \tau]$是一个弱斜配对, 并且$\tau$可逆,则在某个条件下弱双交叉积$B\bowtie_\tau H$是一个弱双代数. 如果$(B,H, \sigma)$是弱相关Long双代数, 并且$\sigma$可逆,则弱双交叉积$B^{OP}\bowtie_\sigma H$可以被构造. 它的乘法是:$(x\otimes h)(y\otimes g)=\Sigma\sigma(y_1, h_1)y_2x\otimes h_2g\sigma^{-1}(y_3, h_3),$ 特别地, 如果$(B, H,\sigma)$是相关Long双代数, 则$(B^{OP \bowtie_\sigma H,\beta)$是Long双代数当且仅当对任意$b, d\in B^{OP}; g, \ell\in H$,$\Sigma\sigma^{-1}(b, g_2\ell)\sigma(d, g_1)=\Sigma\sigma^{-1}(b,\ell g_1)\sigma(d, g_2),$ 其中$B$为$H$的子Hopf代数,$\beta$定义为$\beta(b\bowtie_\sigma h\otimes c\bowtie_\sigma g)=\varepsilon_H(h)\varepsilon_{B^{OP}}(c)\sigma^{-1}(b, g).$ 对于Sweedler 4维Hopf代数$H$, 作者给出一个例子说明:此弱双交叉积$(B^{OP}\bowtie_\sigma H, \beta)$不仅是一个Long双代数,而且是一个非可换和非余可换的8维Hopf代数. 最后, 设$B,H$都是弱双代数, $\sigma: B\otimes H\rightarrow k$是一个线性映射, 作者给出了$(B,\sigma,\leftharpoonup, \Delta_B)$是弱相关右$(H, B)$ -重模代数的充分必要条件.  相似文献   

5.
作者引入了非齐型空间上的弱Herz空间,并建立了一类次线性算子在这些空间中的弱型估计. 作为应用, 证明了由Calder\'on-Zygmund算子和$\os$函数生成的交换子在弱Herz空间中的弱型估计,其中$r\ge1$. 并且Orlicz空间$\os$当$r=1$时即为$\rb$空间;当$r>1$时为$\rb$的子空间.  相似文献   

6.
利用耗散不等式研究了切换多项式非线性系统的输入-状态稳定性分析问题,在任意切换信号下,给出了使得切换多项式非线性系统输入-状态稳定的充分条件.采用平方和分解方法来寻找切换多项式非线性系统的输入-状态稳定共同Lyapunov函数.数值算例验证了所提方法的可行性.  相似文献   

7.
王军平  叶秀  张然 《计算数学》2016,38(3):289-308
本文简述弱有限元方法(weak Galerkin finite element met,hods)的数学基本原理和计算机实现.弱有限元方法对间断函数引入广义弱微分,并将其应用于偏微分方程相应的变分形式进行数值求解,而数值解的弱连续性则通过稳定子或光滑子来实现.弱有限元方法针对广义函数而构建,是经典有限元方法的一种自然拓广,且能够弥补经典有限元方法的某些缺憾,也因此在科学与工程计算领域具有广泛的应用前景.  相似文献   

8.
设 $k, m$ 是两个正整数, $a\ ( \ne 0)$是有穷复数. $\mathcal{F}$ 是区域 $D$ 内的一族亚纯函数, $f\in\mathcal{F}$ 的零点重数至少为 $k$, $P$ 是多项式,次数或者 ${\rm deg}\, P\geq3$ 或者 ${\rm deg}\, P=2$ 且 $P$ 只有一个不同的零点.若对于 $\mathcal{F}$ 中的任意两个函数 $f$ 和 $g$, $P(f){({f^{(k)}})^m}$ 与 $P(g){({g^{(k)}})^m}$ 在 $D$ 内 IM 分担 $a$, 则 $\mathcal{F}$ 在 $D$ 内正规.  相似文献   

9.
给出了Banach空间的p-弱近似性质和p-有界弱近似性质的定义,获得了这些性质的一些刻画.利用这些刻画证明了如果一个Banach空间X的对偶空间X~*有p-弱近似性质(或p-有界弱近似性质),则X有p-弱近似性质(或p-有界弱近似性质),在一般情况下反之不成立.  相似文献   

10.
利用稳定化方法讨论拉格朗日乘子法得到的具有弱对称应力的线弹性问题. 用线性元和分片常数分别逼近变分问题的应力和位移. 并通过添加稳定项$G_1(\cdot,\cdot)$, $G_2(\cdot,\cdot)$和$G_3(\cdot,\cdot)$ 使相应混合离散变分问题满足弱BB条件. 接着详细研究了变分问题的解与稳定混合有限元解之间的误差估计,最后用两个数值算例验证理论分析的有效性.  相似文献   

11.
In this paper, a weak Galerkin finite element method for the Oseen equations of incompressible fluid flow is proposed and investigated. This method is based on weak gradient and divergence operators which are designed for the finite element discontinuous functions. Moreover, by choosing the usual polynomials of degree i ≥ 1 for the velocity and polynomials of degree i ? 1 for the pressure and enhancing the polynomials of degree i ? 1 on the interface of a finite element partition for the velocity, this new method has a lot of attractive computational features: more general finite element partitions of arbitrary polygons or polyhedra with certain shape regularity, fewer degrees of freedom and parameter free. Stability and error estimates of optimal order are obtained by defining a weak convection term. Finally, a series of numerical experiments are given to show that this method has good stability and accuracy for the Oseen problem.  相似文献   

12.
The second order elliptic equation, which is also know as the diffusion-convection equation, is of great interest in many branches of physics and industry. In this paper, we use the weak Galerkin finite element method to study the general second order elliptic equation. A weak Galerkin finite element method is proposed and analyzed. This scheme features piecewise polynomials of degree $k\geq 1$ on each element and piecewise polynomials of degree $k-1\geq 0$ on each edge or face of the element. Error estimates of optimal order of convergence rate are established in both discrete $H^1$ and standard $L^2$ norm. The paper also presents some numerical experiments to verify the efficiency of the method.  相似文献   

13.
This paper introduces a weak Galerkin (WG) finite element method for the Stokes equations in the primal velocity-pressure formulation. This WG method is equipped with stable finite elements consisting of usual polynomials of degree k≥1 for the velocity and polynomials of degree k?1 for the pressure, both are discontinuous. The velocity element is enhanced by polynomials of degree k?1 on the interface of the finite element partition. All the finite element functions are discontinuous for which the usual gradient and divergence operators are implemented as distributions in properly-defined spaces. Optimal-order error estimates are established for the corresponding numerical approximation in various norms. It must be emphasized that the WG finite element method is designed on finite element partitions consisting of arbitrary shape of polygons or polyhedra which are shape regular.  相似文献   

14.
In this article, a new weak Galerkin mixed finite element method is introduced and analyzed for the Helmholtz equation with large wave numbers. The stability and well‐posedness of the method are established for any wave number k without mesh size constraint. Allowing the use of discontinuous approximating functions makes weak Galerkin mixed method highly flexible in term of little restrictions on approximations and meshes. In the weak Galerkin mixed finite element formulation, approximation functions can be piecewise polynomials with different degrees on different elements and meshes can consist elements with different shapes. Suboptimal order error estimates in both discrete H1 and L2 norms are established for the weak Galerkin mixed finite element solutions. Numerical examples are tested to support the theory.  相似文献   

15.
Weak Galerkin finite element method is introduced for solving wave equation with interface on weak Galerkin finite element space $(\mathcal{P}_k(K), \mathcal{P}_{k−1}(∂K), [\mathcal{P}_{k−1}(K)]^2).$ Optimal order a priori error estimates for both space-discrete scheme and implicit fully discrete scheme are derived in $L^∞(L^2)$ norm. This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes. Finite element algorithm presented here can contribute to a variety of hyperbolic problems where physical domain consists of heterogeneous media.  相似文献   

16.
In this paper, we introduce numerical schemes and their analysis based on weak Galerkin finite element framework for solving 2‐D reaction–diffusion systems. Weak Galerkin finite element method (WGFEM) for partial differential equations relies on the concept of weak functions and weak gradients, in which differential operators are approximated by weak forms through the Green's theorem. This method allows the use of totally discontinuous functions in the approximation space. In the current work, the WGFEM solves reaction–diffusion systems to find unknown concentrations (u, v) in element interiors and boundaries in the weak Galerkin finite element space WG(P0, P0, RT0) . The WGFEM is used to approximate the spatial variables and the time discretization is made by the backward Euler method. For reaction–diffusion systems, stability analysis and error bounds for semi‐discrete and fully discrete schemes are proved. Accuracy and efficiency of the proposed method successfully tested on several numerical examples and obtained results satisfy the well‐known result that for small values of diffusion coefficient, the steady state solution converges to equilibrium point. Acquired numerical results asserted the efficiency of the proposed scheme.  相似文献   

17.
This paper considers weak Galerkin finite element approximations on polygonal/polyhedral meshes for a quasistatic Maxwell viscoelastic model. The spatial discretization uses piecewise polynomials of degree $k (k ≥ 1)$ for the stress approximation, degree $k+1$ for the velocity approximation, and degree $k$ for the numerical trace of velocity on the inter-element boundaries. The temporal discretization in the fully discrete method adopts a backward Euler difference scheme. We show the existence and uniqueness of the semi-discrete and fully discrete solutions, and derive optimal a priori error estimates. Numerical examples are provided to support the theoretical analysis.  相似文献   

18.
This paper presents a new numerical method for the compressible Navier–Stokes equations governing the flow of an ideal isentropic gas. To approximate the continuity equation, the method utilizes a discontinuous Galerkin discretization on piecewise constants and a basic upwind flux. For the momentum equation, the method is a new combined discontinuous Galerkin and finite element method approximating the velocity in the Crouzeix–Raviart finite element space. While the diffusion operator is discretized in a standard fashion, the convection and time-derivative are discretized using discontinuous Galerkin on the element average velocity and a Lax–Friedrich type flux. Our main result is convergence of the method to a global weak solution as discretization parameters go to zero. The convergence analysis constitutes a numerical version of the existence analysis of Lions and Feireisl.  相似文献   

19.
This article concerns numerical approximation of a parabolic interface problem with general $L^2$ initial value. The problem is discretized by a finite element method with a quasi-uniform triangulation of the domain fitting the interface, with piecewise linear approximation to the interface. The semi-discrete finite element problem is furthermore discretized in time by the $k$-step backward difference formula with $ k=1,\ldots,6 $. To maintain high-order convergence in time for possibly nonsmooth $L^2$ initial value, we modify the standard backward difference formula at the first $k-1$ time levels by using a method recently developed for fractional evolution equations. An error bound of $\mathcal{O}(t_n^{-k}\tau^k+t_n^{-1}h^2|\log h|)$ is established for the fully discrete finite element method for general $L^2$ initial data.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号