• Computers and Concrete
• /
• v.30 no.4
• /
• pp.289-299
• /
• 2022

Based on the random cracking theory, the cylinder RVE model of reinforced concrete is established and the damage process is divided into three stages as the evolution of the cracks. The stress distribution along longitude direction of the concrete and the steel bar in the cylinder model are derived. The equivalent elastic modulus of the RVE are derived and the user-defined field variable subroutine (USDFLD) for the equivalent elastic modulus is well integrated into the ABAQUS. Regarding the tensile rebars and the concrete surrounding the rebars as the equivalent homogeneous transversely isotropic material, and the FEM analysis for the reinforced concrete beams is conducted with the USDFLD subroutine. Considering the concrete cracking and interfacial debonding, the macroscopic damage process of the reinforced concrete beam under four-point bending loading in the simulation. The volume fraction of rebar and the cracking degree are mainly discussed to reveal their influence on the macro-performance and they are calibrated with experimental results. Comparing with the bending experiment performed with 8 reinforced concrete beams, the bending stiffness of the second stage and the ultimate load simulated are in good agreement with the experimental values, which verifies the effectiveness and the accuracy of the improved finite element method for reinforced concrete beam.

• Structural Engineering and Mechanics
• /
• v.72 no.6
• /
• pp.713-722
• /
• 2019

Micromechanics is a technique for the analysis of composites or heterogeneous materials which focuses on the components of the intended structure. Each one of the components can exhibit isotropic behavior, but the microstructure characteristics of the heterogeneous material result in the anisotropic behavior of the structure. In this research, the general mechanical properties of a 3D anisotropic and heterogeneous Representative Volume Element (RVE), have been determined by applying periodic boundary conditions (PBCs), using the Asymptotic Homogenization Theory (AHT) and strain energy. In order to use the homogenization theory and apply the periodic boundary conditions, the ABAQUS scripting interface (ASI) has been used along with the Python programming language. The results have been compared with those of the Homogeneous Boundary Conditions method, which leads to an overestimation of the effective mechanical properties. According to the results, applying homogenous boundary conditions results in a 33% and 13% increase in the shear moduli G₂₃ and G₁₂, respectively. In polymeric composites, the fibers have linear and brittle behavior, while the resin exhibits a non-linear behavior. Therefore, the nonlinear effects of resin on the mechanical properties of the composite material is studied using a user-defined subroutine in Fortran (USDFLD). The non-linear shear stress-strain behavior of unidirectional composite laminates has been obtained. Results indicate that at arbitrary constant stress as 80 MPa in-plane shear modulus, G₁₂, experienced a 47%, 41% and 31% reduction at the fiber volume fraction of 30%, 50% and 70%, compared to the linear assumption. The results of this study are in good agreement with the analytical and experimental results available in the literature.

• Structural Engineering and Mechanics
• /
• v.71 no.6
• /
• pp.699-712
• /
• 2019

The reinforced concrete lining of hydraulic pressure tunnels tends to crack under high inner water pressure (IWP), which results in the inner water exosmosis along cracks and involves typical hydro-mechanical interaction. This study aims at the development, validation and application of an indirect-coupled method to simulate the lining cracking process. Based on the concrete damage plasticity (CDP) model, the utility routine GETVRM and the user subroutine USDFLD in the finite element code ABAQUS is employed to calculate and adjust the secondary hydraulic conductivity according to the material damage and the plastic volume strain. The friction-contact method (FCM) is introduced to track the lining-rock interface behavior. Compared with the traditional node-shared method (NSM) model, the FCM model is more feasible to simulate the lining cracking process. The number of cracks and the reinforcement stress can be significantly reduced, which matches well with the observed results in engineering practices. Moreover, the damage evolution of reinforced concrete lining can be effectively slowed down. This numerical method provides an insight into the cracking process of reinforced concrete lining in hydraulic pressure tunnels.

• Advances in Computational Design
• /
• v.9 no.1
• /
• pp.39-52
• /
• 2024

If the governing differential equation arising from engineering problems is treated as an analytic, continuous and derivable function, it can be expanded by one point as a series of finite numbers. For the function to be zero for each value of its domain, the coefficients of each term of the same power must be zero. This results in a recursive relationship which, after applying the natural conditions or the boundary conditions, makes it possible to obtain the values of the derivatives of the function with acceptable accuracy. The elastoplastic analysis of an inhomogeneous thick sphere of metallic materials with linear variation of the modulus of elasticity, yield stress and Poisson's ratio as a function of radius subjected to internal pressure is presented. The Beltrami-Michell equation is established by combining equilibrium, compatibility and constitutive equations. Assuming axisymmetric conditions, the spherical coordinate parameters can be used as principal stress axes. Since there is no analytical solution, the natural boundary conditions are applied and the governing equations are solved using a proposed new method. The maximum effective stress of the von Mises yield criterion occurs at the inner surface; therefore, the negative sign of the linear yield stress gradation parameter should be considered to calculate the optimal yield pressure. The numerical examples are performed and the plots of the numerical results are presented. The validation of the numerical results is observed by modeling the elastoplastic heterogeneous thick sphere as a pressurized multilayer composite reservoir in Abaqus software. The subroutine USDFLD was additionally written to model the continuous gradation of the material.