Finite Element Simulation of NiTiNb Shape Memory Alloy Pipe-Joint Subjected to Coupled Transformation and Plastic Deformation

*e assembling process of Ni47Ti44Nb9 alloy pipe joints considering the phase transformation and plasticity was numerically simulated for the first time with a developed constitutive model. *e simulated process was based on the experimental material parameters, which were determined with the experimental tensile results of Ni47Ti44Nb9 shape memory alloy (SMA) and steel bars.*e results showed that, after assembly, theMises stress distributed uniformly along the longitudinal direction of the NiTiNb joint, but nonuniformly along the radial direction. *e maximum σeq does not appear at the inner wall of the joints due to the coupling effect of the plastic deformation and the recoverable transformation. *e contact pressure distributed uniformly along the circumferential direction, but nonuniformly along the longitudinal direction. *e sizes of the SMA joint and the pipe should be properly matched to ensure contact during the stage of the rapid reverse phase transformation to obtain stable connection performance. *e pull-out force was also computed, and the results were in good agreement with the experimental results. *e results obtained can provide available information for the optimization of the design parameters of the high-performance SMA pipe-joint, such as inner diameter and assembly clearance.


Introduction
Due to their two unique properties, the shape memory effect (SME), and pseudoelasticity (PE), the shape memory alloy (SMA) has drawn significant attention in a broad range of commercial applications, such as automotive, aerospace, robotics, and biomedical domains [1,2]. e pipe-joint is manufactured with SMA even widely used as the oil pipe connection in aircraft [3]. Among the existing SMAs, NiTi binary alloy exhibits excellent mechanical properties and good SME (with about 8% recoverable strain). While the NiTi SMA pipe-joint should be stored and transported in liquid nitrogen to avoid the unexpected recovery of the predeformation, in contrast, NiTiNb SMA exhibits a much larger transformation temperature hysteresis (∼150°C) compared with conventional NiTi binary alloy (∼30°C) [4], and the pipe joints made of NiTiNb alloys can be stored and transported at room temperature, which is greatly convenient for engineering applications.
To date, the literature referring to experiments and numerical simulations of NiTiNb SMA is quite limited. Zhang et al. [5] analyzed the constitutional phases of Ni 47 Ti 44 Nb 9 with a scanning electron microscope and found that it is composed of three phases: TiNi matrix, β-Nb particles, which are softer than the TiNi matrix and can easily be deformed plastically, and compound Ti3(Ni, Nb)2.
increase of Ms is attributed to the micromechanism of Nb solution into NiTi alloy [6]. Jiang et al. [7] also investigated the influence of Nb addition on phase transformation, microstructure, and mechanical properties of equiatomic NiTi SMA. ey found the Nb addition has a slight influence on the transformation hysteresis. It can be attributed to that the NiTiNb SMA is not subjected to plastic deformation and the β-Nb precipitate phase is unable to relax the elastic strain energy in martensitic interfaces. Chen et al. investigated the thermomechanical behavior of Ni 47 Ti 44 Nb 9 (at.%) SMA under different thermomechanical loading histories. eir experimental results show that the equivalent stress-strain behavior and the activated martensite variants depend strongly on the thermomechanical loading history [8]. Jiang et al. found that severe plastic deformation based on equal channel angular pressing plays an important role in enhancing the transformation hysteresis of Ni 47 Ti 44 Nb 9 (at.%) SMA [9]. e effects of prestrain and heat treatment on shape recovery stress and phase transformation temperature were experimentally investigated by Uchida et al. [10], wherein NiTiNb SMAs with various Nb content (6, 9, and 12 mol%). e recovery stresses of NiTiNb alloys increase with the increase of prestrain and decrease after reaching the maximum recovery stress of around 450 to 500 MPa at about 9% prestrain. e effect of low Nb content (3.5 at.%, 4.5 at.%, 5.0 at.%) on the mechanical properties of NiTiNb SMA has been investigated in detail by differential scanning calorimetry (DSC), tensile test, and scanning electron microscope (SEM) [11]. e experimental results show that with the increase of Nb content, the yield strength increases, and the elongation decreases gradually.
Recently, the metal joint industry is developing rapidly, including the product method [12], the combination approach [13], and new processing technology [14]. While the literature on SMA joints is relatively few, Tabesh et al. [15] designed a commercial NiTiNb coupler and tested it. ey found that the thermal-induced contact pressure depends on several factors such as the dimensions and properties of the coupler-pipe system. ey also developed coupler-ring finite element model with the constitutive equations established by Machado and Lagoudas [16]. is 3-dimensional phenomenological constitutive model focuses on the generation and recovery of transformation strains that occur because of martensitic transformation. e Gibbs energy is selected to be the thermodynamic potential in the constitutive model. e comparison and evaluation of the shape memory effect under different expansion temperatures and expansion strains were performed. It is reported that the largest strain for the NiTiNb SMA joints appears on the inner wall, and there exists a larger strain gradient along the axis [17]. Piotrowski and his coworker did excellent studies on the constitutive model of NiTiNb SMA and the application to tightening rings [18,19]. In their research studies, the Nb precipitates were considered as elastic-plastic inclusions embedded in the NiTi matrix. Heating can induce inverse transformation, and the SMA rings will be tightened automatically on the pipes. e variation of the contact pressure against temperature was measured and compared with the numerical results [18]. In addition, the effects of the ridge thickness, the distance between ridge and ring extremity, and the geometrical defects were simulated [19]. However, the constitutive models adopted in the above excellent work did not consider the plasticity of SMA and its effects, while during predeformation, plastic deformation involved in the NiTi matrix can obviously increase the transformation hysteresis temperature [20]. e plastic deformation plays an important role in thermomechanical behavior during the following assembly process and should be taken into account in the analysis. Great efforts were made in the constitutive model of SMAs [21][22][23]; however, due to the complex microstructure of NiTiNb SMA, few models can accurately describe its specific behaviors. Chen et al. [24] developed a three-dimensional phenomenological model, which is coupled with the transformation and plastic behavior of SMA. In the model, the strain is separated into elastic, thermal, transformational, and plastic strains, and two yield functions are adopted to describe the transformation and plastic deformation, respectively. e model can well describe the behavior of the SMAs subjected to complex thermomechanical loading paths, including shape memory effect, pseudoelasticity, coupled transformation and plastic deformation, and the effect of plastic deformation on the inverse transformation.
In the present study, the assembling process of Ni 47 Ti 44 Nb 9 alloy pipe joints is simulated with the threedimensional SMA model developed by Chen et al. [24]. Two types of the pipe-joint assembly processes are considered: (1) the same steel pipes but the SMA joints are of different outer diameters (the inner diameters of the SMA joints are the same); (2) the same SMA joints but the steel pipes are of different outer diameters (the inner diameters of the steel pipes are the same). e main content of this article contains two parts. First, the material parameters of NiTiNb joints and steel pipes are determined with an experiment, based on which the finite element (FE) model is built. Second, the distribution of Mises stress, contact pressure, and other thermal-mechanical quantities in the pipe-joint systems is analyzed and discussed.

