THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES

Similar documents
Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

One Dimensional Axial Deformations

Application to Plane (rigid) frame structure

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS

829. An adaptive method for inertia force identification in cantilever under moving mass

DUE: WEDS FEB 21ST 2018

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

New Method for Solving Poisson Equation. on Irregular Domains

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

CHAPTER 4. Vector Spaces

Indeterminate pin-jointed frames (trusses)

In this section is given an overview of the common elasticity models.

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

CENTROID (AĞIRLIK MERKEZİ )

Module 3: Element Properties Lecture 1: Natural Coordinates

PART 8. Partial Differential Equations PDEs

Lecture Note 3. Eshelby s Inclusion II

CENTROID (AĞIRLIK MERKEZİ )

Lecture 12: Discrete Laplacian

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Finite Element Modelling of truss/cable structures

Grid Generation around a Cylinder by Complex Potential Functions

Chapter Newton s Method

Difference Equations

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

INDETERMINATE STRUCTURES METHOD OF CONSISTENT DEFORMATIONS (FORCE METHOD)

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

ACTM State Calculus Competition Saturday April 30, 2011

CHAPTER 4 MAX-MIN AVERAGE COMPOSITION METHOD FOR DECISION MAKING USING INTUITIONISTIC FUZZY SETS

CHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

6.3.4 Modified Euler s method of integration

Second Order Analysis

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

Consistency & Convergence

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

Effects of Boundary Conditions on Cross-Ply Laminated Composite Beams

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

The Study of Teaching-learning-based Optimization Algorithm

Plate Theories for Classical and Laminated plates Weak Formulation and Element Calculations

Numerical analysis of buried pipes subjected to lateral soil movements

A new Approach for Solving Linear Ordinary Differential Equations

SMOOTHED PARTICLE HYDRODYNAMICS METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Numerical Heat and Mass Transfer

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

11. Dynamics in Rotating Frames of Reference

Buckling analysis of single-layered FG nanoplates on elastic substrate with uneven porosities and various boundary conditions

Important Instructions to the Examiners:

A new method of boundary treatment in Heat conduction problems with finite element method

Preliminary Design of Moment-Resisting Frames

Ballot Paths Avoiding Depth Zero Patterns

Note: Please use the actual date you accessed this material in your citation.

8.592J: Solutions for Assignment 7 Spring 2005

Numerical Solution of Ordinary Differential Equations

A boundary element method with analytical integration for deformation of inhomogeneous elastic materials

Resource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS

Inductance Calculation for Conductors of Arbitrary Shape

CHAPTER-5 INFORMATION MEASURE OF FUZZY MATRIX AND FUZZY BINARY RELATION

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava

ORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp.

Numerical Simulation of Wave Propagation Using the Shallow Water Equations

PHYS 1443 Section 003 Lecture #17

1 Matrix representations of canonical matrices

x yi In chapter 14, we want to perform inference (i.e. calculate confidence intervals and perform tests of significance) in this setting.

FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

Lecture 8 Modal Analysis

Computers and Mathematics with Applications

Dynamic Programming. Preview. Dynamic Programming. Dynamic Programming. Dynamic Programming (Example: Fibonacci Sequence)

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

FUZZY FINITE ELEMENT METHOD

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

AGC Introduction

DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Solution 1 for USTC class Physics of Quantum Information

DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION

σ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review

THE EFFECT OF BEAM TO COLUMN CONNECTION IN ARC PORTAL FRAME

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

Richard Socher, Henning Peters Elements of Statistical Learning I E[X] = arg min. E[(X b) 2 ]

Δ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions

Flow Induced Vibration

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Global Sensitivity. Tuesday 20 th February, 2018

A Hybrid Variational Iteration Method for Blasius Equation

