矩阵位移法等效结点荷载

合集下载

结构力学二7-矩阵位移法

结构力学二7-矩阵位移法
简记为
e e
e
1
4ie
1
e
1e 2ie 2e 4ie
F k
2ie
---单元刚度方程 其中
k e称作单元刚度矩阵(简称作单刚)
1
ie
e 1
e
F2e
单元刚度矩阵中元素的物理意义
e e 4ie 2ie k k e 11 12 k e e 2 i 4 i k k e 21 22 e
F 4i 2ie
e 1 e e 1
e 2
ie
e 1
e
F2e
F
e 1

e 2
2
e F2e 2ie1e 4ie 2
1e
F
e
2e
F2e
e 1
F1 4ie 2ie 1 2 i 4 i F e 2 2 e
1
1/2
2
M-图(kN· m)
(2)乘大数法 若 i 0 ,则将总刚主对角 元素 kii 乘以大数N.
6kN.m
3kN.m
i1 1 i2 2
2 3
P3
1
4 2 0 1 6 2 12 4 3 2 P 0 4 8 3 3
4 2 0 1 / 2 1 0 F 2 4 1 / 4 1 2 1 / 4 0 8 4 1 / 4 2 3 2 F 4 8 0 1 1 1
q
练习: 求图示结构的等效结点荷载. q
1 2 3 4
1
2

结构力学应用-矩阵位移法

结构力学应用-矩阵位移法

3、集成总刚
(6)定位向量法:对号入座,同号相加 定位向量法:对号入座,
4.综合结点荷载
综合结点荷载 {F}={FD}+{FE} }――直接结点荷载 ①{FD}――直接结点荷载 }――等效结点荷载 ②{FE}――等效结点荷载 (7-1)局部坐标系单元固端力 (7-2)整体坐标系单元固端力 (7-3)单元等效结点荷载。 单元等效结点荷载。
等效原则: 等效原则: ——两种荷载对基本体系产生相同的结点位移。 两种荷载对基本体系产生相同的结点位移 ——两种荷载对基本体系产生相同的结点位移。
矩阵位移法的计算步骤及示例
矩阵位移法计算平面刚架 计算机计算――程序化) 程序化) (计算机计算 程序化
1. 编码、整理原始数据 编码、
(1)整体与局部坐标系 ) (2)结点位移编码 ) 单元编码 (3)原始数据: )原始数据: E 、A i、I i、l i 定位向量{λ} 定位向量 e, αi([ T ]) ])
几点补充说明
1、结点位移分量编号,定位向量 、结点位移分量编号,
——引入支承条件:已知位移约束的方向,编码为零。 引入支承条件:已知位移约束的方向,编码为零。 引入支承条件
2、铰结点处理: 铰结点处理: 铰结点处理
铰结的各杆杆端的转角均为基本未知量 ——分别编码(统一单元,程序简单) 分别编码(统一单元,程序简单) 分别编码
矩阵位移法
矩阵位移法——基本原理与位移法相同 基本原理与位移法相同 矩阵位移法 *数学工具 —— 矩阵运算
1、矩阵知识 矩阵: (1)矩阵:A 方阵: 方阵: 阶方阵A相应的行列式 (2)行列式:n阶方阵 相应的行列式 )行列式: 阶方阵 相应的行列式D 若D=0,A为奇异矩阵 (3)矩阵运算 相等:加减:数乘: 相等:加减:数乘: l aik 乘法: 乘法:Cmn=Aml*Bln,则 cij =

李廉锟《结构力学》(上册)配套题库【课后习题】(矩阵位移法)【圣才出品】

李廉锟《结构力学》(上册)配套题库【课后习题】(矩阵位移法)【圣才出品】

第10章矩阵位移法复习思考题1.矩阵位移法的基本思路是什么?答:矩阵位移法的基本思路:(1)单元分析单元分析是指将结构先分解为有限个较小的单元,即离散化,在较小的范围内分析单元的内力与位移之间的关系,建立单元刚度矩阵或单元柔度矩阵。

(2)整体分析整体分析将将单元分析中的各单元集合成原来的结构,要求各单元满足原结构的几何条件(包括支承条件、结点处的变形连续条件)和平衡条件,建立整个结构的刚度方程或柔度方程,以求解原结构的内力和位移。

(3)支承条件引入支承条件,修改结构原始刚度方程。

(4)求解解算结构刚度方程,求出结点位移,计算各单元杆端力。

2.试述矩阵位移法与传统位移法的异同。

答:矩阵位移法与传统位移法的异同点:(1)相同点传统位移法的基本原理,是以在小变形的基础的结构体系中,内力是可以叠加的,位移也是可以叠加的,而矩阵位移法是按传统位移法的基本原理运用矩阵计算内力和位移的方法。

