有限元分析第二讲
- 格式:pdf
- 大小:402.31 KB
- 文档页数:69
§5-5数值积分1、问题的提出在上一节中对等参元进行单元分析时要进行下列积分: (i) 单元刚度矩阵(ii)体积力的等效结点力(iii)边界力的等效结点力(iv)温升载荷的等效结点力式(5-4-5)~(5-4-8)分别归结为计算以下两种形式的积分对于上述积分仅在单元的形状十分规则的情况下才能得到解析的结果(精确值),一般情况只能用数值积分方法(主要是高斯求积法)求近似值。
虽然数值积分是“被迫“采用的,但后来发现:有选择地控制积分点的个数和位置,可以方便地实现我们的某些特殊意图。
这样一来,数值积分就成为有限元分析的一个重要组成部分,以至本来可以精确积分的三角形单元也常常采用数值积分。
2、数值积分的基本概念任何积分工作取决于三个要素:给定的积分区间,给定的被积函数,具体的积分方法。
下面以一维情况为例介绍数值积分的基本概念 (i) 梯形法函数()x f 在区间(a,b)的积分可以表达为 ()()ini ibax f W dx x f I ∑⎰=≈=1⎰⎰⎰---111111),()(dxdxy x f dx x f 、 [][][][][][][]ηξd d J t B E B tdxdyB E B k T Te det 1111⎰⎰⎰⎰--=={}[][]ηξσd d J t f f N td f f N r y x T y x T eV det 1111⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧=⎰⎰⎰⎰--{}[]{}ηξσγd Jd t B T det 01111T ⎰⎰--={}[]()()⎪⎪⎭⎫ ⎝⎛=⎭⎬⎫⎩⎨⎧=⎰⎰⎰--dy y f dx x f tds q p N r T 1111,ΓΓ(5-4-5) (5-4-8) (5-4-7) (5-4-6)i W :权系数; i x :积分样点;()i x f :积分样点的函数值。
梯形法的求积公式为其中,1--=n ab h ,而a b W ni i -=∑=1(ii) 当被积函数为n-1次多项式P n-1(x )时,则由n 个样点及其样点值(x i , P n-1(x i ),i=1,n )可以精确重构这个多项式,从而可以得到精确解。
第二章 有限元分析基本理论有限元法的基本思路是将一个连续求解区域分割成有限个不重叠且按一定方式相互连接在一起的子域(单元),利用在每一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数。
单元内的场函数通常由未知场函数或其导数在单元各个节点的数值和其插值函数来近似表示。
这样,未知场函数或其导数在各个节点上的数值即成为未知量(自由度)。
根据单元在边界处相互之间的连续性,将各单元的关系式集合成方程组,求出这些未知量,并通过插值函数计算出各个单元内场函数的近似值,从而得到全求解域上的近似解。
有限元将一个连续的无限自由度问题变成离散的有限自由度问题进行求解。
如果将区域划分成很细的网格,也即单元的尺寸变得越来越小,或随着单元自由度的增加及插值函数精度的提高,解的近似程度将不断被改进。
如果单元是满足收敛要求的,近似解最后可收敛于精确解。
2.1 有限元分析的基本概念和计算步骤首先以求解连续梁为例,引出结构有限元分析的一些基本概念和计算步骤。
如图2-1,连续梁承受集中力矩作用。
将结构离散为三个节点,两个单元。
结构中的节点编号为1、2、32.1.1单元分析在有限元分析过程中,第一步是进行结构离散,并对离散单元进行分析,分析的目的是得到单元节点的力与位移的关系。
单元分析的方法有直接法和能量法,本节采用直接法。
从连续梁中取出一个典型单元e ,左边为节点i ,右边为节点j 。
将节点选择在支承点处,单元两端只产生转角位移e i θ、ej θ,顺时针转动为正。
独立的单元杆端内力为弯矩i m 、j m ,顺时针为正。
记:{}e j i eu ⎭⎬⎫⎩⎨⎧=θθ为单元e 的节点位移向量;{}ej i em m f ⎭⎬⎫⎩⎨⎧=为单元e 的杆端力向量。
根据结构力学位移法可得如下平衡方程:⎪⎭⎪⎬⎫+=+=e j e e i e e j ej e e i e e i k k m k k m θθθθ22211211 (2-1)式中:ee e e ee i k k i k k 2412212211====,lEIi e =,EI 、l 分别为单元e 的抗弯刚度和长度。
第二章有限元分析基础有限元分析是一种常用的工程计算方法,在工程学科中被广泛应用。
本章将介绍有限元分析的基本概念和基础知识。
有限元分析是一种数值分析方法,用于求解复杂的物理问题。
它的基本思想是将一个连续的物体或结构离散化为有限数量的基本单元,通过在每个单元上进行计算,最终得到整个物体或结构的行为。
这些基本单元通过节点连接在一起,形成了一个有限元网格。
通过在每个节点上求解方程,可以得到整个物体或结构的应力、变形等相关信息。
在有限元分析中,有三个重要的步骤:建模、离散和求解。
建模是指将实际物体或结构转化为数学模型的过程。
在建模过程中,需要确定物体或结构的几何形状、边界条件和力学性质等。
离散是指将物体或结构划分为有限数量的基本单元。
常用的基本单元有三角形、四边形和六面体等。
离散过程中需要确定每个基本单元的几何属性和材料性质等。
求解是指在离散的基础上,通过求解节点上的方程,得到物体或结构的应力、变形等结果。
求解过程中,需要确定节点的位移和应变等参数。
有限元分析的基本假设是在每个基本单元内,应力和应变满足线性关系。
这意味着在小变形和小位移的情况下,有限元分析是有效的。
此外,为了提高计算精度,通常会增加更多的基本单元。
但是,增加基本单元数量会增加计算复杂度和计算时间。
因此,在实际应用中,需要根据问题的复杂程度和计算资源的限制进行权衡。
有限元分析广泛应用于各个领域,例如结构力学、热传导、电磁场、流体力学等。
在结构力学中,有限元分析可以用于求解静力学和动力学问题。
在热传导中,有限元分析可以用于求解温度分布和热流问题。
在电磁场中,有限元分析可以用于求解电荷和电场分布等。
在流体力学中,有限元分析可以用于求解流速和压力分布等。
总之,有限元分析是一种重要的工程计算方法,可以用于求解各种物理问题。
通过建模、离散和求解等步骤,可以得到物体或结构的应力、变形等结果。
有限元分析在工程学科中有着广泛的应用前景,对于工程设计和优化起着重要作用。
第二章有限元分析基本理论有限元分析是一种数值计算方法,广泛应用于结构分析、流体力学、热传导等工程领域。
它通过将连续的物理问题离散化为有限个简单的子问题,再通过数值方法求解这些子问题,最终得到原始问题的近似解。
有限元分析的基本理论包括三个方面:离散化、加权残差和求解方法。
首先是离散化。
离散化是指将原始的连续问题转化为离散的子问题。
有限元分析中常用的离散化方法是将求解区域分割成有限的子域,称为单元。
每个单元内部的场量(如位移、温度等)可以用其中一种函数近似表示。
离散化的关键是选择适当的单元形状和适量的节点,使得子问题的离散解能够较好地近似原问题的解。
接下来是加权残差方法。
加权残差方法是有限元分析的核心思想,用于构造子问题的弱型方程。
弱型方程是原始问题的一种积分形式,由应力平衡和边界条件推导而来。
在加权残差方法中,我们引入加权函数,将弱型方程乘以权函数,再对整个求解区域进行积分,从而将连续问题转化为离散问题。
通过选择合适的权函数,可以使得该离散问题具有良好的数学特性,比如对称、正定等。
最后是求解方法。
有限元分析的求解方法主要包括直接法和迭代法。
直接法适用于小型问题,通过对离散问题的系数矩阵进行直接求解,得到场量的离散解。
而迭代法适用于大型问题,通过迭代求解线性代数方程组,得到场量的近似解。
迭代法的常用算法有雅可比法、高斯-赛德尔法、共轭梯度法等。
在求解中还需要注意计算误差的控制和收敛性的判定。
除了这三个基本理论,有限元分析还有一些相关的概念和技术。
例如,网格生成用于生成离散化的单元网格;后处理用于对离散解进行可视化和数据分析;材料模型用于描述材料的本构关系。
这些概念和技术在具体的有限元分析应用中,有着重要的作用。
综上所述,有限元分析的基本理论包括离散化、加权残差和求解方法。
离散化将连续问题转化为离散子问题,加权残差方法用于构造子问题的弱型方程,求解方法用于求解离散问题。
掌握这些基本理论,对于理解和应用有限元分析方法具有重要意义。