A beam theory fracture mechanics approach for strength analysis of beams with a hole

Similar documents
COMPARISON OF THERMAL HAND MODELS AND MEASUREMENT SYSTEMS

Resistivity survey at Stora Uppåkra, Sweden

Sound absorption properties of a perforated plate and membrane ceiling element Nilsson, Anders C.; Rasmussen, Birgit

Analytical Studies of the Influence of Regional Groundwater Flow by on the Performance of Borehole Heat Exchangers

Transmuted distributions and extrema of random number of variables

Aalborg Universitet. Lecture 14 - Introduction to experimental work Kramer, Morten Mejlhede. Publication date: 2015

Citation for published version (APA): Burcharth, H. F. (1988). Stress Analysis. Aalborg: Department of Civil Engineering, Aalborg University.

Optimizing Reliability using BECAS - an Open-Source Cross Section Analysis Tool

Lidar calibration What s the problem?

Design possibilities for impact noise insulation in lightweight floors A parameter study

Graphs with given diameter maximizing the spectral radius van Dam, Edwin

Resistivity imaging for mapping of quick clays

Matrix representation of a Neural Network

A Note on Powers in Finite Fields

Modelling Salmonella transfer during grinding of meat

The dipole moment of a wall-charged void in a bulk dielectric

Numerical modeling of the conduction and radiation heating in precision glass moulding

Reflection of sound from finite-size plane and curved surfaces

Radioactivity in the Risø District July-December 2013

Sorption isotherms for textile fabrics and foam used in the indoor environment

Lattice Boltzmann Modeling From the Macro- to the Microscale - An Approximation to the Porous Media in Fuel Cells -

Reducing Uncertainty of Near-shore wind resource Estimates (RUNE) using wind lidars and mesoscale models

On a set of diophantine equations

Geoelectrical and IP Imaging Used for Pre-investigation at a Tunnel Project. Danielsen, Berit Ensted; Arver, Henrik; Karlsson, T; Dahlin, Torleif

Crack Tip Flipping: A New Phenomenon yet to be Resolved in Ductile Plate Tearing

Prediction of airfoil performance at high Reynolds numbers.

Radioactivity in the Risø District January-June 2014

Benchmark of Femlab, Fluent and Ansys

An Equivalent Source Method for Modelling the Lithospheric Magnetic Field Using Satellite and Airborne Magnetic Data

A note on non-periodic tilings of the plane

The Wien Bridge Oscillator Family

Geometry explains the large difference in the elastic properties of fcc and hcp crystals of hard spheres Sushko, N.; van der Schoot, P.P.A.M.

The M/G/1 FIFO queue with several customer classes

Dynamic Effects of Diabatization in Distillation Columns

De Prisco, Renato; Comunian, Michele; Grespan, Franscesco; Piscent, Andrea; Karlsson, Anders

A comparison study of the two-bladed partial pitch turbine during normal operation and an extreme gust conditions

Reflection and absorption coefficients for use in room acoustic simulations

A numerical study of vortex-induced vibrations (viv) in an elastic cantilever

Embankment dam seepage assessment by resistivity monitoring

Simulating wind energy resources with mesoscale models: Intercomparison of stateof-the-art

GIS Readiness Survey 2014 Schrøder, Anne Lise; Hvingel, Line Træholt; Hansen, Henning Sten; Skovdal, Jesper Høi

Realtime 3D stress measurement in curing epoxy packaging

Tilburg University. Strongly Regular Graphs with Maximal Energy Haemers, W. H. Publication date: Link to publication

Aalborg Universitet. Land-use modelling by cellular automata Hansen, Henning Sten. Published in: NORDGI : Nordic Geographic Information

Accelerating interior point methods with GPUs for smart grid systems

Spatial decay rate of speech in open plan offices: the use of D2,S and Lp,A,S,4m as building requirements Wenmaekers, R.H.C.; Hak, C.C.J.M.

Estimation of Peak Wave Stresses in Slender Complex Concrete Armor Units

Kommerciel arbejde med WAsP

Aalborg Universitet. All P3-equipackable graphs Randerath, Bert; Vestergaard, Preben Dahl. Publication date: 2008

Multi (building)physics modeling

Determination of packaging induced 3D stress utilizing a piezocoefficient mapping device

Higher-order spectral modelling of the diffraction force around a vertical circular cylinder

Global Wind Atlas validation and uncertainty

Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)

Geophysical techniques for investigating shallow and deep landslides: results in the framework of PREVIEW project (FP6 UE)

A systematic methodology for controller tuning in wastewater treatment plants

