特来此请教一个有关刚度矩阵的问题

浏览:68496 回答:1
用百度搜了半天,发现这个有限元论坛最为正规了,所以抱着很大的希望来请教一个问题,望高手不吝赐教:
都说刚度矩阵是奇异矩阵,也就是说它的行列式为零,可是当我求等参单元的刚度阵为啥总是特别大(大概为10的80次方级别),而且我特意找了本书把它的二次四边形的刚度矩阵(16×16的)的行列式计算了一下,竟然为10的100多次方级别的。
也不知道是不是我求解方法有问题,我也算了下二次四边形的刚度矩阵,对称性符合了,而且分别把每行每列的1、3、5、7、9、11、13、15和2、4、6、8、10、12、14、16这些项相加,结果也都接近于零,问题就是行列式的值特别大。
我想是不是因为等参单元节点的实际坐标都映射为标准几何体的节点坐标,比如四边形映射为正方形,六面体映射为正方体,又用了型函数来转化,最后还要积分求刚度矩阵,通过了这一系列的变化,致使它的刚度矩阵本来就可以不为零?
希望有热心的朋友能够粘贴上一个等参单元的刚度矩阵,看到底行列式是否为零
邀请回答 我来回答

全部回答

(1)
默认 最新
csyky2007
我是楼主,我弄明白了,呵呵,原来是我用matlab把算出的符号矩阵转换成数值矩阵时存在舍入误差,刚度矩阵每个项本来就很大(由于E一般都很大),误差再经过连乘的放大,就造成行列式的值特别大,不过刚度矩阵还是对的
2008年12月8日
评论 点赞

没解决?试试专家一对一服务

换一批
    App下载
    技术邻APP
    工程师必备
    • 项目客服
    • 培训客服
    • 平台客服

    TOP