Transactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
|
|
- Piers Casey
- 5 years ago
- Views:
Transcription
1 A path-independent integral for the calculation of stress intensity factors in three-dimensional bodies C.M. Bainbridge," M.H. Aliabadi," D.P. Rooke* "Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton SO40 7AA, UK ''Structural Materials Centre, Defence Research Agency, Farnborough GUN 6TD, UK Abstract A path-independent surface integral is derived for the calculation of the stress intensity factor in three dimensions; it is similar to an existing contour integral developed for two dimensional problems. The integral is based upon a type of Green's function, which is used in conjunction with Betti's reciprocal theorem. The proposed integral is used to evaluate stress intensity factors for different benchmark problems, and the results are found to be in good agreement with the values found in previous analyses. 1 Introduction In the field of aerospace engineering it is imperative to be able to characterise a safe working lifetime for the various components. Inevitably, all components must be assumed to contain small micro-mechanical flaws which may initiate small cracks. Various loading regimes may cause these cracks to grow in size, until catastrophic failure occurs at some critical crack size. Linear elastic fracture mechanics postulates that the component lifetime is dependent upon a single quantity, the stress intensity factor (SIF), which is a function of geometry and of loading. This dependence has resulted in much work being done to find the SIF for various different geometries and loadings. A different stress intensity factor may be defined for each of the three basic modes of loading one symmetric, and two anti-symmetric modes denoted by KN where N = /,//or ///. A separate issue, which must be emphasised, is that many important structures are fully three dimensional, and admit no realistic two dimensional simplifications. Therefore, efficient methods to calculate the SIF for three dimensional bodies are of great benefit to the engineer. Traditionally, direct methods of calculating the stress intensity factor, such as from the crack opening displacement, rely heavily on accurate representation of the stress field near the crack front. This often proves difficult to achieve in practice, due to the unbounded nature of the stress field along the crack front. For this reason, interest has developed in using path-independent integrals, such as the J-mtegral developed by Rice [1]. Path-independent integrals offer the advantage that the elastic field values need not be calculated in the immediate
2 48 Boundary Element Technology vicinity of the crack front. However, the J-integral has the disadvantage that the three stress intensity factors are not given directly, but are coupled together. This necessitates the adoption of various special decoupling techniques. The presence of a third dimension further complicates this coupling of the SIFs given by the J-integral. Another approach which has gained currency over the last decade is that of weight functions, and the related topic of fundamental fields, as developed by Bueckner [2]. Fundamental fields are defined by Bueckner [2] to be elastic fields which satisfy the governing equations of elasticity theory, and also have a certain asymptotic behaviour on approaching the crack front. It should be noted that these fundamentalfields,although satisfying the elasticity equations, are unphysical, due to their unbounded strain energy along the crack front. The usual physical elastic field, with the standard asymptotic behaviour along the crack front, is termed a regularfield.unfortunately, there is a limitation to the direct applicability of fundamental fields to finite geometries, arising from the fact that they are defined only for an infinite elastic medium. Adopting the philosophy of both of the above approaches, it is possible to obtain a path-independent integral which uses directly the fundamental fields. A similar approach has been adopted before, in two dimensions, as for example by Stern et al. [4] and also by Zhang and Mikkola [3]. The approach followed here is derived from that used by Zhang and Mikkola in [3], and is a straightforward application of Betti's reciprocal theorem. 2 Derivation of the Path-Independent Integral Consider an arbitrary three dimensional cracked domain O, under two states of loading denoted by B and S in Figure may be a subregion of larger elastic R B State B State S Figure 2.1: The two loading states Figure 2.2: Point forces applied in state S domains R& and R^ which need not be identical, as long as 12 C R& and 12 C R^. State B is the actual physical loading state, consisting of a displacement field u, and a stressfieldresulting from a body force f, and boundary tractions t. Now define state 5 to be an auxiliary loading state in which 12 is considered to be cut from an infinite elastic body R^. The subregion 12 is subject to the restriction of the (Mode I) regular elastic field that acts throughout the infinite body; this regular field consists of stresses a' and displacements u'. Thus, 12
3 Boundary Element Technology 49 may be considered to be subject to boundary conditions u' (displacement), and F (traction), which are those values of u' and t' ((/ = a'n, where n is the normal to the boundary P) on the boundary I\ This regular field will first be used to derive a path independent integral, then later be subject to a limiting procedure, and its limit will be equated to a fundamental field (defined by Bueckner in [2], and known analytically) with a condensation point. The above cut-out procedure eliminates the difficulties caused by the fact that the fundamental field is defined only for an infinite elastic body. Since it is known that fundamentalfieldsmay be obtained by means of a limiting process (see for example [5]), the following procedure may be developed. Considerfirstthe loading of state S to consist of a doublet of equal and opposite opening point forces of magnitude F', one applied perpendicularly to each crack face at distance d from the crack front in a subregion 0 of an infinite elastic space RS, as in Figure 2.2. The displacement field in state S is denoted by u', and the stress field by </ (boundary tractions resulting from the stress field are represented by t'). Application of Betti's reciprocal theorem to the domain H bounded by a surface F, gives when both loading states are applied, F'UB + 1 1'. udf = j f. u'dft + / 1. ii'df, (2.1) r o r where UB is the crack opening displacement in state B at the position where the point force F' was applied in state 5. If d -> 0 and F' -» oo such that = RI = constant, (2.2) then the constant RI is known as Bueckner's field strength (by analogy to a similar constant in 2-D work, see for example [3]), and is defined in [6]. As d -» 0 the displacement UB (strictly UB is the difference in the displacements on the two crack faces) approaches the well-known limiting value on the crack front. Following Hartranft and Sih in [7], this is known to be given by the plane strain value: lim us = lim 2v(r = d,0 =,) = Kr, (2.3) where KI is the mode I stress intensity factor; E is Young's modulus; f/ is Poisson's ratio; and v is the theoretical displacement in the direction perpendicular to the crack face. Also, u' and t' approach the values for a fundamental field with a condensation point, denoted by u* and t* which can easily be derived for certain simple crack front configurations in the infinite cracked space R$. These are given by Bains, Aliabadi and Rooke [8] for the case of an rectilinear crack front in an infinite space, and by Sham and Zhou [9] for a penny shaped crack embedded in an infinite space. Without loss of generality, Bueckner's fundamental field strength may be chosen to be unity, that is,
4 50 Boundary Element Technology This gives, on neglecting body forces in state J9, (t.u* - t*. u)df. (2.5) ^ ' As the surface F is chosen arbitrarily, the above integral should be path independent; the above integral may then be evaluated for a given problem by any convenient method of quadrature. Note that, as implied above, the fields u* and t* should be chosen to be suitable for the problem at hand. For a problem in a finite body with a rectilinear crack front, it is natural to choose those fields derived for a rectilinear crack front. 3 Test Problems 3.1 Model I : Rectangular bar with a straight crack front The first test problem was that of a planar edge crack with a rectilinear crack front in a rectangular bar, which has a well established solution. In this test, the aim was to establish the validity of the integral for a simple well-known problem before extending the technique to more complex cases. The problem considered is that of a rectangular bar of width 6, thickness t and height 2/&, with a planar rectilinear edge crack of length a lying in the (z,z)-plane, as shown in Figure 3.1. The bar is subject to uniaxial tensile stress a along the y axis. In the example presented here the following ratios hold: a/6 = 0.5, t/b = 1.0 and h/b = 3.0. The numerical solution to the elastic boundary value problem was obtained using the commercial boundary element package BEASY [10]. The surface T was modelled using internal points, which were distributed to form a BEM-like discretisation of F into elements, and the stress and displacement field values were calculated at these internal points. The integral in equation (2.5) was then calculated using Gaussian quadrature over this discretisation of the surface F. The surface F was taken to be a sphere of radius r centred on the point P, the mid-point of the crack front; thus the SIF is obtained for that point. Results obtained using the above formulation are compared with the plane strain value obtained by Keer and Freedman [11], and are shown graphically in Figure 3.3. Note that the normalised SIF, Ki/(cr^/ a) is plotted. As can be seen, the agreement is acceptable, the difference being less than 2% for 0.1 < r/a < 0.8. It should be noted that for r/a > 0.6 there appears to be a slight divergence of the calculated stress intensity factor. This is probably due to the fact that as r/a increases, the contour F includes contributions from the surface- breaking region of the crack, where the stress intensity factor does not exist in the traditional sense. The source of the error for r/a < 0.1 is likely to be caused by the difficulty of calculatingfieldvalues accurately for internal points which lie close to the boundary.
5 Boundary Element Technology Model II: Penny-shaped crack in a large cylindrical bar This problem consists of a cylindrical bar of cross-sectional radius R and height 2/i, containing a penny-shaped crack of radius a lying in the (x,y) plane, as shown in Figure 3.2. The bar is subject to uniaxial tensile stress a directed a -y Figure 3.1: Bar with straight fronted crack Figure 3.2: Bar with pennyshaped crack along the z axis. In order to approximate a penny shaped crack in an infinite bar, the following ratios were taken, a/r = 0.1 and h/r = 3. The surface T was taken to be a sphere of radius r centred on the point P of the crack front, which lies on the positive x axis. Results obtained for the above model are again compared in Figure 3.4 with the results given by Cartwright and Rooke in [12]. As can be seen, the agreement is acceptable, differing by less than 2% with Reference [12] for r/a > Model III: Elliptical crack in a large cylindrical bar The problem consists of a cylindrical bar of radius R and height 2/t, containing an elliptical crack of semi-major axis a, and semi-minor axis 6. The bar is subject to uniaxial tensile stress a along the z axis. The following ratios were used for computation: a/r = Q.l,h/R = 3. Two aspect ratios were considered for the ellipse: b/a = 0.5 and b/a = 0.9. The contour surface F was taken to be a sphere, centred on the various points P,Q etc. of the crack front, where the SIF was to be calculated. The ratio used for the integration contour was r/a = 0.1. A technique known as approximation of fundamental fields was used to model the fundamental field of an elliptical crack using that of a penny shaped crack. The penny shaped crack chosen was that which had the same curvature as the elliptical crack at the point where the SIF was required. Results obtained for the above model for both b/a = 0.5 and b/a = 0.9 are given in Figure 3.6, compared with results from Rooke and Cartwright [12].
6 52 Boundary Element Technology 2.85 CO 1.08 ***** 192 nodes ##### 788 node* Plane Strain 1.04 ** ft) Current Work Rooke * Cartwright oi.oo r/a ratio i.oo r/a ratio i.oo Figure 3.3: Normalised Mode I SIF for straight fronted crack Figure 3.4: Normalised Mode I SIF for penny shaped crack The results are seen to be acceptable, being within 2% of Reference [12]. 4 Discussion and Conclusions 4.1 Sources of error It may be noted that generally, as r/a» 0, the solution becomes less accurate. The most likely explanation of this fact is that for small values of r/a, some of the internal points at which u and t are calculated lie very close to part of the external boundary (namely the crack surface). This gives rise to numerical difficulties in calculating the solution at these internal points. 4.2 Possible improvements in accuracy To utilise the technique of approximation of fundamental fields, it is desirable to have small r/a radius contour surfaces. Various options are available to improve the accuracy of the BEM internal point solution on these small radius contour surfaces. The most obvious approach is to use afinerbem mesh for the stress analysis, which allows more accurate internal point results for points somewhat closer to the boundary. However, this causes a large increase in the number of degrees of freedom of the BEM problem, and hence a large increase in solution time. An alternative approach is to use more sophisticated integration techniques in the calculation of elastic field values for those internal points that lie close to the boundary, for example the work of Cruse and Aithal [14], and also Mi and Aliabadi [15].
7 Boundary Element Technology i so GO QO Parametric angle of ellipse Figure 3.5: Elliptical crack in a large cylindrical bar Figure 3.6: Variation of SIF with angle along crack front for elliptical crack Acknowledgement This work was sponsored by the Defence Research Agency, Farnborough, U.K. British Crown Copyright 1995/DRA. Published with the permission of the Controller of Her Britannic Majesty's Stationery Office. References [1] J. R. Rice (1968), A path independent integral and the approximate analysis of strain concentration by notches and cracks. Trans. ASME, J. Appl. Mech., 35, [2] H. F. Bueckner (1987), Weight functions and fundamental fields for the penny-shaped and the half-plane crack in three-space, Int. J. Solids Structures, 23, [3] Z. L. Zhang and T. P. J. Mikkola (1992), A simple path-independent integral for calculating mixed mode stress intensity factors. Fatigue Fract. Engng Mater. Struct., 15, [4] M. Stern, E. B. Becker and R. S. Dunham (1976), A contour integral computation of mixed mode stress intensity factors, Int. J. Fract., 12, [5] H. F. Bueckner (1989), Observations on weight functions, Engng Anal. Boundary Elements, 6, 3-18
8 54 Boundary Element Technology [6]R. Bains, M. H. Aliabadi, and D. P. Rooke, Fracture mechanics weight functions in three dimensions: subtraction of fundamental fields, Int. J. Num. Meth. Engng., 35, [7] R. J. Hartranft and G. C. Sih (1977), Stress singularity for a crack with an arbitrarily curved crack front, Engng Frac. Mech., 9, [8] R. Bains, M. H. Aliabadi, and D. P. Rooke (1992), Fracture mechanics weight functions in three dimensions: subtraction of fundamental fields, Int. J. Num. Meth. Engng, 35, [9] T.-L. Sham and Zhou, Y. (1989), Computation of three-dimensional weight functions for circular and elliptical cracks, Int. J. Fract., 41, [10] BEASY, Boundary Element Analysis System Software, Computational Mechanics, Southampton, U.K. [11] L. M. Keer and J. M. Freedman (1973), Tensile strip with edge cracks, Int. J. Engng Sci., 11, [12] D. P. Rooke and D. J. Cartwright (1976), Compendium of Stress Intensity Factors, HMSO. [13] R. S. Bains, M. H. Aliabadi, D. P. Rooke (1992), Stress intensity factor weight functions for 3-D cracks infinitegeometries, in Advances in Boundary Element Methods for Fracture Mechanics, Eds. M. H. Aliabadi, C. A. Brebbia, Computational Mechanics Publications, Southampton, U.K. [14] T. A. Cruse and R. A. Aithal (1993) Non-singular boundary integral equation implementation, Int. J. Num. Meth. Engng., 36, [15] Y. Mi and M. H. Aliabadi (1994), Semi-analytical integration method for 3D near-hypersingular integrals, in Proc!(?** Int. Con}, on Boundary Element Methods, pp , Ed. C. A. Brebbia, Computational Mechanics Publications, Southampton, U.K
A modified quarter point element for fracture analysis of cracks
ndian Journal of Engineering & Materials Sciences Vol. 14, February 007, pp. 31-38 A modified quarter point element for fracture analysis of cracks Sayantan Paul & B N Rao* Structural Engineering Division,
More informationApplication of a non-local failure criterion to a crack in heterogeneous media S. Bavaglia*, S.E. Mikhailov*
Application of a non-local failure criterion to a crack in heterogeneous media S. Bavaglia*, S.E. Mikhailov* University of Perugia, Italy Email: mic@unipg.it ^Wessex Institute of Technology, Ashurst Lodge,
More informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
More informationCritical applied stresses for a crack initiation from a sharp V-notch
Focussed on: Fracture and Structural Integrity related Issues Critical applied stresses for a crack initiation from a sharp V-notch L. Náhlík, P. Hutař Institute of Physics of Materials, Academy of Sciences
More informationStress intensity factors under combined tension and torsion loadings
Indian Journal of Engineering & Materials Sciences Vol. 19, February 01, pp. 5-16 Stress intensity factors under combined tension and torsion loadings A E Ismail a *, A Ariffin b, S Abdullah b & M J Ghazali
More informationThermal load-induced notch stress intensity factors derived from averaged strain energy density
Available online at www.sciencedirect.com Draft ScienceDirect Draft Draft Structural Integrity Procedia 00 (2016) 000 000 www.elsevier.com/locate/procedia 21st European Conference on Fracture, ECF21, 20-24
More informationPROPAGATION OF CURVED CRACKS IN HOMOGENEOUS AND GRADED MATERIALS
PROPAGATION OF CURVED CRACKS IN HOMOGENEOUS AND GRADED MATERIALS Abstract Matthew T. Tilbrook, Robert J. Moon and Mark Hoffman School of Materials Science and Engineering University of New South Wales,
More informationCALCULATION OF FRACTURE MECHANICS PARAMETERS FOR AN ARBITRARY THREE-DIMENSIONAL CRACK, BY THE EQUIVALENT DOMAIN INTEGRAL METHOD 1
CALCULATION OF FRACTURE MECHANICS PARAMETERS FOR AN ARBITRARY THREE-DIMENSIONAL CRACK, BY THE EQUIVALENT DOMAIN INTEGRAL METHOD 1 G. P. NIKISHKOV 2 and S. N. ATLURI 3 Center for the Advancement of Computational
More informationAnalysis of Square-shaped Crack in Layered Halfspace Subject to Uniform Loading over Rectangular Surface Area
Copyright 2015 Tech Science Press CMES, vol.109-110, no.1, pp.55-80, 2015 Analysis of Square-shaped Crack in Layered Halfspace Subject to Uniform Loading over Rectangular Surface Area H. T. Xiao 1,2,3,
More informationModelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach
Modelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach Giang Dinh Nguyen, Chi Thanh Nguyen, Vinh Phu Nguyen School of Civil, Environmental and
More informationSize effect in the strength of concrete structures
Sādhanā Vol. 27 Part 4 August 2002 pp. 449 459. Printed in India Size effect in the strength of concrete structures B L KARIHALOO and Q Z XIAO Division of Civil Engineering School of Engineering Cardiff
More informationULTRASONIC REFLECTION BY A PLANAR DISTRIBUTION OF SURFACE BREAKING CRACKS
ULTRASONIC REFLECTION BY A PLANAR DISTRIBUTION OF SURFACE BREAKING CRACKS A. S. Cheng Center for QEFP, Northwestern University Evanston, IL 60208-3020 INTRODUCTION A number of researchers have demonstrated
More informationAn accelerated predictor-corrector scheme for 3D crack growth simulations
An accelerated predictor-corrector scheme for 3D crack growth simulations W. Weber 1 and G. Kuhn 2 1,2 1 Institute of Applied Mechanics, University of Erlangen-Nuremberg Egerlandstraße 5, 91058 Erlangen,
More informationINFLUENCE OF THE LOCATION AND CRACK ANGLE ON THE MAGNITUDE OF STRESS INTENSITY FACTORS MODE I AND II UNDER UNIAXIAL TENSION STRESSES
INFLUENCE OF THE LOCATION AND CRACK ANGLE ON THE MAGNITUDE OF STRESS INTENSITY FACTORS MODE I AND II UNDER UNIAXIAL TENSION STRESSES Najah Rustum Mohsin Southern Technical University, Technical Institute-Nasiriya,
More informationOn elastic compliances of irregularly shaped cracks
International Journal of Fracture 114: 245 257, 2002. 2002 Kluwer Academic Publishers. Printed in the Netherlands. On elastic compliances of irregularly shaped cracks I. SEVOSTIANOV 1 and M. KACHANOV 2
More informationDEVELOPMENT OF TEST GUIDANCE FOR COMPACT TENSION FRACTURE TOUGHNESS SPECIMENS CONTAINING NOTCHES INSTEAD OF FATIGUE PRE-CRACKS
Transactions, SMiRT-23 Division II, Paper ID 287 Fracture Mechanics and Structural Integrity DEVELOPMENT OF TEST GUIDANCE FOR COMPACT TENSION FRACTURE TOUGHNESS SPECIMENS CONTAINING NOTCHES INSTEAD OF
More informationL e c t u r e. D r. S a s s a n M o h a s s e b
The Scaled smteam@gmx.ch Boundary Finite www.erdbebenschutz.ch Element Method Lecture A L e c t u r e A1 D r. S a s s a n M o h a s s e b V i s i t i n g P r o f e s s o r M. I. T. C a m b r i d g e December
More informationTHE TRACTION BOUNDARY CONTOUR METHOD FOR LINEAR ELASTICITY
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng. 46, 1883}1895 (1999) THE TRACTION BOUNDARY CONTOUR METHOD FOR LINEAR ELASTICITY ZHOU SHENJIE*, CAO ZHIYUAN AND SUN
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationFinite element analysis of longitudinal debonding between fibre and matrix interface
Indian Journal of Engineering & Materials Sciences Vol. 11, February 2004, pp. 43-48 Finite element analysis of longitudinal debonding between fibre and matrix interface K Aslantaş & S Taşgetiren Department
More informationLaboratory 4 Bending Test of Materials
Department of Materials and Metallurgical Engineering Bangladesh University of Engineering Technology, Dhaka MME 222 Materials Testing Sessional.50 Credits Laboratory 4 Bending Test of Materials. Objective
More informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Analysis of singular and near-singular integrals by the continuation approach D. Rosen & D.E. Cormack Department of Chemical Engineering and Applied Chemistry, University of Toronto, Toronto M5S 1A4, Canada
More informationExtended Boundary Element Method approach for Direct and Accurate Evaluation of Stress Intensity Factors. Ibrahim Alatawi
Extended Boundary Element Method approach for Direct and Accurate Evaluation of Stress Intensity Factors Ibrahim Alatawi Thesis submitted towards the degree of Doctor of Philosophy Mechanics Research Group
More informationStress and Displacement Analysis of a Rectangular Plate with Central Elliptical Hole
Stress and Displacement Analysis of a Rectangular Plate with Central Elliptical Hole Dheeraj Gunwant, J. P. Singh mailto.dheerajgunwant@gmail.com, jitenderpal2007@gmail.com, AIT, Rampur Abstract- A static
More informationCrack Tip Plastic Zone under Mode I Loading and the Non-singular T zz -stress
Crack Tip Plastic Zone under Mode Loading and the Non-singular T -stress Yu.G. Matvienko Mechanical Engineering Research nstitute of the Russian Academy of Sciences Email: ygmatvienko@gmail.com Abstract:
More informationTransactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN
Unitary weight functions for semi-infinite crack in the half-plane K.L. Molski Institute of Cranes and Heavy Machinery Engineering, Warsaw University of Technology, Narbutta 85, 02-524 Warsaw, Poland ABSTRACT
More informationIntegral equations for crack systems in a slightly heterogeneous elastic medium
Boundary Elements and Other Mesh Reduction Methods XXXII 65 Integral equations for crack systems in a slightly heterogeneous elastic medium A. N. Galybin & S. M. Aizikovich Wessex Institute of Technology,
More informationJ. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, Bratislava, Slovakia
Evaluation of fracture parameters for functionally gradient materials J. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, 842 20 Bratislava, Slovakia
More informationStress intensity factor analysis for an interface crack between dissimilar isotropic materials
Stress intensity factor analysis for an interface crack between dissimilar isotropic materials under thermal stress T. Ikeda* & C. T. Sun* I Chemical Engineering Group, Department of Materials Process
More informationEfficient 2-parameter fracture assessments of cracked shell structures
Efficient 2-parameter fracture assessments of cracked shell structures B. Skallerud, K.R. Jayadevan, C. Thaulow, E. Berg*, K. Holthe Faculty of Engineering Science The Norwegian University of Science and
More informationElastic and Elastic-Plastic Behaviour of a Crack in a Residual Stress Field
Residual Stresses 2016: IC-10 Elastic and Elastic-Plastic Behaviour of a Crack in a Residual Stress Field Guiyi Wu a, Chris Aird b, David Smith and Martyn Pavier c* Department of Mechanical Engineering,
More informationInternational Journal of Innovative Research in Science, Engineering and Technology Vol. 2, Issue 7, July 2013
ANALYSIS AND MITIGATION OF STRESS CONCENTRATION FACTOR OF A RECTANGULAR ISOTROPIC AND ORTHOTROPIC PLATE WITH CENTRAL CIRCULAR HOLE SUBJECTED TO INPLANE STATIC LOADING BY DESIGN OPTIMIZATION Shubhrata Nagpal
More informationElastic behaviour of an edge dislocation near a sharp crack emanating from a semi-elliptical blunt crack
Chin. Phys. B Vol. 19, No. 1 010 01610 Elastic behaviour of an edge dislocation near a sharp crack emanating from a semi-elliptical blunt crack Fang Qi-Hong 方棋洪, Song Hao-Peng 宋豪鹏, and Liu You-Wen 刘又文
More informationCRACK INITIATION CRITERIA FOR SINGULAR STRESS CONCENTRATIONS Part I: A Universal Assessment of Singular Stress Concentrations
Engineering MECHANICS, Vol. 14, 2007, No. 6, p. 399 408 399 CRACK INITIATION CRITERIA FOR SINGULAR STRESS CONCENTRATIONS Part I: A Universal Assessment of Singular Stress Concentrations Zdeněk Knésl, Jan
More informationFRACTURE MECHANICS FOR MEMBRANES
FRACTURE MECHANICS FOR MEMBRANES Chong Li, Rogelio Espinosa and Per Ståhle Solid Mechanics, Malmö University SE 205 06 Malmö, Sweden chong.li@ts.mah.se Abstract During fracture of membranes loading often
More informationThree-dimensional thermo-mechanical analysis of layered elasticlplastic solids
Three-dimensional thermo-mechanical analysis of layered elasticlplastic solids W. Peng & Y.-T. Hsia Seagate Technology, U.S.A. Abstract At high temperature a layered solid undergoes intense thermal loading
More informationLinear Elastic Fracture Mechanics
Measure what is measurable, and make measurable what is not so. - Galileo GALILEI Linear Elastic Fracture Mechanics Krishnaswamy Ravi-Chandar Lecture presented at the University of Pierre and Marie Curie
More informationA novel approach to predict the growth rate of short cracks under multiaxial loadings
A novel approach to predict the growth rate of short cracks under multiaxial loadings F. Brugier 1&2, S. Pommier 1, R. de Moura Pinho 2, C. Mary 2 and D. Soria 2 1 LMT-Cachan, ENS Cachan / CNRS / UPMC
More informationFinite Element Investigation on the Stress State at Crack Tip by Using EPFM Parameters
Finite Element Investigation on the Stress State at Crack Tip by Using EPFM Parameters FRANCESCO CAPUTO, ALESSANDRO DE LUCA, GIUSEPPE LAMANNA 1, ALESSANDRO SOPRANO Department of Industrial and Information
More informationStress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension, compression or shear
Xi an 2-25 th August 217 Stress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension compression or shear Heyin Qi 1 Mingming Chen 2 Yonghong Duan 3 Daxu
More informationSurface crack subject to mixed mode loading
Engineering Fracture Mechanics 65 (2000) 1±14 www.elsevier.com/locate/engfracmech Surface crack subject to mixed mode loading M.Y. He a, J.W. Hutchinson b, * a Materials Engineering Department, University
More informationSchur decomposition in the scaled boundary finite element method in elastostatics
IOP Conference Series: Materials Science and Engineering Schur decomposition in the scaled boundary finite element method in elastostatics o cite this article: M Li et al 010 IOP Conf. Ser.: Mater. Sci.
More informationBIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS
BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS H. Kumazawa and T. Takatoya Airframes and Structures Group, Japan Aerospace Exploration Agency 6-13-1, Ohsawa, Mitaka,
More informationA CRITERION OF TENSILE FAILURE FOR HYPERELASTIC MATERIALS AND ITS APPLICATION TO VISCOELASTIC-VISCOPLASTIC MATERIALS
MTS ADHESIVES PROGRAMME 1996-1999 PERFORMANCE OF ADHESIVE JOINTS Project: PAJ1; Failure Criteria and their Application to Visco-Elastic/Visco-Plastic Materials Report 2 A CRITERION OF TENSILE FAILURE FOR
More informationTensile behaviour of anti-symmetric CFRP composite
Available online at www.sciencedirect.com Procedia Engineering 1 (211) 1865 187 ICM11 Tensile behaviour of anti-symmetric CFRP composite K. J. Wong a,b, *, X. J. Gong a, S. Aivazzadeh a, M. N. Tamin b
More informationElastic Crack Interaction Limit of Two Interacting Edge Cracks in Finite Body
Elastic Crack Interaction Limit of Two Interacting Edge Cracks in Finite Body R. Daud, M.A. Rojan Division of Applied Mechanics, School of Mechatronic Engineering, Pauh Putra Campus, Universiti Malaysia
More informationStresses and Displacements in Functionally Graded Materials of Semi-Infinite Extent Induced by Rectangular Loadings
Materials 2012, 5, 210-226; doi:10.3390/ma5020210 Article OPEN ACCESS materials ISSN 1996-1944 www.mdpi.com/journal/materials Stresses and Displacements in Functionally Graded Materials of Semi-Infinite
More informationTransactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
A direct formulation for three-dimensional elastoplastic boundary elements A. P. Cisilino*,**, M. H. Aliabadi* and J. L. Otegui** *Wessex Institute of Technology, As hurst Lodge, As hurst, Southampton
More informationComputer Methods in Applied Mechanics and Engineering, 66, (1988).
1 FINITE ELEMENT ANALYSIS OF STRESS CONCENTRATION AROUND A BLUNT CRACK IN A COSSERAT ELASTIC SOLID S. Nakamura and R.S. Lakes adapted from Computer Methods in Applied Mechanics and Engineering, 66, 257-266
More informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
Elastic-plastic model of crack growth under fatigue using the boundary element method M. Scibetta, O. Pensis LTAS Fracture Mechanics, University ofliege, B-4000 Liege, Belgium Abstract Life of mechanic
More informationConservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body
Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference
More informationTransactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN
Strain energy release rate calculation from 3-D singularity finite elements M.M. Abdel Wahab, G. de Roeck Department of Civil Engineering, Katholieke Universiteit, Leuven, B-3001 Heverlee, Belgium ABSTRACT
More informationDETERMINATION OF STRESS INTENSITY FACTORS ALONG CRACKED SURFACES IN PIPING ELBOWS STRUCTURES E.M.M.FONSECA*, F.Q.MELO**, R.A.F.
1 th PORTUGUESE CONFERENCE ON FRACTURE - 6 DETERMINATION OF STRESS INTENSITY FACTORS ALONG CRACKED SURFACES IN PIPING ELBOWS STRUCTURES E.M.M.FONSECA*, F.Q.MELO**, R.A.F.VALENTE** *CENUME IDMEC Instituto
More informationSingular asymptotic solution along an elliptical edge for the Laplace equation in 3-D
Singular asymptotic solution along an elliptical edge for the Laplace equation in 3-D Samuel SHANNON a, Victor PERON b, Zohar YOSIBASH a a Pearlstone Center for Aeronautical Engineering Studies, Dept.
More informationPROPAGATION OF A MODE-I CRACK UNDER THE IRWIN AND KHRISTIANOVICH BARENBLATT CRITERIA
Materials Science, Vol. 39, No. 3, 3 PROPAGATION OF A MODE-I CRACK UNDER THE IRWIN AND KHRISTIANOVICH BARENBLATT CRITERIA I. I. Argatov, M. Bach, and V. A. Kovtunenko UDC 539.3 We study the problem of
More informationA FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE
A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE Daxu Zhang and D. R. Hayhurst School of Mechanical, Aerospace and Civil Engineering,
More informationDetermination of Stress Intensity Factor for a Crack Emanating From a Rivet Hole and Approaching Another in Curved Sheet
International OPEN ACCESS Journal Of Modern Engineering Research (IJMER) Determination of Stress Intensity Factor for a Crack Emanating From a Rivet Hole and Approaching Another in Curved Sheet Raghavendra.
More information7.4 The Elementary Beam Theory
7.4 The Elementary Beam Theory In this section, problems involving long and slender beams are addressed. s with pressure vessels, the geometry of the beam, and the specific type of loading which will be
More informationMAE 323: Chapter 6. Structural Models
Common element types for structural analyis: oplane stress/strain, Axisymmetric obeam, truss,spring oplate/shell elements o3d solid ospecial: Usually used for contact or other constraints What you need
More informationFracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load
Fracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load Sriram Chintapalli 1, S.Srilakshmi 1 1 Dept. of Mech. Engg., P. V. P. Siddhartha Institute of Technology.
More informationFailure Modes and Bearing Capacity Estimation for Strip Foundations in C-ɸ Soils: A Numerical Study
Failure Modes and Bearing Capacity Estimation for Strip Foundations in C-ɸ Soils: A Numerical Study Paul Akagwu, Aaron Aboshio International Science Index, Civil and Environmental Engineering waset.org/publication/10005793
More informationSIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION Theodore V. Gortsas
More informationInfluence of impact velocity on transition time for V-notched Charpy specimen*
[ 溶接学会論文集第 35 巻第 2 号 p. 80s-84s (2017)] Influence of impact velocity on transition time for V-notched Charpy specimen* by Yasuhito Takashima** and Fumiyoshi Minami** This study investigated the influence
More informationACCURACY OF APPROXIMATION OF STRESS FIELD IN CRACKED BODIES FOR FAILURE ZONE EXTENT ESTIMATION
VIII International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-8 J.G.M. Van Mier, G. Ruiz, C. Andrade, R.C. Yu and X.X. Zhang (Eds) ACCURACY OF APPROXIMATION OF STRESS
More informationINVESTIGATION ON THE STRESS CONCENTRATION IN METALLIC FLAT PLATES DUE TO HOLES WITH DIFFERENT CONFIGURATIONS
International Journal of Mechanical Engineering and Technology (IJMET) Volume 8, Issue 7, July 2017, pp. 1718 1725, Article ID: IJMET_08_07_189 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=8&itype=7
More informationTransactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN
A computational method for the analysis of viscoelastic structures containing defects G. Ghazlan," C. Petit," S. Caperaa* " Civil Engineering Laboratory, University of Limoges, 19300 Egletons, France &
More informationFinite element modelling of infinitely wide Angle-ply FRP. laminates
www.ijaser.com 2012 by the authors Licensee IJASER- Under Creative Commons License 3.0 editorial@ijaser.com Research article ISSN 2277 9442 Finite element modelling of infinitely wide Angle-ply FRP laminates
More informationFatigue Damage Development in a Steel Based MMC
Fatigue Damage Development in a Steel Based MMC V. Tvergaard 1,T.O/ rts Pedersen 1 Abstract: The development of fatigue damage in a toolsteel metal matrix discontinuously reinforced with TiC particulates
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
A Boundary Element Method for analysis of bone remodelling M. Martinez, M.H. Aliabadi, H. Power PF&M&x TW/fwfe qf T^cAWogy, As/mr Ashurst, Southampton SO40 7AA, UK Abstract A boundary element formulation
More informationEFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES
EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES OF WOOD OR PLYWOOD CONSIDERED AS ORTHOTROPIC MATERIALS Information Revied and Reaffirmed March 1956 No. 1510 EFFECT OF ELLIPTIC
More informationTransactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN
Evaluation of approximate methods for elastic-plastic analysis of notched bodies S. S0rb0, G. Harkegard The Norwegian Institute of Technology, Trondheim and Power Generation, Baden, Switzerland ABB ABSTRACT
More informationFracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach
University of Liège Aerospace & Mechanical Engineering Fracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach Ludovic Noels Computational & Multiscale Mechanics of
More informationEffective stress assessment at rectangular rounded lateral notches
Focussed on characterization of crack tip fields Effective stress assessment at rectangular rounded lateral notches Enrico Maggiolini, Roberto Tovo, Paolo Livieri University of Ferrara Enrico.maggiolini@unife.it,
More informationPGroupN background theory
12/12/03 Dr Francesco Basile, Geomarc Ltd PGroupN background theory Estimation of the deformations and load distributions in a group of piles generally requires the use of computer-based methods of analysis.
More informationAn exact reanalysis algorithm using incremental Cholesky factorization and its application to crack growth modeling
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 01; 91:158 14 Published online 5 June 01 in Wiley Online Library (wileyonlinelibrary.com). DOI: 10.100/nme.4 SHORT
More informationStress magnification due to stretching and bending of thin ligaments between voids
Arch Appl Mech 6) DOI 1.17/s419-6-8-6 ORIGINAL V. A. Lubarda X. Markenscoff Stress magnification due to stretching and bending of thin ligaments between voids Received: 8 December 5 Springer-Verlag 6 Abstract
More informationJ-Integral Evaluation of Surface Cracks in Round Bar under Mode III Loadings
Research Journal of Applied Sciences, Engineering and Technology 7(10): 1985-1993, 2014 ISSN: 2040-7459; e-issn: 2040-7467 Maxwell Scientific Organization, 2014 Submitted: June 17, 2013 Accepted: June
More informationExamination in Damage Mechanics and Life Analysis (TMHL61) LiTH Part 1
Part 1 1. (1p) Define the Kronecker delta and explain its use. The Kronecker delta δ ij is defined as δ ij = 0 if i j 1 if i = j and it is used in tensor equations to include (δ ij = 1) or "sort out" (δ
More informationFINITE GEOMETRY CORRECTION FACTORS FOR THE STRESS FIELD AND STRESS INTENSITIES AT TRASNVERSE FILLET WELDS
Riggenbach i FINITE GEOMETRY CORRECTION FACTORS FOR THE STRESS FIELD AND STRESS INTENSITIES AT TRASNVERSE FILLET WELDS by KANE RYAN RIGGENBACH Submitted in partial fulfillment of the requirements For the
More information9 MECHANICAL PROPERTIES OF SOLIDS
9 MECHANICAL PROPERTIES OF SOLIDS Deforming force Deforming force is the force which changes the shape or size of a body. Restoring force Restoring force is the internal force developed inside the body
More informationAssessment of boundaryelement method for modelling the structural response of a pavement
Assessment of boundaryelement method for modelling the structural response of a pavement (PhD Student, Lecturer DEC of FCTUC) (Full Professor - IST) Pavements: Materials, design and performance a, LNEC,
More informationFAILURE ASSESSMENT DIAGRAM ASSESSMENTS OF LARGE-SCALE CRACKED STRAIGHT PIPES AND ELBOWS
Transactions, SMiRT-23, Paper ID 093 FAILURE ASSESSMENT DIAGRAM ASSESSMENTS OF LARGE-SCALE CRACKED STRAIGHT PIPES AND ELBOWS R A Ainsworth 1, M Gintalas 1, M K Sahu 2, J Chattopadhyay 2 and B K Dutta 2
More informationULTIMATE STRENGTH OF SQUARE PLATE WITH RECTANGULAR OPENING UNDER AXIAL COMPRESSION
Journal of Naval Architecture and Marine Engineering June, 2007 http://jname.8m.net ULTIMATE STRENGTH OF SQUARE PLATE WITH RECTANGULAR OPENING UNDER AXIAL COMPRESSION M. Suneel Kumar 1*, P. Alagusundaramoorthy
More informationKeyword Construction Science, Normal Stress, Curved Beams, Hook, Neutral Axis.
International Journal of Innovative Research in Engineering & Management (IJIREM) ISSN: 2350-0557, Volume-3, Issue-1, January-2016 Theoretical Demonstration of the Normal Stress Formula in the Hook Msc.
More informationTransactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN
On the Computation of Elastic Elastic Rolling Contact using Adaptive Finite Element Techniques B. Zastrau^, U. Nackenhorst*,J. Jarewski^ ^Institute of Mechanics and Informatics, Technical University Dresden,
More informationANALYSIS OF SHELL-LIKE STRUCTURES BY THE BOUNDARY ELEMENT METHOD BASED ON 3-D ELASTICITY: FORMULATION AND VERIFICATION
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, VOL. 41, 541 558 (1998) ANALYSIS OF SHELL-LIKE STRUCTURES BY THE BOUNDARY ELEMENT METHOD BASED ON 3-D ELASTICITY: FORMULATION AND VERIFICATION
More informationStress intensity factors determination for an inclined central crack on a plate subjected to uniform tensile loading using FE analysis
Stress intensity factors determination for an inclined central crack on a plate subjected to uniform tensile loading using FE analysis Paulo C. M. Azevedo IDMEC and Faculdade de Engenharia da Universidade
More informationSIZE EFFECTS IN THE COMPRESSIVE CRUSHING OF HONEYCOMBS
43rd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Con 22-25 April 2002, Denver, Colorado SIZE EFFECTS IN THE COMPRESSIVE CRUSHING OF HONEYCOMBS Erik C. Mellquistand Anthony M.
More informationOn the Path-Dependence of the J-Integral Near a Stationary Crack in an Elastic-Plastic Material
Dorinamaria Carka Chad M. Landis e-mail: landis@mail.utexas.edu Department of Aerospace Engineering and Engineering Mechanics, University of Texas at Austin, 10 East 4th Street, C0600 Austin, TX 7871-035
More information1 Static Plastic Behaviour of Beams
1 Static Plastic Behaviour of Beams 1.1 Introduction Many ductile materials which are used in engineering practice have a considerable reserve capacity beyond the initial yield condition. The uniaxial
More informationMaterial Subjected to Uniaxial Tens
NAOSITE: Nagasaki University's Ac Title Author(s) Citation Experimental and Numerical Research Material Subjected to Uniaxial Tens Yang, Lei; Jiang, Yujing; Li, Shuca Journal of Geotechnical and Geoenvi
More informationEFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A SPHERE AGAINST A RIGID FLAT A FINITE ELEMENT STUDY
Proceedings of the International Conference on Mechanical Engineering 2009 (ICME2009) 26-28 December 2009, Dhaka, Bangladesh ICME09- EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A
More informationAbstract. Introduction
Structural vibration isolation by rows of piles S.E. Kattis", D. Polyzos*, D.E. Beskos* "Department of Mechanical Engineering, University ofpatras, G^-26 j 00 f6zrr^ Greece ^Department of Civil Engineering,
More informationAnalysis of Open Circular Reinforced Concrete Shell Roof with Curved Edge Simply Supported
Proceedings of IOE Graduate Conference, 2017 Volume: 5 ISSN: 2350-8914 (Online), 2350-8906 (Print) Analysis of Open Circular Reinforced Concrete Shell Roof with Curved Edge Simply Supported Abhishek Adhikari
More informationImproved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method
Center for Turbulence Research Annual Research Briefs 2006 313 Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method By Y. Khalighi AND D. J. Bodony 1. Motivation
More informationNumerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach
Numerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach S. Stefanizzi GEODATA SpA, Turin, Italy G. Barla Department of Structural and Geotechnical Engineering,
More informationINTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011
Interlaminar failure analysis of FRP cross ply laminate with elliptical cutout Venkateswara Rao.S 1, Sd. Abdul Kalam 1, Srilakshmi.S 1, Bala Krishna Murthy.V 2 1 Mechanical Engineering Department, P. V.
More informationBack Analysis of Measured Displacements of Tunnels
Rock Mechanics and Rock Engineering 16, 173--180 (1983) Rock Mechanics and Rock Engineering 9 by Springer-Verlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe
More informationExperimental validation of a fracture plane model for multiaxial cyclic loading
Experimental validation of a fracture plane model for multiaxial cyclic loading A. Carpinteri, R. Brighenti, A. Spagnoli Dipartimento di Ingegneria Civile, Universita di Parma Viale delle Scienze, 431
More informationModeling of Interphases in Fiber-Reinforced Composites Under Transverse Loading Using the Boundary Element Method
Y. J. Liu 1 Assistant Professor, e-mail: Yijun.Liu@uc.edu Mem. ASME N. Xu Graduate Student Department of Mechanical, Industrial, and Nuclear Engineering, P.O. Box 210072, University of Cincinnati, Cincinnati,
More information