Validation of Boundary Layer Winds from WRF Mesoscale Forecasts over Denmark

Tectonic Thinking Beim, Anne; Bech-Danielsen, Claus; Bundgaard, Charlotte; Madsen, Ulrik Stylsvig

HIF-1 can act as a tumor suppressor gene in murine Acute Myeloid Leukemia.

Optimal acid digestion for multi-element analysis of different waste matrices

Measurements of the angular distribution of diffuse irradiance

Wake decay constant for the infinite wind turbine array Application of asymptotic speed deficit concept to existing engineering wake model

Impact of the interfaces for wind and wave modeling - interpretation using COAWST, SAR and point measurements

Online algorithms for parallel job scheduling and strip packing Hurink, J.L.; Paulus, J.J.

On machine measurements of electrode wear in micro EDM milling

Development of a Dainish infrastructure for spatial information (DAISI) Brande-Lavridsen, Hanne; Jensen, Bent Hulegaard

Citation for pulished version (APA): Henning, M. A., & Yeo, A. (2016). Transversals in 4-uniform hypergraphs. Journal of Combinatorics, 23(3).

Phase Equilibrium in Amino Acid Salt Systems for CO2 Capture

The colpitts oscillator family

Aalborg Universitet. Published in: IEEE International Conference on Acoustics, Speech, and Signal Processing. Creative Commons License Unspecified

Design of Robust AMB Controllers for Rotors Subjected to Varying and Uncertain Seal Forces

The IEA Annex 20 Two-Dimensional Benchmark Test for CFD Predictions

On the Einstein-Stern model of rotational heat capacities

Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)

Generating the optimal magnetic field for magnetic refrigeration

Alkali Release from Typical Danish Aggregates to Potential ASR Reactive Concrete

Minimum analysis time in capillary gas chromatography. Vacuum- versus atmospheric-outlet column operation Leclercq, P.A.; Cramers, C.A.M.G.

On the Unique D1 Equilibrium in the Stackelberg Model with Asymmetric Information Janssen, M.C.W.; Maasland, E.

Solar radiation and thermal performance of solar collectors for Denmark

Room Allocation Optimisation at the Technical University of Denmark

Steady State Comparisons HAWC2 v12.2 vs HAWCStab2 v2.12

On numerical evaluation of two-dimensional phase integrals

Application of mesoscale models with wind farm parametrisations in EERA-DTOC

Comparisons of winds from satellite SAR, WRF and SCADA to characterize coastal gradients and wind farm wake effects at Anholt wind farm

Organ doses determined using a RANDO phantom for different radionuclide depositions and photon energies

Tilburg University. Two-Dimensional Minimax Latin Hypercube Designs van Dam, Edwin. Document version: Publisher's PDF, also known as Version of record

The calculation of stable cutting conditions in turning

Thermodynamics of monomer partitioning in polymer latices: effect of molar volume of the monomers

Characterization of microporous silica-based membranes by calorimetric analysis Boffa, Vittorio; Yue, Yuanzheng

Vindkraftmeteorologi - metoder, målinger og data, kort, WAsP, WaSPEngineering,

Benchmarking of optimization methods for topology optimization problems

Notes on cooperative research in the measurement of Gottwein-temperature Veenstra, P.C.

Diffusion of sulfuric acid in protective organic coatings

Strongly 2-connected orientations of graphs

Exposure Buildup Factors for Heavy Metal Oxide Glass: A Radiation Shield

Validation and comparison of numerical wind atlas methods: the South African example

Aalborg Universitet. Seawall Overtopping Tests sea erosion study : Sohar Cornice Frigaard, Peter. Publication date: 1990

Convergence analysis for a class of LDPC convolutional codes on the erasure channel

Analysis of Biaxially Stressed Bridge Deck Plates

Polydiagnostic study on a surfatron plasma at atmospheric pressure

Transcription:

A beam theory fracture mechanics approach for strength analysis of beams with a hole Danielsson, Henrik; Gustafsson, Per-Johan Published in: International Network on Timber Engineering Research - Proceedings Meeting 48 Published: 5-8- Document Version Publisher's PDF, also known as Version of record Link to publication Citation for published version (APA): Danielsson, H., & Gustafsson, P-J. (5). A beam theory fracture mechanics approach for strength analysis of beams with a hole. In International Network on Timber Engineering Research - Proceedings Meeting 48 [INTER / 48-9-] Karlsruhe, Germany. General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal Take down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. L UNDUNI VERS I TY PO Box7 L und +4646

Download date:. Jul. 8

