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

四边形八节点等参元matlab程序 悬臂钢梁，尺寸如图一所示；v0.3。h1，E2.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 at short e FP1bjd.txt,rt; 打开数据文件 读入限制数据 EfscanfFP1,f,1; 弹性模量 vfscanfFP1,f,1; 泊松比 hfscanfFP1,f,1; 厚度 NELEMfscanfFP1,d,1; 单元数 NPOINfscanfFP1,d,1; 总结点数 NNODEfscanfFP1,d,1; 单元节点数 NFPOINfscanfFP1,d,1; 受力结点数 NVFIXfscanfFP1,d,1; 约束结点个数 LNODSfscanfFP1,f,[NNODE,NELEM]; 单元定义 单元结点号（逆时针） COORDfscanfFP1,f,[2,NPOIN]; 结点号 x,y坐标整体坐标下 FPOINfscanfFP1,f,[3,NFPOIN]; 节点力结点号、X方向力向右正，Y方向力向上正 FIXEDfscanfFP1,d,[3,NVFIX]; 约束信息数组n,3 n受约束节点数目, n,1约束点号 n,2与n,3分别为约束点x方向和y方向的约束状况,受约束为1否则为0 ******************************************************************* ******************************************************************* 平面应力问题的求解 ******************************************************************* ******************************************************************* 刚度矩阵的生成 计算刚度矩阵，并对约束条件进行处理 Kezeros2*NNODE,2*NNODE; 单元刚度矩阵并清零 HKzeros2*NPOIN,2*NPOIN; 张成总刚矩阵并清零 调用子程序 生成单元刚度矩阵 for m1NELEM m为单元号 KeKE,v,h,... COORDLNODSm,1,1,COORDLNODSm,1,2,... COORDLNODSm,3,1,COORDLNODSm,3,2,... COORDLNODSm,5,1,COORDLNODSm,5,2,... COORDLNODSm,7,1,COORDLNODSm,7,2; 调用单元刚度矩阵 aLNODSm,; 临时向量,用来记录当前单元的节点编号 对总刚度矩阵的处理 for j18 for k18 HKaj*2-1aj*2,ak*2-1ak*2HKaj*2-1aj*2,ak*2-1ak*2... Kej*2-1j*2,k*2-1k*2; end end end 对荷载向量进行处理 FORCEzeros2*NPOIN,1; 张成总荷载向量并清零 for i1NFPOIN b1FPOINi,1*2-1;b2FPOIN