The Effect of Hole Geometry on the Nonlinear Nanomechanics of γ-Graphyne Structures: A Finite Element Analysis

1. Introduction

The introduction of monolayer graphene in 2004 was a scientific breakthrough in the field of materials due to its combination of superior mechanical, electromagnetic, optical, and thermal properties [1,2]. In recent years, efforts have been made to discover materials and structures at the nanoscale that are akin to graphene or possess superior physical and structural attributes compared to graphene. Among the most promising recently identified graphene-like nanomaterials is the family of graphyne materials, initially examined theoretically by Balaban and colleagues [3], and subsequently by Narita and Nagai through first-principles analysis [4]. Graphyne represents a novel category of two-dimensional (2D) substances that have recently garnered interest owing to their distinct physical and chemical characteristics. Comprising sp- and sp²-bonded carbon atoms, these substances exhibit a honeycomb configuration akin to graphene, but with extra atoms or functional groups affixed to the carbon atoms. The presence of these additional atoms or groups gives rise to a wide range of properties, including enhanced electrical conductivity, mechanical strength, and chemical stability [5,6,7].

Despite the successful synthesis of graphdiyne by Li et al. in 2010 [8], there is limited experimental data available on the mechanical properties of graphyne structures in the literature. According to Ivanovskii [9], various types of graphyne and graphdiynes were thoroughly investigated for basic structural and mechanical properties, which were determined by experiments or evaluated by ab initio theory. Additionally, as reported by Zhang et al. [10], the Raman spectra of graphyne and graphdiyne were examined in a systematic manner, and the behavior of these materials under mechanical stress and strain was also investigated. A study by Xiao et al. [11] used atomic force microscopy experiments to determine the thermomechanical properties of graphdiyne films.

The effective design of advanced applications using nanostructures requires a thorough understanding of their properties and an ability to predict them using reliable theoretical or computational tools. Numerous molecular dynamics (MD) simulations and ab initio studies have been carried out on structures like graphene, carbon nanotubes, and fullerenes to comprehend a broad spectrum of properties [12,13,14] and behaviors [15,16,17]. Additionally, the finite element method has been used to provide a more comprehensive approach to the computational modeling of nanostructures. With regards to nanostructures, several forms of finite elements, such as beams [18,19,20], spring-like elements [21,22,23], shells [24,25], and bar elements [26,27], have been formulated to analyze the mechanical performance of these materials.

Atomistic finite element models simulate interatomic interactions of atoms using appropriate finite element formulations. Couto and Silvestre [28] introduced a FEA modeling approach to estimate the elastic behavior of graphyne. They conducted uniaxial and biaxial tensile simulations using computational models of zigzag and armchair configurations and compared their results with equivalent molecular dynamics and density functional theory models finding that they were in good agreement. Lee et al. [29] applied an atomistic FE model to analyze the sensitivity of various graphyne-based resonators under different boundary conditions and size configurations. They compared the results with those of graphene-based resonators and concluded that graphyne derives higher sensitivity. Georgantzinos et al. [30,31] investigated the elastic and elastoplastic mechanical properties of γ-graphyne family structures utilizing FE analysis. They conducted various parametric studies, concentrating on the size-dependency of the sheet structure on the mechanical properties, and generalized their results using regression models. Galhofo and Silvestre [32] derived stress-strain curves to examine both the monotonic and hysteretic behavior of γ-graphyne sheets using the same geometric configurations as in [28]. Siampanis et al. [33] examined the stress-strain behavior of graphyne structures under shear loading and obtained higher shear strength and stiffness compared to other structures of this family.

The available studies on the mechanical properties of graphyne structures by numerical methods assume ideal sheet structures without any defects or holes. However, in the case of graphene, a particular interest in the modeling of the structure lies in the assessment of the influence of nano-holes on the mechanical beahvior. This is in line with the investigation of nanoporous graphene or graphene sheet with holes that have been studied due to its good properties for gas absorption [34], water desalination [35], energy storage [36,37], and nanoelectronics applications [38]. Genoese et al. [39] introduced a new interatomic potential, the Dreading potential in which they added damping functions, to investigate the in-plane and out-plane tensile behavior of monolayer graphene sheets. In their study, they examined the effect of the diameter of a hole located in the middle of the sheet on the stress and strain distributions. Another study on the effect of holes on graphene sheets can be found in [40] where Muraru et al. developed a new software tool to generate graphene-based molecular models. Yoon et al. [41] examined the mechanical resonance properties of porous graphene sheets computationally and experimentally. However, graphyne structures fall behind in terms of published papers regarding the role of holes and their effect on mechanical properties. Recently, Lee et al. [42] used an atomistic FE model to examine the effect of holes on Young’s moduli of porous γ-graphyne sheets under tensile loading, concluding that an increase in porosity decreases Young’s modulus.