Abeamtheoryfracturemechanics approachforstrengthanalysisof beamswithahole HenrikDanielsson,DivisionofStructuralMechanics,LundUniversity,Sweden PerJohanGustafsson,DivisionofStructuralMechanics,LundUniversity,Sweden Keywords:Glulam,Hole,FractureMechanics,StrengthAnalysis. Introduction INTER / 48-9 - Aholeinamemberconstitutesasuddenchangeinthecrosssectionandconcen tratedperpendiculartograintensilestressandshearstressappearinthevicinityof thehole,seefigure..thisstresssituationmayforrelativelylowexternalloads leadtocrackpropagationinthefiberdirection.lookingatthedesignapproachesfor beamswithaholeintimberengineeringdesigncodesoverthelastdecades,itcanbe seenthattheissuehasbeentreatedinmanydifferentways.thetheoreticalback groundsonwhichthedesignapproachesarebasedshowfundamentaldifferences andtherearealsomajordiscrepanciesbetweenthestrengthpredictionsaccordingto thedifferentcodesaswellasbetweentestsandpredictionsaccordingtocodes (Aicher&Höfflin4,Höfflin5,Danielsson&Gustafsson8,Danielsson& Gustafsson). Thereisatthemomentnofullyaccepteddesignmethodbasedonacompletelyra tionalmechanicalbackground.thereareforexamplenodesignequationsforbeams withaholeinthecontemporaryversionofeurocode5(4).however,design equationsarefoundinthegermanandaustriannationalannexestoec5.thisde signapproachoriginatesfromtheworkpresentedbykolbandepple(985),although simplificationsandempiricalmodificationshavebeenaddedovertime. + - - + 335

Forapplicationswithconcentratedperpendiculartograintensilestressandshear stress,conventionalmaximumstressfailurecriteriaareseldomofanyusebutin steadstrengthanalysisbymeansoffracturemechanicsismorerelevant.theec5 designequationsforendnotchedbeamsanddoweltypeconnectionsloadedatan angletograinareexamplesoffracturemechanicsbaseddesignequations. Concerningthecaseofendnotchedbeams,anequationfortheenergyreleaserate duringcrackextensionandthecorrespondingbeamstrengthwasderivedbycompli anceanalysisusingbeamtheoryalmost3yearsago(gustafsson988).initsec5 implementation,thegustafssondesigncriterionisexpressedasacomparisonofa nominalshearstressandthereducedshearstrength,althoughthedecisivematerial propertiesfromthefracturemechanicsapproacharefractureenergyandstiffness. Inspiredbytherelativesimplicityofthedesignequationforendnotchedbeams,a similarbutgeneralizedmethodwasoutlinedbygustafsson(5).thisgeneralized methodisinthepresentpaperindetaildevelopedandpresentedforthecaseof beamswithahole.theapplicationofabeamtheoryfracturemechanicsapproach forbeamswithaholeismorecomplexthanforendnotchedbeamsdueto: thesignificantinfluenceofshearmakesitnecessarytoconsidertherespective contributionsofmodesiandiitothetotalenergyreleaserate, thenormalforceactingonthecrosssectionpartsaboveandbelowtheholemust beconsideredand thecrosssectionforcesandmomentsactingonthepartsaboveandbelowthe holearestaticallyindeterminate. Inthispaper,abeamtheoryfracturemechanicsapproachforstrengthanalysisof beamswithaholeispresented.thisapproachisbasedonassumptionsofanortho tropicmaterialbehavior,abeammodelaccordingtotimoshenkobeamtheoryanda mixedmodefracturecriterionbasedonlinearelasticfracturemechanics(lefm). Thestrengthanalysisapproachisevaluatedbymeansofacomparisonbetweenthe oreticallypredictedbeamstrengthsandstrengthsfoundfromexperimentaltests.for thispurpose,testresultsforbeamswithcircularholes(höfflin5,aicher&höfflin 6)areusedaswellastestresultsofbeamswithsquareholes(Danielsson8). Thetestsofbeamswithsquareholesarefurtheralsopresentedanddiscussedin (Danielsson&Gustafsson8)and(Danielsson&Gustafsson). Furtherevaluationiscarriedoutbymeansofcomparisontootherstrengthanalysis methods,includingbothcodetypemethodsandmoregeneralmethodsbasedonlin earandnonlinearfracturemechanicsapproachescarriedoutusingdand3dfinite element(fe)analysis. 336

