matlab求刚度,求整体刚度矩阵matlab程序

班级：机自11-13 学号：03111323

求整体刚度矩阵。

部分分析思路：

各单元信息

单元编号

—

①

②

③

④

整体编码 1,2,3 2,4,5 5,3,2 3,5,6

局部编码

i, j,m

i, j,m

i, j,m

i, j,m

以整编码体表示 的单元刚度矩阵

K ⑴

K ⑴

K ⑴

K

11 K

12 K

13 K ⑴

K ⑴

K ⑴

K

21 K

22 K

23 K ⑴

K ⑴

K ⑴

K

31 K

32 K

33 K ⑵

K ⑵

K ⑵

K

22 K

24

K

25

K ⑵

K ⑵

K ⑵

K

42 K

44 K

45 K ⑵

K ⑵

K ⑵

K

52 K

54 K

55 K ⑶

K ⑶

K ⑶

K

55 K

53 K

52

K ⑶

K ⑶

K ⑶

K

35 K

33 K

32 K ⑶

K ⑶

K ⑶

K

25 K

23 K

22 K ⑷ K ⑷

K ⑷

K

33 K

35

K

36

K ⑷

K ⑷

K ⑷

K

53 K

55 K

56

K ⑷

K ⑷

K ⑷

K

63 K

65 K

66

整体刚度矩阵

姓名：吴佳侣

例：如图所示有限元模型，弹性模量为

E =1，厚度为t =1，为简化计算取

，

K 12 K 13 K 14 K 15 K 16

K 22 K 23 K 24 K 25 K 26 K 32 K 33 K 34 K 35 K 36 K 42 K 43 K 44 K 45 K 46 K 52 K 53 K 54 K 55 K 56 K 62 K 63 K 64 K 65 K 66