有限元大作业matlab课程设计例
子
2020年4月19日
1
文档仅供参考
有 限 元 大 作 业 程 序 设 计
学校：天津大学
院系：建筑工程与力学学院专业：01级工程力学 姓名：刘秀 学号：\\\\\\\\\\\\\\\\\\\\\\ 指导老师：
2020年4月19日
2
文档仅供参考
连续体平面问题的有限元程序分析
[题目]：
如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界上受正向分布压力，
p?1kNm，同时在沿对角线
y轴上受一对
集中压力，载荷为2KN，若取板厚t?1，泊松比v?0。
2k1kN/2k
[分析过程]：
由于连续平板的对称性，只需要取其在第一象限的四分之一部分参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。
3
2020年4月19日
文档仅供参考
[程序原理及实现]：
用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN，经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。 该程序的特点如下：
问题类型：可用于计算弹性力学平面问题和平面应变问题 单元类型：采用常应变三角形单元 位移模式：用用线性位移模式
载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷 材料性质：弹性体由单一的均匀材料组成
约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至
少应有对三个自由度的独立约束
方程求解：针对半带宽刚度方程的Gauss消元法
4
2020年4月19日
1kN/