Theoryandmodelformulation Thebasicconceptofthepresentapproachforanalysisofbeamswithaholerelates torectangular/squareholesandbyassumptionalsocircularholesareanalyzed.cross sectionforcesandmomentsinthebeampartsaboveandbelowtheholearedeter minedbyequilibriumconsideringkinematicassumptionsaccordingtotimoshenko beamtheory.aninfinitesimallyshortpartofthebeamattheendoftheholeisthen considered.theforcesandmomentactingacrossahorizontalsectionoftheinfinites imallyshortpartofthebeam,dividingitintoonepartaboveandonepartbelowan assumedcrack,aredeterminedbyequilibrium.theenergyreleaseratesformodesi andiiarethenobtainedbyusingthemethodofworkofcrackclosurewithconsidera tiontothedeformationsoftheinfinitesimalpartsaboveandbelowtheassumed crack.thebeamstrengthwithrespecttocrackingisthenfoundbyusingamixed modefracturecriterion. Inordertofacilitateaconvenientformulation,whichisconsistentforcircularand rectangular/squareholeswithorwithoutroundedcorners,anumberofassumptions andsimplificationsareintroduced.thesearepartlybasedonengineeringapproxima tionsandpartlybasedonexperiencefromexperimentaltests.. Basicassumptionsanddefinitionofgeometry Definitionsandnotationfortheparametersusedtodescribethebeamandthehole geometryarefoundinfigure..acircularholewithdiametermayformallybe regardedasarectangular/squareholewith= =.Theparameterdefinesthe holeplacementwithrespecttobeamheightandisgivenbythecoordinateofthe centreofthehole.,andrefertovaluesattheedgeoftheholeonitsrighthand side,althoughthefigurebelowmayseemtoindicateanotherlocation. Theformulationofthepresentapproachisgeneralinthesensethatitallowsforany combinationofcrosssectionalforces,and.dependingonloadingconditions andthesignsof,and,differentareasintheholevicinityareexposedtoper pendiculartograintensilestressandmayexperiencecracking,seefigure..the mostcommoncaseforpracticalengineeringpurposesishoweverbelievedtobe transversalloadinggivingacombinationofshearforceandbendingmomentonly. a INTER / 48-9 - h/ h/ y x r s h u h d h l M V N h b 337

Thelocationofthepointofcrackinitiationisassumedtobeknownaprioriandcrack propagationisassumedtooccurintheparalleltograindirection,whichisassumed tobealignedwiththebeamlengthdirection(thedirection).thechosenpredefined cracklocationisbasedonresultsfromtheexperimentaltestsmentionedaboveof beamsloadedinbending,havingeitheracircularorasquarehole.thetestswith squareholeshavingroundedcorners(. )showedthatcrackingcommonly startedatthetoprightcornerofthehole.crackinitiationcommonlyoccurredwithin theroundedpartoftheholecorner.thetestswithcircularholesshowedcrackinitia tionattheholeperipheryatanangleofapproximately45 withrespecttothebeam centralaxis,forholesplacedcentricallywithrespecttothebeamheight.basedon thesefindingsfromexperimentaltests,thelocationofthecrackisinthepresentap proachassumedtobeatanangleof45 fromthebeamlengthdirectionaccordingto Figure.a),yieldingadistancefromtheupperhorizontalholeedgetothecrack planeof.3. Thegeometryofthebeampartsaboveandbelowtheholeareinthecalculations simplifiedinthesensethattheseareassumedtobeprismatic.asillustratedinfigure.b)e),thereareanumberoffeasibleinterpretationsoftheholegeometrywhich resultinthissimplification.theinterpretationusedforallanalysespresentedhereis accordingtofigure.b).interpretationofholegeometryaccordingtofigures. c)d)orsomeothersimplificationmayhoweverpossiblyyieldequalorbetteragree mentbetweenmodelstrengthpredictionsandexperimentallyfoundbeamstrength. a) b) r 45 o.3r assumed crack locaon r 45 o c) d) e) 338

. Elementcrosssectionalforcesandmoments Thecrosssectionalforces, and (=,)inthebeampartsabove()andbe low()theholearedeterminedbyequilibriumconsiderationsofthestaticallyinde terminatesystemshowninfigure.3.kinematicassumptionsaccordingtotimo shenkobeamtheoryareusedforbeamelementsand.theverticalcrosssections onthetwosidesoftheholeareassumedtoremainplaneatloading.withthesimpli ficationgiveninfigure.b),thebeampartsaboveandbelowtheholeareassumed tohaveconstantheightsaccordingto () () Toaccountfortheelasticclampingofthebeamelementsinanapproximateway, theseareinthecalculationsgivenalengthslightlylongerthantheactuallengthof thehole,=+.5,basedonnumericalresultspresentedbypetersson(974). L e h M e N h V Accordingtothekinematicassumptions,theinfinitesimallythincrosssectiontothe rightoftheholeremainsplanepriortocrackinitiationandpropagation.atcrack propagationhowever,itissplitintotwocrosssectionsaccordingtofigure.4.the crosssectionsectionalforces, and forcingthecrosssectiontoremain planeintheuncrackedstate maybedeterminedfromequilibriumgiving (3) (4) (5) where, and aretheresultingforcesandthemomentfromthenormal andshearstressesactingontherightsideofthecrosssectionpartabovethecrack planeaccordingto (6) (7) (8) 339

