Elastic-plastic Model Research Articles

Influence of Material Degradation on Deformation of Paraglider during Flight.

The aim of this article is to determine experimentally and numerically the influence of material degradation on the deformation of a paraglider during flight. The presented method regards numerical modeling of pressure distribution over the wing and its effect on paraglider behavior; the considerations are preceded by experiments on three types of Polyamide 6.6 paraglider fabrics, subjected and not subjected to thermal, UV and flexing degradation. Scanning electron microscope (SEM) records allowed to determine the structural characteristics of the analyzed samples. Air permeability and mechanical tests are the input data for the computational simulations. When a pressure drop of 200 Pa is applied, all the analyzed samples are impermeable, except for those damaged by flexing. Thus, flexing damage has the greatest influence on the air permeability change among all considered aging factors. Aging caused by UV radiation has the greatest influence on mechanical properties. No major influence of thermal ageing on the mechanical properties of the considered samples is observed. Safety factors of the considered materials not subjected to degradation range between 3.94 and 6.00. Safety factor of fabric no. 1 subjected to the UV degradation is equal to 1.33; this result does not secure a safe usage of the considered material. The methodology described in this research can help to predict paraglider covering materials' behavior in flight; it assumes many cases, i.e., applying a new material or the material at any point of its life cycle. Thus, the practical implications of this model supported by numerical methods may result in saving time and cost in producing prototypes, as well as potentially assessing the safety of used wings. Future research activity can introduce the application of different elastic-plastic damage models to determine the paraglider behavior during collapse.

ПРОГНОЗУВАННЯ НЕСУЧОЇ СПРОМОЖНОСТІ ПЛИТНОГО ФУНДАМЕНТА ЗА ЧИСЛОВИМ МЕТОДОМ ГРАНИЧНИХ ЕЛЕМЕНТІВ

Construction is one of the leading branches of the national economy in the historical aspect of its development. The first design task is to determine the strength of building structures. Therefore, the study of the stress-strain state and related calculations are the most important in construction. The purpose of studying soil mechanics and foundation construction techniques is the calculation and construction of structures on or in soil. The main task is the construction of structures with a sufficient degree of reliability.The selection of an adequate theoretical model remains the main problem of soil mechanics. Indeed, the deformation of the dispersed granular material of the soil takes place during the mutual sliding of the grains, the rheology of the soil is complex, as evidenced by a large amount of experimental material. Today, the path of development of soil mechanics is related to the study of problems within the framework of the elastic-plastic dilatation model and the improvement of this model based on experiments. A mathematical model of a technical object at the micro level is a system of differential equations in partial derivatives, the exact solution of which can be obtained only in a few partial cases, therefore a discrete model is built using numerical methods that use the Poisson idea that the behavior of a complex model can be represented by the behavior of its individual component elements. The intensive development and widespread use of computers significantly brought fundamental mathematical problems closer to applied ones, and strengthened their mutual influence.The emergence of a new, powerful and general method of research - a numerical experiment, more than ever before closely connected the physical content of the problem, its mathematical formulation, numerical methods of calculation and modern computers. The work uses the numerical method of boundary elements. A promising way to develop foundations and foundation structures is to use the ratios of the theory of plastic flow, and the level of development of soil mechanics significantly affects the economy and reliability of the decisions made.

A coupled mechanical-electrochemical phase-field formulation for understanding the evolution of lithiated-silicon sponge

The unexpected formation of lithiated-silicon sponge during cyclic charge of lithium-ion batteries with silicon anodes has been observed in recent reports. The severe hollowing of silicon electrodes leads to an irreversible deformation and the generation of new surfaces, which greatly affects the structural integrity and cycling stability of the cells. However, the understanding of lithiated-silicon sponge still remains as a great challenge. In this work, a coupled mechanical-electrochemical model based on the phase-field method is proposed, which can describe the void evolution and the surface instability in Si electrodes. The model combines a description of the lithium diffusion, vacancy nucleation and annihilation, and an elastic-plastic constitutive model with concentration-dependent material properties. The simulation results reveal that the void within/on the silicon electrode tends to grow toward/along the electrode surface and toward other voids, leading to the electrode surface disruption and void merging, as well as the intra-electrode tunnel formation and electrolyte intrusion. Importantly, it is found that mechanical pressure can not only suppress the growth of a single void but also prevent the linkage and merging between multiple voids. The proposed theoretical framework and the findings of this work provide an in-depth understanding of the complicated mechanical-electrochemical coupling behavior of Si electrodes.

