Mathematica Programs:

Frame Analysis Program

I made a program that analyses frames for the course of Computational Mechanics taught by Prof. Corigliano. Initially
the program only analyzed static cases, but I modified it after the course on structural dynamics to analyze structures
under dynamic loading using modal analysis. Shear deformation is accounted for and the user can impose fixed
displacements, elastic restraints, rigid offsets. The results show the internal stresses diagrams.
The benchmark below was presented during the Computational Mechanics course, the results were compared with
SAP2000:
The program shows the internal forces diagram as shown below

The user can select a certain member to read the exact values of the internal forces

Concerning the dynamic aspect of the program, the program was used for a project on the course of earthquake
design. A simple frame was analyzed, the obtained modal results were very close to the ones listed in the project (the
error was around 3%)

Mode Mode 3 Mode 4

4 8 12 16 4 8 12 16

13 23 14 24 15 25 16 13 23 14 24 15 25 16

3 7 11 15 3 7 11 15

9 20 10 21 11 22 12 9 20 10 21 11 22 12
Out[3992]= Mode Shape
2 6 10 14 2 6 10 14

5 17 6 18 7 19 8 5 17 6 18 7 19 8

1 5 9 13 1 5 9 13

1 2 3 4 1 2 3 4

Modal Frequency Hz 4.27 14.76

4 8 12 16 4 8 12 16

13 23 14 24 15 25 16 13 23 14 24 15 25 16

3 7 11 15 3 7 11 15

9 20 10 21 11 22 12 9 20 10 21 11 22 12
Out[3991]= Mode Shape
2 6 10 14 2 6 10 14

5 17 6 18 7 19 8 5 17 6 18 7 19 8

1 5 9 13 1 5 9 13

1 2 3 4 1 2 3 4

Modal Frequency Hz 30.38 49.14

The structure was then analyzed under the effect of an earthquake and the displacements were computed at each
story level.

4th Level
3rd Level
Out[9177]=
2nd Level
1st Level
Reinforced Concrete Multispan Beam Design Per the ACI Code:

This program was done for my Bachelor of Sciences graduation project. The user chooses the number of spans and
the dimension of the concrete beam along with the dead and live loads acting on the system and the program will
compute the ultimate shear and moment envelopes and compute the required steel reinforcement area and distribute
the steel.

Figure 1 Beam with the imposed loadings

Figure 3 Suggested steel distribution

Column Under Biaxial Bending

This application computes the moments and the axial load acting on a concrete cross section given the ortientation
and depth of the neutral axis. The user can choos the cross section properties (dimension, number of steel bars and
their diameters) and the program will compute the results.
The figure below shows the input on the right and the ouptut on the left.
 Neutral Axis
a
265.

30

Concrete Properties

Pn Mny Mnx F'c

Steel Properties
Out[1]= Number Of Bars In The X Direction
Number Of Bars In The Y Direction

Bar Diameter 28.7

Section Properties
Width
750
Depth
400
Concrete Cover
40

Another program using nearly the same code was created to compute the M-P surface for a concrete cross section,
given its properties. A similar one was done for uniaxial bending.
The figure below shows the M-P surface for a certain column cross section.
1000
200
400300

4000

2000

2000

0
200
400
The figure below shows a M-P surface in 2D. The cross section properties can be modified by the controllers on the
right.

Concrete Properties

C30 30 4 Bars f 'c25 MPa F'c

Steel Properties
Number Of Bars In The X Direction
Number Of Bars In The Y Direction
Ties Diameter 8 10 12

Section Properties
Width
300
Depth
300
Concrete Cover
40

Load
Axial Force
100
Applied Moment
20
Ties Diameter mm 8
Concrete Cover mm 40