首页/文章/ 详情

中厚板边缘效应

4月前浏览3046

在    中厚板理论中,存在一类特殊的难题--边缘效应现象。边缘效应是指:在靠近板的某种边界的狭小区域内,转角与内力场发生剧烈的梯度变化,其中自由边界和软简支边界附近发生的边界效应现象最为强烈。能否准确捕捉这一剧烈变化对于一些工程结构的设计,比如超大型浮式结构,具有特殊的重要意义。但是对于这一问题,由于传统有限元法的固有局限性--位移型有限元应力精度相对较低,根本无法有效计算,只能通过不断加密网格或者自适应网格加密来处理,然而,即使采用极端精细的网格,也不一定能够准确反映边界效应及给定面力的边界条件。以开洞剪力墙的内力分析为例,说明这一问题。

▲图1

如图1所示为一个开门洞剪力墙的局部模型。其中,楼层高度为    ,取    宽局部进行分析,在该部分开了一个    高,    宽的门洞,墙的厚度取    、    和    三种情况。材料弹性模量    ,泊松比    。在所分析板结构角点作用一垂直于板的集中荷载    ,在同样网格密度下,用    所带的板/壳单元、厚板/壳单元以及实体单元进行分析。为比较方便,分析    截面右侧中点    的三个内力合力:剪力    、弯矩    和扭矩    。计算结果如表1所示。

▲表1 不同厚度的开洞剪力墙内力计算结果比较

表1中,几种模型都取    作为加密单元边长。    薄板/壳单元和厚板/壳单元的分析结果基本一致,但与其实体单元在同样单元尺寸情况下计算的内力合力的结果有很大的差异。其中,在墙厚    的情况下,弯矩存在    的误差,扭矩存在    的误差,剪力存在    的误差;在板厚为    的情况下,弯矩的误差为    ,扭矩误差为    ,剪力误差为    ;在板厚为    的情况下,弯矩的误差为    ,扭矩误差为    ,剪力误差为    


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

有限元| 刚架临界荷载的新算法

针对前一篇有限元 | 有限元法计算刚架的临界荷载中计算刚架临界荷载的不足之处,提出新的算法。梁单元新的弹性刚度矩阵和几何刚度矩阵▲图1图1所示为杆件单元上的任一微段 , 是其失稳之前的位置, 为轴压力达到临界值时可能出现的分支平衡位置。设微段由平衡位置 发生无限小的横向虚位移至 。弹性稳定分析是二阶分析,虚位移应该在变形后。 和 分别为虚位移发生前、后微段的长度。微段的轴向虚应变可表示为 由弧长公式 (2)代入(1),并略去高阶微量,得由横向虚位移产生的轴向虚应变为 在梁单元的4个节点自由度中, 和 的单位不同,现在用 将单位统一起来。其中 是单元长度。梁的挠度 其中 于是 两种坐标系的映射 由(9)(10)得 梁横截面虚应变 其中 ▲图2如图2所示,单元轴向荷载所作的虚功为 式中单元的轴向荷载 以受拉为正。(3)代入(4)可得 内力虚功 由虚功原理 可得 记 则有 上式是一个求广义特征值问题。经积分计算后可得 算例▲图3已知简支梁的长度为 ,划分2个单元,求其临界荷载。 考虑边界条件之后,整体弹性刚度矩阵 整体几何刚度矩阵 由 得 其中 ,在MATLAB中求特征值与特征向量的代码如下A = [2,3,1,0; 3,12,0,-3; 1,0,4,1; 0,-3,1,2];B = [4/3, 1, -1/3 0; 1, 24, 0, -1; -1/3, 0, 8/3, -1/3; 0, -1, -1/3, 4/3];[X,D] = eig(A,B)% D是对角矩阵,对角线上的元素是特征值。 D中最小的值是0.1243,则 结果与材料力学相同。特征值0.1243对应得特征向量为 则 便是屈曲模态。来源:数值分析与有限元编程

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