整体分析及总体刚度矩阵的性质

合集下载

最新7.4-单元刚度矩阵组装及整体分析

最新7.4-单元刚度矩阵组装及整体分析

根据全结构的平衡方程可知,总体刚度矩阵是由单元刚度矩阵集合而成的个结构的计算模型分成个单元,那么总体刚度矩阵可由各个单元的刚度矩阵组装而成,即是由每个单元的刚度矩阵的每个系数按其脚标编号“对号入座”叠加而成的将总体坐标轴分别用表示,对某单元有式中,和分别是局部坐标系和总体坐标系下的单元结点位移向量是该单元在总体坐标系下的单元刚度矩阵.将单元结点的局部编号换成总体编号,其中右上角的上标表示第单元所累加上的子矩阵具有相同的下标,的那些子矩阵的累加总体刚度矩阵第行的非零子矩阵是由与结点相联系的那,从环绕点各单元移置而来的结点载荷为式中表示对环绕结点的所有单元求和,环绕结点的各单元施加于结点的结点力为.因此,结点的平衡方程可表示为得到以结点位移表示的结点的平衡方程,为整体刚度矩阵,为全部结点位移组成的向量,为全部结点载荷组成的向量式中,是总体坐标系下的结点载荷向量,为坐标转换阵.构是处于自由状态,在结点载荷的作用下,结构可以产生任意的刚体位移的条件下,仍不能通过平衡方程惟一地解出结点位移.约束的种类包括使某些自由度上位移为零,,或给定其位移值,还有给定支承刚为了理解这个方法,我们把方程分块如下:其中,假设是给定的结点位移;是无约束的(自由)结点位移因而是已知的结点力;其中,不是奇异的,因而可以解方程(一旦知道了,求得未知结点力.殊情况下,我们可以删除对应于的各行和各列(即删行删列法),故可把方程简写为由于全部给定的结点位移通常都不能在位移向量的开始或终了,故分块法的编号方法是很麻烦因此,为了引入给定的边界条件,可以采用下述等价的方法如果把给定为,则载荷向量为结点自由度总数中对应于的行和列为零,而对角线元素为)在载荷向量中引入规定的值,即对全部规定的结点位移均应反复运用上述过程(步骤(置大数法的思路是:在总体刚度矩阵中,把指定位移所对应的行和列的对角元素乘上一个很大的数,如,此行其他元素保持不变,同时把该行对应的载荷项也相应地用来代替,这里为指定位移,于是原平衡方程组变为除第行外,其他各行仍保持原来的平衡特性,而第个方程式展开为由于上式中的比其他项的系数大得多,求和后可略去其小量,则上式变为即.边上有,若结构的总体坐标系为为斜支座的局部坐标系(见图对于边界结点,须限定方向位移,为此,将边界结点的位移及载荷都变换到局部坐标轴系设轴与斜支座的轴夹角为,逆时针为正,其中,.)中第行左右两边前乘以上式的系数矩阵仍然是对称的,而且此方程中结点位沿轴表示,这样,限定方向的位移异性,解这个线性代数方程组可求出结位移.阶线性代数方程,需进行次消元行元素作为主元行,为主元,对第行元素()的消元公式为式中等的上角码(次消元后的系数矩阵和载荷阵分别记为及.式表时第我们把消元最后结果记为,为上当回代求解时,已经解得总体刚度平衡方程中,,是单位上三角矩阵,.记,则.由其中第一个方程解得,再由第二个方程解得,向上回代,可得,由得依此类推可求得.由平衡方程组解出位移后,从中分离出各单元的结点位移,再通过方程)等计算各单元的应变、应力和结点力等内力。

土木工程-结构力学-重点分析

土木工程-结构力学-重点分析

学习目标1、理解矩阵位移法的内容2、掌握单元分析3、掌握整体分析4、掌握内力计算的原理5、掌握单元荷载处理6. 掌握桁架分析矩阵位移法矩阵位移法以传统的位移法为理论基础;以矩阵作为数学表达形式;以计算机作为计算工具三位一体解决各种杆系结构受力、变形等问题。

采用矩阵进行运算,公式紧凑,形式统一,便于使计算过程规格化和程序化。

适应计算机自动化计算的要求。

矩阵位移法结构力学传统方法与结构矩阵分析方法,二者同源而有别:在原理上同源,在作法上有别前者在“手算”的年代形成,后者则着眼于“电算”,计算手段的不同,引起计算方法的差异。

与传统的力法、位移法相对应,在结构矩阵分析中也有矩阵力法和矩阵位移法,或称柔度法与刚度法。

矩阵位移法由于具有易于实现计算过程程序化的优点而广为流传。

矩阵位移法1、矩阵位移法的基本思路a、方法的选择b、基本假设和基本原理线弹性、小变形。

满足叠加原理、功能原理c、正负号规定杆端内力、杆端位移、结点位移和结点力规定当与坐标轴正方向一致时为正;矩阵位移法1、矩阵位移法的基本思路原结构--离散--单元分析--整合2、离散(单元划分)为了减少基本未知量的数目,跨间集中荷载作用点可不作为结点,但要计算跨间荷载的等效结点荷载;跨间结点也可不作为结点,但要推导相应的单元刚度矩阵,编程序麻烦。

矩阵位移法 {}[]{}{}ee ef F k F δ=+单元分析的目的: 建立单元刚度方程单元分析的方法:利用形常数获得刚度系数,形成刚度矩阵; 利用载常数(固端力)叠加获得等效结点力。

单元分析如何操作:按自然位置选每跨为一个单元,支座处作为结点;分别给单元和结点编号;以结点位移作为基本未知量。

l li 2 i1 M 1 M2 M 3单元分析刚度矩阵的物理意义:•单元刚度矩阵是杆端力与杆端位移之物理关系的转换矩阵;•矩阵的阶数与杆端位移分量数相等;•系数kij 表示第j 个单位位移分量引起的第i 个杆端力分量数值的大小;•单元刚度矩阵具有对称性kij =kji 。

[工学]第4章 平面问题的有限元法-3刚度矩阵

[工学]第4章 平面问题的有限元法-3刚度矩阵
* 1 1 * 2 * 3 3
* T
F
T
* * * * * x x y * * y z z xy xy yz yz zx zx
({ } )
T
e T
R
e
(f)
而单元内的应力在虚应变上所做的功为
tdxdy
(g)
这里我们假定单元的厚度t为常量。把(d)式及(4-16) 式代入上式,并将提到积分号的前面,则有
({ } )
e T
B D B
T
e
tdxdy
根据虚位移原理,由(f)和(h)式可得到单元的虚功方程 即 e T e e T e T ({ } ) R ({ } ) B D B tdxdy 注意到虚位移是任意的,所以等式两边与相乘的项应该相等, 即得
R
e
B D Btdxdy
T
e

k B D B tdxdy
e T
(4-24) (4-25)
则有
R k
e e
e
上式就是表征单元的节点力和节点位移之间关系的刚 度方程,[k]e就是单元刚度矩阵。如果单元的材料是均质的 ,那么矩阵 [D] 中的元素就是常量,并且对于三角形常应 变单元,[B]矩阵中的元素也是常量。当单元的厚度也是常 量时,因 dxdy ,所以式(4-24)可简写为
1 2 4 7 11 3 5 8 6 9 10 15
12
13
14
图 4-6 a
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1 15
2
3
4
5

abaqus提取整体刚度矩阵

abaqus提取整体刚度矩阵

abaqus提取整体刚度矩阵abaqus是一款常用的有限元分析软件,主要用于模拟和分析各种工程结构的力学行为。

在abaqus中,整体刚度矩阵是一个重要的概念,它能够描述结构在受力作用下的刚度特性。

本文将深入探讨abaqus 如何提取整体刚度矩阵,并分享对该概念的观点和理解。

一、整体刚度矩阵的概念整体刚度矩阵是指在有限元分析中,将结构划分成若干个离散的单元后,通过单元刚度矩阵的叠加得到的描述结构整体刚度特性的矩阵。

整体刚度矩阵反映了结构在受力作用下的刚度响应,是进行结构力学分析的重要工具。

二、abaqus提取整体刚度矩阵的方法在abaqus中,提取整体刚度矩阵的方法主要有以下步骤:1. 创建有限元模型:需要在abaqus中创建一个准确表达所研究结构的有限元模型。

这包括定义结构的几何形状、材料性质以及边界条件等。

2. 定义材料属性:在有限元分析中,材料的力学性质对整体刚度矩阵具有重要影响。

在abaqus中需要明确定义结构中所使用的材料的力学性质,包括弹性模量、泊松比等。

3. 定义加载条件:接下来,需要定义结构在受力作用下的加载条件。

这可以是施加在结构上的力或约束条件等。

4. 进行力学分析:有了有限元模型、材料属性和加载条件后,就可以进行力学分析。

在abaqus中,通常使用有限元方法求解结构的响应,得到结构的位移和应力等。

5. 提取整体刚度矩阵:通过分析结果,abaqus提供了方便的工具来提取整体刚度矩阵。

用户可以在abaqus的后处理模块中选择相应的输出选项来得到整体刚度矩阵的结果。

三、对整体刚度矩阵的理解整体刚度矩阵是结构力学分析中的一个关键概念,对于研究和理解结构的强度和刚度特性具有重要意义。

整体刚度矩阵可以用来计算结构在受力作用下的位移、应力和应变等响应,进而评估结构的安全性和可靠性。

从数学角度看,整体刚度矩阵是由单元刚度矩阵叠加得到的。

单元刚度矩阵描述了单个有限元单元在特定边界条件下的刚度特性。

7.4 单元刚度矩阵组装及整体分析

7.4 单元刚度矩阵组装及整体分析

7.4 单元刚度矩阵组装及整体分析7.4.1 单刚组装形成总刚根据全结构的平衡方程可知,总体刚度矩阵是由单元刚度矩阵集合而成的.如果一个结构的计算模型分成个单元,那么总体刚度矩阵可由各个单元的刚度矩阵组装而成,即[K]是由每个单元的刚度矩阵的每个系数按其脚标编号“对号入座”叠加而成的.这种叠加要求在同一总体坐标系下进行.如果各单元的刚度矩阵是在单元局部坐标下建立的,就必须要把它们转换到统一的结构(总体)坐标系.将总体坐标轴分别用表示,对某单元有式中,和分别是局部坐标系和总体坐标系下的单元结点位移向量;[T]为坐标转换阵,仅与两个坐标系的夹角有关,这样就有是该单元在总体坐标系下的单元刚度矩阵.以后如不特别强调,总体坐标系下的各种物理参数均不加顶上的横杠.下面就通过简单的例子来说明如何形成总体刚度矩阵.设有一个简单的平面结构,选取6个结点,划分为4个单元.单元及结点编号如图3-27所示.每个结点有两个自由度.总体刚度矩阵的组装过程可分为下面几步:图7-27(1)按单元局部编号顺序形成单元刚度矩阵.图7-27中所示的单元③,结点的局部编号顺序为.形成的单元刚度矩阵以子矩阵的形式给出是(2)将单元结点的局部编号换成总体编号,相应的把单元刚度矩阵中的子矩阵的下标也换成总体编号.对下图3-27所示单元③的刚度矩阵转换成总体编号后为(3)将转换后的单元刚度矩阵的各子矩阵,投放到总体刚度矩阵的对应位置上.单元③的各子矩阵投放后情况如下:(4)将所有的单元都执行上述的1,2,3步,便可得到总体刚度矩阵,如式(3-9).其中右上角的上标表示第单元所累加上的子矩阵.(3-9)(5)从式(3-9)可看出,总体刚度矩阵中的子矩阵AB是单元刚度矩阵的子矩阵转换成总体编号后具有相同的下标,的那些子矩阵的累加.总体刚度矩阵第行的非零子矩阵是由与结点相联系的那些单元的子矩阵向这行投放所构成的.7.4.2 结点平衡方程我们首先用结构力学方法建立结点平衡方程.连续介质用有限元法离散以后,取出其中任意一个结点,从环绕点各单元移置而来的结点载荷为式中表示对环绕结点的所有单元求和,环绕结点的各单元施加于结点的结点力为.因此,结点的平衡方程可表示为(3-10)以[K]代入平衡方程,得到以结点位移表示的结点的平衡方程,对于每个结点,都可列出平衡方程,于是得到整个结构的平衡方程组如下:式中,[K]为整体刚度矩阵,为全部结点位移组成的向量,为全部结点载荷组成的向量.当然,如果各点的载荷向量也是在单元局部坐标下建立的,在合成以前,也应把它们转换到统一的结构(总体)坐标系下,即式中,是总体坐标系下的结点载荷向量,为坐标转换阵.7.4.3 位移边界条件在有限元法对结构进行整体分析时,建立了整体刚度矩阵[K],也得到了结构的刚度平衡方程,即.结构刚度方程的求解相当于总刚[K]求逆的过程.但是,从数学上看,未经处理的总刚是对称、半正定的奇异矩阵,它的行列式值为零,不能立即求逆.从物理意义看,在进行整体分析时,结构是处于自由状态,在结点载荷的作用下,结构可以产生任意的刚体位移.所以,在已知结点载荷的条件下,仍不能通过平衡方程惟一地解出结点位移.为了使问题可解,必须对结构加以足够的位移约束,也就是应用位移边界条件.首先要通过施加适当的约束,消除结构的钢体位移,再根据问题要求设定其他已知位移.所以,处理位移边界条件在有限元分析步骤中十分重要.约束的种类包括使某些自由度上位移为零,,或给定其位移值,还有给定支承刚度等,本书涉及前两种.处理约束的方法,常用的有删行删列法、分块法、置大数法和置“1”法等,下面分别予以介绍.1、删行删列法若结构的某些结点位移值为零时(即与刚性支座连接点的位移),则可将总体刚度矩阵中相应的行列、删行删列划掉,然后将矩阵压缩即可求解.这种方法的优点是道理简单.如果删去的行列很多,则总体刚度矩阵的阶数可大大缩小.通常用人工计算时常采用该方法.若用计算机算题,在程序编制上必带来麻烦,因为刚度矩阵压缩以后,刚度矩阵中各元素的下标必全改变.因而一般计算机算题不太采用.2.分块法为了理解这个方法,我们把方程分块如下:(3-11)其中,假设是给定的结点位移;是无约束的(自由)结点位移.因而是已知的结点力;是未知的结点力.方程(3-11)可以写为即(3-12)和(3-13)其中,不是奇异的,因而可以解方程(3-12)得出(3-14)一旦知道了,就可以由方程(3-13)求得未知结点力.在全部给定的结点自由度都等于零的特殊情况下,我们可以删除对应于的各行和各列(即删行删列法),故可把方程简写为(3-15)3.置“1”法由于全部给定的结点位移通常都不能在位移向量的开始或终了,故分块法的编号方法是很麻烦的.因此,为了引入给定的边界条件,可以采用下述等价的方法.可以把方程(3-12)和(3-13)合在一起写为(3-16)在实际计算中,方程(3-16)所示的过程可以在不重新排列所述方程的情况下用下述分块的方法为进行.步骤(1)如果把给定为,则载荷向量P可以修改为为结点自由度总数.步骤(2)除对角线元素以外,使[K]中对应于的行和列为零,而对角线元素为1,即步骤(3)在载荷向量中引入规定的值,即对全部规定的结点位移均应反复运用上述过程(步骤(1)到(3)).应当指出,由于这个过程保持了方程的对称性,因此,[K]可以按带状存储,而且几乎不会增加编制程序的工作量.4.置大数法置大数法的思路是:在总体刚度矩阵中,把指定位移所对应的行和列的对角元素乘上一个很大的数,如,此行其他元素保持不变,同时把该行对应的载荷项也相应地用来代替,这里为指定位移,于是原平衡方程组变为除第行外,其他各行仍保持原来的平衡特性,而第个方程式展开为由于上式中的比其他项的系数大得多,求和后可略去其小量,则上式变为即.这样就用近似方程组代替原方程组,得到近似满足边界条件的解.当指定位移为零时,只要将对角元素乘上一个大数,而相应的载荷项经证明可以不置零.删行删列法适用于指定零位移点,而置大数法适用于给定位移(包括零位移).5.斜支座的处理对于简单的约束情况(如限定某些结点位移为零或取得给定数值),可以用前述置大数法处理.有的结构在直角坐标系内建立了位移方程组,但在某个斜边上受有法向约束.如图3-28所示正方形固支板,受均布横向载荷,对此,可利用对称性而只计算其1/8,如图中ABC部分,其中AC为固支边,按对称性,AB边上有,但在BC边上应限定绕BC的转用等于零.为处理此类斜边上的约束,须对斜边上的结点做坐标变换.若结构的总体坐标系为为斜支座的局部坐标系(见图3-29).对于边界结点,须限定方向位移,为此,将边界结点的位移及载荷都变换到局部坐标轴系.设轴与斜支座的轴夹角为,逆时针为正,图7-28 图7-29 则依据第二单中坐标转换关系有其中,.或写成(3-17)与位移关系相同有(3-18)将上两式带入结构刚度方程有(3-19)这样把位移到列阵中凡是斜支座的结点位移矢量都用局部坐标表示了.将式(3-19)中第行左右两边前乘以(3-20)由上式可见:凡是边界点的斜支座,在刚度方程中对应于斜支座的位移和载荷向量均可直接斜支座的局部坐标值,总刚度距阵中的相应行列需作相应的变换.上式的系数矩阵仍然是对称的,而且此方程中结点位沿轴表示,这样,限定方向的位移就很方便了.实际计算中,并不需要建立结构总的位移方程组后再进坐标变换.而可以在形成单元刚度矩阵和结点载荷之后,就对斜支座点进行坐标变换,把变换后的单元刚度矩阵和结点载荷叠加入总刚度矩阵和总载荷的相应位置,最后叠加形成的也就是方程组(3-20),即需要处理的结点,应该在单元计算中完成坐标变换后再叠加,当结构有不同的斜边约束时,都可以这样处理,只不过对不同边上的结点,应按不同的方向余弦矩阵变换就是了.7.4.4 总刚度平衡方程的求解应用有限元法,最终都是归结为解总体刚度平衡方程,它实际上是以总体刚度矩阵为系数矩阵的大型线性代数方程组.通过对结构施加位移边界条件,消除了结构的刚体位移,从而消除总体刚度矩阵的奇异性,解这个线性代数方程组可求出结位移.我们已知,总体刚度矩阵具有大型、对称、稀疏、带状分布、正定、主元占优势的特点,稀疏表示将对称消元法进一步改造,使之适合总刚的等带宽二维存储.(4)因子化法(三角分解)又称Cholesky分解,适合一维变带宽存储总刚.这上方法储效率高,计算速度快,应用较为普遍.此外,还有一种方法,叫做波前法.波前法实际上也是一种改进的高斯消去法.它建立一个称为“波前”的空间,各单元刚度系数依次进入波前.一旦与某自由度有关的所有单元的刚度系数全部装入,便可将相应的变量消去.经过消元的方程的系数随即退出波前,存放在计算机的外存中.这样就可腾出空间装入新的刚度系数.所以,波前法不需要生成完整的总刚,而是边组装边消元,“成熟”一个消去一个.消元完成后,全部系数都已存储在计算机的外存或缓冲区中.回代时将各方程的系数按“先出后入”的顺序调入内存求解.由此可见,这种方法是利用计算机充裕的外存资源,以多耗取机时来缓解内存不足的矛盾,以便适应较大规模的问题.随着计算机技术的发展,内存资源不断扩大,对具有稀疏、带状性质的有限元刚度方程,这种以时间换取空间的办法得不偿失.另一方面,波前法的阐述和程序设计比较复杂,且对多种单元并存的结构使用不便.所以,本书不拟介绍波前法.本书第九章将详细讨论适合整体存储总刚的高斯消去法和适合一维变带宽存储的因子化法以及有关的程序设计问题,以下仅列出这两种方法的梗概.1、高斯消去法高斯循序消去法的一般公式:对于n阶线性代数方程,需进行次消元.采用循序消去时,第m次消元以m-1次消元后的m行元素作为主元行,为主元,对第行元素()的消元公式为(3-21)式中等的上角码(m),表示该元素是经过第m次消元后得到的结果.同样,可以把经过m次消元后的系数矩阵和载荷阵分别记为及.式表时第m 次消元是在经m-1次消元的基础上进行的.消元过程中,主元及被消元素的位置可见图3-30(a).图中阴影部分已完成消元过程的元素,主元行以下的矩阵为待消部分.在进行第m次时,1-m行元素的消元过程已经完成,其中的元素就是消元最后得到的上三角阵中的元素. m行发下的元素消元过程尚未结束,连同m行元素在内构成一个待消的方阵.消元共需进行n-1次.消元完成后,即可回代求解.我们把消元最后结果记为,为上三角阵,回代公式可写作(3-22)回代过程自后向前进行.当回代求解时,已经解得.回代示意图见图3-30(b),阴影部分为已求得解答的部分.图7-30 高斯消去法2.三角分解法总体刚度平衡方程中,[K]是对称、正定矩阵,因而可做如下分解(3-23)其中,则是单位上三角矩阵,.代入整本结构平衡方程记,则.即由向下回代.由其中第一个方程解得,再由第二个方程解得,……,依此类推可求得{Y}.又由向上回代,可得,由得依此类推可求得.由上述过程可见,三角分解法求解线性代数方程组的关键是对系数矩阵进行三角分解.7.4.5 求解内力由平衡方程组解出位移后,从中分离出各单元的结点位移,再通过方程(3-3)、(3-4)和(3-6)等计算各单元的应变、应力和结点力等内力。

有限元总体刚度矩阵的特点

有限元总体刚度矩阵的特点

有限元总体刚度矩阵的特点
以下是 8 条关于有限元总体刚度矩阵特点的内容:
1. 有限元总体刚度矩阵那可是具有对称性啊!就好像你照镜子,左边和右边是一样的呢!比如说在分析一个结构体的时候,它的内力和变形的关系也是对称的呀。

2. 嘿,它还有稀疏性哦!这就像一片大森林里,树与树之间有很多空隙呢!比如说对于一个大规模的结构,很多元素之间其实并没有直接的联系,这就导致矩阵中有很多零元素呢。

3. 有限元总体刚度矩阵还具有带状分布特点呢!就好比一条有规律的彩带一样。

像是在分析那种规则形状的结构,你就能明显看到这种带状的特征呀。

4. 哇塞,它的奇异性也不能忽视呀!这就好像一个独特的存在。

举个例子,当结构处于某些特殊状态时,它的表现就会特别明显呢。

5. 有限元总体刚度矩阵是可分块的哟!就好像把一个大拼图分成小块一样。

比如在处理复杂结构的时候,我们就可以把它分块来更好地理解和计算呀。

6. 它还具有正定性能哦!这就如同一个总是积极向上的小伙伴。

像是一个稳定的结构,它的总体刚度矩阵就会表现出这种正定性呢。

7. 有限元总体刚度矩阵对结构的描述超级精确啊!简直就像给结构拍了一张超级清晰的照片。

比如在设计一个重要的工程结构时,它的作用可太大了。

8. 有限元总体刚度矩阵的这些特点真的很重要好不好!这可是我们理解和分析结构的关键呀!没有它,很多复杂的问题都没法解决呢!
我的观点结论:有限元总体刚度矩阵的这些特点让我们能够更深入、准确地对各种结构进行分析和理解,在工程和科学研究中有着不可或缺的地位。

7.4单元刚度矩阵组装及整体分析报告材料

7.4单元刚度矩阵组装及整体分析报告材料

7.4 单元刚度矩阵组装及整体分析7.4.1 单刚组装形成总刚根据全结构的平衡方程可知,总体刚度矩阵是由单元刚度矩阵集合而成的.如果一个结构的计算模型分成个单元,那么总体刚度矩阵可由各个单元的刚度矩阵组装而成,即[K]是由每个单元的刚度矩阵的每个系数按其脚标编号“对号入座”叠加而成的.这种叠加要求在同一总体坐标系下进行.如果各单元的刚度矩阵是在单元局部坐标下建立的,就必须要把它们转换到统一的结构(总体)坐标系.将总体坐标轴分别用表示,对某单元有式中,和分别是局部坐标系和总体坐标系下的单元结点位移向量;[T]为坐标转换阵,仅与两个坐标系的夹角有关,这样就有是该单元在总体坐标系下的单元刚度矩阵.以后如不特别强调,总体坐标系下的各种物理参数均不加顶上的横杠.下面就通过简单的例子来说明如何形成总体刚度矩阵.设有一个简单的平面结构,选取6个结点,划分为4个单元.单元及结点编号如图3-27所示.每个结点有两个自由度.总体刚度矩阵的组装过程可分为下面几步:图7-27(1)按单元局部编号顺序形成单元刚度矩阵.图7-27中所示的单元③,结点的局部编号顺序为.形成的单元刚度矩阵以子矩阵的形式给出是(2)将单元结点的局部编号换成总体编号,相应的把单元刚度矩阵中的子矩阵的下标也换成总体编号.对下图3-27所示单元③的刚度矩阵转换成总体编号后为(3)将转换后的单元刚度矩阵的各子矩阵,投放到总体刚度矩阵的对应位置上.单元③的各子矩阵投放后情况如下:(4)将所有的单元都执行上述的1,2,3步,便可得到总体刚度矩阵,如式(3-9).其中右上角的上标表示第单元所累加上的子矩阵.(3-9)(5)从式(3-9)可看出,总体刚度矩阵中的子矩阵AB是单元刚度矩阵的子矩阵转换成总体编号后具有相同的下标,的那些子矩阵的累加.总体刚度矩阵第行的非零子矩阵是由与结点相联系的那些单元的子矩阵向这行投放所构成的.7.4.2 结点平衡方程我们首先用结构力学方法建立结点平衡方程.连续介质用有限元法离散以后,取出其中任意一个结点,从环绕点各单元移置而来的结点载荷为式中表示对环绕结点的所有单元求和,环绕结点的各单元施加于结点的结点力为.因此,结点的平衡方程可表示为(3-10)以[K]代入平衡方程,得到以结点位移表示的结点的平衡方程,对于每个结点,都可列出平衡方程,于是得到整个结构的平衡方程组如下:式中,[K]为整体刚度矩阵,为全部结点位移组成的向量,为全部结点载荷组成的向量.当然,如果各点的载荷向量也是在单元局部坐标下建立的,在合成以前,也应把它们转换到统一的结构(总体)坐标系下,即式中,是总体坐标系下的结点载荷向量,为坐标转换阵.7.4.3 位移边界条件在有限元法对结构进行整体分析时,建立了整体刚度矩阵[K],也得到了结构的刚度平衡方程,即.结构刚度方程的求解相当于总刚[K]求逆的过程.但是,从数学上看,未经处理的总刚是对称、半正定的奇异矩阵,它的行列式值为零,不能立即求逆.从物理意义看,在进行整体分析时,结构是处于自由状态,在结点载荷的作用下,结构可以产生任意的刚体位移.所以,在已知结点载荷的条件下,仍不能通过平衡方程惟一地解出结点位移.为了使问题可解,必须对结构加以足够的位移约束,也就是应用位移边界条件.首先要通过施加适当的约束,消除结构的钢体位移,再根据问题要求设定其他已知位移.所以,处理位移边界条件在有限元分析步骤中十分重要.约束的种类包括使某些自由度上位移为零,,或给定其位移值,还有给定支承刚度等,本书涉及前两种.处理约束的方法,常用的有删行删列法、分块法、置大数法和置“1”法等,下面分别予以介绍.1、删行删列法若结构的某些结点位移值为零时(即与刚性支座连接点的位移),则可将总体刚度矩阵中相应的行列、删行删列划掉,然后将矩阵压缩即可求解.这种方法的优点是道理简单.如果删去的行列很多,则总体刚度矩阵的阶数可大大缩小.通常用人工计算时常采用该方法.若用计算机算题,在程序编制上必带来麻烦,因为刚度矩阵压缩以后,刚度矩阵中各元素的下标必全改变.因而一般计算机算题不太采用.2.分块法为了理解这个方法,我们把方程分块如下:(3-11)其中,假设是给定的结点位移;是无约束的(自由)结点位移.因而是已知的结点力;是未知的结点力.方程(3-11)可以写为即(3-12)和(3-13)其中,不是奇异的,因而可以解方程(3-12)得出(3-14)一旦知道了,就可以由方程(3-13)求得未知结点力.在全部给定的结点自由度都等于零的特殊情况下,我们可以删除对应于的各行和各列(即删行删列法),故可把方程简写为(3-15)3.置“1”法由于全部给定的结点位移通常都不能在位移向量的开始或终了,故分块法的编号方法是很麻烦的.因此,为了引入给定的边界条件,可以采用下述等价的方法.可以把方程(3-12)和(3-13)合在一起写为(3-16)在实际计算中,方程(3-16)所示的过程可以在不重新排列所述方程的情况下用下述分块的方法为进行.步骤(1)如果把给定为,则载荷向量P可以修改为为结点自由度总数.步骤(2)除对角线元素以外,使[K]中对应于的行和列为零,而对角线元素为1,即步骤(3)在载荷向量中引入规定的值,即对全部规定的结点位移均应反复运用上述过程(步骤(1)到(3)).应当指出,由于这个过程保持了方程的对称性,因此,[K]可以按带状存储,而且几乎不会增加编制程序的工作量.4.置大数法置大数法的思路是:在总体刚度矩阵中,把指定位移所对应的行和列的对角元素乘上一个很大的数,如,此行其他元素保持不变,同时把该行对应的载荷项也相应地用来代替,这里为指定位移,于是原平衡方程组变为除第行外,其他各行仍保持原来的平衡特性,而第个方程式展开为由于上式中的比其他项的系数大得多,求和后可略去其小量,则上式变为即.这样就用近似方程组代替原方程组,得到近似满足边界条件的解.当指定位移为零时,只要将对角元素乘上一个大数,而相应的载荷项经证明可以不置零.删行删列法适用于指定零位移点,而置大数法适用于给定位移(包括零位移).5.斜支座的处理对于简单的约束情况(如限定某些结点位移为零或取得给定数值),可以用前述置大数法处理.有的结构在直角坐标系内建立了位移方程组,但在某个斜边上受有法向约束.如图3-28所示正方形固支板,受均布横向载荷,对此,可利用对称性而只计算其1/8,如图中ABC部分,其中AC为固支边,按对称性,AB边上有,但在BC边上应限定绕BC的转用等于零.为处理此类斜边上的约束,须对斜边上的结点做坐标变换.若结构的总体坐标系为为斜支座的局部坐标系(见图3-29).对于边界结点,须限定方向位移,为此,将边界结点的位移及载荷都变换到局部坐标轴系.设轴与斜支座的轴夹角为,逆时针为正,图7-28 图7-29 则依据第二单中坐标转换关系有其中,.或写成(3-17)与位移关系相同有(3-18)将上两式带入结构刚度方程有(3-19)这样把位移到列阵中凡是斜支座的结点位移矢量都用局部坐标表示了.将式(3-19)中第行左右两边前乘以(3-20)由上式可见:凡是边界点的斜支座,在刚度方程中对应于斜支座的位移和载荷向量均可直接斜支座的局部坐标值,总刚度距阵中的相应行列需作相应的变换.上式的系数矩阵仍然是对称的,而且此方程中结点位沿轴表示,这样,限定方向的位移就很方便了.实际计算中,并不需要建立结构总的位移方程组后再进坐标变换.而可以在形成单元刚度矩阵和结点载荷之后,就对斜支座点进行坐标变换,把变换后的单元刚度矩阵和结点载荷叠加入总刚度矩阵和总载荷的相应位置,最后叠加形成的也就是方程组(3-20),即需要处理的结点,应该在单元计算中完成坐标变换后再叠加,当结构有不同的斜边约束时,都可以这样处理,只不过对不同边上的结点,应按不同的方向余弦矩阵变换就是了.7.4.4 总刚度平衡方程的求解应用有限元法,最终都是归结为解总体刚度平衡方程,它实际上是以总体刚度矩阵为系数矩阵的大型线性代数方程组.通过对结构施加位移边界条件,消除了结构的刚体位移,从而消除总体刚度矩阵的奇异性,解这个线性代数方程组可求出结位移.我们已知,总体刚度矩阵具有大型、对称、稀疏、带状分布、正定、主元占优势的特点,稀疏表示将对称消元法进一步改造,使之适合总刚的等带宽二维存储.(4)因子化法(三角分解)又称Cholesky分解,适合一维变带宽存储总刚.这上方法储效率高,计算速度快,应用较为普遍.此外,还有一种方法,叫做波前法.波前法实际上也是一种改进的高斯消去法.它建立一个称为“波前”的空间,各单元刚度系数依次进入波前.一旦与某自由度有关的所有单元的刚度系数全部装入,便可将相应的变量消去.经过消元的方程的系数随即退出波前,存放在计算机的外存中.这样就可腾出空间装入新的刚度系数.所以,波前法不需要生成完整的总刚,而是边组装边消元,“成熟”一个消去一个.消元完成后,全部系数都已存储在计算机的外存或缓冲区中.回代时将各方程的系数按“先出后入”的顺序调入内存求解.由此可见,这种方法是利用计算机充裕的外存资源,以多耗取机时来缓解内存不足的矛盾,以便适应较大规模的问题.随着计算机技术的发展,内存资源不断扩大,对具有稀疏、带状性质的有限元刚度方程,这种以时间换取空间的办法得不偿失.另一方面,波前法的阐述和程序设计比较复杂,且对多种单元并存的结构使用不便.所以,本书不拟介绍波前法.本书第九章将详细讨论适合整体存储总刚的高斯消去法和适合一维变带宽存储的因子化法以及有关的程序设计问题,以下仅列出这两种方法的梗概.1、高斯消去法高斯循序消去法的一般公式:对于n阶线性代数方程,需进行次消元.采用循序消去时,第m次消元以m-1次消元后的m行元素作为主元行,为主元,对第行元素()的消元公式为(3-21)式中等的上角码(m),表示该元素是经过第m次消元后得到的结果.同样,可以把经过m次消元后的系数矩阵和载荷阵分别记为及.式表时第m 次消元是在经m-1次消元的基础上进行的.消元过程中,主元及被消元素的位置可见图3-30(a).图中阴影部分已完成消元过程的元素,主元行以下的矩阵为待消部分.在进行第m次时,1-m行元素的消元过程已经完成,其中的元素就是消元最后得到的上三角阵中的元素. m行发下的元素消元过程尚未结束,连同m行元素在内构成一个待消的方阵.消元共需进行n-1次.消元完成后,即可回代求解.我们把消元最后结果记为,为上三角阵,回代公式可写作(3-22)回代过程自后向前进行.当回代求解时,已经解得.回代示意图见图3-30(b),阴影部分为已求得解答的部分.图7-30 高斯消去法2.三角分解法总体刚度平衡方程中,[K]是对称、正定矩阵,因而可做如下分解(3-23)其中,则是单位上三角矩阵,.代入整本结构平衡方程记,则.即由向下回代.由其中第一个方程解得,再由第二个方程解得,……,依此类推可求得{Y}.又由向上回代,可得,由得依此类推可求得.由上述过程可见,三角分解法求解线性代数方程组的关键是对系数矩阵进行三角分解.7.4.5 求解内力由平衡方程组解出位移后,从中分离出各单元的结点位移,再通过方程(3-3)、(3-4)和(3-6)等计算各单元的应变、应力和结点力等内力。

有限单元法部分课后题答案

有限单元法部分课后题答案

1.1 有限单元法中“离散”的含义是什么?有限单元法是如何将具有无限自由度的连续介质问题转变成有限自由度问题的?位移有限元法的标准化程式是怎样的?(1)离散的含义即将结构离散化,即用假想的线或面将连续体分割成数目有限的单元,并在其上设定有限个节点;用这些单元组成的单元集合体代替原来的连续体,而场函数的节点值将成为问题的基本未知量。

(2)给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律,即通过插值以单元节点位移表示单元内任意点的位移。

因节点位移个数是有限的,故无限自由度问题被转变成了有限自由度问题。

(3)有限元法的标准化程式:结构或区域离散,单元分析,整体分析,数值求解。

1.3 单元刚度矩阵和整体刚度矩阵各有哪些性质?各自的物理意义是什么?两者有何区别?单元刚度矩阵的性质:对称性、奇异性(单元刚度矩阵的行列式为零)。

整体刚度矩阵的性质:对称性、奇异性、稀疏性。

单元Kij物理意义Kij 即单元节点位移向量中第j个自由度发生单位位移而其他位移分量为零时,在第j个自由度方向引起的节点力。

整体刚度矩阵K 中每一列元素的物理意义是:要迫使结构的某节点位移自由度发生单位位移,而其他节点位移都保持为零的变形状态,在所有个节点上需要施加的节点荷载。

2.2什么叫应变能?什么叫外力势能?试叙述势能变分原理和最小势能原理,并回答下述问题:势能变分原理代表什么控制方程和边界条件?其中附加了哪些条件?(1)在外力作用下,物体内部将产生应力σ和应变ε,外力所做的功将以变形能的形式储存起来,这种能量称为应变能。

(2)外力势能就是外力功的负值。

(3)势能变分原理可叙述如下:在所有满足边界条件的协调位移中,那些满足静力平衡条件的位移使物体势能泛函取驻值,即势能的变分为零δ∏p=δ Uε+δV=0此即变分方程。

对于线性弹性体,势能取最小值,即δ2∏P=δ2Uε+δ2V≥0此时的势能变分原理就是著名的最小势能原理。

有限元法基础重点归纳(精)

有限元法基础重点归纳(精)
29、常应变三角形单元:当单元确定后。矩阵B是常量,单元中任一点的应变分量也是常量的单元。
30、有限元法的任务:建立和求解整个弹性体的节点位移和节点力之间的关系的平衡方程。31、单元刚度矩阵:表达了单元节点位移与节点力之间的转换关系。
32、单元刚度矩阵的性质:①单元刚度矩阵中每个元素有明确的物理意义②K e是对称矩阵③K e的每一行或每一列元素之和为零,因此K e为奇异矩阵④K e不随单元的平行移动或作n π角度的转动而改变。33、刚度集成法集成规律:①先对每个单元求出其单元刚度矩阵K e ,而且以分块形式按节点编号顺序排列②将单元刚度矩阵扩大阶数为2n*2n ,并将单元刚度矩阵中的子块按局部码与总码的对应关系,搬到扩大后的矩阵中,形成单元贡献矩阵K e。③将所有单元贡献矩阵同一位置上的分块矩阵简单叠加成总体刚度矩阵中的一个子矩阵,各行各列都按以上步骤即形成总体刚度矩阵K。34、整体刚度矩阵的性质:①整体刚度矩阵是对称矩阵②整体刚度矩阵中每一元素的物理意义:整体刚度矩阵的第一列元素代表使第一个节点在x方向有一单元位移,而其余节点位移皆为零时必须在节点上施加的里。对于K的其余各列也有类似意义③整体刚度矩阵K的主对角线上的元素总是正的④整体刚度矩阵K是一个稀疏阵⑤整体刚度矩阵K是一个奇异阵。35、带形矩阵:整体刚度矩阵K的非零元素分布在以主对角线为中心的斜带形区域内的矩阵。
γxy
=E 1−μ
2∗
1−μ2
γxy
42、制造位移函数:{u (x,y =α1+α2x +α3y
v (x,y =α4+α5x +α6y
43、等参单元精度比四边形单元高,四边形精度比三角形精度高。
44、轴对称问题:很多工程物件,它们的几何形状承受的载荷以及约束条件都对称于其一固定轴,这即为对称轴,此时载荷作用下的位移、应变和应力也对称于该对称轴的问题。45、等参数单元:优点:①形状方位任意,适应性好,精度高,容易构造高阶单元②具有统一形式,规律性强,采用数值积分算,程序处理方便③高阶等参单元精度高,描述复杂边界,形状能力强,所需单元少。缺点:①单元各方向尺寸要尽量接近②单元边界不能过于曲折,不能有拐点折点,尽量接近直线或抛物线③边之间夹角要尽量接近直角④单元形状不能过度畸变,边中节点不能过于偏离中间。46、有限元法基础理论:弹性力学,材料力学

弹性力学整体刚度矩阵的特点与存储方法

弹性力学整体刚度矩阵的特点与存储方法

元计算秉承中国科学院数学与系统科学研究院有限元自动生成核心技术(曾获中科院科技进 步二等奖、国家科技进步二等奖),通过自身不懈的努力与完善,形成一系列具有高度前瞻性和 创造性的产品。
元计算产品适用范围广泛,目前有国内外专业客户300余家,涉及美、加、日、韩、澳、德、 新等国,遍布石油化工、土木建筑、电磁电子、国防军工、装备制造、航空航天……等多个领域。
弹性力学整体刚度矩阵的特点 与存储方法
用有限元方法分析复杂工程问题时,结点的数目比较多,整体刚度矩阵的阶数通常也是很高的。那么, 是否在进行计算时要保存整体刚度矩阵的全部元素?能否根据整体刚度矩阵的特点提高计算效率。 整体刚度矩阵具有以下几个显著的特点:对称性,稀疏性,非零系数带形分布。 1)对称性 由单元刚度矩阵的对称性和整体刚度矩阵的集成规则,可知整体刚度矩阵必为对称矩阵。利用对称性, 只保存整体矩阵上三角部分的系数即可。 2)稀疏性 单元刚度矩阵的多数元素为零,非零元素的个数只占较小的部分。如图3-26 所示的结构,结点2 只和 通过单元联接的1、3、4、5 结点相关,结点5 只和通过单元联接的2、3、4、6、8、9 结点相关。由单元 刚度矩阵的物理意义和整体刚度矩阵的形成方式可知,相关结点2、3、4、6、8、9 及结点5 本身产生位移 时,才使结点5 产生结点力,其余结点产生位移时不在该结点处引起结点力。在用分块形式表示的整体矩 阵中,与相关结点对应的分块 矩阵具有非零的元素,其它位置上的分块矩阵的元素为零,如图3-27 所示。
二维等带宽存储 设整体刚度矩阵[K]为一个n n的矩阵,最大半带宽为d。利用带形矩阵的特点和对称性, 只需要保存以d 为固定带宽的上半带的元素,称为二维等带宽存储。进行存储时,把整体刚度 矩阵[K]每行中的上半带元素取出,保存在另一个矩阵[K*]的对应行中,得到一个n x d 矩阵[K*]。 把元素在[K]矩阵中的行、列编码记为r、s,在矩阵[K*]中的行、列编码记为r*、s*,对应关 系如下: r*=r,s*=s-r+1

刚度矩阵的性质和存储

刚度矩阵的性质和存储

! ! !
k ji ! kn1i
! ! !
k jj ! kn1 j
! ! !
k jn1 !
⎥ ⎥ ⎥
kn1n1 ⎥⎦
过虚功原理得到证明。
6
3、结构刚度矩阵主对角线上的元素恒为正值
由性质(1)可知,任一主对角线上元素kii是使 节点位移Δi为一单位位移,其它节点位移为零时必须 在第i号位移方向施加的力Pi。它的方向自然应与位移 方向相同,因而是正值。
⎢ ⎢
"
"
"
"
"
"
"
"
"
⎥ ⎥
⎢0 ! kii ! kij ! kim ! 0⎥ i
⎢ ⎢
"
"
"
"
"
"
"
"
"
⎥ ⎥
[k ]2n×2n = ⎢0 ! k ji ! k jj ! k jm ! 0⎥ j
⎢ ⎢
"
"
"
"
"
"
"
"
"
⎥ ⎥
⎢0
⎢ ⎢
"
! !
kmi "
! "
kmj "
! "
kmm "
! "
0⎥
"
个节点的节点数大得多,因而,结构刚度矩阵中很大
t
一部分元素是零,即所谓的稀疏矩阵。
9
5、结构刚度矩阵是一个奇异矩阵
从单元刚度矩阵的奇异性讨论中知,处于静力 平衡状态的无约束单元可以发生任意的刚体位移。 与单元刚度矩阵是奇异矩阵的理由一样,无约束结 构的结构刚度矩阵[K]也是奇异矩阵,即[K]的行列 式为零。

刚架的整体刚度矩阵[详细]

刚架的整体刚度矩阵[详细]
第9章 矩阵位移法
§9-1 概述 §9-2 单元刚度矩阵(局部坐标系) §9-3 单元刚度矩阵(整体坐标系) §9-4 连续梁的整体刚度矩阵 §9-5 刚架的整体刚度矩阵 §9-6 结构整体结点荷载 §9-7 计算步骤和算例
▲ 竖向杆件坐标变换的简化技巧 §9-8 忽略轴向变形时刚架的整体分析 §9-9 桁架及组合结构的整体分析
0 30 100 0 30 50 3

104 ×300 0
0 300 0
0
0
3 0 12 30 0 12 30 0
0
30 50
0
30
100
4
1
3
单元② 900
k② T T k ② T
0 1 0
1 0 0
0
T
0
01
0
1 0
0
1 0 0
0 0 1
1 2 300 0
0
12 30
0
12
30
2
104
×
0 300
30 0
100 0
0 300
30 0
50 3
0
0
0 12 30 0 12 30 0
0
30 50
0
30
100
4
1230 00
12 0 30 12 0 30 1
0
300
0
0
300
0
2
104
×
30 12
0 0
100 30 30 12
0 0
解:1)编号、建立坐标如图所示。 2)写出各单元局部坐标下的 刚度矩阵
1(1,2,3) ①

2 y (0,0,0)

整体及总体刚度矩阵的性质概述

整体及总体刚度矩阵的性质概述

整体及总体刚度矩阵的性质概述整体及总体刚度矩阵是一个方阵,其尺寸等于结构体系中自由度的个数。

整体刚度矩阵可以表示为K = [k11, k12, ..., k1n; k21,k22, ..., k2n; ..., kn1, kn2, ..., knn],其中ki,j表示结构体系中第i个自由度在受到第j个自由度作用时的刚度系数。

总体刚度矩阵是整体刚度矩阵的一种特殊形式,在结构的分析与计算中较为常见。

总体刚度矩阵一般通过将各个单元的刚度矩阵逐个组合得到。

总体刚度矩阵包含了结构的所有自由度,反应了整个结构在受力作用下的刚度特性。

1. 对称性:整体及总体刚度矩阵是对称矩阵,即kij = kji。

这是由于结构在平衡状态下受力成立的一个基本条件。

对称性使得计算和分析过程更加简化,可以减少计算量。

2.正定性:整体及总体刚度矩阵是正定矩阵,即对于任意非零的向量v,v^TKv>0。

正定性保证了整体及总体刚度矩阵的特征值均为正数,即不存在零特征值。

这意味着结构不会出现无穷大的位移和变形,具有稳定性和可靠性。

3.奇异性:整体及总体刚度矩阵是奇异矩阵的条件是存在零特征值。

如果结构体系有刚度为零的单元或自由度,则整体及总体刚度矩阵的秩将小于其自由度的个数,从而成为奇异矩阵。

奇异性代表结构的不稳定性,需要进行特殊处理或修正。

4.加法性:整体及总体刚度矩阵具有加法性,即当结构被分解成若干个结构单元(子结构)时,每个子结构的刚度矩阵加和得到整个结构体系的刚度矩阵。

这使得结构计算和分析可以被分解和简化,提高了效率。

5.可逆性:整体及总体刚度矩阵是可逆矩阵,即存在逆矩阵K^(-1),使得K·K^(-1)=K^(-1)·K=I。

逆矩阵的存在保证了结构计算的唯一性,可以通过刚度矩阵求解结构的位移和反力。

6.非线性性:整体及总体刚度矩阵的计算涉及到结构的几何非线性和材料非线性。

当结构存在较大的变形和应力非线性时,刚度矩阵的计算需要进行迭代,并考虑材料的非线性特性。

有限元 整体分析

有限元 整体分析
K r s,K rs,K r s,K rs 表示结点S(S=i,j,m)在水平方向、垂直方向产
生单位位移时,在结点r(r=i,j,m)上分别所要施加的水平结点
力和垂直结点力的大小。例如 K i j 表示结点j在垂直方向产生 单位位移时,在结点i所需要施加的水平结点力的大小。
单元刚度矩阵的性质:
1)对称性: K 是对称矩阵
元素K的脚码,标有“-”的表示水平方向,没有标“-” 的表示垂直方向。
单元刚度矩阵的物理意义:
Ur