因此矩阵位移法和传统位移法的基本原理在实质上是一致的。

(2)不同点①矩阵位移法中一般考虑杆件轴向变形的影响,传统位移法忽略杆件的轴向变形;②矩阵位移法一般在计算机上进行计算,可以解决大型复杂问题;传统位移法的计算手段一般是手算,只用来解决简单问题。

3.矩阵位移法中,杆端力、杆端位移和结点力、结点位移的正负号是如何规定的?答:杆端力沿局部坐标系的、的正方向为正,杆端弯矩逆时针为正;杆端位移的正负同杆端力和弯矩。

结点力沿整体坐标系x、y的正方向为正,结点力偶逆时针为正;结点位移的正负同结点力和力偶。

4.为何用矩阵位移法分析时,要建立两种坐标系?答:因为单元刚度矩阵是建立在杆件的局部坐标系上的,但对于整体结构,各单元的局部坐标系可能不尽相同,在研究结构的几何条件和平衡条件时,需要选定一个统一的坐标系即为整体坐标系,另外按局部坐标系建立的单元刚度矩阵可以通过坐标转换到整体坐标系中,从而得到整体坐标系中的单元刚度矩阵。

故建立两种坐标系使矩阵位移法的思路更清晰,物理意义更明确,且不会影响计算结果。

第九节矩阵位移法

第九节矩阵位移法

(2 =1)
0
6EI l2 2EI l
0
6EI
l2 4EI
l
e
…(9-4)
F e k ee
…(9-5)
即为一般单元的刚度方程。其中 k e 称为局部坐标系中的单
元刚度矩阵。
2、一般单元刚度矩阵的性质
(1)单元刚度系数的意义
单元刚度矩阵中的每个元素称为单元刚度系数 kij ,其物理
意义表示由于单位杆端位移引起的杆端力。
( v1e
v2e
)
Fye1
6EI l2
(1e
2e )
12EI l3
( v1e
v2e )
Fye2
6EI l2
(1e
e 2
)
12 l
EI
3
( v1e
v2e
)
Fx1 M1
1
v1
Fy1
u1
…(9-2)
e
1
M2
Fx2
2 Fy2
v2
u2
2
式(9-1) 、(9-2)即为局部坐标系下平面刚架一般单元的单元刚度方
ke T Tk eT
F e kee
即为单元e在整体坐标中的单元刚度方程 其中 k e为整体坐标系的单元刚度矩阵,和 k e 同阶,且具有类似的性质。
§9-4 结构的整体刚度矩阵
作用在结构上的荷载与结构的结点位移, 也存在一一对应的关系,即为结构的整体刚 度方程。结构的整体刚度方程反映了结点荷 载和结构位移之间的关系,其实质就是位移 法的基本方程。求解方法一种是传统位移法, 另一种是直接刚度法。
l
Fxe1
EA l
u1e
EA l
u2e

结构力学龙第六版课后答案详解第八章

结构力学龙第六版课后答案详解第八章

结构力学龙第六版课后答案详解第八章
1.在矩阵位移法中,结构在等效结点荷载作用下的内力,与结构在原有荷载作用下的内力相同。

()
【答案】错
【解析】结构在等效结点荷载与在原有荷载作用下的内力不相同,还可能有固端力作用。

2.当荷载作用在静定多跨梁的基本部分时,对附属部分是否引起内力?当荷载作用在静定多跨梁附属部分时,对基本部分是否引起内力?
答案:(1)静定多跨梁的基本部分能独立承受荷载并维持平衡;因此,当荷载只作用于基本部分时,对附属部分不引起内力。

(2)静定多跨梁的附属部分依赖于基本部分才能维持平衡。

当荷载作用在附属部分时,附属部分的支座反力将通过连接较作用于基本部分;因此,当荷载作用在静定多跨梁附属部分时,对基本部分将引起内力。

矩阵位移法

矩阵位移法

k22坐k11标局k01成部1k029坐200标时kk20与32,3 整局k0体12部45 单k0k20514
0 k26 k26
To 47
k e ke
刚和有何整k关体3k3系单33 ?刚k0k间454535
k35 00
k3k6 36
0 k56
对称对称
kk5544
kk65k66 66
F e FEe k e e
单元杆端位移矩阵
e 1
2
3
4
T e
单元刚度矩阵(应熟记)
12 6l 12 6l
k
e
EI l3
6l
12
4l 2 6l
6l 12
2l
2
6l
6l 2l 2 6l 4l 2
是转角位移方程的矩阵表示
单元等效结点荷载矩阵
根据单跨梁的载常数,可得
向上满跨均布荷载 q 作用
(F FE )e k e e F e FEe k e e
连续梁单元需要 进行坐标转换吗?
连续梁的局部坐标与整 体坐标一致,所以不需 要转换。
第一种做法
桁架单元如何
进行坐标转换? T
力的转换
T
F1
F2
F3
F4
T
cos
0
位移的转换
sin
0
0
cos
0 T F1
sin F2
1 2
3. 坐标转换问题
在搞清单元特性后,像位移法一样,需将单 元拼装回去。在结点处位移自动满足协调条件 的基础上,令全部结点平衡,即可建立求解位 移的方程,这是下一节将讨论的内容。
除连续梁外,一般结构单元不全同方位, 为保证协调和平衡,应将杆端位移和杆端力 都转换成统一的,对整体坐标的量,因此要 先解决坐标转换问题。下面先讨论自由式梁 单元的转换问题。

