韩志强-有限元程序设计.doc

-
. z.
基于Matlab语言按平面三角形单元划分的构造有限元程序设计
专 业： 建筑与土木工程
班 级： 建工研12-2
**：志强
上数值数据为程序运行前输入的初始数据，存为“*t〞文本格式，同时必须放在Matlab工作目录下，路径不对程序不能自动读取指定初始文件，运行出错。初始数据文本文件输入格式如下列图：
Matlab语言程序源代码：
程序中变量说明
-
. z.
NNODE 单元节点数
NPION 总结点数
NELEM 单元数
NVFI* 受约束边界点数
FI*ED 约束信息数组
NFORCE 节点力数
FORCE 节点力数组
COORD 构造节点坐标数组
LNODS 单元定义数组
YOUNG 弹性模量
POISS 泊松比
THICK 厚度
B 单元应变矩阵(3*6)
D 单元弹性矩阵(3*3)
S 单元应力矩阵(3*6)
A 单元面积
ESTIF 单元刚度矩阵
ASTIF 总体刚度矩阵
-
. z.
ASLOD 总体荷载向量
ASDISP 节点位移向量
ELEDISP 单元节点位移向量
STRESS 单元应力
FP1 数据文件指针
-
. z.
Matlab语言程序代码
%******************************************************************************
%初始化及数据调用
format short e %设定输出类型
clear %去除存变量
FP1=fopen('*t','rt'); %翻开输入数据文件，读入控制数据
NELEM=fscanf(FP1,'%d',1), %单元个数〔单元编码总数〕
NPION=fscanf(FP1,'%d',1), %结点个数〔结点编码总数〕
NVFI*=fscanf(FP1,'%d',1) %受约束边界点数
NFORCE=fscanf(FP1,'%d',1), %结点荷载个数
YOUNG=fscanf(FP1,'%e',1), %弹性模量
POISS=fscanf(FP1,'%f',1), %泊松比
THICK=fscanf(FP1,'%f',1) %厚度
LNODS=fsc