S i, m j,
(K r s us K rsvs )(r i,j, m)
Vr

S i, m j,
(K r s us K rsvs )(r i,j, m)
单元刚度矩阵的每一个元素都有明显的物理意义。
3-3 单元刚度矩阵
根据虚功原理,有
F

e
[B] [D][B]tdxdy δ
T
e
K
e

[B] [D][B]tdxdy
T
F K δ 建立了单元的结点力与结点位移之间的关系, Ke 称为单 元刚度矩阵。它是6*6矩阵,其元素表示该单元的各结点沿坐标 方向发生单位位移时引起的结点力,它决定于该单元的形状、 大小、方位和弹性常数,而与单元的位置无关,即不随单元或 坐标轴的平行移动而改变。
2)奇异性: Ke 是奇异矩阵, K
e
e
0
单元刚度矩阵所有奇数行的对应元素之和为零,所有偶 数行的对应元素之和也为零。由此可见,单元刚度矩阵各列元 素的总和为零。由对称性可知,各行元素的总和也为零。
单元刚度矩阵的性质: 例题:求下图所示单元的刚度矩阵,设 1、求[B]

有限元复习题

有限元复习题

有限元法基本原理复习资料1、线性弹性力学中一般哪些基本假设 ?什么是理想弹性体?2、线弹性材料物体内任意一点,一定存在三个相互垂直的主应力σ1、σ2、σ3,假设材料的柏松比为μ,弹性模量为E,则三个应变ε1、ε2、ε3可以表达为:3、弹性力学基本方程的导出,可从三方面分析:通过平衡微分方程建立了应力、体力和面力之间的关系。

