Spreadsheet for SPT to Bearing Capacity using Various Equations

Various Geotechnical Engineers have worked out on finding bearing capacity from SPT N value.  You can download its LibreOffice Calc Spreadsheet or Excel Spreadsheet.  If you have any suggestion or opinion, then, use comment box below. The formulas are shown in below figure:

SPT to Bearing Capacity
SPT to Bearing Capacity

Analysis of Isolated Footing using Elastic Foundation (Springs) in Abaqus

  1.  Choose Units, i.e., MKS or FPS. For this problem, we will use N, m and seconds as units.
  2. Make Rectangular footing in Sketch in Parts. Give it thickness as we want to model it as Solid Plate. Take footing dimension as 2x2x0.5
  3. Define Material as Concrete density (23.60) (note: correct this value when you practice it to 23.60e3), Elasticity (20e6) and Poisson’s ratio (0.2).
  4. We need to partition the plate so that we get point for application of load in centre. Tools>Partition in Part Module.
  5. Assign instance in assembly and then define interaction of Elastic Foundation in interaction module. Give value of Spring as 1. Note that the purpose of taking simplified value of load, dimension and spring is to get easy to interpret value of deflection to verify calculations.
  6. Mesh the footing from mesh module.
  7. Define load application step and apply concentrated load of 4 in negative 3-direction.
  8. Create Job and run the analysis.
  9. Open Result and check maximum deflection. As applied pressure is 4/4 = 1 and resisting stiffness per area is also 1 so from F = K*x, the resulting deflection will be 1 units. Again note that such values are pre-decided to have an easy cross-check on calculations.
  10. Check plots and deflection contour. And save the database cae file.

CAE file of abaqus can be downloaded from: CAE file

How to Download Accelerogram for given Earthquake Event?

Following are steps for downloading acceleration time history for any earthquake record:

  1. Open http://service.iris.edu/irisws/timeseries/docs/1/builder/
  2. Check station info from http://ds.iris.edu/wilber3/find_event
  3. You can find more about parameters of timeseries inputs from find event page of wilber3
  4. Network = II, Station = AAK, Location = 00, Channel = BH1, Start time = 2015-10-26 09:11:18, End time = 2015-10-26 09:21:18, Correction = Perform Instrument Correction, Frequency Limits = 0.001-0.02-0.02-9. Check differentiate, taper = 0.1, taper window = cosine, output = plot
  5. The url link is generated below, press to open in new tab. A time history plot will be opened in new window, which, you can save as image.
  6. Change the output type to ASCII: 2 column format, and click the below link to download the data.
  7. You can also change to audio to hear calculated sound of earthquake.
  8. Change channel to BH2 and then BHZ to download other components.
  9. Open the text file in spread sheet to process data.
Timeseries Data form
Timeseries Data form


Below is video explaining above steps: FB page of Engineering Programs

How to Download Seismogram for Given Earthquake Event?

Following are steps for downloading seismogram from IRIS Wilber3 for any earthquake event from any station:

  1. Open http://ds.iris.edu/wilber3/find_event
  2. You can view recent requests to download or Load Event Data.

    Find Earthquake Event and Station Info at Wilber 03
    Find Earthquake Event and Station Info at Wilber 03
  3. Select Past 12 months, M4.0+ in combo box of Load Event Data and Specify start and end date as 2015-10-26
  4. Scroll down to select event at date (utc) 2015-10-26 09:09:42 and click it to open.

    Station Details for Given Event
    Station Details for Given Event
  5. You can check coordinates, date, depth, and Related event page. Select the station for which you want to download data. Use map to zoom for stations.
  6. In the list below, click the station for which you want to download data.
  7. Select AAK. A popup will appear showing seismogram, phase arrival, time range. You can also select the instrument to preview its data. You can change time range to customize the seismograph. Close this pop up.

    Station Details for given Earthquake Event Showing Seismograms
    Station Details for given Earthquake Event Showing Seismograms
  8. With station AAK selected, press Show Record Section in side pan on right to see detailed seismograms.

    Show Record Section shows recorded seismograms at different stations.
    Show Record Section shows recorded seismograms at different stations.
  9. Press Request Data button, give your desired time range, data output format (SAAC or miniseed), bundle type, and fill the request info along with your email to let know you when the download file is ready.

    Request Data form helps in downloading earthquake record.
    Request Data form helps in downloading earthquake record.
  10. Press Submit button and check save these settings as default.

