有限元理论与技术-习题-有限元法.(优选)

有限元理论与技术-习题-有限元法.(优选)
有限元理论与技术-习题-有限元法.(优选)

填空题:

1、利用有限单元法求解弹性力学问题时,简单来说包含结构离散化、单元分析、整体分析三个主要步骤。

2、有限单元法首先将连续体变换成为离散化结构,然后再用结构力学位移法进行求解。其具体步骤分为单元分析和整体分析两部分。

3、每个单元的位移一般总是包含着两部分:一部分是由本单元的形变引起的,另一部分是由于其他单元发生了形变而连带引起的。

4、每个单元的应变一般总是包含着两部分:一部分是与该单元中各点的位置坐标有关的,是各点不相同的,即所谓变量应变;另一部分是与位置坐标无关的,是各点相同的,即所谓常量应变。

5、为了能从有限单元法得出正确的解答,位移模式必须能反映单元的刚体位移和常量应变,还应当尽可能反映相邻单元的位移连续性。

6、为了使得单元内部的位移保持连续,必须把位移模式取为坐标的单值连续函数,为了使得相邻单元的位移保持连续,就不仅要使它们在公共结点处具有相同的位移时,也能在整个公共边界上具有相同的位移。

7、在有限单元法中,单元的形函数N i在i 结点N i= 1 ;在其他结点N i= 0 及∑N i= 1 。

8、为了提高有限单元法分析的精度,一般可以采用两种方法:一是将单元的尺寸减小,以便较好地反映位移和应力变化情况;二是采用包含更高次项的位移模式,使位移和应力的精度提高。

9、在有限单元法中,结点力是指结点对单元的作用力。(√)

10、在平面三结点三角形单元的公共边界上应变和应力均有突变。(√)

11、形函数N i(xi,yi)= __(i=j)

N i(xi,yi)= __(i≠j)

简答题:

1、有限元分析的基本思路

答:首先,将物体或求解域离散为有限个互不重叠仅通过节点互相连接的子域(即单元),原始边界条件也被转化为节点上的边界条件,此过程称为离散化。

其次,在单元内,选择简单近似函数来分片逼近未知的求解函数,即分片近似。具体做法是在单元上选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各

变量或其导数的节点值与所选用的插值函数组成的线性表达式,这是有限元法的创意和精华所在。而整体区域上的解函数就是这些单元上的简单近似函数的组合。

最后,基于与原问题数学模型(基本方程和边界条件)等效的变分原理或加权残值法,建立有限元方程(即刚度方程),从而将微分方程转化为一组变量或其导数的节点值为未知量的代数方程组。从而借助矩阵表示和计算机求解代数方程组得到原问题的近似解。

2、简述有限元法求解中离散处理所遵循原则。

答:(1)几何逼真,(2)受力真实,(3)计算准确,(4)计算量少,(5)单元编号遵循右手准则(相邻单元编号差值最小)。

3、 针对附图所示的有限元结构,组集出整体刚度矩阵K 。(单元刚度矩阵用K e 表示,单元刚度矩阵元素用

e

ij

K 表示)。

单元刚度矩阵:

1 2 3

2 4 3

??????????=3332

31

232221131211k k k k k k k k k K ①

321

K ②=??????????3334

32

434442

232422

k k k k k k k k k 3

42

5 4

4 6 3 2

5 4

K ③=??????????3336

34

636664434644k k k k k k k k k 364 K ④=??????????4445

42

545552242522

k k k k k k k k k 4

52

单元刚度贡献矩阵:

1 2 3 4 5 6 1 2 3 4 5 6

K ①

=??

?

???

??

?

?

?????

?????0000

000000000000000000000k 33

32

31

232221

131211①①①①①①①①

①k k k k k k k k

65

4321

K ②

=???

???

???

???

?????

???000

000000000000000000000

0444342

3433

32242322k k k k k k k k k ②

6

5432

1

1 2 3 4 5 6 1 2 3 4 5 6

K ③

=?????????

?

?????

?????6664

63

4644433634330

00

00000

00000000000

0000000k k k k k k k k k 65

4321

K ④

=?????????

???

?????

???00

000000000

000000000000

0555452454442

252422k k k k k k k k k ④

6

54321

整体刚度矩阵:

1 2 3 4 5 6

K =?????????

??????

?????++++++++++++③③③④

③④④③②③②④②③②③②①②①①④④②②①④②①①①①

①6664

63

555452464544

444443

4342

42343433

3333323231

252424232322

222221

1312110

0000000000

k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 6

5432

1

3、用有限元法分析实际工程问题有哪些基本步骤?需要注意什么问题?

答:1.结构的离散化,2 单元分析2.1 选择位移函数2.2 载荷等效2.3 单元刚度矩阵3 整体分析3.1 集成等效节点载荷3.2 集成整体刚度 矩阵3.3 约束边界条件

1)建立实际工程问题的计算模型 利用几何、载荷的对称性简化模型

建立等效模型

2)选择适当的分析工具

侧重考虑以下几个方面:

多物理场耦合问题

大变形

网格重划分

3)前处理(Preprocessing)

建立几何模型(Geometric Modeling,自下而上,或基本单元组合)

有限单元划分(Meshing)与网格控制

4)求解(Solution)

给定约束(Constraint)和载荷(Load)

求解方法选择

计算参数设定

5)后处理(Postprocessing)

后处理的目的在于分析计算模型是否合理,提出结论。

用可视化方法(等值线、等值面、色块图)分析计算结果,包括位移、应力、应变、温度等;

最大最小值分析;

特殊部位分析。

4、在有限单元法中,位移模式应满足哪些基本条件。

答:(1)位移模式必须包含单元刚体位移;

(2)位移模式必须包含单元的常应变;

(3)位移模式在单元内要连续,且唯一在相邻单元之间要协调。

5、简述有限单元法结构刚度矩阵的特点。

答:(1)对称性;(2)奇异性;(3)主对角元恒正;

(4)稀疏性;(5)非零元素带状分布。

6、简述有限单元法中单元刚度矩阵的性质。

答:(1)单元刚度矩阵为对称矩阵;

(2)单元刚度矩阵为奇异矩阵;

