平面桁架单元刚度矩阵推导

2020年02月19日 星期三

这里只讲杆单元刚度矩阵的推导,基于前提是已经根据势能最小原理得到了局部坐标系下单元刚度矩阵表达式,即对于弹性杆单元,通用单元刚度矩阵表达式为……

​ 这里只讲杆单元刚度矩阵的推导,基于前提是已经根据势能最小原理得到了局部坐标系下单元刚度矩阵表达式,即对于弹性杆单元,通用单元刚度矩阵表达式为:

1 基本信息

——沿杆单元方向建立局部坐标,起点为0,端点为杆单元长度

——沿杆单元方向位移

对于其他符号意义可参见于另一文:有限元推导基本符号

参考书籍:辛科维奇《有限元方法》第五版曾攀等译

2 形函数

​ 采用两节点拉格朗日多项式形式的插值函数(关于插值函数可参见于《有限元方法》第五版第八章) 即这里采用的这里形函数为:

3 等参单元

​ 由于形状规则,对于上述坐标选取下积分依旧方便计算,故理论上可以不采用单元映射方法(更多可见于可参见于《有限元方法》第五版第九章)。不过这里通用起见,采用等参单元,即另取局部坐标,满足: 相应形函数为: 雅可比矩阵为:

4 局部坐标系下单元刚度矩阵

这里: 则: 故: 由于积分容易得到,故不采用高斯积分方法而直接求解定积分: 如果扩展到y向,则有:

5 整体坐标系下单元刚度矩阵

存在转换矩阵 (这里符号与前述雅可比矩阵相同,但矩阵并不相同)使得:

故: 则: 相应有整体坐标系下单元刚度矩阵为:

6 整体刚度矩阵

​ 每一个单元通过节点转换矩阵集成入整体刚度矩阵中去,即: 相应的如果有节点荷载向量为,则得到最终有限元方程为: 上述杆单元刚度矩阵的推导可以认为是最简单的了。对于实际的复杂单元,上述求解单元刚度矩阵过程中需要进行积分。由于一般情形下无法求出不定积分,故通常采用高斯积分形式。

7 释疑

​ 对于梁单元亦可采用上述方法进行推导,只不过形函数、应变微分算子、物理方程等不相同。

​ 上问推导可能存在的一个困惑是怎么会出现两个局部坐标系。其原因在于本文推导时为讲究方法通用性应用了等参变换,实际上对于本例而言是没有必要的,不进行等参变换依旧能得到相同结果。由于采用了等参变换,第一个局部坐标就相当于第二个局部坐标的整体坐标了。


2020年2月19日 星期三 天气多云转阴

精选博客

word的按样式输出

在编写一个标准,写的时候呢就会把条文正文和条文说明写在一块,方便编写和对照查看(如下图所示)。但是呢最终的成文必然是两者需要分开的。如果条文很多的话就会非常麻烦。可是word又没有图层功能,可以快速筛选出两者。于是尝试采用下python的相关扩展库来实现这一功能。……

继 续 阅 读

高斯积分相关

对于一般的积分表达式,是无法求出其定积分形式的。对于有限元问题而言,主要就是单元刚度矩阵的问题。(当然对于梁单元而言,是可以给出单元刚度矩阵显式表达式的)。因此得用到数值积分方式。而这其中,如果积分具有如下形式,则可以应用高斯积分。……

继 续 阅 读

底层倾覆力矩的问题

抗规中对于框架-剪力墙结构中基于底层框架倾覆力矩占总倾覆力矩的比值进行分类设计。并在条文说明中给出了计算底层框架倾覆力矩的公式,但这个公式是否合理?本文给出相关推导并进行适当讨论……

继 续 阅 读