矩阵位移法 单元刚度矩阵 整体刚度矩阵

小好

想问问大家，在矩阵位移法中，怎么用编程实现由刚度矩阵生成整体刚度矩阵？比如说有N个单元，M个节点。单元刚度矩阵是6×6，整体刚度矩阵是3M×3M。请大家指教啊？真的懵了

小好

顶起来，希望大家知道的给我解释下啊

pasuka

参考王勖成的有限单元法一书

小好

谢谢楼上的，我刚下载完了，马上看看

caoer

这个问题很简单，但是文字很难描述。
看看书把，先看1维问题，再2维，3维。
总之总体n个节点， 总体刚度矩阵 就是 n＊维数 x n＊维数

风随影动

我也学习学习

jiangsuwubing

有限元书上都有的

whujsh

很简单，随便找一个含刚架的有限元程序设计的课本看一下就可

refeihc

我再补充一下吧。假设单元刚度矩阵[ke]是6X6的，记为
      [ke]=[k(i,j)],  i=1,2,...,6; j=1,2,...,6
总刚阵[K]是(3M)X(3M)的，记为
      [ZK]=[K(i,j)],  i=1,2,...,3M; j=1,2,...,3M
再设单元节点自由度依次是n1、n2、n3、n4、n5、n6，那么就要把单元刚度阵中的各个元素加到总刚阵中去，对应关系是
   k(1,1) ---> K(n1,n1),    k(1,2) ---> K(n1,n2),    ... ,        k(1,6) ---> K(n1,n6)
     k(2,1) ---> K(n2,n1),    k(2,2) ---> K(n2,n2),    ... ,        k(2,6) ---> K(n2,n6)
                                                     ......
     k(6,1) ---> K(n6,n1),    k(6,2) ---> K(n6,n2),    ... ,        k(6,6) ---> K(n6,n6)

a04051127

单刚节点都对应有总刚编号 对号加到总刚度阵就行

lovebbtyf

想问问大家，在矩阵位移法中，怎么用编程实现由刚度矩阵生成整体刚度矩阵？比如说有N个单元，M个节点。单元刚度矩阵是6×6，整体刚度矩阵是3M×3M。请大家指教啊？真的懵了
小好 发表于 2010-1-12 16:27

我也在研究从单元刚度矩阵，组装为整体刚度矩阵的程序，如果你找到了能应用的程序，麻烦也教我一下喽。