Finite element analysis of reinforcement for gas insulated switchgear

Considering the burgeoning requirement for reinforcing gas-insulated metal-enclosed switchgear, a three-dimensional elastic-plastic finite element model was established to examine the stress distribution before and after reinforcement. The results showed that, under an internal pressure of 3.8 MPa, the maximum stress of the defective sample surpassed the yield strength of the material, leading to its destruction. On the other hand, the maximum stress value of the reinforced sample, with a gel and a reinforcing plate, decreased to 145 MPa and was lower than the yield strength of the material. Moreover, the results of the hydrostatic test further demonstrated that the gel was effective in reducing the stress concentration at the defect within a certain range. In particular, the results indicated that the gel was capable of decreasing the maximum stress value by 33 MPa, thereby safeguarding the structural integrity of the switchgear.

Which asperity scales matter for true contact area? A multi-scale and statistical investigation

The true contact area between two surfaces is only a small fraction of the apparent macroscopic contact area; it governs many interfacial properties such as friction and contact resistance and depends sensitively on roughness. However, for real-world multi-scale surface topography, it is not clear which size scales of roughness govern the true contact area. This study investigates true contact area for a real-world surface that has been characterized across all scales from Angstroms to centimeters. Elastic and elastic-plastic contact is investigated using both a multiscale framework and a statistical roughness model. The multiscale method is a rough-surface contact-modeling technique based on Archard's stacked scales from a spectrum of the surfaces, which has shown promise when compared to previous experimental and numerical results. In contrast, statistical models assume that the asperities follow a defined height distribution and are in contact when taller than the mean surface separation. The results show that even the smallest scales can have a significant influence on the contact area, especially when the contact is elastic. However, when the contact is elastic-plastic, the influence of smaller scales can be limited depending on the character of the roughness. For self-similar, fractal-like roughness across some scales, the pressure tends to saturate at those scales. This work also explores the inclusion of scale-dependent yield strength. Both the multiscale and statistical models predict that the inclusion of scale-dependent strength causes the predicted contact area of the elastic-plastic models to come into closer agreement with that of the elastic model, especially when a wider range of size scales are included. In addition, both types of models predict that below a certain scale, smaller asperities flatten under contact pressure and will no longer influence the predicted contact area. Taken together, this work helps to guide the accurate modeling of rough-surface contact, and provides insights into which scales can be modified to improve performance in manufactured components.

Multi-objective optimization of sawblade multi-spot pressure tensioning process based on backpropagation neural network and genetic algorithm

ABSTRACT Resources scarcity has created strict requirements for efficient wood processing. Multi-spot pressure tensioning can effectively stabilize the sawblade and reduce cutting losses, but it is difficult to apply in practice due to the complexity of the process. In this study, the elastic-plastic solid model of multi-spot pressure tensioning was established, which was simplified to a thermal expansion shell model, and the mapping relationship between the two models was determined. The feasibility of the mapping method was verified by experiments. The backpropagation neural network (BPNN) was trained on a database composed of 8160 working conditions and the relative error was found to be less than 5%. Indenter displacement, pressing point quantity, and indenter radius were shown to change the degree of tension to varying extent. The pressing-point-distribution radius determines the development direction of sawblade performance. Increasing the number of pressing point circles can expand the adjustment range for sawblade performance. Both sawblade performance and tensioning energy consumption can be optimized by the genetic algorithm (GA). The optimal process parameters for different applications can be obtained. The combination of finite element method, BPNN, and GA can effectively optimize the multi-spot pressure tensioning process to improve the sawblade performance.

Multifactorial analysis of a gateroad stability at goaf interface during longwall coal mining – A case study