Below is video tutorial of above steps:

Compiling and Using CMU Sphinx in GNU Debian Linux

CMU Sphinx, also called Sphinx in short, is the general term to describe a group of speech recognition systems developed at Carnegie Mellon University. These include a series of speech recognizers (Sphinx 2 – 4) and an acoustic model trainer (SphinxTrain). Here are brief steps for compilation and usage of CMU Sphinx in Debian.

  1. Download sphinxbase and pocketsphinx from github:
  2. Install prerequisite packages  as:
  3.   Add write access to /usr/local for current user as:
  4. Export following paths:
  5. Change directory to sphinx base and compile it:
  6. In case of any update:
  7. Then change directory to pocketsphinx and compile it:
  8. In case of any update:
  9. To test run below command and wait for ready. Say a word in microphone and it will guess and type that word:

Recording Voice from Mic in Matlab

Below is sample code for recording voice from mic in matlab:



RampCalc0.1 – Parametric input for ramp in FE.

This program calculate data of beam, col and slab for ramp. It creates CSV file that can be used in sap2k or any FE package. It requires gnuoctave/matlab to work. You can download it from here:

Scripts Info:

  1. spiral.m = calculate and plot ramp in octave/matlab
  2. beams.m  = calculate data of beam for ramp to be used in FE
  3. cols.m    = calculate data of columns for ramp to be used in FE
  4. slabs.m   = calculate data of triangular slab for ramp to be used in FE

Background Theory:

Circle plot requires x = r*sin(t) and y = r*cos(t) where t increments from 0 to any given radians. If we increment z from 0 to given height h with some increment, and plot with x and y of a circle,  we will get a 3D spiral curve. Two such spiral curve can be used to calculate data point for ramp. Below is the sample code. You can run this by copy and paste directly in gnu octave or matlab or  create script out of it.

Ramp plot in octave
Ramp plot in octave

By using, beams, cols and slabs script function, csv data will be created that can be formatted as per any finite element code. Here I formatted it for sap2K.

Ramp view in sap2k
Ramp view in sap2k

A spreadsheet with sample sap2k format is attached with this code. Below is video tutorial explaining usage of this code:


Civil Quiz Engine

You can run this quiz engine from below Quiz Main Page link.

The purpose of this project is listed in following points:

  1. Preparation for Competition and Qualification Tests
  2. Clarify Civil Engineering Concepts.
  3. Hacking Civil Engineering Knowledge for better understanding and Learning.

You can run this quiz engine from Quiz Main Page

Lectures on Construction Project Planning, Management and Engineering Economics Spring 15

Followings are the lecture series on Project Planning, Management and Engineering Economics for Spring 2015. Main pdf lectures from 03 to 10 are provided by Dr. Prof. Syed Muhammad Ali. While the rest of lectures, further readings and matlab m files are provided by Me.

Please note that these are in preparation. Most of the content collected from various sources over internet.

Construction Project Planning & Management:

Lecture 1: Intro

Lecture 2: WBS, Kaizen

Lecture 03a: Project Management Process Group

Lecture 03b: Project Manager’s Role

Lecture 04a: Planning Commission

Lecture 04b: Project Life Cycle

Lecture 05: Project Appraisal

Lecture 06: Scheduling Resources

Lecture 07 & 08: ISO9000, Auditing Process, and ABM 

Lecture 08: Further Readings 01: Manual for Development Projects

Lecture 08: Further Readings 02: Approval of Development Projects

Lecture 08: Further Readings 02a: Project Approving Authorities in Pakistan

Lecture 08: Further Readings 03: Recent News about increase in CDWP limit