Experiment and Finite Element Model
e transformation temperatures of the Ni 47 Ti 44 Nb 9 alloy used in the present work were obtained by DSC [8], as shown in Figure 1. e initial four transformation temperatures are determined from the curves, as summarized in Table 1. It was reported that there is a characteristic deformation temperature (Ms + 30°C) and a maximum tensile strain (about 16%), at which the transformation hysteresis can effectively increase and the strain recovery ratio remains significantly high [25]. In general, NiTiNb joints are expanded at Ms + 30°C before assembly, and the circumferential tensile stress is the main stress formed during this process. e corresponding tensile strain can reach ε � 16%. In this study, Ni 47 Ti 44 Nb 9 alloy was used as the pipe-joint material, and the stress-strain curve (with ε max � 0.16) during tensile deformation at −60°C and subsequent heating recovery were used to determine material parameters.
Material composition and the size of test samples are similar to the literature [4], as shown in Figures 2(a) and 2(b). e uniaxial tensile tests of Ni 47 Ti 44 Nb 9 alloy under −60°C were conducted on the MTS-858 electrohydraulic servo system, with axial loading capacity ±15 kN and torque ±100 N-m (see Figure 2(c)). e stain-control mode was adopted and the loading/unloading rates were chosen as ±1.3 × 10 4 s −1 . As shown in Figure 2(c), the Epsilon (3542-010M-025-LHT) extensometer, with an accuracy grade of 0.5 (ISO9513) and temperature range of −270°C to 200°C, was used to measure the tensile strain. e red box is the installation diagram. A GDW-50 temperature chamber was used to provide the testing temperature ranging from −70°C to 150°C, with fluctuation of ±0.5°C and nonuniformity of ±2°C (see Figure 2(d)). e heating/ cooling rates were controlled ±2°C/min to avoid the inhomogeneous distribution of temperature within the testing specimen. e stress-strain (σ 33 − ε 33 ) curve of loading and unloading is shown in Figure 3. e tests were repeated three times with virtually the same results. is repeating practice was emphasized by Krahmer et al. [26] and was used to check the stability of material response. A typical "elasticitytransformation-strengthen" loading curve can be found and during unloading, no obvious recoverable strain occurred. e following constraint heating result is shown in Figure 4(a) and the stress-free heating result is shown in Figure 4(b). For the case of Figure 4(a), during the heating process, the ε 33 was kept unchanged by the MTS machine. One can find that the stress σ 33 decreases with the increase of temperature (T) from −60°C to about 50°C; this is due to thermal expansion of the alloy. Heating leads to a tendency of elongation along the sample's axial direction, while ε 33 was kept unchanged by the MTS machine at this time. us, the compressure occurred and σ 33 decreased. When the temperature was heated to 50°C, the stress σ 33 increased with the increase of temperature (T). It  ] Advances in Materials Science and Engineering 3 can be attributed to the reverse phase transformation that occurred in the heating process and the sample tended to shorten. Similarly, by the constraint of collets, stretching stress appeared and σ 33 increased. For the case of Figure 4(b), it can be found that the residual strain decreases slowly before the temperature is around 50°C. Between 50°C to 75°C, the residual strain decreases rapidly, which is attributed to the reverse phase transformation. When the temperature exceeds 75°C, the rate of decrease goes back to the initial level. Comparing with Figure 1, one can find that the predeformation increases the As from −68.29°C to about 50°C. e mechanical property of the steel pipe to be connected was also tested at room temperature, as shown in Figure 5. In the initial stage of uniaxial tension, the steel exhibited elastic deformation, while the stress raised to about 450 MPa; obvious plastic deformation can be found. e steel exhibited an elastoplastic feature of linear strain-hardening.

