[FEM][有限元][编程][Matlab][Code by myself] 平面4节点线性“完全积分”单元(CPS4)

  (业余时间 和小伙伴一起写写有限元程序 ) 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所; 2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 平面4节点线性“完全积分”单元,用Abaqus中的命名规则,即CPS4。 问题描述( Problem Description) XZ平面内的悬臂梁,悬臂长度2.0m,梁高0.5m,梁宽0.2m。梁左端嵌固,右端受到-z方向的集中力1000kN。材料弹性模量E=200000MPa,材料泊松比为0.3。 MATLAB编程实现该悬臂梁的弹性静力分析,采用Q4平面单元进行模拟,并将MATLAB的计算结果与SAP2000、midas Gen分析结果进行对比。 MATLAB 编程 disp stress SAP2000 …

[FEM][有限元][编程][Matlab][Code by myself] Constant Stress Triangle Element 三角形常应变单元(CST)

(  业余时间 和小伙伴一起写写有限元程序 ) 程序作者 ( Author ) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1.广州容柏生建筑结构设计事务所;2.华南理工大学建筑设计研究院 基本概念 ( Concept ) 三节点三角形单元有三个节点和三个直边,单元内应力为常数,因此称为常应变三角形单元,即Constant Stress Triangle Element(常应力三角形单元),简称CST单元。 问题描述( Problem Description) 平面内的悬臂梁,悬臂长度2.0m,梁高0.5m,梁宽0.2m。梁左端嵌固,右端受到-z方向的集中力1000kN。材料弹性模量E=200000MPa,材料泊松比为0.3。 基于MATLAB编程实现该悬臂梁的弹性静力分析,采用CST平面单元进行模拟,并将MATLAB的计算结果与SAP2000、Midas Gen分析结果进行对比。 MATLAB 编程 mesh …

[Midas2Abaqus] Verification Example Static-31: Elliptic membrane under uniformly distributed load [ Midas Gen 验证实例 Static-31: 椭圆薄膜受均布荷载]

【1】Problem Description ( Midas Gen’s Verification Examples Static-31:  Elliptic membrane under uniformly distributed load) An elliptical membrane of thickness, t and geometrical properties is subjected to a uniformly distributed outward pressure. …

[Abaqus] Plate with Hole Stress Analysis [带孔平板应力分析]

利用ABAQUS进行弹性圆孔薄板的应力分析,复习弹性力学知识。 例子(Example) 如图所示为一带圆孔的平面薄板,圆孔的直径相对板的尺寸较小,板受两侧均匀平面拉力作用。根据弹性力学理论,圆孔的A点和B点存在应力集中,当板接近无限大时,A点应力的绝对值约为施加的水平应力绝对值的3倍,B点应力的绝对值约为施加的水平应力的绝对值。以下采用ABAQUS对该问题进行分析。 有限元模拟 (FEM Analysis) 算例中的平面板厚度较薄,且为平面受力,可简化为平面应力问题,由于板两个方向具有对称性,因此可以等效为图右的1/4模型的分析。 采用ABAQUS进行分析,单元采用最简单的常应力应变三角形单元(Constant Strain Triangle)。基本步骤如下: (1)创建部件(Create Part) 模型空间选为 2D Planar, 类型为 Deformable,基本特性为 Shell,绘制如下部件轮廓。 (2)定义材料 (Define Material) 定义弹性材料,输入弹性模量和泊松比 (3)定义截面 (Define Section) 定义截面,属性为固体和各向同性,并指定厚度。 (4)划分部件网格 (Mesh Part) 给网格布种,这里按边布种,并指定相应的偏心参数,使靠近圆孔附近的网格更密。指定单元属性为三角形,单元阶次为Linear,相应的单元名字为CPS3,对应为三角形平面应力单元,划分网格: …