《结构力学习题集》(下)-矩阵位移法习题及答案 (2)

《结构力学习题集》(下)-矩阵位移法习题及答案 (2)

第七章 矩阵位移法一、就是非题1、单元刚度矩阵反映了该单元杆端位移与杆端力之间得关系。

2、单元刚度矩阵均具有对称性与奇异性。

3、局部坐标系与整体坐标系之间得坐标变换矩阵T 就是正交矩阵。

4、结构刚度矩阵反映了结构结点位移与荷载之间得关系。

5、用 矩 阵 位 移 法 计 算 连 续 梁 时 无 需 对 单 元 刚 度 矩 阵 作 坐 标 变 换。

6、结 构 刚 度 矩 阵 就是 对 称 矩 阵 ,即 有K i j = K j i ,这 可 由 位 移 互 等 定 理 得 到 证 明 。

7、结构刚度方程矩阵形式为:,它就是整个结构所应满足得变形条件。

8、在直接刚度法得先处理法中,定位向量得物理意义就是变形连续条件与位移边界条件。

9、等效结点荷载数值等于汇交于该结点所有固端力得代数与。

10、矩阵位移法中,等效结点荷载得“等效原则”就是指与非结点荷载得结点位移相等。

11、矩阵位移法既能计算超静定结构,也能计算静定结构。

二、选择题1、已知图示刚架各杆EI = 常数,当只考虑弯曲变形,且各杆单元类型相同时,采用先处理法进行结点位移编号,其正确编号就是:(0,1,2)(0,0,0)(0,0,0)(0,1,3)(0,0,0)(1,2,0)(0,0,0)(0,0,3)(1,0,2)(0,0,0)(0,0,0)(1,0,3)(0,0,0)(0,1,2)(0,0,0)(0,3,4)A.B.C.D.21341234123412342、平面杆件结构一般情况下得单元刚度矩阵,就其性质而言,就是:A.非对称、奇异矩阵;B.对称、奇异矩阵;C.对称、非奇异矩阵;D.非对称、非奇异矩阵。

3、单元i j 在图示两种坐标系中得刚度矩阵相比:A.完全相同;B.第2、3、5、6行(列)等值异号;C.第2、5行(列)等值异号;D.第3、6行(列)等值异号。

4、矩阵位移法中,结构得原始刚度方程就是表示下列两组量值之间得相互关系:A.杆端力与结点位移;B.杆端力与结点力;C.结点力与结点位移;D.结点位移与杆端力。

矩阵位移法程序化解题方法

矩阵位移法程序化解题方法

矩阵位移法解法步骤解:1)、单元及结点位移分量统一编码单元及结点位移分量编码、整体坐标系如图所示,局部坐标系横轴正向在各单元上标出。

注:编结点位移分量总码时,后处理法和先处理法有区别:采用后处理法编码时暂不考虑边界条件对支座处位移分量的限制,皆视为一般情形处理;采用先处理法时,对已知为零的位移分量总是以零编码。

对于连接于铰结点的杆端编码时,线位移采用同码,而角位移异码。

2)、形成局部坐标中单元刚度矩阵 k e:首先,计算各单元杆件的几何特征:⋯ ⋯各单元的单元刚度矩阵如下:单元①: ⋯ ⋯3)、形成整体坐标中单元刚度矩阵:(计算公式: k e = T T ke T ) 整体坐标系中的各单元刚度矩阵转换如下:单元①: ⋯ ⋯4)、集成整体刚度矩阵 K (单元集成法或直接刚度法):首先,由各单元的局部码与总码的对应关系写出各单元的定位向量如下:λ e = ⋯ ⋯ T其次,将各单元刚度矩阵 k e 按其定位向量 λ e 在整体刚度矩阵 K 中定位并累加 得整体刚度矩阵如下:K =(⋯ ⋯)5)、计算综合等效结点荷载向量 F P :①、计算局部坐标系中各杆件单元的固端力向量:F P e =(F N1F ,F Q1F ,M 1F ,F N2F ,F Q2F ,M 2F )T ②、转换整体坐标系中各杆件单元的固端力向量:{F P }e =(F x1F ,F y1F ,M 1F ,F x2F ,F y2F ,M 2F )T ③、将各杆件单元的固端力反其指向,并按其定位向量 λ e 在综合等效结点荷载向量 F P 定位并累加,得综合等效结点荷载向量如下:F P = ⋯ ⋯ T6)、计入边界条件条件,写出刚度方程并解之:刚度方程: K Δ = F P采用后处理法时,对已知为零的结点位移,在整体刚度矩阵 K 中将其所对应行列的主元素记为1,其余都变为零,然后写出刚度方程,解之。

