Rješenje čvora rešetke spajalom velikog promjera

Similar documents
Truss girder joint with a large diameter mechanical fastener

EXPERIMENTAL ANALYSIS OF COMBINED ACTION OF BENDING, SHEAR AND TORSION ON TIMBER BEAMS

Shear Modulus and Shear Strength Evaluation of Solid Wood by a Modified ISO Square-Plate Twist Method

RELIABILITY OF GLULAM BEAMS SUBJECTED TO BENDING POUZDANOST LIJEPLJENIH LAMELIRANIH NOSAČA NA SAVIJANJE

EXPERIMENTAL ANALYSIS OF THE STRENGTH OF A POLYMER PRODUCED FROM RECYCLED MATERIAL

COMPARISON OF MODULY OF ELASTICITY OBTAINED BY NON- DESTRUCTIVE AND DESTRUCTIVE TESTS OF TIMBER SAMPLES

PARAMETRIC ANALYSIS OF LARGE DIAMETER MECHANICAL FASTENER IN REINFORCED GLULAM

Mode I Critical Stress Intensity Factor of Medium- Density Fiberboard Obtained by Single-Edge- Notched Bending Test

ANALYTICAL AND NUMERICAL PREDICTION OF SPRINGBACK IN SHEET METAL BENDING

Algoritam za množenje ulančanih matrica. Alen Kosanović Prirodoslovno-matematički fakultet Matematički odsjek

TEORIJA SKUPOVA Zadaci

Projektovanje paralelnih algoritama II

BENDING-SHEAR INTERACTION OF LONGITUDINALLY STIFFENED GIRDERS

NAPREDNI FIZIČKI PRAKTIKUM 1 studij Matematika i fizika; smjer nastavnički MJERENJE MALIH OTPORA

MATHEMATICAL MODELING OF DIE LOAD IN THE PROCESS OF CROSS TUBE HYDROFORMING

DIJAGRAMI ZA ODABIR POPREČNOG PRESJEKA NOSAČA OD DRVA ZA RAZLIČITE PROTUPOŽARNE OTPORNOSTI

LINEAR AND NON-LINEAR ANALYSES OF CABLE-STAYED STEEL FRAME SUBJECTED TO SEISMIC ACTIONS

Red veze za benzen. Slika 1.

NAPREDNI FIZIČKI PRAKTIKUM II studij Geofizika MODUL ELASTIČNOSTI

Metode praćenja planova

Dijagram moment savijanja zakrivljenost za armiranobetonske grede

Mechanical analysis of timber connection using 3D finite element model

Čelični plošni elementi opterećeni u svojoj ravnini: faktori izbočivanja i kritična naprezanja

ON DERIVATING OF AN ELASTIC STABILITY MATRIX FOR A TRANSVERSELY CRACKED BEAM COLUMN BASED ON TAYLOR EXPANSION

PROBABILISTIC ASSESSMENT OF CALCULATION RESISTANCE MODELS OF COMPOSITE SECTION WITH PARTIAL SHEAR INTERACTION

STRESS OF ANGLE SECTION SUBJECTED TO TRANSVERSAL LOADING ACTING OUT OF THE SHEAR CENTER

ANALYSIS OF THE RELIABILITY OF THE "ALTERNATOR- ALTERNATOR BELT" SYSTEM

INVESTIGATION OF UPSETTING OF CYLINDER BY CONICAL DIES

Uvod u relacione baze podataka

DIPLOMSKI RAD Bojan Mihaljević

The Prediction of. Key words: LD converter, slopping, acoustic pressure, Fourier transformation, prediction, evaluation

STRESSES AND DEFORMABILITY OF ROCK MASS UPON OPEN PIT EXPLOITATION OF DIMENSION STONE. Siniša DUNDA, Petar HRŽENJAK and Trpimir KUJUNDŽIĆ

INNOVATION OF PARAMETER α cc FOR DESIGN RESISTANCE OF HIGH-STRENGTH CONCRETE COLUMNS

Mjerenje snage. Na kraju sata student treba biti u stanju: Spojevi za jednofazno izmjenično mjerenje snage. Ak. god. 2008/2009

FTN Novi Sad Katedra za motore i vozila. Drumska vozila Uputstvo za izradu vučnog proračuna motornog vozila. 1. Ulazni podaci IZVOR:

A NEW THREE-DIMENSIONAL CHAOTIC SYSTEM WITHOUT EQUILIBRIUM POINTS, ITS DYNAMICAL ANALYSES AND ELECTRONIC CIRCUIT APPLICATION

FIBER KONAČNI ELEMENT U NELINEARNOJ ANALIZI KVADRATNIH SPREGNUTIH CFT STUBOVA FIBER FINITE ELEMENT IN NONLINEAR ANALYSIS OF SQUARE CFT COLUMNS

DESIGN AND CALCULATION OF RING SPRINGS AS SPRING ELEMENTS OF THE WAGON BUFFER UDC : Jovan Nešović

PRIPADNOST RJEŠENJA KVADRATNE JEDNAČINE DANOM INTERVALU

Numerical model for determining fire behaviour of structures

Modelling of Experimental Tests of Wooden Specimens from Scots Pine (Pinus sylvestris) with the Help of Anisotropic Plasticity Material Model

Steel plate elements loaded in their plane buckling factors and critical stresses

STRUCTURAL ANALYSIS OF NORTH ADRIATIC FIXED OFFSHORE PLATFORM

Oracle Spatial Koordinatni sustavi, projekcije i transformacije. Dalibor Kušić, mag. ing. listopad 2010.

A COMPARATIVE EVALUATION OF SOME SOLUTION METHODS IN FREE VIBRATION ANALYSIS OF ELASTICALLY SUPPORTED BEAMS 5

Analiza ab elemenata u scenariju progresivnog urušavanja metodom konačnih elemenata

Mathcad sa algoritmima

SVEUČILIŠTE JOSIPA JURJA STROSSMAYERA U OSIJEKU GRAĐEVINSKI FAKULTET OSIJEK

ESTIMATION OF UNIAXIAL COMPRESSIVE AND TENSILE STRENGTH OF ROCK MATERIAL FROM GYPSUM DEPOSITS IN THE KNIN AREA

Discrete dome model for St. Jacob cathedral in Šibenik

Dynamic analysis of 2-D and 3-D quasi-brittle solids and structures by D/BEM

KLASIFIKACIJA NAIVNI BAJES. NIKOLA MILIKIĆ URL:

PROJEKT NOSIVE ARMIRANOBETONSKE KONSTRUKCIJE OBITELJSKE KUĆE

