A Unified Empirical Equation for Determining the Mechanical Properties of Porous NiTi Alloy: From Nanoporosity to Microporosity

The mechanical characteristics of a monolithic (non-porous) crystalline or amorphous material are described by a well-defined set of quantities. It is possible to change the mechanical properties by introducing porosity into this material: as a rule, the strength values decrease with the introduction of porosity. Thus, porosity can be considered as an additional degree of freedom that can be used to influence the hardness, strength and plasticity of the material. In the present work, using porous crystalline NiTi as an example, it is shown that the mechanical characteristics such as the Young's modulus, the yield strength, the ultimate tensile strength, etc. demonstrate a pronounced dependence on the average linear size $\bar{l}$ of the pores. For the first time, an empirical equation is proposed that correctly reproduces the dependence of the mechanical characteristics on the porosity $\phi$ and on the average linear size $\bar{l}$ of the pores in a wide range of sizes: from nano-sized pores to pores of a few hundred microns in size. This equation correctly takes into account the limit case corresponding to the monolithic material. The obtained results can be used directly to solve applied problems associated with the design of materials with the necessary combination of physical and mechanical characteristics, in particular, porous metallic biomaterials.


Introduction
Porous materials have a unique combination of physical and mechanical properties that is not found in their monolithic analogues [1,2,3,4,5].For example, the presence of percolating voids and a significant specific surface area typical of porous materials provide their lower density, lower thermal conductivity and high catalytic activity compared to non-porous analogues [6].The potential functional applicability of porous materials is largely determined by their mechanical properties M [7,8,9,10,11].Here the quantity M denotes a set of properties such as the Young's modulus E, the yield stress σ y , the shear modulus K s , the bulk modulus G b , etc., i.e.M = {E, σ y , K s , G b , . . .}. ( It is well known that the mechanical properties M of the monolithic crystalline or amorphous material are mainly determined by thermodynamic conditions: the pressure p and the temperature T [12,13,14,15]. Let us assume that there is a bulk porous material whose characteristic dimensions -length, width and height -are comparable.In this case, the linear size of the material can be estimated as 3 √ V , where V is the volume of the material.Obviously, such specific porous materials as porous rods, membranes, porous thin films are not included in this category.The porosity of the bulk material can be quite unambiguously characterized by two quantities: the pore size distribution function P (l) and the interpore partition size distribution function P (g).Other characteristics of the porous structure -such as the most probably pore shape, the type of pores, the average pore surface area, etc. -should be directly or implicitly related to P (l) and P (g).On the other hand, it is convenient to characterize the porous structure in terms of scalar quantities.Therefore, the functions P (l) and P (g) can be replaced by the pair of quantities with some loss of uniqueness characterization: l = l P (l) dl P (l) dl and ḡ = g P (g) dg P (g) dg , where l is the average linear size of the pores, ḡ is the average linear size of the interpore partitions.
In practice, both l and ḡ can be estimated by statistical analysis of "patterns" (pore contours) on thin sections of the material.In addition, these two quantities correlate with a well-known quantity such as the total porosity of the material ϕ = V 0 /V , where V 0 is the volume of empty space in the material, the total volume of which is V .The porosity ϕ is dimensionless and ϕ ∈ (0; 1); for ϕ = 0 we have the case of the monolithic material.It is therefore clear that it is convenient to use the porosity ϕ instead of one of the two quantities l and ḡ to characterize a porous material [16,17,18,19].
In fact, one has two physical quantities: (i) the average linear size of a pore and (ii) the average linear pore size for a whole system.First, the average linear size of a pore is defined as the average of all possible distances between opposite walls of a considered single pore.The distance between opposite walls means the length of the segment that is common to the perpendiculars drawn from the opposite walls.Second, the average linear pore size l for a whole system is just an average characteristic for all the pores of the material.It can be simply defined as the normalized first moment of the pore size distribution function P ( l) [see Eq. ( 2)].In the general case, we have l ∈ (0; Here, l0 is the quantity characterising the average interatomic distance in the monolithic material. The value of this quantity can approximate the size of point defects.If we define the average size of the atoms forming the alloy as d, then the value of l0 can be estimated as l0 ≃ (1/3)d.For l = l0 we have the monolithic material; and for l → 3 √ V there is no material as such.The range of possible values of l is limited by the value of the porosity ϕ according to the following inequality: On the ( l, ϕ)-diagram, a porous material can only be realized within a strictly defined range of values of these quantities, as shown in Figure 1.
In the case of a porous material with non-zero values of l and ϕ, an arbitrary mechanical characteristic from the set M [see relation (1)] will be a function of l and ϕ, i.e.
Almost all strength properties of the material decrease with increasing l and ϕ [20,21,22,23].In general, the following conditions can be formulated which characterize the decreasing strength properties with increasing porosity and/or pore size.
The condition (5a) is taken into account by almost all known models.So, the example of such a model is the Bert's power law equation for porous materials with elliptically shaped pores [24,25,26]: where K = 0.75 + 1.25(b/c).Here, the parameter c is the average linear size of the pores along the applied stress; b is the average linear size of the pores perpendicular to the direction of the applied stress; M 0 is the mechanical properties of the monolithic material; ϕ m is the maximum possible porosity for the material.Equation ( 6) is a modification of the linear Rossi expression [27], and this expression is mainly used to describe the dependence of elastic properties on porosity.In Eq. ( 6), the degree function is required for a more accurate description of the mechanical properties at porosity ϕ > 0.2.In addition, in the case of an isotropic porous system with isolated pores whose sizes have a random distribution, the dependence of mechanical properties on porosity can be described by Mackenzie's semi-empirical equation [28]: Here, the parameters d and g characterize the shape of the closed pores.At the same time, Eq. (7) incorrectly describes the quantity M(ϕ) in the case of a porous system with open pores, due to the difficulties associated with the determination of the linear dimensions of the pores.There are empirical model that do not explicitly take into account the morphology of a porous system.The parameters in such models are fitted through approximation of empirical data for a sample with a specific synthesis protocol, which makes the application of such models more universal.An example of such a model is the Ryshkewitch-Duckwoth exponential equation [29,30]: where the exponent B is the empirical constant.It is important to note that this equation correctly describes the behaviour of mechanical properties at porosity values up to 0.6.Bal'shin power law equation is suitable for reproducing the mechanical properties of porous materials over a wide range of porosities [31]: where the exponent n is the parameter that depends on the type of material and on the deformation protocol of the material.Thus, all the above equations take into account the fact that the introduction of porosity into the material leads to a significant decrease in the elastic, plastic and strength properties according to the exponential or power law.
Let the porosity ϕ be fixed.The average linear size l of the pores can vary over a fairly wide range, and the maximum possible pore size lmax cannot exceed the size of the material.For such a fixed porosity ϕ, the maximum possible pore size lmax can be realized in a material with a single large pore; and a porous material with the minimum possible pore size lmin is manifested by the presence of a certain number of pores with an extremely small non-zero size, i.e. lmin ̸ = 0.The aim of the present work is to find out how the mechanical strength properties of the porous material depend on the pore size.It is assumed that a fixed porosity (ϕ =const) can be provided by both millimeter-sized pores and extremely small pores of the order of tens of nanometres.This problem is solved using the example of such a well-known metal alloy as NiTi under normal conditions (temperature T = 300 K and pressure p = 1 atm).
The NiTi alloy is considered as a promising material for use in engineering and medicine due to its unique combination of physical properties, including shape memory effect, superelasticity, high corrosion resistance and biocompatibility [18,32,33].Materials based on this alloy have found wide application as actuators in microelectromechanical devices [34], in the aerospace industry [35,36], in the manufacture of medical implants [37,38].Porous NiTi has a large specific surface area, which allows it to be used to treat or replace damaged human organs due to its ability to penetrate and implant biological tissue into the pore space, where this tissue develops naturally [39].At the same time, porous NiTi, like most other porous materials, has a lower resistance to deformation than its monolithic analogue [16].It should be noted that the introduction of porosity into NiTi alloy can change its microstructure limiting phase transitions associated with austenitic-martensitic transitions [33].The presence of defects and free surface due to porosity as well as the formation of a percolating branched structure by the crystalline matrix can limit the superelasticity and shape memory effect [40,41].In the present work, we do not study the effect of porosity on the austenitic-martensitic phase transitions because the porosity of the considered porous NiTi samples is fixed.
The NiTi alloy is an intermetallic compound whose crystal structure has been identified as the high-temperature B2 cubic phase (austenite).The NiTi alloy can undergo a martensitic transition to the low-temperature orthorhombic B19' structure (martensite) with monoclinic distortion through an intermediate rhombohedral R-phase.This transition can be induced by heat treatment and forces such as tension, compression and bending.The structural transition occurs on cooling through the scheme B2→R→B19', while the reverse martensitic transition B19'→R→B2 is realized on heating [39].Due to the austenitic-martensitic phase transitions, NiTi alloy has a shape memory effect and superelasticity -the ability to fully recover its original shape after removal of an external load, which can significantly exceed the yield strength of the material [42,43].
At normal atmospheric pressure, NiTi has the relatively high liquidus temperature T l ≃ 1580 K, which allows this alloy to be classified as heat resistant [44].Therefore, NiTi is widely used in the aerospace, automotive, electronics and biomedical industries, for example, in the manufacture of self-sealing thermomechanical joints, thermally sensitive mechanical actuators, robotic elements, electronic device actuators, instruments and implants for cranio-cerebral and maxillofacial surgery, traumatology, orthopaedics, etc. [4,45].

Materials and Methods
To  (I) Monolithic material.This case can be considered as the limit case of a porous material, where the average linear pore size l approximates the linear sizes of the point defects, i.e. l ≃ l0 ≃ 0.09 nm.
The mechanical properties of monolithic NiTi depend on the type of crystal lattice.For example, for B2 phase monolithic NiTi the Young's modulus E, which characterizes the ability of the material to resist tensile stress under elastic deformation, is E ≃ 75 ± 10 GPa, whereas E for the B19' phase is E ≃ 32 ± 9 GPa [see Table 1].The yield strength, which determines the limit of elastic behavior and the onset of plastic deformation on the stress-strain curve, is σ y ≃ 0.8 ± 0.1 GPa for austenite and σ y ≃ 0.1 ± 0.03 GPa for martensite.The ultimate strength, which is the stress at which the material begins to fracture, is σ f ≃ 1.9 ± 0.3 GPa [47,48,49].The ultimate strain ϵ f , which corresponds to ultimate strength, is ϵ f > 15 % and depends on the quality of the raw material as well as on the production protocol of the crystalline alloy (i.e. on the degree of mixing of the alloy components, on the heating temperature, on the cooling rate of the melt).
(II) Nanoporous alloy (ϕ = 0.55, l ≃ 4.4 nm).According to the classical definition [50], the nanoporous material is a porous material with characteristic pore sizes of less than 100 nm.In this case, we propose to consider NiTi with the porosity ϕ = 0.55, where the average linear pore size is l ≃ 4.4 nm.This value of the quantity l follows from the pore distribution P (l) over linear sizes obtained by molecular dynamics simulations of the nanoporous alloy and shown in Figure 3.The obtained distribution P (l) is the result of averaging over five independent molecular dynamics simulations performed from different initial configurations under identical thermodynamic conditions.Nonequilibrium molecular dynamics simulations of this nanoporous alloy have been used to determine stress-strain curves and to identify key mechanical properties.Figure 4(a) shows the stress-strain curves for uniaxial tensile strain obtained for different porous systems from independent molecular dynamics calculations.These stress-strain curves vary widely due to the inhomogeneous structure of the porous systems.Therefore, these curves are averaged and the most probable stress-strain curve is found [in Figure 4(a) this curve is shown in red].The Young's modulus E is determined as the slope of the linear part of this curve in the low strain region corresponding to the elastic regime.
The yield strength (or conventional yield strength) σ y is defined as the stress corresponding to 0.2% plastic deformation.The ultimate strength σ f is the maximum stress on the averaged stressstrain curve.Thus, for uniaxial tensile strain, we find that nanoporous NiTi is characterized by the Young's modulus E ≃ 17 ± 3 GPa, the yield strength σ y ≃ 0.38 ± 0.15 GPa and the ultimate strength σ f ≃ 1.1 ± 0.3 GPa.The limit value of the deformation ϵ f , at which the destruction of the nanoporous alloy is initiated, is ϵ f ∼ 13.4 %, which is comparable to the tensile limit of monolithic NiTi.Details of the molecular dynamics calculations associated with the preparation of nanoporous NiTi and the determination of mechanical properties are given in "Appendix: Experimental and simulation details".(III) Porous alloy with micron-sized pores (ϕ = 0.55, l ≃ 78 µm).In this case, we propose to consider NiTi with the porosity ϕ = 0.55 and whose characteristic pore size is l ≃ 78 µm.This alloy was prepared by powder sintering (experimental details are given in the section "Appendix: Experimental and simulation details").The value l ≃ 78 µm was obtained from the pore size distribution P (l), which, in turn, was determined using the stereometric method of random section planes.The pore size distribution shown in Figure 3 is the result of averaging over different thin sections.The obtained samples show the percolating porous structure, where the voids are in the form of channels of approximately the same width, mainly due to the peculiarities of sintering the mixture of powdered raw material.The mechanical properties of porous NiTi with micronsized pores are determined from the stress-strain curve obtained by averaging the experimental results for different samples under identical thermodynamic and deformation conditions.Figure 4(b) shows that the transition region from elastic to plastic deformation is blurred at strains < 0.5%.This is mainly due to the fact that structural transitions between the austenitic B2-phase and the martensitic B19'-phase occur at the tensile strain, which causes the pseudoplastic behavior of the porous material at strains from 0.1 to 0.3%.As a result, it was found that for porous NiTi alloys the Young's modulus is E ≃ 5.85 ± 0.8 GPa, the yield strength is σ y ≃ (4.2 ± 3) × 10 −3 GPa and the ultimate strength is σ f ≃ (17.8 ± 8) × 10 −3 GPa [see Table 1].Such relatively low strength characteristics in the case of samples with micron-sized pores can be related to the quality of the initial powdered raw material and the powder sintering conditions [51].Namely, in the case of porous NiTi, the use of powder metallurgy methods can lead to the formation of defects inherent in the crystal structure, as well as the formation of nano-sized cracks at the places of powder sintering, which significantly reduce the strength characteristics of the porous alloy [52,53].

Results and discussion
Let us define the general conditions for an arbitrary mechanical property from the set M as a quantity depending on the average linear size l and on the porosity ϕ.First, it must satisfy the conditions (5) characterising that the value of M decreases with increasing l and ϕ.Second, in the case of the monolithic material, we have l = l0 , ϕ = 0 and M( l, ϕ) = M 0 .
In addition, if there is only one void or an extremely small number of voids with the linear size l, which is much smaller than the size of the material 3 √ V , then the values of the mechanical properties will also approach the values of the monolithic material: Furthermore, we can formulate the following boundary condition: which corresponds to the monolithic material with a set of extremely small voids.And, finally, for the limit situation, where the material is completely absent, we have For the ϕ-dependence of an arbitrary mechanical characteristic from the set M, we propose to use the empirical Bal'shin equation ( 9), which ensures the fulfillment of the conditions (10a) and (10d).On the other hand, we obtain that the l-dependence of an arbitrary mechanical characteristic is reproduced by the power law of the form where the exponent α depends on the corresponding quantity from the set M. The expression (11) correctly reproduces the values of the mechanical properties for the cases considered above -the monolithic alloy, the nanoporous alloy and the alloy with micron-sized pores, as shown in Figure 5.
Then, taking into account the conditions ( 5) and ( 10) as well as Eqs.( 9) and ( 11), we find the following general equation on the complete fulfilment of the condition which is performed for the full set of mechanical properties M = {E, σ y , σ f }.In Eq. ( 12), the quantity V is the volume of the system.We are not aware of any models that produce an equation for the mechanical characteristic as a function of two parameters: porosity ϕ and average linear pore size l, as this is accounted for in the proposed equation (12).The closest in its physical meaning to our model ( 12) is the Knudsen model with equation which is applied to porous ceramics with polycrystalline matrix [54].Equation ( 13) allows us to approximate the dependence of the strength properties taking into account the porosity ϕ and the parameter G characterising the size of crystalline grains in the material.Here, K, a and b are empirical parameters.Equation ( 13) is given as the product of two contributions, where the first contribution is defined as a power-law dependence on the grain size G, and the second contribution is defined as an exponentially decaying dependence on the porosity ϕ.However, it should be noted that the Knudsen model does not apply to porous alloys, but to porous ceramics [55,56].The Bal'shin model with equation ( 9) does not take into account the dependence of M on the average linear pore size l.Nevertheless, it can be generalised to the following form: Comparing ( 12) and ( 14), we can see that the generalization of the Bal'shin model represents only a special case accounted by equation ( 12), where M 0 ( l0 / l) γϕ has a weak dependence on ϕ.
Equation ( 12) has the following features.First of all, this equation reproduces the nature of the change in mechanical properties with the change in the porosity ϕ and in the average linear pore size l.In this case, the equation depends on the range of values of l and ϕ defined by relation (12b).
Here, the quantity M 0 is the value of the corresponding mechanical property for the case of the monolithic alloy, and the quantity l0 is the average interatomic distance, as mentioned above.In the case of monolithic NiTi alloy, we have l0 ≃ 0.09 nm.Equation ( 12) contains three dimensionless parameters n, γ and λ which take positive values.The parameter n is identical to the parameter in the Bal'shin equation, and in the case of NiTi it takes the value n = 3 for all considered mechanical properties.Furthermore, the parameters γ and λ characterize how the mechanical property of the material decreases with increasing void size: the larger the value of the parameters γ and λ, the faster the value of the mechanical property decreases.In this case, the parameter γ corrects the value of the mechanical characteristic mainly in the limit of large values of the average linear pore size l.The parameter λ is responsible for changing the mechanical characteristic in the limit of small values of l.It is noteworthy that in the limit of large values of l, and also when γ → 0 and λ >> 0, Eq. ( 12) tends to the Bal'shin equation ( 9).In the case of NiTi, we find that γ = 0.028 ± 0.002 and λ = 0.25 ± 0.03 for the Young's modulus, while γ = 0.65 ± 0.04 and λ = 0 for the yield strength σ y and the ultimate strength σ f .Figure 6 shows the Young's modulus E, the yield stress σ y and the ultimate strength σ f as functions of the porosity ϕ and the average linear pore size l calculated for NiTi alloy using Eq. ( 12).It can be seen from this figure that the values of the mechanical properties approach the strength parameters of monolithic NiTi in the limit of small ϕ and l.On the other hand, when the pore sizes become comparable to the size of the system, i.e. l − l0 = 3 √ ϕV , the theoretical limit is reached, which satisfies the condition (3), under which the material does not physically exist.At the same time, the change in the mechanical properties is more pronounced at the porosity ϕ ∈ (0, 0.6] and/or at the average linear pore sizes l < 10 nm.The same result is obtained from the analysis of Eq. ( 12) with the parameters n, γ and λ defined for NiTi.Thus, the introduction of insignificant porosity into an almost homogeneous monolithic material leads to a significant reduction in its mechanical properties.In the limit of high porosity (i.e. at ϕ > 0.6) and the average linear pore size corresponding to a highly porous material with micron or millimeter sized pores, the mechanical properties are weakly dependent on ϕ and l.Such porosity conditions correspond to a material in which the porous structure is percolating, and the increase in porosity is mainly due to the thinning of the interpore partitions.
The mean relative error (MRE) between the calculated values of the mechanical properties and the result of Eq. ( 12) is ≈ 4.3 % in the case of the Young's modulus E [see Table 2].For the yield strength σ y and the ultimate strength σ f , the MRE takes the value ∼ 28 %, which is to be expected since the values of these quantities can vary over a wide range [48,49,57].In this case, the relative error can be reduced by analyzing a large set of empirical and simulation data for NiTi samples with different porosity and with different average linear pore size.In addition, given the recent advances in the application of the artificial intelligence methodology to solving physical problems, it is expected that the problem of increasing the amount of data for the mechanical properties can be solved using machine learning methods [58,59].
predict the strength properties of porous alloys with different compositions and porosity parameters if sufficient empirical and simulation data are available.
The results obtained in the present study allow us to approach the solution of actual problems of fundamental character, which may include: (i) understanding the mechanisms by which crack formation and growth occur in porous metal alloys; (ii) understanding how porous structure parameters, such as pore size distribution, interpore partition size distribution, pore shape, pore wall surface area, porosity, etc., affect the elastic and plastic properties and strength of porous metal alloys; (iii) explanation of the influence of structural heterogeneity arising due to the non-sphericity of pores, their orientation and distribution in the system on the mechanical properties of the porous alloy.The solution of these problems is only possible by combining experimental results for macroscopic mechanical properties with the results of molecular dynamics simulations of non-equilibrium processes on microscopic scales.

Acknowledgment
The work was supported by the Russian Science Foundation (project No.  Instron, USA) at a strain rate not higher than 7.5 × 10 −6 m/sec.The surface of the pore walls, the microstructure and the fracture fractograms of the samples were studied by using Quanta 3D scanning electron microscopes (manufacturer: FEI Company, USA) in the secondary electron mode in high vacuum and at the accelerating voltages from 20 to 30 kV and the beam sizes 5-20 nm.

Computational details and methods
Molecular dynamics simulation of crystalline NiTi has been carried out for the sample with the B2 crystal lattice consisting 12 577 atoms of Ni and 12 577 atoms of Ti.The atoms have been placed in a cubic system with the linear size L ≃ 9 nm.The interaction between Ni and Ti atoms is determined using the second nearest neighbour modified embedded atom method (2NN MEAM) [60,61,62]: This potential reproduces the physical properties of NiTi for a wide thermodynamic region, as has been shown previously in the works [63,64,65].In Eq. ( 15), E is the total energy of the system; F i (ρ i ) is the embedding function for the atom i within the background electron density ρ i ; the pair potential ϕ ij (r ij ) and the screening function S ij are evaluated at the distance r ij between atoms i and j; the cutoff radius of this potential is 5.0 Å for the considered binary system.
Porous NiTi samples were formed by cutting out the voids through by removing some of the atoms from the monolithic system.Only those atoms that were located inside imaginary ellipsoids were removed from the system.The positions and spatial orientations of these ellipsoids were determined randomly.At the same time, the ellipsoids can intersect with each other to form open voids.The values of the semi-axes of the ellipsoids were set randomly and were chosen from the range [2.0; 6.0] nm.Changing the semi-axes of the ellipsoids in the considered range makes it possible to obtain an open-type porous structure with the porosity ϕ = 0.55, which coincides with the porosity of the experimental samples.After cutting out the voids, the sample was brought to a state of thermodynamic equilibrium at the temperature T = 300 K and the pressure 1 atm for the time 60 ps.This time is sufficient to minimise the free energy of the void surfaces and to obtain a stable porous structure.The molecular dynamics simulation was performed using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) software (Sandia National Laboratories, USA) [66].The 3D visualisation of the simulation results and the analysis of the porous structure were performed using the OVITO software (OVITO GmbH, Germany) [46].
The resulting porous samples were subjected to uniaxial tensile stress along the OX axis at the constant rate 1.0 • 10 10 s −1 .It is important to note that such ultra-high strain rates are typical for the molecular dynamics simulation method used to study the process of material destruction under high dynamic loading [67,68].This rate differs from the experimental strain rate, which can lead to overestimation of mechanical properties in the case of molecular dynamics simulation, for example, as shown in Refs.[69,70].The value of the stress σ xx was determined using the Irwin-Kirkwood equation [71]: Here, V is the volume of the system; m i is the mass of the ith atom; ϑ xi is the x velocity component of the ith atom; F ijx is the x component of the force between the particles with the labels i and j; r ijx is the x component of the distance between the particles i and j.The tensile strain ϵ was determined by the relation ϵ(t) = (L x (t) − L (init) x )/L (init) x , where, L (init) x is the length of the simulation box before deformation, and L x (t) is this length at time t after the start of deformation.

