连续梁的矩阵位移法优秀课件
合集下载
结构力学-矩阵位移法-PPT
a11 AB a21
当p=l时才能相乘
a12 b11 a22 b21
a12 a22
共形
b11 a11 BA b21 a21
非共形
(2)不具有交换律,即 AB BA
6、转置矩阵 将一个阶矩阵的行和列依次互换,所得的阶矩 阵称之为原矩阵的转置矩阵,如:
任意矩阵与单位矩阵相乘仍等于原矩阵,即 AI =A IA =A
10、逆矩阵
在矩阵运算中,没有矩阵的直接除法,
除法运算由矩阵求逆来完成。例如,若
AB = C
则
B=A 1 C
-
此处A-1 称为矩阵A的逆矩阵。
一个矩阵的逆矩阵由以下关系式定义:A A 1 = A 1 A =I
矩阵求逆时必须满足两个条件: (1)矩阵是一个方阵。 (2)矩阵的行列式不为零,即矩阵是非奇异矩阵(行列 式为零的矩阵称为奇异矩阵)。
矩阵位移法(刚度法):
结点力
P
F
(物理条件)
结点位移
(几何条件)
(平衡条件)
杆端力
杆端位移
r11 z1 r12 z 2 L r1i zi R1p 0 r21 z1 r22 z 2 L r2i zi R2p 0 r31 z1 r32 z 2 L r3i zi R3p 0
结构力学
STRUCTURE MECHANICS
第十章
矩阵位移法
知识点:
• • • • 矩阵位移法的基本要点 常见单元单元刚度矩阵的建立 单元刚度矩阵的坐标变换 矩阵位移法计算连续梁和刚架
教学基本要求:
掌握矩阵位移法的基本要点;
理解各种常见单元杆端位移和杆端力的对应 关系,理解单刚矩阵的建立方法及过程,能正确 写出常见单元的单刚方程;理解坐标变化的意义 及方法。 掌握前处理法计算连续梁和不考虑轴线变形 的刚架,结合刚架理解后处理法的基本思想。
当p=l时才能相乘
a12 b11 a22 b21
a12 a22
共形
b11 a11 BA b21 a21
非共形
(2)不具有交换律,即 AB BA
6、转置矩阵 将一个阶矩阵的行和列依次互换,所得的阶矩 阵称之为原矩阵的转置矩阵,如:
任意矩阵与单位矩阵相乘仍等于原矩阵,即 AI =A IA =A
10、逆矩阵
在矩阵运算中,没有矩阵的直接除法,
除法运算由矩阵求逆来完成。例如,若
AB = C
则
B=A 1 C
-
此处A-1 称为矩阵A的逆矩阵。
一个矩阵的逆矩阵由以下关系式定义:A A 1 = A 1 A =I
矩阵求逆时必须满足两个条件: (1)矩阵是一个方阵。 (2)矩阵的行列式不为零,即矩阵是非奇异矩阵(行列 式为零的矩阵称为奇异矩阵)。
矩阵位移法(刚度法):
结点力
P
F
(物理条件)
结点位移
(几何条件)
(平衡条件)
杆端力
杆端位移
r11 z1 r12 z 2 L r1i zi R1p 0 r21 z1 r22 z 2 L r2i zi R2p 0 r31 z1 r32 z 2 L r3i zi R3p 0
结构力学
STRUCTURE MECHANICS
第十章
矩阵位移法
知识点:
• • • • 矩阵位移法的基本要点 常见单元单元刚度矩阵的建立 单元刚度矩阵的坐标变换 矩阵位移法计算连续梁和刚架
教学基本要求:
掌握矩阵位移法的基本要点;
理解各种常见单元杆端位移和杆端力的对应 关系,理解单刚矩阵的建立方法及过程,能正确 写出常见单元的单刚方程;理解坐标变化的意义 及方法。 掌握前处理法计算连续梁和不考虑轴线变形 的刚架,结合刚架理解后处理法的基本思想。
结构力学第五版第十章矩阵位移法ppt课件
k12
p3 k31 k32 k33 3
k112
简记为 P k---结构刚度方程
k21 k31
k 211 =1 k22
1
k
1 22
1
k32
2
k
2 21
k --结构刚度矩阵(总刚)
k11 k111 k21 k211
k31 0
k13 k121
k23 k33
=1
3
k12 k112 k22 k212 k121 k32 k221
1 2 3
6 3 P3
3 (P3 01 4 2 ) /(8 N ) 3 0
六.非结点荷载
(1).等效结点荷载
PE
PPEE12
PE3
PE1
PE 2
PE 3
---结构等效结点荷载
“等效”是指等效结点荷载引起的结点 位移与非结点荷载引起的结点位移相同
(2).等效结点荷载的计算
1
4
6/ 1.5
8
1.5 1 1
3
2
2
EI1 6 EI 2 24
4m 4m 12m
1
2
1
2
EI1 6
8m
34
3
2
3
1
2
k 2
4
24 4
/12
4 1 2 8 2 3
34
12
k
3
3 1.5
1.5 1 3
3
2
4
3 1.5 0 0
k 1.5 11
4
0
0 4 11 1.5
0
0
1.5
1 2
Fq
2
ql 2 /12 ql2 /12
矩阵位移法ppt课件
e
i
u j , Fxj
e
vi , Fyi
0 0 0 0 EA l 0 EA l 0
v j , Fyj
0 0 0 0
ui v i u j v j
e
杆端力向量
单刚矩阵
杆端位移向量
22
坐标变换
上述单刚方程是在单元坐标系下建立的,单元按结点平衡拼装成结构之前, 由于结构中单元的方位一般不全相同,因此,应将杆端位移和杆端力都转换 成统一的、对整体坐标的量,这是同一矢量在不同坐标系中的变换问题,简 称为坐标变换。 坐标变换矩阵 —单元系与结构系的关系: F x yj y 结构系 x 轴沿逆时针转至单元系 y x 轴所转过的角度记为 。 Fxj j Fyi —单元系下的杆端力 Fxj e —结构系下的杆端力 Fyj i Fxi e e Fxi Fxie Fxi cos Fyi sin Fyi 表示为 e e e F F sin F cos yi xi yi x o e e Fxje Fxj cos Fyj sin 矩阵形式 e e Fyje Fxj sin Fyj cos
原始总刚度方程无法直接求解。因为原始总刚度方程表示结构全部结点的平 衡方程,结构的结点分为两大类:有约束的支座结点(如图示结构的1、2结 点)和无约束的内部结点(如图示结构的3结点);前类结点已知结点位移而 未知结点力(支座反力),后类结点已知结点力而未知结点位移。 M 要求支座反力必须先求后类结点的结点位移,故,原 P 3 始总刚度方程必须考虑边界位移条件修正为结构刚度 I, A I, A 方程后才能求解。 2 根据支座位移边界条件的处理方式不同,矩阵位移法 1 可分为先处理法和后处 理法。
i
u j , Fxj
e
vi , Fyi
0 0 0 0 EA l 0 EA l 0
v j , Fyj
0 0 0 0
ui v i u j v j
e
杆端力向量
单刚矩阵
杆端位移向量
22
坐标变换
上述单刚方程是在单元坐标系下建立的,单元按结点平衡拼装成结构之前, 由于结构中单元的方位一般不全相同,因此,应将杆端位移和杆端力都转换 成统一的、对整体坐标的量,这是同一矢量在不同坐标系中的变换问题,简 称为坐标变换。 坐标变换矩阵 —单元系与结构系的关系: F x yj y 结构系 x 轴沿逆时针转至单元系 y x 轴所转过的角度记为 。 Fxj j Fyi —单元系下的杆端力 Fxj e —结构系下的杆端力 Fyj i Fxi e e Fxi Fxie Fxi cos Fyi sin Fyi 表示为 e e e F F sin F cos yi xi yi x o e e Fxje Fxj cos Fyj sin 矩阵形式 e e Fyje Fxj sin Fyj cos
原始总刚度方程无法直接求解。因为原始总刚度方程表示结构全部结点的平 衡方程,结构的结点分为两大类:有约束的支座结点(如图示结构的1、2结 点)和无约束的内部结点(如图示结构的3结点);前类结点已知结点位移而 未知结点力(支座反力),后类结点已知结点力而未知结点位移。 M 要求支座反力必须先求后类结点的结点位移,故,原 P 3 始总刚度方程必须考虑边界位移条件修正为结构刚度 I, A I, A 方程后才能求解。 2 根据支座位移边界条件的处理方式不同,矩阵位移法 1 可分为先处理法和后处 理法。
《矩阵位移法》课件
实际工程案例分析
总结词
为了验证矩阵位移法的有效性,可以通过实际工程案例 进行分析。通过与实验结果的对比,可以评估方法的精 度和可靠性。
详细描述
选取具有代表性的实际工程案例,如高层建筑、大跨度 桥梁等,利用矩阵位移法进行计算,并将结果与实验数 据进行对比。通过对比分析,可以评估矩阵位移法的精 度和可靠性,为该方法在实际工程中的应用提供依据。 同时,也可以针对不同工程案例的特点,对矩阵位移法 进行优化和改进,提高其适用性和计算效率。
05
矩阵位移法的优缺点
优点
精确度高
矩阵位移法基于严格的数学推导,能 够精确地计算出结构的位移和内力, 尤其适用于复杂结构的分析。
适用性强
矩阵位移法可以处理多种类型的载荷 ,包括静载、动载以及温度载荷等, 适用范围广泛。
便于计算机化
矩阵位移法的计算过程可以通过计算 机程序实现,便于进行大规模的结构 分析。
多尺度方法
将矩阵位移法应用于多尺度问题 ,考虑不同尺度之间的相互作用 和影响,为复杂系统提供更准确 的模拟结果。
THANKS
感谢观看ts
目录
• 引言 • 矩阵位移法的基本概念 • 矩阵位移法的实施步骤 • 矩阵位移法的应用实例 • 矩阵位移法的优缺点 • 未来展望与研究方向
01
引言
什么是矩阵位移法
矩阵位移法是一种数值分析方法,用 于求解线性方程组和解决各种数值计 算问题。
它通过将原问题转化为矩阵形式,利 用矩阵运算来求解未知数,具有高效 、精确和灵活的特点。
并行计算
利用并行计算技术,将计算任务分解为多个子任务,同时运行在多 个处理器上,加快计算速度。
智能优化
结合人工智能和机器学习技术,自动调整算法参数,实现自适应优 化,提高算法的效率和稳定性。
矩阵位移法PPT课件
9
将一个在荷载作用下的连续结构剖分成若干个各自独立 的单元,原结构可以看成是由各单元在连接点(称结点) 连接而成的体系——化整为零
为了减少基本未知量的数目,跨 间集中荷载作用点可不作为结点, 但要计算跨间荷载的等效结点荷 载;跨间结点也可不作为结点, 但要推导相应的单元刚度矩阵, 编程序麻烦。
10
l
1 2E
l3 6EI
l2
I
6EI l2
2EI
l
1 2E I
l3
6EI l2
6EI l2
4EI
l
23
(2)只考虑轴向变形的桁架单元
由于
vie
v
e j
ie
e j
0
, 可将式
(9—5)中删去第2、3、5、6行(列),则
EA
k
e
l
M
e j
0
0
0
12EI l3 6EI l2
0
12 l
EI
3
6EI
l2
0
6EI l2 4EI l
0
6EI l2
2EI
l
EA l 0
0 EA l 0
0
0
12EI l3
6EI l2
0
12EI l3
6EI l2
0
6EI l2 2EI l
0
6EI l2
第二节 单元刚度矩阵
1. 一般单元杆端力和杆端位移的表示方法 图9-1所示平面刚架中的一等截面直杆单元e。设杆件除弯曲变形 外,还有轴向变形。杆件两端各有三个位移分量(两个移动、一个 转动),杆件共有六个杆端位移分量,这是平面杆系结构单元的一 般情况,故称为一般单元。单元的两端采用局部编码i和j。现以i 点为原点,以从i向j的方向为轴的正方向,并以轴正向逆时针转 过90°为的正方向。这样的坐标系称为单元的局部坐标系。字母、 上面的一横是局部坐标系的标志。i端、j端分别称为单元的始端
将一个在荷载作用下的连续结构剖分成若干个各自独立 的单元,原结构可以看成是由各单元在连接点(称结点) 连接而成的体系——化整为零
为了减少基本未知量的数目,跨 间集中荷载作用点可不作为结点, 但要计算跨间荷载的等效结点荷 载;跨间结点也可不作为结点, 但要推导相应的单元刚度矩阵, 编程序麻烦。
10
l
1 2E
l3 6EI
l2
I
6EI l2
2EI
l
1 2E I
l3
6EI l2
6EI l2
4EI
l
23
(2)只考虑轴向变形的桁架单元
由于
vie
v
e j
ie
e j
0
, 可将式
(9—5)中删去第2、3、5、6行(列),则
EA
k
e
l
M
e j
0
0
0
12EI l3 6EI l2
0
12 l
EI
3
6EI
l2
0
6EI l2 4EI l
0
6EI l2
2EI
l
EA l 0
0 EA l 0
0
0
12EI l3
6EI l2
0
12EI l3
6EI l2
0
6EI l2 2EI l
0
6EI l2
第二节 单元刚度矩阵
1. 一般单元杆端力和杆端位移的表示方法 图9-1所示平面刚架中的一等截面直杆单元e。设杆件除弯曲变形 外,还有轴向变形。杆件两端各有三个位移分量(两个移动、一个 转动),杆件共有六个杆端位移分量,这是平面杆系结构单元的一 般情况,故称为一般单元。单元的两端采用局部编码i和j。现以i 点为原点,以从i向j的方向为轴的正方向,并以轴正向逆时针转 过90°为的正方向。这样的坐标系称为单元的局部坐标系。字母、 上面的一横是局部坐标系的标志。i端、j端分别称为单元的始端
【实用】矩阵位移法PPT文档
局部码总码
由变于形建 连筑续工,程中刚(架和1连)续梁结构较1 多,故这里将只1介绍先处理法。
(1) 0
0
(2) 2 (3) 3 (4) 0 (5) 0 (6) 4
2
3
0
0
4
(2) 0 (3) 0 (4) 1 (5) 2 (6) 3
0
0
1
2 3
1 2
[k] 1 = 3
在进行整体分析的时候,必须要考虑支承边界条件,而这一条 件可以在形成整体刚度方程之前或之后处理,因而形成了先处理法 和后处理法两种矩阵位移法。
后处理法是先不考虑支承条件,将所有6×6的单元刚度方程一 并组集成整体刚度方程。由于还未考虑支承条件,故整体刚度方程 一定是一个奇异方程,整体刚度矩阵一定是一个奇异矩阵,在只有 引入支承边界条件后,才能消除这种奇异性,方程才可求解。后处 理法,整体刚度矩阵物理意义明确,易于修改边界条件,程序简单 ;但后处理法整体刚度矩阵较大,占用计算机内存较多,因此后处 理法对于结点多、支座约束少、必须考虑轴向变形的结构,得到广 泛应用。
先处理法是在进行整体分析前考虑支承边界条件,也就是说对 于单元刚度方程,不必把位移已知的行和对应的单元刚度矩阵的列 组集到总体刚度方程中去。这样做的好处是,最终形成的结构刚度 方程阶数小,不用再修正,即可直接求解。
先处理法特别适用于有铰结点的结构、支承结点较多、通常不 考虑轴向变形的刚架结构以及甚至连剪力都不考虑的连续梁结构的 求解。由于建筑工程中刚架和连续梁结构较多,故这里将只介绍先 处理法。实际上,两种方法由单刚组集总刚的原理是一样的,只是 后处理法待总刚生成后,再引入边界条件加以修正。
与连续梁相比: 到总体刚度方程中去。
三、单元集成过程
连续梁的矩阵位移法
第三章 连续梁的矩阵位移法
§ 3.1 概述 § 3.2 连续梁的单元刚度矩阵 § 3.3 整体刚度矩阵 § 3.4 非结点荷载的处理 § 3.5 连续梁的矩阵位移法举例
§ 3.1 概述
一、结构矩阵分析方法
结构矩阵分析方法的广泛应用是近年来结构力学最重要的
发展之一,这与计算机技术的迅速发展有直接的关系。它是以 传统的结构力学作为理论基础,以矩阵作为数学表述形式,以 电子计算机作为计算手段的三位一体的方法。
MM3223
(2)
------称为单元杆端力列阵。
(1)
1 2
(1)
(2)
2 3
(
2
)
------称为单元杆端位移列阵。
§ 3.3 整体刚度矩阵
将方程组也用矩阵表示:
4i1
2i1
0
1
M1
2i1 0
4i1 4i2 2i2
2i2 4i2
32
M2 M3
简写为: K F ------称为整体刚度方程
有非结点荷载作用时的单元杆端力,可以由两部分叠加而 得:一部分是结点受有约束、各杆件为固端梁情况下的杆端力 (固端力),另部分是综合结点荷载作用下的杆端力,即
F(e) Ff(e)k(e) (e)
§ 3.5 直接刚度法的解题步骤和算例
直接刚度法中后处理作法的解题步骤: 1.对各单元和结点进行编号 2.计算整体坐标系的单元刚度矩阵。 3.将各单元刚度矩阵的子块“对号入座”形成整体刚度矩阵。 4.计算总的荷载列阵,建立整体刚度方程。 5.引入支承条件,修改整体刚度矩阵和整体刚度方程。 6.解整体刚度方程求各结点位移。 7.计算各单元的杆端力,并进一步求各单元的其它内力。 。 8.校核。
§ 3.1 概述 § 3.2 连续梁的单元刚度矩阵 § 3.3 整体刚度矩阵 § 3.4 非结点荷载的处理 § 3.5 连续梁的矩阵位移法举例
§ 3.1 概述
一、结构矩阵分析方法
结构矩阵分析方法的广泛应用是近年来结构力学最重要的
发展之一,这与计算机技术的迅速发展有直接的关系。它是以 传统的结构力学作为理论基础,以矩阵作为数学表述形式,以 电子计算机作为计算手段的三位一体的方法。
MM3223
(2)
------称为单元杆端力列阵。
(1)
1 2
(1)
(2)
2 3
(
2
)
------称为单元杆端位移列阵。
§ 3.3 整体刚度矩阵
将方程组也用矩阵表示:
4i1
2i1
0
1
M1
2i1 0
4i1 4i2 2i2
2i2 4i2
32
M2 M3
简写为: K F ------称为整体刚度方程
有非结点荷载作用时的单元杆端力,可以由两部分叠加而 得:一部分是结点受有约束、各杆件为固端梁情况下的杆端力 (固端力),另部分是综合结点荷载作用下的杆端力,即
F(e) Ff(e)k(e) (e)
§ 3.5 直接刚度法的解题步骤和算例
直接刚度法中后处理作法的解题步骤: 1.对各单元和结点进行编号 2.计算整体坐标系的单元刚度矩阵。 3.将各单元刚度矩阵的子块“对号入座”形成整体刚度矩阵。 4.计算总的荷载列阵,建立整体刚度方程。 5.引入支承条件,修改整体刚度矩阵和整体刚度方程。 6.解整体刚度方程求各结点位移。 7.计算各单元的杆端力,并进一步求各单元的其它内力。 。 8.校核。
第三章 连续梁的矩阵位移法
§ 3.1 概述
一、结构矩阵分析方法 结构矩阵分析方法
结构矩阵分析方法的广泛应用是近年来结构力学最重要的
发展之一,这与计算机技术的迅速发展有直接的关系。 发展之一,这与计算机技术的迅速发展有直接的关系。它是以 传统的结构力学作为理论基础,以矩阵作为数学表述形式,以 传统的结构力学作为理论基础,以矩阵作为数学表述形式, 电子计算机作为计算手段的三位一体的方法。 电子计算机作为计算手段的三位一体的方法。
------为结点力 为结点力 荷载) (荷载)列阵
------称为整体刚度矩阵 称为整体刚度矩阵
结构刚度矩阵 的性质: 的性质:
1、对称性:结构刚度矩阵是一个对称矩阵,即位于主对角线 、对称性:结构刚度矩阵是一个对称矩阵, 两边对称位置的两个元素是相等的。 两边对称位置的两个元素是相等的。 2、由于连续梁结构为几何不变体系,因此其整体刚度矩阵为 、由于连续梁结构为几何不变体系, 非奇异矩阵。 非奇异矩阵。 3、结构刚度矩阵是一带状矩阵。 、结构刚度矩阵是一带状矩阵。
a)
1
q (x)
b)
Fqe1
q(x )
1 2
c)
Fqe2
Fqe1
Fqe2
2
δ1
=
δ2
=
2
+
1
1、在施加荷载之前先在结点处各加上一个刚臂用以限制结 、 点角位移,这样,单元即成为固端梁,而后施加荷载。 点角位移,这样,单元即成为固端梁,而后施加荷载。由于荷 载作用,在各杆端将产生固端剪力和固端弯矩。 载作用,在各杆端将产生固端剪力和固端弯矩。 2、在原结构的结点处分别施加与约束反力数值相等、方 、在原结构的结点处分别施加与约束反力数值相等、
结构矩阵分析方法的基本思想是:把整个结构看作是由若 基本思想是
结构力学教学课件09矩阵位移法ppt
所在行、列的副元素以及同行 的未知结点荷载改为0
满足边界条件3 0,
保持矩阵原有阶数和对称性
上节课内容概述
✓边界支承条件的处理; ✓非节点荷载的移置; ✓连续梁的矩阵分析; ✓坐标变换
静力等效原则 移到邻近结点
仅有结点荷载 作用的结构
假想约束 固定各结点
M1F, j
&
M
F 2
,
j
矩阵位移法 分析
0
0.4
69
0.625
0
0.469
1.25
(1) 各杆在局部坐标系中的 单元刚度矩阵
3.0
0
k
(2)
106
0 3.0
0
0
0 0.12 0.3
0 0.12
0.3
0 0.3 1.0 0 0.3 0.5
3.0 0 0 3.0 0 0
0 0.12 0.3
0 0.12 0.3
单元①
1
1
2 2
单元②
1
2
2 3
刚架的整体刚度矩阵,对号入座
k
k11 k21
(1) (1)
k (1) 12
k22 (1) k11 (2) k (2)
21
k12
(2)
k22
(2)
3.75 0
0 3.75 0
0
0
0
0
0
0.234 0.469
0 0.234 0.469
0.6 0.8 0 0
0
0
T
0 0
0
0
01 0
0 0
0 0 cos sin 0
0 0 sin cos 0
00 0
满足边界条件3 0,
保持矩阵原有阶数和对称性
上节课内容概述
✓边界支承条件的处理; ✓非节点荷载的移置; ✓连续梁的矩阵分析; ✓坐标变换
静力等效原则 移到邻近结点
仅有结点荷载 作用的结构
假想约束 固定各结点
M1F, j
&
M
F 2
,
j
矩阵位移法 分析
0
0.4
69
0.625
0
0.469
1.25
(1) 各杆在局部坐标系中的 单元刚度矩阵
3.0
0
k
(2)
106
0 3.0
0
0
0 0.12 0.3
0 0.12
0.3
0 0.3 1.0 0 0.3 0.5
3.0 0 0 3.0 0 0
0 0.12 0.3
0 0.12 0.3
单元①
1
1
2 2
单元②
1
2
2 3
刚架的整体刚度矩阵,对号入座
k
k11 k21
(1) (1)
k (1) 12
k22 (1) k11 (2) k (2)
21
k12
(2)
k22
(2)
3.75 0
0 3.75 0
0
0
0
0
0
0.234 0.469
0 0.234 0.469
0.6 0.8 0 0
0
0
T
0 0
0
0
01 0
0 0
0 0 cos sin 0
0 0 sin cos 0
00 0
第八章矩阵位移法-1-PPT精品
F yi
F
(e)
Fi
(e)
F
j
F3 F4
Mi F xj
F5
F6
F yj
M j
8-1 概述
32
四.坐标系选择
(e)
F
Fi
(e)
F2
F
yi
δ j 3 u j
F j F3 Fxj
4 v j
F4 Fyj
8-1 概述
29
刚架单元
局部坐标系:
(e)
δ
δ i (e) δ j
δ1 (e) δ 2 δ 3 δ 4
u i (e) vi
i
(e)
F
u j
(e)
F i F j
• 计算机计算:速度快,不怕重复。不具 有灵活性。所以处理格式最好统一,单 元类型要少。
8-1 概述
4
3.矩阵位移法是有限元法的雏形
结构矩阵分析也称为杆件结构的有限元法。
4.有限元法分析的基本思路:
离散:把整体拆开,分解成若干个单元; 集合:将单元按一定的条件集合成整体
把结构计算转化为单元的分析和集合问题。
8-1 概述
5
一.结构的离散化
单元:用杆的轴线代表。 结点:各杆轴线之间的交点。
按照自然数的顺序,对所有结点和 单元进行编号。
8-1 概述
6
例1
8-1 概述
7
例2
8-1 概述
8
例3
8-1 概述
9
结构的离散化要点
1.结点整体码:对结点以数字顺序编号, 相关结点编号的最大差值应尽可能小。
2.单元编号: 对单元按一定顺序以数字 编号,习惯上用①、②等来标记。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
总结为:“化整为零,积零为整”
§ 3.2 连续梁的单元刚度矩阵
y
M1,1
x1
① i1
M2,2
2
② i2
M3,3
3
M ,(1) (1)
M1
12 1
1
M (1) 12
① M ,(1) (1) 21 2
M2
② M ,(2) (2) 23 2
M ,(2) (2) 32 3
M3
i1
M (1) 21
2
二、结构矩阵分析方法的分类
与传统的力法、位移法和混合法对应,也有矩阵力法、矩 阵位移法和矩阵混合法。矩阵位移法具有易于实现计算过程程 序化的优点而被广泛应用,我们主要介绍矩阵位移法。
矩阵位移法又分为刚度法和直接刚度法。两者的基本原理 并无本质的区别,只是在形成所谓整体刚度矩阵时使用的方法
不同。直接刚度法比较简便得多,因此得到广泛的应用。这里 就只介绍矩阵位移法中的直接刚度法。
三、矩阵位移法的基本思路
矩阵位移法的作法同上所述:是先把结构拆散成有限数目
的杆件单元进行单元分析而后进行整体分析也就是将这些单元 再集合一起,使其满足平衡条件和位移连续条件恢复为原结构。
基本思路及过程
矩阵位移法分析问题的过程是,首先进行离散化和单元 分析,然后进行整体分析,考虑单元的集合得出基本方程组 ,通过解线性方程组求出结构的位移并求出结构的内力。
0(e)
M3 M3(22) M3 0
4i112i12M10
2i11(4i14i2)22i23M2 0(f)
2i224i23M30
即为位移法 方程
引入矩阵形式(式a、b)可写为:
M M1221 (1) 2 4ii1 1 4 2ii1 1 1 2 (1)g M M3 22 3 (2) 2 4ii2 2 4 2ii2 2 3 2 (2)g
分析过程:
1.对结构的结点和单元进行编号;
2.进行结构的离散化:将结构拆成两个杆件单元①和②;
3.进行单元分析:建立单元刚度矩阵;
4.进行整体分析:将离散化的各单元重新集合,满足原结
构的平衡条件和位移连续条件,而得到整体刚度方程。我们利 用已求得的各单元刚度矩阵形成整体刚度矩阵。形成整体刚度 矩阵的方法,以直接刚度法最为常用。
(5)根据边界条件修改结构原始刚度矩阵计算自由结点位移。 (6)计算在结构坐标系中由杆端位移产生的杆端力;再计算单 元在局部坐标系中的杆端力。
(7)计算支座反力。 (8)校核。
§ 3.4 非结点荷载的处理
以上关于矩阵位移法的讨论,是说结构的结点位移作为基 本未知量。在讨论中,我们只考虑了作用结点荷载的情况。由 此所得到的矩阵位移法基本方程,即整体刚度方程,表述了结 点位移和给点荷裁的关系。而实际上,不论是恒载还是活载常 常是作用在杆件单元上的均布荷载、分布荷载或集中荷载。对 于这种非结点荷载的处理,一种方法是,不论均布或分布荷载 都适当地改用若干集中荷载加以代替,并把集中荷载的作用点 也看作结点。这样处理的结果是,加多了单元和结点位移,从 而增加了计算工作量。另一种则是目前通用的处理方法,即采
4i1 2i1
0
K 2i1
4i1 4i2
2i2
0
2i2
4i2
------称为整体刚度矩阵
1 2
3
------为结点 位移列阵
F
M
M
1
2
M 3
------为结点力 (荷载)列阵
结构刚度矩阵 的性质:
1、对称性:结构刚度矩阵是一个对称矩阵,即位于主对角线 两边对称位置的两个元素是相等的。
2
2 (c)
(2) 3
3
M M M
(1) 12
(1) 21
(2) 23
4i1 1 2i1 2 2i1 1 4i1 2 4i2 2 2i2
3
( d
)
M
(2) 32
2i2 2
4
i
2
3
由结点平衡条件:
再将(d)式代入,得:
M1 M2
M1(12) M2(11)
M1 0 M2(23) M2
M (2) 23
i2
M (2) 32
3
单元 ①:
M(1) 12
M(1) 21
4i11(1) 2i11(1)
42ii1122((11))(a)
由位移连续条件得:
单元 ②:
M(2) 23
M(2) 32
4i22(2) 2i22(2)
42ii2233((22))(b)
(1) 1
(1) 2
1 (2)
连续梁的矩阵位移法优秀课件
§ 3.1 概述
一、结构矩阵分析方法
结构矩阵分析方法的广泛应用是近年来结构力学最重要的
发展之一,这与计算机技术的迅速发展有直接的关系。它是以 传统的结构力学作为理论基础,以矩阵作为数学表述形式,以 电子计算机作为计算手段的三位一体的方法。
结构矩阵分析方法的基本思想是:把整个结构看作是由若
简写为: F(1) k(1)δ(1) (h) F(2) k(2)δ(2) (i)
------称为单元刚度方程
其中:
k (1)
4i1 2i1
2i1
4i1
j
k (2)
4i2 2i2
பைடு நூலகம்
------称为单元刚度矩阵。
2i2
4i2
k
矩阵中的各元素称为单元刚度影响系数。
F (1)
MM1221
(1)
F (2)
干单个杆件(称为单元)所组成的集合体。
1)单元分析:在进行分析时,首先把结构拆散成有限数目
的杆件单元(结构的离散化),写出各单元杆端的力与位移两 者的关系式;
2)整体分析:即将这些单元再集合一起,使其满足平衡
条件和位移连续条件也就是保证离散化了的杆件单元重新集合 后仍恢复为原结构;
3)解方程组:求出结构的结点位移和内力。
2、由于连续梁结构为几何不变体系,因此其整体刚度矩阵为 非奇异矩阵。
3、结构刚度矩阵是一带状矩阵。
*
\
*
\
0
*
*
*
\
\
\
*
*
*
\
\
\
*
*
*
\
\
\
*
*
*
0
\
\
*
*
综上所述,可将直接刚度法的解算步骤归纳如下:
(1)将结点和单元进行编号;选择结构坐标系和局部坐标系。 (2)把所有结点力沿结构坐标系分解;建立结点位移列向量和 结点力列向量(两者的分量要一一对应)。 (3) 计算结构坐标系中各单元刚度矩阵的四个子块。 (4)将各单元刚度矩阵的四个子块,按其两个下标在结构原始 刚度矩阵中“对号入座”。
MM3223
(2)
------称为单元杆端力列阵。
(1)
1 2
(1)
(2)
2 3
(
2
)
------称为单元杆端位移列阵。
§ 3.3 整体刚度矩阵
将方程组也用矩阵表示:
4i1
2i1
0
1
M1
2i1 0
4i1 4i2 2i2
2i2 4i2
32
M2 M3
简写为: K F ------称为整体刚度方程
§ 3.2 连续梁的单元刚度矩阵
y
M1,1
x1
① i1
M2,2
2
② i2
M3,3
3
M ,(1) (1)
M1
12 1
1
M (1) 12
① M ,(1) (1) 21 2
M2
② M ,(2) (2) 23 2
M ,(2) (2) 32 3
M3
i1
M (1) 21
2
二、结构矩阵分析方法的分类
与传统的力法、位移法和混合法对应,也有矩阵力法、矩 阵位移法和矩阵混合法。矩阵位移法具有易于实现计算过程程 序化的优点而被广泛应用,我们主要介绍矩阵位移法。
矩阵位移法又分为刚度法和直接刚度法。两者的基本原理 并无本质的区别,只是在形成所谓整体刚度矩阵时使用的方法
不同。直接刚度法比较简便得多,因此得到广泛的应用。这里 就只介绍矩阵位移法中的直接刚度法。
三、矩阵位移法的基本思路
矩阵位移法的作法同上所述:是先把结构拆散成有限数目
的杆件单元进行单元分析而后进行整体分析也就是将这些单元 再集合一起,使其满足平衡条件和位移连续条件恢复为原结构。
基本思路及过程
矩阵位移法分析问题的过程是,首先进行离散化和单元 分析,然后进行整体分析,考虑单元的集合得出基本方程组 ,通过解线性方程组求出结构的位移并求出结构的内力。
0(e)
M3 M3(22) M3 0
4i112i12M10
2i11(4i14i2)22i23M2 0(f)
2i224i23M30
即为位移法 方程
引入矩阵形式(式a、b)可写为:
M M1221 (1) 2 4ii1 1 4 2ii1 1 1 2 (1)g M M3 22 3 (2) 2 4ii2 2 4 2ii2 2 3 2 (2)g
分析过程:
1.对结构的结点和单元进行编号;
2.进行结构的离散化:将结构拆成两个杆件单元①和②;
3.进行单元分析:建立单元刚度矩阵;
4.进行整体分析:将离散化的各单元重新集合,满足原结
构的平衡条件和位移连续条件,而得到整体刚度方程。我们利 用已求得的各单元刚度矩阵形成整体刚度矩阵。形成整体刚度 矩阵的方法,以直接刚度法最为常用。
(5)根据边界条件修改结构原始刚度矩阵计算自由结点位移。 (6)计算在结构坐标系中由杆端位移产生的杆端力;再计算单 元在局部坐标系中的杆端力。
(7)计算支座反力。 (8)校核。
§ 3.4 非结点荷载的处理
以上关于矩阵位移法的讨论,是说结构的结点位移作为基 本未知量。在讨论中,我们只考虑了作用结点荷载的情况。由 此所得到的矩阵位移法基本方程,即整体刚度方程,表述了结 点位移和给点荷裁的关系。而实际上,不论是恒载还是活载常 常是作用在杆件单元上的均布荷载、分布荷载或集中荷载。对 于这种非结点荷载的处理,一种方法是,不论均布或分布荷载 都适当地改用若干集中荷载加以代替,并把集中荷载的作用点 也看作结点。这样处理的结果是,加多了单元和结点位移,从 而增加了计算工作量。另一种则是目前通用的处理方法,即采
4i1 2i1
0
K 2i1
4i1 4i2
2i2
0
2i2
4i2
------称为整体刚度矩阵
1 2
3
------为结点 位移列阵
F
M
M
1
2
M 3
------为结点力 (荷载)列阵
结构刚度矩阵 的性质:
1、对称性:结构刚度矩阵是一个对称矩阵,即位于主对角线 两边对称位置的两个元素是相等的。
2
2 (c)
(2) 3
3
M M M
(1) 12
(1) 21
(2) 23
4i1 1 2i1 2 2i1 1 4i1 2 4i2 2 2i2
3
( d
)
M
(2) 32
2i2 2
4
i
2
3
由结点平衡条件:
再将(d)式代入,得:
M1 M2
M1(12) M2(11)
M1 0 M2(23) M2
M (2) 23
i2
M (2) 32
3
单元 ①:
M(1) 12
M(1) 21
4i11(1) 2i11(1)
42ii1122((11))(a)
由位移连续条件得:
单元 ②:
M(2) 23
M(2) 32
4i22(2) 2i22(2)
42ii2233((22))(b)
(1) 1
(1) 2
1 (2)
连续梁的矩阵位移法优秀课件
§ 3.1 概述
一、结构矩阵分析方法
结构矩阵分析方法的广泛应用是近年来结构力学最重要的
发展之一,这与计算机技术的迅速发展有直接的关系。它是以 传统的结构力学作为理论基础,以矩阵作为数学表述形式,以 电子计算机作为计算手段的三位一体的方法。
结构矩阵分析方法的基本思想是:把整个结构看作是由若
简写为: F(1) k(1)δ(1) (h) F(2) k(2)δ(2) (i)
------称为单元刚度方程
其中:
k (1)
4i1 2i1
2i1
4i1
j
k (2)
4i2 2i2
பைடு நூலகம்
------称为单元刚度矩阵。
2i2
4i2
k
矩阵中的各元素称为单元刚度影响系数。
F (1)
MM1221
(1)
F (2)
干单个杆件(称为单元)所组成的集合体。
1)单元分析:在进行分析时,首先把结构拆散成有限数目
的杆件单元(结构的离散化),写出各单元杆端的力与位移两 者的关系式;
2)整体分析:即将这些单元再集合一起,使其满足平衡
条件和位移连续条件也就是保证离散化了的杆件单元重新集合 后仍恢复为原结构;
3)解方程组:求出结构的结点位移和内力。
2、由于连续梁结构为几何不变体系,因此其整体刚度矩阵为 非奇异矩阵。
3、结构刚度矩阵是一带状矩阵。
*
\
*
\
0
*
*
*
\
\
\
*
*
*
\
\
\
*
*
*
\
\
\
*
*
*
0
\
\
*
*
综上所述,可将直接刚度法的解算步骤归纳如下:
(1)将结点和单元进行编号;选择结构坐标系和局部坐标系。 (2)把所有结点力沿结构坐标系分解;建立结点位移列向量和 结点力列向量(两者的分量要一一对应)。 (3) 计算结构坐标系中各单元刚度矩阵的四个子块。 (4)将各单元刚度矩阵的四个子块,按其两个下标在结构原始 刚度矩阵中“对号入座”。
MM3223
(2)
------称为单元杆端力列阵。
(1)
1 2
(1)
(2)
2 3
(
2
)
------称为单元杆端位移列阵。
§ 3.3 整体刚度矩阵
将方程组也用矩阵表示:
4i1
2i1
0
1
M1
2i1 0
4i1 4i2 2i2
2i2 4i2
32
M2 M3
简写为: K F ------称为整体刚度方程