通过几何方程建立了应变、位移和边界位移之间的关系。

通过物理方程建立了应变与应力之间的关系。

4、写出并理解弹性力学的基本方程。

a.平衡微分方程:b.几何方程:1. 平面问题中的几何方程:2. 空间问题的几何方程:c、物理方程:或者:为体积应变即:简写成:{σ}=[D]{ε} 式中[D]称为弹性矩阵,它完全由弹性常数E 和μ 决定。

4、请表述如图所示边界条件:5、如图所示的线弹性材料可以归结为:平面应力问题,其不为零的应力分量有:6、如图所示的线弹性材料可以归结为:平面应变问题,其不为零的应变分量有:εx ,εy,γxy7、描述并理解平面问题的基本方程平面应力问题和平面应变问题都只有8 个独立的未知量,它们只是x 和y 的函数,因此统称平面问题。

1. 平面问题的平衡微分方程2. 平面问题中的几何方程:3. a.平面应力问题中的物理方程:记作{σ}=[D]{ε} 其中[D]为弹性矩阵。

b.平面应变问题中的物理方程:记作{σ}=[D]{ε} 其中[D]为弹性矩阵。

比较两种平面问题的弹性矩阵,可以发现,将平面应力问题物理方程中的弹性常数E 、μ 换成就可得到平面应变问题物理方程。

