A Numerical Analysis of Riveted Lap Joint Containing Multiple-site Damage

Similar documents
An experimental investigation on crack paths and fatigue behaviour of riveted lap joints in aircraft fuselage

Subject Index. Bearing damage, in bolted, composite-metal

AMTS STANDARD WORKSHOP PRACTICE. Bond Design

SIMULATION AND EXPERIMENTAL WORK OF SINGLE LAP BOLTED JOINT TESTED IN BENDING

AN INNOVATIVE FEA METHODOLOGY FOR MODELING FASTENERS

EFFECT OF ADHEREND SHAPE ON STRESS CONCENTRATION REDUCTION OF ADHESIVELY BONDED SINGLE LAP JOINT

Finite Element Analysis of Multi-Fastened Bolted Joint Connecting Composite Components in Aircraft Structures

Disclaimer for FAA Research Publication

INFLUENCE OF SOME MODIFICATIONS OF LOCAL GEOMETRY ON THE STRESS STATES IN ADHESIVE BONDED LAP JOINTS

Structural Strength of Lapped Cold-Formed Steel Z-Shaped Purlin Connections with Vertical Slotted Holes

STRESS ANALYSIS OF RIVETED LAP JOINT USING FINITE ELEMENT METHOD

EXPERIMENTAL INVESTIGATION OF FATIGUE BEHAVIOUR IN COMPOSITE BOLTED JOINTS

Fastener Modeling for Joining Parts Modeled by Shell and Solid Elements

INFLUENCE OF PILES ON LOAD- SETTLEMENT BEHAVIOUR OF RAFT FOUNDATION

Maximizing the Fatigue Crack Response in Surface Eddy Current Inspections of Aircraft Structures

Failure of Engineering Materials & Structures. Code 34. Bolted Joint s Relaxation Behavior: A FEA Study. Muhammad Abid and Saad Hussain

Finite Element Modeling of Early Stage Self-loosening of Bolted Joints Haoliang Xu 1, a, Lihua Yang 1, b,, Lie Yu 1,2, c

Study of an interference fit fastener assembly

An experimental study on fatigue crack growth in lap joints with multiple fastener holes

THE INFLUENCE OF GEOMETRIC PARAMETERS AND MECHANICAL PROPERTIES OF ADHESIVE ON STRESS ANALYSIS IN ADHESIVELY BONDED ALUMINUM SINGLE LAP JOINT

Glued laminated timber beams repair.

ANALYSIS OF LATERAL STIFFNESS FOR INFILLED FRAME WITH OPENING

3-D Finite Element Analysis of Bolted Joint Using Helical Thread Model

1. Enumerate the most commonly used engineering materials and state some important properties and their engineering applications.

Effect of Bolt Layout on the Mechanical Behavior of Four Bolted Shear Joint

UC San Diego UC San Diego Electronic Theses and Dissertations

A Numerical study on proper mode and frequency selection for riveted lap joints inspection using Lamb waves.

An Alternative Formulation for Determining Stiffness of Members with Bolted Connections

LS-DYNA USED TO ANALYZE THE MANUFACTURING OF THIN WALLED CANS AUTHOR: CORRESPONDENCE: ABSTRACT

Modeling Multi-Bolted Systems

Stress Analysis of Flanged Joint Using Finite Element Method

Two basic types of single

The influence of gouge defects on failure pressure of steel pipes

Moment-Resisting Connections In Laminated Veneer Lumber (LVL) Frames

Module 10 : Improvement of rock mass responses. Content

Experimental Evaluation of Metal Composite Multi Bolt Radial Joint on Laminate Level, under uni Axial Tensile Loading

ON ACHIEVING AN OPTIMAL RIVETED LAP JOINT DESIGN FOR FIBRE METAL LAMINATES

Shear Stress Analysis of Single Chain Riveted Lap Joint

DESIGN OPTIMISATION OF 3D WOVEN T-JOINT REINFORCEMENTS

An Investigation of Optimal Pitch Selection to Reduce Self-Loosening of Threaded Fastener under Transverse Loading