Analiza čelične i aluminijske konstrukcije na lokaciji Dugopolje-visina rešetke 2,36m

Theoretical and Numerical Approach in Determining the Thermal and Stress Loads in Train Disc Brakes

Application of Reduced Stiffness of Complex Laminate in Finite Elements for Chair Analysis

Effect of foundation in dynamic analysis of concrete gravity dams

Philippe Jodin. Original scientific paper UDC: :519.6 Paper received:

Yu.G. Matvienko. The paper was presented at the Twelfth Meeting New Trends in Fatigue and Fracture (NT2F12) Brasov, Romania, May, 2012

ANALIZA I MODELIRANJE SPREGNUTIH SISTEMA DRVO BETON: PRIMENA NA PRORAČUN MOSTOVSKE KONSTRUKCIJE PREMA EN

FIZIKALNA KOZMOLOGIJA VII. VRLO RANI SVEMIR & INFLACIJA

Optimizacija Niza Čerenkovljevih teleskopa (CTA) pomoću Monte Carlo simulacija

Static inelastic analysis of steel frames with flexible connections

VELOCITY PROFILES AT THE OUTLET OF THE DIFFERENT DESIGNED DIES FOR ALUMINIUM EXTRUSION

ZANIMLJIV NAČIN IZRAČUNAVANJA NEKIH GRANIČNIH VRIJEDNOSTI FUNKCIJA. Šefket Arslanagić, Sarajevo, BiH

QUARRY STABILITY ANALYSIS FOR COMPLEX SLIP SURFACES USING THE MATHSLOPE METHOD

USPOREDNI PRORAČUN ARMIRANOBETONSKE PLOČE PREMA GRANIČNOM STANJU NOSIVOSTI (GSN) I UPORABLJIVOSTI (GSU)

SVEUČILIŠTE JOSIPA JURJA STROSSMAYERA U OSIJEKU GRAĐEVINSKI FAKULTET OSIJEK DIPLOMSKI RAD

ANALYTICAL, NUMERICAL AND EXPERIMENTAL STRESS ASSESSMENT OF THE SPHERICAL TANK WITH LARGE VOLUME

Using the Energy Balance Method in Estimation of Overhead Transmission Line Aeolian Vibrations


ODREĐIVANJE DINAMIČKOG ODZIVA MEHANIČKOG SUSTAVA METODOM RUNGE-KUTTA

SOUND SOURCE INFLUENCE TO THE ROOM ACOUSTICS QUALITY MEASUREMENT

EXPERIMENTAL ANALYSIS OF AN ORIGINAL TYPE OF STEEL SPACE TRUSS NODE JOINT UDC : : (045)

STRUCTURAL VEHICLE IMPACT LOADING UDC =111. Dragoslav Stojić #, Stefan Conić

Termodinamika. FIZIKA PSS-GRAD 29. studenog Copyright 2015 John Wiley & Sons, Inc. All rights reserved.

COMPARISON OF LINEAR SEAKEEPING TOOLS FOR CONTAINERSHIPS USPOREDBA PROGRAMSKIH ALATA ZA LINEARNU ANALIZU POMORSTVENOSTI KONTEJNERSKIH BRODOVA

Kontrolni uređaji s vremenskom odgodom za rasvjetu i klimu

Fajl koji je korišćen može se naći na

EFFECT OF LAYER THICKNESS, DEPOSITION ANGLE, AND INFILL ON MAXIMUM FLEXURAL FORCE IN FDM-BUILT SPECIMENS

SVEUČILIŠTE U ZAGREBU FAKULTET STROJARSTVA I BRODOGRADNJE DIPLOMSKI RAD. Mateja Bičak. Zagreb, 2014.

Alternativna metoda za analizu izvijanja lameliranih kompozitnih greda

Eksperimentalno ispitivanje injektiranog prstenastog preklopnog spoja uslijed vlačnog opterećenja

Frost Formation Phenomenon in a Fin-and-Tube Heat Exchanger

STATISTICAL ANALYSIS OF WET AND DRY SPELLS IN CROATIA BY THE BINARY DARMA (1,1) MODEL

Slika 1. Slika 2. Da ne bismo stalno izbacivali elemente iz skupa, mi ćemo napraviti još jedan niz markirano, gde će

Geometrijski smisao rješenja sustava od tri linearne jednadžbe s tri nepoznanice

APPROPRIATENESS OF GENETIC ALGORITHM USE FOR DISASSEMBLY SEQUENCE OPTIMIZATION

Periodi i oblici titranja uobičajenih okvirnih AB građevina

Đorđe Đorđević, Dušan Petković, Darko Živković. University of Niš, The Faculty of Civil Engineering and Architecture, Serbia

REVIEW OF GAMMA FUNCTIONS IN ACCUMULATED FATIGUE DAMAGE ASSESSMENT OF SHIP STRUCTURES

Primjena tehnika mekog računalstva za predviđanje svojstva zrakopropusnosti tkanina

PARAMETRIC OPTIMIZATION OF EDM USING MULTI-RESPONSE SIGNAL-TO- NOISE RATIO TECHNIQUE

ANALYTICAL AND NUMERICAL STRESS ANALYSIS OF THE ROTARY KILN RING

CLINICAL. Neodoljiva ponuda iz Ivoclar Vivadenta PROLJEĆE LJETO. Ponuda traje od: ili do isteka zaliha

NUMERIČKA ANALIZA PROCESA ŠIRENJA PUKOTINA KONSTRUKCIJA

ATOMSKA APSORP SORPCIJSKA TROSKOP

Ground vibrations level characterization through the geological strength index (GSI)

Elastic - plastic analysis of crack on bimaterial interface

Transcription:

UDK 64.07.3.00.5:64.078.8 Građevinar 0/03 Primljen / Reeived: 0.6.03. Ispravljen / Correted: 5.0.03. Prihvaćen / Aepted: 5.0.03. Dostupno online / Available online: 0..03. Rješenje čvora rešetke spajalom velikog promjera Autori: Krunoslav Pavković Miljenko Haiman Mladen Meštrović Vlatka Rajčić Rješenje čvora rešetke spajalom velikog promjera Izvorni znanstveni rad Dr.s. Krunoslav Pavković dipl.ing.građ. Tehničko Veleučilište u Zagrebu Graditeljski odjel krunoslav.pavkovi@tvz.hr U radu je prikazan priključak spajalom velikog promjera za čvorove rešetkastih nosača velikih raspona. Eksperimentalnim istraživanjima i metodom konačnih elemenata provedena je analiza otpornosti i modula popustljivosti priključka. Eksperimentalna ispitivanja provedena su na četiri rešetkasta nosača s predloženim priključkom. U programskom paketu Abaqus/CAE s UMAT podrutinom prikazana je nelinearna analiza modela rešetke. Dobiveni rezultati prikazani u radu obrađeni su prema EN 5 a pokazuju da predloženi priključak ima veliku otpornost i modul popustljivosti primjenjiv za rešetke velikih raspona te da izraz dan u EC5 za karakteristični pritisak po omotaču rupe treba mijenjati. Ključne riječi: čvor rešetke nelinearna analiza spajalo velikog promjera ulijepljene šipke UMAT podrutina Original sientifi paper Krunoslav Pavković Miljenko Haiman Mladen Meštrović Vlatka Rajčić Truss girder joint with a large diameter mehanial fastener Prof.dr.s. Miljenko Haiman dipl.ing.građ. Sveučilište u Zagrebu Arhitektonski fakultet miljenko.haiman@arhitekt.hr A joint with a large diameter fastener for onneting large span truss girder nodes is presented in the paper. The joint strength and slip modulus were analysed by means of experimental testing and the finite element method. The experimental testing was onduted on four truss girders with the proposed joint. The nonlinear truss girder model analysis was onduted using the Abaqus/CAE software with the UMAT subroutine. The results obtained and analysed aording to EN 5 show that the proposed joint has a large strength and slip modulus and that it an be used for large-span truss girders. The results also show that the equation for embedding ompression strength aording to EC5 should be modified. Key words: truss girder joint non-linear analysis large diameter fastener glued-in rods UMAT subroutine Prof.dr.s. Mladen Meštrović dipl.ing.građ. Sveučilište u Zagrebu Građevinski fakultet mestar@grad.hr Prof.dr.s. Vlatka Rajčić dipl.ing.građ. Sveučilište u Zagrebu Građevinski fakultet vraji@grad.hr Wissenshaftliher Originalbeitrag Krunoslav Pavković Miljenko Haiman Mladen Meštrović Vlatka Rajčić Fahwerkträgerknoten mit Befestigungselementen großen Durhmessers Fahwerkträgerknoten mit Befestigungselementen großen Durhmessers Zusammenfassung In dieser Arbeit ist ein Knotenpunkt mit Befestigungselementen großen Durhmessers für Fahwerkträger bedeutender Spannweiten dargestellt. Festigkeit und Steifigkeit der Verbindung sind durh experimentelle Versuhe und numerishe FEM-Analysen ermittelt. Die Versuhe sind an vier Trägern mit der vorgeshlagenen Verbindung abgeshlossen während nihtlineare Analysen mit dem Programm Abaqus/CAE im UMAT Unterprogramm durhgeführt sind. Die erhaltenen Resultate zeigen dass der Knoten ausreihende Festigkeit und Steifigkeit vorweist und für Träger großer Spannweiten angewandt werden kann. Außerdem wird auf die notwendige Modifizierung der Gleihung für Einbettungsdruk nah EC5 hingewiesen. Shlüsselwörter: eingeklebte Stäbe Verbindung großen Durhmessers nihtlineare Analyse Fahwerkträgerknoten UMAT Unterpogramm GRAĐEVINAR 65 (03) 0 869-878 869

Građevinar 0/03 Krunoslav Pavković Miljenko Haiman Mladen Meštrović Vlatka Rajčić. Uvod Prilikom projektiranja konstrukija velikih raspona statički sustav rešetkastih nosača svojim malim utroškom materijala i jednostavnošću logičan je i ekonomski isplativiji u odnosu na punostijene lamelirane nosače. Međutim drveni rešetkasti nosači u našim uvjetima se vrlo rijetko koriste za svladavanje velikih raspona tj. kao alternativa se gotovo uvijek koriste lamelirani nosači. Glavni razlog za ovakvo stanje u praksi je problem rješavanja čvorova rešetke lobi drvne industrije te negativan stav arhitekata prema rešetkastim nosačima. Uobičajeno je i najčešće rješenje za spajanje drvenih rešetkastih nosača čeličnim čvornim limovima i štapastim spajalima koji zahtijevaju veliku površinu elemenata rešetke kako bi se sila adekvatno prenijela unutar čvora uz dostatnu duktilnost. Sve to često rezultira povećanjem elementa rešetke kako bi se potreban broj spajala mogao razmjestiti na dostatnoj površini što povećava potrošnju drva i čelika koja prelazi granie ekonomske isplativosti. Postoje i drugi razlozi zbog kojih se izbjegavaju rešetkasti nosači za svladavanje velikih raspona a to su: smanjena vatrootpornost u području spojeva i otežano zadavanje nadvišenja što je nužno kod raspona većih od 300 metara. Zbog uvjerenja da je jednostavnije utrošiti i više od 50 % više lameliranog drva za punostijeni nosač nego rješavati problem čvorova rešetke i vatrootpornosti nameće se pitanje može li se stanje promijeniti uvođenjem novih sustava spajanja. Cilj rada je pronaći rješenje priključka za rešetke velikog raspona koje će omogućavati brzu montažu zadavanje nadvišenja u montaži te povećanu vatrootpornost. Naime postoji jedno razvijeno rješenje za rešetke velikih raspona koje je još u fazi istraživanja SP-Trätek Borås Švedska a bazira se na spojnom sredstvu velikog promjera []. Autori ovog konepta bavili su se također mišlju da su dosadašnja rješenja kompliirana i preskupa te da je potrebno novo rješenje primjenjivo za rešetke velikih raspona... Konept priključka Na slii. prikazan je konept priključka koji je tema rada. Konept je osmišljen primjenom spojnog sredstva velikog promjera. Glavni element priključka je čelična ijev velikog promjera koji se određuje prema potrebnoj otpornosti priključka a ugrađuje se u rupu izbušenu u pojasu okomito na os elementa. Viji su elementi priključka koji povezuju ijev s elementom koji unosi silu u priključak a promjer i broj vijaka ovisi o veličini sile koju priključak treba preuzeti. Drveni element u koji se unosi sila (pojas rešetke) na slii. širine je 00 mm međutim njegova se širina određuje prema potrebnoj otpornosti pojasa rešetke. Pojas rešetke konipiran je tako da se njegova širina može jednostavno mijenjati brojem lamela ovisno o potrebnoj otpornosti a visina se određuje širinom lamela. Druga polovina priključka sastoji se od elementa koji unosi silu u spoj. Unos sile u spoj može biti riješen pomoću drvenog elementa s ulijepljenom šipkom [] (kako je prikazano na slii.) ili direktno iz čelične šipke. Šipka se na svome kraju neovisno o tome je li ulijepljena ili samostalna montira u čahuru za nastavljanje. Slika. Prikaz priključka rešetke.. Duktilnost modul popustljivosti i otpornost spojeva Duktilnost je karakteristika materijala da ima izražene velike deformaije u plastičnom području prije loma a može se reći da je to nužna karakteristika za konstrukije koje se nalaze u seizmički aktivnim područjima. Duktilnost (D) se izražava kao odnos između pomaka u točkama sile popuštanja (d y ) i čvrstoće (d max ) [3]: D = δ max δ y Drvene konstrukije prilikom potresnog opterećenja trebaju pružiti dostatnu razinu duktilnosti kako bi iz deformaija u plastičnom području proizašla potrebna disipaija energije statičkog sustava. CEN (European Committee for Standardisation) propisuje postupak za određivanje mehaničkih karakteristika spojeva u drvenim konstrukijama kroz sljedeće norme: -- za ispitivanje spojeva statičkom silom u normi EN 689 -- za ispitivanje spojeva iklički promjenjivom silom u normi EN 5. Međutim za određivanje duktilnosti spojeva opterećenih statičkom silom navedeni propisi ne daju posebne smjernie [4]. S obzirom na to da se ispituje spoj opterećen statičkom silom za određivanje duktilnosti primijeniti će se smjernie iz norme EN 5. Značajna karakteristika svakog spoja su: iniijalni modul popustljivosti spoja K u i modul popustljivosti spoja K s. Modul () 870 GRAĐEVINAR 65 (03) 0 869-878