(3)单元刚度矩阵主对角线元素恒为正值;

(4)单元刚度矩阵仅与单元本身有关。

7、简述有限单元法中形函数的性质。

答:(1)形函数N i在节点i处的值为1,在其他两个节点j,m处的值为0;

(2)在单元上任意一点处,3个形函数的和都等于1。

8、简述有限元法中选取单元位移函数(多项式)的一般原则。

答:一般原则:

(1)广义坐标的个数应该与结点自由度数相等;

(2)选取多项式时,常数项和坐标的一次项必须完备;

(3)多项式的选取应由低阶到高阶;

(4)尽量选取完全多项式以提高单元的精度。

9、要保证有限单元法计算结果的收敛性,位移函数必须满足那些条件?

答:(1)位移函数必须能够反映单元的常量应变;

(2)位移函数必须能够反映单元的刚性位移;

(3)位移函数在单元内部必须是连续函数;

(4)位移函数必须保证相邻单元间唯一协调。

10、为了保证有限单元法解答的收敛性,位移模式应满足哪些条件?

答:为了保证有限单元法解答的收敛性,位移模式应满足下列条件:

(1)位移模式必须能反映单元的刚体位移;

(2)位移模式必须能反映单元的常量应变;

(3)位移模式应尽可能反映位移的连续性。

11、以三节点三角形单元为例,简述有限单元法求解离散化结构的具体步骤。

答:(1)取三角形单元的结点位移为基本未知量。

(2)应用插值公式,由单元的结点位移求出单元的位移函数。

(3)应用几何方程,由单元的位移函数求出单元的应变。

(4)应用物理方程,由单元的应变求出单元的应力。

(5)应用虚功方程,由单元的应力出单元的结点力。

(6)应用虚功方程,将单元中的各种外力荷载向结点移置,求出单元的结点荷载。

(7)列出各结点的平衡方程,组成整个结构的平衡方程组。

12、在有限单元法中,为什么要求位移模式必须能反映单元的刚体位移?

答:每个单元的位移一般总是包含着两部分:一部分是由本单元的形变引起的,另一部分是

本单元的形变无关的,即刚体位移,它是由于其他单元发生了形变而连带引起的。甚至在弹性体的某些部位,例如在靠近悬臂梁的自由端处,单元的形变很小,单元的位移主要是由于其他单元发生形变而引起的刚体位移。因此,为了正确反映单元的位移形态,位移模式必须能反映该单元的刚体位移。

13、在有限单元法中,为什么要求位移模式必须能反映单元的常量应变?

答:每个单元的应变一般总是包含着两部分:一部分是与该单元中各点的位置坐标有关的,是各点不相同的,即所谓变量应变;另一部分是与位置坐标无关的,是各点相同的,即所谓常量应变。而且,当单元的尺寸较小时,单元中各点的应变趋于相等,也就是单元的应变趋于均匀,因而常量应变就成为应变的主要部分。因此,为了正确反映单元的形变状态,位移模式必须能反映该单元的常量应变。

14、采用有限单元法怎么样求解弹性力学问题(基本思路和基础步骤)?

答:基本思想:

根据近似分割和能量极值原理,把求解区域离散为有限个单元的组合,研究每个单元的特性,组装各单元,通过变分原理,把问题化成线性代数方程组求解。

基本步骤:

(1)将结构进行离散化,包括单元划分、结点编号、单元编号、结点坐标计算、位移约束条件确定;

(2)等效结点力的计算;

(3)刚度矩阵的计算(先逐个计算单元刚度,再组装成整体刚度矩阵);

(4)建立整体平衡方程,引入约束条件,求解结点位移;

(5)应力计算。

15、简述有限单元法常分析的问题。

答:(1)静力学分析,(2)模态分析,(3)动力学分析,(4)热力学分析,(5)其他:接触分析、压杆稳定性分析、结构-流体耦合分析等。

16、有限元法的单元分析中如何导出单元刚度矩阵?

17、在有限元方法中如何形成整体刚度矩阵?如何进行载荷的等效处理?

18、单元节点数目、位移函数、计算精度三者之间存在何种关系?

19、什么是等参变换单元?为什么要采用等参变化?

答:(1)是在取位移模式和进行图形变换时所取变换参数相同,叫做等参变换单元

(2) 等参变换是为了可以把任意四边形和正方形单元联系起来。由于等参变换的采用使等参

单元的刚度、质量、阻尼、荷载等特性矩阵的计算仍在前面所表示单元的规则域内进行,因此不管各个积分形式的矩阵表示的被积函数如何复杂,仍然可以方便地采用标准化的数值积分方法计算。

20、如何采用有限元法求解动力学方程?

21、为什么说单元位移函数是有限元法的关键?单元位移函数收敛的条件是什么?

最新文件---------------- 仅供参考--------------------已改成-----------word文本

--------------------- 方便更改

赠人玫瑰,手留余香。

有限元基础知识归纳