Design, Fabrication, and Validation of an Ultra-Lightweight Membrane Mirror

Design of a Piezoelectric-based Structural Health Monitoring System for Damage Detection in Composite Materials

Model Correlation of Dynamic Non-linear Bearing Behavior in a Generator

Dynamics of High-speed Machining of Aerospace Structures using Finite-element Analysis

A training course delivered at a company s facility by Matrix Engineering, an approved provider of Bolt Science Training

Nonlinear behavior of Reinforced Concrete Infilled Frames using ATENA 2D

Mechanical behavior of fiberglass reinforced timber joints

Strain Measurements with the Digital Image Correlation System Vic-2D

Heat-Mechanics Interaction Behavior of Laminated Rubber Bearings under Large and Cyclic Lateral Deformation

In-Situ Damage Detection of Composites Structures using Lamb Wave Methods

DETECTION AND SIZING OF SHORT FATIGUE CRACKS EMANATING FROM RIVET HOLES O. Kwon 1 and J.C. Kim 1 1 Inha University, Inchon, Korea

A training course delivered at a company s facility by Matrix Engineering, an approved provider of Bolt Science Training

EXPERIMENTAL INVESTIGATION ON LASER BENDING OF METAL SHEETS USING PARABOLIC IRRADIATIONS

Dowel. Design. Performance-Based World of Concrete Official Show Issue. Lift-truck design changes require a new look at joint durability

Slip and yield resistance of friction type of high strength bolted connections with over-sized holes

EXPERIMENTAL ANALYSIS OF BOLT LOOSENING DYNAMICS CHARACTERISTIC IN A BEAM BY IMPACT TESTING

DESIGN, ANALYSIS AND OPTIMIZATION OF CURVE ATTACHMENT ON COMPOSITE HYBRID LAP JOINT

The Behaviour Of Round Timber Sections Notched Over The Support On The Tension Face. Justin Dewey

Indian Journal of Engineering An International Journal ISSN EISSN Discovery Publication. All Rights Reserved

DESIGN OF MACHINE MEMBERS-I

A Study on Effect of Sizing Bolt Hole in Single-Lap Connection Using FEA

What happens to bolt tension in large joints?, Fasteners, Vol. 20, No. 3, 1965, Publication No. 286

FINITE ELEMENT SIMULATIONS OF THE EFFECT OF FRICTION COEFFICIENT IN FRETTING WEAR

DEVELOPMENT OF A NOVEL TOOL FOR SHEET METAL SPINNING OPERATION

Design and Evaluation of a Bolted Joint for a Discrete Carbon-Epoxy Rod-Reinforced Hat Section

Quantitative Crack Depth Study in Homogeneous Plates Using Simulated Lamb Waves.

Preliminary Design and Analysis of an In-plane PRSEUS Joint

Anti-check bolts as means of repair for damaged split ring connections

MANUFACTURING TECHNOLOGY

THE INFLUENCE OF GOUGE DEFECTS ON FAILURE PRESSURE OF STEEL PIPES

MODELLING AND EXPERIMENTS FOR THE DEVELOPMENT OF A GUIDED WAVE LIQUID LEVEL SENSOR

EFFECT OF YARN CROSS-SECTIONAL SHAPES AND CRIMP ON THE MECHANICAL PROPERTIES OF 3D WOVEN COMPOSITES

VIBRATIONAL MODES OF THICK CYLINDERS OF FINITE LENGTH

Keywords: Bracing bracket connection, local deformation, selective pallet racks, shear stiffness, spine bracings.

A finite element stress analysis of aircraft bolted joints loaded in tension

NUMERICAL ANALYSIS OF SCREW ANCHOR FOR CONCRETE

LABORATORY EVALUATION OF CONCRETE-FILLED GFRP DOWELS IN JOINTED CONCRETE PAVEMENTS

Structural Bolting. Notice the Grade 5 has a much smaller head configuration and a shorter shank then the grade A325 structural bolt.

Dowel connections in laminated strand lumber

EDDY CURRENT INSPECTION FOR DEEP CRACK DETECTION AROUND FASTENER HOLES IN AIRPLANE MULTI-LAYERED STRUCTURES

