[FEM][有限元][编程][Matlab][Code by myself] 2D Timoshenko梁单元

(  有空和小伙伴一起写写有限元程序 ) 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所;2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 欧拉梁单元基于一定的假设(Kirchhoff假设),在梁的高度远小于其跨度的时候,可以忽略梁的横向剪切变形,此时采用欧拉梁单元进行模拟,能够得到较为满意的结果。但对于跨高比较小的深梁,梁的剪切变形将引起附加挠度,使得原来垂直于轴线的截面在变形后将不再与轴线垂直,且发生翘曲。此时需采用能够考虑横向剪切变形的梁单元进行模拟。 考虑剪切变形修正的经典梁单元和Timoshenko梁单元是两种较为常用的能够考虑梁剪切变形的梁单元,但这两种梁单元仍假定原来垂直于中面的截面在梁变形后仍保持为平面。 在列式方面,Timoshenko梁的基本特点是将挠度和截面转角分别插值。Timoshenko梁单元应用广泛,也容易据此推广到板壳单元 问题描述( Problem Description) 一榀XZ平面内的刚架结构,结构几何信息如图所示;节点1、5处为固定支座,节点4处受到+x方向P=200kN的集中力作用。结构中各杆件采用相同的材料,弹性模量E=30000MPa,梁、柱截面面积分别为0.08m2和0.16m2,梁、柱截面惯性矩分别为0.0128/12m4和0.0256/12m4,梁柱抗剪面积分别为0.0667m2和0.1333m2,材料泊松比为0.2。 基于MATLAB编程实现该框架结构的弹性静力分析,所有构件采用剪切修正梁单元进行模拟,并将基于MATLAB编程计算的结果与Abaqus分析结果进行对比。 MATLAB 编程 Abaqus 位移结果对比: 支座反力对比: 可以看出,MATLAB编程计算结果是和Abaqus计算结果是完全一致的。 …