The presence of holes in a graphyne structure can play a practical role in several areas such as electronics, energy storage, and catalysis. In electronics, the holes can be used to create electronic devices with high electron mobility, such as field-effect transistors. In energy storage, graphyne with holes can be used as a high-capacity electrode material in lithium-ion batteries. The holes can also be used as active sites for catalytic reactions in chemical synthesis, such as hydrogenation reactions. The size and shape of the holes can also be used to control the behavior of the graphyne, such as its strength, electronic conductivity, and chemical reactivity. Additionally, by controlling the size and shape of the holes, it may be possible to create new types of graphyne structures with unique properties that can be used in a variety of applications. As a recent addition to the family of carbon allotropes, graphyne structures have not been extensively studied in comparison to other carbon nanostructures. In this study, we investigate the effect of hole geometry on the mechanical nonlinear response of γ-graphyne structures, using an efficient nonlinear finite element method (FEM). Specifically, we examine the mechanical properties of graphyne, graphdiyne, graphyne-3, and graphyne-4 structures with circular, triangle, and rhombus hole configurations. The implementation of the FEA model is done under large deformations, and it is an atomistic spring-based finite element model, in which a unique spring is used to simulate the interatomic interactions between carbon atoms.

2. Computational Model

2.1. Geometric Characteristics

Graphyne consists of three types of bonds: aromatic, single, and triple. The repeatable unit cell of the structure is the single-triple-single bond interconnection. The number of the unit-cell triple bonds defines the graphyne variation, i.e., one bond yields graphyne, two bonds yield graphdiyne, three bonds yield graphyne-3 and four bonds yield graphyne-4. Figure 1 shows the structural configuration of a graphyne sheet along with the bond variations.

Figure 2 depicts the geometric failures, which are the circular hole, the square-type hole and rhombus-type hole. The diameter of the circular hole is denoted as D, the length of the side of the square-type hole is A, and the lengths of the rhombus-type hole are B and C. The angle at all corners of the rhombus, in the case studied, is 90 degrees. The holes are located at the geometrical center of the structure.

The effects of the geometric parameters of the holes are indirectly calculated via the calculation of the area of the hole. The dimensions of all γ- graphyne models were approximately 10 nm × 10 nm. The investigations are carried out for five different aspect ratios, including 5%, 10%, 15%, 20% and 25%. The hole aspect ratio AR is obtained by dividing the area of the hole (Af) by the total area of each nanosheet (At).

$$AR=\frac{Af}{At}$$

(1)

This approach allows for a more generalized analysis of the impact of hole geometry on the mechanical properties of the graphyne structures, as it accounts for variations in hole size and shape by normalizing the hole area to the total area of the nanosheet.

2.2. Force Field Description

The aggregate potential energy of a graphyne molecular arrangement can be depicted as the cumulative energies from interatomic interactions. Considering that nonbonded interactions are negligible for the molecular systems under investigation, and only in-plane attributes are being assessed, the potential energy of a γ-graphyne configuration is articulated as the cumulative energies of bond stretching and bond angle bending:

$$U_{\mathrm{tot}}=\sum U_r^i+\sum U_{\theta }^i,$$

(2)

where $U_r^i$ represents the bond stretching term and $U_{\theta }^i$ denotes the bond angle bending term, respectively. The parameter i denotes the kind of bond, assuming the values s, a, and t corresponding to the single, aromatic, and triple bond, respectively. According to the Morse potential field [31], the potential energy terms can be written as

$$U_r^i=D_e^i\left\{\left[1-e^{-{{\rm B}}^i\left(r^i-r_0^i\right)}\right){]}^2-1\right\}.$$

