平面三角形单元常应变单元matlab程序的编制.pdf

平面三角形单元常应变单元matlab程序的编制.pdf

  1. 1、本文档共12页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
三角形常应变单元程序的编制与使用 有限元法是求解微分方程边值问题的一种通用数值方法,该方法是一种基于 变分法(或变分里兹法)而发展起来的求解微分方程的数值计算方法,以计算机 为手段,采用分片近似,进而逼近整体的研究思想求解物理问题。 有限元分析的基本步骤可归纳为三大步:结构离散、单元分析和整体分析。 对于平面问题,结构离散常用的网格形状有三角形、矩形、任意四边形,以三个 顶点为节点的三角形单元是最简单的平面单元,它较矩形或四边形对曲边边界有 更好的适应性,而矩形或四边形单元较三节点三 角形有更高的计算精度。 开始 Matlab 语言是进行矩阵运算的强大工具,因 此,用 Matlab 语言编写有限元中平面问题的程序 输入初始数据 有优越性。本章将详细介绍如何利用 Matlab 语言 编制三角形常应变单元的计算程序,程序流程图 生成单刚集成总刚 见图 1。 有限元法中三节点三角形分析结构的步骤如 施加约束信息 下: 1)整理原始数据,如材料性质、荷载条件、 约束条件等,离散结构并进行单元编码、 生成荷载向量 结点编码、结点位移编码、选取坐标系。 2 )单元分析,建立单元刚度矩阵。 边界条件处理 3 )整体分析,建立总刚矩阵。 4 )建立整体结构的等效节点荷载和总荷载矩 计算结点位移 阵 5 )边界条件处理。 6 )解方程,求出节点位移。 计算单元应力 7 )求出各单元的单元应力。 8 )计算结果整理。计算结果整理包括位移和 计算结果整理 应力两个方面;位移计算结果一般不需要 特别的处理,利用计算出的节点位移分 结束 量,就可画出结构任意方向的位移云图; 而应力解的 误差表现在单元内部不满足平衡方程,单 图 1 程序流程图 - 1 元 与单元边界处应力一般不连续,在边界上应 力解一般与力的边界条件不相符合。 1.1 程序说明 %******************************************************************* % 三角形常应变单元求解结构主程序 %*******************************************************************  功能:运用有限元法中三角形常应变单元解平面问题的计算主程序。  基本思想:单元结点按右手法则顺序编号。  荷载类型:可计算结点荷载。  说明:主程序的作用是通过赋值语句、读取和写入文件、函数调用等完成算 法的全过程,即实现程序流程图的程序表达。 % 1 程序准备 format short e %设定输出类型 clear all %清除所有已定义变量 clc %清屏  说明: format short e - 设定

文档评论(0)

wx171113 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档