NiTiNb Constitutive Model.
e constitutive model developed by Chen et al. [24] was used, in which the effects of plastic deformation were considered. As the effects of plastic deformation of Nb phase and martensite on the reverse phase transformation temperature are similar, we make no distinction between the two plasticities and make them equivalent to a unified plastic evolution equation. e brief introduction of the model is described as follows. e model assumes that the free energy V is the function of elastic strain ε e , transformation strain ε t , tensor α and accumulated plastic strain p associated with plasticity-induced hardening, and temperature T: where where ε and ε p are the total strain and plastic strain. ermal strain ε θ can be expressed as follows: in which θ is the coefficient tensor of thermal expansion and T 0 is a reference temperature. Considering the interaction of transformation and plasticity, V can be further expressed as follows: where Φ e ε e , T � 1 2ρ where C is the standard isotropic elasticity tensor, C p and C t the are material parameters, and ρ and C v are the density and specific heat, respectively, which are assumed invariant during phase transformation and plastic region. τ M (T) can be defined as follows: where B represents the sensitivity of the stress with respect to temperature T and M 0 is the reference temperature below Advances in Materials Science and Engineering which no twinned martensite is observed [27], and it is usually regarded as the martensite transformation finished temperature M f . ‖ε t ‖ denotes the Euclidean norm of ε t of rank two. e function Γ(ε t ) is similar to the additional indicator function adopted by Frémond and Nedjar [28], which satisfies the following phenomenological physical constraints: in which c is the Lagrange multiplier and ε L is the critical maximal value of phase transformation under saturated state; it is a given constant in the process of calculation. e effects of plastic deformation on phase transformation is considered by Q(p), and it can be defined as follows: where Q 0 and ξ are the material parameters. e function Q(p) is based on the fact that there is a close relationship between the inverse transformation start temperature As and the permanent strain [20]. Such a function is used to describe the phenomenon that the influence of plastic deformation on the phase transformation will reach a certain limit. e drag stress of plastic is defined as follows: in which R 0 , ζ, and Y 0 are the material parameters. e above constitutive equations can be rewritten in incremental ones and calculated via the return mapping scheme (an implicit algorithm). e corresponding user material subroutine UMAT was developed and embedded into FE code ABAQUS, with which the installation process of NiTiNb pipe joints can be simulated. A detailed model and algorithm can be found in [24].