h N M V crack locaon h N M V V M N N s M s V s central axis.3 Relativedisplacementsatcrackpropagation Atcrackinitiationandpropagationinthebeamlengthdirection,theoriginalbeam crosssectionissplitintotwoseparatecrosssections.thetworespectivecrosssec tionsarestillassumedtoremainplanebutsincetheyareseparated,theymaybede formedindifferentways.therelativedisplacementsandandtherelativerota tionofthetwocrosssections,aboveandbelowthecrackplane,areillustratedin Figure.5.Therelativedisplacement,relatedto,correspondstomodeIcrack deformation.therelativedisplacementandtherelativerotationareinfluenced byboth and withcorrespondingtomodeiicrackdeformationandassumed tocontributetomodeicrackdeformation. dx aer crack propagaon before crack propagaon -v -θ h crack locaon -u h-h Therelationsbetweentherelativedisplacementsandandtherelativerotation andthecrosssectionsectionalforcesandmoment, and maybeexpressedas (9) wherethedisplacementandforcevectorsaregivenby () () andthecompliancematrixisgivenby 34

() Thecomponentsofthecompliancematrixaredeterminedbasedonthecompliance ofthecrosssectionsofinfinitesimallengthaboveandbelowthecrackplane.consid eringkinematicassumptionsaccordingtotimoshenkobeamtheorygives (3) (4) (5) (6) where istheparalleltograinstiffnessand istheshearstiffness..4 Crackclosurework,energyreleaserateandfracturecriterion AccordingtoLEFM,seee.g.textbook(Hellan985),crackpropagationisgovernedby crackpropagationcriteriawhichmaybeformulatedbasedon,forexample,thecon ceptofstrainenergyreleaserateorstressintensityfactors.thestrainenergyrelease ratecan,accordingtolefm,alsobecalculatedintermsof,i.e. theworkrequiredtocompletelycloseapropagatedcrack.thecrackclosurework mayforthepresentapplicationbeexpressedas (7) where and refertothecrackclosureworkrelatedtomodeiandiirespec tively.indices,andrefertotherelativedisplacementsandtherotationillus tratedinfigure.5above.accordingtolefm,thestressintensityfactors and areproportionaltotheappliedload.theprincipleofsuperpositionmaybeusedfor multipleloadcasesgivingcontributionsinthesamemodeofdeformation.forthe presentapplicationthismeansthatthemodeistressintensityfactormaybeex pressedas (8) TherelationshipbetweenthemodeIandmodeIIstressintensityfactors andthe correspondingenergyreleaserates (=I,II)isgivenby (9) where (=I,II)isameasureofthestiffnessofthematerialwithrespecttothecor respondingmodeofdeformation(gustafsson). 34

Theenergyreleaserateatcrackextensionoveranarea=isequaltothecor respondingcrackclosureworkgivingforthepresentcase () and () () (3) (4) Inordertodeterminethebeamstrengthwithrespecttocracking,acrackpropaga tioncriterionneedstobechosen.inthepresentapplications,withingeneralmixed modeofloading,thefollowinginteractioncriterionhasbeenused (5) where and aretheenergyreleaseratesinmodeiandiigivenaboveandwhere and arethecorrespondingcriticalenergyreleaserates(orfractureenergies). ResultspresentedinSections4and5,relatingtoanalysisofbeamswithaholeusing thepresentapproach,arebasedon==.5andmaterialpropertyparametersac cordingtotable.thechoiceofcrackpropagationcriterionisnotobvious,andother criteriacouldpossiblybemoresuitable.incaseofnegativevaluesof and/or, thenegativecontributionisignoredandfractureinpuremodeioriiisconsidered. Parameter Notation Value Paralleltograinstiffness MPa Shearstiffness 6MPa Fractureenergy,modeI 3Nm/m Fractureenergy,modeII 9Nm/m 3 Comparisontoendnotchbeamequation Thepresentanalysisapproachforbeamswithaholehasmanyfeaturesincommon withthelefmbaseddesignapproachfornotchedbeams,originallypresentedby Gustafsson(988)andwithmodificationsincludedinEC5.Thepresentapproachfor strengthanalysisofbeamswithaholemayalsobeusedforanalysisofnotched beams,sincethisisessentiallyonlyaspecialcaseintermsofthegeneralgeometry illustratedinfigure.. Asmentionedintheintroduction,therearehoweversomedifferencesinthegeneral formulationbetweenthetwoapproaches.fortheendnotchedbeamapproach,no 34

