本文作者:钢结构设计

有限元分析的节点和单元(有限元分析的节点和单元有哪些)

钢结构设计 2个月前 ( 06-23 ) 4154 抢沙发
有限元法是一种有效解决数学问题的解题方法。其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,单元上所作用的力等效到节点上,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,就是用叉值函数来近似代替 ,借助于变分原理或加权余量法,将微分方程离散求解。
本文目录

什么是有限元

有限元法是一种有效解决数学问题的解题方法。其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,单元上所作用的力等效到节点上,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,就是用叉值函数来近似代替 ,借助于变分原理或加权余量法,将微分方程离散求解。

有限元中壳单元和实体单元有什么本质区别

从本质上来说,所有单元都是基于变分原理的,通过插值函数将节点位移与单元内部应变联系都一起,再通过虚功原理或变分原理得到单元刚度与恢复力的计算公式。壳单元与实体元的节点数、节点自由度数、单元位移模式(插值函数)不一样,所以导致他们的单刚和恢复力表现形式不一样,而这就是有限元实体单元和壳单元二者的本质区别1,单元自由度不同;2,单元数量不同;3,基本理论不同,4,应用对象也不同。1、单元自由度。有限元中实体单元只有三个自由度,而壳单元有六个自由度。

有限元方法与有限差分到底有什么区别

1.1 概念
有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
1.2 差分格式
(1)从格式的精度来划分,有一阶格式、二阶格式和高阶格式。
(2)从差分的空间形式来考虑,可分为中心格式和逆风格式。
(3)考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
1.3 构造差分的方法
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
2. FEM
2.1 概述
有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。
2.2 原理
有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。
根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。
(1)从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法;
(2)从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格;
(3)从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。
不同的组合同样构成不同的有限元计算格式。
对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的*方误差最小;在配置法中,先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。
有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等。常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比。在二维有限元中,三角形单元应用的最早,近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元,常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。
2.3 基本原理与解题步骤
对于有限元方法,其基本思路和解题步骤可归纳为:
(1)建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点。
(2)区域单元剖分,根据求解区域的形状及实际问题的物理特点,将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点序号和相应的边界值。
(3)确定单元基函数,根据单元中节点数目及对近似解精度的要求,选择满足一定插值条件的插值函数作为单元基函数。有限元方法中的基函数是在单元中选取的,由于各单元具有规则的几何形状,在选取基函数时可遵循一定的法则。
(4)单元分析:将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近;再将近似函数代入积分方程,并对单元区域进行积分,可获得含有待定系数(即单元中各节点的参数值)的代数方程组,称为单元有限元方程。
(5)总体合成:在得出单元有限元方程之后,将区域中所有单元有限元方程按一定法则进行累加,形成总体有限元方程。
(6)边界条件的处理:一般边界条件有三种形式,分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。对于自然边界条件,一般在积分表达式中可自动得到满足。对于本质边界条件和混合边界条件,需按一定法则对总体有限元方程进行修正满足。
(7)解有限元方程:根据边界条件修正的总体有限元方程组,是含所有待定未知量的封闭方程组,采用适当的数值计算方法求解,可求得各节点的函数值。
3. 有限体积法
有限体积法(FiniteVolumeMethod)又称为控制体积法。其基本思路是:将计算区域划分为一系列不重复的控制体积,并使每个网格点周围有一个控制体积;将待解的微分方程对每一个控制体积积分,便得出一组离散方程。其中的未知数是网格点上的因变量的数值。为了求出控制体积的积分,必须假定值在网格点之间的变化规律,即假设值的分段的分布的分布剖面。从积分区域的选取方法看来,有限体积法属于加权剩余法中的子区域法;从未知解的近似方法看来,有限体积法属于采用局部近似的离散方法。简言之,子区域法属于有限体积发的基本方法。有限体积法的基本思路易于理解,并能得出直接的物理解释。离散方程的物理意义,就是因变量在有限大小的控制体积中的守恒原理,如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。限体积法得出的离散方程,要求因变量的积分守恒对任意一组控制体积都得到满足,对整个计算区域,自然也得到满足。这是有限体积法吸引人的优点。有一些离散方法,例如有限差分法,仅当网格极其细密时,离散方程才满足积分守恒;而有限体积法即使在粗网格情况下,也显示出准确的积分守恒。就离散方法而言,有限体积法可视作有限单元法和有限差分法的中间物。有限单元法必须假定值在网格点之间的变化规律(既插值函数),并将其作为近似解。有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。有限体积法只寻求的结点值,这与有限差分法相类似;但有限体积法在寻求控制体积的积分时,必须假定值在网格点之间的分布,这又与有限单元法相类似。在有限体积法中,插值函数只用于计算控制体积的积分,得出离散方程之后,便可忘掉插值函数;如果需要的话,可以对微分方程中不同的项采取不同的插值函数。
4. 比较分析
有限差分法(FDM):直观,理论成熟,精度可眩但是不规则区域处理繁琐,虽然网格生成可以使FDM应用于不规则区域,但是对区域的连续性等要求较严。使用FDM的好处在于易于编程,易于并行。
有限元方法(FEM):适合处理复杂区域,精度可眩缺憾在于内存和计算量巨大。并行不如FDM和FVM直观。不过FEM的并行是当前和将来应用的一个不错的方向。
有限容积法:适于流体计算,可以应用于不规则网格,适于并行。但是精度基本上只能是二阶了。FVM的优势正逐渐显现出来,FVM在应力应变,高频电磁场方面的特殊的优点正在被人重视。
比较一下:
有限容积法和有限差分法:一个区别就是有限容积法的截差是不定的(跟取的相邻点有关,积分方法离散方程),而有限差分就可以直接知道截差(微分方法离散方程)。有限容积法和有限差分法最本质的区别是,前者是根据积分方程推导出来的(即对每个控制体积分),后者直接根据微分方程推导出来,所以前者的精度不但取决于积分时的精度,还取决与对导数处理的精度,一般有限容积法总体的精度为二阶,因为积分的精度限制,当然有限容积法对于守恒型方程导出的离散方程可以保持守恒型;而后者直接由微分方程导出,不涉及积分过程,各种导数的微分借助Taylor展开,直接写出离散方程,当然不一定有守恒性,精度也和有限容积法不一样,一般有限差分法可以使精度更高一些。
当然二者有联系,有时导出的形式一样,但是概念上是不一样的。
至于有限容积法和有限元相比,有限元在复杂区域的适应性对有限容积是毫无优势可言的,至于有限容积的守恒性,物理概念明显的这些特点,有限元是没有的。目前有限容积在精度方面与有限元法有些差距。
有限元方法比有限差分优越的方面主要在能适应不规则区域,但是这只是指的是传统意义上的有限差分,现在发展的一些有限差分已经能适应不规则区域。对于椭圆型方程,如果区域规则,传统有限差分和有限元都能解,在求解效率,这里主要指编程负责度和收敛快慢、内存需要,肯定有限差分有优势。