Rješenje čvora rešetke spajalom velikog promjera Građevinar 0/03 popustljivosti može se odrediti iz nagiba prava određenog sa 0 % i 40 % maksimalne sile te odgovarajućeg pomaka na 0 % maksimalne (d 0 ) i 40 % maksimalne sile (d 04 ) prema izrazu: F F K = 04 max 0 S δ δ 04 0 max Iniijalni modul popustljivosti spoja određuje se iz nagiba prava određenog s početnim pomakom i pomakom pri sili od 40 % maksimalne sile F max (d 04 ). Razlika između ova dva modula određena je nepovratnim deformaijama uzrokovanim dosjedanjem spojnog sredstva. Budući da se u promatranom priključku radi o spajalu velikog promjera mjerodavan mod otkazivanja je uslijed dosezanja tlačne čvrstoće po omotaču rupe drvenog elementa pri čemu spajalo ostaje nedeformirano. Računska otpornost (R d ) za mehanizam otkazivanja drva na tlak po omotaču rupe dana je izrazom [5]: R d = f had t d (3) gdje su: f had - računska otpornost na tlak po omotaču rupe za silu pod kutom α u odnosu na vlakana t - debljina drvenog elementa za koji se određuje računska otpornost d - promjer spajala. Karakteristična otpornost po omotaču rupe spajala određuje se eksperimentalnim ispitivanjima prema normi EN 383 [6]. Prema spomenutoj normi karakteristična čvrstoća na tlak po omotaču rupe dobiva se utiskivanjem spajala određenog promjera u drveni element u kojem je prethodno napravljena rupa odgovarajućeg promjera. Utiskivanje spajala provodi se pod određenim kutom u odnosu na vlakana a karakteristična čvrstoća po omotaču rupe dobiva se iz sile (F max5 ) koja odgovara pomaku od 50 mm.. Priprema istraživanja.. Geometrija i korišteni materijali rešetkastih nosača za eksperimentalna ispitivanja Za potrebe eksperimentalnih ispitivanja proizvedena su četiri modela rešetkastih nosača s predloženim novim koneptom priključka osnog raspona 60 m i osne udaljenosti pojaseva 0 m. Rešetke su izvedene s tri polja odnosno sa četiri vertikalne tlačne ispune i vlačnim dijagonalnim ispunama u prvom i zadnjem polju (slika ). Elementi za rešetkaste nosače proizvedeni su od lamela debljine 40 mm. Svaki element sastojao se od pet lamela a nakon blanjanja imao je poprečni presjek 0 00 [m]. Za lijepljenje lamela upotrijebljeno je melaminsko ljepilo "Caso adhesives MUF system 47/56" proizvedeno u Caso adhesives AB Stokholm Švedska. Proizvedeni lamelirani () elementi prema mehaničkim karakteristikama primijenjenog drva za lamele mogu se svrstati u klasu GL4h prema EC5 normi. Kako bi se dobio uvid u ijeli konept priključka izvedene su ulijepljene šipke dubine uljepljivanja 350 m u vertikalnim elementima ispune rešetke. U sredini drvenog elementa izbušena je rupa 0 mm većega promjera od nazivnog promjera navojne šipke kako je prikazano u radu [7]. Navojne šipke odabrane su kvalitete 8.8. a za njihovo uljepljivanje upotrijebljeno je epoksidno univerzalno ljepilo proizvođača KGK d.o.o. iz Karlova. Ljepilo se vodi pod komerijalnim nazivom "Epoon '88". Kvaliteta primijenjenih materijala za čelične elemente odabrana je prema naprezanjima dobivenim preliminarnim numeričkim modelima. Za izradu podložnih profiliranih pločia čahura za nastavljanje i ijevi odabrana je kvaliteta čelika oznake 4CrMo4 f u =600 N/mm prema EN-007-. Za elemente spoja koji su manje izloženi naprezanjima (prsteni) upotrijebljen je čelik S35JRCu prema normi EN 007-. Za vijke u spojevima korišteni su inbus viji prema DIN-u 9 kvalitete.9. a za unos sile u spoj iz čeličnih dijagonala korištena su dva vijka M0 kvalitete.9. te za konstrukijsko povezivanje vertikale na pojas rešetke jedan vijak M kvalitete.9. Za vlačne dijagonale upotrijebljene su čelične šipke promjera 0 mm umjesto drvenih elemenata s ulijepljenim šipkama kako je prikazano na slii. Razlog zamjene drvenih dijagonala čeličnima jest jednostavnije mjerenje i kontrola naprezanja u čeličnim dijagonalama te manja proizvodna ijena rešetki za eksperimentalna ispitivanja. Vlačne dijagonale izvedene su od kvalitete čelika 4CrMo4 f u =600 N/mm... Postupak eksperimentalnih ispitivanja Oslanjanje nosača provedeno je pomoću čeličnih ploča dimenzija 00 000 [mm] na osnom razmaku 6000 m. Te su ploče bile položene na čelične profile sa zaobljenjem na vrhu te su na taj način simulirani zglobni osloni. Za vrijeme ispitivanja gornji je pojas rešetke bočno pridržavan konstrukijom sastavljenom od čeličnih ijevi (slika.) kako ne bi došlo do gubitka ukupne stabilnosti nosača. Proes ispitivanja sastojao se od nanošenja konentriranog opterećenja u trećinama raspona s konstantnim prirastom. Unos opterećenja provodio se uz pomoć dva hidraulična ilindra međusobno spojena tako da je pritisak u oba ilindra bio uvijek isti. Ispitivanje rešetkastih nosača provodilo se uz relaksaiju ukupne sile na 600 kn uslijed čega je sila vraćena na 00 kn a nakon toga kontinuirano povećavana sve do loma rešetki. Unesena sila mjerena je uređajima koji su postavljeni između hidrauličnih ilindara i gornjeg pojasa rešetke što je vidljivo na slii. Za vrijeme ispitivanja mjerena su naprezanja u čeličnim dijagonalama pomoću tenzometara s oznakama od T do T4. Mjerenje deformaija rešetke provodilo se u trećinama raspona mjernim uređajima s oznakama L500/03 i L500/04. GRAĐEVINAR 65 (03) 0 869-878 87

