四边形四节点等参元matlab 程序 2013 年12 月28 日 — 1 — 悬臂钢梁，尺寸如图一所示；v =0.3。h=1，E=2.1e11. 图一 悬臂钢梁 图二 单元划分与结点编号 四边形四节点等参元matlab 程序 2013 年12 月28 日 — 2 — 附录： %---------------四边形四节点等参元 matlab 计算程序---------------------------- % 2013 年 % 13 级建筑与土木工程 Bru ce % 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.tx t','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):约束点号 四边形四节点等参元matlab 程序 2013 年12 月28 日 — 3 — %(n,2)与(n,3)分别为约束点x 方向和 y 方向的约束情况,受约束为 1 否则为 0 % 刚度矩阵的生成 %计算刚度矩阵，并对约束条件进行处理 Ke=zeros(2*NNODE,2*NNODE); % 单元刚度矩阵并清零 HK=zeros(2*NPOIN,2*NPOIN); % 张成总刚矩阵并清零 %调用子程序 生成单元刚度矩阵 for m=1:NELEM %m 为单元号 Ke=K(E,v,h,... COORD(LNODS(m,1),1),COORD(LNODS(m,1),2),... C...