有限元知识点归纳 1.、有限元解的特点、原因? 答:有限元解一般偏小,即位移解下限性 原因:单元原是连续体的一部分,具有无限多个自由度。在假定了单元的位移函数后,自由度限制为只有以节点位移表示的有限自由度,即位移函数对单元的变形进行了约束和限制,使单元的刚度较实际连续体加强了,因此,连续体的整体刚度随之增加,离散后的刚度较实际的刚度K为大,因此求得的位移近似解总体上将小于精确解。 2、形函数收敛准则(写出某种单元的形函数,并讨论收敛性)P49 (1)在节点i处N i=1,其它节点N i=0; (2)在单元之间,必须使由其定义的未知量连续; (3)应包含完全一次多项式; (4)应满足∑Ni=1 以上条件是使单元满足收敛条件所必须得。可以推证,由满足以上条件的形函数所建单元是完备协调的单元,所以一定是收敛的。 4、等参元的概念、特点、用时注意什么?(王勖成P131) 答:等参元—为了将局部坐标中几何形状规则的单元转换成总体(笛卡尔)坐标中的几何形状扭曲的单元,以满足对一般形状求解域进行离散化的需要,必须建立一个坐标变换。即: 为建立上述的变换,最方便的方法是将上式表示成插值函数的形式,即: 其中m是用以进行坐标变换的单元节点数,xi,yi,zi是这些结点在总体(笛卡尔)坐标内的坐标值,Ni’称为形状函数,实际上它也是局部坐标表示的插值函数。称前者为母单元,后者为子单元。 还可以看到坐标变换关系式和函数插值表示式:在形式上是相同的。如果坐标变换和函数插值采用相同的结点,并且采用相同的插值函数,即m=n,Ni’=Ni,则称这种变换为等参变换。 5、单元离散?P42 答:离散化既是将连续体用假想的线或面分割成有限个部分,各部分之间用有限个点相连。每个部分称为一个单元,连接点称为结点。对于平面问题,最简单、最常用的离散方式是将其分解成有限个三角形单元,单元之间在三角形顶点上相连。这种单元称为常应变三角形单元。常用的单元离散有三节点三角形单元、六节点三角形单元、四节点四边形单元、八节点四边形单元以及等参元。 6、数值积分,阶次选择的基本要求? 答:通常是选用高斯积分 积分阶次的选择—采用数值积分代替精确积分时,积分阶数的选取应适当,因为它直接影响计算精度,计算工作量。选择时主要从两方面考虑。一是要保证积分的精度,不损失收敛性;二是要避免引起结构总刚度矩阵的奇异性,导致计算的失败。

有限元理论基础

有限元理论基础

有限元理论基础 2.1 数值模拟技术 2.1.1数值模拟技术简介 在工程技术领域中许多力学问题和场问题,实质上就是在一定的边界条件下求解一些微分方程。对于少数简单问题,人们可以通过建立它们的微分方程与边界约束求出该问题的解析解。但是对于比较复杂的数学方程问题以及不规则的边界条件通过激吻戏法往往难以求解,而需要借助各种数值模拟方法活的相应的工程数值解,这就是所谓的数值模拟技术。 在实际工程领域中,用数值模拟技术可以对复杂的工程结构进行受力和响应分析,这样可以在设计或者加工前预知实体结构工作状态下的大概情况。 目前在工程实际应用中,常用的数值求解方法有:有限单元法、有限差分法、边界元等但从实用性和使用范围来说,有限单元法则是随着计算机技术的发展而被广泛应用的一种行之有效的数值计算方法。 2.2.2 有限元法 有限元法是一种基于能量原理的数值计算

方法,是解决工程实际问题的一种有效的数值计 算工具。它是里茨法的另一种表示形式,它可应用里茨法分析的所有弹性理论。 限元法是处理连续的结构体离散或有限个单元集合,也就是将连续的求解域离散为一定数量的单元集合体。且每个单元都具有一定的节点,相邻单元通过节点相互连续,同时使用等效节点力代替作用于单元上的力和选定场函数的节点值作为基本未知量。并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律:进而利用力学中的某些变分原理去建立用以求解节点未知量的有限元法方程,从而将一个连续域中的无限自由度问题化为离散域中的有限自由度问题。求解后,可利用解出的节点值和设定的插值函数确定整个单元集体上的场函数。有限元求解问题中的单元分析:t t t a k F= 式中::t F单元节点作用力。 t K:单元刚度矩阵。 t a:单元节点位移。 通过单元分析确定单元刚度矩阵,建立单元节点作用力和单元为伊关系。有限元求解问题时建立 的结构整体平衡方程:P KU=

有限元理论方法

关于有限元分析法及其应用举例 摘要:本文主要介绍有限元分析法,作为现代设计理论与方法的一种,已经在 众多领域普遍使用。介绍了它的起源和国内外发展现状。阐述了有限元法的基 本思想和设计方法。并从实际出发,例举了有限元法的一个简单应用———啤 酒瓶的应力分析和优化,表明了利用有限元分析法的众多优点。随着计算机的 发展,基于有限元分析方法的软件开发越来越多。本文也在其软件开发方面进 行阐述,并简单介绍了一下主流软件的发展情况和使用范围。并就这一领域的 未来发展趋势进行阐述。 关键词:有限元分析法软件啤酒瓶 Abstract:This thesis mainly introduces the finite element analysis, as a modern design theory and methods used widely in in most respects. And this paper introduces its origins and development in world. It also expounds the basic thinking and approach of FEM..Proceed from the actual situation,this text holds the a simple application of finite-element method———the analysis and optimized of an beer bottle and indicate the the numerous benefits of finite element analysis .As computers mature and based on the finite element analysis of the software development is growing. This article introduces its application in the software development aspects as well, and briefly states the development and scope of the mainstream software. And it’s also prospect future development tendency in this area . Key: Finite Element Analysis Software Beer bottle 0 绪论 有限元法(Finite Element Method,FEM),是计算力学中的一种重要的方法,它是20世纪50年代末60年代初兴起的应用数学、现代力学及计算机科学相互渗透、综合利用的边缘科学。有限元法最初应用在工程科学技术中,用于模拟并且解决工程力学、热学、电磁学等物理问题。对于过去用解析方法无法求解的问题和边界条件及结构形状都不规则的复杂问题,有限元法则是一种有效的分析方法。有限元法的基本思想是先将研究对象的连续求解区域离散为一组有限个且按一定方式相互联结在一起的单元组合体。由于单元能按不同的联结方式进行组合,且单元本身又可以有不同形状,因此可以模拟成不同几何形状的求解小区域;

有限元分析基本理论问答 基础理论知识