STRUCTURAL TIMBER DESIGN

NON-LINEAR CONNECTION MODELS IN TIMBER ENGINEERING

Experimental and numerical investigation of tube sinking of rectangular tubes from round section

Fatigue and Fretting Studies of Gas Compressor Blade Roots

Wire and tube Drawing

1/2/2016. Lecture Slides. Screws, Fasteners, and the Design of Nonpermanent Joints. Reasons for Non-permanent Fasteners

Module 2 WAVE PROPAGATION (Lectures 7 to 9)

Riveted Lap Joints in Aircraft Fuselage

Welded connections Welded connections are basically the same design in AISI as in AISC. Minor differences are present and outlined below.

EXPERIMENTAL INVESTIGATION OF ADHESIVE STRENGTHS OF ADHESIVELY BONDED JOINTS

Enhancing the low frequency vibration reduction performance of plates with embedded Acoustic Black Holes

Corso di Studi di Fabbricazione

THIN-WALLED HOLLOW BOLTS

Behaviour of fibre reinforced composite beams with mechanical joints

WITHDRAWAL AND LATERAL STRENGTH OF THREADED NAILS

Hanger bolts and solar fasteners in sandwich panels

NOTICE. The above identified patent application is available for licensing. Requests for information should be addressed to:

Experimental And FE Analysis Of Eccentric Loaded Symmetrical And Unsymmetrical Bolted Joint With Bolt Pretension

Influence of Lubrication and Draw Bead in Hemispherical Cup Forming

Tension Perpendicular to Grain Strength of Wood, Laminated Veneer Lumber (LVL), and Cross-Banded LVL (LVL-C)

Transcription:

Appl. Math. Inf. Sci. 7, No. 2L, 699-704 (2013) 699 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/072l46 A Numerical Analysis of Riveted Lap Joint Containing Multiple-site Damage Dazhao YU Naval Aeronautical and Astronautical University, Yantai, Shandong, 264001, P. R. China Received: 7 Oct. 2012, Revised: 26 Jan. 2013, Accepted: 28 Jan. 2013 Published online: 1 Jun. 2013 Abstract: The evaluation of stress intensity factors (SIFs) for cracks in mechanically fastened joints is one of the central issues in a damage tolerance analysis. Accurate stress intensity solutions may be difficult to determine due to geometric complexity along with variations in fastener load transfer and fastener interference and material thickness because of corrosion. Detailed finite element models that include specific aspects of fastener for lap joints analysis, but such representations are often impractical for large lap joints involving many fasteners. In order to reduce the number of degrees of freedom in a particular model, a methodology was implemented that efficiently depicts mechanical fasteners in lap joints using finite elements. Uncorroded and corroded lap joints with three crack scenarios were studied by the computationally efficient model. The effects of pillowing corrosion and fastener interference were also included in the model. The results show that the effect of MSD and thinning of the material is to increase substantially stress intensity factors (SIFs) values compare to that of a single crack without corrosion. For a given cyclic stress range, SIF decreases with increasing rivet interference level. This is particularly true for shorter crack lengths. SIF values for the longer crack are not sensitive to the length of shorter crack on the opposite side of the rivet hole. The effects of corrosion pillowing cause SIF values on the faying surface to be larger than those on the opposing surface, and the crack front distorts from the straight front. Moreover, stresses on the faying surface may exceed the yield strength of the material. Keywords: lap joint, stress intensity factor, multiple site damage, fastener interference 1 Introduction Multiple site damage (MSD) cracks have been of great concern since the Aloha accident in April 1988 [1-7]. The coalescence of MSD cracks may lead to sudden fracture and adversely affect the damage tolerance of structures. This is because the anticipated failure load will be significantly underestimated if consideration is taken solely of the fracture resistance of a single crack. MSD is widespread fatigue damage with the unusual feature that small fatigue cracks appear at multiple rivet locations along the lap joint. One central issue in quantitatively evaluate the effect of MSD cracks on the structure integrity is evaluating stress intensity factors (SIFs) for MSD cracks at rivet holes in mechanically riveted joints. Accurate stress intensity solutions may be difficult to determine due to geometric complexity along with variations in rivet load transfer and rivet interference and material thickness because of corrosion, especially the interaction of MSD cracks. Detailed finite element models of the structure and associated rivets can be used to estimate the state of stress in lap joint. However, it is often impractical to explicitly model specific details of rivets geometry for a significant number of rivets. In order to reduce the number of degrees of freedom (DOFs) in a particular model, individual rivets connecting adjacent layered structures may be idealized using simple spring elements [3]. Of course, the adopted spring force-displacement relationship has a substantial effect on the load transfer in a given joint as well as on the SIF solutions for any cracks in the joint [8]. Cope and Lacy verified this approach by considering different combinations of explicitly modeled fasteners and spring-rivet fasteners in single shear lap joint models with three rows of fasteners (cf. Figure.1.1 for three such models). Based on previous investigation, this study investigates the accuracy of the efficient modeling methods to determine SIFs for riveted lap joints with MSD. The influence of rivet interference and corrosion damage on SIFs solutions was also investigated. Corresponding author e-mail: ytyudazhao@yahoo.com.cn

