Abstract — Earthquake Induced Liquefaction
experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Analytical & Numerical Modeling Methodology
Earthquake-inducedliquefaction of soils poses a serious georisk in geotechnical designs, construction and the application of geotechnical structures around the world. In this study, the applicability of three soft computing models for liquefaction classification, a topic of significant importance within the fields of geotechnical and earthquake engineering has been evaluated. Twelve input parameters are used to classify the liquefaction potential for 234 data sets collected from an earthquake-induced liquefaction prone granular material environment. For developing the SVM_Poly, SVM_RBK models, an extensive number of trials were conducted using various combinations of C and d for polynomial kernels and C and ∂ for radial basis function kernel-based support vector machines (SVMs) utilizing user-defined parameters. In the same way, several experiments were conducted with a fixed value of C and ∂ kernel specific parameters in order to determine an appropriate value of error-insensitive zone (∋).Similarly, for the random forest classifier (RFC) model, the number of variables used (m) and the number of trees to be grown (k) are two user-defined parameters. These optimum values of m and k parameters are fixed using trial and error process and the same fixed values. The best model was developed as evidence from the confusion matrixes and statistical indicators. The calculated values of confusion matrixes and statistical indicators for training and testing shows that an accuracy of 0.89 indicates the model is correct in its predictions 89% of the time. A sensitivity of 0.85 signifies the model correctly identifies 85% of actual positive instances, while a specificity of 0.94 implies correct identification of 94% of actual negative instances. A precision of 0.94 suggests that when the model predicts a positive instance, it is correct 94% of the time. The Phi Correlation Coefficient, with a value of 0.82, indicates a strong positive correlation between predicted and actual values.Furthermore, the model exhibits a Mean Absolute Error (MAE) of 0.2351, reflecting a relatively low average error in predictions. The Root Mean Squared Error (RMSE) value of 0.3115 indicates better accuracy in predicting the target variable.Finally, all the developed models exhibit promising performance across various evaluation metrics, with low error measures (MAE and RMSE), high accuracy, and strong performance in correctly identifying both positive and negative instances, as evidenced by sensitivity and specificity. The high precision and Phi Correlation Coefficient further affirm the reliability and accuracy of the model’s predictions. However, among the three models FRC model is the best for classifying the RFC models using a comprehensive set of seismic and soil parameters to accurately classify earthquake-induced liquefaction potential, with the RFC model demonstrating superior predictive performance.
Finite Element Results & Code Compliance
Liquefaction of granular materials refers to a phenomenon where a saturated or partially saturated granular material, such as sand or silt, temporarily loses its strength and behaves like a liquid under stress, often as a result of seismic activity, such as an earthquake 1 . Liquefaction occurs when the pore water pressure within the soil increases, reducing the effective stress that holds the soil particles together 2 . Sands and silts with low cohesion are most susceptible due to the voids between soil particles which are filled with water, which can increase pressure during shaking 3 .Vibrations from an earthquake cause the soil to compact, which pushes water into the voids, leading to increased pore water pressure 4 .As pore water pressure rises, the ability of the particles to support weight (effective stress) diminishes 5 .Structures can settle, tilt, or collapse due to the loss of support from the liquefied soil 6 . Soil may behave like a fluid, causing landslides or flow of material 7 . Ground displacement can cause large, horizontal movements and excess water and sand may be ejected to the surface, forming mounds 8 . Techniques like vibro-compaction or dynamic compaction increase the density of the soil, reducing the likelihood of liquefaction 5 . Installing drainage systems helps reduce pore water pressure by allowing water to escape more easily 3 .
Scope & Practical Applications
- High-rise building and spatial frame analysis under seismic & wind loading
- Finite element continuum modeling in ABAQUS for non-linear material plasticity
- Classical lamination theory (CLT) and failure prediction for composite structures
- Spatial hazard mapping, landslide susceptibility, and catchment terrain analysis in QGIS
- Detailed bar bending schedules (BBS), structural drawings, and calculation booklets
3D Frame Analysis
Spatial modeling in ETABS & STAAD.Pro
Non-Linear FEA
Constitutive damage models in ABAQUS
QGIS Geomatics
DEM terrain modeling & catchment maps
Stiffness Solver
Direct matrix assembly in MATLAB
Code Compliance
IS 456, IS 1893, ACI 318, Eurocodes
Stress Envelopes
Critical axial, shear & moment vectors
Engineering Software & Technologies
Built using industry-standard engineering analysis platforms, numerical code solvers, and geomatics tools.
| Software Platform | Analytical Function & Role | Engineering Tier |
|---|---|---|
| ETABS v21 | Multi-story building analysis, response spectrum, P-Delta, wind & seismic lateral loads | Structural Core |
| STAAD.Pro Connect | 3D space frame analysis, steel connection design, foundation reaction extraction | Steel & Concrete |
| ABAQUS Standard | Non-linear continuum mechanics, concrete damaged plasticity, composite lamina failure | Non-Linear FEA |
| QGIS Desktop | Digital Elevation Model (DEM) hydrological mapping, slope stability raster overlays | Spatial Geomatics |
| MATLAB R2024b | Global stiffness matrix assembly, modal eigenvalue extraction, composite ABD matrix calculation | Mathematical Solver |
Analytical Engineering Workflow
Systematic integration from CAD architectural geometry and GIS terrain mapping to finite element simulation and code validation.
Geometry & BIM
(Concrete/Steel/FRP)
Spatial Analysis
Non-Linear Mesh
& Foundation Coupling
& Detailing
Verified Structural Diagrams & FEA Stress Contours
Advantages of the Earthquake Induced Liquefaction Project
Why this project delivers superior academic grades and provides high technical evaluation standing.
Multi-Software Cross-Validation
Integrates ETABS, STAAD.Pro, and ABAQUS modeling to verify results across both linear and non-linear finite element solvers.
Rigorous Code Compliance
Fully compliant with IS 456, IS 1893 (Seismic Criteria), IS 800 (Steel Design), ACI 318, and Eurocode provisions.
Non-Linear Material Plasticity
Captures crack propagation, yield surfaces, and post-peak softening behavior using ABAQUS Concrete Damaged Plasticity.
Composite Lamination Analysis
Evaluates Classical Lamination Theory (CLT), ABD matrices, and Hashin/Tsai-Wu failure envelopes for high-strength FRP laminates.
Spatial Geomatics Coupling
Combines high-resolution QGIS Digital Elevation Models with slope stability analysis for accurate terrain hazard zoning.
Deterministic Matrix Solvers
Custom MATLAB scripts assemble global stiffness matrices [K] and extract natural frequencies without relying purely on black-box GUIs.
Second-Order P-Delta Verification
Accurately assesses geometric non-linearity and sway amplification induced by gravity and lateral wind/seismic combinations.
Detailed Structural Detailing
Delivers complete Bar Bending Schedules (BBS), ductile detailing drawings per IS 13920, and steel connection drawings.
Direct Industry Relevance
Builds competencies demanded by top structural design consultancies, EPC contractors, and infrastructure developers.
Complete Turnkey Package
Includes full model files (.EDB, .STD, .CAE), MATLAB code, calculation sheets, IEEE-format report, and viva defense presentation.
Sample MATLAB Code, FEA Output & Member Schedules
Representative excerpts illustrating the finite element stiffness assembly, command solver output, and structural capacity envelopes.
[*] Initializing ETABS API Engine v21 & STAAD.Pro Analytical Space Frame... [✓] Dead Load (DL), Live Load (LL), Wind (WL), and Seismic (EQX/EQY) Loaded ──────────────────────────────────────────────────────────────────────── [MODAL ANALYSIS] Mode 1 (Translational X): T1 = 0.842 sec (Mass Part: 81.4%) [MODAL ANALYSIS] Mode 2 (Translational Y): T2 = 0.791 sec (Mass Part: 82.8%) [MODAL ANALYSIS] Mode 3 (Torsional Theta): T3 = 0.512 sec (T/Tx = 0.608 < 0.90 OK) [DRIFT CHECK] Max Inter-Story Drift: Δ = 0.0018 < 0.0040 (IS 1893 Clause 7.11.1 PASSED) [P-DELTA EFFECT] Maximum Stability Index θ: 0.041 < 0.10 (Secondary Moments Negligible) ──────────────────────────────────────────────────────────────────────── [✓] Design Envelope Complete: Zero Overstressed Members across 1.5(DL + LL) & 1.2(DL + LL ± EL)
[*] Generating critical member capacity envelope across factored combinations...
| Member Tag | Section Designation | Axial Pu (kN) | Bending Mz (kNm) | Shear Vy (kN) | Capacity Ratio | IS Code Status |
|---|---|---|---|---|---|---|
| Col C1-G1 | 450 x 450 mm M40 | 2,410.8 | 184.2 | 84.5 | 0.742 | PASSED |
| Col C2-G1 | 500 x 500 mm M40 | 3,180.4 | 218.6 | 112.4 | 0.814 | PASSED |
| Beam B101 | 300 x 600 mm M35 | 42.5 | 264.8 | 168.2 | 0.688 | PASSED |
| Beam B102 | 300 x 600 mm M35 | 38.1 | 295.4 | 182.1 | 0.762 | PASSED |
[✓] Structural Analysis Complete: Zero overstressed members. Design capacity factor > 1.0 throughout.
Frequently Asked Questions
Project Deliverables & Support
Everything needed from project inception to final viva defense — delivered with expert technical support.
Model Files
Input Decks
Scripts
Report
Slides
Drawings
Sheets
Q&A
Support
Extension