悬臂钢梁，尺寸如图一所示；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 %******************************************************************* %******************************************************************* %========...