有限元离散基本原理

对于不同物理性质和数学模型的问题,有限元求解法的基本步骤是相同的,只是具体公式推导和运算求解不同。有限元求解问题的基本步骤通常为:第一步:问题及求解域定义:根据实际问题近似确定求解域的物理性质和几何区域。第二步:求解域离散化:将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域,习惯上称为有限元网络划分。显然单元越小(网络越细)则离散域的近似程度越好,计算结果也越精确,但计算量及误差都将增大,因此求解域的离散化是有限元法的核心技术之一。第三步:确定状态变量及控制方法:一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示,为适合有限元求解,通常将微分方程化为等价的泛函形式。第四步:单元推导:对单元构造一个适合的近似解,即推导有限单元的列式,其中包括选择合理的单元坐标系,建立单元试函数,以某种方法给出单元各状态变量的离散关系,从而形成单元矩阵(结构力学中称刚度阵或柔度阵)。为保证问题求解的收敛性,单元推导有许多原则要遵循。 对工程应用而言,重要的是应注意每一种单元的解题性能与约束。例如,单元形状应以规则为好,畸形时不仅精度低,而且有缺秩的危险,将导致无法求解。第五步:总装求解:将单元总装形成离散域的总矩阵方程(联合方程组),反映对近似求解域的离散域的要求,即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行,状态变量及其导数(可能的话)连续性建立在结点处。第六步:联立方程组求解和结果解释:有限元法最终导致联立方程组。联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量,将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。

矩阵法和有限元法的区别

