Fastener Modeling for Joining Parts Modeled by Shell and Solid Elements

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

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

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

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

IDEA Connections. User guide

Modeling Multi-Bolted Systems

TUTORIAL 4: Combined Axial and Bending Problem Sketch Path Sweep Initial Project Space Setup Static Structural ANSYS

Stress Analysis of Flanged Joint Using Finite Element Method

Prying of a Large Span Base Plate Undergoing a Moment Load Applied by a Round Pier

Development of a Numerical Technique for the Static Analysis of Bolted Joints by the FEM

IDEA Connection 8. User guide. IDEA Connection user guide

Corporate Subscription. NAFEMS reference library at the click of a button

Abaqus CAE (ver. 6.9) Contact Tutorial

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

Korean standards of visual grading and establishing allowable properties of softwood structural lumber

ANALYSIS OF LATERAL STIFFNESS FOR INFILLED FRAME WITH OPENING

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

Instruction Manual for installing

Instruction Manual for installing

Finite Element Study of Using Concrete Tie Beams to Reduce Differential Settlement Between Footings

Creasing and Folding are critical steps in the box forming process. Dipl. Ing. Lukas Pescoller

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

AN INNOVATIVE FEA METHODOLOGY FOR MODELING FASTENERS

Monopile as Part of Aeroelastic Wind Turbine Simulation Code

IJSRD - International Journal for Scientific Research & Development Vol. 4, Issue 05, 2016 ISSN (online):

Finite Element Analysis per ASME B31.3

INFLUENCE OF PILES ON LOAD- SETTLEMENT BEHAVIOUR OF RAFT FOUNDATION

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

Status of Coil Structural Design and Magnetic-Structural Analysis

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

Abaqus Beam Tutorial (ver. 6.12)

AMTS STANDARD WORKSHOP PRACTICE. Bond Design

Isotropic. Anisotropic. Orthogonal. Isotropic orthogonal

Quasi-static Contact Mechanics Problem

Monitoring The Machine Elements In Lathe Using Vibration Signals

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

GENERAL GUIDELINES FOR APPLICATION OF THE EXTENDED SUBTRACTION METHOD IN SASSI SOIL-STRUCTURE INTERACTION ANALYSIS

EFFECTS OF GEOMETRY ON MECHANICAL BEHAVIOR OF DOVETAIL CONNECTION

Note: Conditions where bending loads are imposed on the bolt e.g. non-parallel bolting surfaces, should be avoided.

Abaqus/CAE (ver. 6.14*) Plate/Shell Tutorial

IDEA Connection 8. User guide. IDEA Connection user guide

The predictive model for strength of inclined screws as shear connection in timber-concrete composite floor

REVIEW OF THREADED FASTENERS LOOSENING AND ITS EFFECTS

Structural and Thermal Analysis of Bolted joint of Coiler Drum in Steckel Mill using Finite Element Method

Glued laminated timber beams repair.

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

UC San Diego UC San Diego Electronic Theses and Dissertations

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

Modal vibration control of submarine hulls

An Alternative Formulation for Determining Stiffness of Members with Bolted Connections

Effect of Yarn Twist on Young s Modulus of Fully-green Composites Reinforced with Ramie Woven Fabrics ABSTRACT

NON-LINEAR CONNECTION MODELS IN TIMBER ENGINEERING

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

Module 10 : Improvement of rock mass responses. Content

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

EFFECT OF HOLE CLEARANCE ON BOLT LOADS IN PULTRUDED GRP TENSION JOINTS

Mechanical joints. Major diameter Mean diameter Minor diameter Pitch p chamfer. Root Crest. Thread angle 2a. Dr. Salah Gasim Ahmed YIC 1

Thread Definitions -The terminology of screw threads

Advance Steel. Tutorial

Bolt Material Types and Grades 1- Bolts made of carbon steel and alloy steel: 4.6, 4.8, 5.6, 5.8, 6.8, 8.8, 10.9 Nuts made of carbon steel and alloy

Effect of Pile Bending Stiffness on Static Lateral Behavior of a Short Monopile in Dry Sand

CH # 8. Two rectangular metal pieces, the aim is to join them

RlGIDITY AND STRENGTH OF WALL FRAMES BRACED WlTH METAL STRAPPING

Should We Upgrade to Phase 2 6.0?