(3)

$${\mathrm{U}}_{\mathsf{\theta }}^{i,\mathrm{j}}=\frac{1}{2}{k}_{\theta }^{i,\mathrm{j}}{({\theta }^{i,\mathsf{\beta }}-{\theta }_0^{i,\mathrm{j}})}^2[1+k_{\mathrm{sextic}}^{i,\mathrm{j}}\left({\theta }^{i,\mathrm{j}}-{\theta }_0^{i,\mathrm{j}}{)}^4\right]$$

(4)

In Equation (2), $r^i$ and $r_0^i$ represent the changed and original bond lengths, respectively. $D_e^i$ and ${{\rm B}}^i$ are constants that vary based on the type of bond (s, a, t). In Equation (3), ${\theta }^{i,\mathrm{j}}$ and ${\theta }_0^{i,\mathrm{j}}$ are the angles between the deformed and original state of two adjacent bonds i and j, respectively. $k_{\theta }^{i,\mathrm{j}}$and $k_{\mathrm{sextic}}^{i,\mathrm{j}}$ are constants that control the linear and nonlinear force-deformation relationship of bond angle bending interactions. It is important to note that nonbonded interatomic interactions and resulting nonlocality typically have minimal impact on the overall mechanical behavior of graphene-like nanostructures and can be omitted in the mathematical formulation.

These formulas can be differentiated to mimic the force field inside the nanostructure using a mix of translational and rotational springs. The longitudinal rigidity of the translational springs is determined by differentiating Equation (2) concerning the variation in bond length ($\Delta r^i=r^i-r_0^i$), and the necessary rigidity of the rotational spring can be obtained by differentiating Equation (3) concerning the variation in the bending angle ($\Delta {\theta }^{i,\mathrm{j}}={\theta }^{i,\mathrm{j}}-{\theta }_0^{i,\mathrm{j}}$),

$$F_r^i\left(\Delta r^i\right)=\frac{\partial U_r^i}{\partial \left(\Delta r^i\right)}=2B^i{D}_e^i\left(1-e^{-B^i\Delta r^i}\right)e^{-B^i\Delta r^i},$$

(5)

$$M_{\theta }^{i,j}\left(\Delta {\theta }^{i,j}\right)=\frac{\partial U_{\theta }^{i,j}}{\partial {(\Delta {\theta }^{i,j})}^2}=k_{\theta }^{i,j}\Delta {\theta }^{i,j}\left[1+4\stackrel{\approx 0}{\stackrel{�}{k_{\mathrm{sextic}}^{i,j}{(\Delta {\theta }^{i,j})}^3}}\right]\approx k_{\theta }^{i,j}\Delta {\theta }^{i,j},$$

(6)

where the developed axial force between two bonded atoms $F_r^i\left(\Delta r^i\right)$ is a result of a change in their inter-atomic distance, while a bending moment $M_{\theta }^{i,j}\left(\Delta {\theta }^{i,j}\right)$ arises from a change in the angle between two linked bonds.

2.3. Finite Element Model

The simulation of potential energies resulting from bond length and bending angle variations is consistent with the finite element formulation outlined in [30,32]. This is achieved by utilizing two-nodes, spring-based line elements that connect two bonded atoms. The longitudinal rigidity of the translational springs is determined by differentiating equation (4) concerning the variation in bond length, while the rotational spring is derived by differentiating equation (5) with respect to bending angle variation.

$$\frac{\partial F_r^i}{\partial \left(\Delta r^i\right)}=2B^i{D}_e^i\left(2-e^{-B^i\Delta r^i}\right)e^{-2B^i\Delta r^i},\frac{\partial M_{\theta }^{i,j}}{\partial {(\Delta {\theta }^{i,j})}^2}\approx k_{\theta }^{i,j}.$$

(7)

The stiffness matrix of the finite elements and their force-displacement behavior are deduced by using a local Cartesian coordinate system ($\overline{x},\overline{y}$). This results in a 2x2 elemental matrix that incorporates two fundemantal coefficients in the $\overline{x}$ and $\overline{y}$ directions:

$$k^{\mathrm{el}}=\left[\begin{array}{cccc}{k}_{\overline{x}}^{\mathrm{el}}& 0& -k_{\overline{x}}^{\mathrm{el}}& 0\\ 0& k_{\overline{y}}^{\mathrm{el}}& 0& -k_{\overline{y}}^{\mathrm{el}}\\ -k_{\overline{x}}^{\mathrm{el}}& 0& k_{\overline{x}}^{\mathrm{el}}& 0\\ 0& -k_{\overline{y}}^{\mathrm{el}}& 0& k_{\overline{y}}^{\mathrm{el}}\end{array}\right].$$

(8)

where the symbolization $\mathrm{el}$ of the element gets distinct string values: sas, ast, sts, and$\mathrm{tst}$indicating that the current formulation necessitates the implementation of four kinds of spring-like elements with varying longitudinal and transverse stiffness coefficients, as illustrated in [30,32]. For example, the notation 'sts' corresponds to a spring-based element used to describe a triple bond located between two single bonds. The nonlinear axial stiffness coefficients of these elements can be determined using the subsequent equation.

$$k_{\overline{x}}^{\mathrm{el}}=\{\begin{array}{c}2{(B^{\mathrm{s}})}^2{D}_e^{\mathrm{s}}(2-{\mathrm{e}}^{B^{\mathrm{s}}\Delta \overline{x}}){\mathrm{e}}^{-2B^{\mathrm{s}}\Delta \overline{x}},\mathrm{el}=\mathrm{ast},\mathrm{tst}\\ 2{(B^{\mathrm{a}})}^2{D}_e^{\mathrm{a}}(2-{\mathrm{e}}^{B^{\mathrm{a}}\Delta \overline{x}}){\mathrm{e}}^{-2B^{\mathrm{a}}\Delta \overline{x}},\mathrm{el}=\mathrm{sas}\\ 2{(B^{\mathrm{t}})}^2{D}_e^{\mathrm{t}}(2-{\mathrm{e}}^{B^{\mathrm{t}}\Delta \overline{x}}){\mathrm{e}}^{-2B^{\mathrm{t}}\Delta \overline{x}},\mathrm{el}=\mathrm{sts}\end{array},$$

(9)

$$k^{\mathrm{el}}=\left[\begin{array}{cccc}{k}_{\overline{x}}^{\mathrm{el}}& 0& -k_{\overline{x}}^{\mathrm{el}}& 0\\ 0& k_{\overline{y}}^{\mathrm{el}}& 0& -k_{\overline{y}}^{\mathrm{el}}\\ -k_{\overline{x}}^{\mathrm{el}}& 0& k_{\overline{x}}^{\mathrm{el}}& 0\\ 0& -k_{\overline{y}}^{\mathrm{el}}& 0& k_{\overline{y}}^{\mathrm{el}}\end{array}\right].$$

(10)

The notation el of the element is assigned one of four different string values: sas, ast, sts, and tst. This indicates that the current formulation necessitates the implementation of four types of spring-like elements, each with distinct longitudinal and transverse stiffness coefficients, as outlined in [21,31]. For example, the notation sts represents a spring-like element used to describe a triple bond between two single bonds. The nonlinear axial stiffness coefficients of these elements can be determined using the following equation.

$$k_{\overline{x}}^{\mathrm{el}}=\{\begin{array}{c}2{(B^{\mathrm{s}})}^2{D}_e^{\mathrm{s}}(2-{\mathrm{e}}^{B^{\mathrm{s}}\Delta \overline{x}}){\mathrm{e}}^{-2B^{\mathrm{s}}\Delta \overline{x}},\mathrm{el}=\mathrm{ast},\mathrm{tst}\\ 2{(B^{\mathrm{a}})}^2{D}_e^{\mathrm{a}}(2-{\mathrm{e}}^{B^{\mathrm{a}}\Delta \overline{x}}){\mathrm{e}}^{-2B^{\mathrm{a}}\Delta \overline{x}},\mathrm{el}=\mathrm{sas}\\ 2{(B^{\mathrm{t}})}^2{D}_e^{\mathrm{t}}(2-{\mathrm{e}}^{B^{\mathrm{t}}\Delta \overline{x}}){\mathrm{e}}^{-2B^{\mathrm{t}}\Delta \overline{x}},\mathrm{el}=\mathrm{sts}\end{array}.$$

(11)

In accordance with the simplified modeling technique, detailed in [42], the bond angle bending variation can be effectively addressed by utilizing an suitable rigidity coefficient in the$\overline{y}$-direction for each element $\mathrm{el}=\mathrm{sas},\mathrm{ast},\mathrm{sts},\mathrm{tst}$. This is provided by:

$$k_{\overline{x}}^{\mathrm{el}}=\{\begin{array}{c}2{(B^{\mathrm{s}})}^2{D}_e^{\mathrm{s}}(2-{\mathrm{e}}^{B^{\mathrm{s}}\Delta \overline{x}}){\mathrm{e}}^{-2B^{\mathrm{s}}\Delta \overline{x}},\mathrm{el}=\mathrm{ast},\mathrm{tst}\\ 2{(B^{\mathrm{a}})}^2{D}_e^{\mathrm{a}}(2-{\mathrm{e}}^{B^{\mathrm{a}}\Delta \overline{x}}){\mathrm{e}}^{-2B^{\mathrm{a}}\Delta \overline{x}},\mathrm{el}=\mathrm{sas}\\ 2{(B^{\mathrm{t}})}^2{D}_e^{\mathrm{t}}(2-{\mathrm{e}}^{B^{\mathrm{t}}\Delta \overline{x}}){\mathrm{e}}^{-2B^{\mathrm{t}}\Delta \overline{x}},\mathrm{el}=\mathrm{sts}\end{array}.$$

(12)

To investigate the elastoplastic behavior of a nanostructure, the equilibrium equation for each nonlinear spring-based element can be written as:

$$k^{\mathrm{el}}{u}^{\mathrm{el}}=f^{\mathrm{el}}.$$

(13)

The equilibrium equation for each hypothetical nonlinear spring-based element, where $u^{\mathrm{el}}$ and $f^{\mathrm{el}}$ represent the respective force and displacement vectors for each element. These vectors have the following notation if $i$ and $j$ are the two nodes of the finite element:

$$u^{\mathrm{el}}={\left[\begin{array}{cccc}{u}_{\overline{x}i}^{\mathrm{el}}& u_{\overline{y}i}^{\mathrm{el}}& u_{\overline{x}j}^{\mathrm{el}}& u_{\overline{y}j}^{\mathrm{el}}\end{array}\right]}^{\mathrm{T}},$$

(14)

$$f^{\mathrm{el}}={\left[\begin{array}{cccc}{f}_{\overline{x}i}^{\mathrm{el}}& f_{\overline{y}i}^{\mathrm{el}}& f_{\overline{x}j}^{\mathrm{el}}& f_{\overline{y}j}^{\mathrm{el}}\end{array}\right]}^{\mathrm{T}}.$$

(15)

The system of nonlinear equations can be assembled into its final form, through the transformation of the global coordinate system to the elemental stiffness equation for each finite element, in accordance with the constraints of nodal connectivity:

$$K^{\mathrm{el}}(U^{\mathrm{el}})U^{\mathrm{el}}=F^{\mathrm{el}},$$

(16)

where $F^{\mathrm{el}}$, $U^{\mathrm{el}}$, and $K^{\mathrm{el}}{(U)}^{\mathrm{el}}$is the assembled force vector, assembled displacement vector, and assembled deformation-dependent stiffness matrix, respectively.

To simulate the nonlinear mechanical tensile tests, appropriate supports and loads are applied as boundary conditions. The next step is to use an established incremental-iterative method based on the Newton-Raphson algorithm to address the nonlinearity of the global stiffness matrix numerically and determine the response of the structure as a result.

3. Results and Discussion

The proposed finite element model has been previously validated through comparisons with other available results in the open literature [31] regarding pristine graphyn structures. However, a number of comparisons were done between the outcomes of this computational approach and related data taken from various literature sources, as provided in Table 1, in order to validate the numerical calculation of the key mechanical properties of γ-graphynes. This comparison serves as further verification of the validity and accuracy of the proposed model. These sources employed different approaches to extract the data. The results indicate a reasonable agreement with the findings of previous studies and demonstrate that the model can be used to predict the mechanical properties of graphyne structures with holes with a high degree of accuracy, providing confidence in the accuracy of the computational scheme, and thus it can be useful for the design and optimization of graphyne-based materials and devices.

The behavior of the elastic modulus of graphyne structures (graphyne, graphdiyne, graphyne-3, graphyne-4) with increasing aspect ratio of the hole is an important aspect of understanding the mechanical properties of these materials. Figure 3 provides insights on this behavior illustrating the change in elastic modulus for different types of holes. Figure 3a,b show the behavior of the elastic modulus in the x and y axis for the case of a central circular hole for the four structures. Figure 3b,c depict the case of a square hole, while Figure 3c,d describe the case of a rhombus-type hole. The elastic modulus is seen to decrease as the size of the hole increases in all circumstances. Additionally, the rate of decrease in the elastic modulus with an increase in hole size is similar for all cases. This suggests that the graphyne structures have a relatively consistent response to the introduction of holes in their structure. However, small fluctuations in the elastic modulus may be attributed to the non-symmetric geometric effects caused by the hole creation, depending on the number of atoms and bonds cut off. As the length of the acetylene series increases, the elastic modulus is generally lower in all cases. This highlights the importance of the number of atoms and bonds in determining the mechanical properties of graphyne structures. Close values are observed concerning the elastic modulus in the y-direction between the graphyne and graphdiyne specifically, indicating that these two structures may have similar mechanical properties in certain directions.

The mechanical tensile strength of graphyne structures (graphyne, graphdiyne, graphyne-3, graphyne-4) seems to be considerably affected by the aspect ratio of the hole. Figure 4 illustrates the change in ultimate strength for different types of holes. Figure 4a,b show the ultimate strength in the x and y axis for the case of a central circular hole for the four structures. Figure 4b,c depict the case of a square hole, while Figure 4c,d describe the case of a rhombus-type hole. As the size of the hole increases, the structure can experience a significant reduction in strength, with losses of up to almost 50% observed for aspect ratios near 0.25. The rate of decrease in the ultimate strength with an increase in hole size is similar for all cases. Small fluctuations in the ultimate strength may be observed for the same reasons as previously.

As the length of the acetylene series increases, the ultimate strength is generally lower in all cases. This emphasizes the significance of the number of atoms and bonds in establishing the mechanical properties of graphyne structures.

In Figure 5, we observe the variation in fracture strain for different types of holes. Figure 5a,b present the fracture strain in the x and y axis for the case of a central circular hole for the four structures. Figure 5b,c depict the case of a square type- hole, while Figure 5c,d describe the case of a rhombus-type hole. It is noticeable that in all cases, the rate of decrease in the fracture strain with an increase in hole size is similar for all cases. Slight variations in the fracture strain can be attributed to non-symmetric geometric effects caused by the hole creation, depending on the number of atoms and bonds cut off. Similarly, here, as the length of the acetylene series increases, the fracture strain is generally lower in all cases.

We conducted a nonlinear regression analysis to extract analytical equations that capture the hole-dependent variations in mechanical properties of the four types of graphyne with central holes, aiming to enable a more systematic and simplified prediction of their mechanical properties. The following fitting 3-parameter function of hole aspect ratio was used for all the attempted property approximations:

$$Y/Y_0=a+b{e}^{-ARc},$$

(17)

where Y is the suggested mathematical model to fit the data computed for all examined properties E_x, E_y, σu_x, σu_y, εf_x, and εf_y whereas, a, b, and c are constant parameters that must be found using a different regression procedure for each mechanical property variation. Y₀ is the corresponding property of pristine graphyne structure. Eq. (20) from [31] is a valuable tool for predicting the mechanical properties of graphyne structures, by considering the effect of size variations. This equation can be used to estimate Y_o, an important parameter for understanding the mechanical behavior of these structures. By utilizing Eq. (17), which relates Y_o to the corresponding mechanical properties Y, predictions of any size can be made, providing a comprehensive understanding of the mechanical response of graphyne structures with holes.

By incorporating the FEM results and performing a regression analysis using Eq. (17), the optimized values of the parameters can be calculated. Appendix A presents the regression curves for each graphyne structure type, hole type, and mechanical property in Figure A1, Figure A2 and Figure A3. These figures also include the average curve obtained by considering a single equation for all graphyne structure types. Table 2 provides the parameter values for the unique equation applied to all graphyne structure types. For a more accurate prediction, Table A1 contains the values of the parameters for Eq. (17) that were obtained from the fitting process for mechanical properties data of graphyne and each hole type. The fitting process for the mechanical properties data based on Eq. (17) yielded the parameter values shown in Table A2, regarding the graphdiyne. Table A3 and A4 provide the parameter values obtained from fitting Eq. (17) to the mechanical properties data for each hole type in graphyne-3 and graphyne-4, respectively.

The nonlinear FEA method was implemented under large deformations to examine the mechanical properties of γ-graphyne structures with circular, triangle, and rhombus hole configurations. The study used an atomistic spring-based finite element model, simulating the interatomic interactions between carbon atoms. This approach provided a comprehensive understanding of the mechanical response of graphyne structures with holes, enabling a more systematic and simplified prediction of their mechanical properties. In addition to the FEA used in this study, the nonlinear response of γ-graphyne structures can also be investigated using several other robust numerical modeling methods. For instance, the Finite Difference Method [46], Bezier Multi-Step Method [47], and Differential Quadrature Method [48] can provide alternative approaches to understanding the nonlinear response of γ-graphyne structures and can be considered for future studies or comparisons.

The study methodically investigates the nonlinear mechanical response of γ-graphyne structures, including variations such as graphyne, graphdiyne, graphyne-3, and graphyne-4, under the influence of different hole geometries (circular, square, and rhombus) and aspect ratios (5%, 10%, 15%, 20%, 25%). The investigation reveals that hole geometry significantly impacts the elastic behavior, deformation, and strength of the structures, with the elastic modulus, ultimate strength, and fracture strain decreasing as hole size increases. This is attributed to the variations in stress distribution and deformation patterns caused by the holes. The study also presents analytical equations derived from regression analysis to capture hole-dependent variations in mechanical properties, providing a valuable tool for systematic and simplified prediction of mechanical properties. This comprehensive analysis underscores the importance of considering hole geometries in the design and optimization of graphyne-based materials and devices, ultimately aiding in the development of materials with desired mechanical properties for specific applications.

The study of hole geometry, encompassing circular, triangular, and rhombus configurations, is pivotal for comprehending the nonlinear mechanical behavior of γ-graphyne structures. Different hole geometries lead to variations in the stress distribution, deformation patterns, and ultimately, the mechanical properties of the material. Specifically, the hole geometry influences the initiation and propagation of cracks, affecting the material’s strength and failure mechanisms. This understanding is crucial for optimizing the design and fabrication of γ-graphyne structures, as it can lead to enhanced mechanical properties, making the material more suitable for various applications, such as in nanoelectronics and nanomechanical systems.

4. Conclusions

In this study, we investigated the effect of hole geometry on the mechanical nonlinear response of γ-graphyne structures, including graphyne, graphdiyne, graphyne-3, and graphyne-4 structures, under large deformations using a finite element analysis (FEA) method. We found that the size and shape of the holes significantly impacted the structures’ elastic behavior, deformation, and strength, highlighting the importance of considering hole geometries in the design and fabrication of graphyne-based materials and devices. The results indicated that the size and shape of the holes significantly impacted the structures’ elastic behavior, deformation, and strength. Specifically, the elastic modulus of the graphyne structures decreased with an increase in the size of the hole, and the ultimate strength and fracture strain experienced a similar rate of decrease with an increase in hole size. Small fluctuations in these properties were observed due to non-symmetric geometric effects caused by the hole creation. Our study provides a valuable tool for predicting the mechanical properties of graphyne structures with holes by incorporating the effect of size variations, enabling a more systematic and simplified prediction of their mechanical properties. The analytical equations extracted from the regression analysis, which capture the hole-dependent variations in mechanical properties, were presented. This study highlights the importance of considering the hole geometries in the design and fabrication of these materials, providing a comprehensive understanding of the mechanical response of graphyne structures with holes, which can be useful for the design and optimization of graphyne-based materials and devices. Overall, our findings can aid in the optimization of the design of graphyne structures for specific mechanical applications.