8、结构的分类与基本特征(1) 按结构在空间的位置分结构可分为平面结构和空间结构两大类(2) 按结构元件的几何特征分① 杆系结构:梁、拱、桁架、刚架、桁构结构等 。

② 板壳结构③ 实体结构实体结构的长、宽、高三个尺寸都很大,具有同一量级。

④ 混合结构 (3) 按结构自由度分① 静定结构——自由度为零的几何不变结构。

弹性力学整体刚度矩阵的特点与存储方法

弹性力学整体刚度矩阵的特点与存储方法

性的过程。
02
整体刚度矩阵可以用于描述结构的动态特性,如固有
频率、模态振型等。
03
通过动态分析,可以预测结构的振动、冲击和稳定性
等行为,为结构的优化设计和安全评估提供依据。
06 结论与展望
研究结论
整体刚度矩阵具有对称性
整体刚度矩阵是弹性力学分析中的重要矩阵,它具有对称性,即矩阵的主对角线元素和副对角线元素相等。这种对称 性反映了弹性力学中应力和应变之间的对称关系。
整体刚度矩阵的存储方法
由于整体刚度矩阵具有对称性,可以采用压缩存储的方法来减小存储空间。具体而言,可以只存储矩阵的上三角部分 或下三角部分,然后通过对称性来获取另一部分的元素。这样可以大大减少存储空间,提高计算效率。
整体刚度矩阵的稳定性
整体刚度矩阵的稳定性对于弹性力学分析的精度和稳定性至关重要。在实际应用中,需要采用稳定的方 法来计算和存储整体刚度矩阵,以避免误差的积累和失真。
弹性力学整体刚度矩阵的特点与存 储方法
contents
目录
• 引言 • 弹性力学基础 • 整体刚度矩阵的特点 • 整体刚度矩阵的存储方法 • 整体刚度矩阵的应用 • 结论与展望
01 引言
背景介绍
弹性力学是研究物体在受力作用下的 变形和内力的学科,整体刚度矩阵是 弹性力学中描述物体刚度的核心工具 。
在结构优化设计中的应用
01
结构优化设计是指在满足一定约束条件下,寻找使 某一目标函数达到最优的结构设计方案。
02
整体刚度矩阵可以用于描述结构的刚度特性,进而 用于结构优化设计的计算和分析。
03
通过调整结构的设计变量,可以改变整体刚度矩阵, 从而得到最优的结构设计方案。
在动态分析中的应用