Optimizing the Natural Frequencies of Beams via Notch Stamping

AutoCAD Inventor - Solid Modeling, Stress and Dynamic Analysis

Dowel connections in laminated strand lumber

FAST VISUALISATION OF SAFETY MARGINS OF THE W7-X PLASMA VESSEL

Connection and performance of two-way CLT plates

Wear Analysis of Multi Point Milling Cutter using FEA

MOULDABILITY OF ANGLE INTERLOCK FABRICS

Vertex Detector Mechanics

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

Design of Machine Elements I Prof. G. Chakraborty Department of Mechanical Engineering Indian Institute of Technology Kharagpur

IDEA CSS 7 General cross-section

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

CHAPTER 2 ELECTROMAGNETIC FORCE AND DEFORMATION

Eurocode EN Eurocode 3: 3 Design of steel structures. Part 1-1: General rules and rules for buildings

NUMERICAL ANALYSIS OF SCREW ANCHOR FOR CONCRETE

CATIA V5 Workbook Release V5-6R2013

Acoustic Emission Signals versus Propagation Direction for Hybrid Composite Layup with Large Stiffness Differences versus Direction

Studies on free vibration of FRP aircraft Instruments panel boards

Stress Analysis of T-Flange Bolted Joint with a Simplified Spring and Beam Model

STRUCTURAL TIMBER DESIGN

Bhagwan mahavir college of Engineering & Technology, Surat.

Comparative Evaluation of Resistance Made Simple Shear Connection with Bolts and With Welding

Beam & Header Technical Guide. LP SolidStart LVL. 2900F b -2.0E. U.S. Technical Guide U.S. TECHNICAL GUIDE

APA Performance Rated Rim Boards

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

Combined stress analysis of mitered spline furniture joints under diagonal loading

RECENTLY DESIGNED BOW-STRING RAILWAY BRIDGES IN SLOVAKIA

Fatigue crack propagation in uniaxial loading and bending fatigue in 20 khz testing. Mohamed Sadek PhD Student Karlstad university

Modeling Basic Mechanical Components #1 Tie-Wrap Clip

Load-carrying capacity of timber frame diaphragms with unidirectional support

A study of Vibration Analysis for Gearbox Casing Using Finite Element Analysis

STUDY AND ANALYSIS OF ANGULAR TORQUING OF ENGINE CYLINDER-HEAD BOLTS USING TORQUE-TO-YIELD BOLTS: A CASE STUDY

Introduction to Engineering Design

ABSTRACT II. OBJECTIVE I. INTRODUCTION III. METHODOLOGY

Fatigue Analysis of VMC 450 Spindle

MODELLING OF CONCRETE PAVEMENT DOWEL-SLAB INTERACTION

Transcription:

2007-08 Fastener Modeling for Joining Parts Modeled by Shell and Solid Elements Aleander Rutman, Chris Boshers Spirit AeroSystems Larry Pearce, John Parady MSC.Software Corporation 2007 Americas Virtual Product Development Conference October 11-12, 2007 Detroit, MI Abstract The continued trend in finite element modeling of structures is toward the increased use of solid elements. For FE analysis of components containing two or more mechanically attached parts, it is often necessary to model each fastener to adequately represent the load transfer in the joint. It is desirable to be able to rapidly prepare model of each fastener to represent the actual joint behavior in a FEM containing solid elements, shell elements, or both. The FEM representation of the fastener should capture the shear and bending stiffness of the fastener and the bearing stiffness at the fastener joined part interaction. The stiffness of a fastener joint is a function not only of the fastener geometry and material, but also the properties of the joined structural members. A set of rigid elements should be used to ensure proper interaction of joined parts. The FEM fastener representation previously proposed meets this requirement quite well, and can be used to quickly prepare fastener representations when the structural members are modeled using shell elements. The work presented in this paper etends the fastener FEM formulation previously developed to enable its use with solid elements. This paper describes a procedure to automate the creation of nodes, elements, MPCs, and constraints using MSC.Nastran to represent individual fasteners to develop mechanical connections between structural members modeled by shell elements, solid elements, or both. The procedures described in this work can be used with any number of joint members, and with a variety of element types. 1