Purpose. Creating a generalized algorithm to account for factors (coal seam thickness, enclosed rock mechanical properties, the dimension and bearing capacity of artificial support patterns) causing a gateroad state under the effect of longwall face and goaf. Methods. The assessment of the gateroad stability is based on numerical simulation of the rock stress-strain state (SSS). The finite element method is used to find out the changes in the SSS of surrounding rocks at various stages of longwall mining. The elastic-plastic constitutive model and Hoek-Brown failure criterion implemented in codes RS2 and RS3 (Rocscience) are applied to determine rock displacements dependently on the coal seam thickness, enclosed rock strength, width and strength of artificial support (a packwall comprised of hardening mixture “BI-lining”). To specify the mechanical properties of the packwall material a series of experimental tests were conducted. A computational experiment dealing with 81 combinations of affecting factors was carried out to estimate the roof slag and floor heaving in the gateroad behind the longwall face. A group method of data handling (GMDH ) is employed to generalize the relationships between rock displacements and affecting factors. Findings. The roof-to-floor closure in the gateroad has been determined at the intersection with the longwall face and goaf dependently on the coal seam thickness, enclosed rock strength, width of the packwall, and strength of hardening material. It is revealed that the support material gains the strength value of 30 MPa on the 3rd day from its beginning to use which is fully corresponding to the requirements of protective element bearing capacity. The possibility of using untreated mine water to liquefy the mixture is proved, that allows simplifying and optimizing the solute mixing and pumping technology. Originality. This study contributes to improving the understanding of the factors that influence the stability of underground mining operations and highlights the importance of utilizing numerical simulations in optimizing mining designs. The impact of each factor on the resulting variable (decrease in cross-section of gate road by height) based on the combinatorial algorithm of structural identification of the model is estimated as follows: the packwall width is 48%, the thickness of coal seam is 25%, the strength of enclosing rocks is 23%, and the strength of the packwall material is 4%. Practical implications. The findings provide stakeholders with a technique to determine reasonable parameters for support and protective systems, and the predictive model developed can be used to mitigate potential instability issues in longwall mining excavations. The results have implications under similar geological settings and can be valuable for mine design and optimization in other regions.

Numerical Investigation on Reduction Characteristics of Welding Deformation According to Standoff Distance in Heat Sink Welding Process

Welding with a trailing heat sink process is proposed as a useful technique to reduce welding deformation without a significant loss of manufacturing efficiency. In particular, this method is effective than the existing mechanical or thermal correction method in thin plate steel structures with low stiffness. In this study, a three-dimensional thermal elastic-plastic finite element model using a 3-mm thick 304L stainless steel sheet was developed to simulate heat sink welding. The reduction characteristics of welding deformation according to the standoff distance between the welding heat source and the heat sink zone was investigated in the numerical simulation. The results confirmed that the longitudinal and angular deformations were significantly reduced compared with those by conventional welding. The longitudinal deformation was reduced by approximately 60%, from 5.50 mm to 2.22 mm, and the angular deformation was reduced by approximately 70%, from 1.68 mm to 0.51 mm. The overall results obtained in the numerical analysis indicated that welding deformation significantly decreased as the standoff distance was shorter owing to the temperature cooling effect of the welds.

ПРОТИЗСУВНА СПОРУДА ПІДВИЩЕНОЇ НЕСУЧОЇ ЗДАТНОСТІ ТА ЗНИЖЕНОЇ МАТЕРІАЛОМІСТКОСТІ

Recently, there has been active construction development of territories with difficult engineering and geological conditions. Construction on landslide and landslide-prone slopes leads to increased loads on individual areas, which can cause activation of old or the development of new landslides. As a result of the occurrence and activation of landslides, unacceptable deformations are formed, which lead to a violation of the stability of the slope, as well as the destruction of buildings and structures. Therefore, when performing construction works on slopes, it is necessary to provide measures for engineering protection of territories, buildings and structures located on or near slopes. On many sites, the placement of slopes and the erection of supporting prisms, buttresses and counterbanks is impossible. This is due to the need to observe the limits of alienation, the presence of residential and industrial buildings. In this case, the most rational means is the installation of pile retaining structures. They are optimal for shear strength of up to 15÷20 m. If it is necessary to contain deep landslides, bored piles are exposed to large bending moments, which leads to an increase in the number of piles or their diameters. As a result, the volume of drilling work and consumption of materials increases significantly. Therefore, issues of rationalization of retaining anti-slide structures are of particular practical interest. The paper calculation models were developed and numerical modeling of the interaction of elements of two structures of retaining anti-slide structures with the soil massif was performed. Numerical modeling was performed in the licensed geotechnical software complex Midas GTS NX, which implements the finite element method. The behavior of the soil in the calculation schemes was described by the Mohr-Coulomb elastic-plastic model. As a result of calculations of the stress-strain state of the elements of the retaining structures compatible with the soil, the coefficients of the reserve of stability were determined by the reduction method and the value of the internal forces (bending moments and transverse forces) arising in the elements of the structures. The reinforcement of the retaining elements was selected for the obtained internal forces using the "SCAD" software complex. On the basis of a comparison of the costs of the main materials necessary for the construction of structures, the economic feasibility of using a landslide protection retaining structure in the form of buttresses made of bored piles has been proven.