Figure 1 :
Figure 1: Diagram demonstrating the possible range of l and ϕ for which a porous material is realized.Here, the average linear size l of the pores is measured in the reduced unit 3 √ V and takes values in the range ( l− l0 )/ 3 √ V ∈ (0; 1); V is the volume of the material.The thick red line indicates the theoretical limit, which results from equation l − l0 = 3 √ ϕV .A porous material can only be realized in the coloured region, where the condition l − l0 ≤ 3 √ ϕV is satisfied.Limit values ϕ = 0 and ( l − l0 )/ 3 √ V = 0 correspond to the monolithic material, while ϕ = 1 and solve the problem formulated above, we propose to consider monolithic NiTi alloy without pores as well as porous NiTi with micron-sized pores prepared by powder sintering method [Figure 2(a)] and nanoporous NiTi obtained through simulation [Figure 2(b)].The following conditions are considered here:

Figure 2 :
Figure 2: (a) Thin section of porous NiTi with micron-sized pores obtained by powder sintering.The pores and the solid matrix are identified by scanning electron microscopy through analysis of the thin sections of porous samples, where the pores are black because they are filled with impregnating resin and the solid frame is grey in color.(b) Snapshot of nanoporous NiTi obtained by simulation using OVITO software [46].The solid matrix and the pores are indicated by arrows.

