<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ENG</journal-id><journal-title-group><journal-title>Engineering</journal-title></journal-title-group><issn pub-type="epub">1947-3931</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/eng.2024.164009</article-id><article-id pub-id-type="publisher-id">ENG-132895</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Numerical Optimization by Finite Element Method of Stainless Steel/Glass-Epoxy Composite Bolted Joint under Tension and Compression
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Christian</surname><given-names>Schmitt</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Arnaud</surname><given-names>Kremeur</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Pawel</surname><given-names>Lipinski</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Julien</surname><given-names>Capelle</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>National Engineering School of Metz, University of Lorraine, Metz, France</addr-line></aff><pub-date pub-type="epub"><day>29</day><month>04</month><year>2024</year></pub-date><volume>16</volume><issue>04</issue><fpage>102</fpage><lpage>122</lpage><history><date date-type="received"><day>5,</day>	<month>December</month>	<year>2023</year></date><date date-type="rev-recd"><day>27,</day>	<month>April</month>	<year>2024</year>	</date><date date-type="accepted"><day>30,</day>	<month>April</month>	<year>2024</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The aim of this study was to optimize the geometry and the design of metallic/composite single bolted joints subjected to tension-compression loading. For this purpose, it was necessary to evaluate the stress state in each component of the bolted join. The multi-material assembly was based on the principle of double lap bolted joint. It was composed of a symmetrical balanced woven glass-epoxy composite material plate fastened to two stainless sheets using a stainless pre-stressed bolt. In order to optimize the design and the geometry of the assembly, ten configurations were proposed and studied: a classical simple bolted joint, two joints with an insert (a BigHead&lt;sup&gt;R&lt;/sup&gt; insert and a stair one) embedded in the composite, two &amp;#8220;waved&amp;#8221; solutions, three symmetrical configurations composed of a succession of metallic and composites layers, without a sleeve, with one and with two sleeves, and two non-symmetrical constituted of metallic and composites layers associated with a stair-insert (one with a sleeve and one without). A tridimensional Finite Element Method (FEM) was used to model each configuration mentioned above. The FE models taked into account the different materials, the effects of contact between the different sheets of the assembly and the pre-stress in the bolt. The stress state was analyzed in the composite part. The concept of stress concentration factor was used in order to evaluate the stress increase in the highly stressed regions and to compare the ten configurations studied. For this purpose, three stress concentration factors were defined: one for a monotonic loading in tension, another for a monotonic loading in compression, and the third for a tension-compression cyclic loading. The results of the FEM computations showed that the use of alternative metallic and composite layers associated with two sleeves gived low values of stress concentration factors, smaller than 1.4. In this case, there was no contact between the bolt and the composite part and the most stressed region was not the vicinity of the hole but the end of the longest layers of the metallic inserts.
 