采用先处理法时,由于在进行位移分量编码时已考虑边界条件,因而无须再计入,只写出刚度方程求解即可。

矩阵位移法等效结点荷载

矩阵位移法等效结点荷载
将各单元的杆端力反号转移到结点上把结点相关单元同一方向的力相加后以整体坐标系为标准判断正负号一致为正反之为负即得到该结点等效结点荷载列根据整体位移码将等效结点荷载列阵与直接结点荷载列阵元素对应相加即得到结构的综合结点荷载列阵
第九章 矩阵位移法
9.4 等效结点荷载
非结点荷载的处理
有限元分析的重要一步是把一个连续的结构看成是 由离散单元在结点处连接拼装而成。 由离散单元在结点处连接拼装而成。而把作用在结构上 的荷载统统考虑作用在结点上。 的荷载统统考虑作用在结点上。然而无论是恒载还是活 常常是分布作用在单元上。 载,常常是分布作用在单元上。 对这种非结点荷载的处理: 对这种非结点荷载的处理: 方法一:把分布荷载改用若干集中荷载代替, 方法一:把分布荷载改用若干集中荷载代替,并把集 中荷载的作用点选作结点; 中荷载的作用点选作结点; 方法二:等效结点荷载法。 方法二:等效结点荷载法。
FP l FP l 2 − FP l 2 8
q

1(1,2,3)


ql 2 12
ql 2
没有结 点位移
0
+
0
ql 2 ql − 2 −
F l FP l − P F l − ql − P 2 2 FP l 12 8 8 ql 0 只有结 2 ql 点位移 2 ql 2
2
ql − 2

12
0
ql 2 − 12
(1)
FP 2(4,5,6) ② 3(7,8,9) q ①
1(1,2,3)
{F }
F e
F F F = F F F
F 1 F 2 F 3 F 4 F 5 F 6

e
{F }

矩阵位移法(载荷与边界)

矩阵位移法(载荷与边界)
e 1 e 2 T
---单元等效结点荷载
e P t 0 FP 1 e 0 t P FP 2 e T t FP 1 e T t FP 2
T
P
q
2 1 0 2 1 0 3
T
3 1 0 3 0 0 0
6 1 6 1 6 48 6 24 1 k 1 6 1 6 6 24 6 48
T
3(1,0,2)
4(1,0,3)
2
1
1(0,0,0)
3
2(0,0,0)
Fp2 x

0,0,0 4
5
0,0,0
y
2 FP 2 5 0 6 3 0 7
三、非结点载荷的计算
e
FP 4
FP 6 FP 3
FP1
FP 5
FP 2
FP1 FP 2 e FP FP 6 ---单元固端力
e
e
Pe T e FP
1
2
6 4 8 48
24
3
K
1 6
1 24 2 48 3
P Pj FP
10 Pj 0 0
FP14
q
3(1,0,2)
4(1,0,3)
2 1
1(0,0,0)
3
2(0,0,0)
F
1 P3
F
1 P2
FP11
30 1 60 2 1 FP 3 30 60 4
解方程
2 6 6
6 96 24

【实用】矩阵位移法PPT文档

【实用】矩阵位移法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的单元刚度方程一 并组集成整体刚度方程。由于还未考虑支承条件,故整体刚度方程 一定是一个奇异方程,整体刚度矩阵一定是一个奇异矩阵,在只有 引入支承边界条件后,才能消除这种奇异性,方程才可求解。后处 理法,整体刚度矩阵物理意义明确,易于修改边界条件,程序简单 ;但后处理法整体刚度矩阵较大,占用计算机内存较多,因此后处 理法对于结点多、支座约束少、必须考虑轴向变形的结构,得到广 泛应用。
先处理法是在进行整体分析前考虑支承边界条件,也就是说对 于单元刚度方程,不必把位移已知的行和对应的单元刚度矩阵的列 组集到总体刚度方程中去。这样做的好处是,最终形成的结构刚度 方程阶数小,不用再修正,即可直接求解。
先处理法特别适用于有铰结点的结构、支承结点较多、通常不 考虑轴向变形的刚架结构以及甚至连剪力都不考虑的连续梁结构的 求解。由于建筑工程中刚架和连续梁结构较多,故这里将只介绍先 处理法。实际上,两种方法由单刚组集总刚的原理是一样的,只是 后处理法待总刚生成后,再引入边界条件加以修正。
与连续梁相比: 到总体刚度方程中去。
三、单元集成过程