The practice of arranging a pit deeper than the level of the soles of the foundations of the surrounding buildings

A rational constructive-technological decision to arrange a fence of a deeper pit for a new building than the depth of laying strip foundations of the surrounding buildings has been substantiated and tested in practice. The decision to arrange the fence was tested based on the results of complex experimental and theoretical studies under the conditions of a wet loess massif. At same time, regulatory restrictions on additional subsidence of these buildings foundations are provided. It has been proven that the additional subsidence of existing objects did not exceed the permissible values according to the norms. This happened due to the arrangement of the pit fence, the phasing of soil extraction in the initial stages under of the soil berm protection. The technology provided for the construction of a pit fence with sheet piling (I-beam No. 30) with a step of 1 m, with timber lagging between them. Then they installed the capping beam, struts, rakers and gradually raised the floor and outer wall of the parking lot of the new building.
 The results of long-term geodetic observations of the foundations deformations of buildings located in the zone of the pit influence of the new building at various construction stages of its enclosure are described.
 Sufficient convergence of the simulation results in the plane condition using the FEM and field studies was determined (the relative error did not exceed 20%). An elastic-plastic soil model of the "basement - foundations of an existing building - pit enclosure" system was used for modeling.
 The article describes examples of the results of FEM modeling of soil massif deformations at various stages of pit fence construction.
 It was determined that during construction the largest horizontal movements of the pit fence vary from 8 to 23 mm in the area of existing buildings.The maximum vertical movements of the bases foundations of these buildings amounted to 10 mm, which does not exceed the permissible values according to the norms.
 Prospects for further research on assessing the level of reliability of the pit enclosure by determining the variability parameters of the components of the system "foundation - foundation of the existing building - pit enclosure" have been determined.

Application of continuum damage mechanics in fatigue assessment of A516 steel specimens considering residual stresses

The present study evaluates the effects of tensile residual stresses on fatigue behaviour of ASTM A516 pressure vessel steel specimens. In this regard, fatigue specimens containing a smooth notch region are studied. A damage coupled elastic-plastic constitutive material model is developed to investigate effects of residual stresses on fatigue life. Isotropic and kinematic hardening parameters of the material are experimentally determined. Additionally, the damage model parameters are specified from fatigue experiments conducted on standard specimens in different stress ratios. Tensile residual stresses are introduced into the smooth notched specimens through employing four-point-bending method. The hole drilling technique is utilized to measure the residual stress magnitudes. The results indicate that about two third of fatigue life of the specimens is decreased due to the existed tensile residual stresses. Furthermore, fatigue crack initiation and propagation behaviour in the specimens containing residual stress are studied based on the developed continuum damage model.

Stress evolution in plastically deformed austenitic and ferritic steels determined using angle- and energy-dispersive diffraction

In the presented research, the intergranular elastic interaction and the second-order plastic incompatibility stress in textured ferritic and austenitic steels were investigated by means of diffraction. The lattice strains were measured inside the samples by the multiple reflection method using high energy X-rays diffraction during uniaxial in situ tensile tests. Comparing experiment with various models of intergranular interaction, it was found that the Eshelby-Kröner model correctly approximates the X-ray stress factors (XSFs) for different reflections hkl and scattering vector orientations.The verified XSFs were used to investigate the evolution of the first and second-order stresses in both austenitic and ferritic steels. It was shown that considering only the elastic anisotropy, the non-linearity of sin2ψ plots cannot be explained by crystallographic texture. Therefore, a more advanced method based on elastic-plastic self-consistent modeling (EPSC) is required for the analysis. Using such methodology the non-linearities of cos2φ plots were explained, and the evolutions of the first and second-order stresses were determined. It was found that plastic deformation of about 1–2% can completely exchange the state of second-order plastic incompatibility stresses.

A new strategy for predicting fracture of U‐notched specimens made of Al‐6061‐T6 and Al‐5083 using extended finite element method

