机械仿真基于matlab二维钢桁架分析与设计【含Matlab源码 2103期】.zip

clear; clc; close all; %输入 %-------------------------------------------------------------------------- S=24; %跨度 N=12; %托架数（偶数） H1=1; %边缘高度 H2=3.7; %跨中高度 DL=-4000; %上弦静载荷 LL=-4000; %上弦活载 %节点坐标和自由度 %------------------------------------------ for i=1:1:N+1 Node{i}.x=S*(i-1)/N; Node{i}.y=0; Node{i}.DOFx=2*i-1; Node{i}.DOFy=2*i; end for i=N+2:1:1.5*N+2 Node{i}.x=S*(i-N-2)/N; Node{i}.y=H1+(H2-H1)*(i-N-2)/(N/2); Node{i}.DOFx=2*i-1; Node{i}.DOFy=2*i; end for i=1.5*N+3:1:2*N+2 Node{i}.x=S*(i-N-2)/N; Node{i}.y=H1+(H2-H1)*(2*N+2-i)/(N/2); Node{i}.DOFx=2*i-1; Node{i}.DOFy=2*i; end %单元连通性、截面、弹性模量和屈服强度 %-------------------------------------------------------------------------- for i=1:1:N %下弦 Element{i}.Nodes=[i i+1]; Element{i}.A=1472e-6; Element{i}.rmin=0.0237; Element{i}.E=200e9; Element{i}.Fy=250e6; end for i=1:1:N %上和弦 Element{N+i}.Nodes=[N+1+i N+2+i]; Element{N+i}.A=1472e-6; Element{N+i}.rmin=0.0237; Element{N+i}.E=200e9; Element{N+i}.Fy=250e6; end for i=1:1:N/2 %对角线成员 Element{2*N+i}.Nodes=[i+1 N+1+i]; Element{2*N+i}.A=1472e-6; Element{2*N+i}.rmin=0.0237; Element{2*N+i}.E=200e9; Element{2*N+i}.Fy=250e6; end for i=N/2+1:1:N %对角线成员 Element{2*N+i}.Nodes=[i N+2+i]; Element{2*N+i}.A=1472e-6; Element{2*N+i}.rmin=0.0237; Element{2*N+i}.E=200e9; Element{2*N+i}.Fy=250e6; end for i=1:1:N+1 %垂直成员 Element{3*N+i}.Nodes=[i N+1+i]; Element{3*N+i}.A=1472e-6; Element{3*N+i}.rmin=0.0237; Element{3*N+i}.E=200e9; Element{3*N+i}.Fy=250e6; end %支撑、支撑位移和作用力 %------------------------------------------------- Support=[1 2 2*N+2]; U_s=[0 0 0]'; Force=zeros(4*N+4,1); Force(2*N+4)=(DL+LL)/2; Force(4*N+4)=(DL+LL)/2; Force(2*N+6:2:4*N+2)=DL+LL; %加工 %-------------------------------------------------------------------------- [Node,Element,F_s]=Analysis(Node,Element,Support,U_s,Force); %检查 AISC 代码 %-------------------------------------------------------------------------- Element=Check_Code(Element); %输出 %-------------------------------------------------------------------------- Mag_Factor=20; Plot_Results(Node,Element,Mag_Factor,0); Plot_Results(Node,Element,Mag_Factor,1); title('拉紧'); Plot_Results(Node,Element,Mag_Factor,2); title('压力'); Plot_Results(Node,Element,Mag_Factor,3); title('轴向力'); Plot_Results(Node,Element,Mag_Factor,4); title('设计比例');