有限元方法基础教程第三版答案第二单元

有限元方法基础教程第三版答案第二单元

有限元方法基础教程第三版答案第二单元1. 诉述有限元法的定义答:有限元法是近似求解一般连续场问题的数值方法2. 有限元法的基本思想是什么答:首先,将表示结构的连续离散为若干个子域,单元之间通过其边界上的节点连接成组合体。

其次,用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。

3. 有限元法的分类和基本步骤有哪些答:分类:位移法、力法、混合法;步骤:结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移。

4. 有限元法有哪些优缺点答:优点:有限元法可以模拟各种几何形状复杂的结构,得出其近似解;通过计算机程序,可以广泛地应用于各种场合;可以从其他CAD软件中导入建好的模型;数学处理比较方便,对复杂形状的结构也能适用;有限元法和优化设计方法相结合,以便发挥各自的优点。

缺点:有限元计算,尤其是复杂问题的分析计算,所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。

对无限求解域问题没有较好的处理办法。

尽管现有的有限元软件多数使用了网络自适应技术,但在具体应用时,采用什么类型的单元、多大的网络密度等都要完全依赖适用者的经验。

5. 梁单元和平面钢架结构单元的自由度由什么确定答:由每个节点位移分量的总和确定6. 简述单元刚度矩阵的性质和矩阵元素的物理意义答:单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵单元刚度矩阵中元素aml的物理意义为单元第L个节点位移分量等于1,其他节点位移分量等于0时,对应的第m个节点力分量。

