Scribd
Upload
45 views

Assignment 1 Ce627: Structure Design Lab: Satyam Kumar PHD, Roll: 164046010 Indian Institute of Technology, Bombay

This document contains the input files for two structural analysis problems solved using a MATLAB program. Problem 1 defines a 4 node frame with 5 members, applied loads, and restraints. The MATLAB output shows the generated stiffness matrix and results of the analysis. Problem 2 defines a similar but different 4 node frame model, with different properties, loads, and restraints. The MATLAB program is then used to analyze this second frame model.

Uploaded by

satyam
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
45 views

Assignment 1 Ce627: Structure Design Lab: Satyam Kumar PHD, Roll: 164046010 Indian Institute of Technology, Bombay

This document contains the input files for two structural analysis problems solved using a MATLAB program. Problem 1 defines a 4 node frame with 5 members, applied loads, and restraints. The MATLAB output shows the generated stiffness matrix and results of the analysis. Problem 2 defines a similar but different 4 node frame model, with different properties, loads, and restraints. The MATLAB program is then used to analyze this second frame model.

Uploaded by

satyam
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 32

Assignment 1

CE627: Structure design lab

Submitted by

Satyam kumar

PhD , Roll: 164046010

Indian Institute of Technology, Bombay


INDIAN INSTITUTE OF TECHNOLOGY, BOMBAY
MUMBAI 400076

Figure 1: Problem 1 Frame

Problem 1: Matlab input file

% FILE: MDL_input

% Simple program for linear elastic structural analysis. The program


% 'MDL-main' requires user input of the model in this. Please follow the
% same input file format if entering a different structural model. This
% file is called from the main program and the file name serves as input
% into the main program.

% INPUT number of nodes and number of members

nnod = 4; %number of nodes


nmem = 5; %number of elements
nodedof = 3; %number of degrees of freedom/node (x-y-rot, plane frame)
%assumed to be 3 in many locations in program
fdof = nnod*nodedof %number of possible dofs (size of full KG matrix)

% INPUT node coordinates (node#, x-y coordinates and restraints x-y-rot


% for restraints: 0-free 1-fixed)
% node xcor ycor xres yres zres
ND = [ 1 0.0 8000.0 1 0 1;
2 6000.0 8000.0 0 0 1;
3 10000.0 0.0 1 1 1
4 0 0.0 1 1 1];

% INPUT member data (member #, node ni, node nj, release at ni, nj,
% Area, Modulus of Elasticiy, Moment of Inertia)
% member ni nj pin-i pin-j Area E Iz
MD = [1 4 1 1 1 4000 70 1 ;
2 1 2 1 1 4000 70 1 ;
3 3 2 1 1 4000 70 1 ;
4 3 1 1 1 4000 70 1 ;
5 4 2 1 1 4000 70 1 ];