Parameters Identification.
Based on the experimental results, the parameters of the NiTiNb SMA constitutive model can be identified. e four phase transition temperatures are determined from the DSC curves. e elasticity modulus E, the starting stress of martensite transformation (σ t s ) and plasticity (σ p s ), and the maximum strain of martensite transformation (ε L ) can be determined from the uniaxial tension. e B, R 0 , and Q 0 can be determined from the constrained heating test. e relative parameters of the Ni 47 Ti 44 Nb 9 alloy and the steel are summarized in Table 1.

Finite Element Model.
e commercially available finite element code ABAQUS ® is used for simulation. e UMAT subroutine based on the constitutive model is developed and embedded in ABAQUS. A kind of aeronautics oil pipe joints is simulated, which is 3 mm in inner diameter and 11 mm in length. As the contact force mainly affects the contact area between the joint and pipe, to reduce computation time, we take the total length of two pipes as 13 mm. To compare the effect of the thickness and assembly interference on the mechanical property during assembling process, two groups of pipe-joint systems were simulated: (1) steel pipes of same size and the NiTiNb joints with the same inner diameter but different external diameters; (2) NiTiNb joints of the same size and steel pipes with the same inner diameter but different external diameters. e specific sizes of the pipe-joint and steel are listed in Table 2.
In general, we expand the NiTiNb joint to ε � 16% at T � −60°C, heat the temperature to room temperature, and then assemble the pipes and joint, followed by heating the joint to a prescribed high temperature (e.g., 100°C). Whole assembling process can easily be resolved into the five steps, as shown in Table 3.
Considering the symmetry of the system, we only calculated a quarter of the system using an axisymmetric model, as shown in Table 3. In the FE model of the pipe-joint system, 8-node linear brick elements were adopted to mesh the NiTiNb joint, and 8-node linear brick reduced integration elements were adopted to mesh the steel pipe. Since contact behavior between the outer surface of the steel pipe and the inner surface of the NiTiNb joint is surface contact forms during the assembly and drawing process, the "penalty" contact pair form of tangential behavior is added. As the elastic modulus of the steel pipe is greater than that of the NiTiNb joint, the outer surface of the steel pipe is defined as the active contact surface, while the inner surface of the NiTiNb joint is defined as the passive contact surface, with a friction coefficient of 0.2. To avoid the nonconvergence caused by rapid strain recovery in the NiTiNb joint during heating, the maximum time increment is set as 0.001 (the step duration is 1.0).

Model Checking.
e uniaxial tension at −60°C and subsequent constrained heating recovery are simulated firstly to verify the correctness of the constitutive model. e results of the calculation are compared with experimental results, as shown in Figures 6 and 7. e tensile stress-strain (σ 33 − ε 33 ) curve calculated with the constitutive model agrees reasonably with the experimental results. However, since the phenomenological model in this study does not consider the local instability and softening of the material caused by the martensite nucleus, the initial stress drop of the phase transformation cannot be described. For the plastic deformation stage, the von Mises stress of the finite element simulation results is slightly larger than the experimental results. As indicated above, the yield function defined in the present work is a convex function. erefore, the numerical simulation results after entering the plastic deformation stage are convex compared with the experimental results. e measured residual strain after unloading is less than that of numerical simulation, because the phenomenological model does not consider the local evolution characteristics of material properties during unloading, such as the local reverse phase transformation and reorientation, so it shows stronger linear elasticity compared with the experimental results. However, for the engineering application of the SMA pipejoint, the above deficiencies have less impact on its service performance. e comparison between the experimental and numerical results shows that the model can well describe the main characteristics of the deformation in tension and the subsequent thermal recovery property. It indicates that the model can reasonably be used to simulate the behavior of NiTiNb SMA joints.