distinctionbetweenmodesiandiiismadeintermsofthecrackpropagationcrite rionandthematerialischaracterizedbythemodeifractureenergyonly.another differencerelatestotheeffectofelasticclampingofthereducedbeampartsatthe sectionwherethenotchorholecornerislocated,whichisconsideredpartlydiffer entinthetwoapproaches.theexpressionforthenominalbeamstrengthwithre specttocrackingatarightanglednotchderivedbygustafsson(988)reads (6) where isthecrackshearforce, isthenetcrosssectionareaatthereduced crosssectionandwhereandaredefinedinfigure3.. Illustrationsofthepredictedbeamstrengthsasinfluencedbythenormalizedbeam heightandthenormalizednotchlengthaccordingtothetwoapproachesare showninfigure3..thecomparisonisbasedonabeamofheight=5mm,with = =3Nm/m andvaluesoftheexponentinequation5as==.for thischoice,i.e.for = and==,isthemixedmodecrackpropagationcrite rionthesameforbothcalculations.moreorlessslightdifferencesinpredicted strengthsaretobeexpectedduetothedifferencesinthemodelformulationswith respecttotheinfluenceoftheelasticclampingofthereducedcrosssection..5.5 αh.5h..4.6.8 Normalized beam height α h 4 Comparisontoexperimentaltests..4.6.8 Normalized notch length β Examplesoftheoreticallypredictedstrengthsforbeamswithaholearepresentedin Figures4.and4.,wherealsoresultsofexperimentaltestsonbeamswithcircular holes(höfflin5,aicherandhöfflin6)andwithsquareholes(danielsson8) areshown.theoreticalbeamstrengthsrefertothecrackpropagationloadpredicted accordingthetheorypresentedinsection,usingequation5with==.5and withmaterialdataaccordingtotable.thenetcrosssectionareaisdefinedas =( )andthestrengthfromexperimentaltestscorrespondstotheloadatthe instantofcrackpropagationacrosstheentirebeamwidth..5.5.5h βh h 343

Beamandholegeometriesarefortheseillustrationschosenbasedonavailabilityof testresults,yieldingpartlydifferentgeometriesforthecaseofcircularandsquare holesrespectively.resultspresentedinfigure4.relatetotheinfluenceofbeam height,holesizeandbendingmomenttoshearforceratioforholescentricallyplaced withrespecttobeamheight(=).resultspresentedinfigure4.relatetotheinflu enceofholecornerradiusandholeplacementwithrespecttobeamheight. 3.5.5.5 Test results a = hd = r =.3h M/(Vh) =.65..4.6.8 Beam height h [m] 3.5.5.5 Test results h = 45 mm M/(Vh) =.65 ±.5..4.6.8 Normalized hole height hd/h [ ] 3.5.5.5 Test results a = hd = h/3 r =.hd M/(Vh) =.7..4.6.8 Beam height h [m] 3.5.5.5 Test results a = hd = h/3 r =.hd h = 63 mm M(Vh) =.7..4.6.8 Normalized hole height hd/h [ ] V 3.5.5.5 Test results a = hd =.3h h = 45 mm 4 6 8 Moment to shear force ratio M/(Vh) [ ] 4 6 Moment to shear force ratio M/(Vh) [ ] 3.5.5.5 Test results a = hd = h/3 r =.hd h = 63 mm 344

