Excel For Civil Engineers Preview PDF

PREVIEW ONLY
ANI
A INTRODDUC
CTIOONT
TO:
E CELFO
EXC ORC
CIVILEN
NGINEE
ERS

by:Gu
untharPangarribuan
FilesIInclud. t TranslaationR
Direct Revision
n1.1
20042007
This book is intended to introduce a beginner level of using Microsoft Excel
in civil engineering practices

Translated and re-produced in Nov, 2015.
An Introduction to: Excel for Civil Engineers 2

http://exceljunior.blogspot.co.id
TABLE OF CONTENTS
CHAPTER 1: BASICS OF EXCEL ................................................................................. 6

1.1 Worksheet and Workbook .................................................................................... 6
1.2 Data Type ............................................................................................................. 7
1.3 Formula .............................................................................................................. 10
1.4 Built-In Function .................................................................................................. 12
1.5 Array Formula ..................................................................................................... 14
1.6 Data Formatting .................................................................................................. 16
1.7 Error Message .................................................................................................... 17
1.8 Printing ............................................................................................................... 17
1.9 Making Charts .................................................................................................... 18
1.10 Engineering Drawing .......................................................................................... 21
1.11 Visual Basic for Application ................................................................................ 31
1.11.1 Creating Macro ...................................................................................... 32
1.11.2 Recording Macro .................................................................................... 35
1.11.3 Procedure .............................................................................................. 37
1.11.4 Running Macro ....................................................................................... 39
1.11.5 VBA Dictionary ....................................................................................... 40
CHAPTER 2: EXCEL FUNCTIONS .............................................................................. 41
2.1 Math and Trigonometry Functions ...................................................................... 41
2.2 Logical Functions ................................................................................................ 44
2.3 Lookup Functions ............................................................................................... 47
2.4 Text Functions .................................................................................................... 49
2.5 Data Analysis Functions ..................................................................................... 52
2.5.1 Linear Regression .................................................................................. 52
2.5.2 Polynomial Regression .......................................................................... 61
2.5.3 Interpolation ........................................................................................... 63
2.5.4 Statistical Data ....................................................................................... 68
2.5.5 Circular Reference ................................................................................. 75

CHAPTER 3: CREATING MACRO .............................................................................. 79
3.1 Function Procedure ............................................................................................ 79
3.2 Sub Procedure .................................................................................................... 84
3.3 Control Structures ............................................................................................... 86
3.3.1 Looping .................................................................................................. 86
3.3.2 Branching ............................................................................................... 89
3.4 User Defined Function Problems ........................................................................ 93
3.5 Structure of Program ........................................................................................ 108
3.5.1 Input Output Form ................................................................................ 109
3.5.2 Work With Modules .............................................................................. 110
3.5.3 Tips ...................................................................................................... 112
3.6 Chart Macro ...................................................................................................... 114
3.7 Manipulation on Program Steps ....................................................................... 120
CHAPTER 4: MATRIX PROGRAM ............................................................................ 124
4.1 Matrix Definition ................................................................................................ 124
4.1.1 Types of Matrix .................................................................................... 125
4.1.2 Matrix Operation .................................................................................. 127
4.2 Program for Matrix Operations ......................................................................... 137
4.3 Matrix Method for Structural Analysis ............................................................... 145
4.3.1 Upper Structure .................................................................................... 145
4.3.2 Sub Structure ....................................................................................... 147
CHAPTER 5: NUMERICAL METHOD ........................................................................ 148
5.1 Numerical Integration........................................................................................ 148
5.2 Numerical Differentiation .................................................................................. 152
CHAPTER 6: PROGRAM FOR 2D FRAME STRUCTURE ANALYSIS ..................... 157
6.1 Case Example .................................................................................................. 157
6.2 Sign Convention for Diagram ............................................................................ 173
6.3 Application ........................................................................................................ 175
CHAPTER 7: PROGRAM FOR 2D TRUSS STRUCTURE ANALYSIS ...................... 176
7.1 Case Example .................................................................................................. 177

7.2 Application ........................................................................................................ 189
CHAPTER 8: BEAM ON ELASTIC FOUNDATION .................................................... 193
8.1 Case Example .................................................................................................. 195
8.1 Application ........................................................................................................ 199
CHAPTER 9: LATERALLY LOADED STRUCTURE ................................................. 202
9.1 Case Example .................................................................................................. 202
9.2 Application ........................................................................................................ 209
CHAPTER 10: ONE DIMENSIONAL CONSOLIDATION ........................................... 212
10.1 Application 1 ..................................................................................................... 218
10.2 Application 2 ..................................................................................................... 219
CHAPTER 11: AUTOCAD SCRIPT FILE ................................................................... 222
11.1 Creating Scripts in Worksheet .......................................................................... 222
11.2 Creating Scripts in VBA .................................................................................... 239

Chapter 1:
BASICS OF EXCEL
The work performed by Excel is basically the job of entering data which is then processed to
obtain the desired results. It is, in principle, the same as entering data into an electronic
calculator. However, the data entered here consists of various types and coupled with existing
Excel facilities makes it possible to present the appearance of numbers, text, associated tables,
graphs and a database. This makes Excel becomes well integrated to create a reporting text
and the problem analysis as well.
Each version of Excel to be developed to always make changes and additions of new facilities,
while is still maintaining compatibility with previous versions. But the changes do not alter the
basic features of this software as a worksheet for computing applications. In this book, we will
be working with Microsoft Excel 2007.
1.1 Worksheet and Workbook

When opening Excel, by default Book1 is the name of the first workbook. This workbook
consists of 3 worksheets named Sheet1, Sheet2 and Sheet3. Excel worksheet is also referred
to as a spreadsheet that is the sheet for processing text and numbers.
Figure 1.1 shows the elements of a workbook. At the top of the page there is a title bar displays
the workbook name. Underneath, there is a Ribbon, new interface intoduced in Excel 2007,
which is a navigation tool replaces the menu and tool bars in earlier Excel version as a tool of
access to Excel commands. All commands are grouped and placed into tabs for particular
purposes, thus a tab contains groups of commands.
Worksheet is divided into rows and columns. In Excel 2007, the number of columns and rows
have been improved from previous versions. Columns are from A to Z, then AA, AB to XFD
(16,384 columns), while the rows start from 1 to 1,048,576. The intersection point of column and
row forms a cell as a place to fill data. Each cell has an address referred to by a column and a
row, for example cell B4 is a cell in column B and 4th row. The address of the selected cell can
be seen in the Name Box below the Ribbon.

A collection of cells is called a cell range, forming an array extending horizontally or vertically.
For example range A1: A5 is a collection of cells from A1 to A5, forming an array of 5 x 1 or the
range A1: C5 form a 5 x 3 array, and so on.
The name of the worksheet can be replaced by other names, by clicking twice on the sheet tab,
press Delete to clear the name, then write the new name; or by right-clicking on the tab and
then on the shortcut menu click Rename to change its name.
1.2 Data Type

The input data in the spreadsheet can be divided into some types that are:
a. Text: alphabet characters and text: A, B, -Z, AB, A2, Computers ...
b. Numbers: numerical data: 1,2,3,0, -1, -2,4,5.85...
c. Date: the date data typically refers to the setting of a computer calendar or formulated
in the calculation.
d. Hours: Data of hours generally are referred to computer time setting or formulated in the
calculation.
e. Formula: mathematical expressions that calculate two or more values produce a new
value.
f. Function: a function that is used for various applications such as calculations, finance,
mathematics and trigonometry, statistics, database, logic and others.

Quick Access Toolbar Tittle Bar e
Close
C
Command Tab Close Window
Command
ds related to the tab
Ribbo
on
Formulla Bar
Name Box
Active Cell
Scroll Bar
Sheet Tab
Figure 1.1: Ex
xcel 2007 Worksheet
W
Table 1.1: Excel Basic Commands
Shortcut
S
Name Descriiption
o Icon
or
Offiice Button Iccon to basicc commandss below:
S
Save a file with
w the nam me of the wo orkbook thatt
Save
iss being done e
S
Store the filee with another name in nto a folderr
a
and directorry, or to a disc. By default file e
e
extension off .xlsx is give
en (free macros). If the
e
Save As
s A
Alt+FA
file contains macro (VBA A) and needss to activate
e
e
each time it is opened, an extension of .xlsm
e
extension mu ust be given.
New O
Open a new workbook
Open O
Open an exissting workbo ook
Print P
Print a workssheet throug gh the printerr device
An Introdduction to: Excel

E for Civiil Engineers 8
http://excceljunior.blog
gspot.co.id
Close Close the wo
C orkbook (Clo
ose Window icon)
Exitt C
Close the pro
ogram and exit
e Excel (iccon Close)
Quiick Access
A
Too
olbar:
Unddo Cancel the la
C ast job
Red
do R
Repeat the la
ast job
Tab
b HOME
Copy the da
C ata of a ce ell or range of cells to o
a
another loca
ation. To select
s the data
d to bee
Copy
c
copied or cut, then clickk on the cell or range off
c
cells where the
t data is th hen press Cttrl+C
M
Move the daata to anothe er location. Commands
s
such as Cut, Copy and Paste can be found in
Cut the pop-up menu by clicking the right r mouse
b
button on th
he highlighte ed cells or byb pressing
C
Ctrl+X
P the data
Put a into a new w cell or ran nge of cells.
P
Paste comm mand is done e after preceeding Cut orr
Paste
C
Copy. To plaace the data, click the ne ew location,
then press Enter
E or by pressing Ctrl+V
D
Display ert dialog box. Used to insert a new
Inse w
Insert > Cells Alt+IE c
cell, row orr column at a the pointter position
(highlighted cell)
c
Insert > Sheet
Alt+IR Insert a new row at the pointer
p positiion
Row
ws
Insert > Sheet
Alt+IC Insert a new column at th
he pointer position
Columns
A
Alt+HO
Format .> Format row, consists of the row height setting,
F
(H
Home >
Row
w Height o hide and to unhide row
to w
F
Format)
Format column, consistts of the co
F olumn width
Alt+HO
A
Format .> s
setting, automatic width (autofit sele
ection), hide
e
(H
Home >
Column Width a show co
and olumns, and showing infformation off
F
Format)
the standard column widdth
Alt+HO
A
Format cell such as numbers, text, font used,
F
Format > Cell (H
Home >
b
borders, ors and workksheet protection
colo
F
Format)

gspot.co.id
Table 1.2 Operators and mathematical relationships
Operator Description
+ Summation
Subtraction
* Multiplication
/ Division
% Percent
^ Exponentiation
Relationships
= Equal to
> Greater than
< Less than
<> Not equal to
>= Greater then or equal to
<= Less than or equal to
1.3 Formula
By definition, a formula is a mathematical expression to calculate the results of two or more
values. Formula can consist of numbers, mathematical operators, functions, reference cell or
range of cells. Cells and cell range are often given a name, for example "A" to B2, or "B" for a
range of D2: D6. Naming cells will be discussed later in this section. All formula always starts
with an equal sign (=). For example, in cell B4 is written the following formula:

=B2+B3 summing the data in cell B2 to cell B3
=B4*B5 multiplying the data in cell B4 to cell B5
=SUM(D2:D4) summing the number of cells D2, D3 and D4
When a formula to be copied, it needs to be considered what is the type of cell in the formula.
Relative cell is the cell that will adjust to its new place when copied. Cell C1 that contain formula
=A1+B1 will be =C2+D2 in Cell E2. The result is as shown in Figure 1.2.
Absolute cell is a cell that does not change the address if copied to another place. The notation
is to add $ before the name of a column or row number, for example: $A$4. Adding $ in the
name of any column or row number will only change one address. Column A in $A4 will remain
when copied, but the row number will adjust to its new location. And vice versa with cell A$4.
Such cells are called semi-absolute cell.
A B C D E F
1 20 30 =A1+B1
2 =C2+D2
3 =$A$1+B1
4 =$A$1+D2
5 =$A1+B1
6 =$A2+D2
7 =A$1+B7
8 =C$1+D8
9
10
Figure 1.2: Copied formula and the results in spreadsheet
Formula stating cell relationships or range of cells, for example, =B4+C30 or =SUM (A4: C20)
are more difficult to read than the mathematical relationship with the more practical variable
name, for example, =x+n or SUM(B). Name of cell or cell range is created by clicking Formula >
Define Name. This name applies to all sheets in a workbook. If it only represents the specific
sheet, for example, Sheet1, the writing is Sheet1!x, i.e. to the variable x in Sheet1. The use of
variables is highly recommended because it will simplify a formula. The names used and the
locations can be seen in Name Manager and you can also create, edit and delete a name.

1.4 Built-In Function
Excel has many built-in functions to build complex formulas, some of them are shown in Table
1.3. Functions such as mathematic and trigonometry or statistic, for example, are the most
common functions used in engineering practices e.g. to produce calculation data sheet in
laboratory soil testing. Before using a function, it is advised to know well the function and its
arguments. The reference could be found the Help menu or press F1.
Table 1.3 Excel Built-in Functions
A. Math and Trigonometry

ABS(num) Returns the absolute value of a number
ACOS(num) Returns the arccosine of a number
ACOSH(num) Returns the inverse hyperbolic cosine of a
number
ATAN(num) Returns the arctangent of a number
EXP(num) Returns e raised to the power of a given number
FACT(num) Returns the factorial of a number
INT(num) Rounds a number down to the nearest integer
LOG(num, base) Returns the logarithm of a number to a specified
base
LOG10(num) Returns the base-10 logarithm of a number
MDETERM(array) Returns the matrix determinant of an array
MINVERS(array) Returns the matrix inverse of an array
MMULT(array1,array2) multiplying of 2 arrays
PI() Value of pi = 3.141592654
RAND() Random value between antara 0 dan 1
SIGN(num) Sign of number. Sign 1 or 0 or -1 if the number

is positive, zero or negative, respectively.
SIN(num) Returns the sine of the given angle
SINH(num) Returns the hyperbolic sine of a number
SQRT(num) Returns square root of a number
SUM (num1,num2,) Add the numbers
SUMPRODUCT(array1,array2) Multiplies corresponding components in the
given arrays
TAN(num) Returns the tangent of a number
TANH(num) Returns the hyperbolic tangent of a number
B. IS Function
ISBLANK(value) TRUE if value is empty
ISLOGICAL(value) TRUE if value is logical value (TRUE or FALSE)
ISNUMBER(value) TRUE if value is a number
ISTEXT(value) TRUE if value is a text
C. Statistic
AVERAGE(num1,num2,) Average value of the numbers
COUNT(value1,value2,) Counts the number of cells that contain numbers
within the list of arguments
COUNTA(value1,value2,) Counts the number of cells that are not empty
within the list of arguments
LINEST(ys,xs,const,stats) Returns the parameters of a linear trend. Const
and stats are logical values (see in Excel Help)
MAX(num1,num2,) Maximum value in a list of arguments
MIN(num1,num2,) Minimum value in a list of arguments
SLOPE(ys,xs) Returns the slope of linear regression line
INTERCEPT(ys,xs) Returns the intercept of the linear regression line
TREND(ys,xs, new xs,const) Returns values along a linear trend. Const is a
logical value specifying whether to set the
constant b = 0 in y = mx + b relationship
D. Lookup and Reference
COLUMNS(reference) Returns the column number of the given
reference
INDEX(array,row_num,column_num) To choose a value from a reference or array
ROWS(reference) Returns the row number of the given reference

TRANSPOSE(array) Returns the transpose of an array
HLOOKUP(value,table,row_in) Search for a value based on the row index in the
table data arranged horizontally
VLOOKUP(value,table,col_in) Search for a value based on the column index in
the table data arranged vertically
E. Logical
AND(logical1,logical2,..) Returns TRUE if all of its arguments are TRUE
IF(log_value,value_if_ true,value_if_false) Specifies a logical test to perform
NOT(logical) Reverses the logic of its argument: NOT(TRUE)
= FALSE
OR(logical1,logical2) TRUE if one of its argument is TRUE
F. Text
CHAR(num) Returns the character specified by the code
number
EXACT(text1,text2) Checks between two text strings and returns
TRUE if they are exactly the same, otherwise
returns FALSE
FIND(text1,text2,start_num) Finds one text value within another with start
number
LEFT(text, num_character) Returns the first character based on the
specified number of characters
LEN(text) Returns the number of characters in a text string
RIGHT(text, num_character) Returns the last character or characters in a text
string, based on the specified number of
characters
TRIM(text) Removes spaces from text except for single
spaces between words
UPPER(text) Converts text to uppercase
1.5 Array Formula

By Excel definition, array formula is a formula that can perform multiple calculations and then
return either a single result or multiple results. Array formulas act on two or more sets of values
known as array arguments. Each array argument must have the same number of rows and
columns. For an example, the built-in LINEST returns two results which are Slope and Y-

Intercept. This function gives linear regression lines that fit with the known x and y-values, as
shown in Figure 1.3.
A B
1 x y
2 -1.0 -5.0
3 -0.5 0.0
4 1.0 5.0
5 2.0 4.0
6 3.0 0.5
7 4.0 -5.0
8 5.0 -12.0
9 Slope Y-intercept 9 Slope Y-intercept
10 -1.204 0.537 10 ={LINEST(B2:B8,A2:A8)} {=LINEST(B2:B8,A2:A8)}
Figure 1.3: Array formula for linear regression and the results
Array formula is entered by pressing Ctrl+Shift+Enter and Excel will automatically insert the
formula in brackets ({}). To display values returned by LINEST function, make LINEST formula
in cell A10 as shown in Figure 1.3, and then select the range A10: B10 > press F2 > Ctrl+Shift+
Enter.
Here are some other examples of array formula:
{=SUM(A2: B2*A3:B3)} is a single result
{=TREND(A2:A6,B2:B6)} is a separately 5 results
Alternatively, the INDEX function can be used to return the values, as shown below. LINEST
has two results indexed sequentially that are Slope and Y-intercept.
9 Slope Y-intercept 9 Slope Y-intercept

10 -1.204 0.537 10 =INDEX(LINEST(B2:B8,A2:A8),1) =INDEX(LINEST(B2:B8,A2:A8),2)
The example of array formula for linear regression using LINEST and TREND function with
associated chart can be seen in Regression.xls file on the CD included in this book.

1.6 Data Formatting
Data and cells formatting can be done through the Home tab > Format > Cells. In the Format
Cells dialog box, there are several sheet tabs which are Number, Alignment, Font, Border,
Fill and Protection. Number formatting in Number tab is divided into several categories such
as general, number, date and so on. For example, for custom category or user preference is
shown in Table 1.4. Alignment is the setting for layout of text in a cell such as vertical and
horizontal alignment, control and orientation of text in a cell. Fonts are the settings that relate to
the character, such as the font face used (Arial, Times New Roman, Tahoma, ...), style (italic,
bold ...), size, color and so on. Border is to make the margins, line types (straight, dotted, thin,
thick), and color. Fill to create a displayed cell with shading (horizontal, diagonal, dot, ...),
shading color and the color of the cell background. Protection is used to select the condition
(locked or not) of data and formulas (hidden or not) in a cell when worksheet protection is
enabled. To protect a worksheet click Review tab > in Changes group click Protect Sheet.
Table 1.4: Custom format for numbers
Format Displayed
0.00 2343.00
0.00E+00 2.34E+03
##0.0+0 2.3E+3
#,##0.00 2,343.00
#,##0.00% 234,300.00%
0.0 m 2343.0 m
Column width and row height can be adjusted by: first way, by dragging the mouse (hold the left
mouse button while moving) when the pointer position at the boundaries of the column and row
headings (there is a sign "J" or "L"), or second way is by right-clicking the mouse on the column
or row heading and from the shortcut menu click Row Height or Colomn Width, then enter
desired value in the textbox. The other way is through Home > Format.

1.7 Error Message
Error message will be appeared if a formula is not working as it should be, for examples,
arguments in a formula are not complete, the data type does not match, or wrong in developing
the logic. Another case, for example is division by 0 (zero) or the width of the column is not wide
enough to display a result. The message started with the # sign followed by the type of error.
Table 1.5 shows some error messages that may appear.
Table 1.5 Error Messages
Messages Remarks
Cell contains number, date or time that is wider
##### than the cell width or minus result in the date and
time format
(Not Available) when a value is not available to a
#N/A
function or formula.
#DIV/0! Error due to division by zero (0).
#NAME? Excel does not recognize text in a formula
#REF! When a reference cell is not valid
The types of arguments used in the function are
#VALUE! wrong. The argument can be a numerical value,
text, name, label, reference cell, cell range and
function..
1.8 Printing
Before printing the worksheet to paper, you can set the area by blocking the area that will be
printed (press Shift + "J" or "L" or drag the mouse). Afterward, click the Page Layout tab >
Print Area > Set Print Area. Printing area can be seen in Office Button > Print > Print
Preview or to press the shortcut key Ctrl+F2. Printing then performed through > Print > OK
or pressing Ctrl+P.