Transcription:

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES W. C. Lao Department of Cvl Engneerng, Feng Cha Unverst 00 Wen Hwa Rd, Tachung, Tawan SUMMARY: The ndentaton etween clndrcal ndentor and angle-pl compostes s nvestgated. An exact Green s functon for the surface dsplacement of an angle-pl eam s derved usng Jones s dsplacement formulaton. B matchng the dsplacements of the top surface of the angle-pl lamnate and those of the ndentor wthn the contact regon, we can fnd the contact stress dstruton, the magntude of contact force and the ndentaton. In order to speed up the formulaton of fundamental deflecton mode shape of Green s functon, method of ntal functon (MIF) s adopted to fnd the dsplacements n the top and ottom surface of an angle-pl lamnate. In the MIF, the order of sstem equaton s onl 3 3, whch wll reduce the computaton tme. Numercal results for angle-pl lamnates sujected to clndrcal ndentaton are evaluated to verf the applcalt of the proposed method. KEYWORDS: angle-pl lamnates, ndentaton, pont matchng method, method of ntal functons INTRODUCTION The contact prolems of composte lamnate are mportant n the analss of foregn oject mpact and flexural measurement of composte materals. Several authors have studed the contact ehavors etween a clndrcal ndentor and orthotropc eams[,,3]. Sanar solved the

contact prolem of orthotropc eam the superposton of half space soluton and eam soluton[]. B applng the exact dsplacement soluton for clndrcal endng of orthotropc eams, the exact Green s functon s formulated n solvng the contact prolems[3]. However the ndentaton etween clndrcal ndentor and angle-pl compostes s seldom nvestgated. In ths stud, an exact Green s functon for the surface dsplacement of an angle-pl eam s derved usng Jones s[4] dsplacement formulaton. Ths Green s functon satsfes the dsplacement and stress contnut condton at the nterface, and the prescred stress condton at the top and ottom surfaces. B matchng the dsplacements of the top surface of the angle-pl lamnate and those of the ndentor wthn the contact regon, we can fnd the contact stress dstruton, the magntude of contact force and the ndentaton. In constructng the exact Green s functon, we have to solve a matrx of order 6N 6N, wth M tmes, where N s the pl numer and M s the numer of asc flexural modes. If the pl numer ncreases, the computaton tme wll ncrease drastcall. In order to speed up the calculaton, method of ntal functon (MIF)[5] s adopted to fnd the dsplacements n the top and ottom surface of an angle-pl lamnate. In the MIF, the order of sstem equaton s onl 3 3, whch wll reduce the computaton tme. Numercal examples are llustrated to demonstrate the contact ehavors of angle-pl lamnates. CYLINDRICAL BENDING OF ANGLE-PLY LAMINATES In order to stud the contact ehavor of angle-pl lamnates, the clndrcal endng prolems of an ansotropc lamnate s nvestgated frst whch wll provde the Green s functon for contact prolem. Assumng plane stan condton n the x-drecton, an ansotropc eam has wdth and thcness h n the -drecton s shown n Fgure. The stran feld s defned as u v w ε x 0, ε, ε x v w u u γ +, γ x, γ x () The consttutve relaton for a generall ansotropc materal s xx c c c c 0 3 0 0 6 c c c c ε 3 0 0 6 c c c c ε 3 3 33 0 0 36 () 0 0 0 c c γ 44 45 0 x 0 0 0 c c 0 γ 45 55 x x c c c 0 0 c γ 6 6 36 66 x

In whch c j s the elastc constant for each laer. The eam s smpl supported along 0 and, wth the followng oundar condtons w(0,)w(,)0 (3a) ( 0, ) (, ) 0 (3) x ( 0, ) (, ) 0 (3c) x (, 0) (, 0) (, 0) 0 (3d) x ( h, ) q ( ), ( h, ) ( h, ) 0 (3e) x If the nerta force s neglected, the equlrum equaton can e wrtten as + 0 (4a) x, x, + 0 (4),, + 0 (4c),, The Dsplacement felds for ths prolem can e expressed as nπ u (, ) Un ( ) cos nπ v(, ) Vn ( ) cos nπ w (, ) Wn ( ) sn The aove canddates satsf Eqns. (3a), (3) and (3c). When susttute Eqn. (5) nto Eqn.(4) and Eqn.(), the followng equatons can e derved. '' '' pc66 U( ) + c55 U ( ) pc6 V( ) + c45v ( ) + pc ( 36 + c45) W' 0 (6) '' '' pc6 U( ) + c45 U ( ) pc V( ) + c44 V ( ) + pc ( 3 + c44 ) W' 0 (7) '' pc ( 45 + c36) U'( ) pc ( 44 + c3) V'( ) + c33 W( ) pc44 W ( ) 0 (8) In whch pnπ /, the suscrpt n of U, V and W are omtted for smplct. In order to solve the U (), V () and W (), we let [4] U ( ) U V ( ) V 0 0 W ( ) W e 0 e λ λ e λ a c55, a c44, a3 c33 p c66, p c, 3 p c44 (0) c c45, d pc6, e pc ( 36 + c45), e pc ( 3 + c44 ) B susttutng Eqns. (9) and (0) nto Eqns. (6)-(8), the nontrval solutons exst f the followng relaton holds λ a + λ c d λe det λ c d λ a + λe 0 () λe λe λ a + 3 3 (5) (9)

