四节点matlab,四边形四节点等参元matlab程序.doc

四边形四节点等参元matlab程序

悬臂钢梁，尺寸如图一所示；v=0.3。h=1，E=2.1e11.

图一 悬臂钢梁

图二 单元划分与结点编号

附录：

%---------------四边形四节点等参元 matlab计算程序----------------------------

% 2013年

% 13级建筑与土木工程 Bruce

% B15-405

% 本程序只输出了节点位移、单元中心点的应力

%*******************************************************************

% 变量说明

% E v h

% 弹性模量 泊松比 厚度

% NPOIN NELEM NVFIX NNODE NFPOIN

% 总结点数 , 单元数, 约束结点个数, 单元节点数 ,受力结点数

% COORD LNODS

% 结构节点整体坐标数组, 单元定义数组,

% FPOIN FORCE FIXED

% 结点力数组， 总体荷载向量, 约束信息数组

% HK DISP

% 总体刚度矩阵,结点位移向量

%******************************

clear all

format short e

FP1=fopen(' sjd.txt','rt'); %打开数据文件

%%读入控制数据

E=fscanf(FP1,'%f',1); %弹性模量

v=fscanf(FP1,'%f',1