7. 有限元法基本方程中的每一项的意义是什么P14 答:Q——整个结构的节点载荷列阵(外载荷、约束力);整个结构的节点位移列阵;结构的整体刚度矩阵,又称总刚度矩阵。

8. 位移边界条件和载荷边界条件的意义是什么答:由于刚度矩阵的线性相关性不能得到解,引入边界条件,使整体刚度矩阵求的唯一解。

9. 简述整体刚度矩阵的性质和特点P14 答:对称性;奇异性;稀疏性;对角线上的元素恒为正。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1、对称性。 只存贮矩阵的上三角部分,节省近一半的存贮容量。
2、稀疏性。 矩阵的绝大多数元素都是零,非零元素只占一小部分。
整体刚度矩阵的特点
2、稀疏性。 矩阵的绝大多数元素都是零,非零元素只占一小部分。
节点5只与周围的六个节点 (2、3、4、6、8、9)用三角形 单元相连,它们是5的相关节 点。只有当这七个相关节点产 生位移时,才使该节点产生节 点力,其余节点发生位移时并 不在该节点处引起节点力。因 此,在矩阵[K]中,第5行的非 零子块只有七个(即与相关节 点对应的七个子块)。
整体刚度矩阵的特点
2、稀疏性。
1 2 3 4 5 6 7 8 9 10
一般,一个节点的相关结
1
点不会超过九个,如果网格中
2
有200个节点,则一行中非零
3
子块的个数与该行的子块总数
4
相比不大于9/200,即在5%以
5
下,如果网格的节点个数越多,
6
则刚度矩阵的稀疏性就越突出。
7
利用矩阵[K]的稀疏性,
8
可设法只存贮非零元素,从而
9
可大量地节省存贮容量。
10
整体刚度矩阵的特点
3、带形分布规律。
上图中,矩阵[K]的非零元素分布在以对角线为中心 的带形区域内,称为带形矩阵。在半个带形区域中(包括对角 线元素在内),每行具有的元素个数叫做半带宽,用d表示。 半带宽的一般计算公式是:
半带宽 d = ( 相邻结点码的最大差值 + 1 ) * 2 上图中相邻节点码的最大差值为4,故d=(4+1)*2=10
m2 4
i2 j3 5 m4
i4 6
[K jj](1)
[K jm ](1)
[K ji ](1)
[K mm ](1) [K jj](2) [K ii ](3)
[K mi ](1) [K im ](3)
[K jm ](2)
[K ii ](1)
[K mm](3) [K jj](4)
[K ji ](2) [K ij ](3)
利用带形矩阵的特点并利用对称性,可只存贮上半带的 元素,叫半带存贮。
整体刚度矩阵的特点
图(a)中的矩阵[K]为n行n列矩阵,半带宽为d。半带存贮
时从[K]中取出上半带元素,按图(b)中的矩阵 [K]* 的排列方
式进行存贮,即将上半部斜带换成竖带。存贮量n*d,存贮量 与[K]中元素总数之比为d/n,d值越小,则存贮量约省。
整体刚度矩阵的形式
2、刚度矩阵的集成规则:
1)在整体离散结构变形后,应保 证各单元在节点处仍然协调地相互 连接,即在该节点处所有单元在该 节点上有相同位移,
2
1
i
2
i
L
in i
1
4



