杆件结构的有限元法
第六章杆系结构
第六章杆件系统结构有限元法杆件系统是由几何特征为长度比横梁面的两个尺寸大很多的杆件连接而成的结构体系。
起重机械和运输机械的动臂、汽车的车架、钢结构等，都是由金属的杆件组成的。
杆件系统的有限元法在机械、建筑、航空、造船等各个工程领域得到了广泛的应用。
若杆件之间由铰相连，并且外载荷都作用在铰节点上，则该体系称为桁架。
有限元中将桁架的单元称为杆单元，即桁架是由仅承受轴向拉压的杆单元的集合。
如果杆件之间是由刚性连接，则该体系是刚架，刚架的单元称为梁单元。
梁单元可以承受轴力、弯矩、剪力及扭矩的作用。
第一节等截面梁单元平面刚架结构——所有杆件的轴线以及所有外力作用线都位于同一平面内，并且各杆件都能在此平面内产生平面弯曲，从而结构的各个节点位移都将发生在这个平面内。
一、结构离散化原则：杆件的交叉点、边界点、集中力作用点、位移约束点、分布力突变的位置都要布置成节点，而不同横截面的分界面和不同材料的分界面都要成为单元的分界面。
平面桁架对于桁架结构，因每个杆件都是一个二力杆，故每个杆件可设置成一个单元。
平面桁架结构每个节点有2个自由度，分别是u 和v ，每个单元有4个自由度。
最大半带宽B=(2+1)×2=6。
一维单元和二维单元的混合应用：左边部分是平面问题的二维板件结构（黑线部分），右面框架部分是一维杆件结构（红线部分）。
xy采用平面4节点四边形单元模拟二维板件，用平面杆单元单元模拟一维杆件结构。
离散化后，共有37个节点，32个单元，其中4节点四边形单元16个，杆单元单元16个。
因为平面4节点四边形单元和平面杆单元单元每个节点都有2个自由度，4节点四边形单元的刚度矩阵是8×8，平面杆单元的刚度矩阵是4×4。
整体刚度矩阵刚[]k 的维数是227474n n ⨯=⨯。
其中部分总刚子块为[](1)(2)(3)(4)777777777722k k k k k ⨯⎡⎤⎡⎤⎡⎤⎡⎤=+++⎣⎦⎣⎦⎣⎦⎣⎦(4)(6)(19)11,1111,1111,1111,1122k k k k ⨯⎡⎤⎡⎤⎡⎤⎡⎤=++⎣⎦⎣⎦⎣⎦⎣⎦最大半带宽B=［(8-2) +1］×2=14。
龙驭球《结构力学》笔记和课后习题(含真题)详解(矩阵位移法)【圣才出品】
第9章 矩阵位移法9.1 复习笔记一、矩阵位移法的基本思路矩阵位移法又称为杆件结构的有限元法。
分析的两个基本步骤：（1）单元分析；（2）整体分析。
单元分析：建立杆端力与杆端位移间的刚度方程，形成单元刚度矩阵。
整体分析：将单元合成整体，按照刚度集成规则形成整体刚度矩阵，建立位移基本方程。
二、单元刚度矩阵（局部坐标系）进行单元分析，推导单元刚度方程和单元刚度矩阵。
单元刚度方程是指由单元杆端位移求单元杆端力的一组方程，可以用“”表示，由位移求力称为“正问题”。
相应的由力求位移称为“反问题”。
正问题的解是唯一的确定的，但是反问题则可能无解，如果有解也非唯一解。
当外部荷载为不平衡力系时，反问题无解；当外荷载为平衡力系时，反问题有解但是因为杆件除本身变形外还可有任意刚体位移，此时反问题的解不唯一。
本书暂不考虑反问题的求解。
1．一般单元图9-1所示为平面刚架中的一个等截面直杆单元.单元的两个端点采用局部编码1和2，由端点1到端点2的方向规定为杆轴的正方向，在图中用箭头标明。
F →∆e图9-1图中采用坐标系，其中轴与杆轴重合。
这坐标系称为单元坐标系或者局部坐标系。
字母、的上面都画了一横，作为局部坐标系的标志。
推导单元刚度方程时，有以下几点需要注意：重新规定正负号规则、讨论杆件单元的一般情况、采用矩阵表示形式。
在局部坐标系中，图9-2所示的位移、力分量方向为正方向。
图9-2杆件性质：长度l ，截面面积A ，截面惯性矩I ，弹性模量E ；杆端位移u 、v 、θ。
根据杆端位移可以推导出下面两组刚度方程：（9-1）x y x x y（9-2）将上述六个刚度方程列成矩阵形式：（9-3）其中就是局部坐标系下单元刚度矩阵，即为（9-4）2．单元刚度矩阵的性质 （1）单元刚度系数的意义e e ek F∆=eK代表单元杆端第j 个位移分量等于1时所引起的第i 个杆端力分量。
（2）是对称矩阵，即。
（3）一般单元的是奇异矩阵，即，因此不存在逆矩阵。
杆件结构的有限元法
第一篇 有限元法
第二章 杆件结构的有限元法
当结构长度尺寸比两个截面方向的尺 寸大得多时，这类结构称为杆件。工程中 常见得轴、支柱、螺栓、加强肋以及各类 型钢等都属于杆件。
杆件结构可分为珩杆和梁两种。
和其他结构采用铰连接的杆称为珩杆。珩杆的连接处可以自由转动， 因此这类结构只承受拉压作用，内部应力为拉压应力。影响应力的 几何因素主要是截面面积，与截面形状无关。 和其他结构采用固定连接的杆称为梁。链的连接处不能自由转动， 因此梁不仅能够承受拉压，而且能承受弯曲和扭转作用。这类杆件 的内部应力状态比较复杂，应力大小和分布不仅与截面大小有关， 而且与截面形状和方位有很大关系。 建立有限元模型时，这两类杆件结构可用相应的杆单元和梁单元离散。
Ke 1 kkaa
ka
ka
中的元素在总刚度矩阵中应在位置第1行、第2行的第1列，第2列
k k
1 11
1 21
k
1 12
k
1 22
0
0
0 0 0
第2个单元的节点号为2和3，则单元刚度矩阵叠加到总刚度矩阵 的第2行、第3行的第2列、第3列元素上
0 0 0
0
k
2 22
k
2 23
0
k
2－3 杆件系统的有限元法
一、铰支杆系统的有限元计算格式 上面求解弹簧系统的有限元方法可以直接用力求解受轴向力的杆件系统。 均质等截面铰支杆，刚度值可由材料力学中力与变形的关系中获得
AE F1 L u1
k AE L
均质等截面铰支杆的力－位移方程可写为
F F12ALE11 11uu12
坐标变换
由杆件组成的机构体系称为杆系，如起重机、桥梁等。 由珩杆组成的杆系称为珩架，由梁组成的杆系称为刚架。
2_杆系结构有限元分析1
( x) Nii N j j
x x N 1 , N 其中 i 为形函数。 j l l
由材料力学扭转可知
d dN e e M GI p GI p θ GI p B θ dx dx
其中 B
dN 1 1 dx l l
§1-2 扭转杆单元
e
外力势能 V u
e
e T
fe
e
1 e T e e e T 总势能 U V u K u u f e 2
e e
§1-1 拉(压)杆单元
1 e T e e e T U V u K u u f e 2
e e e
根据最小势能原理，势能泛函取驻值的必要条件
空间杆单元坐标变换矩阵
0 T 0
单元在两个坐标系中刚度矩阵转换关系同样有
K e T T K ' T
e
矩阵中仅仅包含有坐标的倾角，仅平行移动坐标轴，刚度矩阵 中元素值不变，矩阵的阶数也不改变。
§1-2 扭转杆单元
结点位移向量θe i , j
T
结点力向量
平衡关系
杆单元结点力向量
f U i
e
Uj
T
单元在外力和内力作用下处于平衡状态，反映单元平衡状态 的关系式就是刚度方程。下面利用最小势能原理推导单元的 刚度方程。 最小势能原理：在满足连续条件和边界条件的位移中，满足 平衡条件的位移其总势能最小，反之亦然。 单元总势能
e U e V e
M e Mi , M j
T
杆件发生自由扭转时，待求位移是截面的扭转角 ( x) 在局部坐标系中，每一个点将具有一个基本未知位移，最简单 的单元位移函数可以设为
第五章杆系结构的有限元法
第五章 杆系结构的有限元法 5.1 引言杆系结构是工程中应用较为广泛的结构体系，包括平面或空间形式的梁、桁架、刚架、拱等。
其组成形式虽然复杂多样，但用计算机进行分析时却较为简单。
杆系结构中的每个杆件都是一个明显的单元。
杆件的两个端点自然形成有限元法的节点，杆件与杆件之间则用节点相连接。
显然，只要建立起杆件两端位移与杆端力之间的关系，则整体平衡方程的建立与前几章完全相同。
杆端位移与杆端力之间的关系，可用多种方法建立，包括前面几章一直采用的虚功原理，但是采用材料力学、结构力学的某些结论，不仅物理概念清晰、直观，而且推导过程简单明了。
因此，本章将采用这种方法进行单元分析。
至于整体平衡方程的建立，则和前面几章所讲的方法一样，即借助于单位定位向量，利用单元集成法进行。
5.2 平面桁架的有限元分析平面桁架在计算上有以下几个特点： 1． 杆件的每个节点仅有两个线位移； 2． 杆件之间的连接为理想铰，即在节点处各杆件可相对自由转动，且杆件轴线交于一点。
3． 外载荷均为作用于节点的集中力。
由于以上特点，所以在理论上各杆件只产生轴向拉、压力，截面应力分布均匀，材料可得到充分利用，因此桁架结构往往用于大跨结构。
5.2.1 局部坐标系下的单元刚度矩阵从平面桁架中任取一根杆件作为单元，称作桁架单元，单元长为L ，横截面面积为A ，图5.1。
两端节点分别用i 和j 表示，规定从i 到j 的连线方向为局部坐标x 轴，垂直于x 的方向为y 轴。
图5.1由于桁架中各杆只产生轴向力和轴向变形，所以节点i 和j 只发生沿x 方向的位移，用i u 和j u 表示，相应的杆端轴力分别用xi F 和xj F 表示。
由虎克定律可推得)()()(j i i j xj j i xi u u L EA u u L EA F u u LEAF --=-=-=将这两个式子写成矩阵形式，就是e j i exj xi u u L EA LEA L EA L EA F F ⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧ (5.1)显然，在局部坐标系下，i 、j 两节点沿y 轴方向的位移0==j i v v ，在y 轴方向的节点力0==yj yi F F 。
有限元法(杆系)
Fjy
FFji Fj
s in cos s in
s in
0 0
0 0 0
0
cos s in
或 F(e) T F (e) （1）
Fiy
i
Fi i
Fix
拉压杆单元
0 Fi e
0 0 0
0 Fj 0
F jy
j
j
uiy ui
uix
u jy
y
Fj
F jx uj
u jx
2)
叠加形成总刚度矩阵,求位移
2sin2
0
sin2 EA sin cos
l
0
0
sin2
sin cos
0 2 cos2 1 sin cos
cos2 0 1
sin cos cos2
sin2 sin cos
sin2 sin cos
0 0 0 0
sin cos cos2 sin cos cos2
• 用单元节点位移表示单元内部位移
第 i 个单元中的位移用所包含的结点位移来表示：
u(x)
ui
ui1 ui Li
(x
xi )
（1- 1）
其中 u i 为第 i 结点的位移， xi 为第 i 结点的坐标。
第 i 个单元的应变为 i ，应力为 i ，内力为 N i ：
i
du dx
ui1 ui Li
x
在局部坐标下，轴向力与轴向位移的关系：
(e)
Fi
1 0 1 0ui e
0
Fj
0
EA
0
0
l 1 0
0
0
0 1 0
0 0 0
有限元(第二章-杆单元部分)tg
−
1 2 1 2 1 2 1 − 2
−
1 2 1 2 1 2 1 − 2
1 2 1 − 2 1 − 2 1 2
按节点号叠加得6×6阶总刚度矩阵
−1 1 0 0 1 0 1 − 1 0 1 + 2 2 [K ] = 0 0 − 1 2 2 0 0 − 1 2 2 1 0 −1 2 2 0 0 1 − 2 2 1 2 2 1 2 2 1 − 2 2 0 0 0 −1 1 1 − 2 2 2 2 1 1 2 2 2 2 1 1 − 2 2 2 2 1 1 − 1+ 2 2 2 2
2-10 刚度矩阵元素的带状分布
【例】对图（a）中结构分别采用图（b）、图 （c）两种编号方式以观察其刚度矩阵的带宽。
对于图(b)、(c) 编号方式的结构，总刚度矩阵 的非零元素分布分别如下图(a)、(b) 所示。
[K ]
e
λ2 AE λµ = L − λ2 − λµ
λµ µ2 − λµ − µ2
Fx1 1 Fy1 AE 0 = L − 1 F x2 Fy 2 0
即：
0 − 1 0 u1 0 0 0 v1 0 1 0 u 2 0 0 0 v 2
{F }= [K e ]{δ }
求各杆单元的λ和μ的值。Φ角是按 逆时针从x轴正向转到单元ij方向的
三杆受力桁架
单元⑴ 单元⑵ 单元⑶
ϕ = 0 o , λ = 1, µ = 0 ϕ = 90 o , λ = 0 , µ = 1 ϕ = 135 o , λ = −
1 1 ,µ = 2 2
单元刚度矩阵分别为
杆结构 分析的有限元方法(有限元)
杆单元形状函数
杆单元刚度矩阵
平面问题中的坐标变换
梁结构分析的有限元方法
梁：承受横向荷载和弯矩的杆件。
梁的主要变形为挠度v
横截面变形前后都垂直于杆变形前的轴线x轴
中性层变形=0
纯弯曲没有剪力，只有弯矩
梁截面的惯性矩
杆结构分析的有限元方法
杆：承受轴向荷载的杆件
最基本的承力结构件：杆、梁
弹簧--简单的承受轴力的结构件
有限元方法中，每一个处理步骤都是标准化和规范化的，
因而可以在计算机上通过编程来自动实现。
F=kδ
k--刚性系数
位移的绝对变化量/杆件的伸长量δ=u2—u1
应力某截面上单位面积上的内力/内力的分布集度
应变相对伸长量单位长度的伸长量
杆单元的特性是节点位移及节点力的方向都是沿轴线方向。
杆结构的力学分析
铰接的杆结构----杆只受轴力-----杆件拉伸问题---可自然离散
两端为铰接的杆件只承受轴力。
各个单元研究（基于局部坐标系的表达)
各个单元研究
离散单元的集合、组装
杆单元及坐标变换
自由度：描述物体位置状态的每个独立变量。
对于杆单元，其节点位移有两个自由度。
杆件屈曲 有限元
杆件屈曲是指在受力作用下，杆件发生稳定失稳的现象，通常包括弹性屈曲和塑性屈曲两种情况。
有限元方法是一种数值分析方法，可以用来模拟和分析结构在受力作用下的变形、应力分布等问题。
在研究杆件屈曲问题时，有限元方法也可以被广泛应用。
在有限元分析中，研究杆件屈曲通常包括以下步骤：
1. 建立模型：根据实际情况和要求，建立杆件的有限元模型，包括定义几何形状、材料力学性质、约束条件等。
2. 网格划分：将杆件模型离散为有限个单元，通常采用三角形或四边形等简单几何形状的单元，并确保单元之间的连接和边界条件设置正确。
3. 施加载荷：在模型中施加适当的载荷或边界条件，模拟实际工程中的受力情况。
4. 求解：通过有限元求解器对模型进行计算，得到杆件在受力下的应力、变形等信息。
5. 后处理：分析求解结果，确定杆件是否处于屈曲状态，了解屈曲位置和形式，评估结构的稳定性。
在研究杆件屈曲问题时，有限元分析可以帮助工程师更好地理解结构在复杂加载条件下的行为，预测结构的稳定性和安全性，优化设计方案，减少实验成本和时间。
通过有限元模拟，可以有效地探索杆件的屈曲特性，为工程实践提供重要参考和支持。
杆件系统有限单元法
（3）单元应力场的表达 由弹性力学中物理方程有：
σ e ( x ) = E eε e ( x ) = E e B e ( x ) ⋅ δ e = S e ( x ) ⋅ δ e
其中Se为单元的应力函数矩阵：
⎡ E S ( x) = E B ( x) = ⎢ − ⎣ l
e e e
e
E ⎤ ⎥ l ⎦
平面梁单元的节点位移δe和节点力Fe为：
δ =⎡ ⎣ui vi θi u j v j θ j ⎤ ⎦
e e
T
F =⎡ ⎣ FNi FQi M i FNj FQj M j ⎤ ⎦
相应的刚度方程为：
T
K e ⋅δ e = F e
将杆单元刚度矩阵与纯弯梁单元刚度矩阵进行组 合，可得到平面梁单元的刚度矩阵：
可以写出节点位移向量和节点力向量：
δ =⎡ ⎣ui u j ⎤ ⎦
e
e
T
T ⎡ ⎤ F = ⎣ FNi FNj ⎦
（1）单元位移模式的表达 由于每个节点只有一个轴向位移，即一个单元共有 两个自由度，因此可假设该单元的位移模式为具有 两个待定系数的函数模式：
u ( x ) = a 0 + a1 x
e
第三章
杆件结构的有限元分析 (FEA)
在杆件系统中根据单元受力的特点，我们可以 把它们分成两大类：杆和梁。为了以后描述的 方便，我们把两端铰接，只受轴向力的基本结 构称为杆单元，而受轴向力和弯矩、扭矩、剪 力共同作用的基本结构称为梁单元。
3.1 平面杆单元
局部坐标系中的杆单元描述
设有一任意的杆单元如图所示，i 和j 为单元的两 个结点，x 为该单元的局部坐标，其原点设在单 元的i 结点。设两个结点在x 方向的位移为 u i 和 u j ，它们的正方向如图3-1 所示，与它们相应的 结点力 FN δ e
《有限元理论与数值方法》第三讲-杆、梁结构有限元分析
Finite Element Theory and Numerical Method
一、杆、梁的物理力学模型
拉压杆单元如图3-6所示，已知等直杆件杆长为 l 横截面面积为 A 材料弹性模量为 E 所受轴向分布载荷集度为 p(x) 杆端位移分别为 u1 u2
杆端力分别记为 F1 F2
1、建立位移场
F1, u1 xa
1
a p(x)
2 F2 , u2
x
设局部坐标系下杆中任意点a的坐标为 xa
因为只有两个边界条件 u1
形函数具有如下性质： 1)本端为1，它端为0 2)单元内任意一点总和为1
N1(0) 1
N1(1) 0
N2 (0) 0 N2 (1) 1
N1() N2 () 1
2、应变分析
du dx
dN dx
ue
dN1 dx
B为应变矩阵或者几何矩阵。
dN2 dx
u
e
1 l
1 l
ue
[B1
B2 ]ue Bue
图示所示桁架 l 2m
EA 1.2106 kN
试求1-2杆和1-4杆单元的局部坐标单元 刚度矩阵
1-2杆：抗拉刚度 EA / l 6106 kN/m
F1 10N 3
1
F2 20N 4
2
ke1
EA l
1 1
1
1
6
105
1 1
1
1
kN
/
m
1-4杆：抗拉刚度 EA /( 2l) 4.24264 105 kN/m
有限单元法课件第四章 杆件系统的有限元法
(a)
(b)
由杆件组成的结构体系称为杆系,如起重机,桥梁等。
由桁杆组成的杆系称为桁架。
由梁组成的杆系成为刚架。
若杆系和作用力均位于同一平面内,则称为平面桁架 或平面刚架,否则称为空间桁架或空间刚架。
由于杆件结构采用一维单元进行离散,所以杆系的网 格划分容易用半自动方法实现。当采用自动网格划 分方法时,杆系的几何模型是由杆件轴线构成的线框 模型。
R
e P
RiP R jP
R
lP
R
R
e F
RiF R jF
Rlx Rly NlT l R l
lF T l
Px dx (l i, j ) Py
e T
Bj dx
kii k ji
kij k jj
其中矩阵元素为
kst D Bt dx B as 0 EA 0 at 0 0 0 bs dx 0 EI 0 bt ct 0 cs 0 0 EAas at dx 0 EIb b EIb c s t s t EIcs bt EIcs ct 0
e
du dx e x 2 B Bi q x d v dx 2
Bj q
e
其中
ai 0 0 Bi 0 b c i i a j 0 0 Bj 0 b c j j 1 12 6 ai a j bi b j 3 x 2 l l l 4 6 2 6 ci 2 x cj 2 x l l l l
杆件结构的有限元法PPT课件
2 2
K e
EA
2
L 2 2
2
EA k e k e
L k e
k
e
其中：k e
2
2
2
2
第28页/共33页
求解整体坐标系下结构受力与位移方程组：
F K
可得到各节点位移，从而可以求出每根杆的 受力，简单推导可得：
pij
EA L
,
ij
单元1：FF12
ka ka
单元2：FF32
kb kb
ka ka
uu12
kb kb
uu32
第12页/共33页
(2)由于整个系统有3个节点，扩充上述方程为3阶：
F1 F2
ka ka
ka ka
00uu21
F3 0 0 0u3
F1 F2
kb kb
** **
行
2j-1 2j
** **
** **
第30页/共33页
刚度矩阵的性质： （1）对称性——关于主对角线对称； （2）稀疏性——大量0元素； （3）带状分布——非0元素在主对角线两侧 呈带状分布。 所以可以对总体刚度矩阵进行压缩存储。方法 是：找出所有各行中非0元素所占最宽一行， 以离对角线最远的元素为基准画一条平行于主 对角线的带子，称为其带宽，方法称为等带宽 存储。由于对称性，带宽的一半称为半带宽。
• （1）形成每个单元刚度矩阵； • （2）由各单元的刚度矩阵按节点号叠加
整个系统的刚度矩阵；
• （3）引入约束条件； • （4）以节点位移为未知量求解线性方程
组
• （5）用每个单元的力-位移关系求的单元
第18页/共33页
第三节 杆件系统的有限元法 简单拉（压）杆的受力特点为作用在直杆 上的外力（体力、面力）合力的作用线一定与 杆的轴线重合，如图所示。
结构有限元法(绪论)
有限单元法的应用 有限单元法在应用上已远远超过了原来的范围。
它已由弹性力学平面问题扩展到空间问题和板壳问题， 能对原子能反应堆、拱坝、飞机、船体、涡轮叶片等 复杂结构进行应力分析；它已出平衡问题扩展到稳定 问题与动力问题，由弹性问题扩展到弹塑性与粘弹性 问题，由结构的应力分析扩展到结构的优化设计。除 此，它在流体力学、热传导、磁场、建筑声学、 生物力学等等方面部有不同程度的应用。
近几十年来，随着电子计算机的高速化和普遍化， 有限元继续不断地向更加广阔、更加深入的方面发展。
有限单元法的发展借助于两个重要工具，在理论 推导方面，采用了矩阵方法，在实际计算中，采用了 电子计算机。有限元、矩阵、计算机是三位一体的。 由于有了现代化的、先进的计算工具，使得有限单元 法近年来以惊人的速度骤然崛起。
又有效的数值方法。
有限单元法的发展历史 有限单元法最初是在五十年代作为处理固体力学
问题的方法出现的。 追溯历史，早在一九四三年，库兰特(courant)已
应用了“单元”概念。在一九五六年，特纳(Turner) 等人把刚架位移法的解题思路，推广应用于弹性力学 平面问题。他们把连续体划分成一个个三角形的和矩 形的单元，单元中位移函数首先采用了近似表达式， 推导了单元刚度矩阵，建立了单元结点位移与结点力 之间的单元刚度方程。
即物体在引起形变的外力被除去以后，能够完 全恢复其原来的形状，这种性质称为“弹性”。如 果材料又服从虎克定律，即外力与变形之间的关系 成正比，这种弹性就叫做“线性弹性”。在这一假 定下的物体只能发生线性弹性变形。
（2）假设物体是连续的 这假设认为整个物体的体积都被组成这个物体
的物质所填满，而不留下任何空隙。这样物体中应 力、应变和位移等等物理量就可看成是连续的，因 而我们就可用坐标的连续函数来表示它们的变化规 律。
杆系结构的有限元法分析
杆系结构的有限元法分析有限元法是一种结构分析方法，常用于分析各种不同类型的结构系统，其中包括杆系结构。
杆系结构是由杆件连接而成的桁架结构，常见于桥梁、塔架和支撑结构等。
利用有限元法进行杆系结构的分析，可以得到结构的位移、应力、应变和刚度等信息，帮助工程师评估结构的稳定性和安全性。
下面将介绍杆系结构的有限元法分析的步骤。
首先，进行前期准备工作。
这包括收集与结构相关的几何信息（如杆件长度、截面形状等）、边界条件（如固定支座、外载荷等）和材料性质（如材料的弹性模量、密度等）。
这些信息将是有限元模型建立所需要的输入参数。
接下来，建立有限元模型。
将杆系结构离散化为一个个的杆单元，采用有限元方法对每个杆单元进行离散近似。
常用的杆单元包括横截面线性杆单元、三节点弯曲杆单元和非线性杆单元等。
然后，确定单元刚度矩阵。
对于横截面线性杆单元，其刚度矩阵可以根据材料性质和几何信息计算得到。
对于弯曲杆单元和非线性杆单元，则需要考虑附加的几何和材料非线性效应。
接着，组装全局刚度矩阵。
将所有杆单元的刚度矩阵按照其关联的节点自由度进行组装。
在组装过程中，需要考虑杆单元之间的关联关系，确保刚度矩阵的正确性和完整性。
然后，应用边界条件。
根据实际情况，将已知的边界条件（如固定支座、已知位移等）施加到全局刚度矩阵中。
这将改变全局刚度矩阵的特征值和特征向量，从而影响结构的响应。
接下来，求解结构的位移和应力。
通过求解结构的整体刚度方程以及施加的边界条件，可以得到结构的位移解向量和应力解向量。
位移解向量描述了结构的变形情况，而应力解向量体现了结构的应力分布情况。
最后，进行后处理。
在得到位移和应力解后，可以计算结构的应变分布、变形形态以及额外的设计指标。
通过这些结果，可以对结构的性能进行评估，以便优化设计。
综上所述，杆系结构的有限元法分析包括前期准备、建立有限元模型、确定单元刚度矩阵、组装全局刚度矩阵、应用边界条件、求解结构的位移和应力以及后处理等步骤。
有限单元法 第2章 杆系结构的有限元法分析
义 & 可以进一步求得单元刚度矩阵为 )
( & # 0# ( $’ $ % 8 . ! 1 # $ ’ 0# # 同时 & 我们可以根据式 $ % 求出等 效 结 点 荷 载 矩 阵 ’ 这 里 要 指 出 的 是 ) 分 布 荷 载 ! .$
! # !! !
! # $! !
! 第 ! 章 ! 杆系结构的有限元法分析 # #! ! """""""""""""""""""""""""""""""""""""""""""
不适定的 " 第九步 # 求解方程组 " 计算结构的整体结点位移列阵 ## 并 进一步 计算各 单元 的应力 分量及主应力 $ 主向 " 第十步 # 求单元内力 # 对计算成果进行整理 $ 分析 # 用表格 $ 图线标示出所需的位移 及应力 " 大型商业软件 % 如 )* + , + 等 & 一般都具有强大的后处理功能 # 能够 由计算 机自 动绘制彩色云图 # 制作图线 $ 表格乃至动画显示 "
矩阵 ’ $ %进行应力 ( 应变分析 ’ 根据材料力学中应变的定义 & 有 ) ! # # $’ 2 + 2 $ ( ( ( ( $’ $’ $’ . 0 ! ! . " 3 3 .% ". . ! ! ! !! "# ’ ’ 2 # 2 #
第三讲 杆件结构有限元分析
l
0
AE
l du d u dx f x udx 0 dx dx
其中E表示弹性模量，A表示横截面积，方程左端得到单元的刚度矩阵。
建立有限元模型
现考虑一个由5个长度相同(le=1m)横截面积不同的杆件构成的一维杆件，各杆弹性模量都为 E=1.0e10pa，A1=0.5m2，A2=0.4 m2，A3=0.3 m2，A4=0.2 m2，A5=0.1 m2，如图1所示，右端给定位移 u右=0.1，左端固定位移u左，分析杆件内位移分布：
根据虚功原理，方程两边乘以虚位移δu，平衡方程可以写为：
其弱形式为：
l
0
[
d ( A x ) f ( x)] udx 0 dx
l
0
A x
l d u dx f x udx Pj u j 0 dx j
基本方程的最终弱形式
其中，右端最后一项可以看作是节点力情况，所以可以不单独列出，同时 x E 所以上式可以继续写为：
网格尺寸设置
网格划分信息
网格划分
选择calculate → calculate，在弹出的对话框，点击OK，保存，前处理完毕。
工程求解
点击工具栏中“求解计算”按钮，完成模型的求解计算。
后处理
点击工具栏中的“后处理”按钮进入GID，查看计算结果，如下图所示。
结果分析： 本章针对一个变截面一维杆件，通过理论分析和ELAB1.0软件实现两种方式来分析，一方面对有限元
几何模型
将其划分为五个单元六个节点，即每根杆件作为一个单元，每个单元的节点关系如下图所示：
单元拓扑关系
确定杆单元的形函数
考虑其中一个杆单元，其两个端点分别为节点1，节点2，基本变量为节点位移u1，u2:：
有限元方法第三章杆系结构有限元
应用实例
某大型桥梁的稳定性分析
采用杆系结构有限元对某大型桥梁进行稳定性分析，评估其在不同载 荷下的变形和承载能力。
高层建筑的抗震性能研究
利用杆系结构有限元模拟高层建筑的抗震性能，分析地震作用下结构 的响应和破坏模式。
汽车悬挂系统的优化设计
通过杆系结构有限元模拟汽车悬挂系统的运动和受力情况，优化悬挂 参数以提高车辆行驶的稳定性和舒适性。
有限元方法第三章杆系结 构有限元
• 引言 • 杆系结构有限元的基本概念 • 杆系结构有限元的建模方法 • 杆系结构有限元的求解方法 • 杆系结构有限元的应用案例 • 结论与展望
01
引言
目的和背景
杆系结构是工程中常见的一种结构形式，广泛应用于桥梁、 建筑、机械等领域。由于其具有复杂的几何形状和受力特性 ，因此需要采用有限元方法进行数值分析。
THANKS
感谢观看
04
杆系结构有限元的求解方法
求解步骤
确定边界条件
根据实际情况，确定杆系结构 的边界条件，如固定、自由、 受压等。
求解线性方程组
将所有单元的平衡方程组合成 一个线性方程组，然后使用数 值方法求解该线性方程组。
建立离散模型
首先将杆系结构离散化为若干 个小的单元，每个单元具有一 定的物理属性。
应用力学平衡方程
杆系结构有限元的优缺点
优点
能够处理复杂的几何形状和边界条件， 适用于大规模问题求解，计算精度可 调，可模拟复杂的结构和场。
缺点
需要针对不同的问题建立不同的模型， 计算量大，需要较高的计算机资源， 对于非线性问题求解较为困难。
03
杆系结构有限元的建模方法
建模步骤
确定研究问题
杆件结构有限元分析
1.0e10 0.1 1 1 u6 1 1 u 1 5 u3 u2
1.0e10 0.2 1 1 u5 1 1 u 1 4
1.0e10 0.3 1 1 u4 1 1 u 1 3 u2 u1
u( x) N1u1 N2u2
通常来说形函数需满足以下条件：
• 在单元内，一阶导数必须存在；
• 在单元之间的连接处，位移必须连续；
推导该问题的单元刚度矩阵表达式
对于基本方程的最终弱形式：
l du d u dx f x udx 0 dx dx
l
0
AE
在一个单元内： u( x) N1u1 N2u2 因此：
K
(4)
1.0e10 0.4 1 1 u3 1 1 u 1 2
K
(5)
1.0e10 0.5 1 1 u2 1 1 u 1 1
确定该问题总体刚度矩阵
将得到的各个单元刚度矩阵按节点编号进行组装，可以形成整体刚度矩阵，同时将所有节点荷载也 进行组装。 刚度矩阵： K K (1) K (2) K (3) K (4) K (5)
le x2 x1
由上面的关系式可得到：
1 1 u1 dN 2 du dN1 u1 u2 dx dx dx le le u2
对δu取与u相同的形函数，并将上面的关系式带入基本方程的最终弱形式中可得：
1 le l du d u AE dx u u 1 2 0 dx dx 1 l e 1 AE le 1 l 1 e u1 u2 AEle 1 le l e 1 u1 le 1 1 d le u2 2
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
可编辑ppt
5
由杆件组成的机构体系称为杆系，如起重机、桥梁等。 由珩杆组成的杆系称为珩架，由梁组成的杆系称为刚架。
可编辑ppt
6
奥运会场馆
鸟巢
空间立体网架
可编辑ppt
7
2-1 引 言
工程中最简单的结构可以认为是铰支的杆件。它的性质完全类似于弹簧。
弹簧系统力F与弹簧伸长量 （位移）之间关
系由胡克定律有
FF12
k11 k21
k12 k22
kk1233uu12
F3 k31 k32 k33u3
利用线弹性系统的叠加原理，找出3×3阶刚度矩阵各元素 的表达式
节点1处的合力 节点2处的合力
节点3处的合力
F1kau1 F2kau1 F30
kau2 kau2kbu2 kbu2
0 kbu3 kbu3
ka
Kka
0
刚度应采用一个矩阵来表示，即 K ，同理，各点的位移也应采用一个
矩阵来表示，即 ，再加上矩阵 F ，就构成了
FK K 称为对应于施加存系统上各节点力的刚度矩阵。
可编辑ppt
9
问题： 1、复杂结构其刚度矩阵是多少阶的？ 2、如何求出？ 3、为什么着重讨论系统的刚度矩阵？
系统的整体刚度矩阵－求出所受外力作 用下各杆件节点处的位移－计算各杆件的 受力和应力
ka ka kb
kb
0 kb kb
对成、奇异矩阵
（2－8）
可编辑ppt
14
用同样的方法可以求解具有更多个弹簧 的串连系统，推导过程乏味。
知道单个弹簧的刚度矩阵－－直接叠加 出多个串联系统的总刚度矩阵。
可编辑ppt
15
知道单个弹簧单元的刚度矩阵，直接叠加出总刚度矩阵
对整个系统来说有3个节点，将上述方程扩大成3阶方程：
F1c
ka F2c kb
u3,F3c
3) 只允许节点3有位移u3，类似于情况1），有
F 3 c k b u 3 ,F 2 c F 3 c k b u 3
u1＝0
u2＝0
由于节点1、2无位移，有
(c)
可编辑pptF1c 0
13
组合弹簧的刚度矩阵
4) 合成。对整个系统来说有3个节点，每个节点只有一个 方向的位移。因此方程式应用如下形式：
的结构，如图（c），叠加结果为：
A A‘
(c) B B‘
作用于节点1上的合力 作用于节点2上的合力
F1 F1a F1b F2 F2a F2b
刚度矩阵
F1 F2
k k
k k
uu12
（2－5）
Ke
k k
k k （2－6）
对成、可奇编异辑矩pp阵t
12
二、组合弹簧的刚度矩阵
u1,F1
u u
1 2
F3 0
0 0 u 3
F1 F2
F3
0
0
0
0 kb kb
0 k kb
b
u u u
1 2 3
矩阵扩大办法
单元数量增多时，相应扩大后的矩阵 就相当大，扩大后的非零元素在矩阵 的什么位置，概念上就不很清楚了。
按矩阵相加原理将两式叠加，
可编辑ppt
11
1)只有节点1可以变形,点2固定
F1a ku1
由力的平衡有
F1a
u1
k
F2a
F1a F2a 0 F2a F1a ku1
2)只有节点2可以变形,点1固定
F 2bk2 uF 1b
A A‘ (a)
F1b
k
u2=0
u2
F2b
u1=0 F1 u1
B B‘
(b)
k
u2 F2
3)根据线弹性系统的叠加原理，叠加1) 2)两种情况，就得到与原始问题一样
F3a 0
2) 只允许节点2有位移u2，这时由于位移的连续性，每个 弹簧在节点2要求有相同的位移，即，弹簧1-2的伸长量与
弹簧2-3的缩短量相等。对弹簧1-2 有拉力kau2，对弹簧 2-3 有压力kbu2
F 2bkakbu2
分别对两弹簧求静力平衡，有 F 1 b kau 2,F 3 b kb u 2
u2,F2
ka
kb
u3,F3
u1,F1a F1b
1
2
3
ka
F2a kb
F3a
u2＝0
(a)
3
u3＝0
ka u2,F2b kb
F3b
u1＝0
u3＝0
(b)
1) 只允许节点1有位移u1，力F1a与位移u1之间的关系
F1a kau1 考虑弹簧1-2，由静力平衡条件有
F 2aF 1akau1 由于u1＝ u2＝0，没有力作用于节点3，因此，
有限元理论与应用
可编辑ppt
1
第一篇 有限元法
可编辑ppt
Hale Waihona Puke 2第一篇 有限元法
第二章 杆件结构的有限元法
可编辑ppt
3
当结构长度尺寸比两个截面方向的尺 寸大得多时，这类结构称为杆件。工程中 常见得轴、支柱、螺栓、加强肋以及各类 型钢等都属于杆件。
可编辑ppt
4
杆件结构可分为珩杆和梁两种。
和其他结构采用铰连接的杆称为珩杆。珩杆的连接处可以自由转动， 因此这类结构只承受拉压作用，内部应力为拉压应力。影响应力的 几何因素主要是截面面积，与截面形状无关。 和其他结构采用固定连接的杆称为梁。链的连接处不能自由转动， 因此梁不仅能够承受拉压，而且能承受弯曲和扭转作用。这类杆件 的内部应力状态比较复杂，应力大小和分布不仅与截面大小有关， 而且与截面形状和方位有很大关系。 建立有限元模型时，这两类杆件结构可用相应的杆单元和梁单元离散。
F F 1 2 k k a a k k a a u u 1 2 F F 3 2 k k b b k k b b u u 3 2
整个系统有3个节点（位移），将上述方程扩大成3阶方程，
F1 F2
ka
k
a
ka ka
0 0
Fk (4—1)
式中k为弹簧的刚度，是弹簧的固有参数。它对应于
力－位移图中F- 关系直线的斜率。
当k和力F已知时，可由下式求出弹簧伸长量
1F
k
弹簧力－位移间关系
可编辑ppt
8
当处理比较复杂的铰支杆系统时，要确定系统在力P的作用下，节点B、 C、D和E处的变形。以便计算各杆件的内应力及各杆所受的轴向力，可 假设整个杆件系统也具有像式(4—1)中k值一样的刚度，这样在力P的作 用下各点的位移就可以用类似式(4—1)的公式计算了，不过．这时的系统
可编辑ppt
10
2-2 弹簧系统的刚度矩阵
一、单个弹簧的刚度矩阵
u1,F1
k
u2,F2
弹簧的作用力向量为
F F
1 2
位移向量为
u u
1 2
FF12kk1211 kk1222uu12
从而这个弹簧的刚度矩阵是2x 2阶的。
为求出它们，将图2—4所示弹簧系统看作两个简单的系统，然后合成。