bài tập tin học ứng dụng

Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu I Phần 1: Bài tập matlab 1.Bài tập a Đoạn chương trình matlab njoint=9;% Number of nodes (include fulcrum) %Node number xcord= [0 0 6 10 10 10];% x-cordinate of node n (m) ycord= [0 8 8];% y-cordinate of node n (m) xsupport= [1 0 0 0];% Support condition 0: free; 1: constraint ysupport= [1 0 0 0];% Support condition 0: free; 1: constraint zsupport= [1 0 0 0];% Support condition 0: free; 1: constraint xload= 1e3*[0 20 20 0 0 0];%Concentrated force at node in x direction (N) yload= [0 0 0 0 0];%Concentrated force at node in y direction (N) zload= [0 0 0 0 0];%Moment at node around z direction (Nm) % Information of element nmemb=10; % Number of elements % Element number 10 joint1= [1 6];% The number of node that the element m start from joint2= [2 9];% The number of node that the element m end to E= 2e11*[1 1 1 1 1 1];% Young modulus of material of element m A= 0.25*[1 1 1 1 1 1];% Area section of element m (m2) I= 0.5^4/12*[1 1 1 1 1 1];% Moment of inertia of element m (m4) % -caculation -skk = zeros(3*njoint,3*njoint);% The stiffness matrix of the whole structure displ = zeros(3*njoint,1);% The displacement vector load = zeros(3*njoint,1);% The force vector support = zeros(3*njoint,1);% The support condition axial1 = zeros(nmemb,1);% The axial force at start point of elements axial2 = zeros(nmemb,1);% The axial force at end point of elements shear1 = zeros(nmemb,1);% The shear force at start point of elements SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu shear2 = zeros(nmemb,1);% The shear force at end point of elements moment1 = zeros(nmemb,1);% The moment at start point of elements moment2 = zeros(nmemb,1);% The moment at end point of elements % Prepare for calculaton for m=1:njoint load(3*m-2:3*m) =[xload(m); yload(m); zload(m)];% Define the force vector support(3*m-2:3*m)=[xsupport(m);ysupport(m);zsupport(m)];% Define the support condition end for m=1:nmemb i=joint1(m);% Node number j=joint2(m);% Node number delx=xcord(j)-xcord(i); dely=ycord(j)-ycord(i); L(m)=sqrt(delx^2+dely^2);% Length of element m cs(m)=delx/L(m);% cosin inclination of element m sn(m)=dely/L(m);% sin inclination of element m end % Define the stiffness matrix of whole structure skk for m=1:nmemb al=L(m); ei=E(m)*I(m); ea=E(m)*A(m); % Define the stiffness of element kff (3X3) k22=[ea/al 0;0 12*ei/al^3 -6*ei/al^2;0 -6*ei/al^2 4*ei/al]; % Define the transfomation matrix T(6x6) T(1:3,1:3)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; T(4:6,4:6)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; % Define the balance Matrix H (3x3) H=[1 0;0 0;0 al 1]; % Define conection matrix [C] (6x3) C=[-H;eye(3)]; % Calculate the element stiffness matrix (6x6) sk=T*C*k22*C'*T'; % The degree-of-freedom number of both edge nodes of element m i=3*joint1(m);% The degre-of-freedom of node j=3*joint2(m);% The degre-of-freedom of node SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu n=[i-2 i-1 i j-2 j-1 j]; % Input the element stifness matrix to structure stifness matrix skk(n,n)=skk(n,n)+sk; end % Define the reaction force and displacement indf=find(support==0);% Find the index of unknown displacement inds=find(support);% Find the index of known displacement % Define part matrix of rearrangement structure stiffess matrix kff=skk(indf,indf); kfs=skk(indf,inds); ksf=skk(inds,indf); kss=skk(inds,inds); % Define unknown displacement displ(indf)= inv(kff)*(load(indf)-kfs*displ(inds)) % Define reaction force load(inds)=ksf*displ(indf)+kss*displ(inds) % Define the element stress for m=1:nmemb al=L(m); ei=E(m)*I(m); ea=E(m)*A(m); % Define the stiffness of element kff (3X3) k22=[ea/al 0;0 12*ei/al^3 -6*ei/al^2;0 -6*ei/al^2 4*ei/al]; % Define the transfomation matrix T(6x6) T(1:3,1:3)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; T(4:6,4:6)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; % Define the balance Matrix H (3x3) H=[1 0;0 0;0 al 1]; % Define conection matrix [C] (6x3) C=[-H;eye(3)]; % The degree-of-freedom number of both edge nodes of element m i=3*joint1(m);% The degre-of-freedom of node j=3*joint2(m);% The degre-of-freedom of node n=[i-2 i-1 i j-2 j-1 j]; force=C*k22*C'*T'*displ(n); axial1(m)=force(1); SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu axial2(m)=force(4); shear1(m)=force(2); shear2(m)=force(5); moment1(m)=force(3); moment2(m)=force(6); end axial1 axial2 shear1 shear2 moment1 moment2 b Kết xuất từ chương trình displ = 1.0e-003 * 0 0.1102 0.0008 -0.0248 0.2064 0.0011 -0.0127 0 0.1087 0.0008 -0.0150 0.2045 0.0011 -0.0074 0 0.1082 -0.0016 -0.0214 0.2040 -0.0022 -0.0114 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu load = 1.0e+004 * -1.1852 -1.0389 3.0151 2.0000 0 2.0000 0 -1.5383 -1.0047 3.4668 0 0 0 -1.2765 2.0436 3.1106 0 0 0 axial1 = 1.0e+004 * -1.0389 -0.3490 -1.0047 -0.3205 2.0436 0.6696 1.2282 0.6858 1.5866 0.5907 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu axial2 = 1.0e+004 * 1.0389 0.3490 1.0047 0.3205 -2.0436 -0.6696 -1.2282 -0.6858 -1.5866 -0.5907 shear1 = 1.0e+004 * 1.1852 0.4134 1.5383 0.9959 1.2765 0.5907 -0.6898 -1.3740 -0.3490 -0.6696 shear2 = 1.0e+004 * -1.1852 -0.4134 -1.5383 -0.9959 -1.2765 -0.5907 0.6898 1.3740 0.3490 0.6696 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu moment1 = 1.0e+004 * 3.0151 0.5137 3.4668 1.7939 3.1106 0.9199 -2.2392 -2.5806 -1.1400 -1.2354 moment2 = 1.0e+004 * 1.7255 1.1400 2.6864 2.1896 1.9955 1.4429 -1.8998 -2.9154 -0.9543 -1.4429 Đơn vị : N,m 20kN 10 6696 N 12282 6858 20436 3205 15866 10047 10389 20kN 3490 5907 4 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang Bài tập tin học ứng dụng 9543 11400 12354 17255 19955 9199 34668 31106 25806 22392 26864 17939 13740 M 12765 6898 15383 11852 5137 18998 Q 29154 5907 9959 4134 6696 14429 21896 11400 3490 14429 GVHD: Th.S Bùi Quang Hiếu 30151 11852 15383 30151 34668 10389 10047 12765 31106 20436 Hình I.1: Sơ đồ kết cấu biểu đồ nội lực, phản lực 2.Bài tập a Đoạn chương trình matlab njoint=8;% Number of nodes % Node number xcord= [0 0 (m) ycord= [0 (m) xsupport= [1 0 free; 1: constraint ysupport= [1 0 free; 1: constraint zsupport= [1 0 free; 1: constraint xload= 1e3*[0 20 20 node in x direction (N) yload= [0 0 node in y direction (N) (include fulcrum) 6 10 10];% x-cordinate of node n 8];% y-cordinate of node n 0 0];% Support condition 0: 0 0];% Support condition 0: 0 0];% Support condition 0: 0 0];%Concentrated force at 0 0];%Concentrated force at SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu zload= [0 0 0 0 0];%Moment at node around z direction (Nm) % Information of element nmemb=8;% Number of element % Element number joint1= [1 6];% The number of node that the from element m start from joint2= [2 8];% The number of node that the from element m end to E= 2e11*[1 1 1 1 1];% Young modulus of material of element m A= 0.25*[1 1 1 1 1];% Area section of element m (m2) I= 0.5^4/12*[1 1 1 1 1];% Moment of inertia of element m (m4) % -Calculation skk = zeros(3*njoint,3*njoint);% The stiffness matrix of the whole structure displ = zeros(3*njoint,1);% The displacement vector load = zeros(3*njoint,1);% The force vector support = zeros(3*njoint,1);% The support condition axial1 = elements axial2 = elements shear1 = elements shear2 = elements moment1 = elements moment2 = zeros(nmemb,1);% The axial force at start point of zeros(nmemb,1);% The axial force at end point of zeros(nmemb,1);% The shear force at start point of zeros(nmemb,1);% The shear force at end point of zeros(nmemb,1);% The moment at start point of zeros(nmemb,1);% The moment at end point of elements % Prepare for calculaton for m=1:njoint load(3*m-2:3*m) =[xload(m); yload(m); zload(m)];% Define the force vector support(3*m-2:3*m)=[xsupport(m);ysupport(m);zsupport(m)];% Define the support condition end for m=1:nmemb i=joint1(m);% Node number j=joint2(m);% Node number delx=xcord(j)-xcord(i); dely=ycord(j)-ycord(i); L(m)=sqrt(delx^2+dely^2);% Length of element m cs(m)=delx/L(m);% cosin inclination of element m SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu sn(m)=dely/L(m);% sin inclination of element m end % Define the stiffness matrix of whole structure skk for m=1:nmemb al=L(m); ei=E(m)*I(m); ea=E(m)*A(m); % Define the stiffness of element kff (3X3) k22=[ea/al 0;0 12*ei/al^3 -6*ei/al^2;0 -6*ei/al^2 4*ei/al]; % Define the transfomation matrix T(6x6) T(1:3,1:3)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; T(4:6,4:6)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; % Define the balance Matrix H (3x3) H=[1 0;0 0;0 al 1]; % Define conection matrix [C] (6x3) C=[-H;eye(3)]; % Calculate the element stiffness matrix (6x6) sk=T*C*k22*C'*T'; % The degree-of-freedom number of both edge nodes of element m i=3*joint1(m);% The degre-of-freedom of node j=3*joint2(m);% The degre-of-freedom of node n=[i-2 i-1 i j-2 j-1 j]; % Input the element stifness matrix to structure stifness matrix skk(n,n)=skk(n,n)+sk; end % Define the reaction force and displacement indf=find(support==0);% Find the index of unknown displacement inds=find(support);% Find the index of known displacement % Define part matrix of rearrangement structure stiffess matrix kff=skk(indf,indf); kfs=skk(indf,inds); ksf=skk(inds,indf); kss=skk(inds,inds); % Define unknown displacement displ(indf)= inv(kff)*(load(indf)-kfs*displ(inds)) % Define reaction force load(inds)=ksf*displ(indf)+kss*displ(inds) SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 10 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu 0 axial1 = 1.0e+004 * -1.5867 -0.4246 0.7284 -0.4337 0.8583 0.7703 1.4474 0.5405 axial2 = 1.0e+004 * 1.5867 0.4246 -0.7284 0.4337 -0.8583 -0.7703 -1.4474 -0.5405 shear1 = 1.0e+004 * 1.7823 0.5526 1.6772 0.9069 0.5405 -1.1621 -0.4246 -0.8583 shear2 = 1.0e+004 * -1.7823 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 13 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu -0.5526 -1.6772 -0.9069 -0.5405 1.1621 0.4246 0.8583 moment1 = 1.0e+004 * 4.4152 0.7335 4.2630 1.0792 2.3684 -3.4476 -1.4770 -1.4777 moment2 = 1.0e+004 * 2.7141 1.4770 2.4458 2.5483 1.9554 -3.5250 -1.0706 -1.9554 Đơn vị : N,m 20kN 8 6 7703 7284 15867 N 4337 14474 8583 20kN 4246 5405 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 14 Bài tập tin học ứng dụng 10706 35250 27141 24458 M 10792 16772 5405 Q 7335 19554 14777 8583 34476 9069 5526 17823 11621 25483 14770 14770 4246 19554 GVHD: Th.S Bùi Quang Hiếu 44152 16772 17823 44152 15867 42630 5405 42630 7284 23684 23684 8583 Hình I.2: Sơ đồ kết cấu biểu đồ nội lực, phản lực 3.Bài tập a Đoạn chương trình matlab njoint=9;% Number of nodes (include fulcrum) %Node number xcord= [0 0 6 10 10 10];% x-cordinate of node n (m) ycord= [0 8 8];% y-cordinate of node n (m) xsupport= [1 0 0 0];% Support condition 0: free; 1: constraint ysupport= [1 0 0 0];% Support condition 0: free; 1: constraint zsupport= [0 0 0 0 0];% Support condition 0: free; 1: constraint xload= 1e3*[0 20 20 0 0 0];%Concentrated force at node in x direction (N) yload= [0 0 0 0 0];%Concentrated force at node in y direction (N) zload= [0 0 0 0 0];%Moment at node around z direction (Nm) SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 15 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu % Information of element nmemb=10; % Number of elements % Element number 10 joint1= [1 6];% The number of node that the element m start from joint2= [2 9];% The number of node that the element m end to E= 2e11*[1 1 1 1 1 1];% Young modulus of material of element m A= 0.25*[1 1 1 1 1 1];% Area section of element m (m2) I= 0.5^4/12*[1 1 1 1 1 1];% Moment of inertia of element m (m4) % -caculation -skk = zeros(3*njoint,3*njoint);% The stiffness matrix of the whole structure displ = zeros(3*njoint,1);% The displacement vector load = zeros(3*njoint,1);% The force vector support = zeros(3*njoint,1);% The support condition axial1 = elements axial2 = elements shear1 = elements shear2 = elements moment1 = elements moment2 = zeros(nmemb,1);% The axial force at start point of zeros(nmemb,1);% The axial force at end point of zeros(nmemb,1);% The shear force at start point of zeros(nmemb,1);% The shear force at end point of zeros(nmemb,1);% The moment at start point of zeros(nmemb,1);% The moment at end point of elements % Prepare for calculaton for m=1:njoint load(3*m-2:3*m) =[xload(m); yload(m); zload(m)];% Define the force vector support(3*m-2:3*m)=[xsupport(m);ysupport(m);zsupport(m)];% Define the support condition end for m=1:nmemb i=joint1(m);% Node number j=joint2(m);% Node number delx=xcord(j)-xcord(i); dely=ycord(j)-ycord(i); L(m)=sqrt(delx^2+dely^2);% Length of element m cs(m)=delx/L(m);% cosin inclination of element m sn(m)=dely/L(m);% sin inclination of element m end SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 16 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu % Define the stiffness matrix of whole structure skk for m=1:nmemb al=L(m); ei=E(m)*I(m); ea=E(m)*A(m); % Define the stiffness of element kff (3X3) k22=[ea/al 0;0 12*ei/al^3 -6*ei/al^2;0 -6*ei/al^2 4*ei/al]; % Define the transfomation matrix T(6x6) T(1:3,1:3)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; T(4:6,4:6)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; % Define the balance Matrix H (3x3) H=[1 0;0 0;0 al 1]; % Define conection matrix [C] (6x3) C=[-H;eye(3)]; % Calculate the element stiffness matrix (6x6) sk=T*C*k22*C'*T'; % The degree-of-freedom number of both edge nodes of element m i=3*joint1(m);% The degre-of-freedom of node j=3*joint2(m);% The degre-of-freedom of node n=[i-2 i-1 i j-2 j-1 j]; % Input the element stifness matrix to structure stifness matrix skk(n,n)=skk(n,n)+sk; end % Define the reaction force and displacement indf=find(support==0);% Find the index of unknown displacement inds=find(support);% Find the index of known displacement % Define part matrix of rearrangement structure stiffess matrix kff=skk(indf,indf); kfs=skk(indf,inds); ksf=skk(inds,indf); kss=skk(inds,inds); % Define unknown displacement displ(indf)= inv(kff)*(load(indf)-kfs*displ(inds)) % Define reaction force load(inds)=ksf*displ(indf)+kss*displ(inds) % Define the element stress SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 17 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu for m=1:nmemb al=L(m); ei=E(m)*I(m); ea=E(m)*A(m); % Define the stiffness of element kff (3X3) k22=[ea/al 0;0 12*ei/al^3 -6*ei/al^2;0 -6*ei/al^2 4*ei/al]; % Define the transfomation matrix T(6x6) T(1:3,1:3)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; T(4:6,4:6)=[cs(m) -sn(m) 0; sn(m) cs(m) 0;0 1]; % Define the balance Matrix H (3x3) H=[1 0;0 0;0 al 1]; % Define conection matrix [C] (6x3) C=[-H;eye(3)]; % The degree-of-freedom number of both edge nodes of element m i=3*joint1(m);% The degre-of-freedom of node j=3*joint2(m);% The degre-of-freedom of node n=[i-2 i-1 i j-2 j-1 j]; force=C*k22*C'*T'*displ(n); axial1(m)=force(1); axial2(m)=force(4); shear1(m)=force(2); shear2(m)=force(5); moment1(m)=force(3); moment2(m)=force(6); end axial1 axial2 shear1 shear2 moment1 moment2 b kết xuất từ chương trình displ = 1.0e-003 * 0 -0.1366 0.4281 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 18 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu 0.0014 -0.0479 0.5654 0.0017 -0.0144 0 -0.1467 0.4267 0.0013 -0.0267 0.5633 0.0017 -0.0105 0 -0.1393 0.4262 -0.0027 -0.0409 0.5629 -0.0034 -0.0137 load = 1.0e+004 * -1.1553 -1.7272 2.0000 0 2.0000 0 -1.5631 -1.6819 0 0 0 -1.2817 3.4091 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 19 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu 0 0 0 axial1 = 1.0e+004 * -1.7272 -0.4336 -1.6819 -0.4154 3.4091 0.8490 1.0939 0.7452 1.7509 0.5364 axial2 = 1.0e+004 * 1.7272 0.4336 1.6819 0.4154 -3.4091 -0.8490 -1.0939 -0.7452 -1.7509 -0.5364 shear1 = 1.0e+004 * 1.1553 0.2491 1.5631 1.2144 1.2817 0.5364 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 20 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu -1.2936 -2.5601 -0.4336 -0.8490 shear2 = 1.0e+004 * -1.1553 -0.2491 -1.5631 -1.2144 -1.2817 -0.5364 1.2936 2.5601 0.4336 0.8490 moment1 = 1.0e+004 * -0.0000 -0.3720 -0.0000 2.0093 -0.0000 0.3649 -4.2490 -4.7489 -1.3686 -1.6152 moment2 = 1.0e+004 * 4.6211 1.3686 6.2523 2.8484 5.1267 1.7809 -3.5127 -5.4915 -1.2332 SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 21 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu -1.7809 Đơn vị : N,m 20kN 10 4154 4336 10939 7452 62523 47489 42490 20093 M 51267 3649 12817 15631 11553 11553 17272 54915 35127 5364 Q 3720 46211 25601 17809 16152 8490 12144 2491 28484 13686 4336 17809 12332 13686 12936 N 34091 16819 17272 20kN 4 8490 5364 17509 15631 16819 12817 34091 Hình I.3: Sơ đồ kết cấu biểu đồ nội lực, phản lực SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 22 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu II Phần 2: 1.Bài tập sap a.Khai báo vật liệu SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 23 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu b.Khai báo tiết diện SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 24 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu c.Sơ đồ kết cấu biểu đồ nội lực sap Đơn vị: N,m -6695.57 3205.17 Q 12765.20 15383.06 5136.86 17255.50 25806.18 13739.99 17939.43 30151.47 11851.74 30151.47 10388.71 15383.06 34667.98 10046.85 14428.58 21896.21 12354 18997.52 11399.86 14428.58 29153.79 6695.57 6898.31 11851.74 9543 9958.91 4134.18 3490.40 -12282.44 10046.85 N -6858.29 10388.71 3490.40 -15865.82 11399.86 22392.36 -5906.91 -20436.56 10 20kN 5906.91 20kN 26864.27 19954.73 9199.06 34667.98 M 31106.05 12765.20 31106.05 20436.56 Hình II.1: Sơ đồ kết cấu, biểu đồ nội lực, phản lực tập SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 25 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu Đơn vị: N,m 20kN 8 9068.85 14770.07 25483.33 5404.78 Q 7335.39 27140.51 44152.28 15867.00 16772.02 42630.35 7284.11 19554.46 24457.73 M 10792.07 44152.28 17823.20 19554.46 10706.24 8582.29 16772.02 11620.95 17823.20 5526.36 4246.05 35249.79 14777.09 -7703.17 15867.00 14770.07 -7284.11 N -8582.89 4246.05 34475.90 -14473.64 20kN 4336.84 -5404.78 42630.35 23683.82 5404.78 23684.82 8582.29 Hình II.2: Sơ đồ kết cấu, biểu đồ nội lực, phản lực tập SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 26 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu Đơn vị:N,m -10938.65 -8490.21 12332.08 17808.93 54915.43 8490.21 13685.67 25601.14 12936.14 Q 11552.67 17272.43 15630.65 16818.92 28484.01 62522.61 20093.08 17808.93 M 51266.71 3648.71 47489.12 42490.27 12816.68 15630.65 3720.41 46210.68 11552.57 16151.92 35126.57 2491.32 -7452.27 12144.27 4336.29 N -34091.35 -5364.41 16818.92 4153.92 17272.43 4336.29 -17508.68 13685.67 4 20kN 10 5364.41 20kN 12816.68 34091.35 Hình II.3: Sơ đồ kết cấu, biểu đồ nội lực, phản lực tập SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 27 ... Hiệp-Lớp:10X1B-Nhóm:68A Trang 22 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu II Phần 2: 1 .Bài tập sap a.Khai báo vật liệu SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 23 Bài tập tin học ứng dụng GVHD: Th.S... 20436.56 Hình II.1: Sơ đồ kết cấu, biểu đồ nội lực, phản lực tập SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 25 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu Đơn vị: N,m 20kN 8 9068.85... 8582.29 Hình II.2: Sơ đồ kết cấu, biểu đồ nội lực, phản lực tập SVTH: Nguyễn Thái Hiệp-Lớp:10X1B-Nhóm:68A Trang 26 Bài tập tin học ứng dụng GVHD: Th.S Bùi Quang Hiếu Đơn vị:N,m -10938.65 -8490.21