AbstractThis research focuses on the numerical investigation of fracture loads of U‐notched specimens made of Al‐6061‐T6 and Al‐5083 under pure opening mode and mixed mode I/II loading conditions. A new methodology is introduced with the aid of the equivalent material concept (EMC) combined with the extended finite element method (XFEM) and linear cohesive zone model (CZM), and the outcomes are compared with a conventional elastic–plastic damage model. The combination of EMC and XFEM result in the numerical EMC‐XFEM model, which provides a fast and accurate method to assess the notch fracture loads and crack growth path. Considering the experimental results, a trial and error loop is performed to calibrate the ductile damage model, which is much more cumbersome than the proposed EMC‐XFEM model. Moreover, the computational cost of elastic–plastic model is much more than the EMC‐XFEM model, although both approaches provide appropriate predictions for fracture of ductile aluminum notches.

Improved inherent strain extraction algorithm based analytical surface modeling for Ti-6Al-4V and SS316L selective laser melting part

Selective laser melting (SLM) is one of the most popular metal additive manufacturing (AM) techniques relied on powder bed and has been widely used in numerous fields. Printing surface modeling is a critical step for product development and quality control of SLM parts as the surface profile has a great impact on the structure and functional performance of printed parts. However, it is still a challenge to predict the printing surface geometry of the SLM part accurately and quickly due to the complex physical processes and thermomechanical behaviors in SLM. To address this barrier, a comprehensive bottom-up modeling approach based on an improved inherent strain method (ISM) is proposed to estimate the surface fluctuation between designed and printed parts and further predict the printing surfaces of SLM parts. The improved inherent strain is extracted from multi-physics analytical calculation considering the thermal influences of in-plane scanning tracks, interlayer re-melting and interlayer strain accumulation. Specifically, temperature filed is obtained from the quasi-static heat analytical model integrating residual temperature involving the heat contribution of in-plane scanning tracks via virtual heat source technique. The interlayer strain accumulation and interlayer re-melting effects are embedded into the inherent strain terms through the layer-lever fluctuation superposition and re-melted depth parameter respectively. Strain compositions are solved through thermo-elastic Green function models and modified elastic-plastic McDowell analytical models, where temperature-dependent material properties are employed. The proposed method has been validated by experimental observations from the literature via two SLM parts using two distinct materials (Ti-6Al-4 V and SS316L).

Fatigue life estimation of stranded copper conductors using the method based on the fatigue damage evolution model

ABSTRACT To consider the effects of plastic damage caused by compaction procedure on the tensile fatigue life of stranded copper conductors, the fatigue life estimation method based on the fatigue damage evolution model was applied to estimate the tensile fatigue life of stranded copper conductors with two different compaction degrees. The macroscopic and mesoscopic elastic-plastic constitutive models of TU2 oxygen-free copper material based on mixed hardening theory and two-scale theory were constructed and the user defined material subroutines (UMAT) for these two models were programmed. The estimation results were compared with fatigue test results and verified the good applicability and accuracy of the applied method on the tensile fatigue life estimation of stranded copper conductors. The effects of the plastic damage caused by the compaction procedure on the tensile fatigue life of stranded copper conductors were obtained.

Effect of L‐T and S‐T orientations on fatigue crack growth in an aluminum 7050‐T7451 plate

AbstractThis study is focused on the anisotropic fatigue crack growth (FCG) behavior of an aluminum AA7050‐T7451 plate. L‐T and S‐T orientations were studied in M(T) samples with W = 50 mm, in mode I loading, with R‐ratio of +0.05. A numerical approach was used, assuming that crack tip plastic strain is the crack driving force. A purely kinematic elastic–plastic model was calibrated using experimental data from low cycle fatigue tests of smooth specimens in L and S orientations. The predicted FCG rates agree well with experimental trends in the Paris' regime, suggesting that cyclic plastic deformation is the main damage mechanism. The numerical model was used to estimate the stress ratio effect for both orientations, which was found to be linked with crack closure variations. However, the closure free predicted trends for both microstructural orientations at R = 0.05 are not overlapped, suggesting an effect of microstructure not linked to crack closure.

Investigation of low cycle fatigue crack propagation behavior of 316H steel at 550℃ based on cyclic response and damage accumulation: experiment and modelling

