第六章 平面刚架程序设计.pptVIP

编码（输入初始控制数据） 结点数NPOIN =4 单元数NELEM =3 约束数NVFIX =5 非零位移数NFLWY =0 结点力数NFPOIN =2 非结点力数NFPRES =2 弹性支撑数NTXZZ =0 弹性模量YOUNG =1 输入初始数据srsj 输入单元连接信息（单元属性） LNODS=fscanf(FP1,%f,[4,NELEM]) 左、右结点号，面积，惯性矩 1, 2, 1.0, 1.0 2, 4, 1.0, 1.0 3, 4, 1.0, 1.0 输入结点坐标 COORD=fscanf(FP1,%f,[2,NPOIN]) x坐标、y坐标 0.0, 0.0 0.0, 2.0 4.0, 0.0 4.0, 2.0 输入初始数据srsj 输入结点荷载信息 FPOIN=fscanf(FP1,%f,[4,NFPOIN])‘ 结点号 X方向力（向右为正） Y方向力（向上为正） M力偶（逆时针为正） 2, 50.0, 0.0, 0.0 4, 0.0, 0.0, -50.0 输入初始数据srsj 输入非结点荷载信息 FPRES=fscanf(FP1,%f,[4,NFPRES])‘ 单元号、荷载类型 荷载大小、距离左端长度 2, 2, -100.0, 2.0 3, 1, -10.0, 0.0 输入初始数据srsj 输入约束信息 FIXED=fscanf(FP1,‘%f’,NVFIX)‘ 约束对应的位移编码 1, 2, 7, 8, 9 输入非零已知位移信息 FLWY=fscanf(FP1,%f,[2,NFLWY]) 非零位移的位移编码和位移值 输入弹性支撑信息 TXZZ=fscanf(FP1,%f,[2,NTXZZ]) 弹性支撑的位移编码和弹性刚度值 计算原始总刚矩阵EK 步骤： 求每一单元、局部坐标系下的、单元刚度矩阵EK 构造函数：ele_EK.m 求每一单元、坐标转换矩阵T 构造函数：zbzh.m 求每一单元、整体坐标系下的、单元刚度矩阵EKT EKT=T’*EK*T； 对号入座，集装总刚HK 计算局部坐标系单元刚度矩阵EKele_EK.M：理论 计算局部坐标系单元刚度矩阵EKele_EK.M：程序 function EK=ele_EK (i, LNODS, COORD,E) NL=LNODS(i,1); %左结点号 NR=LNODS(i,2); %右结点号 dx=COORD(NR,1)-COORD(NL,1); % x 坐标差 dy=COORD(NR,2)-COORD(NL,2); % y 坐标差 L=sqrt(dx^2+dy^2); % 单元长度 c=dx/L; % cos a (与 x 轴夹角余弦） s=dy/L; % sin a (与 x 轴夹角正弦） A=LNODS(i,3); %面积 I=LNODS(i,4); %惯性矩 求每一单元、坐标转换矩阵Tzbzh.m：理论 求每一单元、坐标转换矩阵Tzbzh.m：程序 function T=zbzh(i,LNODS,COORD) NL=LNODS(i,1); %左结点号 NR=LNODS(i,2); %右结点号 dx=COORD(NR,1)-COORD(NL,1); % x 坐标差 dy=COORD(NR,2)-COORD(NL,2);