Calculation of Foundation Deformation under Circular Vertical Loading Using Triangular Axisymmetric Elements

This paper aims to investigate the computational methodology for foundation deformation under circular vertical loading using triangular axisymmetric elements. During the calculation process, it is essential to consider factors such as soil material properties, element material parameters, and element dimensions. The implementation typically involves defining material constitutive models using functions like `materialPropertyMatrix()` and specifying element geometry through nodal coordinates. Prior to computation, preprocessing steps include mesh generation using algorithms like Delaunay triangulation and boundary condition configuration through constraint application functions such as `applyBoundaryConditions()`. Key computational aspects involve solving the system of equations using finite element methods, where the stiffness matrix assembly is performed with functions like `assembleGlobalStiffness()`. The results obtained through this study provide enhanced understanding of foundation deformation mechanisms under circular loading conditions and offer valuable references for related engineering practices, with potential output visualization using `contourPlot()` functions for displacement fields.