</p></abstract><kwd-group><kwd>Bolted Joint</kwd><kwd> Glass-Epoxy Composite</kwd><kwd> Clearance</kwd><kwd> Hybrid Steel-Composite</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>Highlight</title><p>- Ten configurations of stainless steel/woven glass-epoxy bolted joints were modelled using finite element method.</p><p>- Stress concentration factors were evaluated in the composite material under tension and compression loads.</p><p>- The best solution is a hybrid composite bolted joint with alternative stainless steel and glass-epoxy plies and with two sleeves.</p></sec><sec id="s2"><title>1. Introduction</title><p>For the past decades, fiber reinforced polymers have been widely used in engineering fields such as for civil, aeronautic, aerospace, automotive, civil engineering or naval industries thank to their high specific mechanical properties (high strength and high stiffness to weight ratio) and good fatigue resistance. The improvement and optimization of production technologies and processes permit also to reduce the assembly costs because polymer composites components are easy to make in a single sheet.</p><p>Nevertheless, in most of time, engineering structures are made of many components which have to be connected together. Two techniques are usually used to connect structural composite components: mechanical fastening using bolted or riveted joints and adhesive bonding. Especially in aircraft engineering, for security reasons and for improving fatigue properties, hybrid connections based on the combination of these two techniques are widely employed [<xref ref-type="bibr" rid="scirp.132895-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref4">4</xref>] .</p><p>This study being focused on polymer composite bolted joints, thus adhesive bonded joints and bolted/bonded joints are not studied and a particular attention is focused on metallic/composite bolted connections and hybrid joints.</p><p>The main disadvantage of bolted joints comes from the high stressed concentration zone near the fastener hole, which can reduce the stiffness and the ultimate strength due to damages in the composite material such as matrix cracking, fiber failure or plies delamination.</p><p>In literature, a lot of papers treat failure and damage prediction of the bolted joint. Experiments and numerical modelling of single-lap and double-lap, single- and multi-bolted polymer composite joints have been investigated in numerous works [ [<xref ref-type="bibr" rid="scirp.132895-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref8">8</xref>] for instance]. Although, it is possible to reduce the stress level, avoid damage and increase the loads transfer by improving the local geometry and the configuration of the bolted joint. This aspect is rarely treated in literature. Numerous parameters such as joint geometry, bolt diameter, washer diameter, clearance between the bolt and the composite plate, clamping force which induces a pre-stress in the bolt, and stacking sequence can have a great influence on the behavior, the stress state, the stiffness and the ultimate strength of the assembly.</p><p>The effect of washer size and tightening torque on a glass-epoxy bolted joint were studied in [<xref ref-type="bibr" rid="scirp.132895-ref9">9</xref>] . It was found that the strength increases with increasing washer size at a constant tightening torque. When the diameter of the washer decreases, lateral compressive stress in the composite material under the washer occurs and leads to micro-cracks around the hole. The stiffness of a glass-epoxy bolted joint increases with increasing tightening torque thanks to the increasing of contact pressure [<xref ref-type="bibr" rid="scirp.132895-ref9">9</xref>] . For single lap-shear glass-epoxy bolted joints, the failure load increases with the increasing of the clamping torque [<xref ref-type="bibr" rid="scirp.132895-ref10">10</xref>] . Fengrui et al. [<xref ref-type="bibr" rid="scirp.132895-ref11">11</xref>] studied a carbon-epoxy laminate single-lap bolted joint. The authors found that the bolt-hole clearance had a significant influence on the distribution of load while the bolt tightening torque had a negligible effect.</p><p>Muc et al. [<xref ref-type="bibr" rid="scirp.132895-ref12">12</xref>] studied a local reinforcement around the hole of a composite plate using curved fibres as the reinforced material, but the authors pointed out that this technique does not permit to recover the same strength of the composite material without hole, due to the great stress concentration factor in the vicinity of the hole. Yingdan et al. [<xref ref-type="bibr" rid="scirp.132895-ref13">13</xref>] proposed a variable angle tow placement technique to create local reinforcement on composites plates with holes. This technology consists in placing continuous fibers which follow the pre-calculated trajectories of eigenstresses.</p><p>However, these two technologies have not yet been applied to composite bolted joints.</p><p>A hybrid carbon-epoxy composite/titanium alloy assembled with a titanium bolt was studied in [<xref ref-type="bibr" rid="scirp.132895-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref15">15</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref16">16</xref>] . A local reinforcement around the hole was made with alternative carbon-epoxy and titanium alloy layers. It was found a significant increase of the strength of the bolted joint compared to this without metallic layers. A spacecraft payload adaptor was made with CFRP—titanium hybrid laminates and experimentally tested in [<xref ref-type="bibr" rid="scirp.132895-ref17">17</xref>] . The authors pointed out the advantages of this technology. They found a great increase of the load capability of bolted joints, the joint efficiency and the joint robustness. They also pointed out a high potential for steel reinforcement.</p><p>As describes above, most of time, the studied structure is composed of carbon-epoxy plates joined with a steel or titanium bolt. Glass-epoxy composites and hybrid metal/composite structures are less studied.</p><p>The association of glass-epoxy/stainless steel could replace with advantage stainless steel only components in industrial automotive applications such as transport of gas or liquid under pressure in order to reduce the weight of the vehicle and the cost of the transport. In this purpose, the aim of this paper is the study of a single bolted joint under tension-compression made with two stainless steel sheets and a glass-epoxy material assembled with a pre-stress bolt.</p><p>The paper is organized as follows.</p><p>In the “Bolt joint” part, the general principle of the bolted joint studied, its geometry and the size of the sample used and materials chosen are described. Ten different configurations are proposed and defined. In the “Finite element modeling” part, calculi by finite elements method are performed to evaluate the stress state in the composite material. In the “stress analysis” part, in order to describe the increase of the tension-compression stress component in the most stressed zones of the composite material, the concept of stress concentration factor is exploited. All the ten configurations are analyzed and compared each other to select the best suited solution from a mechanical and technological point of view, taking into account the cost of machining of the bolted structure.</p></sec><sec id="s3"><title>2. Bolt Joints</title><sec id="s3_1"><title>2.1. Materials and Specimen</title><p>The bolt joint of this study is based on the principle of the double lap bolted connection. It is designed as a multi-material assembly. In order to evaluate the stress concentration level in the bolted joint area, the model showed in <xref ref-type="fig" rid="fig1">Figure 1</xref> is used and subjected to tension-compression loads. All dimensions are specified in this figure.</p><p>The specimen is composed of (<xref ref-type="fig" rid="fig2">Figure 2</xref>):</p><p>- Two 316 L stainless steel external plates having a length of l<sub>s</sub> = 150 mm, a width of w<sub>s</sub> = 50 mm and a thickness of t<sub>s</sub> = 5 mm. A metallic bead is placed between the two plates (see <xref ref-type="fig" rid="fig2">Figure 2</xref>) in order to allow an efficient fixation of the specimen into the first grip of the testing machine in future experimental tests.</p><p>- A 0/90 balanced woven glass-epoxy composite plate of the same size as the metallic plates, i.e. length l<sub>c</sub> = 150 mm, width w<sub>c</sub> = 50 mm and thickness t<sub>c</sub> = 5 mm. A symmetric stacking sequence with ten layers of same thickness (t<sub>l</sub> = 0.5 mm) is used to avoid undesired warping phenomenon during the tension or compression loading. Two beads are glued at the end of the composite plate in order to produce a balanced tension-compression loading by the second grip of the tensile machine.</p><p>The three plates of the bolt, i.e. two stainless steel sheets and the composite plate, are drilled and jointed together using a stainless steel M6 bolt of strength class A2-80. Its yield strength σ<sub>eb</sub> is equal to 600 MPa and its ultimate strength σ<sub>mb</sub> is equal to 800 MPa. The clearance between the composite material and the bolt is nil and this between the stainless steel plates and the bolt is equal to 0.5 mm. Since the fastener head and the nut washer are in contact with the steel sheets, there is no local crushing of the composite material which could damage it. The diameter of the washer has no influence on possible damage of composite since there is no contact between these two components.</p><p>The main function of the assembly described above is to transfer the load from the metallic part (stainless steel sheets) to the composite part (woven composite plate). For this reason, the analysis of the obtained results is mainly focused on the area around the bolt. Various designs of the bolted joint are studied and improvements are proposed leading to the lowest stress level in the composite material and taking into account technological aspects and cost.</p><p>Each specimen can be decomposed into a “standard part” which is the same for all configurations and an “individual part” which is specific to each configuration studied (<xref ref-type="fig" rid="fig2">Figure 2</xref>). The stress field is analyzed only in this “individual or central” part of the sample.</p></sec><sec id="s3_2"><title>2.2. Configuration 1.1—Simple Bolted Joint</title><p>This technological solution showed in <xref ref-type="fig" rid="fig1">Figure 1</xref> is the simplest one and is treated later on as reference configuration.</p></sec><sec id="s3_3"><title>2.3. Configurations 2.1 and 2.2—Bolted Joint with Aninsert</title><p>A metallic insert embedded in the composite material is introduced to improve the stress transfer through the bolted joint. Two configurations are studied.</p><p>A BigHead<sup>R</sup> insert is composed of a washer welded to a metallic element such as a threated rod, a nail, a screw or a rod [<xref ref-type="bibr" rid="scirp.132895-ref18">18</xref>] . The BigHead<sup>R</sup> insert chosen for this study (<xref ref-type="fig" rid="fig3">Figure 3</xref>) is composed of a washer of diameter 23 mm welded to a M6 nut (configuration 2.1). It can be embedded into the composite material during the manufacturing process. Its role is to increase the binding between the different parts and to allow a better distribution of the load among the concerned plates. However, this kind of joint is non symmetric and applied load can lead to unwanted shear and bending stresses in the sample.</p><p>Configuration 2.2 is a bolted joint with a symmetric insert (<xref ref-type="fig" rid="fig4">Figure 4</xref>(a)) having a circular shape with a diameter of 30 mm and embedded into the woven 0/90 glass-epoxy composite. Each step of the insert is 0.5 mm high, the same as a ply of the composite, except for the central step whose height is equal to 1 mm (<xref ref-type="fig" rid="fig4">Figure 4</xref>(b)).</p></sec><sec id="s3_4"><title>2.4. Configurations 3.1 and 3.2—“Waved” Bolted Joints</title><p>In order to increase the contact surface between the metallic and composite parts, the plates have here the same sinusoidal wave profiles. The load transfer should then be improved. The wave is defined by its half-period T/2 = 8 mm and its amplitude. Two amplitudes h<sub>1</sub> = 2 mm for configuration 3.1 and h<sub>2</sub> = 1 mm for configuration 3.2 are studied (<xref ref-type="fig" rid="fig5">Figure 5</xref>).</p></sec><sec id="s3_5"><title>2.5. Configurations 4.1 to 4.3—Hybrid Composite Bolted Joints</title><p>In order to improve the mechanical characteristics of the bolted joint, a hybrid composite material defined by a succession of metallic and composites layers is studied. The stacking sequence is symmetric, composed of six 0/90 woven glass-epoxy plies and four stainless steel layers. Each layer has the same thickness t<sub>l</sub> = 0.5 mm (<xref ref-type="fig" rid="fig6">Figure 6</xref>).</p><p>In order to avoid high stress concentration in the composite in the vicinity of the end of the metallic plate, a gradual transition of the length of the metallic layers is chosen. The metallic layers near the middle plane of the joint are 5 mm longer than the metallic plates (<xref ref-type="fig" rid="fig7">Figure 7</xref>).</p><p>Three hybrid composite bolted joints are studied (<xref ref-type="fig" rid="fig8">Figure 8</xref>):</p><p>- A simple hybrid composite bolted joint: <xref ref-type="fig" rid="fig8">Figure 8</xref>(a). (configuration 4.1)</p><p>- A hybrid composite bolted joint with a cylindrical sleeve of an external diameter d<sub>1s</sub> = 10 mm separating the bolt and the composite part (configuration 4.2). In this case, the metallic bolt remains in contact only with the metallic insert.</p><p>- A hybrid composite bolted joint with two sleeves and a shoulder (configuration 4.3). The first sleeve has the same external diameter as previously (d<sub>1s</sub> = 10 mm) and its thickness is greater than the thickness of the hybrid composite. The second sleeve has an external diameter d<sub>2s</sub> = 14 mm to avoid the contact between the hybrid composite and the metallic bolt (<xref ref-type="fig" rid="fig8">Figure 8</xref>(c)).</p></sec><sec id="s3_6"><title>2.6. Configurations 5.1 to 5.2—Hybrid Composite Bolted Joints with a Stair-Insert</title><p>For this solution, the hybrid composite material is reinforced with a circular metallic stair-insert embedded in the composite in the vicinity of the bolt.</p><p>Twonon-symmetric solutions are proposed:</p><p>- Hybrid composite bolted joint with a simple stair-insert (configurations 5.1) presented in <xref ref-type="fig" rid="fig9">Figure 9</xref>(a). The insert has a circular shape and is embedded in the hybrid composite. This configuration is similar to 4.2 except that the insert has replaced the sleeve.</p><p>- Hybrid composite bolted joint with insert and sleeve (configurations 5.2) illustrated in <xref ref-type="fig" rid="fig9">Figure 9</xref>(b). Configuration 5.2 is similar to 4.3 except that the insert has replaced the smallest sleeve.</p></sec></sec><sec id="s4"><title>3. Finite Element Modeling</title><p>Tension and compression tests are modelled using the finite element method. The software MSC Marc Mentat<sup>R</sup> is used. All the components of the bolted joints are modelled taken into account the particular constitutive behavior of the different materials, the non-linearity due to the contact (if existing) between the different sheets of the sample as well as the pre-stress in the bolt.</p><sec id="s4_1"><title>3.1. Meshing</title><p>A three dimensional model is built to simulate the bolt joint under tension or</p><p>compression. The specimens have been meshed with 8 nodes hexahedral elements using Hyper Mesh software.</p><p>The steel plates are modelled with 5 elements through the thickness, each element being 1 mm thick. Concerning the composite plate, the thickness of each element is chosen equal to those of each ply, i.e. 0.5 mm.</p><p>The mesh is sufficiently refined in the vicinity of the hole in the composite and steel plates in order to increase the accuracy in local area and predict the stress level correctly. For each configuration, ten tests are performed to ensure the convergence of the numerical solution.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>0 shows two typical examples of meshes of the “individual part” of the models (configurations 4.3 and 5.2).</p></sec><sec id="s4_2"><title>3.2. Mechanical Constitutive Law of the Materials</title><p>The behavior of each material of the bolted joint is assumed to be homogeneous, elastic and linear.</p><p>The stainless steel is supposed isotropic with Young’s modulus E<sub>st</sub> = 195,000 MPa and Poisson’s ratio v<sub>st</sub> = 0.3.</p><p>The balanced woven composite material is supposed orthotropic. Its characteristics are done in <xref ref-type="table" rid="table1">Table 1</xref> with respect to the coordinate system indicated in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p></sec><sec id="s4_3"><title>3.3. Loading and Boundary Conditions</title><p>The tensile strength of the woven composite material is supposed to be equal to 500 MPa. The aim of this paper is to choose the most appropriate bolted joint.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Mechanical properties of woven glass-epoxy composite material</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Young modulus (MPa)</th><th align="center" valign="middle" >Poisson ratio</th><th align="center" valign="middle" >Coulomb modulus (MPa)</th></tr></thead><tr><td align="center" valign="middle" >E<sub>x</sub>: 26,000</td><td align="center" valign="middle" >V<sub>xy</sub>: 0.13</td><td align="center" valign="middle" >G<sub>xy</sub>: 5170</td></tr><tr><td align="center" valign="middle" >E<sub>y</sub>: 26,000</td><td align="center" valign="middle" >V<sub>yz</sub>: 0.34</td><td align="center" valign="middle" >G<sub>yz</sub>: 5050</td></tr><tr><td align="center" valign="middle" >E<sub>z</sub>: 10,750</td><td align="center" valign="middle" >V<sub>zx</sub>: 0.13</td><td align="center" valign="middle" >G<sub>xz</sub>: 5050</td></tr></tbody></table></table-wrap><p>The damage and fracture are not considered during the analysis of obtained results. Consequently, the applied nominal tensile or compressive stress is limited to 70% of the material strength, i.e. is supposed to be σ x x n o m = &#177; 35 0   MPa .</p><p>In order to reduce the computation time, only a quarter or a half of the sampleis modeled depending on the symmetries of the configuration studied, namely:</p><p>- A quarter of the sample is modeled for configurations 1.1, 4.1 to 4.3 and 5.1.</p><p>- A half of the sample is modeled for configurations 2.1, 3.1 and 3.2, 5.2 and 5.3.</p></sec><sec id="s4_4"><title>3.4. Clamping Force</title><p>Every bolted joint, which has to transmit forces between components, can only work with a reliable and sufficient preload or clamping force. Several studies showed that clamping force can affect failure mode and ultimate strength of a bolt joint [ [<xref ref-type="bibr" rid="scirp.132895-ref19">19</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref20">20</xref>] for instance]</p><p>For this reason, a clamping force is taken into account in this research in order to evaluate the stress level in the composite material of the bolted joint. It is generated by a pre-stress in the bolt.</p><p>The pre-stress is usually taken equal to 90 percent of the yield stress of the bolt material. In the current study, the bolt class A2-80 was considered with a yield stress equal to 600 MPa and anultimate stress equal to 800 MPa. Accordingly, the maximal pre-stress in the bolt is then 540 MPa.</p><p>In order to generate the pre-stress, a thermal loading step is used. To simulatea tightening stress on the bolted joint, a negative uniaxial thermal strainis engendered in bolt using the linear thermal expansion of the material:</p><p>ε &#175; &#175; t h = α &#175; &#175;   Δ T (1)</p><p>where ΔT is the negative increment of temperature and</p><p>α &#175; &#175; is a tensor of expansion coefficients such that only the component α 33 is non nil:</p><p>α &#175; &#175; = ( 0 0 0 0 0 0 0 0 α 33 ) (2)</p><p>Such a tensor enables avoiding undesired thermal expansion in transverse directions of the bolt studied.</p><p>The value of ΔT is chosen to obtain the desired pre-stress in the bolt. ΔT = −277 K was necessary to obtain a pre-stress equal to 540 MPa using:</p><p>α 33 = 10 − 5 / K .</p></sec><sec id="s4_5"><title>3.5. Contacts</title><p>The bolted joint is constituted of two types of materials, metallic and composite. Some contacts are managed between the different parts and the materials.</p><p>For a steel/steel contact, a friction coefficient equal to 0.1 has been used. For a steel/glass-epoxy composite, a friction coefficient equal to 0.6 has been applied [<xref ref-type="bibr" rid="scirp.132895-ref21">21</xref>] .</p></sec></sec><sec id="s5"><title>4. Stress Analysis</title><sec id="s5_1"><title>4.1. Stress Discontinuities in the Components of the Bolted Joint</title><p>The stresses are analyzed in the “individual” parts of the bolted joints. The woven compositematerial of the bolted joint described in sub-sections 2.1 to 2.4 is composed of ten plies of same thickness that are supposed to be perfectly bonded one to each other. In the case of a hybrid bolted joint described in sub-sections 2.5 and 2.6, the adjacent layers of stainless steel plies and composite plies are also supposed to be perfectly bonded. Due to the misorientation of the plies constituting the woven glass-epoxy composite material or the difference of elastic properties of the glass-epoxy composite and the stainless steel, some stress tensor components are not necessarily continuous through the interface between two consecutive layers.</p><p>In fact, from a mechanical point of view, the stress vector T → should be continuous through each interface:</p><p>T → = σ &#175; &#175; n → (3)</p><p>where σ &#175; &#175; is the stress tensor and n → the unit vector normal to the interface.</p><p>Let’s consider two adjacent layers denoted (1) and (2) respectively, composed of two different materials, the z direction being perpendicular to the interface separating both layers.</p><p>Since T → ( 1 ) = T → ( 2 ) , one obtains:</p><p>⇒ σ &#175; &#175; ( 1 ) n → = σ &#175; &#175; ( 2 ) n → (4)</p><p>where σ &#175; &#175; ( 1 ) and σ &#175; &#175; ( 2 ) are the stress tensors in materials (1) and (2) respectively and</p><p>n → = { 0 0 1 } (5)</p><p>Consequently, it comes that:</p><p>σ x z ( 1 ) = σ x z ( 2 ) , σ y z ( 1 ) = σ y z ( 2 ) and σ z z ( 1 ) = σ z z ( 2 ) (6)</p><p>It appears that the three components calculated above are continuous while the components σ x x , σ y y and σ x y can be discontinuous.</p><p>In finite element method, stress components are calculated at the Gauss’ points of each element. The number and location of these points depend on the element type used. For hexahedral 8 node-elements, in the case of full integration technology, the stress and strain tensors are evaluated in 8 Gauss’ points. In order to correctly evaluate the stress state in each layer, taking into account the possible discontinuities, the stress state of the given integration point is translated to the nearest node of each element.</p></sec><sec id="s5_2"><title>4.2. Stress Concentration Factors</title><p>The analysis is focused on stress state in the composite part of the bolted joint.</p><p>It is obvious that stresses are increased in the vicinity of the hole when the bolted joint is loaded under tension or compression [<xref ref-type="bibr" rid="scirp.132895-ref22">22</xref>] . Stresses can also be increased in the area of the interface of the composite material and the stainless steel due to the change of material stiffness (see sub-section 4.1).</p><p>The simplest way to evaluate the local stress increase, and to compare the behavior of the different bolted joints described in section 2, is the use of stress concentration factor concept.</p><p>For a monotonic loading, two different stress concentration factors that describe increase of σ x x stress component in the composite material are defined:</p><p>The stress concentration factors for a tension loading denoted K<sub>t</sub> is defined as follows:</p><p>K t = σ x x m a x σ x x n o m where σ x x n o m is equal to 350 MPa (7)</p><p>the stress concentration factors for a compression loading denoted K<sub>t</sub>' is defined as follows:</p><p>K ′ t = σ x x m i n σ x x n o m where σ x x n o m is equal to −350 MPa (8)</p><p>As composite materials can be sensitive to fatigue, the stress concentration factor for a cyclic loading with a ratio R = −1, denoted K Δ σ is defined as follows:</p><p>K Δ σ = σ x x m a x − σ x x m i n 2 σ x x n o m where 2 σ x x n o m is equal to 700 MPa (9)</p></sec></sec><sec id="s6"><title>5. Results and Discussion</title><p>Stress concentration factors K t , K ′ t ’ and K Δ σ in the composite part are calculated using finite element method for each configuration described at section 2. Analyzes take into account the location of maximal stresses which can be the vicinity of the hole in the composite, the vicinity of the circumference of the metallic insert or the vicinity of the end of the biggest metallic ply.</p><p>Figures 12-16 give for each configuration the color maps of σ x x stress component and the values of stress concentration factors K t and K ′ t . The location of the most stressed point is also identified in these figures. <xref ref-type="table" rid="table2">Table 2</xref> and <xref ref-type="table" rid="table3">Table 3</xref> compare the values of stress concentration factors for all configurations.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>1 shows an example of color map and the zones denoted (a) et (b) which can be the most stressed ones depending on the configuration studied.</p><sec id="s6_1"><title>5.1. Simple Bolted Joint</title><p>Obvious results are obtained for the simple bolted joint (configuration 1.1). High concentration factors are identified which are located at the circumference of the hole, in the (x, z) symmetry plane.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Concentration factors for configurations 1.1 to 3.2</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >configuration</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x46.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x47.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x48.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x49.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x50.png" xlink:type="simple"/></inline-formula></th></tr></thead><tr><td align="center" valign="middle" >1.1</td><td align="center" valign="middle" >2.1</td><td align="center" valign="middle" >2.2</td><td align="center" valign="middle" >3.1</td><td align="center" valign="middle" >3.2</td></tr><tr><td align="center" valign="middle" >K<sub>t</sub></td><td align="center" valign="middle" >12.6</td><td align="center" valign="middle" >3.7</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >5.7</td><td align="center" valign="middle" >5.8</td></tr><tr><td align="center" valign="middle" >K<sub>t</sub><sub>’</sub></td><td align="center" valign="middle" >9.7</td><td align="center" valign="middle" >3.6</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >2.7</td><td align="center" valign="middle" >4.1</td></tr><tr><td align="center" valign="middle" >K<sub>Δσ</sub></td><td align="center" valign="middle" >8.9</td><td align="center" valign="middle" >3.6</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >3.1</td><td align="center" valign="middle" >3.1</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Concentration factors for configurations 4.1 to 5.2</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >configuration</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x51.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x52.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x53.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x54.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="/html.scirp.org/file/2-8104435x55.png" xlink:type="simple"/></inline-formula></th></tr></thead><tr><td align="center" valign="middle" >4.1</td><td align="center" valign="middle" >4.2</td><td align="center" valign="middle" >4.3</td><td align="center" valign="middle" >5.1</td><td align="center" valign="middle" >5.2</td></tr><tr><td align="center" valign="middle" >K<sub>t</sub></td><td align="center" valign="middle" >3.6</td><td align="center" valign="middle" >1.3</td><td align="center" valign="middle" >1.3</td><td align="center" valign="middle" >1.3</td><td align="center" valign="middle" >1.3</td></tr><tr><td align="center" valign="middle" >K<sub>t</sub><sub>’</sub></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >1.4</td><td align="center" valign="middle" >1.4</td><td align="center" valign="middle" >1.3</td><td align="center" valign="middle" >1.3</td></tr><tr><td align="center" valign="middle" >K<sub>Δσ</sub></td><td align="center" valign="middle" >2.5</td><td align="center" valign="middle" >1.3</td><td align="center" valign="middle" >1.3</td><td align="center" valign="middle" >1.3</td><td align="center" valign="middle" >1.3</td></tr></tbody></table></table-wrap><p>Very large values of K t equal to 12.6 and K ′ t equal to 17.8 have been obtained, see <xref ref-type="fig" rid="fig1">Figure 1</xref>2. The cyclic concentration factor is close to 9. They are due to the localized transfer of load from the bolt to the composite material. Since there is no clearance between the bolt and the composite plate, the stress concentration factor in compression K ′ t is higher than K t . Such a technological solution is not adapted for the highly loaded joints.</p></sec><sec id="s6_2"><title>5.2. Bolted Joint with an Insert</title><p>Configurations 2.1: The maximal and minimal stresses are located at the circumference of the insert washer. The BigHead<sup>R</sup> insert improves the load transfer. Indeed, the stress concentration factors are significantly smaller than those of simple bolted joint. K t and K ′ t equals 3.7 and 3.6 respectively (<xref ref-type="fig" rid="fig1">Figure 1</xref>3). The cyclic concentration factor K Δ σ is also equal to 3.6. The use of BigHead<sup>R</sup> insert decreases all the concentration factors by nearly three compared to the previous solution. However, the discontinuous variation of the stiffness between the insert and the composite plate still generates the local stress concentration well distinguishable in <xref ref-type="fig" rid="fig1">Figure 1</xref>3.</p><p>Configurations 2.2: In order to obtain a gradual transition between the different parts of the bolted joint, an insert with steps is used around the bolt (<xref ref-type="fig" rid="fig1">Figure 1</xref>4). As in the case of BigHead<sup>R</sup> insert, maximal and minimal stresses are located at the circumference of the insert, in the vicinity of its greatest radius. Stress concentration factors are not so high as for BigHead<sup>R</sup> Insert. K t , K ′ t and K Δ σ are equal to 2 (instead 3.7 and 3.6, see last paragraph).</p></sec><sec id="s6_3"><title>5.3. Waved Bolted Joint</title><p>Configuration 3.1: Theamplitude of the wave is 1 mm and the maximal stress occurs at the vicinity of the hole, in the same way as in the simple bolted joint. This configuration is an improvement of configuration 1.1, but the stress concentration factors are relative high: K t is equal to 5 for a tension loading and the value of K Δ σ is close to 3 (<xref ref-type="fig" rid="fig1">Figure 1</xref>5(a)). In fact, the wave effect is here not sufficient to reduce significantly the stress level in the composite.</p><p>Configuration 3.2: The amplitude of the wave is 2 mm and the stress concentration factors are quite similar to those of configuration 3.1. Nevertheless, the maximal stress is not located at the vicinity of the hole but at the top of the first half-period of the wave (<xref ref-type="fig" rid="fig1">Figure 1</xref>5(b)). This phenomenon is due to the too small radius of the wave.</p><p>Configurations 3.1 and 3.2 have stress concentration factors smaller than those of configuration 1.1 thanks to the waves but seem however to be not satisfactory.</p></sec><sec id="s6_4"><title>5.4. Hybrid Composite Bolted Joint</title><p>Configuration 4.1: A way to reduce the stress level and the stress concentration factor in the vicinity of the hole can be found by adding metallic plies in the composite material [<xref ref-type="bibr" rid="scirp.132895-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref15">15</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref16">16</xref>] [<xref ref-type="bibr" rid="scirp.132895-ref17">17</xref>] . The highest σ x x stress is again located at the vicinity of the hole, here in the first composite ply (<xref ref-type="fig" rid="fig1">Figure 1</xref>6(a)). In fact, from a geometrical point of view, configurations 1.1 and 4.1 are quite similar. Stress concentration factors are not so high for configuration 4.1 thanks to metallic plies that reinforce the composite. K t equals 3.6 and K Δ σ is equal to 2.5.</p><p>In the same way as for configuration 1.1, the absence of clearance between the bolt and the composite material induces K ′ t higher than K t . Configuration 4.1 is better compared to configuration 1.1, but still remains not satisfactory.</p><p>Configuration 4.2: Another way to reduce the stress level in the composite part is obtained by using a sleeve around the bold. In this case, the stress concentration factors are small, the values obtained for K t and K Δ σ are equal to 1.4, which is an acceptable level (<xref ref-type="fig" rid="fig1">Figure 1</xref>6(b)). It should be noted that thanks to the presence of the sleeve, maximal stress in tension and minimal stress in compression don’t occur at the vicinity of the hole because there isn’t any contact between the composite part and the bolt. There are bothlocalized at the end of the longest metallic ply.</p><p>Configuration 4.3: Configuration 4.2 can be modified by adding a second sleeve around the bolt. This corresponds to configuration 4.3. The stress state in the composite part is quite similar to those for configuration 4.2 and the stress concentration factors K t , K ′ t and K Δ σ have approximatively the same values.</p><p>They are at most equal to 1.4 for both configurations 4.2 and 4.3 which is an acceptable value. One can say that these two configurations are satisfactory.</p></sec><sec id="s6_5"><title>5.5. Hybrid Composite Bolted Joint with an Insert</title><p>Configurations 2.1 and 2.2 showed that an insert could improve a simple bolted joint by reducing stress concentration factors. Configuration 4.2 and 4.3 showed that the use of one or two sleeves improves the behavior of a simple hybrid bolted joint (configuration 4.1).</p><p>The idea is here to combine and to take advantage of these two solutions, insert and hybrid composite material, to improve the bolted joint.</p><p>For both configurations 5.1 and 5.2, stress in tension and stress in compression are maximal and minimal at the end of the longest composite ply. The stress concentration factors have the same values as those obtained for configurations 4.2 and 4.3. In fact, for these four configurations, the direct contact between the hybrid composite material and the bolt was avoided. The use of sleeves or stepped inserts does not change the stress state in the composite part: the values of stress concentration factors and the location of maximal stresses remained the same. In the same way as for configurations 4.2 and 4.3, configurations 5.1 and 5.2 are satisfactory.</p></sec></sec><sec id="s7"><title>6. Conclusions</title><p>Among the ten tested configurations, the 4.2, 4.3, 5.1 and 5.2 ones are acceptable with respect to the conventional values of stress concentration factors of less than 2. These configurations are made with hybrid composite bolted joints and sleeve or insert. In fact, thanks to these parts, there is no direct contact between composite material and metallic bolt. The maximal stresses don’t occur in the vicinity of the hole, but at the end of the longest metallic ply. As consequence, the stress concentration factors K t , K ′ t and K Δ σ are at most equal to 1.4. It should be noted that a low stress concentration factor in fatigue K Δ σ should be required for most applications such as aeronautic and automotive ones, because these structures are most of time submitted to cyclic loadings.</p><p>Overall, a bolted joint is intended to be easy and inexpensive to produce. The shape of insert of configurations 5.1 and 5.2 make them expensive to machine. For this reason, the solutions 5.1 and 5.2 are less interesting than the configurations 4.2 and 4.3.</p><p>The stress concentration factors for these two cases have about the same values, however the configuration 4.3 should be preferred since it is less loaded thanks to the shoulder.</p><p>In a next work, the configuration 4.3 will be experimentally tested at room, but also at cryogenic temperatures.</p><p>Indeed, stainless steel/glass-epoxy bolted joints such that of this study can be used as components of cryogenic devices such as fastener of liquefied gas tank or in the design of semi-trailers for transport of liquefied gas [<xref ref-type="bibr" rid="scirp.132895-ref23">23</xref>] .</p></sec><sec id="s8"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s9"><title>Cite this paper</title><p>Schmitt, C., Kremeur, A., Lipinski, P. and Capelle, J. (2024) Numerical Optimization by Finite Element Method of Stainless Steel/Glass-Epoxy Composite Bolted Joint under Tension and Compression. Engineering, 16, 102-122. https://doi.org/10.4236/eng.2024.164009</p></sec></body><back><ref-list><title>References</title><ref id="scirp.132895-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Lim, G.-H., Bodjona, K., Raju, K.P., &lt;i&gt;et al&lt;/i&gt;. (2018) Evolution of Mechanical Properties of Flexible Epoxy Adhesives under Cyclic Loading and Its Effects on Composite Hybrid Bolted/Bonded Joint Design. &lt;i&gt;Composite Structures&lt;/i&gt;, 189, 54-60. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2018.01.049</mixed-citation></ref><ref id="scirp.132895-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Lopez-Cruz, P., Lalibert&amp;#233;, J. and Lessard, L. (2017) Investigation of Bolted/Bonded Composite Joint Behaviour Using Design of Experiments. &lt;i&gt;Composite Structures&lt;/i&gt;,&lt;i&gt; &lt;/i&gt;170, 192-201. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2017.02.084</mixed-citation></ref><ref id="scirp.132895-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Khashaba, U.A. (2018) Static and Fatigue Analysis of Bolted/Bonded Joints Modified with CNTs in CFRP Composites under Hot, Cold and Room Temperatures. &lt;i&gt;Composite Structures&lt;/i&gt;, 194, 279-291. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2018.04.008</mixed-citation></ref><ref id="scirp.132895-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Marannano, G. and Zuccarello, B. (2015) Numerical Experimental Analysis of Hybrid Double Lap Aluminum-CFRP Joints. &lt;i&gt;Composites Part B&lt;/i&gt;:&lt;i&gt; Engineering&lt;/i&gt;, 71, 28-39. &lt;br&gt;https://doi.org/10.1016/j.compositesb.2014.11.025</mixed-citation></ref><ref id="scirp.132895-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Liu, F.R., Lu, X.H., Zhao, L.B., Zhang, J.Y., Xu, J.F. and Hu, N. (2018) Investigation of Bolt Load Redistribution and Its Effect on Failure Prediction in Double-Lap, Multi-Bolt Composite Joints.&lt;i&gt; &lt;/i&gt;&lt;i&gt;Composite Structures&lt;/i&gt;, 202, 397-405. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2018.02.043</mixed-citation></ref><ref id="scirp.132895-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Mandal, B. and Chakrabarti, A. (2018) Numerical Failure Assessment of Multi-Bolt FRP Composite Joints with Varying Sizes and Preloads of Bolts. &lt;i&gt;Composite Stru&lt;/i&gt;&lt;i&gt;c&lt;/i&gt;&lt;i&gt;tures&lt;/i&gt;, 187, 169-178. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2017.12.048</mixed-citation></ref><ref id="scirp.132895-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Nerilli, F. and Vairo, G. (2017) Progressive Damage in Composite Bolted Joints via a Computational Micromechanical Approach. &lt;i&gt;Composites Part B&lt;/i&gt;:&lt;i&gt; Engineering&lt;/i&gt;, 111, 357-371. &lt;br&gt;https://doi.org/10.1016/j.compositesb.2016.11.056</mixed-citation></ref><ref id="scirp.132895-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Hua, X.F., Haris, A., Ridha, M., Tan, V.B.C. and Tay, T.E. (2018) Progressive Failure of Bolted Single-Lap Joints of Woven Fiber-Reinforced Composites. &lt;i&gt;Composite Structures&lt;/i&gt;, 189, 443-454. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2018.01.104</mixed-citation></ref><ref id="scirp.132895-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Khashaba, U.A., Sallam, H.E.M., Al-Shorbagy, A.E. and Seif, M.A. (2006) Effect of Washer Size and Tightening Torque on the Performance of Bolted Joints in Composite Structures. &lt;i&gt;Composite Structures&lt;/i&gt;, 73, 310-317. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2005.02.004</mixed-citation></ref><ref id="scirp.132895-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Choia, J., Hashemini, S.M., Chun, H.-J., Park, J.-C. and Chang, H.S. (2018) Failure Load Prediction of Composite Bolted Joint with Clamping Force. &lt;i&gt;Composite Stru&lt;/i&gt;&lt;i&gt;c&lt;/i&gt;&lt;i&gt;tures&lt;/i&gt;, 189, 247-255. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2018.01.037</mixed-citation></ref><ref id="scirp.132895-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Liu, F.R., Lu, X.H., Zhao, L.B., Zhang, J.Y., &lt;i&gt;et al&lt;/i&gt;. (2018) An Interpretation of the Load Distributions in Highly Torqued Single-Lap Composite Bolted Joints with Bolt-Hole Clearances. &lt;i&gt;Composites Part B&lt;/i&gt;:&lt;i&gt; Engineering&lt;/i&gt;, 138, 194-205. &lt;br&gt;https://doi.org/10.1016/j.compositesb.2017.11.027</mixed-citation></ref><ref id="scirp.132895-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Muc, A. and Ulatowska, A. (2012) Local Fibre Reinforcement of Holes in Composite Multilayered Plates. &lt;i&gt;Composite Structures&lt;/i&gt;, 94, 1413-1419. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2011.11.017</mixed-citation></ref><ref id="scirp.132895-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Zhu, Y.D., Qin, Y.L., Qi, S.J., Xu, H.B., Liu, D. and Yan, C. (2018) Variable Angle Tow Reinforcement Design for Locally Reinforcing an Openhole Composite Plate. &lt;i&gt;Composite Structures&lt;/i&gt;, 202, 162-169. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2018.01.021</mixed-citation></ref><ref id="scirp.132895-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Kolesnikov, B., Herbeck, L. and Fink, A. (2008) CFRP/Titanium Hybrid Material for Improving Composite Bolted Joints. &lt;i&gt;Composite Structures&lt;/i&gt;, 83, 368-380. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2007.05.010</mixed-citation></ref><ref id="scirp.132895-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Camanho, P.P., Fink, A., Obst, A. and Pimenta, S. (2009) Hybrid Titanium-CFRP Laminates for High-Performance Bolted Joints. &lt;i&gt;Composites Part A&lt;/i&gt;:&lt;i&gt; Applied Science and Manufacturing&lt;/i&gt;, 40, 1826-1837. &lt;br&gt;https://doi.org/10.1016/j.compositesa.2009.02.010</mixed-citation></ref><ref id="scirp.132895-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Kolks, G. and Tserpes, K.I. (2014) Efficient Progressive Damage Modeling of Hybrid Composite/Titanium Bolted Joints. &lt;i&gt;Composites Part A&lt;/i&gt;:&lt;i&gt; Applied Science and Man&lt;/i&gt;&lt;i&gt;u&lt;/i&gt;&lt;i&gt;facturing&lt;/i&gt;, 56, 51-63. &lt;br&gt;https://doi.org/10.1016/j.compositesa.2013.09.011</mixed-citation></ref><ref id="scirp.132895-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Fink, A., Camanho, P.P., Andr&amp;#233;s, J.M., Pfeiffer, E. and Obst, A. (2010) Hybrid CFRP/ Titanium Bolted Joints: Performance Assessment and Application to a Spacecraft Payload Adaptor. &lt;i&gt;Composites Science and Technology&lt;/i&gt;, 70, 305-317. &lt;br&gt;https://doi.org/10.1016/j.compscitech.2009.11.002</mixed-citation></ref><ref id="scirp.132895-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">BigHead&amp;#174;. &lt;br&gt;https://www.bighead.co.uk </mixed-citation></ref><ref id="scirp.132895-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Sasikumar, A., Guerrero, J.M., Quintanas-Corominas, A. and Costa, J. (2021) Numerical Study to Understand Thermo-Mechanical Effects on a Composite Aluminium Hybrid Bolted Joint. &lt;i&gt;Composite Structures&lt;/i&gt;, 275, Article ID: 114396. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2021.114396</mixed-citation></ref><ref id="scirp.132895-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Choi, J., Hasheminia, S.M., Chun, H.-J., Park, J.-C. and Chang, H.S. (2018) Failure Load Prediction of Composite Bolted Joint with Clamping Force. &lt;i&gt;Composite Stru&lt;/i&gt;&lt;i&gt;c&lt;/i&gt;&lt;i&gt;tures&lt;/i&gt;, 189, 247-255. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2018.01.037</mixed-citation></ref><ref id="scirp.132895-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Larsen, T., &lt;i&gt;et al&lt;/i&gt;. (2007) Comparison of Friction and Wear for an Epoxy Resin Reinforced by a Glass or a Carbon/Aramid Hybrid Weave. &lt;i&gt;Wear&lt;/i&gt;, 262, 1013-1020. &lt;br&gt;https://doi.org/10.1016/j.wear.2006.10.004</mixed-citation></ref><ref id="scirp.132895-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Cabrera-Gonz&amp;#225;lez, J.A., Vargas-Silva, G. and Barroso, A. (2021) Riveted Joints in Composites, a Practical Tool to Estimate Stresses around the Rivet Hole. &lt;i&gt;Composite Structures&lt;/i&gt;, 263, Article ID: 113735. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2021.113735</mixed-citation></ref><ref id="scirp.132895-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Lee, K.H. and Lee, D.G. (2010) Design of a Hybrid Glass Composite Anchor for LNG Cargo Containment Systems. &lt;i&gt;Composite Structures&lt;/i&gt;, 92, 469-479. &lt;br&gt;https://doi.org/10.1016/j.compstruct.2009.08.031</mixed-citation></ref></ref-list></back></article>