矩阵位移法的计算步骤及示例

矩阵位移法的计算步骤及示例

=
EA l
⎡0.72855 ⎢⎣0.57006
0.57006⎤ 1 2.47855⎥⎦ 2
(5)解算结构刚度方程
28
解算结构刚度方程,求出结点位移
EA ⎡0.72855 l ⎢⎣0.57006
0.57006⎤ 2.47855⎥⎦
⎩⎨⎧uv11
⎫ ⎬ ⎭
=
⎨⎧FP ⎩0
⎫ ⎬ ⎭
Δ1
=
⎩⎨⎧uv11
k(3)
=
⎢ ⎢ ⎢
l(3) 2EI
⎢⎣ l ( 3 )
4
2EI l(3) 4EI l(3)
⎤ ⎥ ⎥ ⎥ ⎥⎦
3 4
(3)集成结构刚度矩阵K
12
由各单元刚度矩阵的上方和右侧的单元定位 向量,集成结构刚度矩阵K,此时结构刚度 矩阵K 为4阶方阵。
1
2
3
⎡ 4EI
2 EI
⎢ ⎢ ⎢
l 2
(1)
EI
k (1) l ( 1 )
4EI + 4EI
0 2 EI
K
=
⎢ l(1) ⎢
⎢0

l(1) l( 2 EI
2
)
k
(2)
4
l EI
(2
+
)
4 EI
l(2)
l(2) l(3)
⎢ ⎢⎣
0
0
2EI k (3)
l(3)
4
0
⎤ ⎥
1

0 ⎥2 ⎥
2 EI l(3)
⎥ ⎥ ⎥
3
4 EI l(3)
⎥ ⎥⎦
4
13
将各杆所需有关数据计算如下:

结构力学——矩阵位移法

结构力学——矩阵位移法

局部坐标系
整体坐标系
18

第三节
单元分析(整体坐标系下的单元分析 )
1、单元坐标转换矩阵 两种坐标系中单元的杆端位移转换关系为:
y
y
u1
u1
1
v1
e
v2 2
u2
x
u2
局部坐标系下的分量
x
e
整体坐标系下的分量
u 1 e v u co s s i n 0 0 1 1 u u 0 0 co s s i n 2 2 v 2 用整体量表示局部量 ?
矩阵位移法是有限元法的雏形,因此结构矩阵分析有时也称为 杆件结构的有限元法。在本章中将使用有限元法中的一些术语 和提法。
6

第一节
矩阵位移法概述
1、矩阵位移法的基本思路 a、方法的选择
位移法与力法之由于选取的基本未知量不同,因此计算次序不同

结构结点力 杆件杆端力

杆件端点位移 结构结点位移
位移法
正问题
力学 模型 解的 性质
e
F
e
反问题 F e

e
将单元视为两端有人为 约束控制的杆件。 e 控制附加约束加以指 定。
将单元视为两端自由的杆 F e 件, 直接加在自由端作 为指定的杆端力。
F 都 为任何值时, 有对应的唯一解,且总 是平衡力系。

第一节
矩阵位移法概述
结构力学传统方法与结构矩阵分析方法,二者同源而有别:
前者在“手算”的年代形成,后者则着眼于“电算”,计算手 段的不同,引起计算方法的差异。
在原理上同源,在作法上有别
与传统的力法、位移法相对应,在结构矩阵分析中也有矩阵力 法和矩阵位移法,或称柔度法与刚度法。矩阵位移法由于具有 易于实现计算过程程序化的优点而广为流传。

《结构力学习题集》(下)-矩阵位移法习题及答案 (2)

《结构力学习题集》(下)-矩阵位移法习题及答案 (2)

第七章 矩阵位移法一、是非题1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。

2、单元刚度矩阵均具有对称性和奇异性。

3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。

4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。

5、用 矩 阵 位 移 法 计 算 连 续 梁 时 无 需 对 单 元 刚 度 矩 阵 作 坐 标 变 换。

6、结 构 刚 度 矩 阵 是 对 称 矩 阵 ,即 有K i j = K j i ,这 可 由 位 移 互 等 定 理 得 到 证 明 。

7、结构刚度方程矩阵形式为:[]{}{}K P ∆=,它是整个结构所应满足的变形条件。

8、在直接刚度法的先处理法中,定位向量的物理意义是变形连续条件和位移边界条件。

9、等效结点荷载数值等于汇交于该结点所有固端力的代数和。

10、矩阵位移法中,等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。