700 Dazhao YU: A Numerical Analysis of Riveted Lap Joint Containing Multiple-site Damage t 1 Adhesive elements difined inside circular mesh consisted of double cracks of equal length emanating from a single rivet near the center of the specimen. The third scenario consisted of double cracks of differing lengths emanating from a single rivet near the center the specimen. Far filed Stress=88.8MPa 1mm Side Straps 1mm Outer Skin t 2 Rivet(spring) element difined at sttsching node Corroded Region 76 Rivet Dia.=4.9mm 1 2 3 4 5 6 7 8 Fig. 1 Lap Joint Models with Various Rivet Model Combinations 12.7 Spring-rivet Elements 1mm Inner Skin 2 Analysis Approach and Tool Carleton University and National Research Council of Canada (NRC) have completed several fatigue tests on uncorroded and corroded lap joint specimens. Figure.2.1 shows the geometry and material [10]. Test specimens have been specially designed to simulate aircraft conditions, in particular the development of MSD and the MSD failure model. The key feature of the concept is the use of bonded side straps to simulate the load transfer from cracked areas to the surrounding structure that occurs on aircraft. This allows MSD cracks to develop and linkup premature failure of the specimen. The specimens were constructed of two 2024-T3 clad aluminum sheets with a thickness of 1mm and a width mm. Sheets were connected to a doubler with a thickness of 4mm. The edges of the spliced sheets were reinforced by straps that were bonded together. In order to verify that if Cope s method suite for lap joint with MSD, a stress analysis was performed to determine crack tip SIF values of MSD lap joint specimens (cf., Figure.2.1). In analyzing the crack growth test data, it was observed that 95% or more of the joint life was controlled by the lead crack behavior. 50% of the occurrences were a single crack (Scenario 1), 10% were double cracks with equal length (Scenario 2), and the remaining 40% were double cracks with different lengths on each side of a rivet (Scenario 3). The thickness loss estimated in the corroded specimens was about 2%. Analysis was only preformed for the lead crack using three different scenarios, as observed in the test data. The first scenario consisted of a single, radial crack emanating from a single rivet near the center of the specimen. The second scenario Fig. 2 NRC Lap Joint Specimen Configuration (unit: mm) The special purpose code FRANC2D/L[11] was used in this study to determine crack tip SIF values of MSD lap joint specimens. FRANC2D/L is a highly interactive finite element program well suited to performing small deformation analyses of two-dimensional layered structures. The mode allows for efficient SIF determination and incorporates an adaptive remeshing algorithm that facilitates mixed-mode crack propagation and incorporates an adaptive remeshing algorithm that facilitates mix-mode crack propagation. The layered capability allows the user to model mechanically fastened and/or adhesively bonded structures, such as planar lap joints and bonded repairs. Individual structural components such as skins, stiffeners, and doublers may be modeled as separate layers in FRANC2D /L, each with its own associated finite element mesh. Overlapping layers may be defined with coincident mid-planes or with an appropriate mid-plane offset to account for joint eccentricity. FRANC2D/L uses a linear bending algorithm to automatically calculate the bending stresses/strains and out-of-plane displacements in addition to the standard in-plane stress, strain, and distributions. 3 FRANC2D/L Model Cope and Lacy modelled single shear lap joints with three rows of mechanical fasteners using different combinations of explicit and spring element

