首页/文章/ 详情

手算Q4单元刚度矩阵

6月前浏览4631


平面四边形等参单元(Q4)的刚度矩阵

由前文可知

 

k是    矩阵,若将    看作函数    ,则    也是    列阵。

 

   是积分点的坐标。将k分块,即

 
 
 
 

如图是一个单元,    。采用4个高斯积分点计算单元刚度矩阵。

四个积分点坐标分别为:

 
 
 
 

积分点权重皆为1

 
 
 
 

于是

附python代码:










import numpy as np
A = np.array([ [-1.5773 , 1.5773 , 0.4227 , -0.4227],[-1.5773 , -0.4227 , 0.4227 , 1.5773] ] )B = np.array([ [3,2],[5,2],[5,4],[3,4] ])C = np.dot(A,B)C = 0.25 *C

同理可得

 
 
 
 

同理可得

 

由此可得

 

其余计算相同。

单元刚度矩阵特点:

1.对称性

2. 奇异性

3. 主对角元素恒正

4. 所有奇数(偶数)行的和为 0

刚度是表示物质形变能力的一个量,也就是说物体抵抗变形的能力,其元素值为单位位移所引起的节点力,与普通弹簧的刚度系数具有同样的物理本质。或者说,是物体产生单位的位移所需要加载的载荷量。刚度矩阵和刚度概念相似,就是把刚度变到了多维 比考虑了在多维的情况下 各个维度的相关性。


★★★★往期相关★★★★

平面四边形等参单元(Q4)的刚度矩阵

数值积分|二元函数的高斯积分

来源:数值分析与有限元编程
pythonUM
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-04-02
最近编辑:6月前
太白金星
本科 慢慢来
获赞 5粉丝 8文章 318课程 0
点赞
收藏
作者推荐

平面四边形等参单元(Q4)的刚度矩阵

坐标系(物理坐标系)下Q4单元的刚度矩阵为 此时应变矩阵 是 的函数 其中 单元刚度矩阵,等效节点荷载,单元应力,应变等物理量是通过 坐标系表达,而在计算时却是在 坐标系下。因此 此时应变矩阵 是 的函数两个坐标系下坐标转换的桥梁为 其中 是 坐标系中单元四个顶点坐标, 是形函数 坐标系中任意函数 在坐标系 的表达式为 ,根据链式求导法则 或者 其中, 是雅可比矩阵 其中 由此可得 将 换成形函数 例如 在 坐标系下 则 k是 矩阵,若将 看作函数 ,则 也是 列阵。 注意 是积分点的坐标。等参数单元(简称等参元)就是对单元几何形状和单元内的参变量函数采用相同数目的节点参数和相同的形函数进行变换而设计出的一种单元类型。优点:由于等参变换的采用使等参单元的刚度、质量、阻尼、荷载等特性矩阵的计算在规则域内进行,因此不管矩阵表示的被积函数如何复杂,仍然可知以方便地采用标准化的数值积分方法计算。也正因为如此,等参元已成为有限元法道中应用最为广泛的单元形式。来源:数值分析与有限元编程

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