有限元分析基础
有限元分析基础
第一讲第一章有限元的基本根念Basic Concepts of the Finite Element Method1.1引言(introduction)有限元(FEM或FEA)是一种获取近似边值问题的计算方法。
边值问题(boundary value problems, 场问题field problem )是一种数学问题(mathematical problems)(在所研究的区域,一些相关变量满足微分方程如物理方程、位移协调方程等且满足特定的区域边界)。
边值问题也称为场问题,场是指我们研究的区域,并代表一种物理模型。
场变量是满足微分方程的相关变量,边界条件代表场变量在场边界上特定的值(物理边界转化为数学边界)。
根据所分析物理问题的不同,场变量包括位移、温度、热量等。
1.2有限元法的基本思路 (how does the finite element methods work)有限元法的基本思路可以归结为:将连续系统分割成有限个分区或单元,对每个单元提出一个近似解,再将所有单元按标准方法组合成一个与原有系统近似的系统。
下面用在自重作用下的等截面直杆来说明有限元法的思路。
等截面直杆在自重作用下的材料力学解答图1.1 受自重作用的等截面直杆图1.2 离散后的直杆受自重作用的等截面直杆如图所示,杆的长度为L ,截面积为A ,弹性模量为E ,单位长度的重量为q ,杆的内力为N 。
试求:杆的位移分布,杆的应变和应力。
)()(x L q x N -=EAdx x L q EA dx x N x dL )()()(-== ⎰-==xx Lx EA q EA dx x N x u 02)2()()( (1))(x L EAq dx du x -==ε )(x L A q E x x -==εσ 等截面直杆在自重作用下的有限元法解答(1)离散化如图1.2所示,将直杆划分成n 个有限段,有限段之间通过一个铰接点连接。
称两段之间的连接点为结点,称每个有限段为单元。
有限元分析基础知识
2000,4
ANSYS单元分类
1. 杆单元,包括二维杆单元和三维杆单元,线性调节 元,主要包括: LINK1,LINK8,LINK10,LINK11,LINK180等。 2. 弹簧阻尼单元,包括COMBIN系列: COMBIN7,COMBIN14,COMBIN37,COMBIN40等。 3. 质量元,MASS21。
ANSYS/Structural求解功能
ANSYS/Structural求解功能
Static -- 结构静力问题(包括线性和非线性问题) Modal -- 模态振动特性计算分析(结构固有频率和振型) Harmonic -- 谐波分析 Transient -- 瞬态分析 Spectrum -- 谱分析 Eigen Buckling -- 特征值屈曲分析(线性) Substructural -- 子结构分析 。。。。。。
2000,4
有限元分析步骤(续)
• 集合所有单元的平衡方程,集合依据的是所有相邻 单元在公共节点 处的位移相等;建立总体的有限元方程组。 • 引入边界条件 • 求解有限元方程组,得到未知节点位移 • 计算单元应力,对不同的单元,对应力的处理还有不同的方法
2000,4
ANSYS文件结构
二进制文件 Jobname.db (数据库文件) Jobname.dbb (备份文件) Jobname.rst (结构分析结果文件) Jobname.rth (热分析结果文件) Jobname.rmg (电磁场分析结果文件) Jobname.rfl (流体分析结果文件) Jobname.tri (三角化刚度矩阵文件) Jobname.emat (单元矩阵文件) Jobname.esav (单元保存文件)
2000,4
简例(续)
有限元分析基础课件第一章
物体离散化 将某个工程结构离散为由各种单元组成的计算模型, 这一步称作单元剖分。 离散后单元于单元之间利用单元的节点相互连接起来; 单元节点的设置、性质、数目等应视问题的性质,描 述变形形态的需要和计算进度而定。 用有限元分析计算所获得的结果只是近似的。如果划 分单元数目非常多而又合理,则所获 得的结果就与实 际情况相符合。
1956年Turener和Clough等用有限元法第一次得 出了平面应力问题的正确答案。 1960年Clough又进一步应用有限元法处理了平面弹 性问题,并提出了有限元法的名称,这才使得有限元 法的理论和应用都得到了迅速发展。 20世纪70年代以后,随着计算机和软件技术的发展 有限元法得到了迅猛的发展。
对于实际的连续结构,任何位置的物体都是相 互连接、相互作用的,而在被离散成有限元模型 后,假设相邻单元除节点外都是不相互连接、不相 互作用的,这一点是不符合实际的,但当单元趋近 无限小、节点无限多时,则这种离散结构将趋近于 实际的连续结构。 有限元法的离散处理的本质就是将原始的无限 自由度的连续体物理系统转换成由有限个节点自由 度组成的离散系统,且当所分割的单元无限小时, 该离散系统完全等价于原始的连续系统。
有限元基础理论
与ANSYS应用
CAD/CAE/CAM:CAD 工具用于产品结构设计,形 成产品的数字化模型,有限元法则用于产品性能的分 析与仿真,帮助设计人员了解产品的物理性能和破坏 的可能原因,分析结构参数对产品性能的影响,对产 品性能进行全面预测和优化;帮助工艺人员对产品的 制造工艺及试验方案进行分析设计。当前,有限元法 在产品开发中的作用,已从传统的零部件分析、校核 设计模式发展为与计算机辅助设计、优化设计、数字 化制造融为一体的综合设计。
增强可视化的前置建模和后置数据处理功能 目前几乎所有的商业化有限元程序系统都有功能很强 的前置建模和后置数据处理模块。使用户能以可视图 形方式直观快速地进行网格自动划分,生成有限元分 析所需数据,并按要求将大量的计算结果整理成变形 图、等值分布云图,便于极值搜索和所需数据的列表 输出。
有限元分析理论基础
有限元理论基础有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
釆用不同的权函数和插值函数形式,便构成不同的有限元方法。
4.加权余量法:是指采用使余量的加权函数为零求得微分方程近似解的方法称为加权余量法。
(Weighted residual method WRM)是一种直接从所需求解的微分方程及边界条件出发,寻求边值问题近似解的数学方法。
加权余量法是求解微分方程近似解的一种有效的方法。
设问题的控制微分方程为:在V域内厶(")-八0 (5.1.1)在S 边界上〃(“)-& = 0 (5.1.2)式中:L、B——分别为微分方程和边界条件中的微分算子;f、g ——为与未知函数u无关的己知函数域值;u——为问题待求的未知函数当弄!J用力u权余•肚法求近丁以解首先在求耳军域上理立一个T式閑数H 一般兵升如下形式:仁土CN=NC(5.1.3)T M式中:c{----------- 彳寺定系数. 也可称为广义坐标;N:--- 取白完备函冬攵*S线.性无关的基函孕攵°由于〃一般只圮彳守求函缨攵U的近1以耳岂因u匕将式(5 1.3) 代入式(5 1 1)牙口式(5 1.2)后将诃•不誉斯兄,昔迅:| R] = L(flb— f在V域内\R B =B(^~g在S 边界上("14)城然 & 、尽反映了r式函竽攵与实解之问的偏差. 它丁门分另U称做内召卩牙口边界余覺。
若在域\'内引入内部权函数硏,在边界S上引入边界权函数W B 则可理立11个消除余甘的条件.一般可农示为:L兀W B1R B dS = 0 (/ = L2.L ,〃) (51-5)• V • S不同的权函数幵;和jr R反映了不同的消除余•眩的准则。
第二章有限元分析基础
自由度 位移 温度 电位 速度,压力 磁位
UX ROTZ UZ ROTX
结构 DOFs
机自学院安全断裂分析研究室
节点和单元
载荷
节点: 空间中的坐标位置,具有一 定自由度,存在相互物理作用。 单元: 一组节点自由度间相互作用 的数值、矩阵描述(称为刚度或系 数矩阵)。单元有线、面或实体以 及二维或三维的单元等种类。 有限元模型由一些简单形状的单元组成,单 元之间通过节点连接,并承受一定载荷。
在某一时刻发生虚位移 * ,虚位移产生虚应变 , 则外力F做的虚功
*
假设结构受到外力F的作用,内部产生应力
,
W
*
T
F
*
在单位体积上,结构的虚变形能为 结构的虚变形能为
T
,则整个
U
V
*
dV
T
根据虚位移原理,有
*
T
F
V
*
dV
1943年,Courant提出有限元法概念 1956年,Turner和Clough第一次用三角形单元离散飞机 机翼,借助有限元法概念研究机翼的强度及刚度 1960年,Clough正式提出有限元法(FEM)
20世纪60年代,我国数学家冯康把FEM总结成凡是椭圆 形偏微分方程都可用FEM求解
u0 (
机自学院安全断裂分析研究室
第二章
分析指导思想
有限元分析基础
化整为零,裁弯取直,以简驭繁,变难为易
历史典故
• 结构分析的有限元方法是由一批学术界和工业界的研究者 在二十世纪五十年代到二十世纪六十年代创立的。 • 有限元分析理论已有100多年的历史,是悬索桥和蒸汽锅 炉进行手算评核的基础。很多著名的大型有限元软件如 NASTRAN、ANSYS、ABAQUS 等。
有限元分析基础
第三章 杆系结构静力分析的有限单元法
3.1.2 坐标系
为了建立结构的平衡条件,对结构进行整体分 析,尚需要建立一个对每个单元都适用的统一坐标系,
即结构坐标系或称之为整体坐标系、总体坐标系。
图3-3 坐标系示意图
精品课件
29
第三章 杆系结构静力分析的有限单元法
3.1.3 向量表示
在有限单元法中力学向量的规定为:当线位移 及相应力与坐标轴方向一致时为正,反之为负;转角位 移和力矩,按右手法则定出的矢量方向若与坐标轴正向 相一致时为正。对于任意方向的力学向量,应分解为沿 坐标轴方向的分量。
精品课件
25
第三章 杆系结构静力分析的有限单元法
3.1 结构离散与向量表示
工程上许多由金属构件所组成的结构,如塔式 桁构支承架、起重机起重臂架、钢结构桥梁、钢结构建 筑等可以归结为杆系结构。杆系结构按各杆轴线及外力 作用线在空间的位置分为平面杆系和空间杆系结构。
杆系结构可以由杆单元、梁单元组成。
(a) Liebherr塔式起重机 起重机
为了保证解的收敛性,选用的位移函数应当满 足下列要求:
a. 单元位移函数的项数,至少应等于单元的 自由度数。它的阶数至少包含常数项和一次项。至于高 次项要选取多少项,则应视单元的类型而定。
精品课件
32
第三章 杆系结构静力分析的有限单元法
b. 单元的刚体位移状态和应变状态应当全部
包含在位移函数中。
c. 单元的位移函数应保证在单元内连续,以
(b) Liebherr履带式
(c) 钢结构桥梁
(d) 埃菲
尔铁塔
精图品3-课1 件杆系结构
26
第三章 杆系结构静力分析的有限单元法
3.1.1 结构离散化
02-01有限元分析基础-理论基础
Kq=f——————(1) 其中:K是整体刚度矩阵;
q是节点位移矩阵; f是载荷矩点位移 解有限元方程Kq=f可得到位移。在根据方
程组的特点来选择合适的计算方法。
通过上述分析了解到,有限元分析的基本 思路是“先离散在组装”,离散为了进行单 元分析,组装为了对整体结构进行分析。
σ=Eε—————(2-4) 将式(2-2)、式(2-3)代入到式(2-4) 后简化得到:
F=(AE/l)Δl—————(2-5) 式(2-5)与弹簧方程F=kx很相似。因此, 受轴向力作用的等截面杆看做一个弹簧,则:
keq=AE/l——————(2-6)
一、有限元分析理论基础
根据上述分析,杆件的截面面积都是在 一个方向上变化的。可以将杆件近似地看做 是由4个弹簧串联起来的模型。
(2)假定一个近似描述单元特性解 为研究典型单元的力学特性,不妨先考虑
横截面积为A、长度为l的杆件在外力F作用下 构件的变形。
杆件的平均应力由下式给出: σ=F/A————(2-2) 杆件的平均正应变ε为
ε=Δl/l————(2-3)
一、有限元分析理论基础
在弹性区域内,应力和应变服从胡克定 律,即:
1.2 定义单元特性 (2)定义单元的力学关系
根据单元的材料、形状、尺寸、节点数目、 位置等参数,找出单元节点力和节点位移的 关系式。 (3)计算等效节点力
物理模型离散化后,假定力是通过节点在 单元间进行传递的,但对于实际连续体,力 是通过单元的公共界面在单元间进行传递。
一、有限元分析理论基础
1.3 组装单元 利用结构中力的平衡条件和边界条件将各
利用以上模型,假定力施加在各节点上。 可根据有图中节点1~节点5的受力情况, 得到各节点上力的静平衡: 节点1:R1-k1(u2-u1)=0 节点2:k1(u2-u1)-k2(u3-u2)=0 节点3:k2(u3-u2)-k3(u4-u3)=0 节点2:k3(u4-u3)-k4(u5-u4)=0 节点2:k4(u5-u4)-P=0
有限元分析基础
有限元分析基础1.1 有限元法的优势在数值分析方法中有限元法是使用最广泛的,因为它对复杂的边界条件、各种几何形状、不同的材料性能有很强的适用性,对力学的各类问题、位势格列问题的计算格式有一定的通用性,这种分析方法还有良好的效率与计算精度。
就力学概念而言,有限单元法的基础还是传统力学分析方法,只不过结构变成了连续体,它是将一个结构力学中的一个结构整体分成了若干个基本结构构件,这些基本结构构件就被称作是“单元”。
有限单元法将一个结构整体看做由有限个相当微小的结构单元体。
进而假设各单元体仅在节点处产生力和位移,这样一个具有无限个自由度的连续体就简化变成了有限个自由度的力学模型,结构分析的方法就可以利用求解。
建立的模型每个节点的节点力和位移求解方式与结构力学方法是完全相同的,这也是精确地,显而易见的是物理模型将单元划分的越小,它将越接近于真实连续体,最终的解也将收敛于精确解。
从数学观点来看,可以利用有限单元法来求解偏微分方程问题的近似解。
很多位势问题和经典连续体的力学问题都是由未知场函数的偏微分方程组与一定的边界条件来表征的,有限元法就是通过变分原理以及分区插值等离散化处理,将这类二次泛函的极值问题形象转化为常见的一维多元线性代数方程,可以便于求解。
有限单元法分析的问题类型决定了采用何种有限元分析主题程序,当然既可以是静力学问题也可以是动力学问题,既可以可以是温度场或者流场问题,也可以是稳态场或者瞬态场问题,线性和非线性的问题等等。
1.2 有限单元法分析实现手段有限元分析软件是实现有限单元法的主要手段。
面向工程的有限元软件在国际上比较通用的有:ANSYS、ABAQUS、ADINA等。
ANSYS作为工程数值模拟软件的代表,是一款具有很多用途的有限元分析类型的软件,随着版本的不同,其分析功能也在完善和扩充当中,它可以灵活的提供结构线性分析和热扰动分析,也可以对一些结构、一定流体、电力问题、电磁场问题及碰撞等问题进行求解。
有限元分析的基本原理
有限元分析的基本原理有限元分析法是一种通用的数值分析技术,它利用有限数目的计算元素来对结构的应力、变形以及失效的可能性进行分析,它简化了复杂的工程结构在实际受力情况下的模拟计算,可以预测出构件的性能、变形和可能失效等。
有限元分析是用数学模型来模拟生活用来模拟工程中结构抗压、抗弯、抗剪、抗疲劳等性能。
有限元分析有三个基本原理:结构变形、力学方程和材料本构方程。
首先,有限元分析的基础原理是结构变形。
结构变形是指在施加外力作用下,受力的结构的空间变形和大小的变化,它是有限元分析的基础,该原理说明了满足力学方程的解决方法如何以有限元的形式出现。
通常情况下,我们会把构件的耦合变形分成很多小的计算元(这些计算元之间有连接约束),减少变形的不确定性,从而提高分析的准确性。
其次,有限元分析的基础原理是力学方程。
满足力学方程条件的解决方案就是有限元分析,也就是把问题分解成很多小的子问题来求解。
力学方程最常见的形式是基于有限元技术的动态和静态结构分析。
动态结构分析是指结构在某个加载下的振动反应,涉及到施加外力、弹性和惯性效应。
静态结构分析则指结构在不同类型外力作用下的变形。
最后,有限元分析的基础原理是材料本构方程。
材料本构方程是指材料受拉力作用而形成变形和应力的关系,它可以用来描述材料在承受外力时的作用。
本构方程有很多不同的形式,最常用的形式是弹性体的本构方程,它说明了当受到外力作用时,材料的拉伸和压缩的反应,从而将其应用于有限元分析技术。
以上就是有限元分析的基本原理,它是构成有限元分析的基础,而且这些基本原理也被广泛应用于工程中对结构性能进行模拟和分析。
有限元分析可以帮助工程师准确地估算出结构在特定加载条件下的变形和应力,也可以帮助他们判断结构在疲劳荷载作用下是否会发生破坏。
有限元分析也可以帮助设计者更好地分析结构在复杂(多变)条件下的性能,以确定结构的最优设计。
所以,有限元分析的基本原理是工程分析的基础,合理的运用可以节约大量的时间和精力,从而达到性能最优的结构设计。
第二章有限元分析基础
第二章有限元分析基础有限元分析是一种常用的工程计算方法,在工程学科中被广泛应用。
本章将介绍有限元分析的基本概念和基础知识。
有限元分析是一种数值分析方法,用于求解复杂的物理问题。
它的基本思想是将一个连续的物体或结构离散化为有限数量的基本单元,通过在每个单元上进行计算,最终得到整个物体或结构的行为。
这些基本单元通过节点连接在一起,形成了一个有限元网格。
通过在每个节点上求解方程,可以得到整个物体或结构的应力、变形等相关信息。
在有限元分析中,有三个重要的步骤:建模、离散和求解。
建模是指将实际物体或结构转化为数学模型的过程。
在建模过程中,需要确定物体或结构的几何形状、边界条件和力学性质等。
离散是指将物体或结构划分为有限数量的基本单元。
常用的基本单元有三角形、四边形和六面体等。
离散过程中需要确定每个基本单元的几何属性和材料性质等。
求解是指在离散的基础上,通过求解节点上的方程,得到物体或结构的应力、变形等结果。
求解过程中,需要确定节点的位移和应变等参数。
有限元分析的基本假设是在每个基本单元内,应力和应变满足线性关系。
这意味着在小变形和小位移的情况下,有限元分析是有效的。
此外,为了提高计算精度,通常会增加更多的基本单元。
但是,增加基本单元数量会增加计算复杂度和计算时间。
因此,在实际应用中,需要根据问题的复杂程度和计算资源的限制进行权衡。
有限元分析广泛应用于各个领域,例如结构力学、热传导、电磁场、流体力学等。
在结构力学中,有限元分析可以用于求解静力学和动力学问题。
在热传导中,有限元分析可以用于求解温度分布和热流问题。
在电磁场中,有限元分析可以用于求解电荷和电场分布等。
在流体力学中,有限元分析可以用于求解流速和压力分布等。
总之,有限元分析是一种重要的工程计算方法,可以用于求解各种物理问题。
通过建模、离散和求解等步骤,可以得到物体或结构的应力、变形等结果。
有限元分析在工程学科中有着广泛的应用前景,对于工程设计和优化起着重要作用。
学习有限元分析需要哪些有限元分析基础知识
学习有限元分析需要哪些有限元分析基础知识?有限元分析具有确保产品设计的安全合理性,同时采用优化设计,找出产品设计最佳方案,降低材料的消耗或成本; 在产品制造或工程施工前预先发现潜在的问题; 模拟各种试验方案,减少试验时间和经费等作用,越来越被应用,越来越的人不断开始学习有限元分析。
对于很多想开始学有限元分析的人都会有这么一个疑问,学习有限元分析需要哪些有限元分析基础知识呢?对于这个问题,看板网根据超过十年的企业和个人有限元分析培训经验,给各位想学习有限元分析的朋友们提点建议。
有限元分析(FEA,Finite Element Analysis)的基本概念是用较简单的问题代替复杂问题后再求解。
它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的(较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。
这个解不是准确解,而是近似解,因为实际问题被较简单的问题所代替。
由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。
有限元分析基础知识主要有,结构强度分析、振动频率分析、谐响应分析、扭曲分析、机构尺寸优化分析、疲劳分析、热力分析、跌落测试、响应谱分析等。
以下是一些建议:1,图书馆或书店都可以买到有限元教材,有的教材讲得深,有的教材讲得浅。
要是想在理论层面往深层次学习,还要学习一些数学基础,比如泛函分析、变分原理,但是,如果不专门研究一般用不了理解那么深刻。
2,要根据你从事的行业而定。
如果做力学有限元分析,起码要懂力学,就要学习力学理论知识,比如弹性力学等;做电磁有限元分析,起码要懂麦克斯韦方程组。
市场上卖的有限元教材一般都是结合力学讲的。
然后你可以学习有限元软件(比如ANSYS、ABAQUS等)解决具体的工程实际问题了。
如果对结构有限元分析感兴趣,应该从材料力学、弹性力学开始。
对应力、应变、平衡方程、本构关系、位移-应变关系等知识有了了解以后,可以学习变分法的知识,。
CAE课有限元分析理论基础
类型。
精度要求
03
根据问题对精度的要求,选择足够高阶的有限元以保证求解精
度。
常用有限元的介绍
四面体有限元
适用于解决三维问题,具有较高的计算效率 和适应性。
壳体有限元
适用于解决薄壁结构问题,能够模拟结构的 弯曲和变形。
六面体有限元
适用于解决二维和三维问题,精度较高但计 算效率较低。
梁有限元
适用于解决细长结构问题,能够模拟结构的 轴向拉伸和弯曲。
CAE课有限元分析理论基础
目 录
• 引言 • 有限元分析的基本原理 • 有限元的分类和选择 • 有限元分析的实现过程 • 有限元分析的应用实例 • 结论与展望
01 引言
目的和背景
目的
有限元分析(FEA)是一种数值分析方法,用于解决复杂的工程问题,如结构 分析、热传导、流体动力学等。本课程旨在使学生掌握有限元分析的基本原理 和应用。
弯曲有限元
适用于解决大变形问题,如结 构动力学、流体动力学等。
非线性有限元
适用于解决非线性问题,如塑 性力学、断裂力学等。
耦合有限元
适用于解决多物理场耦合问题 ,如流体-结构耦合、电磁-热
耦合等。
有限元的选择
问题特性
01
根据问题的物理特性、边界条件和求解精度要求选择合适的有
限元类型。
计算资源
02
考虑计算资源的限制,选择计算效率高、内存占用小的有限元
04 有限元分析的实现过程
建立模型
确定分析对象和边界条件
首先需要明确分析的对象和所受的边界条件, 这是建立有限元模型的基础。
几何建模
根据分析对象的特点,利用CAD软件建立几何 模型。
模型简化
有限元分析的数学基础
3.1 简单问题的解析求解
3.1.1 1D拉压杆问题 一个左端固定的拉杆在其右端承受一外力P,该 拉杆的长度为l,横截面积为A,弹性模量为E, 如图所示。
(1) 基本变量
由于该问题是为沿x方向的一维问题,因此 只有沿x方向的变量,而其它变量为零。即
(2) 基本方程 对原三维问题的所有基本方程进行简化, 只保留沿x方向的方程,有该问题的三大基 本方程和边界条件如下:
∂σ x = 0
①
∂x
εx
=
∂u ∂x
②
③
④ ⑤
(3) 求解 对方程①②③进行直接求解,可得到以下 结果
⑥
其中c和c1为待定常数,由边界条件BC④ 和⑤,可求出⑥中的常数c1=0, 因此,有最后的结果:
⑦
(4) 讨论1 若用经验方法求解(如材料力学的方法), 则需先作平面假设,即假设 为均匀分 布,则可得到
两端力(弯矩)
144
将弯矩以挠度的二阶导数来表示,即
(2) 求解
若用基于dxdy微体所建立的原始方程(即原平面
应力问题中的三大类方程)进行直接求解,比较
麻烦,并且很困难,若用基于以上简化的“特征
建模”方法所得到的基本方程进行直接求解则比
较简单,对本例问题(如为均匀分布),其方程
为:
145
这是一个常微分方程,其解的形式有
146
其中c0……c3为待定系数,可由四个边界条件 BC求出,最后有结果
(3) 讨论 该问题有关能量的物理量计算为:
应变能 147
外力功 势能
148
(1) 基本方程的建立 描述该变形体同样应有三大方程和两类边界 条件,有以下两种方法来建立基本方程。 (a)用弹性力学中dxdy微体建模方法推导三大
第2章有限元分析基础
第2章有限元法基础第1节有限单法的形成一、有限元法的形成在工程技术领域内,经常会遇到两类典型的问题。
其中的第一类问题,可以归结为有限个已知单元体的组合。
例如,材料力学中的连续梁、建筑结构框架和桁架结构。
我们把这类问题,称为离散系统。
尽管离散系统是可解的,但是求解这类复杂的离散系统,要依靠计算机技术;第二类问题,通常可以建立它们应遵循的基本方程,即微分方程和相应的边界条件。
例如弹性力学问题、热传导问题和电磁场问题等。
由于建立基本方程所研究的对象通常是无限小的单元,这类问题称为连续系统。
尽管已经建立了连续系统的基本方程,由于边界条件的限制,通常只能得到少数简单问题的精确解答。
对于许多实际的工程问题,还无法给出精确的解答。
为解决这个困难,工程师们和数学家们提出了许多近似方法。
在寻找连续系统求解方法的过程中,工程师和数学家从两个不同的路线得到了相同的结果,即有限元法。
有限元法的形成可以追溯到二十世纪50年代,来源于固体力学中矩阵结构法的发展和工程师对结构相似性的直觉判断。
从固体力学的角度来看,桁架结构等标准离散系统与人为地分割成有限个分区后的连续系统在结构上存在相似性。
1956年M.J.Turner,R.W.Clough,H.C.Martin,L.J.Topp在纽约举行的航空学会年会上介绍了一种新的计算方法,将矩阵位移法推广到求解平面应力问题。
他们把结构划分成一个个三角形和矩形的“单元”,利用单元中近似位移函数,求得单元节点力与节点位移关系的单元刚度矩阵。
1954—1955年,J.H.Argyris在航空工程杂志上发表了一组能量原理和结构分析论文。
1960年,Clough在他的名为“The finite element in plane stress analysis”的论文中首次提出了有限元(finite element)这一术语。
数学家们则发展了微分方程的近似解法,包括有限差分方法、变分原理和加权余量法。
有限元分析基础-文档资料
2019.8
内容结构
第一章 第二章 第三章 第四章 第五章 第六章 第七章 概述 结构几何构造分析 杆系结构静力分析的有限单元法 平面结构问题的有限单元法 等参元 空间问题的有限单元法 轴对称旋转单元
2
第一章 概述
1.1 有限单元法的概念 1.2 有限单元法基本步骤 1.3 工程实例
21
第二章 结构几何构造分析
② 反对称载荷作用
(a) 变形状态分析
(b) 铲运机工作装置插入工况有限元分析
图1-3 WJD-1.5型电动铲运机
8
第一章 概述
(a) KOMATSU液压挖掘机
(b) 某液压挖掘机动臂限元分析
图1-4 液压挖掘机
9
第一章 概述
图1-5 驾驶室受侧向力应力云图
图1-6 接触问题结构件应力云图
10
第一章 概述
图1-7 液压管路速度场分布云图
18
第二章 结构几何构造分析
对称结构在正对称载荷下,对称轴截面上只能产生 正对称的位移,反对称的位移为零;对称结构在反对称 载荷下,对称轴截面上只有反对称的位移,正对称的位 移为零。
(1) 具有奇数跨的刚架 ① 正对称载荷作用
2.2.3 结构对称性的利用
(a) 对称刚架
(b) 变形状态分析 图2-22对称性利用示意图
图1-8 磨片热应力云图
图1-9 支架自由振动云图
11
第二章 结构几何构造分析
2.1 结构几何构造的必要性 2.2 结构计算基本知识 2.3 结构几何构造分析的自由度与约束2.1 结构几何构造的必要性
结构是用来承受和传递载荷的。如果不计材料的 应变,在其受到任意载荷作用时其形状和位置没有发 生刚体位移时,称之为几何不变结构或几何稳定结构, 反之则称为几何可变结构或几何不稳定结构。几何可 变结构不能承受和传递载荷。对结构进行几何构造分 析也是能够对工程结构作有限单元法分析的必要条件。
有限元分析理论基础
有限元分析概念有限元法:把求解区域看作由许多小的在节点处相互连接的单元(子域)所构成,其模型给出基本方程的分片(子域)近似解,由于单元(子域)可以被分割成各种形状和大小不同的尺寸,所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件有限元模型:它是真实系统理想化的数学抽象。
由一些简单形状的单元组成,单元之间通过节点连接,并承受一定载荷。
有限元分析:是利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。
并利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。
线弹性有限元是以理想弹性体为研究对象的,所考虑的变形建立在小变形假设的基础上。
在这类问题中,材料的应力与应变呈线性关系,满足广义胡克定律;应力与应变也是线性关系,线弹性问题可归结为求解线性方程问题,所以只需要较少的计算时间。
如果采用高效的代数方程组求解方法,也有助于降低有限元分析的时间。
线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。
非线性问题与线弹性问题的区别:1)非线性问题的方程是非线性的,一般需要迭代求解;2)非线性问题不能采用叠加原理;3)非线性问题不总有一致解,有时甚至没有解。
有限元求解非线性问题可分为以下三类:1)材料非线性问题材料的应力和应变是非线性的,但应力与应变却很微小,此时应变与位移呈线性关系,这类问题属于材料的非线性问题。
由于从理论上还不能提供能普遍接受的本构关系,所以,一般材料的应力与应变之间的非线性关系要基于试验数据,有时非线性材料特性可用数学模型进行模拟,尽管这些模型总有他们的局限性。
在工程实际中较为重要的材料非线性问题有:非线性弹性(包括分段线弹性)、弹塑性、粘塑性及蠕变等。
2)几何非线性问题几何非线性问题是由于位移之间存在非线性关系引起的。
当物体的位移较大时,应变与位移的关系是非线性关系。
研究这类问题一般都是假定材料的应力和应变呈线性关系。
它包括大位移大应变及大位移小应变问题。
有限元分析基础
1.什么是等参数单元?(教材)坐标变换和单元内的场函数采用相同数目的节点参数及相同的插值函数,这种变换方法是等参数变换,这种变换方式能满足坐标变换的相容性,采用等参数变换的单元称之为等参数单元。
2.等参数单元的特点、基本条件、划分单元应注意的问题(教材习题)3.应用等参数单元时为什么要采用高斯积分,高斯积分点的数目如何确定?(教材习题)4.薄板弯曲问题的基本假设是什么?(其他参考书)(1)板弯曲钱垂直于中面的法线,在板弯曲后保持为直线,并垂直于弯曲后的中面。
(2)板面各水平层之间相互挤压(3)薄板受垂直于中面的载荷时可以为中间层各点设有平行于板面的位移.5.位移插值必须满足的三个条件:(教材)(1)位移插值函数应能满足单元的刚体位移(2)位移插值函数应能反映常量应变——常应变准则(3)位移插值函数应能保证单元内及相邻单元间位移的连续性——变形协调准则6.什么是轴对称问题?(其他参考书):轴对称物体的形变及应力分布不一定是轴对称的,只有当约束和载荷都对称于旋转轴时,轴对称物体的变形及应力分布才是轴对称的。
我们把满足上述条件的系统应力分析问题称为轴对称问题。
(教材):如果弹性体的几何形状、约束情况以及所受的外力,都是绕某一轴对称的,则弹性体的应力、应变和位移也就对称于这一轴,这种问题称为轴对称问题。
7.刚度矩阵性质(总刚):(1)对称性,关于正对角线对称(2)稀疏性,矩阵中有大量的零元素(3)带状分布,矩阵中非零元素在主对角线两侧呈带状分布10.形函数的性质。
(教材)(1)单元内任一点的三个形函数之和恒等于1,即Ni+Nj+Nm=1.(2)在节点i:Ni=1,Nj=0,Nm=0在节点j:Ni=0,Nj=1,Nm=0在节点m:Ni=0,Nj=0,Nm=111. 有限元法的特点(其他参考书)(1)概念清楚,容易理解(2)适应性强,应用范围广。
(3)有限元法采用矩阵形式表达,便于编制计算机程序,可以充分利用数字计算机的优势。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
(a) 结点载荷处理方式 (b) 等效结点载荷处理方式 图3-2杆系结构离散化示意图
28
第三章 杆系结构静力分析的有限单元法
3.1.2 坐标系
为了建立结构的平衡条件,对结构进行整体分析,
尚需要建立一个对每个单元都适用的统一坐标系,即结
构坐标系或称之为整体坐标系、总体坐标系。
图3-3 坐标系示意图
29
3.1 结构离散与向量表示
工程上许多由金属构件所组成的结构,如塔式桁构 支承架、起重机起重臂架、钢结构桥梁、钢结构建筑等 可以归结为杆系结构。杆系结构按各杆轴线及外力作用 线在空间的位置分为平面杆系和空间杆系结构。 杆系结构可以由杆单元、梁单元组成。
(a) Liebherr塔式起重机
(b) Liebherr履带式起重机
梁单元平面弯曲仅考虑结点的四个位移分 v j , j ,由材料力学知,各截面的转角: i , 量 i, x 故梁单元平面弯曲的位移表达式可分为仅包含四个 3, 2, 4 的多项式 v( x) 1 2 x 3 x 2 4 x 3 待定系数 1,
27
第三章 杆系结构静力分析的有限单元法
c. 变截面杆件可分段处理成多个单元,取各段中点 处的截面近似作为该单元的截面,各单元仍按等截面杆 进行计算。 d. 对曲杆组成的结构,可用多段折线代替,每端折 线为一个单元。如若提高计算精度,也可以在杆件中间 增加结点。 e. 在有限元法计算中,载荷作用到结点上。当结构 有非结点载荷作用时,应该按照静力等效的原则将其
有限元分析基础
内容结构
第一章 第二章 第三章 第四章 第五章 第六章 第七章 概述 结构几何构造分析 杆系结构静力分析的有限单元法 平面结构问题的有限单元法 等参元 空间问题的有限单元法 轴对称旋转单元
2
第一章 概述
1.1 有限单元法的概念 1.2 有限单元法基本步骤 1.3 工程实例
3
第一章 概述
(b) 铲运机工作装置插入工况有限元分析
图1-3 WJD-1.5型电动铲运机
8
第一章 概述
(a) KOMATSU液压挖掘机
(b) 某液压挖掘机动臂限元分析
图1-4 液压挖掘机
9
第一章 概述
图1-5 驾驶室受侧向力应力云图
图1-6 接触问题结构件应力云图
10
第一章 概述
图1-7 液压管路速度场分布云图
(c) 钢结构桥梁 图3-1 杆系结构
(d) 埃菲尔铁塔
26
第三章 杆系结构静力分析的有限单元法
3.1.1 结构离散化
由于杆系结构本身是由真实杆件联接而成,故离散 化比较简单,一般将杆件或者杆件的一段( 一根杆又分 为几个单元 )作为一个单元,杆件与杆件相连接的交点 称为结点。 杆系结构的离散化的要点可参考如下: a. 杆件的转折点、汇交点、自由端、集中载荷作用 点、支承点以及沿杆长截面突变处等均可设置成结点。 这些结点都是根据结构本身特点来确定的。 b. 结构中两个结点间的每一个等截面直杆可以设置 为一个单元。 变换为作用在结点上的等效结点载荷。
30
第三章 杆系结构静力分析的有限单元法
结点位移列向量为
i ui
e
vi i
T
u
j
j
vj j
T
单元e结点位移列向量为
i ui i i j
uj j j
T
结点力向量为
Fi
e
Ui V i
16
第二章 结构几何构造分析
(3) 按结构自由度分 ①静定结构——自由度为零的几何不变结构。其特征: a. 静定结构的内力及支座反力可全部由平衡方程式 求出,并且解答是唯一的。 b. 静定结构的内力及支座反力与材料的性质和截面 特征(几何尺寸,形状)无关。 c. 静定结构上无外载荷作用时,其内力及支座反力 全为零。 d. 若静定结构在载荷作用下, 结构中的某一部分 能不依靠于其它部分, 独立地与载荷保持平衡时,则 其它部分的内力为零。 e. 当将一平衡力系作用于静定结构的一个几何不 变部分时,结构的其余部分都无内力产生。 f. 当静定结构中的一个内部几何不变部分上的载 荷作等效变换时,其余部分的内力不变。 g. 当静定结构中的一个内部儿何不变部分作构造 改变时,其余部分的内力不变。
32
第三章 杆系结构静力分析的有限单元法
b. 单元的刚体位移状态和应变状态应当全部包含在 位移函数中。 c. 单元的位移函数应保证在单元内连续,以及相邻 单元之间的位移协调性。 由单元结点位移,确定待定系数项 当 x 0 时, u ui 当 x l 时, u u j 所以 u u 1 ui 2 j i l 用结点位移表示
结构计算所常用的结点和支座的简化形式: (1)结点: ① 铰结点;② 刚结点;③ 混合结点。 (2)支座: ① 活动铰支座;② 固定铰支座 ; ③ 固定支座 ;④ 定向支座
15
第二章 结构几何构造分析
2.2.2 结构的分类与基本特征
(1) 按结构在空间的位置分 结构可分为平面结构和空间结构两大类 (2) 按结构元件的几何特征分 ① 杆系结构: 梁、拱、桁架、刚架、桁构结构等 。 ② 板壳结构 ③ 实体结构实体结构的长、宽、高三个尺寸都很 大,具有同一量级。 ④ 混合结构
第三章 杆系结构静力分析的有限单元法
3.1.3 向量表示
在有限单元法中力学向量的规定为:当线位移及相
应力与坐标轴方向一致时为正,反之为负;转角位移和 力矩,按右手法则定出的矢量方向若与坐标轴正向相一 致时为正。对于任意方向的力学向量,应分解为沿坐标
轴方向的分量。
(a) 刚架结构示意图
(b) 结点位移和结点力分向量 图3-4 平面刚架分析示意图
Mi
eT
F U
e j
j
Vj
Mj
eT
单元e结点力列向量为
F
e
e Fi e U i Vi F j
Mi U j
Vj
Mj
eT
31
第三章 杆系结构静力分析的有限单元法
3.2 位移函数及单元的刚度矩阵
3.2.1 轴向拉压杆单元的位移的函数
(c) 对称性利用
19
第二章 结构几何构造分析
② 对称刚架承受反对称载荷作用
(a) 对称刚架
(b) 变形状态分析 图2-23 反对称性利用示意图
(c) 反对称性利用
ቤተ መጻሕፍቲ ባይዱ
20
第二章 结构几何构造分析
(2) 具有偶数跨的刚架 ① 正对称载荷作用
(a) 变形状态分析
(b) 对称性利用
图2-24对称性利用示意图
1.1 有限单元法的概念
基本思想:借助于数学和力学知识,利用计算机技术而
解决工程技术问题。
Finite Element Method -_FEM
Finite Element Analysis
4
第一章 概述
三大类型(按其推导方法分):
(1) 直接刚度法(简称直接法): 根据单元的物理意义,建立有关场变量表示的单元 性质方程。 (2) 变分法 直接从求解泛函的极值问题入手,把泛函的极植问 题规划成线性代数方程组,然后求其近似解的一种计算 方法。 (3) 加权余量法 直接从控制方程中得到有限单元方程,是一种近似 解法。
18
第二章 结构几何构造分析
对称结构在正对称载荷下,对称轴截面上只能产生 正对称的位移,反对称的位移为零;对称结构在反对称 载荷下,对称轴截面上只有反对称的位移,正对称的位 移为零。
(1) 具有奇数跨的刚架 ① 正对称载荷作用
2.2.3 结构对称性的利用
(a) 对称刚架
(b) 变形状态分析 图2-22对称性利用示意图
图1-8 磨片热应力云图
图1-9 支架自由振动云图
11
第二章 结构几何构造分析
2.1 结构几何构造的必要性 2.2 结构计算基本知识 2.3 结构几何构造分析的自由度与约束 2.4 自由度计算公式
12
第二章 结构几何构造分析
2.1 结构几何构造的必要性
结构是用来承受和传递载荷的。如果不计材料的 应变,在其受到任意载荷作用时其形状和位置没有发 生刚体位移时,称之为几何不变结构或几何稳定结构, 反之则称为几何可变结构或几何不稳定结构。几何可 变结构不能承受和传递载荷。对结构进行几何构造分 析也是能够对工程结构作有限单元法分析的必要条件。
23
第二章 结构几何构造分析
2.4 自由度计算公式
(1)桁架自由度计算公式 桁架中的结点数为j,杆件数为g,支座链杆数为z, 则桁架的自由度W 为 平面桁架 W 2jgz 空间桁架 W 3j g z (2) 平面混合结构的自由度计算公式 一个平面体系的自由度计算结果,不外下述三种 可能: a. W>0 表明结构缺少必要的约束, 可运动, 故 结构必定是几何可变体系。 b. W=0 表明结构具有保证几何不变所需的最少的 约束数。 c. W<0 表明结构具有多余约束。
5
第一章 概述
1.2 有限单元法基本步骤
(1) 待求解域离散化
(2) 选择插值函数
(3) 形成单元性质的矩阵方程 (4) 形成整体系统的矩阵方程 (5) 约束处理,求解系统方程 (6) 其它参数计算
6
第一章 概述
图1-2 工程问题有限单元法分析流程
7
第一章 概述
1.3 工程实例
(a) 铲运机举升工况测试
24
第三章 杆系结构静力分析的有限单元法
3.1 结构离散与向量表示
3.2 位移函数及单元的刚度矩阵
3.3 坐标变换及单元刚度矩阵 3.4 整体刚度矩阵 3.5 约束处理及求解 3.6 计算示例 3.7 ANSYS桁架结构计算示例 3.8ANSYS刚架结构计算示例