Category: Sample Programming Codes

MATLAB

Modified Newmark Sliding Block (Deformable Sliding Block) For Estimating Seismic Lateral Movement of Slopes

Posted on
 

 Introduction Practitioners often perform Newmark Sliding Block analyses to predict permanent lateral displacements of slopes under earthquake loading. In this method, sliding of a rigid block of soil mass is computed by double integration of acceleration over the length of time during which the yield acceleration (i.e., friction coefficient) is exceeded. For better illustration purposes,… Read More

OPENSEES TCL Code

Beam on Nonlinear Winkler Foundation (BNWF) Model for Soil-Wall/Soil-Pile Interaction Analysis – TCL File

Posted on
 

 The current engineering practice adopts the concept of beam on a nonlinear Winkler foundation (BNWF) in order to perform displacement-based analysis on piles, cantilever walls, and bulkhead systems. This post provides an OpenSees TCL file to develop the BNWF model of a beam supported on a series of nonlinear springs. Similar to other OpenSees TCL… Read More

OPENSEES TCL Code

Single Pile Pushover Analysis Using Continuum Modeling Method – TCL File

Posted on
 

 This Tcl code creates full-scale three-dimensional continuum model of single pile embedded in a layered soil profile. Soil-pile interaction including formation of gap and slippage was simulated by using the BeamContact3D element. The program is capable of performing a lateral pushover analysis on a single pile. The contact elements at the interface of soil and pile… Read More

Excel VBA (Macro)

Macro for computation of acceleration, velocity and displacement response spectra

Posted on
 

 The VBA macro implemented in the enclosed spreadsheet calculates: (1) the time histories of velocity and displacement of the input acceleration time history, and (2) the spectral responses of acceleration (Sa), velocity (Sv), and displacement (Sd). The code is formulated based on the Newmark-bata numerical integration method. Note that C++ format of this code is… Read More

C++ Code

Calculation of pile end bearing capacity and skin friction using CPT data based on the approach of Bustamante and Gianeselli (1982) – C++ Code

Posted on
 

 As you all know one of the substantial advantages of Cone Penetration Test (CPT) is that the test provides a continuous profile of soil response. This is why CPT methods generally give superior predictions of axial pile capacity compared to most conventional methods (Robertson, 2015). Bustamante and Gianeselli (1982) proposed a CPT-based analysis procedure based on… Read More

Excel VBA (Macro)

Macro for calculating pile downdrag force based on Neutral Plane Method

Posted on
 

 Long-term consolidation, seismic dry sand and liquefaction-induced settlements of a soil deposit can impart significant amount of downward load to piles. This load, called downdrag or negative skin friction, is caused by the friction between the pile and the downward moving soil. In my previous post (click here), I explained how to calculate downdrag load… Read More

Excel VBA (Macro)

Macro for calculating pile downdrag force per AASHTO LRFD Bridge Design Manual

Posted on
 

 Long-term consolidation, seismic dry sand and liquefaction-induced settlements of a soil deposit can impart significant amount of downward load to piles. This load, called downdrag or negative skin friction, is caused by the friction between the pile and the downward moving soil. Per AASHTO LRFD Bridge Design guidelines, downdrag force should be determined as follows:… Read More

OPENSEES TCL Code

Continuum model of a single ground improvement element under seismic loading – TCL file

Posted on
 

 OpenSees implementation of a single ground improvement (GI) element in a continuum soil domain is described here. The Unit Cell concept is adopted to develop the ground improvement layout. The unit cell consists of a single GI element located in the center of its equivalent influence zone. The model represents a GI element located on… Read More