Distribution of Mises Stress.
After assembling, the contour maps of Mises equivalent stress, σ eq , in the joints J1, J2, and J3 connected with P3 (Table 2) are shown in Figure 8(a). It can be seen that the σ eq beyond the contact area of the pipe is very small, which implies that the effect of the assembly contact along the axis of pipe can be ignored out of the contact area. e distribution of σ eq in the contact region is quite different among the three cases, and the maximum value of σ eq decreases with the decrease of the outer diameter of the joint. Within the joint, σ eq distributes  Step I

Low temperature
Expanding inner diameter of the joint from 3.0 mm to 3.48 mm (corresponding to the circumferential strain of about 0.16 at its inner wall) at −60°C.

NiTiNb joint
Steel pipe Steel pipe A Y Z Step II

Room temperature
Unloading the expanding pressure, and raise the temperature to 23°C Step III

High temperature
Insert the pipes into the joint, followed by heating the pipe-joint system to 100°C.
Step IV

Room temperature
Cooling the system down to room temperature (23°C). After assembly, add the following step to calculate the pull-out force.
Step V

Room temperature
Pulling out the pipes from the joint. Advances in Materials Science and Engineering uniformly along the longitudinal direction, but nonuniformly along the radial direction. To facilitate the analysis, the contour map of σ eq in the cross-section of the NiTiNb joint is shown in Figure 8(b) and the variation of σ eq against the radius is shown in Figure 9. It can be found that for case J2-P3 and case J3-P3, σ eq increases monotonically from the inner wall to the outer wall, but the changes become insignificant after λR (the radial position) is 2.3 mm. us we do not need to increase the thickness of the NiTiNb SMA joint wall unlimited. While for case J1-P3, σ eq increases firstly until λR � 2.16 mm, and then decreases. It indicates that the distribution of σ eq is closely related to the wall thickness of the joint. An obvious difference is that the maximum σ eq does not appear at the inner wall of the joints, which is different from the previous study [18,19]. is interesting phenomenon is closely related to the coupling effect of the plastic deformation and the recoverable transformation of the NiTiNb joint. Interfered by the external environment, the joint whose maximal equivalent stress distributed on the surface outside or inside will be easier to fail. erefore, appropriate wall thickness and the ratio of expanding diameter should be designed. Figure 10 shows the distribution of σ eq along the radial direction in the pipe. It can be seen that σ eq decreases monotonically with the increase of λR (from inner to outer). Taking J1-P3 as an example, σ eq decreases from 420 MPa at the inner surface to about 300 MPa at the outer surface. On the other hand, with the increase of wall thickness, σ eq in the pipe increases from about 220 MPa (for J3-P3) to 420 MPa (for J1-P3). erefore, when selecting the wall thickness of the NiTiNb SMA joint, the bearing capacity of the steel pipe should be considered to avoid the overlarge restoring force in the assembling process. e von Mises stress, plastic deformation, and phase transformation strain distributed in the NiTiNb SMA joint at the end of different steps are analyzed. For instance, the results of the J1-P3 system are shown in Figure 11. It can be found that, in step I, the Mises stress decreases from inside to outside; the same law can also be found in plastic strain and transformation strain. In the constitutive model of NiTiNb SMA, when the stress is greater than 220 MPa, the material gets into the plastic deformation stage (see Figure 6). erefore, the plastic strain is very small, and the phase transformation strain has not yet reached the critical value in the outer wall of the joint. In step II, the load is removed at −60°C, so the von Mises stress within the joint is very small, while the plastic strain and the phase transformation strain remain unchanged. It is reasonable because the plastic strain is irreversible and the reverse phase transition temperature is raised after the predeformation, so there is no transformation strain to recover. After heating, the transformation strain is partial recovery because the constraint from the steel pipe is existing on the joint. us the transformation strain cannot fully recover. It can also be found that the plastic strain keeps unchanged in the subsequent steps since it comes into being in step I. e above analysis of different variables in the SMA joint coincides with the response of NiTiNb material and deepens the understanding of the deformation characteristics of the SMA joint during the whole assembly process.

