2022（二）矩阵位移法大作业.docxVIP

矩阵位移法大作业 一、计算条件 请用矩阵位移法计算如下图示钢筋混凝土刚架的内力，并绘制内力图（弯矩图，剪力图，轴力图） 说明： 刚架尺寸，查阅计算条件表格； 粱的截面，采用矩形截面，确定方法：梁高按照1/16~1/10的跨度，选三跨梁段中，长度中间数取得梁高。梁的宽度一般为梁高度的1/3~1/2。为便于施工，统一模板模数，均采用50mm的倍数。 柱截面采用方形截面，边长等于粱高； 混凝土强度等级采用C30~C45，各自选择，参照混凝土教材附录选取弹性模量E。 二、计算： 采用MATLAB计算程序完成计算。 熟悉计算程序，读懂命令流，完成问答题。 输入数据：根据程序命令流读取数据的顺序，完成数据文件。 输出数据（将变量矩阵从MATLAB中拷贝到ecxel表中）： 每个单元的局部坐标系单刚矩阵、转换矩阵； 总刚矩阵； 总荷载列阵； 节点位移列阵； 最终内力矩阵FEE 三、回答问题（可以引用程序中的命令流说明） 计算程序如何计算单刚矩阵、转换矩阵、荷载列阵？ 单刚：EK=[E*A/L 0 0 -E*A/L 0 0; 0 12*E*I/L^3 6*E*I/L^2 0 -12*E*I/L^3 6*E*I/L^2; 0 6*E*I/L^2 4*E*I/L 0 -6*E*I/L^2 2*E*I/L; -E*A/L 0 0 E*A/L 0 0; 0 -12*E*I/L^3 -6*E*I/L^2 0 12*E*I/L^3 -6*E*I/L^2; 0 6*E*I/L^2 2*E*I/L 0 -6*E*I/L^2 4*E*I/L]; 转换：T=[ c s 0 0 0 0; -s c 0 0 0 0; 0 0 1 0 0 0; 0 0 0 c s 0; 0 0 0 -s c 0; 0 0 0 0 0 1]; 荷载：for i=1:NFPOIN N1=FPOIN(i,1); N1=N1*3-3; for j=1:3 FORCE(N1+j)=FORCE(N1+j)+FPOIN(i,j+1); 计算程序是如何形成总刚矩阵？ 总刚矩阵程序为 HK=zeros(3*NPOIN,3*NPOIN); for i=1:NELEM EK=ele_EK(i,LNODS,COORD,TY); T=zbzh(i,LNODS,COORD,TY); EKT=T*EK*T; for j=1:2 N1=LNODS(i,j)*3; for k=1:2 N2=LNODS(i,k)*3; HK((N1-2):N1,(N2-2):N2)=HK((N1-2):N1,(N2-2):N2)+EKT(j*3-2:j*3,k*3-2:k*3); 如何引入支承条件？如果支座有已知的非零位移，你认为应该怎样修改语句？ 修复引入相应的支承条件： FIXED=fscanf(FP1,%f,[1,NVFIX]); 修改：在对约束个数进行循环前先进行非零判断再往后运行for循环。 如何计算结点位移？ for i=1:NELEM for j=1:2 N1=LNODS(i,j)*3; EDISP(3*j-2:3*j)=DISP(N1-2:N1); end 如何计算各单元内力？ 程序可为 FEE=zeros(6,NELEM); EDISP=zeros(6,1); for i=1:NELEM for j=1:2 N1=LNODS(i,j)*3; EDISP(3*j-2:3*j)=DISP(N1-2:N1); end EK=ele_EK(i,LNODS,COORD,TY); T=zbzh(i,LNODS,COORD,TY); FE=EK*T*EDISP; for j=1:NFPRES if FPRES(j,1) == i F0=ele_FPRES(j,FPRES,LNODS,COORD); FE=FE+F0; end end FEE(1:6,i)=FE; End; 四、绘图要求 手绘计算简图（尽量按比例绘制），标明外荷载和具体尺寸；并在空白位置标明截面尺寸，混凝土等级等计算参数； 根据最终内力矩阵FEE，绘制结构弯矩、剪力、轴力图； 要求在A4白纸上绘图，图线、数字标注