1. Introduction. With an increase of computer speed along with the amount of available memory, the trend for creation of more detailed models has arisen. It includes not only fine mesh modeling using shell elements but also increased use of solid elements due to etensive application of 3D CAD software packages and the desire of analysts to more realistically represent structural parts and their interaction, including fastener joints. The distribution of loads between the components of a structural assembly depends not only on their dimensions and material properties but also on the stiffness of the fastener joints connecting the components. The accuracy of the finite element analysis is therefore significantly influenced by the representation of fastener joints in the model. A procedure for modeling of fastener joints when the connected parts are represented by shell elements located at the plate mid-planes [1] was developed earlier. However, the etensive use of this procedure during these years by a number of companies has shown the necessity for additional capabilities, including: Modeling of fastener joints when all or some of the connected structural components are represented by solid elements; Ability to obtain tensile loads in fasteners directly from MSC.Nastran output. These requirements are addressed in the presented paper. Other potential capabilities, i.e., the non-linear analysis of models with fastener joints, or the analysis of composite structures, will be addressed later. The work presented in this paper etends the fastener FEM formulation to enable its use with solid elements. The paper describes a procedure for creation of a set of nodes, elements, MPCs, and constraints using MSC.Nastran to represent individual fastener joints. This fastener representation can then be used to model mechanical connections between structural members that are represented by shell elements, solid elements, or any combination of the two. The procedures described in this work can be used with any number of joint members, and can be used with a variety of element types, including quad, tri, tet, and he elements. Higher order elements are also supported. To enable the reading of tensile loads in fasteners directly from MSC.Nastran output, an alternative modeling procedure based on the linear gap techniques [2] was developed. However, this procedure does not consider the effect of the fastener pre-tension and fit. Also, the friction phenomenon in the joint is not considered here. These effects require non-linear analysis and therefore fall beyond the scope of this paper. 2

2. Fastener Modeling Approach. The presented approach to modeling of a fastener joint is based on the definition of each deformation component contributing to joint fleibility, then modeling them using a corresponding finite element. The combination of these elements represents the complete behavior of a fastener joint. Some relative displacements in the fastener joint model are limited to ensure the compatibility of deformations. 2.1. Stiffness of Fastener Joint In a fastener joint (Figure 1) the following stiffness components are considered [1]: Combined plate and fastener translational bearing stiffness; Combined plate and fastener rotational bearing stiffness; Fastener shank shear and bending stiffness. Fastener t p 1 t p 2 t p 3 t p 4 d f 2 nd plate 3 rd plate 4 th plate Figure 1. Fastener Joint. The combined plate i and fastener translational bearing stiffness is S pi bt = i 1 1 + E cpi t E cf where E cp i - compression modulus of plate i material; E cf - compression modulus of fastener material; t - thickness of plate i. p i The combined plate i and fastener rotational bearing stiffness is 3

S bt 1 = 12 1 E cp i t 3 p i i 1 + E cf The bearing stiffness is modeled by elastic elements. A more detailed definition of bearing stiffness is given in [1]. The shear and bending stiffness of a fastener are represented by a beam element. If a plate is composed of a laminated composite material, a reasonable estimate of the bearing translational and rotational stiffness can be obtained by using an average of the apparent in-plane stiffnesses for E cp i. The apparent in-plane stiffnesses for a given composite laminate can be determined using classical lamination theory [3,4]. 2.2. Representation of a Fastener Joint Modeling of a fastener joint is described here in MSC.Nastran terms. Idealization of a plate-fastener system includes the following: Elastic bearing stiffness of a plate and fastener at the contact surface; Bending and shear stiffness of a fastener shank; Compatibility of displacements of a fastener and connected plates in the joint. The presented method creates a plate-fastener system that is illustrated in eamples shown in Figures 2 7: a. Model of a fastener joint connecting three plates modeled by shell elements (Figure 2). b. Model of a fastener joint connecting three plates modeled by single ply of solid elements (Figure 3). c. Model of a fastener joint connecting three plates modeled by two plies of solid elements (Figure 4). Modeling of fastener joints connecting plates represented by high order solid elements is similar to modeling shown in Figure 4. d. Model of a fastener joint connecting three plates modeled as follows (Figure 5): o Plate 1 solid elements; o Plate 2 shell elements; o Plate 3 solid elements. e. Model of a fastener joint connecting three plates modeled as follows (Figure 6): o Plate 1 shell elements; o Plate 2 solid elements; o Plate 3 solid elements. f. Model of a fastener joint connecting three plates modeled as follows (Figure 7): o Plate 1 solid elements; o Plate 2 solid elements; o Plate 3 shell elements. 4