Građevinar 0/03 Krunoslav Pavković Miljenko Haiman Mladen Meštrović Vlatka Rajčić Spojevima koji se nalaze na donjem pojasu rešetke praćena je lokalna deformaija odnosno pomak vrha vijaka koji ulaze u dijagonalu u odnosu na točku koja se nalazi na osi vijka udaljena 000 mm od sredine ijevi. Ukupno četiri mjerna uređaja s oznakama od L00/0 do L00/05 postavljena su na dva promatrana spoja na donjem pojasu rešetke odnosno jedan mjerni uređaj po vijku kojim se unosi sila u spoj iz vlačne dijagonale. Izmjerena naprezanja i pomai praćeni su i zapisivani na osobnom računalu svake sekunde. Slika. Eksperimentalna ispitivanja rešetki Potrebno je spomenuti da su prilikom sastavljanja rešetke u laboratoriju svi viji u spojevima rešetke pritezani kontroliranom silom od 45 kn odnosno momentom pritezanja 70 Nm. To pritezanje vijaka rezultiralo je nadvišenjem u sredini raspona od 0 mm. Za vrijeme ispitivanja izmjerena je vlažnost elemenata rešetki u graniama od 8 % do 37 %. Temperatura u prostoriji za vrijeme ispitivanja druge rešetke bila je 30 C a za vrijeme ostalih ispitivanja 8 C..3. Analiza metodom konačnih elemenata.3.. Geometrija i materijali numeričkog modela rešetke Numerički model rešetkastog nosača napravljen je programskim paketom Abaqus /CAE ver. 6.0 a analiza je provedena sa materijalnom i geometrijskom nelinearnosti. Geometrija numeričkog modela preuzeta je u ijelosti od rešetki za laboratorijska ispitivanja. Za modeliranje drvenih elemenata korišteni su volumni C3D6 konačni elementi određeni preko šest čvorova [8]. Svi čelični elementi spoja te vlačne zatege modelirane su s volumnim C3D8 konačnim elementima određenim preko osam čvorova. Numerički model rešetke u ijelosti je određen sa 68.454 konačna elementa a s obzirom na velik broj konačnih elemenata modelirana je samo jedna polovina rešetke s rubnim uvjetima simetrije kako je prikazano na slii 3. Opterećenje je modelirano u dva koraka; u prvom koraku zadana je sila prednapinjanja vijaka i zatega kao što je provedeno i u eksperimentalnim ispitivanjima a u drugom koraku zadano je opterećenje na rešetku kao pritisak na čeličnu ploču iznad tlačne vertikale kako je prikazano na slii 3. Budući da se radi o nelinearnoj analizi s geometrijskom i materijalnom nelinearnosti opterećenja su zadana u ovisnosti o vremenu. Sila prednapinjanja Slika 3. Prikaz numeričkog modela s rubnim uvjetima i opterećenjima 87 GRAĐEVINAR 65 (03) 0 869-878

Rješenje čvora rešetke spajalom velikog promjera Građevinar 0/03 vijaka zadana je kao "bolt load" dok je prirast sile zadan u ovisnosti o vremenu. Linearni prirast sile prednapinjanja u svakom vijku modeliran je kroz vrijeme od 0 t do 0 t a za opterećenje na rešetku kroz vrijeme od 5 t do 465 t. Grafički prikazi amplituda sile u ovisnosti o vremenu za oba koraka opterećenja dani su na slii 3. Modeliranje mehaničkih karakteristika čelika provedeno je elastoplastično prema dijagramima prikazanim na slii 4. Plastično područje modelirano je pomoću 5 točaka koje su dobivene kao srednja vrijednost iz stvarnih dijagrama a moduli elastičnosti i granie popuštanja zadani su prema vrijednostima danim u tablii. Viji su modelirani bilinearno s graniom popuštanja i modulima elastičnosti prikazanim u tablii. a tangenijalni modul elastičnosti određen je kao 0 % vrijednosti modula elastičnosti. mnogo složeniji od modeliranja čelika te neizvediv u osnovnom programskom paketu Abaqus/CAE za potrebe određivanja karakteristika drva upotrijebljena je UMAT podrutina. Pomoću podrutine određena su dva kriterija popuštanja drva u zoni spoja kriterij Tsai-Wu a za ostalo drvo izvedenia iz Hillovog kriterija popuštanja. Koordinatni sustav kojim se opisuju mehaničke karakteristike drva označen je osima R T i L (redom radijalni transverzalni i longitudinalni smjer). Moduli elastičnosti i moduli posmika za svaki pojedini smjer dobiveni su kao prosječne vrijednosti drva slične zapreminske težine drva (ρ = 44347 kg/m 3 ) koje je upotrijebljeno za proizvodnju rešetke iz dostupne literature [9-4]. Vrijednosti iz navedenih izvora su obrađene i prikazane u tablii. Srednje čvrstoće drva za svaki pojedini smjer usvojene su iz literature [5 6]. Ispitivanja koja su opisana u odabranoj literaturi provedena su na istoj kvaliteti drva zapreminska težina drva je gotovo identična a vrsta i podneblje u kojem je drvo raslo su isti..3.. UMAT podrutina i kriteriji popuštanja drva Slika 4. Radni dijagrami čelika korištenih za elemente spoja Mehaničke karakteristike drva modelirane su elastoplastično i ortotropno. Budući da je ovaj postupak zbog ortotropnosti Podrutina UMAT služila je za određivanje kriterija popuštanja i određivanje mehaničkih karakteristika drva. Prije implementaije podrutine u programski paket ABAQUS kao "user-defined material" trebalo je UMAT podrutinu napisati u programskom jeziku FORTRAN [7]. Jedan od glavnih problema prilikom numeričkog modeliranja drva zbog njegove ortotropnosti i nehomogenosti bio je taj kako odrediti kriterij popuštanja. Zato je kao osnovni kriterij popuštanja u UMAT podrutini primijenjena izvedenia Hillovog kriterija za homogene materijale. Hillov kriterij popuštanja ortotropnog homogenog materijala [8] dan je izrazom: Tablia. Mehaničke karakteristike materijala Materijal Meh. karakteristike Modul elastičnosti [N/mm ] Modul posmika [N/mm ] Poissonov koefiijent Čvrstoća u tlaku [N/mm ] Čvrstoća u vlaku [N/mm ] E L =.340 G LR = 6770 ν LR = 04 f ul = -86 ft ul = 98 Lamelirano drvo GL 4h E T = 4690 G RT = 500 ν RT = 057 f ut = -330 N ft ut = 038 E R = 7460 G LT = 6580 ν LT = 05 f ur = -330 ft ur = 038 - - - f us = 33 f us = 33 Čelik 4CrMo4 E = 00.0540 G = 76.9440 ν = 03 f y =.450 f u =.6350 Čelik S 35JRCu E = 9.690 G = 84.960 ν = 03 f y = 5350 f u = 5958 Vijak kvalitete.9. (modeliran bilinearno) E = 0.0000 G = 84.0000 ν = 03 f y =.0800 - GRAĐEVINAR 65 (03) 0 869-878 873

