An information-analytical system for finite-difference modeling of dynamic stress concentrations in composite materials under impact, tension, and bending loads
Zhuzbayev Serik1, Sarsenbay Magzhan2, Adilova Aknur3, Sabitova Diana4, Khabdolda Bolat5 and Karin Abylay6
PhD Student, Department of Information Systems,L.N. Gumilyov Eurasian National University,Astana,Republic of Kazakhstan2
Senior Lecturer, Department of Information Systems, Faculty of Technology,M.Kh. Dulati Taraz University,Taraz,Republic of Kazakhstan3
Associate Professor, Department of Mathematics, Physics and Informatics,Sh. Ualikhanov Kokshetau University,Kokshetau,Republic of Kazakhstan4
Senior Lecturer, Department of Applied Mathematics and Informatics,Karaganda Buketov University,Karaganda,Republic of Kazakhstan5
Master’s Student, Department of Information Systems,L.N. Gumilyov Eurasian National University,Astana,Republic of Kazakhstan6
Corresponding Author : Sarsenbay Magzhan
Recieved : 08-November-2025; Revised : 19-June-2026; Accepted : 22-June-2026
Abstract
Composite materials (CM) are widely used in aerospace, transportation, and civil engineering due to their high specific strength, stiffness, and wear resistance. However, predicting their dynamic mechanical response under non-stationary loading conditions remains a significant challenge. Stress fields often exhibit singularities at corner points, contact discontinuities, and dissimilar material interfaces, which cannot be accurately resolved using classical homogenization approaches or conventional finite element methods (FEM). This paper presents a second-order accurate computational framework for analyzing the behavior of CM within a unified strip-halfplane configuration. Three mechanical scenarios are investigated: impact resistance (Problem A), tensile strength (Problem B), and bending behavior (Problem C). The proposed framework employs an explicit finite difference method (FDM) based on the method of spatial characteristics and a splitting technique. The numerical scheme is implemented in Python 3.12.3 using NumPy and Matplotlib libraries. A dedicated algorithm is developed to compute stress and velocity fields at singular points and is validated through a three-level grid convergence study. The results indicate that the stress concentration ratio at the singular point P reaches 3.0 times the applied boundary load in Problem A. Under tensile loading in Problem B, the normalized maximum normal stress, σ11k
, attains a value of 0.63 at point P. Furthermore, sensitivity analysis reveals that the extremal stresses in the half-plane material vary by a factor of 3.0 for strip embedment depths of 5h, 10h, and 15h. Grid convergence is confirmed with a root mean square error (RMSE) of 4.3 × 10⁻³ between the medium and fine grids. The observed order of convergence at the singular point P is 1.2, which is consistent with the second-order accuracy of the trapezoidal bicharacteristic integration scheme. These findings provide quantitative reference data for the design and analysis of composite structures containing strip inclusions subjected to impact, tensile, and bending loads.
Keywords
Composite materials, Dynamic stress analysis, Stress singularities, Finite difference method, Grid convergence, Strip-halfplane configuration.
Cite this article
Serik Z, Magzhan S, Aknur A, Diana S, Bolat K, Abylay K. An information-analytical system for finite-difference modeling of dynamic stress concentrations in composite materials under impact, tension, and bending loads. International Journal of Advanced Technology and Engineering Exploration. 2026;13(139):898-920. DOI : 10.19101/IJATEE.2025.121221464