5 Comparisonofstrengthanalysisapproaches Anoverviewoftheratiobetweenexperimentallyfoundbeamcapacityandtheoreti callypredictedcapacityaccordingtosomedifferentapproachesforstrengthanalysis isgiveninfigure5..theexperimentaltestsusedforthiscomparisonconsistof beamswitheithercircularholes(höfflin5,aicher&höfflin6)orsquareholes (Danielsson8)alsousedforcomparisoninSection4. Thestrengthanalysismethodsincludedare:. PFM AProbabilisticFractureMechanicsmethodbasedonDFEanalysisand considerationoffractureductilityaccordingtoageneralizedlefmapproachin combinationwithconsiderationofmaterialvariabilityaccordingtoweibulltheory (Danielsson&Gustafsson).. Weibull ClassicalWeibulltheorybasedonDFEanalysis(Danielsson&Gus tafsson). 3. NLFM Anonlinearfracturemechanicsapproach(acohesivezonemodel)based on3dfeanalysis(danielsson&gustafsson4). 4. Presentapproach seeprevioussections. 5. GlulamHandbook endnotchedbeamanalogyapproachforbeamswithahole foundintheoldversionoftheswedishglulamhandbook(carling)andalso includedinadraftversionofec5(). 6. DINEN995/NA SemiempiricalapproachfoundinGermanandAustrian NationalAnnexestoEC5. 7. Aicheretal DesignapproachbasedonWeibulltheorypresentedbyAicher,Höf flinandreinhardt(7).originallysuggestedtobeusedforcircularholesonly buthereusedalsoforsquareholesassuming= =. 3.5.5 Test results a = hd =.3h h = 45 mm.5 M/(Vh) =.65 s =...3.4.5 Normalized hole corner radius r/hd [ ] Shear strength Vc/Anet [Pa] 3.5.5 Test results, h = 8 mm Test results, h = 63 mm a = hd = h/3.5 r =.hd M/(Vh) =.7.... Hole placement w.r.t beam height s/h [ ] 345

Theoretical capacity / Experimental capacity [ ] 3..5..5..67.5.4.33 ( ) ( ) ( ) FE based methods ( ) ( ) code type methods ( ). PFM. Weibull 3. NLFM 4. Present approach 5. Glulam Handbook 6. DIN EN 995 /NA 7. Aicher et al Forthegeneralapproaches(items4above),experimentallyfoundmeanvaluesare comparedtopredictedstrengthbasedonassumedmeanvaluesofmaterialproperty parameters.forthenlfmapproach(item3above),onlyresultsforthebeamswith squareholesarepresentedsinceanalysisofbeamswithcircularholeswerenotin cludedintheworkpresentedin(danielssonandgustafsson4). Thestrengthaccordingtothethreecodetypeapproaches(items57above)are basedoncharacteristicmaterialstrengthvalues =3.5MPaand =.5MPa, validforstrengthclassgl3haccordingtossen48:3.thepredicted strengthsareforthesemethodscomparedtocharacteristicvaluesoftheexperimen tallyfoundstrengths;see(danielsson&gustafsson)forfurtherinformation. Alltestresultsincludedinthecomparisonrelatetoholeswithlargerheight than allowedtobeusedwithoutreinforcementaccordingtothegermanandaustrianna tionalannexestoec5. 6 Discussion Thepresentedapproachisgeneralinthesensethatitallowsforstrengthanalysisin cludingnotonlyloadingintermsofshearforceandbendingmoment,butalsoaxial force.dependingonrelativeloadlevelsandthesignoftheaxialforce,theassumed locationofthecrack(seefigure.)mayhoweverneedtobeadjusted. 346

Beamswithholesareinpracticaldesignsituationsoftenreinforced.Thereasonfor thisislikelyrelatedtotwodifferentfactors:(i)theactualstrengthreductiondueto theholeand(ii)theuncertaintyandlackofknowledgerelatedtostrengthanalysis anddesignofbeamswithahole.thecrosssectionsectionalforces, and (see Equations35)couldpossiblybeofinterestinrelationtodesignofinternalorexter nalbeamreinforcementintermsoffullythreadedscrews,gluedinrodsorgluedon panels. AnewversionoftheSwedishGlulamHandbookisplannedtobereleasedduring 5,withinwhichtheendnotchbeamanalogyapproachforbeamswithaholewill beremovedinfavourofthedesignapproachedfoundinthegermanandaustrian NationalAnnexestoEC5. 7 Conclusions Fromtheworkpresentedinthispaper,thefollowingconclusionsaredrawn: ThepresentapproachisconsistentwiththeLEFMbaseddesignapproachforend notchedbeamsfoundinec5. Thepresentapproachappearstobeabletocapturethestrongbeamsizeinflu encefoundfromexperimentaltestsofbeamswithaholefairlywell. Thepresentapproachappearsfairlyaccurateincapturingtheabsolutevaluesof thebeamstrengthforbothcircularandsquareholesofvarioussizesandlocations. Althoughbasedonbeamtheoryanalysis,thepresentapproachisnotsuitablefordi rectincorporationintotimberengineeringdesigncodesofpracticeinitscurrent form,becauseofrathercomplexequations.forcertainapplicationsandloadconfigu rations,moreuserfriendlydesignexpressionsmaypossiblybederivedandcould thenserveasabaseforimproveddesignrecommendations.furtherworkregarding verificationandcalibrationishoweverneededbeforethiscanberealized. 8 References AicherS,HöfflinL(4):Newdesignmodelforroundholesinglulambeams.In: Proceedingsof8 th WorldConferenceonTimberEngineering,vol,Lahti,Finland. AicherS,HöfflinL(6):TragfähigkeitundBemessungvonBrettschichtholzträgern mitrundendurchbrüchen SicherheitsrelevanteModifikationenderBemes sungsverfahrennacheurocode5unddin5.materialprüfungsanstalt(otto GrafInstitute),UniversityofStuttgart,Germany. AicherS,HöfflinL,ReinhardtHW(7):RundeDurchbrücheinBiegeträgernaus Brettschichtholz.Teil.TragfähigkeitundBemessung.Bautechnik84():86788. CarlingO():Limträhandbok(Glulamhandbook).SvensktLimträAB,Print&Me diacenterisundsvallab,sweden. 347