The limits of printing page can also be specified through the View tab > Page Break Preview
instead current Normal view. This will show all the data in worksheet to be printed with a page
number. You can also adjust the boundaries of the printing page by dragging the mouse on the
blue boundary line.
Page configuration is then set through the Page Layout tab. Click the Print Titles in the Page
Setup group, and a Page Setup dialog box will appear with several sheets tabs which are
Page, Margins, Header/ Footer, and Sheet. Page is to select the orientation and size of paper,
and the scale of worksheet on paper; Margins is to set the boundaries of the text in the paper;
Header/Footer is to make custom header and footer; and Sheet is to set print area, showing
grids, row and column headings, and the print order.
1.9 Making Charts

A set of data can be well delivered and communicated through a chart image shown the
correlation between the data. In Excel there are many types of chart depictions such as Bar,
Column, Line, Pie, Area and so on, plus a 3-dimensional view. Chart is created through the
Insert tab > Chart, and there will be Charts group as shown in Figure 1.4. For this example, we
will use the data from Figure 1.3. Select Scatter > Scatter with Smooth Lines and Markers
because smooth line between points to be made. Click Select Data on the Ribbon interface to
display the Select Data Source dialog box as shown below.

Enter x and y-values data in the Chart Data Range by clicking a button at the right side of the
input box and then point the pointer to cell A2 and drag the mouse from cell A2 to cell B8. Click
OK and a chart will be created as shown in Figure 1.5. In the chart, the regression line is also
shown with the associated equation. Regression line is created by clicking mouse on the data in
the chart (at point or a line), then right click > Add Tendline > Linear > checked the Display
Equation on Chart. Furthermore, if you want to work with the chart thoroughly, click the graph
area to display Chart Tools tab where commands are grouped in accordance with their task
name, which are Type, Data, Chart Layouts, Chart Styles and Move Chart.
Figure 1.4: Charts group displays chart type options

Figure 1.5: Example of Scatter with smooth lines chart type
Note:
As designed for office applications, Excel offers various types of charts. However, we will not
discuss many different types of charts, and but only focused on Scatter chart type with
straight lines that frequently used in this book.

1.10 Engineering Drawing
Chart type of XY Scatter is suitable for use in Civil engineering practices in giving a drawing
presentation that forms lines or elements of structure. The reason is that every line can stand
alone, so that easily modified and formulated for the intended drawing.
A straight line of XY Scatter chart type is determined by the given coordinates of both ends
required for its input data. If both ends of the line called joint, then a line segment is formed by
the coordinates of two joints. Thus, there will be a series of joint data put into a worksheet table
to create straight lines.
To make it easier to understand the intent and purpose of this discussion, below is given
examples and the followed steps.
Example 1
Create 3 continuous lines drawing through 4 points (1 to 4), as shown in figure below.
The joint coordinates are:

x y
Joint 1 0 0
Joint 2 1 1
Joint 3 2 3
Joint 4 3 4
To create continuous lines as the figure above, it takes the following steps:
1. First step: click on the Insert tab > Scatter > Scatter with Straight Lines and Markers.
2. Chart Area displays nothing because no data on it as shown below. Click Select Data
on the Ribbon interface to display the Select Data Source dialog box.

3. A Select Data Source dialog box is displayed as shown in the figure below (left side).
Click Add to display the Edit Series dialog box (right side). Each series that added
(clicking Add) represents a line that requires a pair of x and y values as inputs. Thus, a
line series shall consist of the coordinates of the both ends, the one to be (x1, y1) and
the other is (x2, y2).
4. Following Step 3 is to enter the line coordinate values i.e. series x and series y-values.
When needed, give a name for the series in the Series name input box.
5. Click OK to end.
The intended joint coordinate and lines for chart depiction can be summarized as below:

Coordinate
Line Joint x1 y1 x2 y2
1 1 2 0.0 0.0 1.0 1.0
2 2 3 1.0 1.0 2.0 3.0
3 3 4 2.0 3.0 3.0 6.0
Thus it can be stated that the first line (no.1) is a line created by given joint 1 and 2 coordinates
(refers to the figure).
In step 4, what we do is to enter x1 and x2 range of values (colored yellow) into the Series X
values input box, while the Series Y values is filled by column y1 and y2 range of values
(colored green). To enter x-coordinates, press Ctrl+left mouse button click at x1 column and
repeat left mouse button click at x2 column to get the range of cells. Do the same way to
enter y-coordinates.
Repeat step 3 to 5 to make 2nd and 3rd line. The result will be shown as below:

Example 2
Draw a simple picture of one floor building as below:
As in Example 1, to simplyfy the portrayal of the drawing it needs to produce joints and lines
coordinate tables. The buiding is composed of 6 joints and 6 straight lines that can be built up
with the following numbering system:

Based on
n the picture
e above, the following tables can be formed:
The jointt coordinates

s for assigning lines on the right tab
ble is obviou
usly the repe
etition of inpu
utting
task of th
he previous joint coordinates on the
e left table. To avoid re
epeated data
a entry manually,
UP function can be used
VLOOKU d by adopting the x, y joint coordinates data of the
t left table
e as a
reference btain x1, y1 and x2, y2
e table to ob 2 coordinatess of the righ e formulas in the
ht table. The
workshee
et are shown
n as follows::
he drawing, do
To get th d the seque
ence steps as
a in Example 1.
e3
Example
Draw the
e following truss structure that conssists of 21 lines and 12 joints. The joint coordin
nates
are tabulated as belo
ow:

gspot.co.id
Joint x y
1 1.0 0.0
2 1.5 2.5
3 2.0 5.0
4 2.5 10.0
5 3.0 15.0
6 3.5 20.0
7 4.0 15.0
8 4.5 10.0
9 5.0 5.0
10 5.5 2.5
11 6.0 0.0
12 3.5 2.5

Before creating the drawing of truss, it is convenient to form lines coordinates as required for
Series X and Series Y input data. The coordinates of truss lines are tabulated as the following:
1 1 2 1.0 0.0 1.5 2.5
2 2 3 1.5 2.5 2.0 5.0
3 3 4 2.0 5.0 2.5 10.0
4 4 5 2.5 10.0 3.0 15.0
5 5 6 3.0 15.0 3.5 20.0
6 6 7 3.5 20.0 4.0 15.0
7 7 8 4.0 15.0 4.5 10.0
8 8 9 4.5 10.0 5.0 5.0
9 9 10 5.0 5.0 5.5 2.5
10 10 11 5.5 2.5 6.0 0.0
11 2 12 1.5 2.5 3.5 2.5
12 12 10 3.5 2.5 5.5 2.5
13 3 12 2.0 5.0 3.5 2.5
14 12 9 3.5 2.5 5.0 5.0
15 3 8 2.0 5.0 4.5 10.0
16 3 9 2.0 5.0 5.0 5.0
17 4 9 2.5 10.0 5.0 5.0
18 4 8 2.5 10.0 4.5 10.0
19 4 7 2.5 10.0 4.0 15.0
20 5 8 3.0 15.0 4.5 10.0
21 5 7 3.0 15.0 4.0 15.0
As the previous examples, the line coordinates data entry used VLOOKUP function. The
formulas writing with VLOOKUP can be seen in Lines.xls file on the CD included in this book.
To obtain the truss drawing, do the sequence steps as in Example 1.
Next example is to draw truss structure subjected to horizontal load so that each joint will be
shifted as far as n unit. Suppose that the displacement at the joints will be from 0 at supports to
maximum +0.5 unit at joint 6 (at the top), so it will be in a range of 0 < n < 0.5 to the right
direction. The table of joints coordinates is then modified and labeled for two conditions, before
and after loading, where n is used in the calculation. It will show as follow:

It is noticed that only the x-values will be changed by adding the horizontal displacement value
of n, from 0 at Joint 1 to a maximum of 0.5 units at the Joint 6. Joint 1 and 11 are fixed supports
so their positions remain in place before and after loading.
After forming the joints coordinates table, next is to create lines coordinates table of Series X
and Y. There are now 42 Series X and Y to be created for the condition of before and after
loading. Table below shows lines coordinates before and after loading:

The reference range of data table in VLOOKUP function is now range A7 to E18. Thus, the
coordinates of the lines before and after loading are respectively referred to column BC and DE
or column index of 2.3, and 4.5. The formulation can be seen in Lines.xls file on the CD included
in this book. Lines.xls contains Example 1, 2 and 3 materials.
The drawing of truss lines before and after the horizontal loading are shown in chart below,
which are in blue and red, respectively.

1.11 Visual Basic for Application
One of the strong points of Excel is its macro capability of using macro programming language
to shorten and simplify the repetitive works (automation task). It can not be separated, however,
from the presence of Microsoft Visual Basic application which is by default embedded in as an
Excel command menu. The Visual Basic application for the repetitive works or macros in an
application program (e.g. Excel or Word) is called Visual Basic for Applications (VBA). In
addition to Microsoft Office, many of non-Microsoft software also use VBA, such as AutoCAD,
CorelDraw, Visio and Norton.

VBA can
n also serve
e as a prog
gramming la
anguage to solve man
ny problemss in science
e and
technolog
gy field, for instance to solve a complex iterattions and the analysis of
o civil strucctures
that are not easily or
o can not be
e done rely on the builtt-in functions and sprea
adsheet stan
ndard
commands. In this re
espect, the analogy
a of using
u macro is therefore a repetitive calculation in an
effort to get
g new outp
put when new input is en
ntered.
Microsoftt Visual Bas

sic as well as
a Excel VBA
A is an Obje
ect Oriented
d Programm
ming, which reads
r
and analyzes objects
s exposed by
b Excel thro
ough object library. The reference to
o these objects is
stored in
n Microsoft Excel Objecct Library (E
EXCEL12.O
OLB). In a programming
p g structure for a
calculatio
on applicatio
on, Excel exxposes the property
p of an object which is the va
alue of inputt data
that to be processed
d by VBA and then retu
urns the ressult to a worrksheet. So, in writing in
nputs
what we have actuallly done is storing
s objecct value into a variable, or
o the value of the varia
able =
value of cell
c property
y of the workksheet (= Ra
ange object, that will be mentioned
m n
next).
1.11.1 Creating Macro

M
To create
e macros, ac
ctivate the Developers
D tab through > Excel Options > Popular > Show
S
Develop e Ribbon, so
per tab in the s it looks likke this:
To go to
o the window e VBA programming language is created, click Visual Bas
w where the sic to
open Vis
sual Basic Editor (VB
BE) window > click the
e Insert me
enu > Modu
ule to show
w the
following window:

gspot.co.id
ogramming language is written in a module th
VBA pro hat has been
n added be
efore by Insert >
Module. Below is an gram code) in a procedu
n example of VBA progrramming language (prog ure in
e:
a module
VBA read
ds, analyzes
s and manip
pulates Exce
el's objects through
t their properties and method
ds. A
VBA cod
de for an ob
bject is orde
erly compossed of the object
o name
e, a period, followed byy the
property or method. The propertty is an attribute of an object
o and iss always acccompanied by
b an
gn (=), while the method is the action
equal sig n carried outt with the ob
bject.
ode created above, Sele

In the co ect and Clea a called method, whereas Selectiion in
arContents are
the second row is an
n object reffers to Rang
ge object in the first row
w. Value on the
t third row
w is a
y of the Ran
property w the value of 10. The
nge object with ed to remove the
e code abovve is intende
existing value
v in the range A1 on
n Sheet1 and replace it with a new value
v of 10. The Worksh
heets
("Sheet1") code refers to active sheet
s name,, in this case
e is Sheet1.

gspot.co.id
To see th
he available
e objects in Excel,
E you can
c access through
t the Object Bro
owser by clicking
the icon on VBE
E toolbar orr press F2. Select
S the Excel
E es in the pull down men
librarie nu as
shown be
elow:
Objects are classifie

ed into Cla
asses (on Class
C windo
ow) and the
eir Memberrs (on Mem
mbers
window) that could be
b a properrty, a functiion or a constant. Nam
mes, Cells, FormulaArra
F ay for
instance,, are membe
ers of the Ra
ange class. Excels obje
ects are arra
anged in a hierarchy;
h ob
bjects
contain other
o objects
s from top le
evel to down
n level. Code below sho
ows an exam
mple of an object
o
hierarchyy from the Application ob
bject to the Range
R objecct:
Applicat
tion.Workb
books("Book
k1.xls").W
Worksheets(
("Sheet1").Range("A1
1")
In the co
ode, Range is
i the property of the Worksheet
W ob
bject that returns a Rang
ge object. It uses
arg) syntax where arg is name of
Range(a o the range
e (A1) to return
r a Ra
ange objectt that
represen
nts cell A1. Moreover,
M W
Worksheet iss the propertty of the Wo
orkbook obje
ect that returns a
Workshe
eet object, Workbook
W is the propertty of the Ap
pplication ob
bject that retturn a Workkbook
object, and the top le
evel of this hierarchy
h is the
t Application object wh
hich is Excel.

gspot.co.id
1.11.2 Recording Macro
The best way to figure out how VBA communicates with Excels objects is to record a macro.
This tool gives a quick introduction of all Excel's object as well as how to write code in VBA. For
instance, we would like to record a macro to delete data on a range of cells, say, in a range of
A1 to E10. We need to create an automation steps to clean the existing data as the program
produces a new one. The steps for recording macro are as follows:
1. Click the Developer tab > Record Macro to display the Record Macro dialog box with
the default name of Macro1 for macro name as shown below:
2. Change macro name to ClearOutput
3. Move the pointer to cell A1 and drag the mouse to the cell E10 to assign the range data
that will be removed.
4. Press Delete to delete the contents of range A1 to E10 and then place the pointer in cell
A1.
5. End the process by clicking Stop Recording.

The recorded macro is stored in the workbook, and to see the code you have to open the Visual
Basic Editor window through Developer tab > Visual Basic. The code is created as the
following:.
Sub ClearOutput ()
'
'Macro ClearOutput
'
'
Range("A1:E10").Select
Selection.ClearContents
End Sub
In the recorded macro above there are VBA statements such as Range (), Select, Selection and
ClearContents. You may a bit familiar with these statements that have been briefly explained in
the previous section.
The recording macro tool is very useful for writing complex code or you may indeed never
imagined before, like how to manipulate a chart object such as to copy, to create a title, to
change the colors, making marker, automation for creating lines, and so on. This can all be
done without the need to understand more about the Excels object hierarchy that such long and
may be a bit confused.

1.11.3 Procedure
VBA Macro is a programming language written in a procedure, for instance, ClearOuput macro
that is written in a Sub procedure. Procedure in VBA is divided into three types as follows:
Sub procedure, procedure that begins with Sub, its name and ended with End Sub
Function procedure, begins with Function, its name and ended with End Function
Event procedure, works when there is a certain event, for example, open a workbook,
click the button, activate the worksheet, and so on.
The syntaxes of above types of procedure are as follows:
1. Sub Procedure
Syntax:
[Private | Public | Friend] [Static] Sub name [(arguments)]

[statements]
[Exit Sub]
[statements]
End Sub
Remarks
Public Optional. Indicates that the Sub procedure is accessible to all other procedures in
all modules. If not defined, the procedure by default is Public
Private Optional. Indicates that the Sub procedure is accessible only to other procedures
in the module where it is declared
Friend Optional. Used only in a class or object module. Can be accessed from other
procedure in all modules in a workbook
Static Optional. Indicates that the Sub procedure's local variables are preserved
between calls.
name Required. Name of the Sub; follows standard variable naming conventions.

arglist Optional. List of variables representing arguments that are passed to the
Sub procedure when it is called. Multiple variables are separated by
commas.
statements Optional. Any group of statements to be executed within the Sub

procedure.
2. Function Procedure
Syntax:
[Private | Public | [Static] Function name [(arguments)][As type]

[statements]
End Function
3. Event Procedure
Syntax:
Sub expression.eventname (parameters)
Remarks:
expression variabel that refers to object
eventname name of event
parameters argumens used in Event procedure
The most frequently used of Event procedure is clicking a button (command button) on the
worksheet. Here is the example code:
Private Sub CommandButton1_Click()

'call other procedure
Call Pro_Sub1
End Sub
The steps and examples of the above procedures will be shown in next chapters.

1.11.4 Running Macro
Macros can be executed in the two following ways:
1. Through the Developer tab > Macros > select the macro name > Run. Function
procedure, Private Sub procedure and Sub procedure written with argument, can not be
run by this way.
2. Clicking Button (Form Control) or the Command Button (Active X control) on the
worksheet.
The second way is commonly used because it is done by only one step. However, it can be
done in two ways, which are as follows:
o Create a button through the Developer tab > Insert > Form Controls > Button > there
is a sign "+" on worksheet, and place it on worksheet by dragging the mouse, the same
way to determine its size as well. The Assign Macro dialog box will be displayed as
below. Enter the name of a macro that to be executed and then click OK.
o Macro is executed by clicking a command button on the worksheet which is associated

with an event procedure. To create command button click the Developer tab> Insert >
Active X Controls > Command Button > there is sign "+" at the pointer, and place it on

the worksheet by dragging the mouse, the same way to determine its size as well. The
results is as shown below:
Double-click the command button to write code in the following procedure:
Private Sub CommandButton1_Click()
End Sub
When a command button is created, the code program is inactive or when the Design Mode is
enabled. Before run macro, click the Design Mode to disable (until it looks unselected).
1.11.5 VBA Dictionary
Statements or words that are used in Excel-VBA are too many because the capability VBA to
support Excel for diverse goals. Some statements are similar to those used by Excel, for
examples, operator and mathematical relationships, some strings functions and math functions.
In this book, all statements related to this book will not be placed in a special section, but it will
be discussed respectively, on each topic in its section.

Chapter 2:
EXCEL FUNCTIONS
Almost identical to a function definition in mathematics, a function in Excel works based on the
given arguments and written as follows:
Result Name Argument/s
y = f (a,b,c,)
Arguments are written after the function name, between parentheses, which are the values used
to perform the operation. The argument can be a numerical value, text, references or range of
cell, name, label or a function.
Excel functions are listed by category such as, math and trigonometry, statistics, finance, logic,
and so on. In addition, Excel also provides macrosheet and Visual Basic for Applications (VBA)
to create a function that is defined by the user (user-defined function).
2.1 Math and Trigonometry Functions

ABS
Returns the absolute value of a number.
Syntax:
ABS(number)
Example:
ABS(4) = 4
ABS(-4) = 4
SQRT(ABS(-81)) = 9

INT
Rounds a number down to the nearest integer
Syntax:
INT(number)
Example:
INT(8.9) = 8
INT(-8.9) = -9
TRUNC
Makes a number to an integer by removing the fractional part of the number. To specify the
precision of the truncation, enter a specified number of digits after the number.
Syntax:
TRUNC (number,number_digits)
Example:
TRUNC(8.9) = 8
TRUNC(-8.9) = -8
TRUNC(PI ()) = 3
TRUNC(PI (), 3) = 3,141
ROUND
Rounds a number to a specified number of digits.
Syntax:
ROUND(number,num_digits)
Example:

ROUND(3.25,1) = 3.3
ROUND(3.247,1) = 3.2
ROUND(3.247,0) = 3
ROUND(32.47,-1) = 30
ROUNDOWN
Rounds a number down to a specified number of digits.
Syntax:
ROUNDDOWN (number,number_digits)
Example:
ROUNDDOWN(3.3,0) = 3
ROUNDDOWN(66.8,0) = 66
ROUNDDOWN(3.14159,3) = 3,141
ROUNDDOWN(-5.24621,2) = -5.24
ROUNDDOWN(1550.24621,-2) = 1500
ROUNDUP
Rounds a number up to the desired number of digits.
Syntax:
ROUNDUP (number, num_digits)
Example:
ROUNDUP (3.3, 0) = 4
ROUNDUP (66.8, 0) = 67
ROUNDUP (3.14159, 3) = 3,142

ROUNDUP (-5.24621, 2) = -5.25
ROUNDUP (1550.24621,-2) = 1600
SUMPRODUCT
The number resulted by multiplying corresponding components of arrays
Syntax:
SUMPRODUCT(array1,array2,)
Example
SUMPRODUCT(1,2,3) = 1 x 2 x 3 = 6
SUMPRODUCT({1,2,3},{4,5,6}) = 1 x 4 + 2 x 5 + 3 x 6 = 36
SUMPRODUCT({1,2,3,hello},{4,5,6,5}) = 36
2.2 Logical Functions

IF
This function has two values of results, which are TRUE and FALSE, with logical test, if met will
do the TRUE value. If the result returns a text, then the text must be between quotes ("text").
Syntax:
IF(logical_test,value_ if true,value_if_ false)
Example 1:
In a class, students which have exam scores below or equal to (<=) 55 will FAIL the exam, and
above the value will PASS the exam.

e 2: Nested IFs with two logical testss
Example
Studentss who FAIL with exam score > 45

5 have cha
ance to take
e a supplem
mentary exa
am to
increase their scores
s. Use the fu
unction to find out the na
ame of stude
ent:
AND
Returns TRUE if alll of its argu

uments are TRUE, and
d returns FA
ALSE if one
e its argume
ent is
FALSE. Generally
G ne
ested with IF
F function.
Syntax:
AND(logiical_1, logical_2, ...)
Example
e 1:
AND(2*2
2=4,2+2=4) Returnss TRUE
100,4<100,102<100)
AND(2<1 Returnss FALSE

gspot.co.id
Example
e 2:
Now, witth the same

e class and
d the same
e students in Example 1, but with
h different exam
e
subjects:: a course score
s can he
elp the stude
ent that receives exam score > 40 and <= 55. The
given forrmula is:
(2 x Exam
m score + 1 x Course sccore) / 3.
COUNTIF
g the number of cells in a range of ce

Counting ells accordin
ng to given criteria.
c
Syntax:
COUNTIF(range,crite
eria)
Example
e1
Looking for
f the numb
ber of studen
nts with speccified value for Example
e 2 above:
COUNTIF(D2: D6,60
0) = 2
COUNTIF(D2: D6,> 60) = 1
OR
ult is TRUE iff one of its arguments

The resu a iss TRUE. Gen
nerally neste
ed the IF fun
nction.
Syntax:

gspot.co.id
OR(logical1,logical2, ...)
Example:
OR(2*2 = 4,2+2=4) Returns TRUE
OR(2<100,4<100,102<100) Returns TRUE
2.3 Lookup Functions

VLOOKUP
Search for a value based on the column index in the table data arranged vertically
Syntax:
VLOOKUP(value,table,col_index,range_lookup)
Example
The experimental result of laboratory test is affected by the size and weight of the used tools.
This example shows how to use VLOOKUP function to read the size and weight of the ring
based on the Ring Calibration table. Suppose it is intended to search the diameter and weight
of ring no. 3 and 7. The table data refers to range A4 to F11, ignoring the column heading.
A B C D E F
RING
1 CALIBRATION
2
Ring Height Diameter Weight Area Volume
3 No. (cm) (cm) (gram) (cm2) (cm3)
4 1 1.950 5.000 32.850 19.635 38.288
5 2 2.000 5.000 32.000 19.635 39.270
6 3 2.000 5.000 33.440 19.635 39.270
7 4 2.000 5.000 34.980 19.635 39.270
8 5 1.950 5.000 33.700 19.635 38.288

9 6 1.925 4.985 32.000 19.517 37.571
10 7 2.000 5.000 33.000 19.635 39.270
11 8 1.975 5.125 34.000 20.629 39.711
VLOOKUP result:
A B C D E F
Ring Diameter Weight
14 No. (cm) (gram)
15 3 5.000 33.440
16 7 5.000 33.000
Formula in cell B15 to C16:
A B C
Ring
Diameter (cm) Weight (gram)
14 No.
15 3 =VLOOKUP(A15,$A$4:$F$11,3) =VLOOKUP(A15,$A$4:$F$11,4)
16 7 =VLOOKUP(A16,$A$4:$F$11,3) =VLOOKUP(A16,$A$4:$F$11,4)
Range_lookup is a logical value that specifies whether VLOOKUP finds an exact match or an
approximate match. If TRUE (or omitted), the values in the first column of table_array must be in
ascending sort order (-2, -1,0,1,2, ..) because VLOOKUP may return incorrect. If FALSE, the
table_array values do not need to be sorted. For TRUE condition, if the lookup_value does not
match with table_array, VLOOKUP returns the next largest value that is less than lookup_value.
Otherwise, if FALSE, VLOOKUP will find an exact match, if it is not found, VLOOKUP returns
#N/A (error).
HLOOKUP
Search for a value based on the row index in the table data arranged horizontally. The
application of HLOOKUP are almost identical to VLOOKUP, in exception that row index is
needed to lookup the value.

Syntax:
HLOOKUP(value,table,row_index,range_lookup)
Example
Transpose (Copy > Paste Special > Transpose) the Ring Calibration table to arrange table
data as shown below. After that, use HLOOKUP the same way as VLOOKUP example above to
search the diameter and weight of ring no. 3 and 7.
Ring No. 1 2 3 4 5 6 7 8
Height (cm) 1.95 2.00 2.00 2.00 1.95 1.93 2.00 1.98
Diameter (cm) 5.00 5.00 5.00 5.00 5.00 4.99 5.00 5.13
Weight (gram) 32.85 32.00 33.44 34.98 33.70 32.00 33.00 34.00
Area (cm2) 19.64 19.64 19.64 19.64 19.64 19.52 19.64 20.63
Volume (cm3) 38.29 39.27 39.27 39.27 38.29 37.57 39.27 39.71
2.4 Text Functions

LEFT, RIGHT, MID
Returns characters in a text string based on the specified number of characters
=LEFT(text,num_chars)
Returns the first character
=RIGHT(text,num_chars)
Returns the last character
=MID(text,start_num,num_chars)
Returns characters starting at specified position
Example

A B
1 ENGINEERING ENG =LEFT(A1,3)
2 745.56 RING =RIGHT(A1,4)
3 ENGINEER =MID(A1,1,8)
4 745.5 =LEFT(A2,5)
CONCATENATE
Combines two or more strings from different cells into one string
Example
A B C D
1
Civil
2 Civil Engineer Engineer
3
4 Point Coordinate
5 x y z
6 2.0 4.0 0.0 2,4,0
7 2.0 4.0 -1.5 2,4,-1.5
8 2.0 4.0 -1.5 (2,4,-1.5)
9
D
1
Civil
2 Engineer =CONCATENATE(A2," ",B2)
3
4 Coordinate
5
6 2,4,0 =CONCATENATE(A6,",",B6,",",C6)
7 2,4,-1.5 =CONCATENATE(A7,",",B7,",",C7)

8 (2,4,-1.5) =CONCATENATE("(",A8,",",B8,",",C8,")")
9
Another way is to use the way of writing as the following:
numeric & "text" & numeric
Example
A B C D
1
Civil
2 Civil Engineer Engineer
3
4 Point Coordinate
5 x y z
6 2.0 4.0 0.0 2,4,0
7 2.0 4.0 -1.5 2,4,-1.5
8 2.0 4.0 -1.5 (2,4,-1.5)
9
D
1
Civil
2 Engineer =A2&" "&B2
3
4 Coordinate
5
6 2,4,0 =A6&","&B6&","&C6
7 2,4,-1.5 =A7&" "&B7&","&C7
8 (2,4,-1.5) ="("&A8&","&B8&","&C8&")"
9

2.5 Data Analysis Functions
Excel has many built-in functions that are used to analyze data obtained from experimental
result. The objective of the analysis is to obtain theoretical parameters that give the best
relationship between theory and experimental result. Three methods will be discussed here
which are: linear regression, polynomial regression and interpolation.
2.5.1 Linear Regression
Linear regression is to determine a straight line that fits or the most closely fits to a number of
points data, providing a linear relationship between two variables. The method used to obtain
the line is called the least squares method. Thus, linear regression consists of a series of points
that fit to a number (n) of points data (xi, yi) written into a straight line equation:
y = Ax + B (2.1)
where:
A = slope of the line,
B = the intersection of the straight line to the Y-axis
The accuracy of the straight line over a number of points data is evaluated by a total deviation
E, which is the sum of squares of the distance e between the points data and the fitted points:
n n
E= e = ( Ax + B y )
i =1
2
i
i =1
i i
2
(2.2)

By making E to minimum, A and B values therefore can be determined to obtain the equation
2.1. Beside E, R2 value is also used for the accuracy measurement by a relationship below:
[y
i =1
i f ( xi )]
2
R2 = 1 2
n
1 n (2.3)
i =1
y y i
2
i
n i =1

R2 value varies from 0 to 1. R2 = 1 is when the regression line coincides with the data. In
polynomial regression, the higher-order of polynomial, the closer R2 value to = 1.
Regression in Excel can be obtained using the TREND function, SLOPE, INTERCEPT and
LINEST or with Trendline, a regression line from data relations in an XY coordinate system.
Trendline is created by the following steps: right-clicking the mouse when the pointer is on one
of the points data in the graph to display a box menu > click Add Trendline to display the
Format Trendline dialog box > select Linear for linear regression.
Example 1: Linear Regression
The soil shear strength determination from a laboratory test obtains the result as shown in
Figure 2.1. In general, the depiction of data generated from the test is rarely a straight line to
show failure envelope, therefore, it needs a fitted straight line to represent the data. Regression
is used to obtain shear strength parameters, which are cohesion (c) and shear angle in () of
the soil mass.
In direct shear test, the failure envelope is a straight line that is expressed into equation:
= c + n tan
Where, shear strength

c cohesion intercept
n normal stress
shear angle or slope of failure envelope

ure 2.1: Regression on Excel Charrt
Figu
Shear sttrength in th
his model iss a linear fu
unction overr the normal stress, wh
here c and tan
t
respectivvely expresses Y-axis intercept and the slope off the line.
In Chart, a regressio
on can be done with Trrendline by the
t following
g steps: righ
ht clicking mouse
m
e pointer is on one of the points da
when the ata on Charrt to displayy a box men
nu. From the
e box
menu se
elect Add Trrendline to display Form
mat Trendliine dialog box
b > Linearr > Backwa
ard to
intercept Y-axis > clic
ck Display Equation
E n Chart. The
on e result is ass shown in Figure
F 2.1.

gspot.co.id
The values of A and B that are shown in line equation, y = Ax + B in direct shear test are equal
to those calculated with the following functions:
A = slope of the line = tan = SLOPE(Y,X)
B = Y-intercept = INTERCEPT(Y,X)
Shear angle therefore can be calculated with ATAN function (in degree):
=ATAN(SLOPE(Y,X))*180/PI()
LINEST and TREND function used and written as follows:
=INDEX(LINEST(Y,X),1) to obtain Slope,
=INDEX(LINEST(Y,X),2) to obtain Y-intercept

In practicces, TREND
D function is entered in a series of worksheet
w formulas to ob
btain fitted points
p
represen
nts a series of
o data (arrayy Y, array X).
X The formu
ula can be written
w as follows:
Shear strrength can also

a ajor 1 and minor 3 principal stressses in the trriaxial
be exprressed in ma
test by plotting (1 3) againsst (1 + 3). The stress condition that fit to a number of points
p
eferred to as
data is re s modified fa
ailure envelope, which exxpressed intto equation:
(1 3) = a + (1 + 3) tan
where a and are the

t modified
d shear stren eters. The parameters c and are then
ngth parame
given thrrough relations below:
= sin1(tan
( )
a
c=
cos
e 2: Logarithm
Example mic Regresssion
It is often
n found that the data re
elationship frrom a test iss presented where the X-axis
X is ma
ade in
log scale
e as shown figure 2.2. Th
he equation 2.1 now beccomes:
y = A ln (x) + B (2.4)
where A and B are constants
The exam
mple data here
h is adop aboratory liquid limit (LL
pted from la L) test, i.e. to
t determine
e soil
moisture content at 25
2 blows. Be
elow are the related informations:

gspot.co.id
Moisture content (Wn) = [(weight of wet soil - weight of dry soil) / (weight of dry soil)] x
100%.
Water content' (Wn') = fitted Wn in the semi-log relationship.
N = blows number
Liquid Limit Sheet
85.00
80.00
Water Content (%)
y = -10.887ln(x) + 105.972
75.00
70.93
70.00
65.00
Data
60.00 Log. Trendline

25
10 100
No. of Blows (N)
Figure 2.2: Regression on Semi-log Liquid Limit Test
Figure 2.2 shows the LL test result and a regression line created through Add Trendline
command. The fitted Wn value (Wn') that represents a number laboratory Wn for each N value
can be obtained through the equation below:
y = -10.887ln(x) + 105.972
Using the above equation, at N = 25 blows, the related fitted Wn = 70.93%.
The calculation sheet of LL test in the worksheet can be made as follows:

A B C D E F G
1 LIQUID LIMIT TEST
2
3
Weight of Can + Can + Water Fitted
4 No. Can (gr) Soil (gr) Dry Soil (gr) Blows Content (%) Wn (%)
Set
5 Wc Ws Wd N Wn Wn
6 1 8.3 52.46 33.29 14 76.71 77.24
7 2 2.8 55.61 32.96 19 75.10 73.92
8 3 2.8 57.14 35.12 29 68.13 69.31
9 4 7.8 52.72 34.62 36 67.49 66.96
Liquid Limit (%) = 70.93
3 Water Fitted
4 Content (%) Wn (%)
5 Wn Wn
6 =100*(C6-D6)/(D6-B6) =TREND(Wn,LOG10(N),LOG10(E6))
7 =100*(C7-D7)/(D7-B7) =TREND(Wn,LOG10(N),LOG10(E7))
Liquid Limit (%) =TREND(Wn,LOG10(N),LOG10(25))
The relationship of Wn 'and N in sheet above is a regression line as pictured in Figure 2.2.
Example 3: Log-log Regression
In the third example, we will examine on the data that is formed in log-log scale in the XY
coordinate system. The given data is the standard sieve size (US) for the particle size
distribution of the granular material as shown in Figure 2.3.

A B C D
1 Sieve Sieve Fitted
2 No. Size Sieve Size
3 (US) (mm) (mm)
4 N S Formula :
5 4 4.750 5.142 =10^TREND(LOG10(S),LOG10(N),LOG10(A5))
8 7 2.800 2.794 Copy the formula down
9 8 2.360 2.415
10 10 2.000 1.894
11 12 1.700 1.552
12 14 1.400 1.312
13 16 1.180 1.134
14 18 1.000 0.998
15 20 0.850 0.889
16 25 0.710 0.697
17 30 0.600 0.572
18 35 0.500 0.483
19 40 0.425 0.418
20 45 0.355 0.367
21 50 0.300 0.328
22 60 0.250 0.268
23 70 0.212 0.227
24 80 0.180 0.196
25 100 0.150 0.154
26 120 0.125 0.126
27 140 0.106 0.107
28 170 0.090 0.086
29 200 0.075 0.072
30 230 0.063 0.062
31 270 0.053 0.052
32 325 0.045 0.043
33 400 0.038 0.034
34 500 0.025 0.027
Figure 2.3: Formula for Log-log data relationship

Regression in log-log data relationship can be obtained using formulas as shown in Figure 2.3,
or by adding Trendline by the following steps: click the Add Trendline > Power > Display
equation on chart and the R-squared value on the graph, the results will be shown as in
Figure 2.4. The general equation of Power Trendline is:
y = 10cXa, c and a = constant (2.5)
Slope and Y-intercept are respectively the constants a and c in equation 2.5, where the
formulas for this series are, respectively:
=SLOPE(LOG (S),LOG (N)) = -1.09022
=INTERCEPT(LOG (S),LOG (N)) = 1.36749
Figure 2.4: Chart from Figure 2.2 data series where TREND is used to
obtained the fitted points

2.5.2 Polynomial Regression
Polynomial regression is used to obtain fitted points of a number of points data which having
non-linear (curve) trend in XY coordinate system. Thus, polynomial regression consists of fitted
points of a number (n) of points data (xi, yi) to be a polynomial equation form:
f (x) = amxm + am-1xm-1 + + a1x + a0 (2.6)
In Equation 2.6, m = a positive integer and the equation is called a polynomial function of order
m. Linear regression is a special form of polynomial regression with m = 1, A = ai and B = ao.
A B
Temp Unit Mass
o
1 ( C) (gr/cm3)
2 X Y
3 4 1.00000
4 16 0.99897
5 17 0.99880
6 18 0.99862
7 19 0.99844
8 20 0.99823
9 21 0.99802
10 22 0.99780
11 23 0.99757
12 24 0.99733
13 25 0.99708
14 26 0.99682
15 27 0.99655
16 28 0.99627
17 29 0.99598
18 30 0.99568
Figure 2.5: Example of two data series for polynomial regression

Example
Figure 2.5 shows two data series from the relationship between water density and temperature,
where X-coordinate is the temperature (oC) and the Y-coordinate is the water density (g/cm3).
By plotting Y-values against X-values, it is apparently visible that the result has a curve trend.
Figure 2.6 shows chart of the plotted data from Figure 2.5 and regression lines created using
Trendline that consist of Linear and Polynomial of order 2 and 3. It is apparently seen that the
polynomial regression is much better than linear regression for this data series. The fitted points
produced from order 2 and 3 regression seemed like really coincides with the data and are only
distinguished by the value of R2, where order 3 (red curve) the value of R2 = 1. Excel 2007
provides the type of polynomial regression to order 6.
Figure 2.6: Regression analysis using Trendline for data series in Figure 2.5

2.5.3 Interpolation
In an engineering calculation, sometimes it only needs a data interpolation without resolve it into
regression analysis. Here is to find y-value that correspond to known x-value a line segment,
based on a number data (xi, yi) tabulation, where xi may continuously be increased or
decreased. Linear interpolation of y-value corresponds to x-value is:
y i +1 y i
y = yi + ( x xi ) (2.7)
xi +1 xi
where ( x xi )( x xi +1 ) 0 and x i x i +1
If y and x-axis has linear-log scale, respectively (semi-log), equation 2.7 becomes:
ln( yi +1 ) ln( yi )
ln( y) = ln( yi ) + ( x xi ) (2.8)
xi +1 xi
where ( x xi )( x xi +1 ) 0 and x i x i +1
Equation 2.8 can also be written:
y x xi y
ln = ln i + 1 (2.9)
yi x i +1 x i yi
Or,
( x x i ) /( x i +1 x i )
y
y = yi i +1 (2.10)
yi
In Excel, the interpolation value can be obtained by TREND function. Data series for this
function is a line segment formed by xi, yi and xi + 1, yi + 1.
Example 1: Linear Interpolation
Given below is two data series of X and Y-array for interpolation:

Table 2.1: Two Data Series for Linear Interpolation
Data X Data Y
0 -1
1 0
2 4
3 5
4 -1
5 -2
6 0
7 3
Plotted in chart:
Worksheets below show formulas to obtain the interpolation of the data in Table 2.1, using
TREND function and Equation 2.7, respectively.
Using TREND function:

Using Equation 2.7:

The interpolation values are in column D, which correspond to new values entered in column C.
The y-value that correspond to a x-value that lies between 0 - 1 in cell C3 is shown in cell D3;
the y-value that correspond to a x-value that lies between 1 - 2 in cell C4 is shown in cell D4,
and so forth.
TREND function and Equation 2.7 can also be used to interpolate x-values from known y-
values, by replacing the variable x with y and vice versa in Equation 2.7, or by switching x and
y-argument in TREND function. The formulas are shown in worksheet above at row 12, columns
C and D.
Example 2: Logarithmic Interpolation
A sands soil sample of 1,000 grams weight is taken and placed in the sieving machine for a
sieve analysis test. The results are presented in Table 2.2, showing sieve size used and percent
finer. The percent finer (or percent passing) is calculated based on cumulative weight of soil
retained on each sieve. Chart of log grain sizes versus percent finer of soil is presented in
Figure 2.7.
Table 2.2: Sieve Analysis Results
A B
1 Sieve Size Finer
2 (mm) %
3 9.50 100.00
4 4.75 97.80
5 2.36 95.48
6 1.18 88.50
7 0.60 70.60
8 0.30 24.10
9 0.15 2.40
10 0.00

The nextt step is to obtain grain size
s distributtion in perce
entage of the
e weight of th
he sample, for
f an
example to obtain grrain size value correspo
onds to 60 percent
p finerr or D60. Th
he correspon
nding
value can be obtain
ned using Eq
quation 2.10
0. Formulass for semilog
g interpolation are show
wn in
workshee
et as follows
s:
n chart:
Plotted in
Fig
gure 2.7: Se
emilog interpolation forr the data of
o Table 2.2

gspot.co.id
For an exercise, a file contains this interpolation example can be seen in Interpolation.xls on the
CD included in this book.
Notes:
There is however more-more flexible way to obtain data interpolation which to resolve
them with a macro (custom function), using Equation 2.7 or 2.10. When entered into
macro, formulas can be shortened and does not depend on a line segment where a
known value is placed. Discussion on the custom function macro will be in the following
chapter in the topic of user defined function.
The interpretation on sieve analysis result includes to determine how well the gradation
of sands soil. Gradation of sands or grained particles (silt to gravel) can be measured by
the coefficient of uniformity (CU) and the curvature coefficient (CC) as defined below.
D60 D 302
CU = C =
D10 , C D 60 D10
The greater the value of CU, the greater grain distribution range, where well-graded
particles usually has CC value between 1 and 3.
2.5.4 Statistical Data
Histogram and Cumulative Distribution
In a statistical data presentation, it is convenient to make the data from observation into
histogram to see the data distribution. In histogram, the observed data plotted against its
frequency distribution and thus, we can get a visual summary and give quick impression on the
observed data.
Example
The following is concrete compressive strength data selected randomly from 40 samples
obtained from the test of characteristics compressive strength of concrete:

Concrete compressive strength test
result (kg/cm2)
376.28 399.30 375.41 370.41
394.43 409.47 393.53 371.43
387.19 374.34 386.30 372.44
390.35 392.29 389.46 390.30
401.27 386.78 390.37 384.82
379.95 361.23 379.08 362.28
380.97 402.90 380.10 400.86
389.33 368.22 388.44 366.35
384.64 399.20 383.76 407.46
396.34 365.36 398.83 363.51
To simplify data distribution, the result data divided into several classes within interval (taken) =
5 kg/cm2, and then tabulated in Table 2.3. The scores from each class of data is called class
frequency, which is represented by height of the bar. To ease data distribution process into
classes can be done by using a function macro or custom function. The discussion on custom
function will be in the next chapter in the topic of user defined function.
Table 2.3: Data distribution to create a histogram
Class
Interval Mid Point Frequency Percent Cumulative
2 2
kg/cm kg/cm N (%) (%)
360-365 362.5 3 7.5 7.5

365-370 367.5 3 7.5 15.0
370-375 372.5 4 10.0 25.0
375-380 377.5 4 10.0 35.0
380-385 382.5 5 12.5 47.5
385-390 387.5 6 15.0 62.5
390-395 392.5 6 15.0 77.5
395-400 397.5 4 10.0 87.5
400-405 402.5 3 7.5 95.0
405-410 407.5 2 5.0 100.0
40 100.0

Histogram chart of Figure 2.8 is made based on Table 2.3 data. To create a histogram take the
following steps: click Insert tab > Column > Column clusterd > Select Data > Add > Series
Values: select range "Percentage" > OK. Click Edit for Horizontal Axis Labels > insert range
"Mid Point".
To create cumulative frequency chart as Figure 2.9, select the type of XY scatter chart with
smoothlines and markers > Select Data > Add > Series X values: select range "Mid Point" >
Y-values: select range "Cumulative".
Figure 2.8: Histogram of concrete compressive strength
Figure 2.9: Cumulative frequency of concrete compressive strength

Range, Mean dan Standard Deviation
The classification or data distribution is generally defined into three statistical parameters, which
are, range, mean and standard deviation. Range is the difference between the lowest and the
highest values of data. From the concrete compressive strength test result, there is a range of
48.24, from the lowest 361.23 to the highest 409.47 in kg/cm2 units. The mean or average is the
sum of the data values divided by the total number of data (n). From the test, the mean is
obtained = 384.87 kg/cm2. If a set of data values are divided into M classes, the mean formula is
written as follows:
m x
i =1
i i
x= (2.11)
n
Where,
mi = frequency of class interval i

xi = mid point of class interval i
Mean is a value that commonly used to represent values of a set of data. In consistent with this,
several measurements are made to state deviation, and one of them is the formulation of mean
square deviation (MSD). The MSD is also called the variance, and the formula is:
m (x
i =1
i i x) 2
(2.12)
n
Another measurement to be used is the standard deviation (S), which is the square root of
MSD:
m (x
i =1
i i x) 2
S= (2.13)
n
Through Equation 2:13, the standard deviation from the test result as in Table 2.3 is calculated
= 12:49 kg/cm2, and with divisor of (n - 1), calculated from the sample, MSD and S is

respectivvely calculatted 394.74 and 16.12 kg/cm2. Th
he standard deviation is
i a measure of
dispersio
on of a set of
o data value
es from its mean;
m where a big or sm
mall of the va
alue indicate
es the
quality off job execution.
Normal Distributio
on
From the
e histogram and measurement of statistical
s parameters ab
bove, it can now be dra
awn a
normal distribution cu
urve. The no
ormal distribution functio
on is,
2
1 x x

1 2 S
f ( x, S , x ) = e (2.14)
S 2
For log-n
normal distrib
bution, it can ed by replacing x with ln((x) in Equatio
n be obtaine on 2.14.
Figure 2..10 and 2:11 mal distribution curve tha

1 show norm at fits with th
he percent frequency
f da
ata in
Table 2.3
3. Figure 2.12 shows the
e formulas used in the presentation of the test re
esult.
10: Normal distribution

Figure 2.1 n curve thatt fits to the data of Tab
ble 2.3

gspot.co.id
100
90
Cumulative frequency (%)

80
70
60
Data
50
Normal
40
30
20
10
0
360 370 380 390 400 410
Concrete compressive strength (kg/cm2)
Figure 2.11: Cumulative normal distribution curve that fits to data of Table 2.3
A B C D E F G
Class
1 Interval Mid Point Frequency Percent Cumulative Normal Normal
2 2
2 kg/cm Kg/cm (%) (%) Cum.
3 k N f
4 360-365 362.5 3 7.5 7.5 3.2 3.7
5 365-370 367.5 3 7.5 15.0 6.1 8.2
6 370-375 372.5 4 10.0 25.0 9.8 16.1
7 375-380 377.5 4 10.0 35.0 13.4 27.8
8 380-385 382.5 5 12.5 47.5 15.7 42.5
9 385-390 387.5 6 15.0 62.5 15.6 58.3
10 390-395 392.5 6 15.0 77.5 13.2 72.9
11 395-400 397.5 4 10.0 87.5 9.6 84.4
12 400-405 402.5 3 7.5 95.0 5.9 92.1
13 405-410 407.5 2 5.0 100.0 3.1 96.5
14 40 100.0
15
16 Mean = 384.87 kg/cm2
17 S= 12.49 kg/cm2

F G
1 Normal Normal
2 Cumulative
3
4 =500/($E$17*SQRT(2*PI()))*EXP(-0.5*((k-$E$16)/$E$17)^2) =NORMDIST(k,$E$16,$E$17,TRUE)*100
5 =500/($E$17*SQRT(2*PI()))*EXP(-0.5*((k-$E$16)/$E$17)^2) =NORMDIST(k,$E$16,$E$17,TRUE)*100
D E F
15
16 Mean = =SUMPRODUCT(k,n)/SUM(n) kg/cm2
17 S= =SQRT(SUMPRODUCT(n,(k-E16)^2)/SUM(n)) kg/cm2
Figure 2.12: Data presentation from the test result to obtain normal
distribution and the formula used in spreadsheet
Formulas:
Mean : =SUMPRODUCT(k,n)/SUM(n)
S : =SQRT(SUMPRODUCT(n,(k-Mean)^2)/SUM(n))
Normal : =500/(S*SQRT(2*PI()))*EXP(-0.5*((k-Mean)/S)^2)
Normal Cumulative : =NORMDIST(k,Mean,S,TRUE)*100
Mean and standard deviation are respectively, calculated using Equation 2:11 and
2:13. The actual scale for normal function (Equation 2.14) is determined from the
interval used, in this case were taken 5 kg/cm2. Multiply then by 100 to get the
depiction in percent (%).
Excel has also the other ways to make distributions, such as Exponential, Weibull,
Gamma and Poisson.

2.5.5 Circular Reference
When formulas are connected with their own cell, either directly or indirectly, then this
condition is called a circular reference. Excel by default can not process a formula that
contains a circular reference until we select Enable Iterative calculation by the
following steps: Office Button > Excel Options > Formulas. Iteration is the repetition
of numerical calculation until a condition is met. The process is executed based on
iterations number entered in the textbox Maximum iterations and Maximum Change
for the maximum change resulting from the last two iterations. Excel will stop counting
when one of the conditions is met, or whichever comes first.
Selecttoallow
iterativecalculation
Cell with a circular reference in a spreadsheet as shown below:
Example 1
A B Result : 1 x iteration A B
1 4 1 4
2 =B2 =(A1+A2)/3 2 0 1.33333

5 x iterations A B
1 4
2 1.97531 1.99177
100 x iterations A B
1 4
2 2 2
In short, it can be explained herein about the iteration method - to describe the way of
Excel executes Example 1. Formulas in cell A2 and B2 are interrelated and have the
form: x = f(x) and calculated by the following equation:
xn+1 = f(xn),
where: xn = predicted value of x at iteration nth
xn + 1 = predicted new value of x at iteration nth + 1
n = number of iterations
Formula in cell sel A2 and B2 can be written: xn+1 = (4 + xn)/3.
At first iteration (n =1), x1 is taken = 0 x2 = (4 + 0)/3 = 1.33333
At n = 2, x2 is taken = 1.33333 x3 =1.77778
At n = 5, x5 = 1.97531 x6 = 1.99177
And so on.
Excel will do this process based on the number of iterations (n) and the maximum change
x n +1 x n that you specify, and will stop if one is met. This analogy is found in the process of
Newton Raphson Method function in Chapter 3, but of course, the method is different.

Application of Excel iteration can also be applied to solve system of linear equations, to obtain
unknown variables of equations. For an example, to obtain values of x, y which satisfy the
equations below:
2x + 3y = 7
3x - 2y = 4
The solution is firstly to re-write each of equation to form x and y:
2x + 3y = 7 x = (7 3y)/2; y = (7 2x)/3
3x 2y = 4 x = (4 + 2y)/3; y = (4 3x)/2
Next, convert the formed x and y into formulas and placed them in cell B2 until C3 as shown in
the worksheet below. In order to perform iteration, the cell must be linked which is in this
operation the formula in cell B3 and C2 are to be interrelated.
10 x iterations
11 x iterations
Values of x and y that satisfy both equations are obtained at iteration 11, which are:
x = 2 and y = 1.
Example 2

A B C 1 x iteration A B C
1 0 10 10 1 0 10 10
2 10 =(B1+B3+A2+C2)/4 =B2 2 10 5 5
3 10 =(B2+B4+A3+C3)/4 =B3 3 10 3.75 3.75
4 0 0 0 4 0 0 0
5 x iterations A B C
1 0 10 10
2 10 8.64128 8.64128
3 10 6.18032 6.18032
4 0 0 0
100 x iterations A B C
1 0 10 10
2 10 8.75 8.75
3 10 6.25 6.25
4 0 0 0
In the cell B2 formula, the initial value for B3 and C2 are taken = 0. The value of B2 will be used
as the initial value in the cell B3 formula, where C3 = 0. The values obtained are then used as a
predicted new value in the next iteration.
In its application, for examples, the Excel iterative calculation or circular reference can be used
in the solution of linear simultaneous equations, slope stability analysis using Bishop method
where the safety factor is on both sides of the equation, and the solution for seepage below dam
using finite difference method.

Chapter 3:
CREATING MACRO
3.1 Function Procedure

In addition to the functions that already exist in Excel, it is often to create a function to perform
calculations which is determined by the user. Such function is called user-defined (UD) function.
This function will be added by Excel into a collection of functions and can be used as well as a
built-in function.
Some UD functions are made by combining formulas and built-in functions. UD function is made
in macrosheet or using VBA. Macrosheet is a sheet where macro is created by XLM language,
the original Excel macro language that is still maintained, at least until the Excel 2007. Figure
3.1 shows the UD functions to obtain cosine of an angle in macrosheet and VBA.
A
1 COSD Function Cosd(x)
2 =RESULT(1) Cons Pi = 3.141592655
3 =ARGUMENT("X",1) Cosd = Cos(Pi * x/180)
4 =RETURN(COS(PI()*X/180)) End Function
Figure 3.1: Function in Macrosheet (left) and VBA (right)
Creating Function in Macrosheet
To insert a macrosheet in worksheet, right-click on the sheet tab to display a command menu,
click Insert > MS. Excel 4.0 Macro > OK, to insert a macrosheet named Macro1 as shown
below:

This sheet is seemin
ngly similar to
t a worksheet, but there are colum
mns that are
e wide enough to
write a function. The
e first row iss to write th
he name of function, fo
or example, COSD in macro
m
example in Figure 3.1.
3 The seccond row is to define the type of ou
utput refers to Table 3.1, for
example, number, te
ext or an arrray, using RESULT
R sta
atement and
d followed by
b the data type.
Next, is to define va
ariables or in
nput data ussing ARGUM
MENT statement, follow
wed by the name
n
a type.
and data
User hass to determin

ne formulas used in the operation of function off the input va
ariables, which is
defined by
b ARGUME
ENT. A varia
able can refe
er a name orr cell referen
nce of worksheet, such as
a "x"
of COSD
D for a cell in
n the worksh
heet. A formula can be given
g a nam
me to simplifyy writing the
e next
formulas. The process in writing
g a function is ended byy RETURN statement,
s w
which returns the
function value.
v Below
w is an exam
mple of macrrosheet maccro to calcula
ate volume of
o a cube:

gspot.co.id
MVOL fu
unction is used in the wo
orksheet as follows:
f
Tab
ble 3.1: Data
a Types in Macrosheet
M t
Vallue Da
ata Type
1 Nu
umber
2 Te
ext
4 Logical
8 Re
eference
16 Errror
64 Arrray
Creating
g Function
n in Visual Basic for Applicatio
A n (VBA)
UD functtion is made
e in a modu
ule that inse
erted through
h Insert tab
b > Module in Visual Basic
B
Editor which
w has be
een shown previously
p in
n Chapter 1. Function procedure
p is begin with word
function, followed by
y its name an
nd ends with
h the word End
E Function
n. The type of a function
n can
be Single, Double or
o all data types
t of Variant, depen
nding on wh
hat it would be produce
ed, a
number or
o an array. Data types in
i VBA are presented
p in Table 3.2.
The selection on these data type

es in addition
n to define a function an
nd argumentts, or variables, it
elated to the
is also re e allocation of computer memory which
w depen
nds on the ra
ange of the data
type (the
e greater the
e value, the greater
g mem
mory allocatiion). If an arrgument is not
n defined, it will
be classiified as a Va
ariant and co
onsumes me
emory more
e than any otther data typ
pe. Defining data
type in a programmin
ng is highly recommende
r ed because it saves com
mputer mem
mory consum
mption
elerates the process of VBA
and acce V in recog
gnizing the data
d type ass well, or in the
t other wo
ord to
acceleratte VBA to ex
xecute a macro.

gspot.co.id
Table 3.1: Data Types in VBA
Data Type Range

Integer -32,768 to 32,767
Long (long integer) -2,147,483,648 to 2,147,483,647
Single (single precision) -3.402823E38 to -1.401298E-45 ( - values)
1.401298E-45 to 3.402823E38 (+ values)
Double (double precision) -1.79769313486231E308 to -4.94065645841247E-324 (-)
4.94065645841247E-324 to 1.79769313486232E308 (+)
Currency -922,337,203,685,477.5808 to 922,337,203,685,477.5807
String (variable length) 0 to approx.2 billion characters
String (fixed length) 1 to approx. 65,400 characters
Byte 0 to 255
Boolean (logika) TRUE or FALSE
Date 1 Januari 100 to 31 Desember 9999
Object Any object reference
Variant (with numbers) Any numeric value up to the range of Double
Variant (with characters) 0 to approx.2 billion characters
Such as Excel, VBA also has many built-in functions that have names similar to Excel built in
functions, say, Cos function in Figure 3.1 or Abs to find the absolute value of a number. Some
Excel built-in functions can also be used in VBA by using the prefix Application. Figure 3.2 is an
example where Excel built-in functions are used in UD functions MLIN to obtain linear
regression parameters. MLIN returns an array result (multiple result) which is the slope of the
line and the intersection with Y-axis (Y-intercept).
Function MLIN(Y, X) As Variant

Slope = Application.Slope(Y, X)
Intercept = Application.Intercept(Y, X)
MLIN = Array(Slope, Intercept)
End Function
Figure 3.1: Excel built-in functions in VBA

Notes:
Module with Option Explicit requires all variables are to be declared. If not, VBA will
show an error message "Variable not defined". To declare a variable use Dim, ReDim,
Private, Public and Static statement.
Example
Option Explicit
Sub OptEx1()
Dim myVar as Integer
MyInt = 5 'error: variable not defined
myVar = 10
End Sub
Dim is to declare a size of variable that already set at the beginning of the program,
while ReDim is to change the size of variable that could be changed during the
execution of program. Variables declared by Dim or Redim at procedure level (between
Sub End Sub statements) are available only to the procedure where the variables are
declared.
Example
Sub OptEx2()
Dim myVar as Integer
Dim myString(1 to 3) as String
Dim n As Integer
n = 3 'an input
ReDim MyInt(n) As Integer

MyInt(1) = 1
MyInt(2) = 3
MyInt(3) = 5
myVar = 10
myString(3) = "Hello"
End Sub

Use Type..End Type with example below to define a user-defned data type.
Type FamilyTree
Name As String
Address As String
Phone As Long
Numsiblings As Byte
End Type
Sub GetRecord()
Dim TheFirst As FamilyTree
TheFirst.Name = "Johan"
TheFirst.Address = "Street 123"
TheFirst.Phone = 1234
TheFirst.Numsiblings = 3
End Sub
3.2 Sub Procedure

Macro that previously used to create a function can also be used to make an engineering
calculation program to solve more complex problem, such as structural analysis. In VBA, such
program is created in the Sub procedure, that is procedure with syntax as briefly discussed in
Chapter 1. In general, the process of input - output data from a Sub procedure is:
1. Read input data from the worksheet
2. To execute data input by VBA
3. Print the results to the worksheet
In Excel-VBA programming, the value of input data is taken from the value of the cell or range of
cells on worksheet or a Range object with value as shown by its Value property. So typing in the
value in the cell in the worksheet with numbers, text, and so on actually is to set the value on
Value property of these cells. Cells property is commonly used to express the Range object with
the syntax Cells (row_index, column_index), i.e. cell (1,1) to return a Range object of cell A1.

Here is an example to set a number to a cell in VBA languange, which sets the value of cell A1
= 10:
Worksheets(Sheet1).Cells(1,1).Value = 10
Value is the default property of the Range object, therefore negligible or means no need to be
written again. If no object qualifier or an objects before period (Worksheets () above), then a
range of cells are on the active sheet, in this case from which the macro runs, for example via a
command button on the worksheet. Use an object qualifier before the range of cells if the data is
not located on the active sheet with the writing of Worksheet ("sheet name").
To complete this topic of discussion, the following is a simple example of Sub procedure named
Pro_Sub1:
Sub Pro_Sub1()
Dim Length As Double
Dim Width As Double
Dim Area As Double
Length = Cells(2, 1)
Width = Cells(3, 1)
'Determine Area
Area = Length * Width
Cells(4, 1) = Area
End Sub
The above procedure is to calculate the area of a rectangle which is the multiplication of the
length and width of its sides. The variables are declared using the Dim statement and the data
type is a number with double precision (Double). Cells property is created without an object
qualifier, therefore, the range of cells is on the active sheet, which the active sheet (opened)
when the code program is executed.
Sign apostrophe (') is a sign where program ignores the comments in it. This sign is used to
insert an explanation. The default color of the text is green.
The process of input-output of Pro_Sub1 procedure can be explained as below:

If the pro
ocedure in th
he above fig
gure is execu
uted, value of
o cell A4 orr Cells (4,1) = 30. To run the
macro, perform
p a seq
quence of stteps as desccribed in Cha
apter 1.11.
In creatin
ng large sca
ale applicatio
on program, normally, the
t program
m is divided into several sub-
programss or procedu
ures Sub and/or with so
ome function
ns that each handles a specific
s taskk then
they calle
ed with Call statement or the name
e of function
n. The advan
ntages of th
his division are
a to
programming become more
make a p m focused, more rea
adable and iff an error occcurred it ca
an be
quickly addressed.
a Chapters
C 6 to
o 10 show th
he exampless of such com
mputer programming.
3.3 C
Control Sttructures
s
A progra
am is create
ed to be able to handle
e any changes either fro
om the data
a input or to
o that
occurred during prog
gram executiion. Therefo
ore, there are
e necessary control structures in ord
der to
ogram rema
make pro ains on its specified
s ste
eps. Such controls
c are looping and
d branching or a
combinattion of both.
3.3.1 Looping
Do While
e...Loop or DoLoop While
W
Repeats a set of sttatements while

w a cond
dition is TR
RUE. If cond
dition becom
mes FALSE
E, the
program will exit the looping and
d move to the
e next code..

gspot.co.id
Syntax:
Do While condition
[statements]
[Exit Do]
[statements]
Loop
Example
Sub exLoop1()
i = 1
Do While i <= 5
Cells(1 + i, 1) = i
i = i + 1
Loop
End Sub
Do UntilLoop or DoLoop Until
Repeats a set of statements until a condition becomes TRUE. If a condition becomes FALSE,
the looping will be discontinued.
Syntax:
Do Until condition
[Statements]
[Exit Do]
[Statements]
Loop
Example

Sub exLoop2()
i = 1
Do Until i >= 5
Cells(1 + i, 1) = i
i = i + 1
Loop
End Sub
ForNext
Repeats a set of statements a given number of times
Syntax:
For counter = start To end [Step step]

[statements]
[Exit For]
[statements]
Next [counter
Example
Sub exLoop3()
j = 1
For i = 1 To 5
Cells(1 + i, 1) = j
j = j + 1
Next i
End Sub
For EachNext
Repeats a set of statements a number of elements in an array. It helps if the number is

unknown. For example, to be used in a function below:
Syntax:
For Each element In array

[statements]
[Exit For]

[statements]
Next [element]
Example
Function exLoop4(Yarray, Xarray) As Variant

Dim n As Integer, x As Single, y As Single
n = 0
x = 0
y = 0
For Each c In Xarray

n = n + 1
x = x + Xarray(n)
y = y + Yarray(n)
Next
a = x * x
b = y * x
exLoop4 = Array(a, b)
End Function
3.3.2 Branching
Select Case
Executes several set of statements depends on the value of an expression
Syntax:
Select Case testexpression

[Case expressionlist-n
[statements-n]]
[Case Else
[elsestatements]]
End Select
Example

Sub exLoop5()
Dim Number
Number = Cells(1, 1) 'a number
Select Case Number
Case 1 To 5
Cells(2, 1) = "Between 1 and 5"
Case 6, 7, 8
Cells(2, 1) = "Between 6 and 8"
Case Else
Cells(2, 1) = "Greater than 8"
End Select
IfThenElse
Executes a set of statements because of certain condition, when it is not appropriate then move
to another state of condition.
Syntax:
If condition (1) Then

statement(1)
ElseIf condition(2) Then
statement(2)
Else
statement(3)
End If
Branching is usually nested in a looping. If the condition is met, the program will exit the looping
with Exit command. If...ThenElse branching can also be combined with And, Or, or Not
statement to state the expression or the relationship between two expressions.
Syntax:

Foror Do Until
statement(1)
If statement(1) And or Or expression(2) Then
statement(2)
Exit For or Exit Do or Exit Function or Exit Sub
End If
Nextor Loop
GotoGosubReturn
This branching is used in a procedure consisting of one or several blocks of statements

(subroutines) that is made for particular task or calculation. The block is then called upon to do
the job with GOSUB statements and returned with RETURN statement. GoTo statement will
branch to a code line, subroutines, or line label, and then continue to the next code.
Example
Sub exBranch()
Dim myNum As Double
myNum = Cells(2, 1)
If myNum >= 0 Then
GoSub MyRoutine
GoTo 10
Else 'means myNum < 0
GoTo 20
End If
MyRoutine:
myNum = Sqr(myNum) 'determine square root of number >=0
Return
10 Cells(3, 1) = myNum
Exit Sub
20 Cells(3, 1) = "#NUM!" 'error message
End Sub

CallSub
Calls a procedure from another procedure
Syntax
[Call] name [argument]
Example 1: Calling Sub procedure from function procedure
Sub Area(x, y, L)
L = x * y
End Sub
Function Volume(a, b, t) As Single

Call Area(a, b, L)
Volume = L * t
End Function
Example 2: Calling Sub procedure from another Sub procedure.
Sub Area(x, y, L)
L = x * y
End Sub
Sub Pro_Vol()
Dim a As Single, b As Single, t As Single
Dim volume As Single
a = Cells(2, 2)
b = Cells(3, 2)
t = Cells(4, 2)
Call Area(a, b, L)
volume = (L * t)
Cells(5, 2) = volume
End Sub

In addition to looping and branching defined above, there is also useful control structure that is
With statement.
With
Execute a series of statements on the same object
Syntax:
With object
[statements]
End With
Example:
With MyObject
.HorizontalAlignment = xlCenter
.VerticalAlignment = xlCenter
.Position = xlLabelPositionRight
.Orientation = xlHorizontal
With .Font
.Name = "Arial"
.FontStyle = "Regular"
.Size = 8
End With
End With
3.4 User Defined Function Problems
Problem 1
Create a macro to determine rectangular area where area = length x width.
Solution
If a and b express length and width respectively, area in macro can be expressed as a function
of (a, b) = a x b. If the name of the function is Area, the expression is:
Area (a, b) = a x b.

To solve this problem, we will use macrosheet languange. Insert a macrosheet by right-clicking
on the worksheet Sheet tab > Insert > MS. Excel 4.0 Macro > OK. Type the macro as follows:
Afterward, we have to define the function name so that it can be recognized by Excel and
incorporated it into Excel functions. The steps are as follows:
1. Click cell A1 where the name of the function (Area) is placed.
2. Click the Formula tab > Define Name to show the New Name dialog box as below.
Enter a macro name in the name text box. Macro name will automatically be referred to
cell A1 with step 1 above.
3. Select Function in the option button, so the function (by default) is categorized as User
Defined or you can choose a category for your function.
4. Click OK.
Now back to the worksheet. Area function is written as follows:

Cells A2 and B2 are
e respectivelly named a and b which
h carried ou
ut by follow step
s No. 1 and
a 2
above through active
e sheet. In th
his step, the
ere is no option (option button)
b for Macro
M catego
ory in
me dialog bo
New Nam ox.
Problem
m2
Find the area betwee

en a straightt line y = mxx + c and the
e X-axis, betw
ween x = a and
a x = b (c 0).
This will be solved us
sing integrall with formula:
Area = (mx + c)dxx

a
In macro
o, area can be
b expressed
d as:
.mb2 + cb) - (ma2 + ca)

function(m,c,a,b) = ( c
Now, bassed on the area

a formulation above, find
f the area
a between y = x + 2, X-a
axis, betwee
en a =
-1, b = 3.. The name of function iss INAREA.
Solution
n
sheet looks like below:

The maccro in macros

gspot.co.id
Problem 2 in chart:
Before applying
a the function, pe
erform step no.1 to 4 as
a in Proble
em 1. The use of function in
workshee
et as follows
s:
Problem
m3
Create a VBA macro

o to assign letter grade for a course
e with score range 0 - 100 as follow
ws, A:
nting scores from 80 - 10
represen 00, B: 65 - 79,
7 C: 51 - 64
4, D: 40 - 50
0 and F: for scores
s below
w 40.
Macro wiill return literral grade witth sign + and
d at given numerical value.

gspot.co.id
Solution
Function VBGrade(g, a, b, c, d) As String

Select Case g
Case Is >= 2 * (100 - a) / 3 + a
VBGrade = "A+"
Case Is >= (100 - a) / 3 + a
VBGrade = "A"
Case Is >= a
VBGrade = "A-"
Case Is >= 2 * (a - b) / 3 + b
VBGrade = "B+"
Case Is >= (a - b) / 3 + b
VBGrade = "B"
Case Is >= b
VBGrade = "B-"
Case Is >= 2 * (b - c) / 3 + c
VBGrade = "C+"
Case Is >= (b - c) / 3 + c
VBGrade = "C"
Case Is >= c
VBGrade = "C-"
Case Is >= 2 * (c - d) / 3 + d
VBGrade = "D+"
Case Is >= (c - d) / 3 + d
VBGrade = "D"
Case Is >= d
VBGrade = "D-"
Case Else
VBGrade = "F"
End Select
End Function
Functions that created in VBA will automatically be entered into Excel functions collection. The
use of the function in worksheet is as follows:

If the teacher wants or course sccoring, it is simply don
s to apply a different standard fo ne by
g value of the
changing e function arrguments.
Problem
m4
Create a VBA macro

o to find vertical stress distribution
d b
below a foundation with
h elastic equ
uation
from Bou
ussinesq - Newmark (19
935). The equation to get the vertical stress () below
b a corn
ner of
a rectang
gular area b x l at depth y as shown in Figure 3.3 is:
1 2MN V V + 1 2 MN
N V
y =0 . + tan 1
4 V + V1
V V V1
where,
M = b/y
N = l/y
V = M2 + N2 + 1
N)2
V1 = (MN
b dan l iss length of re

ectangular siides
If V1> V, then the arc

ctan in the eq
quation need
ds to be add
ded with .

gspot.co.id
Figure
e 3.3: Vertic
cal stress below
b a corn
ner of a recttangular are
ea
Solution
n
Function
n Bous(q, b, l, y) As
A Double
Const pi
i = 3.1415
593
n v, vs, rad, at
Dim m, n,
m = b / y
n = l / y
v = m ^ 2 + n ^ 2 + 1
vs = (m * n) ^ 2
rad = (2 * m * n * Sqr(
(v)) / (v - vs)
If vs
v > v The
en
at = Atn(
(rad) + pi
Else
e
at = Atn(
(rad)
End If
Bous 4 * pi) * (((2 * m * n * Sqr(

s = q / (4 (v) / (v + vs)) * ((v
( + 1) / v))
+ at)
End Func
ction
antage to co
The adva onvert an eq
quation or fo
ormula into a function iss that it grea
atly simplifie
es the
calculatio
on. The abo
ove equation
n with a con
ndition V1 > V certainlyy would be very
v lengthyy and

gspot.co.id
requires a lot attention if built in
nto worksheet formula. It is more easily
e solved
d by typing it into
VBA code.
Now, givven a founda ure of 10 ton/m2.

ation with a size of 3 x 3 m is subjjected to uniform pressu
The vertical stress distribution
d of each inte
erval depth of 0.5 mete
ers below th
he center of
o the
on is shown in the workssheet as belo
foundatio ow:
Boussinesq - Newmark equation is intended to

o find vertica
al stress below a cornerr of a
re
ectangular area,
a therefo oundation off 3 x 3 m is obtained 4 area
ore, below a center a fo
w the same
with e size, with 1.5 x 1.5 m dimension. Furthermore, for the co
ontribution off four
co
orners it sha
all be calcula
ated = 4 y.
D
Division by 0 is not allowed in the calculation, the
erefore the value
v ed in cell A8 shall
entere
be
e an approa
ach = 0 (take
e, for examplle, y = 0.000
01).

gspot.co.id
Chart of stress versus depth is shown as below:
Chart is created through Insert tab> Scatters with smoothlines and markers > Select Data >
Add and then select x and y-data. To modify the chart click once on chart area then click Chart
Tools tab > Chart Layout > Layout 1, and then add chart titles by clicking once on the title,
select Edit Text to type. The depth orientation is made to increase downward by right clicking
the Y-axis > Format Axis > Axis Option > Values in reverse order.
Problem 5
Find root x in equation: ax4 + bx3 + cx2 + dx + c = 0
Solution
This problem will be solved numerically by Newton-Raphson method. The process of calculation
is done in an iterative technique, which root x is gradually approached. Function of equation
depicted in chart in Figure 3.4 and exemplified intercepts the X-axis at the point p and f (x) = 0

or (p, 0). In the first step, it is taken x1-value to get tangent line at y = f(x) which intercepts the X-
axis at the point (x2,0), where x2 is new approach value to p. The next approaching process is
done using x2-value to get x3-value and so on until xn p.
y
y = f(x)
y f(xn) = f(xn)(x xn)
(xn,f(xn))
(p,0) (xn+1,0) x
Thus, by induction it obtains tangent line: y - f(xn) = f '(xn) (x - xn) which intercepts X-axis at point
(xn + 1,0), where:
f ( xn ) )
xn + 1 = xn
f ' ( xn )
Macro
Function NRM(a, b, c, d, e) As Double

Dim funct, deriv, x, xn
Dim n As Integer, i As Integer
n = 1000
i = 1

x = 100
Do
funct = (a
( * x ^ 4)
) + (b * x ^ 3) + (c
c * x ^ 2) _
+ (d * x) + e
deriv = (4
( * a * x ^ 3) + (3
3 * b * x ^ 2) + _
(2 * c * x) + d
xn = x - (funct / deriv)
d
s(xn - x) < 0.00001 Then
If funct = 0 Or Abs
NRM = xn
Exit Function
End If
i = i + 1
x = xn
Loop
p Until i = n
NRM = "Diverg
gen!"
End Func
ction
Function NRM (stands for Newtton-Raphson

n Method) iss now used to find the value of x in the
equation below:
-2.79x3 - 2:11 x2 + 24.51x

2 + 36.6
68 = 0.
This equation form can

c be found
d in the geo
otechnical en
ngineering to
o find the em
mbedment depth
d
of retainiing wall with
h Free-Earth
h Support Method.
M The
e use of the function in worksheet is as
below:

gspot.co.id
Problem 6
Table X is a number of random numeric data from 0 - 100. Create macro to calculate the
amount of data in Table X with a range of values from Y to Z.
A B C D E F G
1 Tabel X
2 77.58 73.27 22.21 70.40 81.40 52.57 98.34
3 68.58 21.86 62.37 37.61 51.24 12.41 18.55
4 60.92 22.26 35.01 23.76 63.34 32.81 83.16
5 98.27 63.85 30.97 34.98 64.44 70.49 17.96
6 71.23 79.63 74.21 100.00 37.95 3.94 95.27
7 10.62 20.69 5.82 74.81 95.09 37.48 0.74
Solution
Macro
Function DNUM(Xarray, Y, Z) As Integer

Dim n As Integer, Num As Integer
Dim c
n = 0
For Each c In Xarray
n = n + 1
If Xarray(n) >= Y And Xarray(n) <= Z Then Num = Num + 1
Next
DNUM = Num
End Function
The use of the function in worksheet is as below:
I J K Formula at Column K
2 0.0 50.0 19 =DNUM(X,I2,J2)
3 50.0 100.0 23 =DNUM(X,I3,J3)
4 0.0 100.0 42 =DNUM(X,I4,J4)

Argument X is an array argument, in this example is the name that has been defined for the
table data in a range A2 to G7. Argument array can be also written in the following manner,
range(A2:G7) or a certain range, for example range(B2:C7), range(A7:G7) and so on. One of
the simplicities of using function in Excel because it is well supported by the spreadsheet to
create variable names, writing function and put the result at once.
Problem 7
Create a macro for linear interpolation example in Chapter 2.5.3. Use Equation 2.7 to solve.
Solution
Macro
Function Linter(Xrange, Yrange, nval) As Double

Dim i As Integer
Dim c, x
x = Xrange(1)
i = 1
For Each c In Xrange
If nval = c Then
Linter = Yrange(i)
Exit Function
'find a new value satisfying the given condition
ElseIf (nval - c) * (nval - x) <= 0 And c <> x Then
'Eq. 2.7:
Linter = Yrange(i - 1) + (Yrange(i) - Yrange(i - 1)) _
* (nval - x) / (c - x)
Exit Function
End If
x = c
i = i + 1
Next
End Function
The use of the function in worksheet is as below:

Interpola
ation in Chartt:
EnterNewValu
ueatthisrange
Yvalue
Example
Problem
m8
Create a macro for semilog

s interrpolation exa apter 2.5.3. Use Equatio
ample in Cha on 2.10 to so
olve.
Solution
n
Macro
Function
n Loginter
r(Xrange, Yrange,
Y nv
val) As Dou
uble
Dim i As Inte
eger
Dim c, x
x = Xrange(1)
i = 1

gspot.co.id
For Each c In
n Xrange
If nval = c Then
Login
nter = Yran
nge(i)
Exit Function
'find a new
n value satisfying
s g the given
n conditio
on
ElseIf (n nd c <> x Then
nval - c) * (nval - x) <= 0 An
'Eq 2.10:
2
Loginter = Yrange(i - 1) * (Y
Yrange(i) / Yrange(i
i -1)) _
^ ((nval
l - x) / (c
c - x))
Exit Function
End If
x = c
i = i + 1
Next
t
End Func
ction
The use of the function in worksh

heet is as be
elow:
ation in Chartt:
Interpola

gspot.co.id
EnterNewValueatthisrange
Notes:
Character "c" and "r" can not be used for variable names in macrosheet, because C and
R stand for column and row, respectively. However, you can define a variable names
such as "c_" or "r_", "cx" or "rx" and so on.
The Problem 1 to 8 can be seen in User Defined Function.xls file on the CD included in
this book.
3.5 Structure of Program

Data processing in Excel-VBA is a sequence of step; read execution print, as described
previously in Section 3.2, which is (rewritten for your you convenience):
1. Read input data from the worksheet
2. To execute data input by VBA
3. Print the results to the worksheet
In large data processing either from the input data or from program execution, a computer
programming (coding) must be planned and managed properly. It is intended in order to make
the program flow becomes simple (not complex), saving coding, and it is easy to solve if an

error is occurred.
o In this section
n, the three of the above steps will be discusse
ed that is ho
ow to
create In
nput Form, Output
O Form and writing
g code in VB
BA module. In
I this topic,, we are going to
create a civil structurral analysis program.
p
3.5.1 Input Outp

put Form
In Excel--VBA progra
amming structure, input--output form
m is worksheet itself and
d can be ma
ade in
the same
e sheet, e.g
g. on Sheet 1. For large
e input data
a, for an exa
ample for sttructural ana
alysis
program as shown in
n Figure 3.4
4, it would be
e better if th
he input data
a extends to
o the right, fo
or an
example, the data entered
e here
e is consistt of numberr of joint, co
oordinates and
a propertie
es of
structure
e elements. By
B creating form widenss on the righ
ht side, the location (rea
ad: rows) off data
output is therefore ca
an be speciffied in the prrogram.
Figurre 3.4: Input-Output Ex

xample Form
m
Figure 3.4 shows that how mucch the input data entere
ed in the forrm does nott change the
e row
numberss for data ou
utput. It thuss obtains a fixed
f input-o
output form and
a does no
ot depend on the
size of th
he input data
a, so that cod
de of read - print steps can
c be specified in the program.
p
Contrary to the inputt data, it is possible

p to exxtend the ou
utput data do
ownward be
ecause end of
o the
ot a boundarry that has to
row is no t be identified for progrram input. Therefore, the
e data prese
ented

gspot.co.id
here following the common way in presenting the results of the structure analysis extending
downward to improve readability.
3.5.2 Work With Modules
A large data processing such as for an application program generally consists of several smaller
sub-programs united into one unit for a goal that is well managed. What is meant by sub-
program is VBA Sub procedure or function procedure given the scope of their respective task
and then placed into a module.
For the convenience, Author distinguishes each procedure into three modules according to their
tasks as follows:
Module 1: Main Module
Module 1 consists of data provision from input-output form (see Figure 3.3), which is mainly
from: geometry data, materials data, determination of node support or fixities and load data.
Module 1 is the main module, which controls the process of read - execution - print.
Module 2: Calculation Process (Analysis)
Procedures on module 2 are associated with the theory of the structure and the analytical
methods used, e.g. building stiffness matrix, generate load matrix and compute nodal
displacement. Module 2 works based on data provided by Module 1.
Module 3: Creating Charts
Procedure in module 3 are associated with the presentation of data in graphical form, e.g. for 1).
Input data: geometry of the structure and loads, and 2). Data output: displacement, moment and
shear distribution. Module 3 works based on data provided by Module 1 and 2.
Figure 3.5 shows in general how Excel-VBA programming structure is depicted in a chart.
Figure 3.6 shows a chart of programming structure in real world; consist of main procedures
(abbreviated to facilitate in writing) that exist in TRUSS2D, program for truss analysis as
presented in Chapter 7.

Figure 3.5:
3 Excel-VB
BA program
mming struc
cture
F
Figure 3.6: Excel-VBA
E p
programmin
ng structure
e for truss analysis
a (TR
RUSS2D)

gspot.co.id
3.5.3 Tips
In design
ning macros, the movem
ment betwee
en workshee
et and Visua
al Basic Edito
or (VBE) win
ndow
is often done
d becaus
se both windows are in
nterconnecte
ed. To open VBE, click on Visual Basic
B
icon Dev
veloper tab and from VBE
V click Miicrosoft Exc
cel icon to go
g back to worksheet,
w or by
pressing shortcut Allt+F11 on both
b windows. Actually, to accelerate the proce
ess of displa
aying
window is to place together in the monito
or screen ass shown in Figure 3.7. Thus, if a VBA
program code is executed, its intteraction to the
t workshe
eet can be viewed simulttaneously.
e the steps to
Here are o put a workksheet and VBE
V in one screen
s and simple
s techn
niques to exe
ecute
line by lin
ne:
1. In
n the worksh
heet window
w click Resto
ore Down iccon in the to
op right corner, then drag the
sh
heet to the left,
l place it approximately half size
e of the mon
nitor screen. Drag the mouse
m
pointer on the
e worksheet window end
ds to adjust its size.
2. N
Next, go to th ndow, in the same way place it on the right sid
he VBE win de of the mo
onitor
sccreen. There
e are now already both these
t window
ws simultaneously.
3. N
Now, go to th
he VBE by clicking
c its window.
w Clickk Close VBA
A Project iccon to close VBA
Project windo
ow on the lefft so only VB
BA module is
i now visiblle (to displayy back, clickk icon
Project Expllorer or press Ctrl+R). To start fro
om the first line press F8
F and presss F8
gain to contiinue to the next
ag n line. If you
y want to jump
j to the desired
d line press Ctrl+F8 at
th
he cursor. Ea
ach line to be executed is marked in
n yellow.
4. Press F9 to select
s progra
am line at which
w executtion will stop
p (breakpoin
nt) when you run
m
macro by clic
cking Run Macro
M (F5). The
T other way
w is to clickk the left ba
ar at program
m line
th
hat you wantt to set a bre
eakpoint untiil red dot app
pears. Click on the red dot
d to removve.
5. The variable value com

mputed durin
ng program
m execution can be se
een its valu
ue by
atttaching a po
ointer on a variable,
v as shown
s in Fig
gure 3.8. This variable re
eads the vallue of

gspot.co.id
the cell of the worksheet. Hold on a second the pointer until the value is displayed, which
is the cell value of the worksheet on the left.
Through step 1 and 2, the movement between worksheet - VBE is simply done by clicking one
window without losing another window. This kind of method is very helpful if you want to run the
program and see the results at once (through step 3 and 4).
The above steps are very helpful in developing flow of program and solve errors as well.
Running step 3 and 4 is actually the process of tracing errors (debugging) by using commands
in the Debug menu and it is also similar to the process of program verification (testing) by
running code line by line. Testing and debugging is a routine activity to build a program that
requires a number of control structures (branching and looping). A condition where a program is
paused for the purpose of developing the code is called break mode.
Figure 3.7: VBA Worksheet window is placed in one screen
One of error type that may occur is when element of an array smaller or larger than the size that
has been assigned to, so when program is running it will show an error message box "subscript
out of range". For example in a looping where the number of process exceeds the subscript of

an element. To get to the location where an error occurs clicks Debug, and Reset to stop the
program.
Figure 3.8: To see a variable value by attaching pointer
3.6 Chart Macro

Excel chart is an object that can be manipulated by VBA through its methods and properties, for
example, in an aim of adding lines (Add series), changing the line color, border color, creating
line styles and setting the thickness (Weight). By creating macros, the job becomes very easy,
fast, and advances to the most difficult level, which would not be done manually.
In Chapter 1 we have discussed how to make lines in the worksheet with formulas and tables. In
this section, it will be created using macro instead of what has done manually. For a case,
Example 3 of Chapter 1.10 will be used here - which is to create truss structure drawing with
joint and coordinate tables as shown below:

Table of joints
Joint x y
1 1.0 0.0
2 1.5 2.5
3 2.0 5.0
4 2.5 10.0
5 3.0 15.0
6 3.5 20.0
7 4.0 15.0
8 4.5 10.0
9 5.0 5.0
10 5.5 2.5
11 6.0 0.0
12 3.5 2.5
Table of Coordinates
1 1 2
2 2 3
3 3 4
4 4 5
5 5 6
6 6 7
7 7 8
8 8 9
9 9 10 Data is programmed
10 10 11
11 2 12
12 12 10
13 3 12
14 12 9
15 3 8
16 3 9
17 4 9
18 4 8
19 4 7
20 5 8
21 5 7

By creating macro, the data of x1, y1, x2, y2 of above table do not need to be formulated
because it will automatically be created by the program. The Macro is shown as follows:
Sub MChart1()
Dim npt, nln, i
npt = 12 'no. of joint
nln = 21 'no. of line
ReDim x(npt), y(npt)
ReDim x1(nln), y1(nln), x2(nln), y2(nln)
'Read joint coord
For i = 1 To npt
x(i) = Cells(4 + i, 2)
y(i) = Cells(4 + i, 3)
Next i
'Read line coordinates
For i = 1 To nln
x1(i) = x(Cells(4 + i, 6))
y1(i) = y(Cells(4 + i, 6))
x2(i) = x(Cells(4 + i, 7))
y2(i) = y(Cells(4 + i, 7))
Next i
ActiveSheet.ChartObjects("Chart 1").Activate
ActiveChart.ChartArea.Select
'creating lines
For i = 1 To nln
With ActiveChart
.SeriesCollection.NewSeries
.SeriesCollection(i).XValues = Array(x1(i), x2(i))
.SeriesCollection(i).Values = Array(y1(i), y2(i))
.SeriesCollection(i).Name = "="""""
End With
ActiveChart.SeriesCollection(i).Select
Next i
End Sub

Now, we
e focus on three
t blockss of stateme
ents enclose
ed by a red
d line, as be
elow. The blocks
descriptio
on are show
wn thereafterr.
Descrription:
B
Block 1 is th
he same to that done using
u VLOO
OKUP formula in Chapte
er 1. This iss the
emptied line coordinates in the shee
et table abovve because it is now pro
ocessed by VBA.
Figure out this, formulas as many ass 21 rows x 4 columns (84
( formulass) for creatin
ng 21
lin
nes be is sim
mply replace
ed with onlyy 6 VBA stattements with
hout limit to any a numb
ber of
lin
nes that will be created, 100 or 200,, 1000, etc.!
B
Block 2 is co
ode to activa
ate a Chart that
t to be manipulated,
m in this case
e is Chart1. Note
th
hat, we will never (everr) create a new
n Chart (C d) because it will make code
Charts.Add
lo
onger. Instea
ad, we will proceed ma ough Insertt > Charts > Scatters)) and
anually (thro

gspot.co.id
format the chart as desired. To identify a chart name right-click on the Chart Area >
Assign Macro, and then see the name in the Macro Name text box.
Block 3 is steps that we do with Add Series in the Select Data Source dialog box for
data entry (see Chapter 1). This VBA looping replaces job for line coordinates data entry
which takes a long time when done manually.
Excel input form for chart macro above as below:
If you want to go further in creating chart it can be continued by setting its properties, such as
line color, marker style and color, or line style (continuous, dashed). The following is code to set
the Chart properties (code enclosed by red lines)
Sub MChart1()
Dim npt, nln, i

npt = 12 'no. of joint
nln = 21 'no. of line
ReDim x(npt), y(npt)
ReDim x1(nln), y1(nln), x2(nln), y2(nln)

'Read joint coordinates
For i = 1 To npt
x(i) = Cells(4 + i, 2)
y(i) = Cells(4 + i, 3)
Next i
'Read line coordinates

For i = 1 To nln
x1(i) = x(Cells(4 + i, 6))
y1(i) = y(Cells(4 + i, 6))
x2(i) = x(Cells(4 + i, 7))
y2(i) = y(Cells(4 + i, 7))
Next i
ActiveSheet.ChartObjects("Chart 1").Activate
ActiveChart.ChartArea.Select
'creating lines
For i = 1 To nln
With ActiveChart
.SeriesCollection.NewSeries
.SeriesCollection(i).XValues = Array(x1(i), x2(i))
.SeriesCollection(i).Values = Array(y1(i), y2(i))
.SeriesCollection(i).Name = "="""""
End With
ActiveChart.SeriesCollection(i).Select
With Selection
.MarkerBackgroundColorIndex = xlNone
.MarkerForegroundColorIndex = xlNone
.MarkerStyle = xlNone
.Smooth = False
End With
With Selection.Border
.ColorIndex = 3 'red
.Weight = xlMedium
.LineStyle = xlContinuous
End With
Next i
End Sub
The above examples is about setting the type and color of lines, which is respectively,
continuous (xlcontinuous) and red (color index = 3). This code uses Border property of the
Selection object, to returns Border object with its properties that are ColorIndex, Weight and
LineStyle. Selection represents SeriesCollection(index) object of ActiveChart object that is an
active chart named Chart 1.

For an exercise, a file contains this interpolation example can be seen in Chart Macro.xls on the
CD included in this book.
3.7 Manipulation on Program Steps

An application program (for specific application) built through a sequence of steps that took
quite long process, start from defining problems, to find solutions, build algorithms, to write and
debug the code. An algorithm is defined as a sequence of logical steps in problem solving.
Algorithm is closely related to the use of the control structures in program that are branching
and looping that have been discussed earlier. In array data basis such as structural analysis, it
is necessary to have little bit of knowledge how to manipulate of program steps or flow of array
data using control structures. For examples, is to read and print data, changing the composition
of elements of an array and create integer variables for storing data.
Example 1
Given a simple example; look at an array data in the worksheet as shown below. Build elements
of array on the left side to be as on the right side, where only the main diagonal value will still
remain while the other = 0.
Here is a code to solve Example 1, starts from reading the data, manipulating the data (code in
red line) and print the results:

Sub ExSTEP1()
Dim i, j, n
ReDim myArray(6, 6) As Integer
'Read input data

For i = 1 To 6
For j = 1 To 6
myArray(i, j) = Cells(2 + i, j)
Next j
Next i
'Manipulate data
For i = 1 To 6
n = myArray(i, i)
For j = 1 To 6
myArray(i, j) = 0
Next j
myArray(i, i) = n
Next i
'Print ouput data

For i = 1 To 6
For j = 1 To 6
Cells(2 + i, 7 + j) = myArray(i, j)
Next j
Next i
End Sub
The code (outlined in red line) is the Authors version. Every programmer must has own steps to
proceed, based on his/her interpretation in the same problem. Similarly, if you are a
programmer or whose accustomed to use Excel-VBA.
Example 2
Create an elements of an array in Example 1 in sequence from bottom to top or from the largest
subscript (myArray (6,6)) to the smallest (myArray (1,1)).

The code
e (The Autho
ors version):
Sub ExST
TEP2()
Dim i, j,
j m, n
ReDim my
yArray(6, 6) As Inte
eger
ReDim rv
vArray(6, 6) As Inte
eger
'Read in
nput data
For i = 1 To 6
For j = 1 To 6
myArray(i
i, j) = Cel
lls(2 + i, j)
Next
t j
Next i
'Manipul
late data
m = 0
n = 0
For i = 6 To 1 St
tep -1
n = n + 1
For j = 6 To 1 Step -1
m = m + 1
rvArray(j
j, i) = myAArray(m, n)
n
Next
t j
m = 0
Next i
'Print data
d
For i = 1 To 6
For j = 1 To 6
Cells(2 + i, 7 + j)
) = rvArra
ay(i, j)
Next
t j
Next i
End Sub

gspot.co.id
Example 3
Assemble the upper triangular array on the left side into the arrangement shown on the right
side, as below.
This example is made as an exercise for readers to make their own code. As previous
examples, the data manipulation in this example is done and completed within two looping. Use
always the Debug menu to verify the program steps. Good luck and hopefully these examples
are helpful to develop your Excel-VBA skill.
Note:
An algorithm or in the (Authors) other words is a manipulation on program steps (flow) is

basic procedure to tell computer how the program will be executed step by step. Although it
may formally a guideline to be followed but everything is gained from experiences and also
exercises. So, the coding more depends on instinct or interpretation of an author on the
problems faced. However, it was still agreed that the fewer steps needed to solve a problem,
the better the program built, meaning the faster execution time by the computer. Thus, it is
not because of the fewer program lines created.
A Note for Beginners:
If you are a beginner in program algorithm, it is fine because your lesson has just started here.
Run the program code of Example 1 and 2 step by step using the Debug menu (for example
by pressing F8) and be focused on the code in the red lines. In break mode, see how the
algorithm for both examples and try to understand the idea behind those code writing. After you
have understood, try to solve Example 3 or any problem you define by yourself, for example that
are developed from Example 1 and 2.

Chapter 4:
MATRIX PROGRAM
The calculation using matrix operation is very suitable and easy to be performed by a computer.
Thus, in line with rapidly growth of computer usage, matrix method is becoming very popular
and widely used instead of the analysis with manual method. It is widely used in mathematical
siences or in civil engineering field.
4.1 Matrix Definition

A matrix is a rectangular array which involves a series of numbers of constituent components or
elements.
Example
a11 a12 . . . . a1n

a
21
a
A = 31

am1 am 2 amn
The notation for a matrix is usually used [] and () or {} for row or column matrix. A matrix can be
written [A] = [aij], where aij is element of the matrix, i = 1,2, .. m and j = 1,2, ... n. Subscript m
indicates the number of rows and n for the number of column, or called has order (m x n). If m =
n, then called a square matrix.

4.1.1 Types of Matrix
1. Row matrix or row vector
where m = 1.
Example: {1 2 3 4}
2. Column matrix or column vector
where n = 1.
1
2

Example: 3

4
3. Square Matrix
Diagonal matrix
All elements = 0, except in the main diagonal.
1 0 0 0
0 2 0 0

0 0 3 0

0 0 0 4
Upper Triangular Matrix
All elements under the main diagonal = 0.
a11 a12 a13 a14 ... a1n

0 a 22 a 23 a 24 ... a 2 n

0 0 a33 a 34 ... a3n

0 0 0 a 44 ... a 4 n
...

0 0 0 0 ... a nn

Lower Triangular Matrix
All elements above the main diagonal = 0.
a11 0 0 0 ... 0
a a 22 0 0 ... 0
21
a31 a32 a33 0 ... 0

a 41 a 42 a 43 a 44 ... 0
... 0

a m1 am2 am3 am4 ... a mn
Scalar matrix
Scalar matrix is a diagonal matrix where the diagonal elements are the same number.
4 0 0 0
0 4 0 0

0 0 4 0

0 0 0 4
Unit matrix
Unit matrix is a scalar matrix in which the value of elements = 1. This matrix is also called
identity matrix and generally denoted by [I]. Multiplying by [I] will return the origin matrix.
1 0 0 0
0 1 0 0

0 0 1 0

0 0 0 1
Band Matrix
Band matrix is a square matrix (n x n) where the non-zero elements are grouped and form a
diagonal elements band.

a11 a12 0 0 ... 0 0
a a 22 0 0 ... 0 0
21
0 0 a33 a34 ... 0 0

0 0 a 43 a 44 ... 0 0
0 0 0 0 ... a n 1,n 1 a n 1,n

0 0 0 0 ... a n ,n 1 a n,n
Symmetric matrix
A matrix [A] is announced symmetric if,
[A]T = [A] or aij = aji
Example:
1 4 5
4 2 6

5 6 3
4.1.2 Matrix Operation
Matrix Addition and Subtraction
If [A] and [B] are two matrices that have equal dimension, then they can be added.
If [A] = [aij] and [B] = [bij] returns [C] = [Cij], it is written:
[C] = [A] + [B]
cij = aij + bij for each i dan j.
Example
2 4 3 0 3 2
[A] = 1 5 4 , [B] = 4 2 3

2 + 0 4 3 3 + 2 2 1 5
[A] + [B] = 1 + 4 5 + 2 4 3 = 5 7 1


Subtract [B] from [A] is equal to add [A] with [-B] or written: [A] - [B] = [A] + [-B].
Matrix Multiplication
Multiplication of two matrices, [A] (m x n) by [B] (n x p) will return [C] of order (m x p). When do
the matrix multiplication the number of columns [A] must be equal to the number of rows [B].
If [A] = [aij], [B] = [bij] and [C] = [Cij] then,
cij = ai1b1j + ai2b2j + + aimbmj
or written as,
m
cij = a
k =1
ik bkj
Example 1
a11 a12
b11 b12
[A] = a 21 a 22 , [B] = b b
a31 a32 21 22
a11 a12
b11 b12
[C] = [A] x [B] = a 21 a 22 b b
a31 a32 21 22
a11b11 + a12 b21 a11b12 + a12 b 22

a 21b12 + a 22 b 22
= a 21b11 + a 22 b 21
a 31b11 + a32 b 21 a31b12 + a32 b 22
Example 2
2 3
3 3
[A] = 4 1 , [B] = 1 5
7 4

2 3
3 3
[C] = [A] x [B] = 4 1 1 5
7 4
2.3 + 3. 1 2.3 + 3.5

= 4.3 + 1. 1 4.3 + 1.5
7.3 + 4. 1 7.3 + 4.5
3 21

= 11 17
17 41
Matrix Inverse
Division is not defined in a matrix operation, for it used the inverse of a matrix. Matrix inversion
-1
is only done on square matrix and written with the notation [ ] . Multiplication by matrix inverse
is a division by the matrix. For example: [A] x [B] = [C] where [A] is a square matrix, the equation
can be written into a division operation: [B] = [A] -1 [C]. [A] -1 is called the inverse of [A].
Matrix Transpose
If [A] is a matrix of order (m x n), then the transpose of [A] is a matrix of order (n x m) or written
as [A]T. Rows and columns of [A] become the columns and rows of [A]T.
[B] = [A]T
bij = aji
Example
1 2 3
[A] =
4 5 6

1 4
5
[A]T = 2
3 6
Some properties associated with matrix transpose:
([A]T)T = [A]
([A]+[B])T = [A]T + [B]T
([A][B])T = [B]T [A]T
System of Linear Equations
Many problems in civil engineering is represented by systems of "n" linear equations with
unknown number "n" and this can be solved using matrix methods. It is often, the number of
equations equal to number of unknowns.
Example: two systems of equations
3x + 2y = 9
5x y = 2
When expressed in matrix form it becomes:
3 2 x 9
5 1 y = 2

In general form, "n" linear equations with "n" unknown number can be expressed as:
a11 a12 a13 ... a1n x1 b1

a a22 a23 ... a2 n x2 b2
21
a31 a32 a33 ... a3n x3 = b3
(4.1)
M M M M M M M
a n1 ... a nn xn bn
an 2 an 3
Simplified into:

[A].{X} = {B} (4.2)
where,
[A] square matrix represents coefficients of Equation 4.2
{X} unknowns column vector
{B} constants column vector
Simply to say that the solution of systems of "n" linear equations in matrix form as expressed in
Equations 4.2 is finding {X}, where [A] and [B] are known values. There are many methods used
out there to solve Equation 4.2, some will briefly be discussed in examples below.
1. Elimination Back Substitution
Example
2x + 3y = 7 (1)
3x 2y = 4 (2)
This is an example of two linear equations with 2 unknowns, x and y.
Elimination
Eliminate x from both equations by subtracting Equation 2 from Equation 1 that has been
given a multiplying factor (Eq. B1):
3 9 21
2x + 3y = 7 (1) x 3x + y= (B1)
2 2 2
3x 2y = 4 (2)
13 13
y=
2 2
y =1
Back Substitution

Substituting y value back into Equation 1, then we have x = 2.
x 2
Thus, solution for the unknowns {X}: =
y 1
2. Directly get the inverse of [A] or [A] -1
From Equation 4.2:
[A]{X} = {B}
Multiply both matrices by [A]-1:
[A]-1[A]{X} = [A]-1{B}, it is known [A]-1[A] = [I]
[I]{X}=[A]-1{B}
{X}=[A]-1{B} (4.3)
Thus, if [A] -1 is known, Equation 4.3 can be solved.
Example
Given two a system of two linear equations:
2x + 3y = 7
3x 2y = 4
In this example the matrix inverse will be found using by Gauss-Jordan method. The
solution is to put the unit matrix [I] next to coefficient matrix [A] such that when the
process is finished [A] become [I] and the previous [I] become matrix inverse [A] -1. The
overall process can be just expressed as:
[A M I ] [ I M [A]-1 ]
The process is done through elementary row operation on [A] and [I] as well to reduce
[A] to [I]. Thus, it is a sequence of row operation to turn elements of [A] = 1 on the main
diagonal and zero elsewhere. For detail see the following example:

Derive matrix coefficient:
2 3
[ A] =
3 2
Put [I] next to [A]:
2 3 1 0
3 2 0 1

Step 1: From above matrix, divide row 1 by 2
3 1 0
1
2 2
3 2 0 1

Step 2: Add row 2 to row 1 times -3
6 1
1 2 2 0
13 3 1
0
2 2
2
Step 3: Multiply row 2 by
13
1
3 0
1 2 6 2
2
0
1 13

26
3
Step 4: Add row 1 to row 2 times
2
8 6
1 0 52 26

0 1 6 2

26 13
So,

8 6
52 26
[A]-1 = 6 2

26 13
Enter [A]-1 into Equation 4.3
{X} = [A]-1.{B}
8 6
x 52 26 7
= 6 2
y 4
26 13
104
x 52
= 26
y
26
So, the solution to this system,
x 2
=
y 1
3. Iterative Method
This method is actually an approach and done by taking initial value of the unknown then
corrected and refined to enter the next value. We have examined the iterative method on
the example above in the previous discussion in Chapter 3. Here is now an example
from a system of three linear equations:
5x + 4y +3z = 12 (4.4)
4x + 7y + 4z = 15 (4.5)
3x + 4y + 4z = 11 (4.6)
Because the iterative process is quite long, so in this example we will be working with
help of a computer through Excel iteration.

Re-write Equation 4.4, 4.5 and 4.6 into forms below:
Equation 4.4: x = 12/5 4y/5 3z/5 = (12 4y 3z)/5
Equation 4.5: y = 15/7 4x/7 4z/7 = (15 4x 4z)/7
Equation 4.6: z = 11/4 3x/4 4y/4 =(11 3x 4y)/4
Next, convert the equations into worksheet formulas and enter in cell B2, C3 and D4.
However, until this stage no iteration to be done because the formulas are still standing
on their own equations.
The iterative process will be carried out after these connections:
The value of y, z of Equation 4.4 = value y, z of Equation 4.5
The value of x, z of Equation 4.5 = value x of Equation 4.4 and z of Equation 4.6
The value of x, y of Equation 4.6 = value x, y of Equation 4.5
In the worksheet, the equations are done as follows:
The iteration steps shown in the worksheet below:

1 x iteration:
5 x iterations:
10 x iterations:
20 x iterations:
50 x iterations:

The iteration process of finding {X} of Equation 4.4 is actually to make x, y, z to the same value
on the worksheet columns B, C and D in satisfying the given equations. It appears that they are
very close in 50th iteration and by rounding obtained:
x 1

y = 1
z 1

4.2 Program for Matrix Operations

Multiplication and inversion of matrix are two main operations in matrix equation found in most
of engineering calculations using matrix method. In this discussion, it will demonstrate how to
create program code for multiplication and inversion of matrix.
Matrix Multiplication
In the following example is given a program to calculate the multiplication of two matrices. The
program is named MMULT1, input data and the results will be done on Sheet1.
Example
Find the multiplication of two matrices in the following equation:
[C] = [A], [B]
where,
2 3
3 3
[ A] = 4 1 , [B ] =
7 4 1 5
We know that the multiplication of a matrix is the multiplication of element row in one matrix by
element column of other matrix. If done manually as in the example of Section, 4.1.2, code to
get [C] is shown as the following:

Sub MMULT1()
Dim MA(3, 2) As Double

Dim MB(2, 2) As Double
Dim MC(1 To 3, 1 To 2) As Double
'===============
'Read Input Data
'===============
'[A] or MA
MA(1, 1) = Cells(6, 2)
MA(2, 1) = Cells(7, 2)
MA(3, 1) = Cells(8, 2)
MA(1, 2) = Cells(6, 3)
MA(2, 2) = Cells(7, 3)
MA(3, 2) = Cells(8, 3)
'[B] or MB
MB(1, 1) = Cells(6, 5)
MB(2, 1) = Cells(7, 5)
MB(1, 2) = Cells(6, 6)
MB(2, 2) = Cells(7, 6)
'=============
'Multiplying
'[C] = [A].[B]
'=============
MC(1, 1) = MA(1, 1) * MB(1, 1) + MA(1, 2) * MB(2, 1)
MC(2, 1) = MA(2, 1) * MB(1, 1) + MA(2, 2) * MB(2, 1)
MC(3, 1) = MA(3, 1) * MB(1, 1) + MA(3, 2) * MB(2, 1)
MC(1, 2) = MA(1, 1) * MB(1, 2) + MA(1, 2) * MB(2, 2)
MC(2, 2) = MA(2, 1) * MB(1, 2) + MA(2, 2) * MB(2, 2)
MC(3, 2) = MA(3, 1) * MB(1, 2) + MA(3, 2) * MB(2, 2)
'==========
'Print [C]:
'==========
Range(Cells(12, 2), Cells(14, 3)).FormulaArray = MC
End Sub
Worksheet form that fits to the above program code is as below:

In this example, [A] has order 3 x 2 and [B] has order 2 x 2 to return [C] of order 3 x 2. Input
data of elemen [A] is made on a range of cells (6,2) to cell (8,3), while element [B] on a cell
range (6.5) to cell (7.6). To run the program, click the command button that has been made in
the worksheet. The results are printed on a range of cells (12, 2) to cell (14, 3), using statement
of FormulaArray = [C], that shown as below:
For a large-sized matrix, the given MMULT1 example is inefficient in term of wasting time and
coding to write code for all matrix (row) operations one by one. Instead, use For-Next looping
for reading input data, to process and print the results. This way saves a lot of code writing.
The code program of above matrix multiplication can be modified by giving the row and column
indices (i, j and k) to [A], [B] and [C], which is done in three looping below:

In a looping above, m and n represents number of rows and columns of [A], respectively, while
n and p represents number of rows and columns of [B]. In the following example, we will use the
modified program code to multiply a 6x6 matrix by a 6x3 matrix. The program is named
MMULT2.
Sub MMULT2()
Dim m, n, p As Integer
m = Cells(2, 5)
n = Cells(3, 5)
p = Cells(4, 5)
ReDim MA(1 To m, 1 To n) As Double
ReDim MB(1 To n, 1 To p) As Double
ReDim MC(1 To m, 1 To p) As Double
'Input data [A]

For i = 1 To m
For j = 1 To n
MA(i, j) = Cells(5 + i, 1 + j)
Next j
Next i
'Input data [B]

For i = 1 To n
For j = 1 To p
MB(i, j) = Cells(5 + i, 9 + j)

Next j
Next i
'Multiplying
'[A](mxn) x [B](nxp) returns
'[C](mxp)
For i = 1 To m
For j = 1 To p
MC(i, j) = 0
For k = 1 To n
MC(i, j) = MC(i, j) + MA(i, k) * MB(k, j)
Next k
Next j
Next i
'Print result
For i = 1 To m
For j = 1 To p
Cells(15 + i, 1 + j) = MC(i, j)
Next j
Next i
End Sub
In this program, the printing of the results that used Array Formula has been replaced by two
orders of For-Next looping. The first and second looping is respectively for printing rows and
columns of [C]. This way is much better in term of automation, because we do not need to
calculate the range of the array that will be generated as it has been known as m x p.
Matrix multiplication can also be done in a worksheet using Excel built-in function MMULT
(array1, array2). For MMULT2 procedure example, a formula can be written as follows:
= MMULT(B6:G11,J6:L11).
To print the result use INDEX function or you can enter the formula above as an array formula.
See about array formula in Chapter 1.5.
For further details, the file contains this example, worksheet formulas and the associated
program code can be seen on CD included in this book.

A B C D E F G H I J K L
1
2 [A] = mxn m= 6
3 [B] = nxp n= 6
4 [C] = mxp p= 3
5
6 [A] = 20 5 0 0 0 0 [B] = 0.25 0 0
7 5 20 0 0 0 0 0.25 1 0
8 0 0 30 10 0 0 -0.50 1 0
9 0 0 10 30 0 0 -0.50 0 1
10 0 0 0 0 12 6 0.25 0 1
11 0 0 0 0 6 12 0.25 0 0
12
13
14 [C] =[A].[B]
15
16 [C] = 6.25 5 0
17 6.25 20 0
18 -20 30 10
19 -20 10 30
20 4.5 0 12
21 4.5 0 6
22
Matrix Inversion
In finding an inverse of matrix, Gauss-Jordan method is among the most popular method to be
used and programmed, with a process that has been described in the previous section. The
steps are an effort to turn a matrix to be searched its inverse (e.g. [A]) becomes an identity
matrix [I] and at the same time, [I] to [A] -1. That is done by elementary row operation as ever
shown in the previous example.

In the program, this can be done by taking one matrix e.g. [C], such that at the end of the
process [C] turns into [A] -1. In the beginning is to divide the first row by element a11 of [A], in an
effort to turn a11 = 1 so that c11 = 1 / a11 and this is analogous to that performed in [I]. If [A] is
n x n matrix, then in the first step the first column of [C] is analogous of that performed in [I],
while the second column to column n is analogous of that performed in [A]. At the end of step 1,
the value of [C] of order 3x3 can be written as follows:
1 a12 a13

a11 a11 a11
[C ] = 21 . a 21
a a a13
a 22 + 12 . a 21 a 23 +
a a11 a11
11
a31 a
a32 + 12 . a31 a33 +
a13
a31
a11 a11 a11
In the second step, the first and second column of [C] is analogous of that performed in [I],
while the third column to column n is analogous of that performed in [A], etc.
Program for finding inverse is created in VBA and the file can be seen on the CD included in this
book. It notes that matrix inversion requires square size (say, m x m) matrix and the determinant
of matrix must not be zero. For the purpose to create worksheet form, the maximum element of
matrix is limited to m = 6, but it can certainly be made more than 6.
Example 1
Find the inverse of [A] = [5].
Solution for Example 1:

Example 2
Find the invers of,
2 3
[ A] =
3 2
Example 3
Find the invers of,
1 3 3
[ A] = 1 4 3
1 3 4
To find inverse of matrix in worksheet can also be used the built-in function MINVERSE. The
way to use it is the same as MMULT.

4.3 Matrix Method for Structural Analysis
Structure can be analyzed as a series of structural elements which are connected to each other
at endpoint called node. Solution of this system of arrangement of elements can be expressed
by simultaneous linear equations that are made in a matrix form. In matrix form, then the
calculation of the structure becomes easier to be solved with the help of a computer.
4.3.1 Upper Structure
The method that will be used for structural analysis is direct stiffness method. A brief
introduction to the basics of this method will be given here; readers are encouraged to read
other references to supplement the discussion on this topic. In short, the aim using of this
method of analysis is to obtain force - displacement relationship expressed as:
{P} = [K]{X} (4.7)
where, {P} vector of external force acting on node
[K] structure stiffness matrix
{X} vector of nodal displacements
By deriving [K], thus it can be obtained force - displacement relationship as expressed by

Equation 4.7. One of the processes here is to derive stiffness matrix of an element, then applied
for other elements and put together to form stiffness matrix of the whole structure. To meet
deformation continuity and compatibility between the elements is by superposing corresponding
force - displacement components.
Furthermore, to simplify the calculation we can partition element of the stiffness matrix at the
nodes that are free to displace (free nodes) and fixed (zero displacement at the supports). If Xb
represents displacements vector at the support, Xf represents displacements vector at free
nodes, Pf and Pb are the external forces vector which correspond to their displacements,
Equation 4.7 can be written as follows:

Pf K ff K fb X f
= (4.8)
Pb K bf K bb X b
If node displacements vector {Xb} = 0, then Equation 4.8 can be written:
{Pf} = [Kff]{Xf} (4.9)
And
{Pb} = [Kbf]{Xf} (4.10)
The displacements vector at free nodes is then obtained using Equation 4.9 that becomes {Xf} =
[Kff]-1{Pf}. If {Xf} is known then the internal forces of the element can be obtained. In this relation,
the internal force of element i is expressed in local coordinate and obtained by the following
equation:
{P}i = {Po} + [K]i{X}i (4.11)
where,
{Po} fixed-end forces: reactions at end of element i due to external load directly acting
on element i
[K]i stiffness matrix of element i in local coordinate system
{X}i deformations matrix = [T]i{Xf}
[T]i transformation matrix of element i
To find the support reactions, it can be directly obtained in the global coordinate of the structure
using Equation 4.10 becomes:
{Rb} + {Pb} = [Kbf]{Xf}
or,
{Rb}= [Kbf]{Xf} {Pb} (4.12)

4.3.2 Sub Structure
The structure analysis in this subject is to obtain a force - displacement relationship of soil -
structure interaction which is modeled by the means of beam on elastic foundation. Generally
it is solved with the finite element method where the system is built by a series of finite elements
connected at endpoint called node.
Some equations are made in order to find force - displacement relationship at the nodes, and
begins by writing the relationship between external forces at the nodes {P} and internal forces of
element {F} in principle of equilibrium.
The relationship between {P} and {F} is expressed as:
{P} = [A]{F} (4.13)
The relationship between node displacements {X} and deformations of element {d}:
{d} = [B]{X}, and from reciprocal theorem it can be shown
{d} = [A]T{X} (4.14)
From elastic behavior of element, the relation between {F} and {d}:
{F} = [S]{d} (4.15)
Subsituting Equation 4.14 into 4.15 to get,
{F} = [S]{d} = [S][A]T{X} (4.16)
Subsituting Equation 4.16 ke 4.13 to get.
{P} = [A]{F} = [A][S][A]T{X} (4.17)
{X} is unknown displacements to be determined, thus
{X} =([A][S][A]T )-1{P} (4.18)
Subsituting {X} back to Equation 4.16 to find internal forces {F}.
Equations 4.16, 4.17 and 4.18 are the fundamental equations in the analysis using finite
element method (J.E. Bowles, 1982).

Chapter 5:
NUMERICAL METHOD
Many problems in science and technology involve the use of integral and differential equations.
If the mathematical equation model is quite simple, then it can be solved by means of analytical
method. However, many of those equations are difficult or may be impossible to be solved
analytically. Instead, it must be solved numerically or using numerical methods.
The result of the numerical method is an approach to its exact value by means of analytically.
Thus, it is just an approximate calculation with estimated accuracy. Numerical methods will
provide many benefits if done by a computer.
5.1 Numerical Integration

The integral of a function f(x) from a to b is written:
f ( x)dx
b
y=
a
has a broad sense as the area bounded by the curve y = f(x) and the X-axis, the between x = a
and x = b.
Analytical solution to integral or exact solution is:
f ( x)dx = [ F ( x) ] = F (a) F (b)

b
y= (5.1)
a
a
where, F(x) is an anti-derivative of f(x) such that F'(x) = f(x).
Numerical integration to calculate area under the curve uses the concept of approach whereby
the area is divided into smaller pieces of area in such that, if combined the total area of the
pieces approaching the exact result. In this example we use trapezoidal method where the area

under the curve is divided into small pieces of area with sides f(xn) and f(xn-1), that resembles a
trapezoid as shown in Figure 5.1.
Figure 5.1: Trapezoidal method for approximating area
Trapezoids are having the same width for x, if there are n intervals (or n trapezoids) from a to
b, then the x = (b - a) / n as shown Figure 5.1. Figure 5.2 shows the comparison between the
area obtained from the total area of the trapezoids and the exact results obtained using
Equation 5.1. The formula to find total area (L) is given by:
x
Ltotal = (( f ( xo ) + f ( x1 )) + ( f ( x1 ) + f ( x2 )) + ... + ( f ( xn 1 ) + f ( xn )) (5.2)
2
x
= ( f ( xo ) + 2 f ( x1 ) + 2 f ( x2 ) + 2 f ( x3 )) + ... + 2 f ( xn 1 ) + f ( xn )) (5.3)
2

Figure 5.2: The comparison of area resulted from trapezoids (left) and
the exact one (right)
Figure 5.1 and 5.2 show, the calculation by means of trapezoidal method is getting closer to the
exact result when taking more trapezoids into account, or by making smaller x.
Furthermore, to complement this discussion we will examine the given example below via VBA.
Example
Find the area bounded by the parabola y = 4x - x2, the X-axis, from x = a and x = b. Find also
the exact value.
Solution
VBA code
Function Trapez(a, b, n)
a = Cells(3, 2)
b = Cells(4, 2)
n = Cells(5, 2)
Dx = (b - a) / n
x = a
For i = 2 To n - 1
x = x + Dx
Sum = Sum + 2 * f(x)
Next i
Sum = Sum + f(a) + f(b)

Trapez = Sum * Dx / 2
End Function

Function f(x)
f = 4 * x - x ^ 2
End Function
Function integf(x)
integf = 2 * x ^ 2 - 1 / 3 * x ^ 3
End Function
In this experiment, we take a = 0 and b = 4 respectively, for n = 10 and 100. The values of a, b
and n are then used as function arguments of Integf and Trapez in the above program. The use
of function in the worksheet and the results are shown below:
In the worksheet, the area calculated using trapezoidal method and its exact value is
respectively expressed as A_num and A_exact. This experiment shows that the smaller width of
trapezoids taken with n = 100, the closer A_num to A_exact, compared with n = 10. However,

the more thorough results requires a longer calculation process and it is certainly suitable done
with the help of computer.
The other methods for approaching area, for examples are by the Simpson method and Gauss
quadrature.
5.2 Numerical Differentiation

Differential equation is any equation that contains derivatives of function.
Example:
dv d 2x
F = m.a = m =m 2 (5.4)
dt dt
d 2v M
= (5.5)
dx 2 EI
u 2u
= cV (5.6)
t z 2
Description:
Equation 5.4 is Newton's second law of motion that states force (F) acting on an object =
mass (m) of the object multiplied by its acceleration (a) where, acceleration (a) is the first
derivative of velocity (v) with respect to time (t) and the second derivative of position (x).
Equation 5.5 shows the relationships between bending moment (M) at any point on the
X-axis and transverse deflection (v) of the beam. EI is flexural rigidity of the beam.
Equation 5.6 is the partial differential equation of one-dimensional consolidation, where
u, t, z are respectively, pore-water pressure, time and the depth of clay layer. Cv is a
coefficient of consolidation where the value is assumed to be constant.
Numerical differentiation solution is an approach to continuous differential form into discrete or

finite model. In this discussion we adopt Taylor approach for it is among the most popular and
widely used in introducing numerical approximation. Taylor theorem for approaching a function
is graphically shown in Figure 5.3.

Figure 5.3: Taylor theorem for approximating a function
The value of the function at point xi+1 located at a distance x from point xi can be approximated
by Taylor series below:
x x 2 x 3 x n
f ( xi +1 ) = f ( xi ) + f ' ( xi ) + f '' ( xi ) + f ''' ( xi ) + ... f n ( x i ) + Rn (5.7)
1! 2! 3! n!
where,
- All derivatives of function through the point (xi, f(xi)).
- Rn: truncation error when the series counted up to n-order.

Taylor series will give an estimation of the function correctly if all derivatives in the series are
taken into account, where n = infinite number. However, in practice, the use of Taylor series is
done by simply taking the first few terms of the series.
Equation 5.7 of the Taylor series can be derived into the discrete form (finite) as shown in
Figure 5.4. The estimation of first derivative of a function f'(x) is approached through several
slope of the line through the value f(x) at x = xi-1 (point A), x = xi (point B) and x = xi+1 (point C)
as in Figure 5.4. From Equation 5.7, with respect to first derivative:
f ( xi +1 ) = f ( x1 ) + f ' ( xi )x + Rn
or,
f f ( xi +1 ) f ( xi )
= f ' ( xi ) = Rn (5.8)
x x
Equation 5.8 is called forward difference formula.
Figure 5.4: Estimation of first derivative of a function
Moreover, backward difference formula is given by:

f f ( xi ) f ( xi 1 )
= f ' ( xi ) = + Rn (5.9)
x x
Central difference formula:
f f ( xi +1 ) f ( xi 1 )
= f ' ( xi ) = + Rn (5.10)
x 2 x
If the function contains two independent variables such as f(x, y), then the derivatives in the
Taylor series will have the following form:
f x f y 2 f x 2 2 f y 2 n f y n
f ( xi +1 , y j +1 ) = f ( x1 , yi ) + + + 2 + 2 + ... n (5.11)
x 1! y 1! x 2! y 2! y n!
Equation 5.11 can be written into the following:
Forward difference formula of first derivative:
f f f i, j
i +1 , j
x x
f f i , j +1 fi , j

y y
Central difference formula of first derivative:
f f f i 1, j
i +1 , j
x 2x
f f f i , j 1
i , j +1
y 2y
Central difference formula of second derivative:
2 f f 2 f i , j + f i +1 , j
i 1, j
x 2
x 2
2 f f 2 f i , j + f i , j +1
i , j 1
y 2
y 2

Subscript i, j on the function f(i, j) is a simplification of f(xi, yj) form. A grid points of function
values in the xy coordinate system is shown in Figure 5.5, it is used to estimate derivative of
partial differential with respect to x and y.
Figure 5.5: Two dimensional grid points of a function
Taylors discretization method whereby derivative of a function is approached as the difference

value of a function between the given independent variable (e.g. x1) and a small increment (e.g.
x1+x), is called the Finite Difference Method (Gilberto E.U, 2004). An approach using finite
difference method implies that the smaller x is used the better the result closer to the exact
value as shown from an approach of slope of the line in Figure 5.4.
In Chapter 10, we will build a program using finite difference method for the one-dimensional
consolidation based on Equation 5.6. In the process, it is necessary to put a multiplier into input
data to divide interval of the independent variables (e.g. x) to be smaller so that the solution is
expected close to exact result. Program is created in VBA and of course, with the advantages of
making chart automation to present the result.

Chapter 6:
PROGRAM FOR 2D FRAME

STRUCTURE ANALYSIS
The matrix program as described in Chapter 4 is now implemented for analysis of structure in
this chapter and the next few chapters. It is expected that these series of references for Excel-
VBA programming will be complementary, as basic knowledge for readers to make an
application program.
Please differentiate between the purpose of programming in Excel and the aim to build an
executable program (using Visual Basic or C++ programming language). Using Excel means
the program as a whole to be convenient and practical for it is created case by case. Each case
is packed in such a way to utilize the facilities on the spreadsheet and then presented as in
creating a technical report. This is consistent with the nature of working with Excel.
6.1 Case Example

Structure of plane frame is composed of 3 members with 4 m long each. The stiffness of
member BC is 2 times of the stiffness of members AB or CD. Point B, and C are the free nodes,
each has 3 displacement components. 1 ton
q = 0.6 ton/m'
B C
0.5 ton
2 EI
4.0 m EI EI 4.0 m
A D
4.0 m

All members of the frame are made of the same linear elastic material with modulus of elasticity,
E = 2.1 x 106 ton/m2. Dimension of member AB = BC = CD = 30 cm.
The given load as the following:
Known external load:
H at node B = 0.5 ton
V at node C = 1 ton
Uniform load, q at member BC = 0.6 tom/m
Uniform load, q is converted to equivalent nodal loads at B and C
(see sign convention in Fig 6.1)
PBC(M) = - PCB(M) = 1/12.q.L2 = -0.8 ton.m.
PB(V) = PC(V) = - .q.L = -1.2 ton.
Sign convention represents the direction of the forces - displacements with reference to either
the global coordinate or the local coordinate. It is set as follows: axial and shear forces follow
the right-hand rule and the moments follow the right-hand screw rule, or in the easy way it is
evident from Figure 6.1, i.e. positive when in the direction of the positive axis.
5 8 Positive direction
B C for forces and
4 7
displacements
6 9
2 11 Y
A 1 D 10 X
Z
3 12
Figure 6.1: Force and displacement components at nodes in global

coordinate system

2 5
1 4
B C
3 6
4
5 2
B C
6 3 1
1
2 A D 5
3 6 4
Figure 6.2: Members force and displacement components in local

coordinate system
5 8
4 7
B C
6 9
5 8
4 7
B C
6 9
2 11
A 1 D 10
3 12
Figure 6.3: Members force and displacement components that conform

to global coordinate system

The process of frame analysis consists of the following steps:
1. Stiffness matrix for plane frame member is derived as the following:
EA EA
L 0 0 0 0
L
12 EI 6 EI 12 EI 6 EI
0 0
L3 L2 L3 L2
6 EI 4 EI 6 EI 2 EI
0 0 2
L2 L L L
EA EA
[K]i= 0 0 0 0 (6.1)
L L
12 EI 6 EI 12 EI 6 EI
0 0 2
L3 L2 L3 L
0 6 EI 2 EI 6 EI 4 EI
0 2
L2 L L L
2. Transformation matrix to convert vector components from local coordinate to global

coordinate is as the following:
For 6 force - displacement vector components of each member,
BC = 0
B C
CD = 270
B C
AB = 90
A D

cos sin 0 0 0 0
sin cos 0 0 0 0

0 0 1 0 0 0

[T]i= 0 0 0 cos sin 0 (6.2)
0 0 0 sin cos 0

0 0 0 0 0 1
3. Member stiffness matrix in the global coordinate system:
[K]s = [T]iT[K]i[T]i (6.3)
Member stiffness matrix is then superposed with other member, which altogether form
global stiffness matrix or called [K] structure. This process is written:
Superpose [K]s [K] structure (6.4)
4. Displacements at the free nodes to be determined by Equation 4.9 that becomes:
{Xf} = [Kff] -1{Pf} (6.5)
5. Internal forces on each member can be found by firstly transforming displacements from
global coordinate back to local coordinate,
{X}i = [T]i{Xf}
After the displacements of member {X}i are obtained, then internal forces on each
member can be found using Equation 4:11:
{P}i = {Po} + [K]i{X}I (6.6)
6. Support reactions are determined using Equation 4.12:
{Rb} = [Kbf]{Xf} {Pb} (6.7)
From above sequence of process then we can create a simple program for frame structure
analysis named FRAME2D.

Programming:
In designing matrix program, it is convenient to mastering design-build in matrix-based data. In

this relation, is how to make matrix variables of Equation 6.3 to 6.7. There are a few things in
coding that actually require a little imagination of a programmer.
For instance is how to assembly [K] structure. For that purpose, it requires to manipulate
program steps, so that [K] members can be placed on their position in [K] structure. The idea is
simply creating an integer variable, namely dindex, to equalize subscripts of [K] member.
Dindex reads index (numbering) of members joint displacements in the global coordinate
system and uses the subscripts to store [K] members in [K] structure and do the superposition
at once.
Program code is as follows:
'Member end-displacements index:

For i = 1 To NM
With Member(i)
dindex(1, i) = 3 * .J1 - 2
dindex(2, i) = 3 * .J1 - 1
dindex(3, i) = 3 * .J1
dindex(4, i) = 3 * .J2 - 2
dindex(5, i) = 3 * .J2 - 1
dindex(6, i) = 3 * .J2
End With
Next i
'Storing members and superposition

For i = 1 to NM 'NM= no. of member
For j = 1 To 6
For n = 1 To 6
S(Idm(j, i), Idm(n, i)) = S(Idm(j, i), Idm(n, i)) + M(j, n, i)
Next n
Next j
Next i
Notes:
Idm variable is the same displacement index as stated by dindex variable, but declared
at a different level. Idm is at public-level, while dindex is at procedure-level. For more
clearly understand how these variables proceeded you can see in FRAME2D code,
provided within this book.
The above code is for completion of the description (6.4).

The other thing is to partition (separation) of the structure stiffness matrix to form [Kff] and [Kbf]
as written in Equation 4.8. To do this, it uses two integer variables namely Ifr to form [Kff] and Ifx
to form [Kbf].
The sequence of steps to partition the structure stiffness matrix is as follows:
1. Creating displacement index at free nodes and supports

'Joint free displacement index, IFr and restraint index, IFx:
n = 1
j = 1
For i = 1 To NP
If Rs(i) = 0 Then IFr(n) = i: n = n + 1
Next i
n = 1
For i = 1 To NS
n = Cells(22, 1 + i)
IFx(3 * i - 2) = 3 * n - 2
IFx(3 * i - 1) = 3 * n - 1
IFx(3 * i) = 3 * n
Next i
2. Sub-matrix [Kff] is obtained as follows:
'submatrix Stiffness, KFF

For i = 1 To DOF
For j = 1 To DOF
SD(i, j) = GS(Idj(i), Idj(j))
Next j
Next i
3 Sub-matrix [Kbf] is obtained as follows:
'submatrix Stiffness, KBF

For i = 1 To DOF
For j = 1 To 3 * NS
SR(j, i) = GS(Irj(j), Idj(i))
Next j
Next i
Notes:
Rs represents joint displacement condition at supports, 1 = fixed, 0 = free.
NS = no. of supports
Variable Idj = Ifr, Irj = Ifx, but they are declared at different level which is at public and
procedure level, respectively. See FRAME2D code provided within this book.

The input and output data in form of FRAME2D program are presented in Figure 6.4. Below are
the results of each step of the analysis.
Member AB
[K] (member stiffness matrix)
47250.00 0.00 0.00 -47250.00 0.00 0.00

0.00 265.78 531.56 0.00 -265.78 531.56
0.00 531.56 1417.50 0.00 -531.56 708.75
-47250.00 0.00 0.00 47250.00 0.00 0.00
0.00 -265.78 -531.56 0.00 265.78 -531.56
0.00 531.56 708.75 0.00 -531.56 1417.50
[T] (member transformation matrix from local to global coordinate system)
0.00 1.00 0.00 0.00 0.00 0.00

-1.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 1.00 0.00
0.00 0.00 0.00 -1.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 1.00
[K]s = [T]iT[K]i[T]i (member stiffness matrix in global coordinate)
265.78 0.00 -531.56 -265.78 0.00 -531.56

0.00 47250.00 0.00 0.00 -47250.00 0.00
-531.56 0.00 1417.50 531.56 0.00 708.75
-265.78 0.00 531.56 265.78 0.00 531.56
0.00 -47250.00 0.00 0.00 47250.00 0.00
-531.56 0.00 708.75 531.56 0.00 1417.50

Member BC
[K]
59850.00 0.00 0.00 -59850.00 0.00 0.00

0.00 540.15 1080.29 0.00 -540.15 1080.29
0.00 1080.29 2880.78 0.00 -1080.29 1440.39
-59850.00 0.00 0.00 59850.00 0.00 0.00
0.00 -540.15 -1080.29 0.00 540.15 -1080.29
0.00 1080.29 1440.39 0.00 -1080.29 2880.78
[T]
1.00 0.00 0.00 0.00 0.00 0.00

0.00 1.00 0.00 0.00 0.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00
0.00 0.00 0.00 1.00 0.00 0.00
0.00 0.00 0.00 0.00 1.00 0.00
0.00 0.00 0.00 0.00 0.00 1.00
[K]s
59850.00 0.00 0.00 -59850.00 0.00 0.00

0.00 540.15 1080.29 0.00 -540.15 1080.29
0.00 1080.29 2880.78 0.00 -1080.29 1440.39
-59850.00 0.00 0.00 59850.00 0.00 0.00
0.00 -540.15 -1080.29 0.00 540.15 -1080.29
0.00 1080.29 1440.39 0.00 -1080.29 2880.78

Member CD
[K]
47250.00 0.00 0.00 -47250.00 0.00 0.00

0.00 265.78 531.56 0.00 -265.78 531.56
0.00 531.56 1417.50 0.00 -531.56 708.75
-47250.00 0.00 0.00 47250.00 0.00 0.00
0.00 -265.78 -531.56 0.00 265.78 -531.56
0.00 531.56 708.75 0.00 -531.56 1417.50
[T]
0.00 -1.00 0.00 0.00 0.00 0.00

1.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 -1.00 0.00
0.00 0.00 0.00 1.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 1.00
[K]s
265.78 0.00 531.56 -265.78 0.00 531.56

0.00 47250.00 0.00 0.00 -47250.00 0.00
531.56 0.00 1417.50 -531.56 0.00 708.75
-265.78 0.00 -531.56 265.78 0.00 -531.56
0.00 -47250.00 0.00 0.00 47250.00 0.00
531.56 0.00 708.75 -531.56 0.00 1417.50

[K] Structure
[K]AB [K]BC [K]CD
1 2 3 4 5 6 7 8 9 10 11 12
1 265.78 0.00 -531.56 -265.78 0.00 -531.56 0.00 0.00 0.00 0.00 0.00 0.00
-
2 0.00 47250.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
47250.00
3 -531.56 0.00 1417.50 531.56 0.00 708.75 0.00 0.00 0.00 0.00 0.00 0.00
-
4 -265.78 0.00 531.56 60115.78 0.00 531.56 0.00 0.00 0.00 0.00 0.00
59850.00
-
5 0.00 0.00 0.00 47790.15 1080.29 0.00 -540.15 1080.29 0.00 0.00 0.00
47250.00
An Introduction to: Excel for Civil Engineers

6 -531.56 0.00 708.75 531.56 1080.29 4298.28 0.00 -1080.29 1440.39 0.00 0.00 0.00
-
7 0.00 0.00 0.00 0.00 0.00 60115.78 0.00 531.56 -265.78 0.00 531.56
59850.00
-
8 0.00 0.00 0.00 0.00 -540.15 -1080.29 0.00 47790.15 -1080.29 0.00 0.00
47250.00
9 0.00 0.00 0.00 0.00 1080.29 1440.39 531.56 -1080.29 4298.28 -531.56 0.00 708.75
10 0.00 0.00 0.00 0.00 0.00 0.00 -265.78 0.00 -531.56 265.78 0.00 -531.56
-
11 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 47250.00 0.00
47250.00
12 0.00 0.00 0.00 0.00 0.00 0.00 531.56 0.00 708.75 -531.56 0.00 1417.50
167
Partition at free nodes [Kff]
4 5 6 7 8 9
4 60115.78 0.00 531.56 -59850.00 0.00 0.00
5 0.00 47790.15 1080.29 0.00 -540.15 1080.29
6 531.56 1080.29 4298.28 0.00 -1080.29 1440.39
7 -59850.00 0.00 0.00 60115.78 0.00 531.56
8 0.00 -540.15 -1080.29 0.00 47790.15 -1080.29
9 0.00 1080.29 1440.39 531.56 -1080.29 4298.28
Partition at support nodes [Kbf]
4 5 6 7 8 9
1 -265.78 0.00 -531.56 0.00 0.00 0.00
2 0.00 -47250.00 0.00 0.00 0.00 0.00
3 531.56 0.00 708.75 0.00 0.00 0.00
10 0.00 0.00 0.00 -265.78 0.00 -531.56
11 0.00 0.00 0.00 0.00 -47250.00 0.00
12 0.00 0.00 0.00 531.56 0.00 708.75
0.00000
0.00000

0.00000

0.00117
0.00002

0.00039

0.00117
Joint displacements at free nodes : {Xf} =
0.00005

0.00017
0.00000

0.00000
0.00000


After getting {Xf} then member deformations are calculated by using relationship:
{X}i = [T]i{Xf}
Member internal forces are calculated by using Equation 6.6:
{P}i = {Po} + [K]i{X}i
Support reactions are calculated by using Equation 6.7:
{Rb} = [Kbf]{Xf} {Pb}
The result of the analysis is presented in FRAME2D input-output form as shown in Figure 6.4
and 6.5 below:
Figure 6.4: FRAME2D input-output form

Figure 6.5: Charts of displacement, diagrams of bending moment,
shear force and axial force
Notes for FRAME2D input and output form:
Data input is blue text.

Scale for presenting output data is divided into 2 cases: true scale for beam members
(continuous beam) and auto scale for frame structure.
NM: number of structural members
NJ: number of joints or nodes
FRAME2D distinguishes 3 types of loading: selftweights by inputting material density
(row 16), nodal loads (row 27 to 30), and uniform loads (row 32 to. 33).
When Plot Geometry button is clicked it shows Chart Window of the geometry of
structure to give geometry verification to that inputted before the program is executed by
clicking RUN Program button. It is expected to minimize kind of errors that caused by
numbering system of the members and joints.

Chart Window command used in Plot Geometry is only compatible with Excel 2003 and
earlier, which is a version where FRAME2D was created (also the next few programs).
When running it in Excel 2007 version, Chart Window (as shown in Figure 6.4) will not be
displayed. However, readers can see it placed before displacement chart in the input-
output form.
Figure 6.5 shows chart of displacement as well as diagrams of moment, shear and axial
forces that are drawn automatically within the program. It is also included maximum
values occurred on nodes. Sign convention used for presenting diagrams of the moments
and internal forces is described in Chapter 6.2.
The above calculation does not show the optimum result as the displacement of the
second member looks like a straight line, contrary to the fact that the member is
subjected to uniform load of 0.6 t/m'. Therefore, it is necessary to expose the result
thoroughly by dividing all members into smaller members.
With the intention of the above point, the frame members is now re-inputted into smaller
members, so that it becomes NM = 12, and NJ = 13. The results are then shown in
Figure 6.6 and 6.7.
The displayed geometry in chart window (only available in Excel 2003) from re-inputted data is
shown below:

Figure 6.6: FRAME2D input-output form Re-inputted

Figure 6.7: Charts of displacement, diagrams of bending moment,
shear force and axial force Re-inputted
The results shown in Figure 6.7 and 6.8 are much better than in Figure 6.5 and 6.6. Thus, it is
convenient to divide into smaller members and for this example at least 4 smaller members with
the same length for each member. Avoid odd divisors (e.g. 5, 7, 9 etc.) so that the results in the
middle of span can be exposed.
6.2 Sign Convention for Diagram

To depict diagrams of moment, shear force and axial force in members, program refers to the
sign convention as follows:

Positive direction
d is as
a shown in the above picture
p and ruled
r as follo
ows:
- Bending mom
ment (M) is given a possitive sign if causing the
e upper fiberrs of membe
ers is
ompressed or
co o to curve upward,
u and is given a negative
n sign
n if to curve downward.
- Shear force (V)

( is given a positive sign
s if rotate
es the memb
ber clockwisse and is givven a
egative sign if rotates th
ne he member counter-clock
c kwise.
- A
Axial force (N
N) is given a positive sig
gn if pulls the
e member apart, and is given a neg
gative
siign if pushes
s the membe
er.
This con
nvention is adopted
a in order
o to havve same inte
erpretation on
o the chartt depiction of
o the
applicatio
on programs
s in the next few chapterrs, thus it is more easily interpreted and also de
esign-
oriented. The progrram automa
atically diffe
erentiates th
he momentss and force
es that occcur in
memberss section, which is blue
e line inside the curve fo
or the posittive direction
n and red fo
or the
e direction. The
negative T example
e below show
ws the coloriing for bending momentt result:

gspot.co.id
6.3 Application
FRAME2D is now applied to the following beam problem:
The beam can be seen to have three main members that are, roll-load segment, load-roll
segment, roll-fixed node, have lengths 3, 3 and 5 meters, respectively. In practice however it
needs to divide the beam into smaller members for exposing the result more accurately.
FRAME2D is then inputted from which each of main members is divided into 4 equal parts. So,
there are 12 members and 13 joints. When the program has run, the charts will be shown as
below.


Excel For Civil Engineers Preview PDF

Uploaded by

Copyright:

Available Formats

Excel For Civil Engineers Preview PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Excel For Civil Engineers Preview PDF

Uploaded by

Copyright:

Available Formats

PREVIEW ONLY

Translated and re-produced in Nov, 2015.

An Introduction to: Excel for Civil Engineers 2

CHAPTER 1: BASICS OF EXCEL ................................................................................. 6

An Introduction to: Excel for Civil Engineers 3

An Introduction to: Excel for Civil Engineers 4

An Introduction to: Excel for Civil Engineers 5

1.1 Worksheet and Workbook

An Introduction to: Excel for Civil Engineers 6

1.2 Data Type

b. Numbers: numerical data: 1,2,3,0, -1, -2,4,5.85...

An Introduction to: Excel for Civil Engineers 7

Table 1.1: Excel Basic Commands

An Introdduction to: Excel

An Introdduction to: Excel

An Introduction to: Excel for Civil Engineers 10

=B4*B5 multiplying the data in cell B4 to cell B5

=SUM(D2:D4) summing the number of cells D2, D3 and D4

An Introduction to: Excel for Civil Engineers 11

Table 1.3 Excel Built-in Functions

A. Math and Trigonometry

An Introduction to: Excel for Civil Engineers 12

An Introduction to: Excel for Civil Engineers 13

1.5 Array Formula

An Introduction to: Excel for Civil Engineers 14

Here are some other examples of array formula:

{=SUM(A2: B2*A3:B3)} is a single result

{=TREND(A2:A6,B2:B6)} is a separately 5 results

9 Slope Y-intercept 9 Slope Y-intercept

An Introduction to: Excel for Civil Engineers 15

Table 1.4: Custom format for numbers

An Introduction to: Excel for Civil Engineers 16

Table 1.5 Error Messages

An Introduction to: Excel for Civil Engineers 17

1.9 Making Charts

An Introduction to: Excel for Civil Engineers 18

Figure 1.4: Charts group displays chart type options

An Introduction to: Excel for Civil Engineers 19

An Introduction to: Excel for Civil Engineers 20

The joint coordinates are:

An Introduction to: Excel for Civil Engineers 21

An Introduction to: Excel for Civil Engineers 22

An Introduction to: Excel for Civil Engineers 23

An Introduction to: Excel for Civil Engineers 24

Draw a simple picture of one floor building as below:

An Introduction to: Excel for Civil Engineers 25

The jointt coordinates

An Introdduction to: Excel

An Introduction to: Excel for Civil Engineers 27

An Introduction to: Excel for Civil Engineers 28

An Introduction to: Excel for Civil Engineers 29

An Introduction to: Excel for Civil Engineers 30

An Introduction to: Excel for Civil Engineers 31

Microsoftt Visual Bas

1.11.1 Creating Macro

An Introdduction to: Excel

ode created above, Sele

An Introdduction to: Excel

Objects are classifie