Contact Pressure Distribution.
e distribution of the contact pressure (σ r ) at the pipe-joint interface is shown in Figure 12. It can be found that the contact pressure distributes uniformly along the circumferential direction. However, it distributes nonuniformly along the longitudinal direction and symmetrically along the middle point. e inhomogeneous stress gradient along the longitudinal direction will bring an inhomogeneous strain on the contact surface, which may increase the surface friction between the pipe and joint. Besides, with the increase of joint thickness, the maximum contact pressure increases, and the stress concentration attenuates. To better understand the result, the longitudinal distribution of σ r within J1, J2, and J3 is shown in Figure 13. For J3, the maximum σ r appears at the  Advances in Materials Science and Engineering end of the contact region, and the minimum σ r appears at 1 mm from the endpoint. For J2, the maximum σ r appears at the endpoint and central area. Such a pressure distribution is advantageous to the seal but is adverse to the effective work of the NiTiNb SMA joint, as the larger the contact pressure at the edge, the easier the destruction. For J1, maximal contact pressure appears at point 2 mm from the endpoint and minimal pressure appears at the endpoint of the contact region. is distribution can not only ensure a good seal but also benefit the effective work of the NiTiNb SMA joint. e thicker the joint wall, the larger the contact pressure. A thicker NiTiNb SMA joint would induce an enclasp force for a better seal effect, which is beneficial for the application of NiTiNb SMA joints.
As contact pressure is one of the key indicators of the pipe-joint systems, the change of contact pressure during the     Figure 9: Distribution of σ eq along the radial direction in the joint at room temperature.

Advances in Materials Science and Engineering
whole assembling process should be concerned. Figure 14 shows the variations of σ r at the pipe-joint interface (point A in Table 3) in the cases J1-P3, J2-P3, and J3-P3 (see Table 2) against T during the assembling process. ere is a similar tendency in the three kinds of joints: with the increase of T, σ r remains zero until T increases to a certain value. is is because the SMA joint is not in contact with the steel pipe. en σ r increases linearly with the increase of T up to T �100°C, which is due to the combined action of thermal explanation and reverse phase transformation. In the following cooling process, T reduces from 100°C to 23°C, a low decrease of σ r can be found. It is attributed to the different   coefficients of thermal expansion of steel pipe and NiTiNb SMA joint. However, the following two main differences can also be found: (1) the temperatures for the onset of inverse transformation are different; i.e., they are about 22°C, 36°C, and 53°C for J1-P3, J2-P3, and J3-P3, respectively; (2) the contact pressures after cooling to room temperature are different, that are about 143 MPa, 126 MPa, and 80 MPa for J1-P3, J2-P3, and J3-P3, respectively. e reason for the differences is that for the same expansion deformation, the thinner the tube wall of the joint, the larger the plastic deformation in the joint. Plastic deformation that occurs during a prestrain process would relax the free energy stored in the material, leading to a higher temperature for the inverse transformation [20]. erefore, the temperature for the onset of inverse transformation increases with the  Figure 12: Contours map of pipe-joint interfacial pressures in J1, J2, and J3 at room temperature.  decrease of the joint thickness. Furthermore, the deformation in the radial direction is nonuniform; the plastic deformation near the inner wall is larger than that near the outer wall, as shown in Figure 11. During heating recovery, the plastic deformation region near the inner wall is subjected to pressure, which may counteract a part of the recovery force. e larger the plastic deformation is, the larger the recovery force would be counteracted, leading to the reduction of the final pressure from J1 to J3. Figure 15 shows the variations of σ r at point A (see the FE model embed in Table 3) against T during heating and cooling processes, in the case of P1 to P5 connected with J1 (see in Table 2). It can be seen that the curve of P1 is remarkably different from that of the others. e contact pressure of P1 is far less than the other four pipes. is difference is very important because the contact pressure directly determines the connection performance of the joint-pipe system. To explain this, the variation of the displacement at point A during stressfree heating after Step I is shown in Figure 16. e heating recovery curve agrees well with the experimental result ( Figure 4(b)). In the initial stage, the displacement recovers slowly. e displacement recovers rapidly at near 20°C and then returns to the initial level. It is interesting to note that the reverse transformation temperatures are about 20°C, which are below the 50°C in Figure 4(b). It can be attributed to the fact that the predeformation in the SMA joint is distributed in gradient while in the test sample, it is distributed uniformly. In Figure 16, the red arrows correspond to the positions of the outer wall of P1 to P5. It can be seen that for P2 to P5, the SMA joint contacts with the steel pipe during the rapid recovery segment, which may account for why the curves of P2 to P5 rise at near 20°C, as shown in Figure 15. Instead, the SMA joint contact with P1 until the rapid reverse phase transformation is over, and the temperature rises to about 50°C. It indicates that the main part of the heating should also occur at T � 20°C, but before T � 50°C the recovery deformation is not sufficient to contact with the pipe. us, no σ r appears. As the joint contacts with the pipe at about T � 50°C, σ r begins to increase. While the main part of the recoverable strain has already occurred, it accounts for the smaller σ r after recovery. At this point, the difference in Figure 15 is explained. From the result, it can be concluded that in order to obtain stable connection performance, the sizes of the SMA joint and the pipe shall be properly matched to ensure contact during the stage of the rapid reverse phase transformation.