a) b) 2 nd plate 2 nd plate 3 rd plate 3 rd plate Figure 2. Model of a Fastener Joint Connecting Three Plates Modeled by s: a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. A description of each of the components of the fastener joint models shown in Figures 2 7 is given in Table 1. The method requires the creation of two additional sets of nodes: Nodes at the plate mid-planes for plates represented by solid elements (gray nodes in Figures 3 7). For plates represented by shell elements, additional nodes are not required because the elements are already located at the plate mid-plane. Nodes for the fastener shank modeling (black nodes in Figures 2 7). This set includes nodes coincident with corresponding nodes at plate mid-planes and nodes located on the intersection of the fastener ais and the outer surfaces of the first and last connected plates. For correct definition of a fastener shear plane and its aial direction, a coordinate system with one of its aes parallel to the fastener ais must be defined in the bulk data. This coordinate system must be used as the analysis coordinate system for both the eisting nodes (white solid element nodes and gray shell element nodes in Figures 2 7) and newly created nodes (gray mid-plane nodes of plates modeled by solid elements and black nodes of a fastener shank in Figures 2-7). 5

a) b) 2 nd plate 2 nd plate 3 rd plate 3 rd plate Figure 3. Model of a Fastener Joint Connecting Three Plates Modeled by s: a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. 2.3. Fastener Shank Modeling The fastener shank is modeled by CBAR or CBEAM elements [5] with corresponding PBAR or PBEAM card for properties definition. For the CBAR or CBEAM elements connectivity, a separate set of nodes (black nodes in Figures 2 7) is created. This set also includes nodes located on the intersection of the fastener ais and the outer surfaces of the first and last connected plates. All CBAR or CBEAM elements representing the same fastener shank reference the same PBAR or PBEAM card [5] with the same cross-sectional area, moments of inertia and torsional constant: 2 A π d 4 4 f π d f π d f = I1 = I 2 = J = 4 64 32 where d f - fastener shank diameter. 6

a) b) Solid Elements Solid Elements Solid Elements 2 nd plate Solid Elements 2 nd plate 3 rd plate 3 rd plate Solid Elements Solid Elements Figure 4. Model of a Fastener Joint Connecting Three Plates Modeled by Two Plies of s: a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. 2.4. Modeling of Interaction Between Fastener Shank and Joined Plates The interaction between a fastener and plates results in bearing deformation of all components of the joint on their surfaces of contact. The combined bearing stiffness of a fastener and connected plates is defined in Section 2.1. The bearing stiffness is presented as a translational stiffness in the direction of aes normal to the fastener ais which define the fastener shear plane, and rotational stiffness about the same aes. For the modeling of the bearing stiffness, both additional sets of grid points mentioned above are used. Each pair of coincident grid points created by a plate mid-plane node for a solid element or shell element node (gray node in Figures 2 7) and the corresponding fastener shank node (black node in Figures 2 7) is connected by a CBUSH element [5] or combination of CELAS2 elements with equal translational stiffness along the aes normal to the fastener ais and equal rotational stiffness about the same aes. The connectivity card CBUSH must be accompanied by a PBUSH card defining the element stiffness. When the plate is modeled by shell elements, the load from the CBUSH element is transferred directly to the shell elements. In the case of the plate modeled by solid elements, load from the plate mid-plane node is transferred to the solid element nodes by a rigid element RBE2. 7

a) b) 2 nd plate 2 nd plate 3 rd plate 3 rd plate Figure 5. Model of a Fastener Joint Connecting Three Plates Modeled as Follows: o Plate 1 s; o Plate 2 s; o Plate 3 s. a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. 2.5. Compatibility of Displacements in the Joint As mentioned above, two modeling procedures are considered here for representation of fastener joints: 1. Modeling using rigid elements only. This approach is recommended when tensile loads in fasteners can be neglected and shear loads on fasteners are the main goal of analysis. In this case load in the direction transverse to the joined plates is transferred through plate contact and the fastener shank is not loaded. 2. Modeling using rigid elements and linear gap techniques. This approach should be applied when the tensile loads in the fasteners are significant and are one of main goals of analysis. However the ecessive use of linear gaps can slow model convergence. The fastener joint model was designed under the following assumptions: No interference of plates under the load; The plate mid-planes stay parallel to each other under load; Planes under the fastener heads stay parallel to the plate mid-planes under load. 8

