9矩阵法.pptVIP

第九章 矩 阵 位 移 法 传统人工手算： 速度低，精度差。 忌繁重的计算工作量。 电子计算机计算： 速度快，精度高。 忌无规律可循。 手算怕繁，电算怕乱！ ?把握共性，抓住本质—— 本是同根生，同根又同源 ?区分差异，理解内涵—— 手算：简化、技巧（怕繁！） 电算：整齐、统一（怕乱！） ?突出重围，理清脉络 防止头晕眼花 二、局部坐标系（单元坐标系、杆件坐标系） 根据单元分析（杆件）与整体（结构）分析的不同需要，采用两种直角坐标系。 局部坐标系以杆轴为x轴，“ 1”为始端，“ 2 ”为终端。1 2为正方向。 局部坐标系下所有量值的正负号规定：杆端位移和杆端力分量与局部坐标方向一致为正。 二、单元刚度矩阵的性质 （1）杆端位移一律用绝对位移。即：除杆端相对位移外，还包含有刚体位移。（请比较位移法） （3）、 [ k ]e是对称矩阵。 由反力互等定理可知,单元刚度矩阵是对称矩阵。 （4）、 [ k ]e6×6中所有元素是（E、A、I、l）的函数，正负号根据单元坐标系的方向而定。 注意：根据单元刚度矩阵，可由{ Δ }e求出{F}e ，且解是唯一的。但不可由{F}e求{Δ}e ，其结果可能无解或非唯一解。这是正反两个问题，不可混淆。 解释：一般单元的单元刚度矩阵之所以为奇异矩阵，是因为计算的单元是两端无任何支承的自由单元。单元本身除弹性变形外，还有任意的刚体位移。 {F}e完全一样，但{Δ}e可以不同。对应于一个平衡力系，可以有多种杆端位移情况。 （6）、单元刚度矩阵可以分块 以平面刚架单元为例，单元 ij 两端点i、j 的位移分量和力的分量可表示为： 一、 单元集成法的力学模型和基本概念 分别考虑每个单元对{F}的单独贡献，整体刚度矩阵由单元直接集成 i1 i2 1 2 ?1 ?2 ?3 F3 {F} 1 = [ F1 1 F2 1 1 ]T F1 1 F2 1 F3 1 令 i2 =0，则 F3 1 =0 [k] = 4i1 2i1 4i1 2i1 1 F1 1 F2 1 = 4i1 2i1 4i1 2i1 ?1 ?2 （a） （b） F1 1 F2 1 F3 1 = 4i1 2i1 4i1 2i1 0 0 0 0 0 ?1 ?2 ?3 1 [K] {?} {F} = 1 [K] = 1 4i1 2i1 4i1 2i1 0 0 0 0 0 单元 1 的贡献矩阵 单元 1 对结点力{F}的贡献 略去其它单元的贡献。 i1 i2 1 2 ?1 ?2 ?3 F1 2 F2 2 F3 2 [k] = 4i2 2i2 4i2 2i2 2 F1 2 F2 2 F3 2 = 4i1 2i1 4i1 2i1 0 0 0 0 0 ?1 ?2 ?3 2 [K] {?} {F} = 2 设 i1 =0，则 F1 2 =0 [K] = 2 4i1 2i1 4i1 2i1 0 0 0 0 0 单元 ? 的贡献矩阵 F3 {F} 2 = [ F1 2 F2 2 2 ]T 单元?对结点力{F}的贡献 略去单元?的贡献。 1 [K] {?} {F} = 1 [K] = 1 4i1 2i1 4i1 2i1 0 0 0 0 0 2 [K] {?} {F} = 2 [K] = 2 4i1 2i1 4i1 2i1 0 0 0 0 0 i1 i2 1 2 1 2 1 2 [K]=（[K] +[K] ）= 1 2 e e [k] [K] [K] e e {F}={F} +{F} =（[K] +[K] ）{?} 1 2 {F}=[K]{?} 整体刚度矩阵为： 单元集成法求整体刚度矩阵步骤： 根据单元?和单元?分别对结点力{F}的贡献，可得整体刚度方程： [k] [K] [K] e e 1 2 [k] = 4i1 2i1 4i1 2i1 1 [K] = 1 4i1 2i1 4i1 2i1 0 0 0 0 0 [k] = 4i2 2i2 4i2 2i2 2 [K] = 2 4i2 2i2 4i2 2i2 0 0 0 0 0 1 2 1 4i1 2i1 4i1 2i1 0 0 0 0 0 2i2 2i2 4i2 [K]= 4i1 2i1 4(i1+i2) 2i1 0 2i2 0 2i2 4i2 4i1