Department of Civil Engineering: Iit Delhi
Department of Civil Engineering: Iit Delhi
IIT DELHI
CEL331 STRUCTURAL ANALYSIS II (2012-13)
PROGRAMMING ASSIGNMENT (A)
Develop a computer program (in C++, FORTRAN, MATLAB or any standard language) to analyse a plane
frame of the type shown in the figure below. The frame has ‘m’ number of bays and ‘n’ number of storeys
and is subjected to vertical loads only. The program should interactively obtain following data from the
user:
The program should automatically generate joint and member numbers and assemble the total structural
stiffness matrix. Those programming in MATLAB should not use existing function for matrix inversion,
but should develop their own. Final output should be in the following form, preferably in the form of an
output file:
Test your program on a test frame consisting of m =5 and n = 3. Check for the equilibrium of external
forces and reactions. Assume columns of 300x300mm and beams 300x600mm cross section. A storey
height of 4m may be considered. Assume M 20 concrete.
SUBMISSION: One report per group should be submitted. It should be typed and should include brief
theory, basic steps executed by the program (in the form of a flow chart) and results of the ‘test’ as
indicated above. The report should be less than 8 pages under all circumstances. Along with report, each
group should submit the source code as well as executable file in a floppy or CD. Last date for submission
is 01 NOVEMBER 2012. Negative marking will be imposed on groups submitting late.
‘n’ storeys
of equal
height
Develop a computer program (in C++, FORTRAN, MATLAB or any standard language) to analyse a plane
frame of the type shown in the figure below. The frame has ‘m’ number of bays and ‘n’ number of storeys
and is subjected to horizontal loads only. The program should interactively obtain following data from the
user:
The program should automatically generate joint and member numbers and assemble the total structural
stiffness matrix. Those programming in MATLAB should not use existing function for matrix inversion,
but should develop their own. Final output should be in the following form, preferably in the form of an
output file:
Test your program on a test frame consisting of m =5 and n = 3. Check for the equilibrium of external
forces and reactions. Assume columns of 300x300mm and beams 300x600mm cross section. A storey
height of 4m may be considered. Assume M 20 concrete.
SUBMISSION: One report per group should be submitted. It should be typed and should include brief
theory, basic steps executed by the program (in the form of a flow chart) and results of the ‘test’ as
indicated above. The report should be less than 8 pages under all circumstances. Along with report, each
group should submit the source code as well as executable file in a floppy or CD. Last date for submission
is 01 NOVEMBER 2012. Negative marking will be imposed on groups submitting late.
‘n’ storeys
of equal
height
2|Page
© Dr. Suresh Bhalla, Department of Civil Engineering, IIT Delhi