In the framework of the unified elastic-plastic constitutive theory, the cyclic response and damage evolution of 316H steel in the low cycle fatigue (LCF) regime was analyzed and modeled. A unified damage-coupled elastic-plastic constitutive model was established based on the results of uniaxial strain-controlled LCF tests, which includes a novel strain-based continuum damage model in the framework of the critical plane theory, an accumulative plastic strain and fatigue damage-related kinematic hardening model, and an isotropic hardening model. The comparison between the simulated and the experimental results indicated that the proposed model is efficient to precisely characterize the cyclic response throughout the whole fatigue process. The model was also used to simulate the fatigue crack growth behavior under different maximum loads and load ratios via finite element (FE) simulation. The effects of the load conditions on the crack growth rate and crack tip parameters were analyzed systematically.

Impact mechanical properties of frozen soil after experiencing freeze–thaw cycles: Experimental and numerical study

When frozen soil serves as the foundation of engineering facilities in cold regions, it often suffers from the combined action of freeze–thaw cycles (FTC) and impact loading, which changes its stability and reliability. This study analyzed frozen soil by conducting FTC experiments with varying FTC numbers and impact loading experiments under varying strain rates. The experimental results revealed that frozen soil exhibited FTC and strain rate effects. The peak stress of frozen soil gradually weakened with a decrease in the strain rate and an increase in the FTC number. The expansion stress caused by the water–ice phase transition during an FTC was investigated, assuming that the expansion stress was comparable to the cyclic load, and the deterioration mechanism of frozen soil under the FTC action was discussed. An expression of FTC damage was derived using thermodynamic and continuous damage mechanical principles. A finite element model of frozen soil, derived from the cohesive zone model, was established by inserting cohesive elements into the common surfaces between solid soil elements. A traditional traction–separation (T–S) model for cohesive elements was improved by considering the influence of the strain rate using the overstress model. An elastic–plastic model was used to model the solid elements. The improved T–S model and elastic–plastic constitutive law were used to implement a numerical simulation via a user subroutine VUMAT integrated in ABAQUS/Explicit. When the FTC damage was considered, the microcrack evolution and impact mechanical behavior of frozen soil under impact loading after experiencing FTC action were observed. Having compared the numerical results with the experimental data, it is concluded that a good experimental–numerical correlation was obtained.

Experimental Study and Numerical Calculation of Projectile Dynamic Engraving Process

The high temperature and pressure gas in the combustion chamber pushes the projectile towards the muzzle during the firing of the gun. The projectile belt is squeezed and deformed by the forcing cone. In order to study the dynamic engraving characteristics of the projectile, a short-barrel gun launcher was designed for experimental research. The projectile dynamics and projectile base pressure curves during the engraving process were obtained. Based on the experimental loading conditions, the J-C elastic-plastic model was used to numerically simulate the dynamic engraving process of the copper elastic belt. The deformation process of the elastic belt and the change characteristics of the resistance after the elastic were studied and the intrusion displacement obtained by the numerical calculation was compared with the experimental results and the results are consistent.

Bond-slip model of anchor-soil interface and analysis of bearing capacity of anchor in frozen soil

It has been recognized that temperature and moisture have obvious influence on the bonding characteristics of anchor-frozen soil interface. Following the principle and method used to explore the interface performance of microelement anchor-soil, this research starts by performing the pullout tests of anchor-soil under different temperature and moisture content conditions. Then, we analyze the characteristics of interface bond stress-displacement and compare the effects of temperature and moisture content on interface mechanical properties. The test results show that the failure mode of anchor changes slightly when soil temperature and moisture content vary, wherein the shear failure which shows brittle characteristics at the anchor-frozen soil interface is dominant. Additionally, we have found that the effect of temperature on the peak bond stress is more significant than moisture. Under different soil temperatures and moisture content, the shape of bond stress-displacement curve of anchor-frozen soil interface is similar. Furthermore, the bond stress-displacement relationship of anchor-frozen soil interface can be fitted by a composite function which is composed of exponential and ideal elastic-plastic model function. Thus, this research recommends to propose a bond-slip model of anchor-frozen soil interface considering temperature effect. Combined with the load transfer method, the force and deformation equation of anchor in frozen soil is derived. Through example analysis, it is discovered that the pullout bearing capacity of anchor increases with the decrease of soil temperature. Moreover, the initial shear stiffness of anchor-soil interface is positively related to the bearing capacity of anchor. At the end, the influence of interfacial residual strength on the bearing capacity of anchor in the later stage of pullout is superior to that in the earlier stage.