Abstract — Climate Change Impact Analysis
and uncertainty in climate projections. Concurrently, opportunities are emerging through advanced monitoring, innovative experimental techniques, computational modeling, climate integration, and reliability-based design. By extending classical bearing capacity models, this study integrates USM to more accurately predict geostructure performance. Analytical insights, supported by case studies, demonstrate the influence of rainfall-induced infiltration on slope stability, shallow foundation capacity, and column-supported embankments. Results reveal that suction enhances soil strength but may diminish rapidly during infiltration, heightening failure risk. The study advocates embedding USM into design codes, modeling frameworks, and early-warning systems to move from reactive to proactive resilience. Bridging theory and practice, it provides a pathway for adapting geotechnical systems to climate variability and ensuring long-term infrastructure durability.
Analytical & Numerical Modeling Methodology
Climate change is reshaping the frequency, intensity, and duration of weather events globally. In many regions, particularly those under tropical climates with residual soils, slope failures are increasingly triggered by climate-driven events such as prolonged droughts, extreme rainfall, and rising temperatures, which are no longer isolated anomalies but part of an increasingly volatile climate system ( Vardon, 2015 ; Yasuhara and Bergado, 2022 ; Costa et al., 2023 ; Kandalai et al., 2023 ; Bridges, 2024 ). For example, climate change has significantly altered rainfall patterns, leading to more intense storms occurring less frequently. This shift, coupled with increased atmospheric temperatures and moisture capacity, results in extreme rainfall events that greatly affect slope stability ( Au, 1998 ; Toll et al., 2008 ; Loveridge et al., 2010 ; Rahimi et al., 2011 ; Scaringi and Loche, 2022 ). Specifically, in tropical regions, intense rainfall infiltration reduces matric suction in residual soils, diminishing shear strength and making slopes more prone to failure. In contrast, arid and semi-arid regions face challenges associated with collapsible or dispersive soils, where wetting after prolonged droughts can induce sudden settlement or erosion. These diverse regional contexts highlight the universal applicability of unsaturated soil mechanics while underscoring the need for tailored design and adaptation strategies. These phenomena significantly affect the behavior and performance of geotechnical infrastructures ( Salimi and Al-Ghamdi, 2020 ; Pantoja Porro et al., 2025 ).
Finite Element Results & Code Compliance
Understanding how these regional processes translate into soil-scale mechanisms is crucial. The degree of infiltration and suction loss is strongly influenced by soil properties such as permeability, water retention characteristics, and hydraulic conductivity. These factors control how quickly pore-water pressures fluctuate during wetting and drying cycles. In climates with pronounced dry and wet seasons, this leads to frequent instability in areas composed of unsaturated soils. At the core of this challenge lies the unsaturated zone of the soil, often overlooked in traditional design but critical in modern geotechnical engineering. This vadose zone, located between the ground surface and the water table, comprises partially saturated soils commonly found beneath slopes, embankments, pavements, and shallow foundations. The hydro-mechanical properties of these unsaturated soils make them highly responsive to climate-driven changes in moisture and temperature, posing risks to stability, deformation, and serviceability ( Jardine, 2020 ; Culligan et al., 2019 ; Rouainia et al., 2020 ; Illés and Nagy, 2022 ; Ng et al., 2024 ; Pham et al., 2024 ).
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
Advantages of the Climate Change Impact Analysis 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