Appl. Math. Inf. Sci. 7, No. 2L, 699-704 (2013) / www.naturalspublishing.com/journals.asp 701 representations of fasteners in order to determine the appropriate level of model refinement necessary foe accurate SIF solutions for cracks emanating from fasteners holes. The results suggest that combinations of explicit and spring element representations of fasteners may be used to develop computationally efficient lap joint models while retaining the essential features of the fastener displacement and load transfer within the joint. In order to reduce the number of DOFs in lap joint model, individual rivets connecting adjacent layered structures were used simple spring element in this study. Corroded Region 76 Rivet Dia.=4.9mm 12.7 Far filed Stress=88.8MPa 1mm Side Straps 1mm Outer Skin 1 2 3 4 5 6 7 8 Spring-rivet Elements 1mm Inner Skin Fig. 3 Mesh of Outer Layer used for FRANC2D/L SIF Calculations The FRANC2D/L model configuration (cf., Figure.3.1) was based on the geometry, material, and loading of the experimental specimens. The model consisted of six-noded triangular plate elements. The inner and outer skins, modelled as discrete layers, were attached only at the rivet locations and along the side straps. Material properties for the adhesive elements and for the spring-rivet elements were selected in order to provide reasonable estimates of the relative rivet displacements and load transfer between the layers. When modelling mechanically riveted lap shear lap joints using a combination of explicit and simplified spring representations of rivets, use of equations(1)-(4) to define the rivet properties ensures that the simplified spring elements will provide approximately the same load transfer and relative displacement between sheets as explicitly modelled rivets. [ E F = 8Eh (1 + ν) + 8 πd 3 ( ) ] h 2 [ ( d B +C + d )] 1 d t 1 t 2 (1) G F = E F 2(1 + ν) k s = πd2 G F 4h (2) (3) k b = 3πd4 E F 64h 3 (4) where E = sheet material elasticity modulus; h = joint eccentricity; ν = Poisson s ratio; t 1 = thickness of sheet in Layer1; t 2 = thickness of sheet in Layer2; k s = shear stiffness; k b = bending stiffness. B = 5.0, C = 0.8, for aluminium rivets. Rivets were modeled using two methods. In the first method, the rivets were modeled as circular, elastic plugs in both layers, and the plugs are attached to the surrounding structure using gap elements. The second method of modeling rivets were modeled by one-dimensional elastic-plastic rivet elements that connect shell finite element nodes in the upper and lower skins. Each rivet element has six degrees of freedom, corresponding to extension, shearing, bending and twisting of the rivet. The stiffness of each degree of freedom was defined by prescribing a force-deflection curve. The axial, flexural, and tensional stiff nesses were computed by assuming that the rivet behaves like a simple elastic rod with a diameter of 4.9mm.The elastic shear stiffness of the rivet element was computed by the equations (3)-(4). Contact interface elements were applied to the rivets in the upper row. Load eccentricity inherent to the single lap joint problem was introduced by offsetting the two skin layers. The effect of material thinning was modeled through a uniform thickness reduction over the corroded region. In addition to the effects of MSD, material thinning due to corrosion damage was also studied. The effect of material thinning was modeled by a uniform reduction in thickness of the upper skin in corroded region. Three levels of corrosion were considered, no corrosion and corrosion equivalent to 2% and 5% of the material thickness. The effect of the transverse pressure to model pillowing was to induce high local bending stresses in the skins around each rivet, thereby producing non-uniform stresses through the thickness of the skin. In order to investigate the effect of rivet interference on the predicted load transfer and SIFs, rivet interference were also modeled by adjusting the stress-displacement curve in the interface elements to produce a residual circumferential compressive stress around the rivet hole. Cracks were introduced into the outer skin in upper rivet row according to three different scenarios.