a) b) 2 nd plate 2 nd plate 3 rd plate 3 rd plate Figure 6. Model of a Fastener Joint Connecting Three Plates Modeled as Follows: o Plate 1 s; o Plate 2 s; o Plate 3 s. a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. These goals are reached by using two sets of RBAR elements and the linear gap techniques. When the first modeling approach (rigid elements only) is used the first set of RBAR elements (black and green arrows in Figures 2 7) satisfies the first assumption. These RBAR elements prevent the movement of plates relative to each other along the fastener ais. The first RBAR element in this set also prevents the fastener movement as a rigid body along its ais. In this case not only is the interference of plates prevented but a gap development in the joint is banned as well. 9

a) b) 2 nd plate 2 nd plate 3 rd plate 3 rd plate Figure 7. Model of a Fastener Joint Connecting Three Plates Modeled as Follows: o Plate 1 s; o Plate 2 s; o Plate 3 s. a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. When the alternative modeling approach (rigid elements and linear gaps) is used, the degree of freedom of RBAR elements in the fastener ais direction is replaced by linear gaps (red arrows in Figures 2 7). The eception is the RBAR element between the upper fastener head and the first plate because it prevents free body movement of the fastener shank. Under this approach, gaps between plates are allowed but interference of plates is banned. The second set of RBAR elements (blue and yellow arrows in Figures 2 7) satisfies the last two assumptions. It forces the plate mid-planes and planes under the fastener heads to remain parallel to each other under load. Two sets of RBAR elements are required when the joined plates are modeled by solid elements because solid elements do not have rotational stiffness at their nodes. 10

Table 1. Components of Fastener Joint Models Component GRID GRID GRID GRID CBAR CBUSH Component Description Fastener shank node (node of fastener shank under the upper head must be constrained against rotation about fastener ais) Node at plate mid-plane Node of solid element involved in the fastener joint modeling Node of solid element not involved in the fastener joint modeling Element for fastener shank modeling Element connecting plate mid-plane node with coincident node of the fastener shank DOF s for Rigid Elements and Linear Gaps for Fastener Ais Directions 1 2 3 2 3 5 6 1 3 4 6 1 2 4 5 RBE2 Rigid element connecting plate mid-plane node with nodes of solid element representing plate 2 3 1 3 1 2 RBAR RBAR Rigid element connecting node of the fastener shank under the upper head (independent node) with the node of shell element representing the first plate (dependent node) Rigid element connecting nodes of shell elements representing two adjacent plates 1 5 6 2 4 6 3 4 5 1 5 6 2 4 6 3 4 5 RBAR RBAR RBAR RBAR Rigid element connecting mid-plane nodes of two adjacent plates Rigid element connecting mid-plane node of the last plate (independent node) with the node of the fastener shank under the lower head (dependent node) Rigid element connecting node of the fastener shank under the upper head (independent node) with mid-plane node of the first plate represented by solid elements (dependent node) Rigid element connecting mid-plane nodes of two adjacent plates 5 6 4 6 4 5 5 6 4 6 4 5 1 4 5 6 2 4 5 6 3 4 5 6 1 4 5 6 2 4 5 6 3 4 5 6 11