11、矩阵位移法既能计算超静定结构,也能计算静定结构。

二、选择题1、已知图示刚架各杆EI = 常数,当只考虑弯曲变形,且各杆单元类型相同时,采用先处理法进行结点位移编号,其正确编号是:2、平面杆件结构一般情况下的单元刚度矩阵[]k 66⨯,就其性质而言,是: A .非对称、奇异矩阵; B .对称、奇异矩阵; C .对称、非奇异矩阵; D .非对称、非奇异矩阵。

3、单元i j 在图示两种坐标系中的刚度矩阵相比:A .完全相同;B .第2、3、5、6行(列)等值异号;C .第2、5行(列)等值异号;D .第3、6行(列)等值异号。

4、矩阵位移法中,结构的原始刚度方程是表示下列两组量值之间的相互关系: A .杆端力与结点位移; B .杆端力与结点力; C .结点力与结点位移; D .结点位移与杆端力 。

5、单 元 刚 度 矩 阵 中 元 素 k ij 的 物 理 意 义 是 :A .当 且 仅 当 δi =1 时 引 起 的 与 δj 相 应 的 杆 端 力 ;B .当 且 仅 当 δj =1时 引 起 的 与 δi 相 应 的 杆 端 力 ;C .当 δj =1时 引 起 的 δi 相 应 的 杆 端 力 ;D .当 δi =1时 引 起 的 与 δj 相 应 的 杆 端 力。

第10章矩阵位移法

第10章矩阵位移法

整体坐标系(结构坐标系):整个结构统一的坐标系。 整体坐标系(结构坐标系):整个结构统一的坐标系。 ):整个结构统一的坐标系
e e e FNi = Fxi cosα + Fyi sin α e e e F i = Fxi sin α + Fyi cosα S e e e FNj = Fxj cosα + Fyj sin α e e F e = Fxj sin α + Fyi cosα Sj
e FNi =
由图a、 , 由图 、d,根据叠加原理可写出
EA e EA e ui uj l l EA e EA e e uj FNj = ui + l l
§10-2 单元刚度矩阵
可写出
12EI e 6EI e 12EI e 6EI e 12EI 6EI 12EI 6EI vi + 2 i 3 v j + 2 j F e = 3 vie 2 ie + 3 vje 2 je Sj l3 l l l l l l l 6EI 4EI e 6EI e 2EI e 6EI 2EI e 6EI e 4EI e Mie = 2 vie + i 2 vj + j M e = 2 vie + i 2 vj + j j l l l l l l l l Fe = Si
F1 F F = 2 F3 F4
Fx1 式中 F1 = Fy1 M 1
Fx 2 F2 = Fy 2 M 2
Fx 3 F3 = Fy 3 M 3
Fx 4 F4 = Fy 4 M 4
结点2、 处 结点外力F 是给定的结点荷载; 结点 、3处:结点外力 2、F3是给定的结点荷载; 支座1、 处 结点外力F 是支座反力, 支座 、4处:结点外力 1、F4是支座反力,如支座有给定结点荷 为结点荷载与支座反力的代数和。 载,则F1、F4为结点荷载与支座反力的代数和。

矩阵位移法

矩阵位移法

原理同源---
(1)以结点位移为基本未知量,
(2)以单元分析为基础(力法计算的 结果单元刚度方程);
(3) 建立平衡方程求出结点位移,
(4) 将结点位移代入单元刚度方 程求得内力
矩 阵 位 移 法
作法有别-(1)矩阵组织数据,矩阵运算;
(2)设计计算机程序(正确);
(3) 原始数据的准备、输入、计算 结果的输出及正确性判别等 特点: 省力;计算速度快;计算结果精度高 ;使用者要力学概念清楚。
1 0 0 1 0 0 8 2 4 i i 2 42 i i 0 2 4 3
修改后的位移 法方程
(6) 解方程
矩 阵 位 移 法
0 1 3.571 2 i 3 12.286 i
(5)引入支承条件修改原始刚度方程
矩 阵 位 移 法
K FP
4 i 2i 0 1 4 2 i 8 i 2i 4 2 i 0 2 4 42 i 3
主1副0法修改后 原始刚度方程
整 体 刚 度 方 程
单元刚度集成法
矩 阵 位 移 法
单元(1)对号 入座
单元刚度集成法 单元(2)对号入 座并累加
矩 阵 位 移 法 单元(3)对号入座
并累加 整体刚度矩阵
连续梁刚度方程
矩 阵 位 移 法
9.5 等效结点荷载向量
矩 阵 位 移 法 加刚臂
去刚臂
(1)加约束求杆端固端弯矩、刚臂约束力矩
矩 阵 位 移 法
(5)集成等效结点荷载向量 形成过程如下:
矩 阵 位 移 法

结构力学 第三十七讲矩阵位移法