Figure 3 :
Figure 3: Pore distributions P (l) over linear sizes obtained for nanoporous NiTi (through simulation) and for porous NiTi with micron-sized pores (by experiment) with the porosity ϕ = 0.55 at the temperature T = 300 K. Inset: image of the thin section obtained by powder sintering, where interpore partitions (light areas) and voids (dark areas) areshown.The linear pore size l is represented as the reduced quantity l/ l, where l is the average linear size of the pores.

Figure 4 :
Figure 4: Stress-strain curves of porous NiTi under uniaxial tension obtained for samples with (a) nano-sized pores and (b) micron-sized pores by means of non-equilibrium molecular dynamics (MD) simulations and experimental studies (Exp.).Insets: thin sections of nanoporous NiTi and porous NiTi with the characteristic pore size l ≃ 78 µm, obtained from MD simulation and from experiment, respectively.The red solid curves represent the averaged results.

Figure 5 :
Figure 5: Dependence of the mechanical properties of NiTi on the reduced average linear pore size log( l/ l0 ).The solid curve is the result of equation M = M 0 ( l0 / l) α , where α ≃ 0.43 for the Young's modulus E and α ≃ 0.22 for the yield strength σy and the ultimate strength σ f .

Figure 6 :
Figure 6: Dependence of the mechanical properties of porous NiTi on the reduced average linear pore size log( l/ l0 ) and on the porosity ϕ: (a) Young's modulus E, (b) yield strength σy and (c) ultimate strength σ f .In all these figures, the curved surface is the result of Eq. (12) obtained with different values of ϕ and l.The red curve is the slice of the M( l, ϕ) surface at the porosity ϕ = 0.55.The round markers indicate the mechanical properties found for porous NiTi with nano-sized pores and for NiTi with micron-sized pores.