702 Dazhao YU: A Numerical Analysis of Riveted Lap Joint Containing Multiple-site Damage 4 FRANC2D/L Analysis Results Results from the FRANC2D/L analysis of the test specimen revealed that a load distribution of 35% in the lead trailing rows and 30% in the center rivet row. The results matched strain gage data reported in the experimental study [8]. This suggests that large lap joint models may be constructed using both explicit and spring element representations of rivets without a serious loss of accuracy. The K-solutions obtained from the FRANC2D/L analysis were recorded and converted to beta factors. The conversion equation is given below: β = K max K res σ re f πc (5) Where K is the stress intensity factor, σ re f is the reference stress based on experimental load level [8], and c is the crack length. To evaluate the effect of MSD cracks, the cracks of three different scenarios were introduced into the outer skin in upper rivet row. Figure.4.1 compares the beta results for the three crack configurations corresponding to Scenarios 1 through 3. The plot shows a significant increase in beta for the case of diametrically opposed cracks compared to that of a single crack, as has been observed in test studies. With increasing, this phenomenon is more apparent. Results for Scenarios2 and 3 show very similar beta values, which indicate that values for the longer crack in Scenario 3 are not sensitive to the length of shorter crack on the opposite side of the hole. Scenario 2 provides an upper bound to beta values. Since Scenario 2 and 3 yielded such close results, Scenario 3 was omitted from most of the subsequent analysis. Without Corrosion and Rivet Interference Values Recorded for crack at H5L 0.0 0.5 Crack Scenario 1 Crack Scenario 2 Crack Scenario 3 Fig. 4 Predictions for Different Crack Scenarios It has been shown that rivet interference significantly impacts crack growth by reducing the initial crack growth rates [3, 9-10]. In order to investigate the effect of rivet interference on SIFs, four different levels of interference (0.00%, 0.25%, 0.5%, 0.78%) were studied. The betas were calculated to assess the effect of various interferences on beta, as shown in Figure.4.2. As expected, for a given, beta decreases with increasing rivet interference level. This is particular true for shorter crack length. The effect of corrosion was modeled through uniform thinning of the material in the corroded region. In this case, one would expect the average stresses in the corroded region, and therefore the beta values, to increase by a factor of approximately t/t corr, where t is the original thickness, and t corr is the corroded panel thickness. The FRANC2D/L prediction yields the expected result, as the beta values for the 2% and 5% corrosion case are approximately 2 and 5 times higher than for the case with no corrosion respectively (cf., Figure.4.3). The results were consistently obtained for all analysis cases which had zero rivet interference. MSD Cracking Scenario 2 Without Corrosion 0.0 0.5 0.00% Interference 0.25% Interference 0.50% Interference 0.78% Interference Fig. 5 Predictions for Different Levels of Rivet Interference 4.5 MSD Cracking Scenario 2 Without Interference 0.0 0.5 No Corrosion 2% Corrosion 5% Corrosion Fig. 6 Predictions for Different Corrosion Level The results also show that the effects of corrosion pillowing cause SIF values on the faying surface to be larger than those on the opposing surface, and the crack front distorts from the straight front. Moreover, stresses

