6平面问题总刚度矩阵的组集.docxVIP

平面问题总刚度矩阵的组集 单元刚度矩阵 T扩展T叠加 三角形单元的单元刚度矩阵 单元结点力与单元结点位 移之间的关系 [Ke]{ e} = {Fe} kn k12 k13 k14 k15 kJ i k k21 k22 k23 k24 k k26 k3i k32 k33 k34 k36 [Ke]= 1 1 k41 k42 k43 k44 k46 k k51 k52 k53 k54 k56 k61 k62 k63 k64 k66 J ￡ 1 1 乍JI Pe} 1 卜 {F 1 Fel FjD 1 J 将单元刚度矩阵按单元结点号分块、命名 r^i ^12 沪21上22 ^13 片 14 % 24 片15切 為5 *26 ] ^^2 褊1七42 *33. ^34 ^43 片 44 七35 *36 ^45十酗 島1片52 % ^54 片55片56 _^61 % %魅4 %抵_ Ke Ke I e _Ke Ke Ke Fie 9个子块! 例：设系统的结点数为6, e为系统中的一个三角形单元， 三个定点的编码分别为 1, 5, 6 Ke Ke Ke 1 丨(Fe Ki1 Ki5 5 1 厂1 K：1 k：5 K56 F5e (K：i K：5 K汕 6】底］ 扩展：将各子块按下标放在总刚度矩阵对应的位置 Ke o o o Kie5 Ki訂 o o o o 0 o! e o o o o 0 o [K]扩」 o o o o 0 o k； o o o k55 0! 1爲 o o o o k66j 叠加扩展后的单元刚度矩阵得原始总刚度矩阵 原始总刚度矩阵：解除所有约束的结构的总刚度矩阵 引入边界条件 69页例题 原始方程 [K]{}={P} [K]―― 12 行、12 列 { }二{ 1 2 3 4 5 （}T ={Ui V U V2 u V3 u W U V U 莎 ={0 v1 0 v2 u3 Vs 0 0 u5 0 u6 0}T {PH{P ￡ P3 P4 P5 P6}T ={Rx -10R2X 000 R4x R4y 0 R5y 0 只6『9 （1）手算：划去方程组的第 1、3、7、8、10、12个方程以 及总刚余下的行中的第 1、3、7、8、10、12列元素，得引 入边界条件后的总刚度矩阵 [K]{「?}珥P} [K] ―― 6 行、6 列 { } ={V1 V2 U3 V3 U5 U6}T {P} ={-10 0 0 0 0 0}T （2）机解：置大数法 将位移为零的自由度所对应的总刚度矩阵的对角线元素乘 以一个大数，载荷对应项置零。 例如，上例中U2 7, U2对应的自由度的序号为 3,贝y将原始 总刚中的k33乘以一个大数，比如1010 ,并将P3赋零（原为R2x ）, 这时方程组的第三个方程成为： 10 k31U1 k32V1 10 k33U2 k34V2 k35U3 k36V3 k311U6 k312V^ = 0 注意，方程左边的各项与第三项相比均可忽略，即上述方程 近似等价于：10l°k33U2=0 而k33不为零，所以这时方程所 表示的意义就是 u2 = 0 置大数后各矩阵的改变: [K]的维数不变，与受约束自由度对应的对角线元素被 改变； { S}没有改变； {P}的维数不变，与受约束自由度对应的元素均为零。 等效结点载荷一非结点载荷的处理 静力等效原则：转换后的结点载荷与原载荷在任意虚位移上 的虚功相等 集中力的等效： { e*} T{P{ M }T{円二{呼⑴区 yM)]T{ P} { M }二[N (xm , yM )]{ *} I {pe}||N(XM yM )]T{P^[Nm]T{P} 例：图示三角形单元ijk的边界jk上作用有沿x方向的集 中力P,力的作用点的坐标为 (0.5,0.75 ),求力P的等效结 点载荷 y ■ 1 0. 5 N 0 “ 解：{ p} = * P ,0 [Nm] = 1 0 Nj 0 Nk 0 1 Ni 0 N j 0 Nk x=0.5 y=0.75 { P?} =[ Nm ]T{ p}=⑴求 Nj(x,y)Ni0Nj { P?} =[ Nm ]T{ p}= ⑴求 Nj(x,y) Ni 0 Nj 0 Nk Ni 0 Nj 0 N k x^.5 -y ^.75 fPl - 讣[Ni 0 Nj 0 Nk l.T 0 x」.5 -yj.75 Nj(x,y)二印 a2X a3y Ni(Xi,yJ = M(0,0) =1 = Ni(Xjy)二 Ni(2,0) =0 二 Ni(Xk,yQ 側(1,0) =0 二 3i = 1 a1 2a2 = 0 = a2--2 a1 a3 = 0 = a3 - -1 二 N j (x, y) =1 —0.5x - y 汕(0.5, 0.75) =1-0.5