After solvng the 6th order characterstc equaton of λ, the fnal form of U n (), V n () and W n () can e elded as 6 * λn n U ne 6 * λn n ( ) V ne () 6 * λn n ( ) W ne U ( ) V W * * * In whch U, V and W are the correspondng egenvectors when solvng Eqn.(). The onl unnowns are n. For a fxed p, there are sx n s n each laer. These constants can e solved from the nterface contnut and the tracton condtons at the top and ottom surfaces. If the lamnate s made of N laers, the contnut of dsplacements and stresses at the nterface requres that + u (, h ) u (, 0) + v ( h, ) v (, 0) + w (, h ) w (, 0) + ( h, ) (, 0 ) (3) + ( h, ) (, 0) + x ( h, ) x (, 0),,,..., N In whch h s the thcness of the -th laer. The tracton condtons on the top and ottom surface mpose that N N (, 0) (, 0) x (, 0) ( h, N) x( h, N ) 0 N ( h, N ) q( ) (4) If the normal stress on the top surface s q( ) o sn p (5) then the coeffcents n can e solved through a set of matrx of order 6N 6N. The vertcal dsplacements w(,h) and w(,0) are of great mportance n the contact prolem. CYLINDRICAL INDENTATION OF SIMPLY SUPPORTED ANGLE-PLY LAMINATES Consder a smpl supported angle-pl lamnate sujected to the ndentaton of a rgd clnder as shown n Fg.. The radus of the ndentor s R. The frcton etween the ndentor and lamnate s neglected. The ndentaton s assumed to e smmetrc aout the mdspan of the eam. The contact length s c. In the contact prolems the dstruton of contact pressure has to e determned. The pont matchng method s adopted to fnd the contact stress dstrutons.

From prevous secton, the dsplacement at the top surface of the eam due to a normal stress o sn( nπ/ ) appled at h s shown to e nπ wh (, ) Wn ( h) sn (6) For an artrar normal stress q() at the top surface, the dsplacement of the eam wll e nπ wh (, ) qw n n( h) sn (7) In whch nπ qn q d ( )sn 0 (8) + c nπ q d ( )sn c From Eqns. (7) and (8), we have nπξ nπ w(, h) Wn ( h) q( ξ )sn dξsn Ω Ω G( ξ) q( ξ ) dξ Where Ω s contact regon, and the Green s functon s defned as[3] nπ nπξ G( ξ ) Wn ( h) sn sn (0) The relatve vertcal dsplacements of an pont w(,h) nsde the contact regon and the center pont w(0.5,h) can e expressed as[] 0.5 ( 0.5) w(, h) w(0.5, h) R R R () [ G( ξ) G(0.5 ξ )] q( ξ ) dξ Ω The contact regon s dscreted nto several segments, wthn each segment the contact stress s assumed to e constant. B the least square method and pont matchng method[,3], we can fnd the contact stress dstrutons. (9) METHOD OF INITIAL FUNCTIONS In constructng the Green s functon for the contact prolems, the vertcal deflecton at top surface W n (h) due to a snusodal stress s needed. If the laer numer ncreases, t wll e tme consumng to fnd the asc deflecton mode. The method of ntal functon[5] s appled to smplf the calculaton of (h). From the stran defntons, the stress-stran relaton and W n equlrum equatons, the followng equaton can e derved[5]:

U V Z 0 A X B 0 Y W In whch c44 / d c45 / d 0 A c / / 45 d c55 d 0 0 U V Z X Y W () (3a) * c66 * B c6 c36 c33 c c * 6 * c c 3 33 c36 c 33 c3 c 33 / c33 (3) d c44c55 c45,,, c c * 3 j3 cj cj,, j,, 3, 6 (4) c33 X, Y, Z x The materal constants for the -th laer should e susttuted n Equatons () and (3). The soluton of Eqn. () s shown to e R( ) e Where 0 B A 0 R(0) R( ) T { U( ), V( ), Z( ), X( ), Y( ), W( )} (6) The state vector R() n the -th laer at artrar heght can e expressed n terms of the state vector at the ntal reference heght of that laer. B a successve multplcaton of the state vector at each laer, W n (h) can e found solvng a set of 6 6 smultaneous equaton. (5)

NUMERICAL EXAMPLES A [± 45] S composte lamnate made of Hexcel graphte/epox s analed as an llustratve example. The eam has a wdth of 0 cm and thcness of cm. Each lamna assumes the followng propertes, E 05.7GPa, E E3 8.57GPa, υ υ3 0.37, υ3 0.306 G G 4.39GPa, G 3.05GPa 3 3 The radus R of the rgd ndentor s 7.6cm. The through-the-thcness stress dstruton s calculated for contact length cmm. Fgure 3 shows the thcness-wse dstruton of the normal stress. The results at three dfferent postons (the contact center, the contact edge and end of eam) are compared. It seems that the normal stress wll have maxmum value at the top surface of the contact center. The nterlamnar stress dstrutons of and x are demonstrated n Fgs. 4 and 5. The solutons satsf the tracton free condtons at oth top and ottom surfaces. Orthotropc solutons are also compared, n whch the effectve elastc constants of ths angle-pl lamnate s calculated usng a 3-D effectve modul theor[6]. The orthotropc results are comparale to the laer-wse solutons. Fgure 6 shows the n-plane stress dstruton of. The relaton of ndentaton and contact force s evaluated n Fg. 7. The ndentaton s xx defned as the dfference of vertcal dsplacement of w(0.5, h) and w(0.5, 0). Fgure 8 s the contact pressure dstruton for dfferent contact length. In order to verf the applcalt of method of ntal functon, results of W n (h) for dfferent n are revealed n Fg. 9. It seems that oth methods agree well for ths case. In some other case, f n s large, the successve multplcaton of state vectors wll cause some numercal nstalt. If ths nstalt can e overcome, the method of ntal functon wll e an effcent tool for contact analss. CONCLUSIONS A method for the analss of ndentaton of angle-pl lamnates sujected to rgd clndrcal ndentor s proposed. Numercal results show that ths method satsfes the gven oundar condtons and tracton condtons. The contact pressure dstruton, magntude of contact force and ndentaton can e evaluated. Ths analtcal results wll compare to the expermental results n the near future.

REFERENCES. Sanar, B. V., Smooth Indentaton of Orthotropc Beams, Compostes Scence and Technolog, 34, (989) 95-.. Sun, C. T. and Sanar, B. V., Smooth Indentaton of an ntall stressed orthotropc eam, Int. J. Solds and Structures,, (985), 6-76. 3. Wu, E., Chao, J. C. and Yen, C. S., Smooth Contact of Orthotropc Lamnates Rgd Clnders, AIAA Journal, 3, 0, (993), 96-9. 4. Jones, A. T., Exact Natural Frequences and Modal Functons for a Thc Off-Axs Lamna, J. of Composte Materals, 5, Oct. (97), 504-50. 5. Sundara Raja Iengar, K. T. and Panda, S. K., Analss of Orthotropc Rectangular Thc Plates, Fre Scence and Technolog, 8 (983), 9-36. 6. Sun, C. T. and Lee, S., Three-dmensonal Effectve Elastc Constants for Thc Lamnates, J. of Composte Materals, 988,, 69-636.