结构力学  第三十七讲矩阵位移法
第十一章 矩阵位移法
第十一章 矩阵位移法
教学内容
教学内容:矩阵位移法基本思想,结构离散化,平面刚架单元 的刚度矩阵(局部坐标系、整体坐标系),坐标转换矩阵,单 元定位向量,单元集结构整体刚度矩阵,等效结点荷载,结构 整体荷载列阵,先处理法。 教学要求: 1、了解矩阵位移法的基本概念, 2、理解一般杆单元局部坐标系下的单元刚度方程;单元刚度 矩阵的性质;连续梁单元的单元刚度矩阵(局部坐标系);平 面刚架单元的单元刚度矩阵(局部坐标系);直接结点荷载; 单元杆端力两种坐标系下的转换关系;坐标转换矩阵及性质; 平面刚架单元整体坐标系下的单元刚度方程;整体刚度矩阵的 性质;结构整体刚度方程;矩阵位移法的计算步骤。 3、掌握用矩阵位移法计算连续梁;用矩阵位移法计算刚架。 重点:连续梁和刚架的矩阵位移法求解。 难点:刚架的矩阵位移法方程的建立。
X ,Y, M满足右手法则。
5
6
6
2 3
3
5
4
1
1
4
2
X
Y
第十一章 矩阵位移法
(5)单元杆端位移:
5
每杆端有:两个线位移(轴线、垂
直轴线)、一个角位移(转角)分量。 2
线位移的正方向与坐标正向正负相同, 3
角位移顺时针为正。
1
1
66
3
5
4
4
2
u1
v1 X1
Y
M1
1
Y1
u2
X
M2
v2
2
Y2 X 2
第十一章 矩阵位移法
(2)单元杆端力向量
1 1
u1
v1
2 2
u2
v2
(1) (e)
(
2)
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

4 5 6 7 8 9
{FE }( 2)
= −[T ] F
T
{ }
F ( 2)
总码 总码
ql 2 0 2 ql 12 = ql 2 0 2 − ql 12
总码
0 F − P 2 FP l − 8 = 0 FP − 2 F l P 8
1 2 3 4 5 6
4 5 6 7 8 9
总码
ql 2 0 ql 2 12 = ql 2 0 ql 2 − 12
总码 总码 4 5 6 1 2 3
1 2 3 0 0 0
1. 等效结点荷载的概念
荷载等效的原则是不改变结构的结点位移情况, 荷载等效的原则是不改变结构的结点位移情况, 原则是不改变结构的结点位移情况 即结构在等效结点荷载作用下的结点位移 结点位移与实际荷载 即结构在等效结点荷载作用下的结点位移与实际荷载 作用下的结点位移相等 相等。 作用下的结点位移相等。
F 2
+
=
F = −F
(C)
E 1
F 1
F = −F
E 2
有非结点荷载作用的单元杆 端力,可以由两部分叠加而得, 端力,可以由两部分叠加而得, 一部分是单元固端力, 一部分是单元固端力,另一部分 是杆端位移产生的杆端力。 是杆端位移产生的杆端力。
练习:用先处理法、 练习:用先处理法、后处理法分别计算结构的 综合结点荷载列阵。 综合结点荷载列阵。
M q 结 点 位 移 等 效
F
E 1
F
E 2
1. 等效结点荷载的概念
M (A) q
=
内 q
FF 1
M (B)
F
F 2
力 相 等
+
F = −F
(C)
E 1
F 1
F = −F
E 2
F 2
结点位 移等效
2. 等效结点荷载的计算
2(4,5,6)
FP ② 3(7,8,9)
FP l ql 2 − 8 12
4 5 6 1 2 3
4 5 6 7 8 9
{FE }
(1)
{FE }
( 2)
ql 2 0 2 ql − 12 ql 2 FP {FE } = − 2 2 F l ql − P 12 8 0 F − P 2 FP l 8
M=FPl FP=ql q
C
FP
D
2(1,2,3)
3(4,5,6)
2(4,5,6)
3(7,8,9)
③ ① ②
4(4,5,7)
③ ①
1(1,2,3)
4(7,8,10)

5(11,12,13)
A
B 1(0,0,0) 5(0,0,0)
ql 12
2 T
{F }
解:1. 局部坐标系下单元固端力
= 0 ql − 2 ql − 12
FP l FP l 2 − FP l 2 8
q

1(1,2,3)


ql 2 12
ql 2
没有结 点位移
0
+
0
ql 2 ql − 2 −
F l FP l − P F l − ql − P 2 2 FP l 12 8 8 ql 0 只有结 2 ql 点位移 2 ql 2
2
ql − 2