矩阵位移法只能用于分析具有已知单元节点力——单元节点位移关系的杆系结构,而不能分析非杆系的连续体结构。两者的区别主要表现为建立的基本原理不同以及求解的侧重点不同。矩阵位移法侧重于根据结构的几何性质和弹性性质建立整体刚度矩阵,根据结构的受载情况形成整体荷载向量,求出各单元的杆端内力,需要叠加固端力;而有限元法侧重于根据已知外荷载的作用和结构材料的性质将结构进行离散化后求解,结点力或是结点应力是直接对位移求导获得的,丢失了固端力一项,致使应力的精度有所下降。通过对*面刚架结构进行求解,得出结构内力,而这一区别在算例结果上得到了验证。

什么是有限元法和有限差分法

有限元法,有限差分法和有限体积法的区别\x0d有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用.该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域.有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组.该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法. 对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式.从差分的空间形式来考虑,可分为中心格式和逆风格式.考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等.目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式.差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定.\x0d构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法.其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度.通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式.\x0d有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解.采用不同的权函数和插值函数形式,便构成不同的有限元方法.有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟.在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成.在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等.根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式.从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等.不同的组合同样构成不同的有限元计算格式.对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的*方误差最小;在配置法中,先在计算域内选取N个配置点.令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0.插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数.有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值.单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等.常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比.在二维有限元中,三角形单元应用的最早,近来四边形等参元的应用也越来越广.对于二维三角形和四边形电源单元,常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等.\x0d对于有限元方法,其基本思路和解题步骤可归纳为\x0d(1)建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点.\x0d(2)区域单元剖分,根据求解区域的形状及实际问题的物理特点,将区域剖分为若干相互连接、不重叠的单元.区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点序号和相应的边界值.\x0d(3)确定单元基函数,根据单元中节点数目及对近似解精度的要求,选择满足一定插值条件的插值函数作为单元基函数.有限元方法中的基函数是在单元中选取的,由于各单元具有规则的几何形状,在选取基函数时可遵循一定的法则.\x0d(4)单元分析:将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近;再将近似函数代入积分方程,并对单元区域进行积分,可获得含有待定系数(即单元中各节点的参数值)的代数方程组,称为单元有限元方程.\x0d(5)总体合成:在得出单元有限元方程之后,将区域中所有单元有限元方程按一定法则进行累加,形成总体有限元方程.\x0d(6)边界条件的处理:一般边界条件有三种形式,分为本质边界条件(狄里克雷边界条件 )、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件).对于自然边界条件,一般在积分表达式中可自动得到满足.对于本质边界条件和混合边界条件,需按一定法则对总体有限元方程进行修正满足.\x0d(7)解有限元方程:根据边界条件修正的总体有限元方程组,是含所有待定未知量的封闭方程组,采用适当的数值计算方法求解,可求得各节点的函数值.\x0d有限体积法(Finite Volume Method)又称为控制体积法.其基本思路是:将计算区域划分为一系列不重复的控制体积,并使每个网格点周围有一个控制体积;将待解的微分方程对每一个控制体积积分,便得出一组离散方程.其中的未知数是网格点上的因变量的数值.为了求出控制体积的积分,必须假定值在网格点之间的变化规律,即假设值的分段的分布的分布剖面.从积分区域的选取方法看来,有限体积法属于加权剩余法中的子区域法;从未知解的近似方法看来,有限体积法属于采用局部近似的离散方法.简言之,子区域法属于有限体积发的基本方法.\x0d有限体积法的基本思路易于理解,并能得出直接的物理解释.离散方程的物理意义,就是因变量在有限大小的控制体积中的守恒原理,如同微分方程表示因变量在无限小的控制体积中的守恒原理一样.限体积法得出的离散方程,要求因变量的积分守恒对任意一组控制体积都得到满足,对整个计算区域,自然也得到满足.这是有限体积法吸引人的优点.有一些离散方法,例如有限差分法,仅当网格极其细密时,离散方程才满足积分守恒;而有限体积法即使在粗网格情况下,也显示出准确的积分守恒.就离散方法而言,有限体积法可视作有限单元法和有限差分法的中间物.有限单元法必须假定值在网格点之间的变化规律(既插值函数),并将其作为近似解.有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化.有限体积法只寻求的结点值,这与有限差分法相类似;但有限体积法在寻求控制体积的积分时,必须假定值在网格点之间的分布,这又与有限单元法相类似.

有限元分析的节点和单元(有限元分析的节点和单元有哪些) 钢结构网架设计

有限元分析的节点和单元有限元分析的节点和单元(有限元分析的节点和单元有哪些)

阅读
分享