Građevinar 0/03 Krunoslav Pavković Miljenko Haiman Mladen Meštrović Vlatka Rajčić ( ) ( ) ( ) = A σr σt B σt σl C σl σr D τrt E τtl F τ LR A = ; B = ; ( fur) ( fut ) ( ful ) ( fut ) ( ful ) ( f ) ur C = ; D = ; ( ful ) ( fur ) ( fut ) ( furt ) = ; F = ( futl ) ( fu LR E ) gdje je s normalno naprezanje t posmično naprezanje a f u čvrstoća u odgovarajućem smjeru. Kako gornji izraz ne uzima u obzir kriterij krtog loma drva u vlaku odnosno ako gornji uvjet nije zadovoljen mehaničke karakteristike za sve smjerove računaju se u plastičnom području. Iz navedenog razloga gornji uvjet ne opisuje dovoljno dobro ponašanje drva u vlaku i krti lom. Mnogo točnije ponašanje bit će ako se izraz rastavi te se za svaki smjer postavi zasebni kriterij [9 0]: σl τrt σr τlr σt = ; TL = ; τ = t t t f f f f f ul ush ur ush ut f ush gdje je f u t vlačna čvrstoća f ush posmična čvrstoća jedinstvena za sve smjerove. Ovim kriterijem za vlak se uzima u obzir interakija normalnog i posmičnog naprezanja za promatranu ravninu a njihovim odvajanjem osigurana je neovisnost tri međusobno ortogonalne ravnine. Kriterij popuštanja u tlaku određen je jednostavnim kriterijem maksimalnih naprezanja. Ako maksimalna naprezanja prelaze graniu popuštanja za jedan smjer tada se za sva tri smjera usvaja plastično područje: (4) (5) Drugi kriterij popuštanja drva u zoni spoja je Tsai-Wuov kriterij a određen je izrazom: FσL FσR F3σT FσLσR F3σLσT F3σTσR F σl Fσ R F3 3 σt 44τLR 55τLT 66τRT F F F = Zanemare li se vrijednosti F F 3 F 3 zbog malog utjeaja izraz (7) može se usvojiti u skraćenom obliku: F σ F σ F σ F σ L R 3 T L F σ R F σ 33 T F τ 44 LR F τ 55 LT F τ = 66 RT F = ; F = ; F3 = ; F = ; F t t = ; t f t f f t f fut fut ful ful fur fur ul ul ur ur F33 = ; F = F = F = t 44 55 66 fut fut f ush Kako bi se dobile zone u radnom dijagramu kriterij popuštanja prikazan izrazom (8) usvojen je na način da je prva grania popuštanja određena ako izraz dosegne vrijednost 00 druga kada je dosegnuta vrijednost 07 i treća kada je dosegnuta vrijednost 0. To je prikazano grafički na slii 5. (7) (8) σl σr σt = ; = ; = (6) f f f ul ur ut gdje je f u tlačna čvrstoća. Potrebno je istaknuti sličnosti između gore usvojenog kriterija popuštanja i Hasinova (980) kriterija popuštanja [] koji je namijenjen za materijale sa znatnom razlikom u mehaničkim karakteristikama u smjeru vlakanaa i okomito na njih. Slika 5. Tsai-Wuov kriterij popuštanja Slika 6. Radni dijagram drva - Hillov kriterij popuštanja Slika 7. Radni dijagram drva - Tsai-Wuov kriterij popuštanja 874 GRAĐEVINAR 65 (03) 0 869-878