i

3
2)整体离散结构各节点应满足平 衡条件。即环绕每个节点的所有单 元作用其上的节点力之和应等于作 用于该节点上的节点载荷Ri,
总码 1
2
34
56
方程为(或节点力表
1
达式):
j2 2 3
m2 4 5
i2 6
[K jj](2)
[K jm](2) [K ji](2)
[K mj](2) [Kij ](2)
[K mm](2)[K mi](2) [Kim ](2) [Kii ](2)
0
F22
0
F42
1
K%
(
2)
Fie Ri
2
e
1
4

Ri
④ i


3
整体刚度矩阵的形式
• 2、整体刚度矩阵的集成方法
• 具阵体集k e成,方然法后是将:其先中对的每每个个单子元块求出k单ij 元送刚到度结矩构刚 度矩阵中的对应位置上去,进行迭加之后即得出结构 刚度矩阵[K]的子块,从而得出结构刚度矩阵[K]。
• 关键是如何找出 k e 中的子块在[K]中的对应位置。
F 2
K
21
K 22
K 23
K 24
K 25
K
26
2
F F
34
K31 K 41
K32 K 42
K33 K 43
K34 K 44
K35 K 45
K36 K 46
3 4
FF56
K51
K61
K52 K62
K53 K63
K54 K64
K55 K65
K56 K66
5 6
其中子向量 i 和 Fi 都是二阶向量,子矩阵 Kij
大刚度矩阵 K%e 。 B、将各单元的扩大刚度矩阵
K%e 迭加,得出结构刚度
矩阵[K]。
2)R R1 L RnT 为节点载荷向量, 1 L nT
为节点位移向量。
局部码
j1
m1, j2 ,i3 i1, m3 , j4
m2
i 2 , j3 , m4
i4
总码
1
2
3
4
5
6
j1 1
m1 j2 2 i3 i1 m3 3 j4
[K mj](3) [K jm ](4)
[K ji ](4)
[K mm](2) [K mi](2)
[K ii ](2) [K jj](3) [K mm](4)
[K mi](4)
[K ii ](4)
2. 整体刚度矩阵的特点
在有限元法中,整体刚度矩阵的阶数通常是很高的,在 解算时常遇到矩阵阶数高和存贮容量有限的矛盾。找到整体 刚度矩阵的特性达到节省存贮容量的途径。
2 34
F52
0
56
整体刚度矩阵的形式

