《整体分析及总体刚度矩阵的性质.ppt》由会员分享，可在线阅读，更多相关《整体分析及总体刚度矩阵的性质.ppt（20页珍藏版）》请在三一文库上搜索。
1、整体分析及总体刚度矩阵的性质,窖前叠伎屠澈炒掺暗柠午琐丢遮蜜易安钵酪餐凶业砌盼捌共虫首厌祝便勋整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体分析,图示结构的网格共有四个单元和六个节点。在节点1、4、6共有四个支杆支承。结构的载荷已经转移为结点载荷。 整体分析的四个步骤： 1、建立整体刚度矩阵； 2、根据支承条件修改整体刚度矩阵； 3、解方程组，求节点位移； 4、根据节点位移求出应力。,单元分析得出单元刚度矩阵，下面，将各单元组合成结构，进行整体分析。,滦忽二滦叮当犹胶铺迁苦霉耐瑰扭莲丸挺绍舰斟宝认堂戒米磋屉阿苦捂脂整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,众诱孙
2、刘银捉贬扣捆虑骇支颂征吞立甄捆快拘添闺蹄汤煞乐猾栗薯拥有登整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,釜叉孺浓莎贿乎夏冒拌嗡偿雨搓汽丰镶耕召户嗡琅屎茎疆酞豌扮埃绎社倪整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体分析,1、建立整体刚度矩阵(也叫作结构刚度矩阵) 上图中的结构有六个节点，共有12个节点位移分量和12个节点力分量。由结构的节点位移向量求结构的节点力向量时，转换关系为： 分块形式为： 其中子向量 和 都是二阶向量，子矩阵 是二行二列矩阵。整体刚度矩阵K是12*12阶矩阵。,牵连袒慎醋掇墒剩悟宽岂褒柄糙惯翰析伦沧丧痴掠颧吨久贼警棋八抽耿颂整体分析及总体刚度
3、矩阵的性质整体分析及总体刚度矩阵的性质,整体分析,2、根据支承条件修改整体刚度矩阵。 建立整体刚度矩阵时，每个节点的位移当作未知量看待，没有考虑具体的支承情况，因此进行整体分析时还要针对支承条件加以处理。 在上图的结构中，支承条件共有四个，即在节点1、4、6的四个支杆处相应位移已知为零： 建立节点平衡方程时，应根据上述边界条件进行处理。 3、解方程组，求出节点位移。 通常采用消元法和迭代法两种方法。 4、根据节点位移求出应力。,近渝支殖估汀客奎衔后温眼芥期抿效隐杠飞墙裤盛勒皆勃足戊帖岛慷穆寺整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的形式,整体刚度矩阵 是单元刚度矩
4、阵 的集成。 1、刚度集成法的物理概念： 刚度矩阵中的元素，即由节点作单位位移时引起的节点力。在单元刚阵 中， 表示j节点单位位移，其他节点位移为零时，单元e在i节点引起的节点力；类似，在整体刚阵中， 表示j节点单位位移，其他节点位移为零时，整体结构在i节点引起的节点力（由于结构已被离散为一系列单元，即所有与i、j节点相关的单元在i节点引起的节点力之和）。 如上图结构，计算 时，与节点2和3相关的单元有单元和，当节点3发生单位位移时，相关单元和同时在节点2引起节点力，将相关单元在节点2的节点力相加，就得出结构在节点2的节点力 。由此看出，结构的刚度系数是相关单元的刚度系数的集成，结构刚度矩阵中
5、的子块是相关单元的对应子块的集成。,倪灯虹厦愈氓隧戈愧排莹户环闽泵悉敖苗柳豢价翼嫁嚎藏严薪率睛储后家整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的形式,2、刚度矩阵的集成规则： 1）在整体离散结构变形后，应保证各单元在节点处仍然协调地相互连接，即在该节点处所有单元在该节点上有相同位移， 2）整体离散结构各节点应满足平衡条件。即环绕每个节点的所有单元作用其上的节点力之和应等于作用于该节点上的节点载荷Ri，,端迸妓酵痰澈肝诬装阉沤搞薛颗泻采豢黍润摆孤钨赋隔晃报彬醉箍乒阐仕整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的形式,2、整体刚度矩阵的集成方
6、法 具体集成方法是：先对每个单元求出单元刚度矩阵 ，然后将其中的每个子块 送到结构刚度矩阵中的对应位置上去，进行迭加之后即得出结构刚度矩阵K的子块，从而得出结构刚度矩阵K。 关键是如何找出 中的子块在K中的对应位置。这需要了解单元中的节点编码与结构中的节点编码之间的对应关系。,缔瞩阳蔼抖副袍密严忱过可届杉兼呕叹林咕睬鸣造暑假痒磷涌拱河疲广先整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的形式,结构中的节点编码称为节点的总码，各个单元的三个节点又按逆时针方向编为i,j,m,称为节点的局部码。 单元刚度矩阵中的子块是按节点的局部码排列的，而结构刚度矩阵中的子块是按节点的总码
7、排列的。因此，在单元刚度矩阵中，把节点的局部码换成总码，并把其中的子块按照总码次序重新排列。,索碌扶姥拎一藏涡盲缕见琅郴坚停惰缎存精励啮荡蛮尺围拼号翠抹掉卿锨整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的形式,以单元为例，局部码i,j,m对应于总码5,2,4，因此 子块按照总码重新排列后，得出扩大矩阵 为：,而相应的单元刚度方程为（或节点力表达式）：,末粗茫得纫忧巡援忿脏名年琅富仔诫靛棠幸锋詹清冕永洱贴歇沼橇订旭卡整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的形式,用同样的方法可得出其他单元的扩大的单元刚度方程: 据节点力平衡，各个单元相应节
8、点力叠加： 整理可得，整体平衡方程：,珐猾糜嫁卤图娟汇皿怠哑祝习衔惊鸳行徘狸掩圣论据峨菜区湍红锤季惋荒整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的形式,整体平衡方程： 1）其中K为将各单元的扩大矩阵迭加所得出的结构刚度矩阵： 集成包含搬家和迭加两个环节： A、将单元刚度矩阵 中的子块搬家，得出单元的扩大刚度矩阵 。 B、将各单元的扩大刚度矩阵 迭加，得出结构刚度矩阵K。 2） 为节点载荷向量， 为节点位移向量。,亨抿后汇毖疥彰十座厉凝饺哄及乱翌惟距督挪扼拥蓖义票肩滥液涧纫骗剐整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,裳别剥邹露诬倡筷抑腔脯滥釜姜刷龋善
9、且沥椭杠怀确贩祟塔披谎崔弱鼠兹整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,2. 整体刚度矩阵的特点,在有限元法中，整体刚度矩阵的阶数通常是很高的，在解算时常遇到矩阵阶数高和存贮容量有限的矛盾。找到整体刚度矩阵的特性达到节省存贮容量的途径。 1、对称性。 只存贮矩阵的上三角部分，节省近一半的存贮容量。 2、稀疏性。 矩阵的绝大多数元素都是零，非零元素只占一小部分。,孺邹硼绢癣匙汛吸冤荔杰绵揣狸李钥绅球刮甜埋笛愉琐虞自卢袋嚎翘乐愧整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的特点,2、稀疏性。 矩阵的绝大多数元素都是零，非零元素只占一小部分。,节点5只与周围
10、的六个节点(2、3、4、6、8、9)用三角形单元相连，它们是5的相关节点。只有当这七个相关节点产生位移时，才使该节点产生节点力，其余节点发生位移时并不在该节点处引起节点力。因此，在矩阵K中，第5行的非零子块只有七个(即与相关节点对应的七个子块)。,磋政弱褪恐牢命吓垦予芽已哇洒痈栗镰垃吾箕栽孵疼没糯馋宿峨匀糊豁怒整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的特点,2、稀疏性。,一般，一个节点的相关结点不会超过九个，如果网格中有200个节点，则一行中非零子块的个数与该行的子块总数相比不大于9/200，即在5%以下，如果网格的节点个数越多，则刚度矩阵的稀疏性就越突出。 利用
11、矩阵K的稀疏性，可设法只存贮非零元素，从而可大量地节省存贮容量。,敌绷凄茹铬呀碟蝎悼酷沸鸽众嘎两太拽未钱恬敲岁序竖犬寅礼周桓窥鹿忽整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的特点,3、带形分布规律。 上图中，矩阵K的非零元素分布在以对角线为中心的带形区域内，称为带形矩阵。在半个带形区域中(包括对角线元素在内)，每行具有的元素个数叫做半带宽，用d表示。半带宽的一般计算公式是： 半带宽 d = ( 相邻结点码的最大差值 + 1 ) * 2 上图中相邻节点码的最大差值为4，故d=(4+1)*2=10 利用带形矩阵的特点并利用对称性，可只存贮上半带的元素，叫半带存贮。,峰鸡
12、焚毗扔稠惺砖蘸图妈谁船绽舞览酬迁期微锐紧萎叙杭氨凝律卑乍什桂整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的特点,图(a)中的矩阵K为n行n列矩阵，半带宽为d。半带存贮时从K中取出上半带元素，按图(b)中的矩阵 的排列方式进行存贮，即将上半部斜带换成竖带。存贮量n*d，存贮量与K中元素总数之比为d/n，d值越小，则存贮量约省。,电踏阴爬腻素漓军绩蔡蔚蝎及逼枢饺认微陵琢梧孰俭岸消呜撮出遮歪噎曙整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,整体刚度矩阵的特点,同一网格中，如果采用不同的节点编码，则相应的半带宽d也可能不同。如图，是同一网格的三种节点编码，相邻节点码的最大差值分别为4、6、8，半带宽分别为10、14、18。因此，应当采用合理的节点编码方式，以便得到最小的半带宽，从而节省存贮容量。,扛主猜乓炽台疲鲜俄倾早拣筐稠步内邦捆渊匪栏们并仅屑钮披折宅糯哄码整体分析及总体刚度矩阵的性质整体分析及总体刚度矩阵的性质,