Rješenje čvora rešetke spajalom velikog promjera Građevinar 0/03 Navedeni kriteriji popuštanja sa svojim uvjetima određeni su da bi se dobile zone u radnom dijagramu drva kako je prikazano na slii 6. i 7..3.3. Metoda konačnih elemenata u UMAT podrutini Nakon usvajanja kriterija popuštanja formulirana je metodom konačnih elemenata podrutina za elastični ortotropni materijal. Zbog složenosti podrutine koja bi opisivala elastoplastično ponašanje materijala primijenjena je formulaija konačnih elemenata za elastično ponašanje materijala. Modifikaija koja je napravljena na podrutini kako bi dala što točnije rezultate u elastoplastičnoj analizi jest u tome što se prema naprezanjima ako nije zadovoljen jedan od kriterija usvaja drugi (plastični) modul elastičnosti. Ovakvo pojednostavljenje je moguće u ovom slučaju jer je unos sile konstantan te se ne promatraju zaostale plastične deformaije. Elastična naprezanja prema metodi konačnih elemenata [] računaju se prema izrazu (9): 3. Analiza rezultata Radni dijagrami promatranih priključaka dobiveni su iz podataka koji su prikupljeni mjerenjem naprezanja u sve četiri dijagonale te mjerenjem pomaka vijaka koji su sastavni dio priključka. Sređivanjem rezultata dobivena je srednja unesena sila te srednji pomak spojnog sredstva za svaki priključak a prikazani su u grafikonu na slii 8. Iz radnih dijagrama koji su dobiveni za priključke rešetki može se prema metodi /6 odrediti sila popuštanja modul popustljivosti i duktilnost priključaka. Također s obzirom na to da se radi o priključima u kojima do popuštanja dolazi prije svega uslijed prekoračenja tlačne čvrstoće po omotaču rupe može se primjenom norme EN 383 [6] odrediti tlačna čvrstoća po omotaču rupe. {s} t = {s} t [M]{e} Dt (9) gdje su: {s} t - vektor naprezanja u vremenu t {s} t - vektor naprezanja u vremenu t [M] - ortotropna matria elastičnih konstanti {e} Dt - vektor deformaija u vremenu t D D D DRRRR D [ M]= D sim LLLL LLRR LLTT RRTT TTTT 0 0 0 0 0 0 0 0 0 DLRLR 0 0 D 0. LTLT {e} = {e L e R e T g RT g LT g LR } T {s} = {s L s R s T t RT t LT t LR } T D RTRT gdje je e normalna deformaija g posmična deformaija. Komponente matrie elastičnih konstanti računaju se prema: D LLLL = E L (-(n RT n TR ))g D RRRR = E R (-(n LT n TL ))g D TTTT = E T (-(n LR n RL ))g D LLRR = D RRLL = E L (n RL (n TL n RT ))g D LLTT = D TTLL = E L (n RL (n RL n TR ))g D RRTT = D TTRR = E R (n TR (n LR n TL ))g D LRLR = G LR D LTLT = G LT D RTRT = G RT γ = ν ν ν ν ν ν ν ν ν LR RL RT TR TL LT RL TR LT Slika 8. Prosječna vrijednost eksperimentalnih rezultata rešetke Da bi se izbjeglo krivo određivanje karakterističnih vrijednosti određena je srednja vrijednost svih rezultata priključaka rešetke što je prikazano na slii 8. (rna linija). Maksimalna sila prema normi EN 5 usvaja se kao manja vrijednost od tri uvjeta: maksimalne sile u trenutku potpunog otkazivanja priključka maksimalne sile do trenutka kada je zabilježen pad sile od 80 % i maksimalne sile pri deformaiji od 300 mm. U obradi dobivenih rezultata na srednju vrijednost dodana je pretpostavka ponašanja priključka u plastičnom području sve do deformaije od 300 mm. Naime analizom eksperimentalnih ispitivanja utvrđeno je da se s velikom sigurnošću može pretpostaviti da će priključak dosegnuti deformaiju od 300 mm prije nego što dođe do otkazivanja nekog njegovog elementa ili pada sile u radnom dijagramu. Također treba napomenuti da je na grafikonu prikazanom na slii 8. sila za sve priključke pa tako i za srednju vrijednost uvećana za 90 kn što je sila prednapinjanja ostvarena prije početka eksperimentalnih ispitivanja. Prema grafikonu prikazanom na slii 8. i određenoj maksimalnoj sili prema gore opisanom postupku određene su karakteristične vrijednosti. Obradom rezultata dobivena je sila popuštanja F y = 8 kn i maksimalna sila F max = 344 kn. Nadalje maksimalna sila pri pomaku od 50 mm GRAĐEVINAR 65 (03) 0 869-878 875

Građevinar 0/03 Krunoslav Pavković Miljenko Haiman Mladen Meštrović Vlatka Rajčić je F max5 = 8563 kn modul popustljivosti K S = 5489 kn/mm a duktilnost D = 808. prikazan je zajedno s radnim dijagramima eksperimentalnih ispitivanja na slii. (rvena linija). Naprezanja dobivena kod sile od 400 kn prikazana su za smjer paralelno s vlakanima S na slii 4. a okomito na vlakana S33 (tangenijalni smjer ako se napravi poveznia s modelom rešetke za eksperimentalna ispitivanja) na slii 5. Naprezanja su prikazana tako da je tumačenje oznaka ograničeno na vlačnu i tlačnu čvrstoću kako bi se lakše uočile zone u kojima naprezanja prelaze čvrstoću drva klase GL 4h. Slika 9. Deformaija srednjeg priključka druge rešetke Slika. Pomai u trećini rešetke dobiveni eksperimentalnim istraživanjem i numeričkim modelom Slika 0. Deformaija rubnog priključka druge rešetke U nastavku su dani zapisi zabilježeni u trećinama lokalno ojačanih rešetki praćeni uz pomoć mjernih uređaja L500/4 i L500/3. Deformaije dobivene numeričkim modelom rešetke prikazane su na slii. u trenutku pune vrijednosti sile prednapinjanja vijaka. Pomai dobiveni numeričkim modelom pokazuju maksimalni pomak od 50 mm na pojasevima rešetke uz rubne uvjete simetrije. Iz dobivenih deformaija razvidno je dobro podudaranje numeričkog modela s nadvišenjem koje je dobiveno prilikom pripreme rešetki za eksperimentalna ispitivanja. Provedenom analizom numeričkog modela rešetke simulirano je ponašanje rešetke do sile od 400 kn. Budući da se radi o vrlo velikim deformaijama naročito u zoni priključaka uslijed prevelikih deformaija konačnih elemenata prekinuta je konvergenija analize prije dosezanja zadane sile od 4500 kn. Pomai u trenutku sile od 400 kn prikazani su na slii 3. Cjelokupni radni dijagram deformaija u odnosu na silu Za naprezanja paralelno s vlakanima može se zaključiti da ne prelaze čvrstoću drva na većim površinama te kao takva nisu od većeg značenja za otpornost i krutost rešetke. Značajna naprezanja za otpornost priključaka i rešetke u jelini dobivena su u smjeru S33. Kao što je vidljivo na slii 3. u zonama oko priključaka pojavljuju se velike površine u kojima naprezanje prelazi čvrstoću okomito na vlakana. To naprezanje uzrokuje velike deformaije priključka odnosno utiskivanje čelične ijevi u pojas rešetke. Ista pojava utiskivanja ijevi je primijećena i prilikom eksperimentalnih ispitivanja. Najveća konentraija naprezanja većih od čvrstoće dobivena je u priključku iznad zglobnog ležaja. U istom priključku prilikom eksperimentalnih istraživanja dolazilo je do velikog utiskivanja čeličnog ležaja u pojas te razdvajanja lamela uslijed prevelikog tlaka okomito na vlakana. Na temelju eksperimentalnih rezultata izveden je izraz za karakterističnu čvrstoću paralelno s vlakanima po omotaču rupe a prikazana je na grafikonu (slika 6.) s još nekoliko izraza za čvrstoću po omotaču rupe u ovisnosti o promjeru spojnog sredstva. Iz grafikona se može vidjeti da izraz u EC5 ima prevelik pad naprezanja u odnosu na porast promjera spojnog sredstva a nešto bolje izraze navode S. Franke i P. Quenneville [3] te K. Sawata i M. Yasumura [4]. 876 GRAĐEVINAR 65 (03) 0 869-878