1. 诉述有限元法的定义 答:有限元法是近似求解一般连续场问题的数值方法 2. 有限元法的基本思想是什么 答:首先,将表示结构的连续离散为若干个子域,单元之间通过其边界上的节点连接成组合体。其次,用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。 3. 有限元法的分类和基本步骤有哪些 答:分类:位移法、力法、混合法;步骤:结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移。 4. 有限元法有哪些优缺点 答:优点:有限元法可以模拟各种几何形状复杂的结构,得出其近似解;通过计算机程序,可以广泛地应用于各种场合;可以从其他CAD软件中导入建好的模型;数学处理比较方便,对复杂形状的结构也能适用;有限元法和优化设计方法相结合,以便发挥各自的优点。 缺点:有限元计算,尤其是复杂问题的分析计算,所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。对无限求解域问题没有较好的处理办法。尽管现有的有限元软件多数使用了网络自适应技术,但在具体应用时,采用什么类型的单元、多大的网络密度等都要完全依赖适用者的经验。 5. ?梁单元和平面钢架结构单元的自由度由什么确定 答:每个节点上有几个节点位移分量,就称每个节点有几个自由度 6. ?简述单元刚度矩阵的性质和矩阵元素的物理意义 答:单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵 单元刚度矩阵中元素aml的物理意义为单元第L个节点位移分量等于1,其他节点位移分量等于0时,对应的第m个节点力分量。 7. 有限元法基本方程中的每一项的意义是什么 答:整个结构的节点载荷列阵(外载荷、约束力),整个结构的节点位移列阵,结构的整体刚度矩阵,又称总刚度矩阵。 8. 位移边界条件和载荷边界条件的意义是什么 答:由于刚度矩阵的线性相关性不能得到解,从而引入边界条件。 9. ?简述整体刚度矩阵的性质和特点 答:对称性;奇异性;稀疏性;对角线上的元素恒为正。 11. 简述整体坐标的概念 答:单元刚度矩阵的坐标变换式把平面刚架的所有单元在局部坐标系X’Y’Z’下的单元刚度矩阵变换到一个统一的坐标系xOy下,这个统一的坐标系xOy称为整体坐标系。 13. 简述平面钢架问题有限元法的基本过程 答:力学模型的确定,结构的离散化,计算载荷的等效节点力,计算各单元的刚度矩阵,组集整体刚度矩阵,施加边界约束条件,求解降价的有限元基本方程,求解单元应力,计算结果的输出。 14. 弹性力学的基本假设是什么。 答:连续性假定,弹性假定,均匀性和各向同性假定,小变形假定,无初应力假定。 15.弹性力学和材料力学相比,其研究方法和对象有什么不同。 答:研究对象:材料力学主要研究杆件,如柱体、梁和轴,在拉压、剪切、弯曲和扭转等作用下的应力、形变和位移。弹性力学研究各种形状的弹性体,除杆件外,还研究平面体、空间体,板和壳等。因此,弹性力学的研究对象要广泛得多。研究方法:弹性力学和材料力学

有限元分析理论基础

有限元分析概念 有限元法:把求解区域看作由许多小的在节点处相互连接的单元(子域)所构成,其模型给出基本方程的分片(子域)近似解,由于单元(子域)可以被分割成各种形状和大小不同的尺寸,所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型:它是真实系统理想化的数学抽象。由一些简单形状的单元组成,单元之间通过节点连接,并承受一定载荷。 有限元分析:是利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。并利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的,所考虑的变形建立在小变形假设的基础上。在这类问题中,材料的应力与应变呈线性关系,满足广义胡克定律;应力与应变也是线性关系,线弹性问题可归结为求解线性方程问题,所以只需要较少的计算时间。如果采用高效的代数方程组求解方法,也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别: 1)非线性问题的方程是非线性的,一般需要迭代求解; 2)非线性问题不能采用叠加原理; 3)非线性问题不总有一致解,有时甚至没有解。 有限元求解非线性问题可分为以下三类:

1)材料非线性问题 材料的应力和应变是非线性的,但应力与应变却很微小,此时应变与位移呈线性关系,这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系,所以,一般材料的应力与应变之间的非线性关系要基于试验数据,有时非线性材料特性可用数学模型进行模拟,尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有:非线性弹性(包括分段线弹性)、弹塑性、粘塑性及蠕变等。 2)几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时,应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题,橡胶部件形成过程为大应变问题。 3)非线性边界问题 在加工、密封、撞击等问题中,接触和摩擦的作用不可忽视,接触边界属于高度非线性边界。 平时遇到的一些接触问题,如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等,当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。

有限元分析理论基础

有限元分析概念 有限元法:把求解区域瞧作由许多小的在节点处相互连接的单元(子域)所构成,其模型给出基本方程的分片(子域)近似解,由于单元(子域)可以被分割成各种形状与大小不同的尺寸,所以它能很好地适应复杂的几何形状、复杂的材料特性与复杂的边界条件 有限元模型:它就是真实系统理想化的数学抽象。由一些简单形状的单元组成,单元之间通过节点连接,并承受一定载荷。 有限元分析:就是利用数学近似的方法对真实物理系统(几何与载荷工况)进行模拟。并利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元就是以理想弹性体为研究对象的,所考虑的变形建立在小变形假设的基础上。在这类问题中,材料的应力与应变呈线性关系,满足广义胡克定律;应力与应变也就是线性关系,线弹性问题可归结为求解线性方程问题,所以只需要较少的计算时间。如果采用高效的代数方程组求解方法,也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别: 1)非线性问题的方程就是非线性的,一般需要迭代求解; 2)非线性问题不能采用叠加原理; 3)非线性问题不总有一致解,有时甚至没有解。 有限元求解非线性问题可分为以下三类:

1)材料非线性问题 材料的应力与应变就是非线性的,但应力与应变却很微小,此时应变与位移呈线性关系,这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系,所以,一般材料的应力与应变之间的非线性关系要基于试验数据,有时非线性材料特性可用数学模型进行模拟,尽管这些模型总有她们的局限性。在工程实际中较为重要的材料非线性问题有:非线性弹性(包括分段线弹性)、弹塑性、粘塑性及蠕变等。 2)几何非线性问题 几何非线性问题就是由于位移之间存在非线性关系引起的。 当物体的位移较大时,应变与位移的关系就是非线性关系。研究这类问题一般都就是假定材料的应力与应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题,橡胶部件形成过程为大应变问题。 3)非线性边界问题 在加工、密封、撞击等问题中,接触与摩擦的作用不可忽视,接触边界属于高度非线性边界。 平时遇到的一些接触问题,如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等,当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。 有限元理论基础

有限元理论与方法