Table 1 (continued). Components of Fastener Joint Models Component RBAR RBAR RBAR RBAR Linear Gap Linear Gap Linear Gap Linear Gap Linear Gap Component Description Rigid element connecting node of the fastener shank under the upper head (independent node) with the coincident node of solid element representing the first plate (dependent node) RBAR element connecting two coincident nodes of solid elements representing two adjacent plates RBAR element connecting the corner node of solid element (independent node) with midplane node of adjacent plate modeled by shell element (dependent node) RBAR element connecting mid-plane node of plate modeled by shell elements (independent node) with the corner node of adjacent plate modeled by solid elements (dependent node) Technique connecting mid-plane nodes of two adjacent plates modeled by shell elements Technique connecting mid-plane node of the last plate modeled by shell elements with the node of the fastener shank under the lower head Technique connecting two coincident nodes of solid elements representing two adjacent plates Technique connecting node of solid element representing the last plate with the coincident node of the fastener shank under the lower head Technique connecting the corner node of plate modeled by solid elements with mid-plane node of adjacent plate modeled by shell elements DOF s for Rigid Elements and Linear Gaps for Fastener Ais Directions 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 If the degrees of freedom of the second set of RBAR elements (blue and yellow arrows) were assigned to the first set of RBAR s (black and green arrows) the last two conditions would not be satisfied because of the absence of rotational stiffness at the solid element nodes. If the degrees of freedom of the first set of RBAR elements (black and green arrows) were assigned to the second set of RBAR elements (blue and yellow arrows) gap or interference development between plate surfaces cannot be prevented. So, when plates are modeled by solid elements the two sets of RBAR elements described above is an absolute necessity for correct formulation of the fastener joint behavior under load. However, when plates are modeled by shell elements, all required degrees of freedom can be assigned to one set of rigid elements (Figure 2). 12

3. Modeling System. A block approach to fastener joint modeling is presented in this section. A model of a fastener joint can be divided into three main parts (Figures 2 7): Interaction of the upper fastener head and fastener shank with the first plate; Interaction of two adjacent plates with each other and with the fastener shank; Interaction of the lower fastener head and fastener shank with the last plate. Each of these model parts can be represented by a few related blocks. Using these blocks a model of a fastener joint connecting any number of joined plates represented by shell or solid elements can be created. 3.1. Interaction of Fastener with the First Plate Interaction of a fastener and the first plate is presented by two blocks (Figure 8) for two possible modeling conditions: The first plate is modeled by shell elements; The first plate is modeled by solid elements. These blocks are the same for both modeling approaches: modeling using rigid elements only and using rigid elements and linear gaps. Rigid elements in these blocks prevent relative rotation of the plate and the fastener head about the aes in the fastener shear plane and relative displacement along the fastener shank. The fastener shank node under the upper head is constrained against rotation about the fastener ais to prevent shank free body rotation. 3.2. Interaction of Two Adjacent Plates and Fastener Shank Interaction of a fastener shank with two adjacent plates is presented by eight blocks (Figures 9 12) for the following modeling conditions: Both adjacent plates are modeled by shell elements (Figure 9); Both adjacent plates are modeled by solid elements (Figure 10); One plate is modeled by shell elements and the adjacent plate is modeled by solid elements (Figure 11); One plate is modeled by solid elements and the adjacent plate is modeled by shell elements (Figure 12). For each of these modeling conditions two separate blocks were developed based on two modeling approaches: modeling using rigid elements only and using rigid elements and linear gaps. Rigid elements keep plate mid-planes parallel to each other, i.e. they prevent relative rotation of the plates about the aes in the fastener shear plane. Also, rigid elements and linear gaps prevent the interference between plates. 13

3.3. Interaction of Fastener with the Last Plate Interaction of a fastener and the last plate is presented by four blocks (Figures 13 and 14) for two possible modeling conditions: The last plate is modeled by shell elements; The last plate is modeled by solid elements. Here also for each of these modeling conditions two separate blocks were developed based on the two modeling approaches described above. Rigid elements in these blocks prevent relative rotation of the plate and the fastener head about aes in the fastener shear plane. The linear gap techniques transfer load from the plate to the fastener head. a) b) Figure 8. Interaction of Fastener and the First Plate: a. The First Plate is Modeled by s; b. The First Plate is Modeled by s. a) b) Plate n Plate n Plate n+1 Plate n+1 Figure 9. Interaction of a Fastener Shank with Two Adjacent Plates Modeled by s: a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. 14

a) b) Plate n Plate n Plate n+1 Plate n+1 Figure 10. Interaction of a Fastener Shank with Two Adjacent Plates Modeled by s: a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. a) b) Plate n Plate n Plate n+1 Plate n+1 Figure 11. Interaction of a Fastener Shank with a Plate Modeled by Shell Elements and Adjacent Plate Modeled by s: a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. 15

a) b) Plate n Plate n Plate n+1 Plate n+1 Figure 12. Interaction of a Fastener Shank with a Plate Modeled by Solid Elements and Adjacent Plate Modeled by s: a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. a) b) Last Plate Last Plate Figure 13. Interaction of Fastener and the Last Plate Modeled by s: c. Model with Rigid Elements Only. d. Model with Rigid Elements and Linear Gaps. 16