DanielssonH(8):Strengthtestsofglulambeamswithquadraticholes Testre port.reporttvsm753,structuralmechanics,lunduniversity,sweden DanielssonH,GustafssonPJ(8):Strengthofglulambeamswithholes Testof quadraticholesandliteraturetestresultcompilation.in:proceedingsofcibw8 Meeting4,StAndrews,Canada,Paperno.CIBW8/44. DanielssonH,GustafssonPJ():Aprobabilisticfracturemechanicsapproachand strengthanalysisofglulambeamswithholes.eurjwoodprod69:4749. DanielssonH,GustafssonPJ(4):Fractureanalysisofgluedlaminatedtimber beamswithaholeusinga3dcohesivezonemodel.engngfractmech4 5:895. DINEN995/NA:38(3):GermanNationalAnnextoEC5 EN48:3(3):Timberstructures Gluedlaminatedtimberandgluedtimber Requirements.CEN. Eurocode5(4):DesignoftimberstructuresPart:General Commonrules andrulesforbuildings.cen.(en995). Eurocode5():DesignoftimberstructuresPart:General Commonrules andrulesforbuildings.finaldraft9.(pren995). GustafssonPJ(988):Astudyofstrengthofnotchedbeams.In:ProceedingsofCIB W8Meeting,Parksville,Canada,Paperno.CIBW8/. GustafssonPJ():Meanstressapproachandinitialcrackapproach.In:AicherS, GustafssonPJ(ed)HallerP,PeterssonH:Fracturemechanicsmodelsforstrength analysisoftimberbeamswithaholeoranotch areportofrilemtc33.report TVSM734,StructuralMechanics,LundUniversity,Sweden. GustafssonPJ(5):Mixedmodeenergyreleaseratebyabeamtheoryappliedto timberbeamswithahole.in:abstractbookof th InternationalConferenceon Fracture,Turin,Italy. HellanK(985):Introductiontofracturemechanics.InternationalEdition.McGraw Hill. HöfflinL(5):RundeDurchbrücheinBrettschichtholzträger Experimentelleund theoretischeuntersuchungen.phdthesis,materialprüfungsanstalt(ottografin stitute),universityofstuttgart,germany. KolbH,EppleA(985):VerstärkungvondurchbrochenenBrettschichtholzbindern. VorschungsvorhabenI.4348,ForschungsundMaterialprüfungsanstaltBaden Württemberg,Germany. PeterssonH(974):Analysisofloadbearingwallsinmultistoreybuildings:stresses anddisplacementscalculatedbyacontinuummethod.phdthesis,chalmersuni versityoftechnology,sweden. ÖNORMB995:4(4):AustrianNationalAnnextoEC5. 348

Discussion The paper was presented by H Danielsson H Blass commented that the stress distribution shown in slide 8 close to the hole is very different from reality. He questioned why the results are so good. H Danielsson agreed but stated that the beam theory is exact for normal stresses but might be different for shear stresses. PJ Gustafsson added that the solution should be exact in terms of energy release rate. As the crack propagated, the normal forces due to the moment should be exact but shear would not be true. BJ Yeh asked how close the hole can be to the support and if there were two holes, how close can they be. H Danielsson responded the distances should be such that the support would not influence the hole. This would also apply to the multiple holes cases. K Malo asked if there would be a limitation to the hole size. H Danielsson said that there would be no limit and the analysis would also work for notched beams. I Smith commented that good results need accurate analysis, real structure, and accurate material properties; when results do not agree, maybe some of these are not working. P Dietsch and H Danielsson discussed the cases of round and square holes and limitations to the model. R Jockwer received clarifications that at this moment there is no recommendation for the size of the hole without reinforcement. 349