第一章 绪论 有限元发展过程: 有限元法在西方起源于收音机和导弹的结构设计,发表这方面文章最早而且最有影响的是西德教授,于1954—1955年间分阶段在《Aircraft Engineering 》上发表上许多有关这方面的论文,并在此基础上写成了《能量原理与结构分析》,此书内容提供了有限元法的理论基础。美国的、 、 和等人于1956年发表了了篇题为《复杂结构的刚度和挠度分析》一文,此文提出了计算复杂结构刚度影响系数的方法,并说明了如何利用计算机进行分析。美国于1960年在一篇介绍平面应力分析的论文中,首先提出了有限元的名字。1965年英国及其合作者解决了将有限元法应用于所有场的问题,使有限元法的应用更加广泛。 有限元法的基本思路: 有限元法的基本思路和基本原理以结构力学中的位移法为基础,把复杂的结构或连续体看成为有限个单元的组合,各单元彼此在节点处连续而组成整体,把连续体分成有限个单元和节点,称之为离散化,先对单元进行特性分析,然后根据各单元在节点处的平衡协调条件建立方程,综合后作整体分析。 这样一分一合,先离散再综合的过程,就把复杂结构或连续体的计算问题转化为简单单元的分析与综合问题。 有限元分析中可采取三种方法: 位移法——取节点位移作为基本未知数 力 法——取节点力作为基本未知数 混合法—— 有限元法分析过程: 1、结构离散化(单元划分) 2、选择位移模式 为了能用节点位移表示单元体的位移、应变和应力,在分析连续体时,必须对单元中位移的分布做出一定的假定,也就是假定位移是坐标的某种简单函数,这种函数称为位移模式或位移函数(形函数)。 {}[]{}e u N δ= (1) 3、分析单元的力学特性 (1)利用几何方程:由位移表达式导出用点位移表示单元应变的关系式 {}[]{} e εδ=B {}ε为单元内任一点的应变列阵 (2) 非线性有限元 线性有限元 几何非线性 材料非线性 有限元

有限元法中的几个基本概念

诚信·公平·开放·共赢 Loyalty Fair Opening Win-win 有限元法中的几个基本概念 有限元法是把要分析的连续体假想地分割成有限个单元所组成的组合体,简称离散化。 这些单元仅在顶角处相互联接,称这些联接点为结点。 离散化的组合体与真实弹性体的区别在于:组合体中单元与单元之间的联接除了结点之外再无任何关联。但是这种联接要满足变形协调条件,即不能出现裂缝,也不允许发生重叠。显然,单元之间只能通过结点来传递内力。 通过结点来传递的内力称为结点力,作用在结点上的荷载称为结点荷载。当连续体受到外力作用发生变形时,组成它的各个单元也将发生变形,因而各个结点要产生不同程度的位移,这种位移称为结点位移。 在有限元中,常以结点位移作为基本未知量。并对每个单元根据分块近似的思想,假设一个简单的函数近似地表示单元内位移的分布规律,再利用力学理论中的变分原理或其他方法,建立结点力与位移之间的力学特性关系,得到一组以结点位移为未知量的代数方程,从而求解结点的位移分量。然后利用插值函数确定单元集合体上的场函数。显然,如果单元满足问题的收敛性要求,那么随着缩小单元的尺寸,增加求解区域内单元的数目,解的近似程度将不断改进,近似解最终将收敛于精确解。 附:FELAC 2.0软件简介 FELAC 2.0采用自定义的有限元语言作为脚本代码语言,它可以使用户以一种类似于数学公式书写和推导的方式,非常自然和简单的表达待解问题的微分方程表达式和算法表达式,并由生成器解释产生完整的并行有限元计算C程序。 FELAC 2.0的目标是通过输入微分方程表达式和算法之后,就可以得到所有有限元计算的程序代码,包含串行程序和并行程序。该系统采用一种语言(有限元语言)和四种技术(对象技术、组件技术、公式库技术生成器技术)开发而成。并且基于FELAC 1.0的用户界面,新版本扩充了工作目录中右键编译功能、命令终端输入功能,并且丰富了文本编辑功能,改善了用户的视觉体验,方便用户快速便捷的对脚本或程序进行编辑、编译与调试。其中并行版在前后处理上进行了相应的改进。

有限元分析中的一些问题

有限元分析的一些基本考虑-—-—-单元形状对于计算精度的影响 笔者发现,在分析复杂问题时,我们所可能出现的错误,竟然是一些很根本的错误,这些根本错误是由于对有限元的基本理论理解不清晰而造成的。 鉴于这个原因,笔者决定对一些基本问题(例如单元形状问题,单元大小问题,应力集中问题等)展开调查,从而形成了一系列文章,本篇文章是这些系列文章中的第一篇. 本篇文章先考虑有限元分析中的第一个基本问题:单元形状问题。 我们知道,单元形状对于有限元分析的结果精度有着重要影响,而对单元形状的衡量又有着诸多指标,为便于探讨,这里首先只讨论第一个最基本的指标:长宽比(四边形单元的最长尺度与最短尺度之比),而且仅考虑平面单元的长宽比对于计算精度的影响。 为此,我们给出一个成熟的算例。该算例是一根悬臂梁,在其端面施加竖直向下的抛物线分布载荷,我们现在考察用不同尺度的单元划分该梁时,对于A点位移的影响。 这五种不同的划分方式,都使用矩形单元,只不过各单元的长宽比不同。 例如第一种(1)AR=1.1,就是长宽比接近1; 第二种(2)AR=1.5,就是长宽比是1。5.其它类推。 第五种(5)AR=24,此时单元的长度是宽度的24倍。 现在我们看看按照这五种单元划分方式对于A点位移的影响,顺便我们也算出了B点的位移,结果见下表.

我们现在仔细查看一下上表,并分析其含义。 我们先考虑第一行,它是第一种单元划分情况,此时每个单元的长宽比是1。1,由此我们计算出A点,B点的垂直位移,可以看到,A点的竖直位移是—1.093英寸,而B点的竖直位移是-0。346英寸。而这两点我们都是可以用弹性力学的方式得到精确解的,其精确解分别是-1。152以及—0。360。这样,我们可以得到此时A点位移误差的百分比是[(—1.093)—(-1。152)]/1。152 =5。2%. 对于其它情况,也采用类似的方式得到A点位移误差的百分比。 从上表可以看出来,随着长宽比的增加,位移误差越来越大,竟然大到56%.因此,如果我们是用长宽比为24的单元进行划分的话,那么我们的结果可以说是完全错误的. 下面按照上表绘制出一张图,该图从形象的角度表达了上表的含义.