用同样的方法可得出其他单
元的扩大的单元刚度方程:
F1e
F2e1 2 Fra bibliotekF3e F4e
K%(e)
34
e=1,2,...4
F5e
F6e
56
• 据节点力平衡,各个单元相 应节点力叠加:
Fie Ri i=1,...6
u1
0,u4
0,v4
0,v6
0
建立节点平衡方程时,应根据上述边界条件进行处理。
3、解方程组,求出节点位移。
通常采用消元法和迭代法两种方法。
4、根据节点位移求出应力。
整体刚度矩阵的形式 整体刚度矩阵 K 是单元刚度矩阵 k e 的集成。
1、刚度集成法的物理概念:
在单元刚刚度阵矩阵中中的k元,e素,即表ki由示ej 节j节点点作单单位位位位移移,时其引他起节的点节位点移力为。
零时,单元e在i节点引起的节点力;类似,在整体刚阵中, 表示j节kij 点 单位位移,其他节点位移为零时,整体结构在i节点 引起的节点力(由于结构已被离散为一系列单元,即所有与i、 j节点相关的单元在i节点引起的节点力之和)。
如上图结构,计算 k23 时,与节点2和3相关的单元有单元 ①和③,当节点3发生单位位移时,相关单元①和③同时在节点 2引起节点力,将相关单元在节点2的节点力相加,就得出结构 在节点2的节点力 k23 k213 k233 。由此看出,结构的刚度 系数是相关单元的刚度系数的集成,结构刚度矩阵中的子块是 相关单元的对应子块的集成。
e
• 整理可得,整体平衡方程: K R
整体刚度矩阵的形式
整体平衡方程:
K R
1)其中[K]为将各单元的扩大矩阵迭加所得出的结构刚度矩
阵:
K K%(1) K%(2) K%(3) K%(4) K%(e)
集成包含搬家和迭加两个环节:
A、将单元刚度矩阵 K e 中的子块搬家,得出单元的扩
4、根据节点位移求出应力。
a
a
整体分析
1、建立整体刚度矩阵(也叫作结构刚度矩阵)
上图中的结构有六个节点,共有12个节点位移分量和12 个节点力分量。由结构的节点位移向量求结构的节点力向量时,
转换关系为: F K
分块形式为:
F 1 K11 K12 K13 K14 K15 K16 1
1
2
3
4
56
7 8 9 10
1
2
3
6
54
7 8 9 10
1
2
9
3 10 8
4 567
这需要了解单元中的节点编码与结构中的节点编码之 间的对应关系。
整体刚度矩阵的形式
结构中的节点编码称为
1
j1
a
节点的总码,各个单元的三 个节点又按逆时针方向编为 i,j,m,称为节点的局部码。

单元刚度矩阵中的子块
2 m1
i1 3
是按节点的局部码排列的,
j2 i3
m3 ③
j4
a
而结构刚度矩阵中的子块是 按节点的总码排列的。因此,
整体分析 及总体刚度矩阵的性质
整体分析
单元分析得出单元刚度矩阵,下面,将各单元组
合成结构,进行整体分析。
图示结构的网格共有四
1 Py1
个单元和六个节点。在节 点1、4、6共有四个支杆支
承。结构的载荷已经转移
a

Py3
2
Px2
a
3
Px3 ③

4

5
为结点载荷。 整体分析的四个步骤:
1、建立整体刚度矩阵; 2、根据支承条件修改整体 刚度矩阵; 6 3、解方程组,求节点位移;
② 4 m2
i2j3 m 4 ④ i4
在单元刚度矩阵中,把节点 的局部码换成总码,并把其
5
6 中的子块按照总码次序重新
排列。
a
a
整体刚度矩阵的形式
以单元②为例,局部码i,j,m对应于总码5,2,4,因此 k (2)
相关文档
最新文档