a) b) Last Plate Last Plate Figure 14. Interaction of Fastener and the Last Plate Modeled by s: a. Model with Rigid Elements Only. b. Model with Rigid Elements and Linear Gaps. 4. Modeling Eamples. The following eamples of fastener joint modeling are presented in this paper: 1. Symmetric double shear joint. Plates are modeled using shell elements (Figures 15 17). For the plate interaction the Rigid Elements Only modeling approach is used. 2. Symmetric double shear joint. Plates are modeled using solid elements (Figures 18 20). The inner plate is modeled by three plies of solid element, and the outer plates are represented by a single ply of solid elements. For the plate interaction the Rigid Elements Only modeling approach is used. 3. Symmetric double shear joint. The inner plate is modeled by shell elements (Figures 21 23), and the outer plates are represented by a single ply of solid elements. For the plate interaction the Rigid Elements Only modeling approach is used. 4. Attachment of a tension clip to a plate. The tension clip is modeled by solid elements and the plate is represented by shell elements (Figures 24 and 25). For the model parts interaction the Rigid Elements and Linear Gaps approach is used. In the first three eamples, load is transferred in the fastener joints by shear in fasteners. Aial loads in the fasteners are negligible. Therefore for modeling of joined plate interaction the Rigid Elements Only approach was used. In the fourth eample, load is transferred by tension in fasteners. In this case interaction of the model parts is represented using the Rigid Elements and Linear Gaps approach. In the first three eamples, under the load, the joined plates slide along each other due to combined plates and fastener translational bearing deformation and the fastener bending and shear deformation. The fastener deformation causes a change of angle between the 17

fastener and the plate or, in other words, their relative rotation. This relative rotation results in non-uniform distribution of bearing stress through the plate thickness. The resultant load transferred through the contact area between the fastener and plate consists of a force in the plate mid-plane and out-of-plane moment. In the structure, the moment is reacted by loads on the plate contact surfaces and does not cause local bending of the plates, i.e. fastener does not guide the plates. The proposed modeling technique takes this phenomenon into account ensuring the plates mid-planes stay parallel to each other under load. Reactions Fasteners Inner Plate Outer Plates Figure 15. Eample 1 Symmetric Double Shear Joint with Plates Modeled by s. Applied Loads Undeformed Model Fastener Deformed Model Fastener Fastener Figure 16. Eample 1. Model Before and After Deformation. 18

Rigid Body Movement Fastener Bending and Shear Deformation Outer Plate Fastener before Deformation Fastener after Deformation Inner Plate Outer Plate Outer Plate Bearing Deformation Inner Plate Bearing Deformation Figure 17. Eample 1. Fastener Joint Deformation. Reactions Fasteners Inner Plate Outer Plates Applied Load (Pressure) Figure 18. Eample 2 Symmetric Double Shear Joint with Plates Modeled by s. 19

Undeformed Model Fastener Deformed Model Fastener Fastener Figure 19. Eample 2. Model Before and After Deformation. Rigid Body Movement Fastener Bending and Shear Deformation Outer Plate Fastener after Deformation Fastener before Deformation Inner Plate Outer Plate Outer Plate Bearing Deformation Inner Plate Bearing Deformation Figure 20. Eample 2. Fastener Joint Deformation. 20

Reactions Fasteners Inner Plate Outer Plates Applied Loads Figure 21. Eample 3 Symmetric Double Shear Joint with Outer Plates Modeled by s and Inner Plate Modeled by Shell Elements. Undeformed Model Fastener Deformed Model Fastener Fastener Figure 22. Eample 3. Model Before and After Deformation. 21

Rigid Body Movement Fastener Bending and Shear Deformation Outer Plate Fastener before Deformation Fastener after Deformation Inner Plate Outer Plate Outer Plate Bearing Deformation Inner Plate Bearing Deformation Figure 23. Eample 3. Fastener Joint Deformation. Applied Loads Fastener Fastener Plate Reactions Reactions Tension Clip Figure 24. Eample 4 Attachment of Tension Clip Modeled by Solid Elements to Plate Modeled by s. 22