有限元分析及应用大作业

有限元分析及应用大作业 作业要求: 1)个人按上机指南步骤至少选择习题中3个习题独立完成,并将计算结果上交; 也可根据自己科研工作给出计算实例。 2)以小组为单位完成有限元分析计算; 3)以小组为单位编写计算分析报告; 4)计算分析报告应包括以下部分: A、问题描述及数学建模; B、有限元建模(单元选择、结点布置及规模、网格划分方案、载荷及边界 条件处理、求解控制) C、计算结果及结果分析(位移分析、应力分析、正确性分析评判) D、多方案计算比较(结点规模增减对精度的影响分析、单元改变对精度的 影响分析、不同网格划分方案对结果的影响分析等) 题一:图示无限长刚性地基上的三角形大坝,受齐顶的水压力作用,试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析,并对以下几种计算方案进行比较: 1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算;(注意ANSYS中用四边形单元退化为三节点三角形单元) 2)分别采用不同数量的三节点常应变单元计算; 3)当选常应变三角单元时,分别采用不同划分方案计算。 解:1.建模: 由于大坝长度>>横截面尺寸,且横截面沿长度方向保持不变,因此可将大坝看作无限长的实体模型,满足平面应变问题的几何条件;对截面进行受力分析,作

用于大坝上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力,满足平面应变问题的载荷条件。因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况P=98000-9800*Y;建立几何模型,进行求解;假设大坝的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.3; 2:有限元建模过程: 2.1 进入ANSYS : 程序→ANSYS APDL 15.0 2.2设置计算类型: ANSYS Main Menu: Preferences →select Structural →OK 2.3选择单元类型: ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182,六节点三角形单元选择Solid Quad 8node 183)→OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window) 2.4定义材料参数: ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK 2.5生成几何模型: 生成特征点: ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS →依次输入四个点的坐标:input:1(0,0),2(10,0),3(1,5),4(0.45,5) →OK 生成坝体截面: ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点,1(0,0),2(6,0),3(0,10) →OK 2.6 网格划分: ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window) 2.7 模型施加约束: 给底边施加x和y方向的约束: ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK 给竖直边施加y方向的分布载荷: ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量;2) 在Result窗口中出现{X},写入所施加的载荷函数: 98000-9800*{Y};3) File>Save(文件扩展名:func) →返回:Parameters →Functions →Read from file:将需要的.func文件打开,参数名取meng,它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边;OK →在下拉列表框中,选择:Existing table →OK →选择需要的载荷为meng参数名→OK 2.8 分析计算: ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load

有限元分析的概念和理论

第五章有限元素方法

§5.1有限元素方法的基本思想 有限元素法是一套求解微分方程的系统化数值计算方法。它比传统解法具有理论完整可靠,物理意义直观明确,适应性强,形式单纯、规范,解题效能强等优点。 从数学上来说, 有限元素方法是基于变分原理。它不象差分法那样直接去解偏微分方程, 而是求解一个泛函取极小值的变分问题。有限元素法是在变分原理的基础上吸收差分格式的思想发展起来的。 采用有限元素法还能使物理特性基本上被保持, 计算精度和收敛性进一步得到保证。 有限元素法优点: - 降低实验所需成本 - 減少試验对象的变异困难 - 方便参数控制 - 可获得实验无法获得的信息

有限元素法基本概念: 元素(element),节点(node),连結元素 有限元素法的基本思想: ?实际的物理問題很难利用单一的微分方程式描述,更无法順利求其解析解. ?有限元素法是将复杂的几何外型結构的物体切割成许多简单的几何形状称之为元素. ?元素与与元素间以“节点”相连. ?由于元素是简单的几何形状,故可以順利地写出元素的物理方程式,並求得节点上的物理量. ?采用內插法求得元素內任意点的物理量.

§5.2二维场的有限元素方法 1. 场域划分的约定 三角形元素。三角形元素越小,场域的分割就越细,计算的精度就会越高。因而在实际应用中是按精度的要求来决定场域内各处三角形元素的大小。 一般规定每个三角形元素的三个边的边长尽量地接近,尽量避免三角形元素具有大的钝角,一般最长的一条边不得大于最短边的三倍。 在分割场域时要求各三角形元素之间只能以顶点相交,即两相邻的三角形元素有两个公共的顶点及一条等长的公共边。不能把一个三角形的顶点取在另一个三角形的边上。 划分时还应当注意要尽量地使由相邻边界节点之间的线段所近似构成的曲线足够光滑。 如果在场域D内有不同的介质,则需要将介质的交面线选为分割线。

有限元法的理论基础