Lecture 08: Further Readings 04: Previous News about increase in CDWP limit

Lecture 08: Further Readings 05: Organization Structure of Ministries in KP Province

Lecture 08: Further Readings 06: Guidelines for Project Management

Lecture 08: Further Readings 07: Handbook of Planning Commission

Lecture 08: Further Readings 08: Critical Chain Method

Lecture 08: Further Readings 09: PCDA

Lecture 08: Further Readings 09a: SWOT Analysis

Lecture 08: Further Readings 10: Time Driven Model for Choosing a Decision

Lecture 08: Further Readings 11: Project Life Cycle in Pakistan

Lecture 08: Further Readings 12: Project Management Videos

Lecture 08: Further Readings 13: Info on ABM

Engineering Economics:

Lecture 09: Engineering Economics Introduction

Lecture 09: Further Readings 01: Engineering Economics Introduction

Lecture 09: Further Readings 02: Engineering Economics Formulae

Lecture 09: Further Readings 03: Net Present Value

Lecture 09: Further Readings 04: Interest

Lecture 09: Further Readings 05: Usury

Lecture 09: Further Readings 06: Salvage Value

Lecture 10: Time Value of Money

Lecture 11: Cost-Benefit Analysis

Lecture 11a: SPI & CPI

Lecture 11b: Break-Even Analysis

Lecture 11c: Internal Rate of Return

Lecture 11: Further Readings 01: ADB on B/C

Lecture 11: Further Readings 02: B/C & IRR on Bridge

Lecture 11: Further Readings 03: IRR

MATLAB m files for Economic Analysis:

Matlab m files for formulae of engineering economics can be found here. You can upload them to octave-online.net to run them. Or you can check it on your desktop using gnu octave or matlab.

Primavera and Vico Constructor Video Tutorials:

Learn project planning with engineering applications using Primavera P3 and BIM – Vico Constructor from Youtube Channel of Engr. Qasim Shinwari at https://www.youtube.com/user/913shinwari

QGIS Lab10 Fall14: QGIS Raster Analysis


This is tenth video tutorial of QGIS. It is about raster analysis. The objective is to learn and use raster analysis in qgis.


  1. Load analysis.qgs
  2. Set project crs to 33S
  3. Add four settelite files and create virtual raster catalog. Raster > Miscellaneous > Build Virtual Raster (Catalog).
  4. Save in exercise_data/residential_development as aerial_photos.vrt
  5. Raster > Projections > Warp (Reproject). to 70300
  6. Raster > Miscellaneous > Merge.
  7. Changing Raster Layer Symbology and using View > Toolbars > Raster
  8. Terrain Analysis (hillshade, Slope and Aspect)
  9. Using the Raster Calculator, aspect@1 <= 90 OR aspect@1 >= 270, save it in aspect_north.tif and also calculate slope_lte2.tif and slope_lte5.tif.
  10. Write criteria for location: ( aspect_north@1 = 1 AND slope_lte5@1 = 1 ) OR slope_lte2@1 = 1 and save it as all_conditions.tif.
  11. Simplifying Raster: all_conditions_sieve.tif and Threshold and Pixel connections values to 8
  12. Open the Raster Calculator again, and build this expression: (all_conditions_sieve@1 <= 0) = 0 and save output as: all_conditions_simple.tif
  13. Raster > Conversion > Polygonize (Raster to Vector). Field name = suitable and save it as: all_terrain.shp
  14. For filer of suitable = 1, save it as suitable_terrain.shp
  15. Raster > Conversion > Rasterize (Vector to Raster) for suitable_terrain.shp, width = 837 and height = 661
  16. Use the Intersect tool (Vector -> Geoprocessing Tools) to create a new vector layer called new_solution.shp which contains only those buildings which intersect the suitable_terrain layer.

Below is the full video lecture. Use full screen for clear content of video. The download link of below video is ogg video or mp4 video. You can play it in vlc player or open it in web browser.

The example files that are used in video can be downloaded from here .


  1. Chapter 08 of QGIS 2.6 Training Manual