四边形四节点等参元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