有限元法的理论基础 有限元法是一种离散化的数值计算方法,对于结构分析而言,它的理论基础是能量原理。能量原理表明,在外力作用下,弹性体的变形、应力和外力之间的关系受能量原理的支配,能量原理与微分方程和定解条件是等价的。下面介绍有限元法中经常使用的虚位移原理和最小势能原理。 1.虚位移原理 虚位移原理又称虚功原理,可以叙述如下:如果物体在发生虚位移之前所受的力系是平衡的(物体内部满足平衡微分方程,物体边界上满足力学边界条件),那么在发生虚位移时,外力在虚位移上所做的虚功等于虚应变能(物体内部应力在虚应变上所做的虚功)。反之,如果物体所受的力系在虚位移(及虚应变)上所做的虚功相等,则它们一定是平衡的。可以看出,虚位移原理等价于平衡微分方程与力学边界条件。所以虚位移原理表述了力系平衡的必要而充分的条件。 虚位移原理不仅可以应用于弹性性力学问题,还可以应用于非线性弹性以及弹塑性等非线性问题。 2.最小势能原理 最小势能原理可以叙述为:弹性体受到外力作用时,在所有满足位移边界条件和变形协调条件的可以位移中,真实位移使系统的总势能取驻值,且为最小值。根据最小势能原理,要求弹性体在外力作用下的位移,可以满足几何方程和位移边界条件且使物体总势能取最小值的条件去寻求答案。最小势能原理仅适用于弹性力学问题。 2.2有限元法求解问题的基本步骤 弹性力学中的有限元法是一种数值计算方法,对于不同物理性质和数学模型的问题,有限元法的基本步骤是相同的,只是具体方式推导和运算求解不同,有限元求解问题的基本步骤如下。 2.2.1问题的分类 求解问题的第一步就是对它进行识别分析,它包含的更深层次的物理问题是什么?比如是静力学还是动力学,是否包含非线性,是否需要迭代求解,要从分析中得等到什么结果等。对这些问题的回答会加深对问题的认识与理解,直接影响到以后的建模与求解方法的选取等。 2.2.2建模 在进行有限元离散化和数值求解之值,我们为分析问题设计计算模型,这一步包括决定哪种特征是所要讨论的重点问题,以便忽略不必要的细节,并决定采用哪种理论或数学公式描述结果的行为。因此,我们可以忽略几何不规则性,把一些载荷看做是集中载荷,并把某些支撑看做是固定的。材料可以理想化为线弹性和各向同性的。根据问题的维数、载荷以及理论化的边界条件,我们能够决定采用梁理论、板弯曲理论、平面弹性理论或者一些其他分析理论描述结构性能。在求解中运用分析理论简化问题,建立问题的模型。 2.2.3连续体离散化 连续体离散化,习惯上称为有限元网络划分,即将连续体划分为有限个具有规则形状的单元的集合,两相邻单元之间只通过若干点相互连接,每个连接点称为节点。单元节点的设置、性质、数目等应视问题的性质、描述变形的需要和计算精度而定,如二维连续体的单元可为三角形、四边形,三维连续体的单元可以是四面体、长方体和六面体等。为合理有效地表示连续体,需要适当选择单元的类型、数目、大小和排列方式。 离散化的模型与原来模型区别在于,单元之间只通过节点相互连接、相互作用,而无其他连接。因此这种连接要满足变形协调条件。离散化是将一个无限多自由度的连续体转化为一个有限多自由度的离散体过程,因此必然引起误差。主要有两类:建模误差和离散化误差。

有限元理论与技术-习题-弹性力学DOC

弹性力学 填空题: 1、连续体力学包括固体力学、流体力学、热力学和电磁场力学,非连续体力学包括量子力学。 2、弹性力学所研究的范围属于固体力学中弹性阶段。 3、弹性力学的基本假定为:连续性、完全弹性、均匀性和各向同性、变形很小、无初应力。 4、连续性假设是指:物体内部由连续介质组成,物体中应力、应变和位移分量为连续的,可用连续函数表示。 5、均匀性和各向同性假设是指:物体内各点和各方向的介质相同,即物理性质相同,物体的弹性常数杨氏模量和泊松比不随坐标和方向的变化而变化。 6、完全弹性假设是指:物体在外载荷作用下发生变形,在外载荷去除后,物体能够完全恢复原形,材料服从胡克定律,即应力与形变成正比。 7、在弹性力学里分析问题,要考虑静力学、几何学和物理学三方面条件,分别建立三套方程为:平衡方程、几何方程和物理方程,三组方程分别表示:应力与载荷关系、应变与位移关系、应力与应变关系。 8、弹性力学研究弹性体由于受外力作用、边界约束或温度改变等原因而发生的应力、形变和位移。 9、在弹性力学中规定,线应变以伸长时为正,缩短时为负,与正应力的正负号规定相适应。 10、在弹性力学中规定,切应变以直角变小时为正,变大时为负,与切应力的正负号规定相适应。 11、物体受外力以后,其内部将发生内力,它的集度称为应力。与物体的形变和材料强度直接有关的,是应力在其作用截面的法线方向和切线方向的分量,也就是正应力和切应力。

12、建立平衡方程时,在正六面微分体的6个面上共有9个应力分量,分别为:,其中正应力为:,剪应力为:,这些应力分量与外载荷共同建立 3 个方程。 13、建立几何方程时,线应变为,角应变为,这些应变与位移共同建立6 个方程。 14、物理方程表示应力与应变的关系,即为胡克定律,其中弹性常数E和μ分别表示材料的杨氏模量和泊松比,物理方程组共包含 6 个方程。 15、平面问题分为平面应力问题和平面应变问题,两者所研究得对象分别为等厚度薄平板和等截面长柱体。 16、平面应力问题和平面应变问题基本方程中:平衡方程和几何方程相同,物理方程不相同。(相同或不相同) 17、表示应力分量与体力分量之间关系的方程为平衡微分方程。 15、边界条件表示边界上位移与约束,或应力与面力之间的关系式。 18、按应力求解平面问题时常采用逆解法和半逆解法。 19、弹性力学中边界条件通常可以分为:位移边界条件、应力边界条件和混合边界条件。 20、弹性力学问题的解法分为解析法、变分法和差分法,就解题方法而言,又分为如下两种方法:位移法和应力法。 21、将平面应力情况下的物理方程中的弹性模量E,泊松比 分别换成及就要得到平面应变情况下相应的物理方程。 22、位移法为物理方程与几何方程联立消除应变分量,得到应力与位移的函数方程式,再与平衡方程联立消除应力,得到载荷与位移的方程式。简答题: 1、在弹性力学中根据什么分别推导出平衡微分方程、几何方程、物理方程,这三个方程分别表示什么关系?

有限元理论基础