12
0
ql 2 − 12
总码
总码
ql 2 − 1 12 ql 2 {FE } = 2 FP 3 − 2 ql 2 FP l 4 − 8 12
3. 综合(自由)结点荷载列阵 综合(自由)
{F} = {FD } + {FE }
综合结点 荷载列阵 直接结点 荷载列阵 等效结点 荷载列阵
1 2 3 0 0 0
4 5 6 1 2 3
{FE }
(1)
= −[T ] F
T
{ }
F (1)
0 − 1 0 = −
1 0 0
0 0 1 0 −1 0
[0]
1 0 0
[0]
总码 总码
{FE }
( 2)
0 0 0 = 0 0 0
(2)
局部坐标系下单元等效结点荷载: 局部坐标系下单元等效结点荷载:
{F }
E
e
=−F
{ }
F e
整体坐标系下单元等效结点荷载: 整体坐标系下单元等效结点荷载:
{FE }
e
= [T ] F E
T
{ } = −[T ] {F
e T
F e
}
总码
ql 0 2 ql − 0 2 ql 2 ql 2 − 12 12 = ql 0 0 2 0 ql 1 − 0 2 ql 2 ql 2 12 − 12
在后处理法中, 在后处理法中,结构的结点荷载列阵等于综合结点荷载列阵 与结构约束反力列阵之和。 与结构约束反力列阵之和。
※结构在实际荷载作用下的单元杆端力
M (A) q
{F }
F
F 2
e
= F
{ } + [k ] {δ }
F e e
e
F
(B)
F 1
M
q
对计算杆端位移而言, 对计算杆端位移而言,等效 结点荷载与原非结点荷载作用效 果等效,由此可以断定, 果等效,由此可以断定,在综合 结点荷载作用下求得的即是杆端 的实际位移。 的实际位移。
{FE }
(1)
{FE }
1(0,0,1)
2(2,3,4)
{FE }(1)
ql 2 0 2 2 F 0 3 − P 3 2 2 ql FP l 4 4 − = 12 {F }( 2) = 8 ql 0 E 0 0 2 FP 0 0 − 2 0 2 F l − ql 1 P 0 12 8
第九章 矩阵位移法
9.4 等效结点荷载
非结点荷载的处理
有限元分析的重要一步是把一个连续的结构看成是 由离散单元在结点处连接拼装而成。 由离散单元在结点处连接拼装而成。而把作用在结构上 的荷载统统考虑作用在结点上。 的荷载统统考虑作用在结点上。然而无论是恒载还是活 常常是分布作用在单元上。 载,常常是分布作用在单元上。 对这种非结点荷载的处理: 对这种非结点荷载的处理: 方法一:把分布荷载改用若干集中荷载代替, 方法一:把分布荷载改用若干集中荷载代替,并把集 中荷载的作用点选作结点; 中荷载的作用点选作结点; 方法二:等效结点荷载法。 方法二:等效结点荷载法。
4 5 6 1 2 3
{FE }(1)
= −[T ] F
T
{ }
F (1)
0 − 1 0 = −
1 0 0
0 0 1 0 −1 0
[0]
1 0 0
[0]
总码
0 0 F F P − P 2 2 FP l FP l 8 − 8 = − = 0 0 FP FP 2 − 2 F l F l − P P 8 8
2
F (1)
0
ql − 2
{F }
F ( 2)
= [0 0 0 0 0 0]
T
{F }
F ( 3)
= 0
FP 2
FP l 8
0
FP 2
F l − P 8
T
2. 整体坐标系下单元等效结点荷载
α 1 = α 2 = −90 , α 3 = 0
o
o
总码
ql 0 2 ql − 0 2 ql 2 ql 2 − 12 12 = ql 0 0 2 0 ql 1 − 0 2 ql 2 ql 2 12 − 12
0
FP l FP l 2
FP l F l − P 2
ql 2 12
8
8
0
0
0
局部坐标系下单元固端力列阵
0 ql − 2 ql 2 − 12 = 0 ql − 2 ql 2 12
在矩阵位移法分析过程中, 在矩阵位移法分析过程中,可用综合结点荷载来 代替原来荷载进行计算。 代替原来荷载进行计算。
F3y F3x 3(2,3) ⑥
F4y
4(4,5)
F4x ⑤ ② F1y ③ ④ F2y F1x ① F2x 1(0,0) 2(1,0) FR1x FR2y FR1y
F1 F2 x F F 2 3x {F } = F3 = F3 y F F 4 4x F5 F4 y
在先处理法中,结构的结点荷载列阵等于综合结点荷载列阵。 在先处理法中,结构的结点荷载列阵等于综合结点荷载列阵。
F3y F3x 3(5,6) ⑥
F4y
4(7,8)
F4x ⑤ F1y② ③ ④ F2y F1x ① F2x 1(1,2) 2(3,4) FR1x FR2y FR1y
相关文档
最新文档