Undeformed Model Deformed Model Tension Clip Tension Clip Plate Plate Fasteners Fasteners Figure 25. Eample 4. Model Before and After Deformation. 5. Patran Utility. To automate the creation of a fastener joint model, the Patran utility was developed. The input form for the utility is shown in Figure 26. This form has following section where the data necessary for modeling of a fastener joint is defined: Element and node ID control. The starting ID of the new nodes and elements can be selected, as it is usually done in majority of input forms. Definition of the group, to which all fastener joint entities will be added. Fastener parameters definition: o Fastener shank diameter. o Grip length. o Fastener material This button opens a subform for a material selection for the fastener. o Modeling of joined plates contact rigid elements or linear gaps. o Definition of symmetry condition [1]: 1.0 - for fasteners not located on the model symmetry planes, 0.5 - for fasteners belonging to one plane of symmetry, 0.25 - for fasteners located on the intersection of symmetry planes. 23

Element and Node ID Control Group to which all the Fastener Joint Entities will be Added Fastener Data: Diameter Grip Length Fastener Material Modeling of Plates Contact Symmetry Condition Fastener Ais and Alignment: Fastener Ais List Node Align Tolerance Nudge Node to Fastener Ais Check Bo Fastener Locations Defined by Point Representing Fastener (One Point per Fastener): List of Points Group: All Nodes or Points in the Group File: Point Locations Imported from a File Fastener Connections (Defines nodes for Connection): All Nodes Nodes in Current Viewport Selected Nodes Nodes of Selected Elements Nodes of Elements Associated to Selected Properties Figure 26. Patran Utility Input Form. 24

Fastener ais definition and nodes alignment: o Ais list Defines eisting coordinate system ID and coordinate ais parallel to the fastener ais, o Nodes align tolerance Nodes within this tolerance will be considered as possible part of the fastener joint, o Nudge node to fastener ais check bo If the bo is checked the selected nodes that do not locate on the fastener ais will be moved (projected) to the ais. Fastener locations Defined by point representing the fastener (one point per fastener). Points can be specified as: o Point list This includes geometric points, nodes, yz locations, etc., o Group All points or nodes in the group will be considered as fastener locations, o File - Point locations will be imported from a file. Fastener connection Defines the nodes that will be considered for connection by fastener joint. Nodes can be specified as: o All nodes - All nodes in the database will by considered, o Current viewport - All nodes in the current viewport will be considered, o Nodes - Only selected nodes will be considered, o Elements - Nodes associated to selected elements will be considered, o Properties - Nodes associated to elements associated to selected property sets. 6. Conclusion. The modeling techniques described above for connected parts represented by shell and solid elements reflects the entire fastener joint behavior including bending, shear, and bearing fleibility and compatibility of displacements in the joint. The modeling approach is presented in terms of MSC.Nastran using CBAR, CBUSH, RBAR, and RBE2 elements and the linear gap techniques; however, it can be applied to another finite element code with similar elements. Acknowledgement The authors would like to acknowledge the contributions and advice of Mr. Robert Bennett of Lockheed Martin. Mr. Bennett's assistance was very valuable in the formulation of fastener modeling techniques using linear gaps. 25

References 1. A. Rutman, A. Viisoreanu, J. A. Parady. Fasteners Modeling for MSC.Nastran Finite Element Analysis. Paper No. 2000-01-5585. Proceedings of the 2000 World Aviation Conference, San Diego, October 2000. 2. L. Proctor. Using Linear Contact in MSC.Nastran and MSC.Patran, 2007, White Paper, The MacNeal-Schwendler Corporation, Los Angeles, CA 3. Jones, R.M., Mechanics of Composite Materials, McGraw-Hill Book Company, New York, 1975 4. Daniel I.M. and Ishai O., Engineering Mechanics of Composite Materials, Oford University Press, New York, 2006 5. MSC.Nastran Quick Reference Guide, 2004, The MacNeal-Schwendler Corporation, Los Angeles, CA Contact Information Aleander Rutman, Ph.D. Larry Pearce Spirit Aerosystems MSC.Sofware (316) 523-7048 (770) 252-6560 aleander.rutman@spiritaero.com larry.pearce@mscsoftware.com Chris Boshers John Parady Spirit Aerosystems MSC.Software (316) 523-0374 (817) 481-4812 15 chris.boshers@spiritaero.com john.parady@mscsoftware.com 26