有限元理论基础 2.1 数值模拟技术 2.1.1数值模拟技术简介 在工程技术领域中许多力学问题和场问题,实质上就是在一定的边界条件下求解一些微分方程。对于少数简单问题,人们可以通过建立它们的微分方程与边界约束求出该问题的解析解。但是对于比较复杂的数学方程问题以及不规则的边界条件通过激吻戏法往往难以求解,而需要借助各种数值模拟方法活的相应的工程数值解,这就是所谓的数值模拟技术。 在实际工程领域中,用数值模拟技术可以对复杂的工程结构进行受力和响应分析,这样可以在设计或者加工前预知实体结构工作状态下的大概情况。 目前在工程实际应用中,常用的数值求解方法有:有限单元法、有限差分法、边界元等但从实用性和使用范围来说,有限单元法则是随着计算机技术的发展而被广泛应用的一种行之有效的数值计算方法。 2.2.2 有限元法 有限元法是一种基于能量原理的数值计算方法,是解决工程实际问题的一种有效的数值计算工具。它是里茨法的另一种表示形式,它可应用里茨法分析的所有弹性理论。 限元法是处理连续的结构体离散或有限个单元集合,也就是将连续的求解域离散为一定数量的单元集合体。且每个单元都具有一定的节点,相邻单元通过节点相互连续,同时使用等效节点力代替作用于单元上的力和选定场函数的节点值作为基本未知量。并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律:进而利用力学中的某些变分原理去建立用以求解节点未知量的有限元法方程,从而将一个连续域中的无限自由度问题化为离散域中的有限自由度问题。求解后,可利用解出的节点值和设定的插值函数确定整个单元集体上的场函数。 有限元求解问题中的单元分析:t t t a k F = 式中::t F 单元节点作用力。 t K :单元刚度矩阵。 t a :单元节点位移。 通过单元分析确定单元刚度矩阵,建立单元节点作用力和单元为伊关系。有限元求解问题时建立的结构整体平衡方程:P KU = 式中:P —结构整体等效点力载荷 K —结构总体刚度矩阵 U —结构节点位移阵列 单元内力的计算:t DBa =σ 式中:D —弹性矩阵 P —应变矩阵 整个结构的有限元分析就是一句上述方程而进行的具体的有限元求解过程如图

有限元分析技术的应用

计算机辅助分析 题目:有限元分析技术的应用 学院:机电工程学院 专业:机械设计制造及其自动化 班级: 姓名: 学号: 年月日

有限元分析技术的应用 摘要 有限元单元法,简称有限元法,是伴随着电子计算机技术的进步而发展起来 的一种新兴数值分析方法,是力学、应用数学与现代计算技术相结合的产物。有 限元法是一种高效能、常用的计算方法。本文主要讲述了有限元的特点、作用、 基本思想、分析步骤,以及有限元的应用,除此之外,也对有限元的应用软件进 和有限元的发展趋势行了简单介绍。 关键词:有限元法,基本思想,应用软件,发展趋势 The application of finite element analysis technology Summary The finite element method, finite element method, is accompanied by advances in computer technology and the development of a new numerical analysis method, is a product of mechanics, applied mathematics and modern technology combine. The finite element method is an efficient computing method, commonly used. This paper mainly describes the characteristics, finite element function, basic thought, analysis steps, and the application of finite element method, in addition, also do a simple introduction on the application software of finite element and finite element development trend. Keywords: finite element method, the basic idea, application, development trend

有限元的基础理论

§1有限元的基础理论 §1-1 概述 有限元法是一种数值计算的近似方法。早在40年代初期就已有人提出,但当时由于没有计算工具而搁置,一直到50年代中期,高速数字电子计算机的出现和发展为有限元法的应用提供了重要的物质条件,才使有限元法得以迅速发展。 有限元法在西方起源于飞机和导弹的结构设计,发表这方面文章最早而且最有影响的是西德的J.H.Argyris教授,于1954–1955年间,他在《Aircraft engineering》上发表了许多有关这方面的论文,并在此基础上写成了《能量原理与结构分析》,此书成为有限元法的理论基础。美国的M.T.Turner,R.W.Clough,H.C.Martin和L.J.Topp等人于1956年发表了一篇题为《复杂结构的刚度和挠度分析》一文,此文提出了计算复杂结构刚度影响系数的方法,说明了如何利用计算机进行分析。美国教授R.W.Clough于1960年在一篇介绍平面应力分析的论文中,首次提出了有限元法的名字。1965年英国的O.C.Zienliewice教授及其合作者解决了将有限元应用于所有场的问题,使有限元法的应用范围更加广泛。 有限元法的优点很多,其中最突出的优点是应用范围广。发展至今,不仅能解决静态的、平面的、最简单的杆系结构,而且还可以解决空间问题、板壳问题、结构的稳定性问题、动力学问题、弹塑性问题和粘弹性问题、疲劳和脆性断裂问题以及结构的优化设计问题。而且不论物体的结构形式和边界条件如何复杂,也不论材料的性质和外载荷的情况如何,原则上都能应用。 §1-2 有限元的基础理论 有限元法的基本思路和基本原则以结构力学中的位移法为基础,把复杂的结构或连续体看成有限个单元的组合,各单元彼此在节点处连接而组成整体。把连续体分成有限个单元和节点,称为离散化。先对单元进行特性分析,然后根据各节点处的平衡和协调条件建立方程,综合后作整体分析。这样一分一合,先离散再综合的过程,就是把复杂结构或连续体的计算问题转化为简单单元的分析与综合的问题。因此,一般的有限元解法包括三个主要步骤:离散化、单元分析、整体分析。 §1-2-1 离散化 一个复杂的弹性体可以看作由无限个质点组成的连续体。为了进行解算,可以将此弹性体简化为有限个单元组成的集合体,这些单元只在有限个节点上铰接,因此,这集合体只具有有限个自由度,这就为解算提供了可能。有无限个质点的连续体转化为有限个单元的集合体,就称为离散化。 §1-2-2 单元分析 单元分析首先要进行单元划分。在工程结构中,一般采用四种类型的基本单元,即标量单元、线单元(杆、梁单元)、面单元和体单元。四种基本单元的若干例子及各单元节点自由度(节点位移)表示在图(1-1)中。而单元划分一般注意下面几点: 一、从有限元本身来看,单元划分的越细,节点布置得越多,计算的结果越精确。但计算时间和计算费用的增加。所以在划分单元时对应兼顾这两个方面。 二、在边界比较曲折,应力比较集中,应力变化较大的地方,单元应划分的细点,而在应力变化平缓处单元划分的大些。单元由小到大应逐渐过渡。 三、对于三角形单元,三条边长应尽量接近,不应出现钝角,以免计算出现较大的偏差。对于矩形单元,长度和宽度也不应相差过大。 四、任意一个三角形单元的角点必须同时也是相邻单元边上的角点,而不能是相邻单元边上的内点。划分其他单元时也应遵循此原则。

相关文档
最新文档