Pull-Out Force.
e pull-out process after assembling is performed at room temperature and the variation of the pull-out force (F) against the pull-out displacement (U) is shown in Figure 17 for J1-P3 system. In the beginning, the frictional resistance increases rapidly, mainly attributed to the static friction generated by contact pressure. e peak is F ≈ 1200 N, which is close to the experimental results (between 1300 N and 2200 N). e difference may result from the existence of the coating in the actual structure, which can not only improve leak tightness of the joint system but also enhances the joint strength by slowing down the stress concentration of the joint [29] and increases the frictional resistance. After the peak, F decreases almost linearly with the increase of U, attributed to the decrease of the remaining contact area.
From the above results and analysis, one can find that the stress is inhomogeneous, whether in the joint, pipe, or on their contact surface.
us employing traditional plane strain assumption or assuming the stress in joint and pipe as constant is inaccurate. A model that utilizes the FE method can compensate for the drawback. It can be obtained that the pull-out force increases with an increase of the NiTiNb SMA joint thickness; at the same time, the sealing property is also improved. erefore, in the design process, some means such as reducing the joint diameter, increasing the wall thickness, and reducing the assembly clearance can be used to improve the contact pressure.

Conclusions
e assembling process of steel pipe-NiTiNb alloy joint systems was simulated with the FE approach, making use of a three-dimensional NiTiNb SMA model that takes into account coupled transformation and plastic deformation. e simulation includes the expansion of the cylindrical cavity at a low-temperature, heating and cooling of the pipe-joint system after assembling, and the following pullout process. Some conclusions can be drawn and listed as follows.
(1) After assembly, the Mises stress within the NiTiNb joint was closely related to the wall thickness of the joint. e maximum Mises equivalent stress did not appear at the inner wall of the joints, and the distribution was closely related to the coupling effects of plastic deformation and the recoverable transformation of the NiTiNb joint. (2) e contact pressure uniformly distributed along ring direction but presented inhomogeneous distribution along the axial direction. e thicker the NiTiNb SMA joint wall, the larger the contact pressure. A thick NiTiNb SMA joint would provide a better seal effect and a larger enclasp force than a thin one. (3) During the heating recovery, the contact pressure began to increase when a certain temperature was reached; then the contact pressure increased linearly and got its maximizing at 100°C. In order to obtain stable connection performance, the sizes of the SMA joint and the pipe shall be properly matched to ensure contact during the stage of the rapid reverse phase transformation. e pull-out force reached its maximizing at the beginning and decreased linearly with the pull displacement.
Data Availability e numerical data and pictures used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest. Advances in Materials Science and Engineering 13