% INPUT Load data at joints (node # correspond to those defined in ND)


% node Px Py M
PD = [ 1 0.0 -400.0 0.0 ;
2 800.0 -400.0 0.0 ;
3 0.0 0.0 0.0 ;
4 0.0 0.0 0.0 ];

% INPUT member load data (member #, constant distributed load, delta temp
% member numbers correspond to MD above)
% Memeber load cases typically need to be predefined withing the program
% and are straightforward to implement by defining fixed end forces.
% They are not implemented here to keep the program simple
% member Wm dT
MDL = [1 0 0;
2 0 0;
3 0 0;
4 0 0;
5 0 0]; % NOT IMPLEMENTED
Matlab output of Problem 1

fdof = 12
nmem = 5

Stiffness matrix:
member =1:

klocal =

35 0 0 -35 0 0
0 0 0 0 0 0
0 0 0 0 0 0
-35 0 0 35 0 0

0 0 0 0 0 0
0 0 0 0 0 0

kglobal =

0 0 0 0 0 0
0 35 0 0 -35 0
0 0 0 0 0 0
0 0 0 0 0 0
0 -35 0 0 35 0
0 0 0 0 0 0
id =

10 11 12 1 2 3

member =2

klocal =
Columns 1 through 6

46.6667 0 0 -46.6667 0 0
0 0 0 0 0 0
0 0 0 0 0 0
-46.6667 0 0 46.6667 0
0 0 0 0 0 0
0 0 0 0 0 0

kglobal =

Columns 1 through 6

46.6667 0 0 -46.6667 0 0

0 0 0 0 0 0
0 0 0 0 0 0
-46.6667 0 0 46.6667 0 0
0 0 0 0 0 0
0 0 0 0 0 0

id = 1 2 3 4 5 6
member =3
klocal =

Columns 1 through 6

31.3050 0 0 -31.3050 0 0

0 0 0 0 0 0

0 0 0 0 0 0
-31.3050 0 0 31.3050 0
0 0 0 0 0 0
0 0 0 0 0 0
kglobal =

Columns 1 through 6

6.2610 -12.5220 0 -6.2610 12.5220 0


-12.5220 25.0440 0 12.5220 -25.0440 0
0 0 0 0 0 0
-6.2610 12.5220 0 6.2610 -12.5220 0

12.5220 -25.0440 0 -12.5220 25.0440 0


0 0 0 0 0 0
id = 7 8 9 4 5 6

member = 4

klocal =
Columns 1 through 6
21.8643 0 0 -21.8643 00
0 0 0 0 0 0
0 0 0 0 0 0

-21.8643 0 0 21.8643 0 0
0 0 0 0 0 0
0 0 0 0 0 0

kglobal =

Columns 1 through 6
13.3319 -10.6655 0 -13.3319 10.6655 0
-10.6655 8.5324 0 10.6655 -8.5324 0
0 0 0 0 0 0
-13.3319 10.6655 0 13.3319 -10.6655 0
10.6655 -8.5324 0 -10.6655 8.5324 0

0 0 0 0 0

id = 7 8 9 1 2 3

member = 5
klocal =

Columns 1 through 6
28 0 0 -28 0 0
0 0 0 0 0 0
0 0 0 0 0 0
-28 0 0 28 0 0
0 0 0 0 0 0

0 0 0 0 0 0
kglobal =
Columns 1 through 6

10.0800 13.4400 0 -10.0800 -13.4400 0

13.4400 17.9200 0 -13.4400 -17.9200 0


0 0 0 0 0 0
-10.0800 -13.4400 0 10.0800 13.4400 0
-13.4400 -17.9200 0 13.4400 17.9200 0
0 0 0 0 0 0

id = 10 11 12 4 5 6

Load and Displacement :


P (Load)= 0
-400
0
800
-400

0
0
0
0
0
0

0
Uf (Displacement)= -9.1886
12.8365
-9.5844
Pr (Reaction)= -501.0366

0
0
-298.3858
479.1707
0
-0.5776

320.8293
0
Figure 2: Problem 2 Frame model

Problem 2: Matlab input file


% FILE: MDL_input

% Simple program for linear elastic structural analysis. The program


% 'MDL-main' requires user input of the model in this. Please follow the
% same input file format if entering a different structural model. This
% file is called from the main program and the file name serves as input
% into the main program.

% INPUT number of nodes and number of members

nnod = 4; %number of nodes


nmem = 5; %number of elements
nodedof = 3 %number of degrees of freedom/node (x-y-rot, plane frame)
%assumed to be 3 in many locations in program
fdof = nnod*nodedof %number of possible dofs (size of full KG matrix)

% INPUT node coordinates (node#, x-y coordinates and restraints x-y-rot


% for restraints: 0-free 1-fixed)
% node xcor ycor xres yres zres
ND = [ 1 0.0 2000.0 1 1 0;
2 5000 2000.0 0 0 0;
3 10000.0 2000 0 1 0;
4 5000 0.0 0 0 1];

% INPUT member data (member #, node ni, node nj, release at ni, nj,
% Area, Modulus of Elasticiy, Moment of Inertia)
% member ni nj pin-i pin-j Area E Iz
MD = [1 1 2 0 0 2000.0 200 200E6 ;
2 2 3 0 0 2000.0 200 200E6 ;
3 4 1 1 1 3000.0 200 1 ;
4 4 2 1 1 8000.0 200 1 ;
5 4 3 1 1 3000.0 200 1 ;];

% INPUT Load data at joints (node # correspond to those defined in ND)


% node Px Py M
PD = [ 1 0.0 0.0 0.0 ;
2 0.0 -100.0 0.0 ;
3 0.0 0.0 0.0 ;
4 0.0 0.0 0.0 ];

% INPUT member load data (member #, constant distributed load, delta temp
% member numbers correspond to MD above)
% Memeber load cases typically need to be predefined withing the program
% and are straightforward to implement by defining fixed end forces.
% They are not implemented here to keep the program simple
% member Wm dT
MDL = [1 0.0 0;
2 0.0 0;
3 0.0 0;
4 0.0 0;
5 0.0 0;]; % NOT IMPLEMENTED
Matlab output of Problem 2
nodedof =

fdof =

12

nmem =

member =

klocal =

1.0e+07 *

Columns 1 through 5

0.0000 0 0 -0.0000 0

0 0.0000 0.0010 0 -0.0000


0 0.0010 3.2000 0 -0.0010
-0.0000 0 0 0.0000 0
0 -0.0000 -0.0010 0 0.0000
0 0.0010 1.6000 0 -0.0010

Column 6

0
0.0010
1.6000

0
-0.0010
3.2000

kglobal =

1.0e+07 *

Columns 1 through 5

0.0000 0 0 -0.0000 0

0 0.0000 0.0010 0 -0.0000


0 0.0010 3.2000 0 -0.0010
-0.0000 0 0 0.0000 0
0 -0.0000 -0.0010 0 0.0000
0 0.0010 1.6000 0 -0.0010

Column 6
0
0.0010
1.6000

0
-0.0010
3.2000

id =

1 2 3 4 5 6

member =

klocal =

1.0e+07 *

Columns 1 through 5

0.0000 0 0 -0.0000 0
0 0.0000 0.0010 0 -0.0000
0 0.0010 3.2000 0 -0.0010

-0.0000 0 0 0.0000 0
0 -0.0000 -0.0010 0 0.0000
0 0.0010 1.6000 0 -0.0010

Column 6

0
0.0010
1.6000
0
-0.0010

3.2000

kglobal =

1.0e+07 *

Columns 1 through 5

0.0000 0 0 -0.0000 0
0 0.0000 0.0010 0 -0.0000
0 0.0010 3.2000 0 -0.0010

-0.0000 0 0 0.0000 0
0 -0.0000 -0.0010 0 0.0000
0 0.0010 1.6000 0 -0.0010

Column 6

0
0.0010
1.6000
0
-0.0010

3.2000

id =

4 5 6 7 8 9

member =

klocal =

Columns 1 through 5

111.4172 0 0 -111.4172 0

0 0 0 0 0
0 0 0 0 0
-111.4172 0 0 111.4172 0
0 0 0 0 0
0 0 0 0 0

Column 6
0
0
0

0
0
0

kglobal =

Columns 1 through 5

96.0493 -38.4197 0 -96.0493 38.4197


-38.4197 15.3679 0 38.4197 -15.3679
0 0 0 0 0

-96.0493 38.4197 0 96.0493 -38.4197


38.4197 -15.3679 0 -38.4197 15.3679
0 0 0 0 0

Column 6

0
0
0
0
0
0
id =

10 11 12 1 2 3

member =

klocal =

800 0 0 -800 0 0
0 0 0 0 0 0
0 0 0 0 0 0

-800 0 0 800 0 0
0 0 0 0 0 0
0 0 0 0 0 0

kglobal =

0 0 0 0 0 0
0 800 0 0 -800 0
0 0 0 0 0 0
0 0 0 0 0 0
0 -800 0 0 800 0

0 0 0 0 0 0
id =

10 11 12 4 5 6

member =

klocal =

Columns 1 through 5

111.4172 0 0 -111.4172 0
0 0 0 0 0
0 0 0 0 0
-111.4172 0 0 111.4172 0
0 0 0 0 0
0 0 0 0 0

Column 6

0
0
0

0
0
0

kglobal =

Columns 1 through 5

96.0493 38.4197 0 -96.0493 -38.4197


38.4197 15.3679 0 -38.4197 -15.3679

0 0 0 0 0
-96.0493 -38.4197 0 96.0493 38.4197
-38.4197 -15.3679 0 38.4197 15.3679
0 0 0 0 0

Column 6

0
0
0
0
0

id =

10 11 12 7 8 9
P=

0
0
0
-100
0
0

0
0
0
0
0

Uf =

-0.0019
-1.3708
-6.3909

0.0000
-2.7415
0.0019
-1.3708
-6.2812
Pr =

-0.0000
50.0000

50.0000
0
17/08/2017

Assignment 1

CE-627: Structural Design Lab

The displacement method is used by most of the structural analysis software programs to obtain
the displacements and reactions in a structure. Solve the following two problems.
Two structures are presented in the following problem. Solve these problems to obtain the
following:
1. Nodal displacements
2. Nodal forces
3. Member forces
Two methods of solutions should be used:
1. Analytical solutions using the matrix methods discussed in the class
2. Matlab implementation of the displacement method
Use appropriate assumptions where necessary. Students have the option to use the uploaded
Matlab program to solve these problems, or they can also write their own Matlab program.
Students should present a comparison of the results obtained using the analytical method and
Matlab implementation.

Problem 1 Problem 2 (use E = 200 GPa)


Final report should be submitted as a pdf file. The report should be professionally prepared.
Date of submission: 27th August, before 11:00 pm.

Notes: This assignment is for a credit of 10%.

1. Grading would be based on following:


Correctness of final results
Presentation and quality of the report
Technical understanding of the concepts
2. Any academic dishonesty would be penalized with zero marks and the person would
further be referred for disciplinary actions.

1
Scanned by CamScanner
Scanned by CamScanner
Scanned by CamScanner
Scanned by CamScanner
Scanned by CamScanner
Scanned by CamScanner
Scanned by CamScanner
Scanned by CamScanner
Comparison of result:

Table 1: Comparison of result for problem 1

Sr.No Reaction and displacement Analytical Matlab


1 D1 -9.195mm -9.1886mm
2 D2 12.836mm 12.8365mm
3 D3 -9.5858mm -9.5844mm
4 A4 -500.835KN -501.0366 KN
5 A7 -298.5858KN -298.3858 KN

6 A8 478.91104KN 479.1707 KN

7 A6 321.0867KN 320.8293KN

8 A5 -0.5537KN -0.5776KN

Table 2: Comparison of result for problem 2

Sr.No Reaction and displacement Analytical Matlab


1 D3 -0.0018156 -0.0019
2 D4 -1.2801mm -1.3708mm
3 D5 -6.051mm -6.3909mm
4 D7 -2.5602869mm -2.7415mm
5 D9 0.0018 0.0019
6 D10 -1.28mm -1.3708mm
7 D11 -5.94 -6.2812mm
8 A2 50KN 50KN
9 A8 50 KN 50KN

You might also like