Ti 3
Ni 4 (7 wt.%).Single sintering was carried out in the electric vacuum furnace SNVE-1.31/16-I4(manufacturer: VARP Vacuum, Russia) for 15 min at the temperature 1528±5 K.The temperature was controlled by means of a tungsten-rhenium alloy thermocouple.The specified temperaturetime regime allows one to obtain experimental samples with an optimal degree of sintering with a high quality of interparticle contacts, a regular porous structure and the lowest possible degree of shrinkage.Sintering was carried out at the pressure 6.65 • 10 −4 Pa with the average heating rate 10 K/min.The porosity of the samples was determined by the weighing method and is equal to ϕ = 0.55.The porous samples have a cylindrical shape with the diameter 11-13 mm and the length 65-80 mm.The sizes of the pores and the interpore partitions were determined by combination of the cutting plane method.The study of physical and mechanical properties was carried out by the tensile method on the universal electromechanical testing machine Instron 3369 (manufacturer:

Table 1 :
[47,48,49] properties and parameters of monolithic NiTi and porous NiTi: porosity ϕ, average linear size l, Young's modulus E, yield strength σy, ultimate strength σ f , strain at the yield strength ϵy and strain at the ultimate strength σ f .The experimental values of the mechanical properties for monolithic NiTi are taken from Refs.[47,48,49], whereas the mechanical properties of porous NiTi with micron-sized pores and nano-sized pores are defined in the present study.

Table 2 :
(12)lts of Eq.(12)for the mechanical properties of porous NiTi with nano-and micron-sized pores as well as the mean relative error (MRE) (|M pred − M calc |/M calc ) × 100%, where M calc are the values of the calculated mechanical properties from Table1, M pred is the result of Eq. (12).