Rješenje čvora rešetke spajalom velikog promjera Građevinar 0/03 Slika. Pomai rešetke nakon uvođenja prednapinjanja Slika 3. Pomai rešetke pri sili od 400 kn Slika 4. Naprezanja S paralelno s vlakanima Slika 5. Naprezanja S33 okomito na vlakana 4. Zaključak Slika 6. Usporedba karakterističnih čvrstoća po omotaču rupe paralelno s vlakanima Iz dobivenih eksperimentalnih rezultata izveden je izraz za karakterističnu čvrstoću paralelno s vlakanima po omotaču rupe. Zaključuje se kako izraz dan u propisima EC5 treba promijeniti i nadopuniti kako bi se omogućilo dimenzioniranje spojeva sa spajalima promjera većega od 30 mm. Nadalje numerički modeli napravljeni programskim paketom Abaqus/CAE i UMAT podrutinom s navedenim Hillovim i Tsai-Wuovim kriterijem popuštanja dali su iznimno dobro podudaranje s eksperimentalnim rezultatima. Ovako složene numeričke analize omogućavaju provedbu parametarske analize na promatranom tipu priključka kojom bi se promatralo ponašanje spajala drugih promjera. GRAĐEVINAR 65 (03) 0 869-878 877

Građevinar 0/03 Krunoslav Pavković Miljenko Haiman Mladen Meštrović Vlatka Rajčić Konačno iz prvih rezultata prikazanih u radu može se zaključiti da je uz daljnja istraživanja na spajalima velikog promjera moguće s ovim tipom priključka projektirati rešetkaste nosače raspona od 300 do 500m. Priključak je prijavljen i zaštićen u Državnom zavodu za intelektualno vlasništvo Republike Hrvatske pod brojem P0049A. Zahvala Autori zahvaljuju Tehničkom veleučilištu u Zagrebu na finanijskoj pomoći bez koje ispitivanja prikazana u radu ne bi bila provedena. Također zahvaljuju na materijalnoj i uslužnoj donaiji tvrtkama: Hoja d.d. iz Slovenije te "Stolariji Grčić" iz Karlova. LITERATURA [] Kobel P.: Modelling of Strengthened Connetions for Large Span Truss Strutures in Department of Strutural Engineering 0 Lund Institute of Tehnology: Sweden. [] Rajčić V. Bjelanović A. Rak M.: Nosivost ulijepljenih navijenih čeličnih šipki u elementima od hrastovine. Građevinar 004. 56(3) p. 55-6. [3] J oahim H. Shadle P.: Dutility aspets of reinfored and nonreinfored timber joints. Engineering Strutures 0 33() p. 308-306. [4] Piazza M. Polastri A. Tomasi R.: Dutility of timber joints under stati and yli loads. Strutures and Buildings 0. 64(sb) p. 79-90. [5] CEN E.C.f.S. Euroode 5: Design of timber strutures. EN995- - Part -: General - Common rules and rules for buildings 004 Brussels: European ommittee for standardization. [6] CEN E.C.f.S. EN 383:007. Timber Strutures-Test methods- Determination of embedment strength and foundation values for dowel type fasteners 007. [7] Tlustohowiz G. Serrano E. Steiger R.: State-of-the-art review on timber onnetions with glued-in steel rods. Materials and Strutures 0. 44(5) p. 997-00. [8] Abaqus Analysis User s Manual in Volume IV: Elamants6.0. [9] Bodig J. Jayne A.B.: Mehanis of Wood and Wood Composites 98 New York: Van Nostrand Reinhold Company. p 7. [0] Carrington H.: The elasti onstants of sprue. Philosophial Magazine 93. 45 p. 055-057. [] Jenkin C.F.: Report on materials used in the onstrution of airraft and airraft engines 90 Her Majesty s Stationery Offie: London. p. 95-3. [] Doyle D.V. Drow J.T. MBurney R.S.: Elasti properties of wood. The Young s Moduli Moduli of Rigidity and Poisson s Ratios of Balsa and Quipo (58). The Elasti Properties of Wood. The Moduli of Rigidity of Sitka Sprue and heir Relations to Moisture Content (58-A). The Elasti roperties of Wood. Young s Moduli and Poisson s Ratio of Sitka Sprue and Their Relations to Moisture Content (58-B). Forest Produts Laboratory 945. 4 5 38. [3] CEN E.C.f.S. EN 94:999 Timber strutures - Glued laminated timber - Strength lasses and determination of harateristi values 999. [4] CEN E.C.f.S. EN 338:003 Strutural timber - Strength lasses 003. [5] Haiman M.: Analiza sigurnosti lameliranih nosača Građevinski fakultet 00 Sveučilište u Zagrebu: Zagreb. [6] Čizmar D.: Robusnost nosivih sustava drvenih rešetkastih nosača Građevinski fakultet 0 Sveučilište u Zagrebu: Zagreb. [7] Abaqus User Subroutines in Referene Manual6.0. [8] Hill R.: A Theory of the yielding and plasti flow of anisotropi metals. Proeedings the royal soiety A 984. 93 p. 8-97. [9] Xu B.H. et al.: Numerial and experimental analyses of multipledowel steel-to-timber joints in tension perpendiular to grain. Engineering Strutures 009. 3 p. 357-367. [0] Xu B.H. Bouhair A. Taazount M.: 3D non-linear finite element modelling of traditional timber onnetions. in World onferene on timber engineering 00. Riva de Garda. [] Hashin Z.: Failure Criteria for Unidiretional Fiber Composites. Journal of Applied Mehanis 980. 47() p. 6. [] Abaqus Analysis User s Manual in Volume 3: Materials6.0. [3] Franke S. Quenneville S.: Bolted and dowelled onnetions in Radiata pine and laminated veneer lumber using the European yield model. 0 Institution of Engineers Australia Australian Journal of Strutural Engineering. [4] Sawata K. Yasumura M.: Determination of embedding strength of wood for dowel-type fasteners The Japan Wood Researh Soiet 3 00. 48 p. 38-46. 878 GRAĐEVINAR 65 (03) 0 869-878