Appl. Math. Inf. Sci. 7, No. 2L, 699-704 (2013) / www.naturalspublishing.com/journals.asp 703 on the faying surface may exceed the yield strength of the material, producing plastic deformation. It should be recognized that linear two-dimensional finite element analysis of riveted lap joints may be ill suited for the case where the out-of-plane displacements become large relative to the sheet thickness. A full three-dimensional analysis may be necessary in order to appropriately account for the influence of bending effects due to joint eccentricity on the load transfer and calculated peak SIFs. This key issue, as well as development of idealized rivet reorientations for use three-dimensional finite element models, will be the focus of further investigations. The last scenario investigated combined all the previously discussed effects: MSD crack, and corrosion. Each of these effects was modeled using the methods previously discussed. Figure.4.3 compares beta for MSD, corrosion, and combined scenario. It is evident that MSD and corrosion damage especial combined scenario significantly increase SIF values compare to that of a single crack. Values Recorded for crack at H5 Single Crack MSD Only 2% Corrosion Only All Effects local stress field is required, such as in the vicinity of cracks emanating from rivet holes. 2. The effect of MSD and corrosion is to increase substantially beta values compare to that of a single crack without corrosion. 3. For a given cyclic stress range, the beta decreases with increase rivet interference level. This is particularly true for shorter crack lengths. When the rivet interference reaches special value, its effect on SIF will not greatly change with interference increasing. 4. values for the longer crack are not sensitive to the length of shorter crack on the opposite side of the hole. 5. The effects of corrosion pillowing cause beta values on the faying surface to be larger than those on the opposing surface, and the crack front distorts from the straight front. Moreover, stresses on the faying surface exceed the yield strength of the material. In this paper, we analyzed WDM optical communication network structure, discussed the add-drop system and proposed a new routing selection method. Routing technology plays an important role in a optical communication network. We discussed the routing design and switching technology for WDM network in optical networks. Routing is the process of selecting paths in a network along which to send network traffic. A good routing selection method can make the network do better performance. Acknowledgements 0.0 0.5 The present work is part of the contribution of National Natural Science Foundation of China(No. 50675221). Thanks for the help. Fig. 7 Combined Scenario: MSD and Corrosion 5 Conclusion A methodology for efficiently depicting rivets in lap joints using finite elements was implemented in order to facilitate numerical analysis of lap joint with MSD. FRANC2D/L model of MSD lap joint was modeled using a combination of explicitly modeled rivets and spring element representations of rivets. SIF solutions were obtained for three scenarios MSD. These results suggest: 1. Combinations of explicit and spring element representations of rivets can be used to develop computationally efficient lap joint models while retaining the essential features of the deformation and load transfer within the joint. Explicit representations of fasteners may be incorporated wherever an accurate depiction of the References [1] D.A. Cope and T. E. Lacy, Stress intensity determination in lap joints with mechanical fasteners, 41st Structures, Structural Dynamics, and Materials Conference on Atlanta. 1, 3-6 (2000). [2] N.C. Bellinger and J.P. Komoroski, Corrosion pillowing stresses in fuselage lap joints. AIAA Journal. 35, 317-320(1997). [3] D.Z. Yu, Study on mechanical properties of bolted joints based on experiment and three-dimensional finite element, Journal of Mechanical Strength.,33, 855-861(2011). [4] Y. Shang, Optimal attack strategies in a dynamic botnet defense model, Applied Mathematics and Computation.6, 29-33(2012). [5] M.T. Alquran, Solitons and periodic solutions to nonlinear partial differential equations by the sine-cosine method, Applied Mathematics and Computation.6, 85-88(2012). [6] S. Furuichi and H. Hino, Mathimatical analyses of 2010 FIFA word cup, Applied Mathematics and Computation.5, 205-219(2011).

704 Dazhao YU: A Numerical Analysis of Riveted Lap Joint Containing Multiple-site Damage [7] N. Chen, Numerical simulation and experimental researches on the LED reliability under temperature loading, Applied Mathematics and Computation.6, 775-779(2012). [8] D.S. Dawicke, E.P. Phillips and D.V. Swensen, Crack Growth from Loaded Rivet Holes, International Symposium on the Structural Integrity of Aging Airplanes. 44 51 (1992). [9] M. James and D. Swenson, FRANC2D/L: A Crack Propagation Simulator for Plane Layered Structures. Version User s Guide, Cornell University. 45-60(2002) [10] J.C. Newman, C.E, Jharris and M.A. James, Fatigue- Life prediction of riveted lap splice joints using smallcrack theory. 19th International Committee on Aeronautical Fatigue Symposium on Edinburgh, Scotland, 18-20(1997) DazhaoYu received the MS degree in Mechanics science from Naval Aeronautical and Astronautical University in 2003, and the PhD degree in Mechanics science from Naval Aeronautical and Astronautical University in 2008. He is currently a professor in Naval Aeronautical and Astronautical University. His research interests are in the areas of the effects of MSD on fatigue life of aircraft.