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

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

图一 悬臂钢梁

图二 单元划分与结点编号

Matlab 输出结果

附录Ⅰ：

有限元ANSYS分析结果

采用PLANE183单元(四边形八节点)单元得出的结构Y向最大位移为-0.216E-04。约等于matlab平面四边形八节点等参元结点Y向最大位移-2.4024E-5。

附录Ⅱ：

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

%——— ———— ———— 主 程 序———— —————

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

% 2012年

% 本程序只能处理集中荷载作用下的情况

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

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

% 变量说明

% E v h

% 弹性模量 泊松比 厚度

% NPOIN NELEM NVFIX NNODE NFPOIN

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

% COORD LNODS

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

% FPOIN FORCE FIXED

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

% HK DISP

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

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

clear all

format short e

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

%%读入控制数据

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

v=fscanf(FP1,%f,1); % 泊松比

h=fscanf(FP1,%f,1); %厚度

NELEM=fscanf(FP1,%d,1); %单元数

NPOIN=fscanf(FP1,%d,1); % 总结点数

NNODE=fscanf(FP1,%d,1); %单元节点数

NFPOIN=fscanf(FP1,%d,1); %受力结点数

NVFIX=fscanf(FP1,%d,1); %约束结点个数

LNODS=fscanf(FP1,%f,[NNODE,NELEM]); % 单元定义： 单元结点号(逆时针)

COORD=fscanf(FP1,%f,[2,NPOIN]); % 结点号 x,y坐标(整体坐标下)

FPOIN=fscanf(FP1,%f,[3,NFPOIN]);

% 节点力：结点号、X方向力(向右正)，Y方向力(向上正)

FIXED=fscanf(FP1,%d,[3,NVFIX]);

%约束信息数组(n,3) n:受约束节点数目, (n,1):约束点号

%(n,2)与(n,3)分别为约束点x方向和y方向的约束情况,受约束为1否则为0

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

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

%========平面应力问题