UNIVERZITET U BEOGRADU  MAŠINSKI FAKULTET                        Vlada M. Gašić    DINAMIČKA INTERAKCIJA   NOSEĆE STRUKTURE I KOLICA   PORTALNIH DIZALICA VISOKIH PERFORMANSI    Doktorska disertacija                                          Beograd, 2012.  UNIVERSITY OF BELGRADE  FACULTY OF MECHANICAL ENGINEERING                        Vlada M. Gašić    DYNAMIC INTERACTION BETWEEN THE  STRUCTURE AND THE TROLLEY OF HIGH  PERFORMANCE GANTRY CRANES     Doctoral Dissertation                                          Belgrade, 2012.        Komisija za ocenu i odbranu disertacije      Mentor:  dr Nenad Zrnić, vanredni profesor      Mašinski fakultet, Univerzitet u Beogradu      Komentor:  dr Srđan Bošnjak, redovni profesor      Mašinski fakultet, Univerzitet u Beogradu          dr Zoran Petković, redovni profesor      Mašinski fakultet, Univerzitet u Beogradu          dr Aleksandar Obradović, redovni profesor      Mašinski fakultet, Univerzitet u Beogradu          dr Milosav Georgijević, redovni profesor      Fakultet tehničkih nauka, Univerzitet u Novom Sadu              Datum odbrane:           god.                        DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA  PORTALNIH DIZALICA VISOKIH PERFORMANSI      REZIME    U  radu  se  analizira dinamičko ponašanje noseće  konstrukcije portalne dizalice usled  dejstva kolica kao pokretnog opterećenja. Dat  je prikaz klasifikacije portalnih dizalica  pri  čemu  su  izdvojene  portalne  dizalice  za  kontejnerske  terminale  sa  svojim  visokih  performansama koje imaju stalnu tendenciju poboljšanja. Prvo je dat koncept primene  analitičkog pristupa  za modeliranje noseće konstrukcije preko  sistema elastičnih  tela  tipa prizmatične  grede  i  razmatranje  slobodnih poprečnih oscilacija. Kao  savremen  i  pre  svega  neophodan,  usvojen  je  kombinovani  pristup  za  istraživanje  naslovnog  problema,  tj.  konačnoelementni  pristup  je  iskorišćen  za  modeliranje  noseće  konstrukcije  portalne  dizalice  a  principi  analitičke  mehanike  su  iskorišćeni  za  modeliranje kolica. Razmatraju se dva najčešća konstrukciona tipa portalne dizalice za  formiranje modela strukture. Kolica su obuhvaćena kroz model pokretne mase, model  pokretnog  oscilatora  i  kroz model  pokretnog  oscilatora  sa  klatnom  koji  predstavlja  originalan model pokretnog opterećenja.  Za  svaki od modela  je utvrđena dinamička  interakcija  između ovih sistema  i postavljeni su matematički modeli koji predstavljaju  sistem diferencijalnih  jednačina drugog  reda  sa promenljivim koeficijentima. Rešenja  su  dobijena  pomoću  originalnih  programa,  na  bazi  metode  direktne  integracije‐ Njumarkove metode.  Identifikacija  i analize odziva  su  izvršene  za dva  realna primera  portalnih  dizalica.  Istražen  je  uticaj  brzine,  ubrzanja/usporenja  i  težine  kolica,  kao  i  uticaj  klaćenja  tereta  i elastične opruge u  sistemu  kolica. Dobijeni  rezultati  se mogu  iskoristiti u početnim fazama konstruisanja portalnih dizalica koje  imaju tendenciju da  ostvare  veoma  visoke  performanse,  u  smislu  ostvarivanja  boljeg  uvida  u  dinamičko  ponašanje.      Ključne reči:    Portalna dizalica, dinamički odziv, pokretno opterećenje, MKE,        pokretna masa, pokretni oscilator, klatno, direktna integracija    Naučna oblast:   Tehničke nauke, mašinstvo    Uža naučna oblast:   Mehanizacija    UDK:      621.874.5:531/534 (043.3)          DYNAMIC INTERACTION BETWEEN THE STRUCTURE AND THE  TROLLEY OF HIGH PERFORMANCE GANTRY CRANES       ABSTRACT    The dynamics of a two‐dimensional gantry crane structure subjected to various types  of moving  load  is  examined  in  this work.  First,  the  classification  of  gantry  cranes  is  suggested and group of gantry cranes at container terminals are distinguished because  of high performances which have tendency to become even better in near future. The  analytical  approach  is  introduced  in  modeling  the  gantry  structure  as  continuous  system with  transverse  vibrations. However, modern  approach,  i.e.  combined  finite  element  and  analytical method  is  adopted  to  solve  the  title  problem.  Two  types  of  structure of gantry cranes are considered. Three types of trolleys are  implemented  in  calculation,  i.e. moving mass, moving  oscillator  and moving  oscillator with  swinging  payload as original model are considered as moving loads acting upon the structure of  the gantry cranes. The interaction between the structure and each moving load model  is  derived  and  the  governing  equations  for  MDOF  systems  are  obtained.  The  postulated  equations,  which  are  second  order  differential  equations  with  time  dependent coefficients, are solved with direct  integration method‐Newmark method.  The  analysis  is  applied  to  two  types  of  gantry  cranes  and  dynamic  responses  are  obtained  for both  the  structure and  the  trolley. There are  studied  factors of moving  loads  such  as  magnitude,  speed,  acceleration,  deceleration  and  factors  within  the  trolley structure such as swinging of the payload and spring stiffness. Numerical results  reveal  that  used  approach  is  useful  and  can  draw  conclusions  for  structural  design  purposes of gantry cranes.      Keywords:   Gantry cranes, Dynamic responses, Moving load, FEA, Moving  mass, Moving oscillator, Swinging payload, Direct integration    Scientific field:   Technical sciences, Mechanical engineering    Scientific discipline:   Dynamics of material handling and conveying machines    i SADRŽAJ 1. UVOD ...................................................................................................................................... 1 1.1 PORTALNE DIZALICE ............................................................................................................... 1 1.2 KLASIFIKACIJA PORTALNIH DIZALICA ........................................................................... 4 1.3 PERFORMANSE PORTALNIH DIZALICA .......................................................................... 9 2. PREDMET I PLAN ISTRAŽIVANJA ............................................................................... 11 2.1 PREDMET PRETHODNIH ISTRAŽIVANJA IZ OBLASTI POKRETNOG OPTEREĆENJA .................................................................................................................................. 11 2.1.1 Problem pokretne sile ............................................................................................... 11 2.1.2 Problem pokretne mase ........................................................................................... 13 2.1.3 Problem pokretnog oscilatora ............................................................................... 15 2.1.4 Problem pokretnog klatna ....................................................................................... 17 2.2 METODOLOGIJA ISTRAŽIVANJA ....................................................................................... 19 2.3 ZAKLJUČNA RAZMATRANJA .............................................................................................. 22 2.4 PLAN ISTRAŽIVANJA .............................................................................................................. 27 2.5 CILJ ISTRAŽIVANJA ................................................................................................................. 28 3. MODELIRANJE NOSEĆIH KONSTRUKCIJA PORTALNIH DIZALICA SISTEMOM ELASTIČNIH TELA ........................................................................................................... 30 4. FORMIRANJE MODELA PORTALNE DIZALICE PRIMENOM KOMBINOVANOG PRISTUPA .......................................................................................................................... 37 4.1 MODEL NOSEĆE KONSTRUKCIJE PORTALNE DIZALICE ........................................ 37 4.1.1 Konačnoelementni model strukture ................................................................... 39 4.1.2 Tip usvojenog konačnog elementa ....................................................................... 43 4.1.3 Matrica krutosti KE modela strukture ............................................................... 47 4.1.4 Matrica inercije KE modela strukture ................................................................ 48 ii 4.1.5 Slobodne neprigušene oscilacije strukture ...................................................... 48 4.1.6 Matrica prigušenja KE modela strukture .......................................................... 50 4.2 MODELI KOLICA ....................................................................................................................... 51 4.2.1 Pokretna masa ...................................................................................................... 53 4.2.2 Pokretni oscilator ............................................................................................... 53 4.2.3 Pokretni oscilator sa klatnom ........................................................................ 54 4.3 MODEL KRETANJA KOLICA ................................................................................................. 54 5. DINAMIČKA INTERAKCIJA NOSEĆE KONSTRUKCIJE I KOLICA......................... 57 5.1 KONAČNOELEMENTNO MODELIRANJE POKRETNOG OPTEREĆENJA STRUKTURE ........................................................................................................................... 61 5.1.1 Prezentacija modeliranja spoljašnjeg opterećenja prema modelu pokretne sile ..................................................................................................................... 66 5.2 MATRIČNA PREZENTACIJA POMERANJA I UBRZANJA TAČKE KONTAKTA STRUKTURE I KOLICA ....................................................................................................... 70 5.3 ODREĐIVANJE INTERAKTIVNIH SILA ZA POSTAVLJENE MODELE KOLICA ..................................................................................................................................... 74 5.3.1 Model pokretne mase ........................................................................................ 76 5.3.2 Model pokretnog oscilatora ........................................................................... 77 5.3.3 Model pokretnog oscilatora sa klatnom.................................................... 78 6. MATEMATIČKI MODELI PORTALNE DIZALICE ...................................................... 81 6.1 MATEMATIČKI MODEL PORTALNE DIZALICE PREMA MODELU POKRETNE MASE ................................................................................................................ 81 6.2 MATEMATIČKI MODEL PORTALNE DIZALICE PREMA MODELU POKRETNOG OSCILATORA .............................................................................................. 83 6.3 MATEMATIČKI MODEL PORTALNE DIZALICE PREMA MODELU POKRETNOG OSCILATORA SA KLATNOM ................................................................ 86 6.4 POSTUPAK REŠAVANJA DIFERENCIJALNIH JEDNAČINA .................................. 90 iii 7. IDENTIFIKACIJA I ANALIZA ODZIVA MODELA ...................................................... 91 7.1 ODZIV MODELA POKRETNE MASE .............................................................................. 93 7.1.A Analiza odziva ‐ model A .......................................................................................... 93 7.1.A.1 Uticaj brzine i ubrzanja kolica ................................................................ 97 7.1.A.2 Uticaj prigušenja u strukturi ................................................................ 101 7.1.A.3 Uticaj inercijalnih efekata pokretne mase ...................................... 103 7.1.B Analiza odziva ‐ model B ....................................................................................... 105 7.1.B.1 Uticaj brzine i ubrzanja kolica ............................................................. 107 7.1.B.2 Uticaj prigušenja u strukturi ................................................................ 111 7.2 ODZIV MODELA POKRETNOG OSCILATORA ............................................................ 111 7.2.A Analiza odziva ‐ model A ....................................................................................... 112 7.2.A.1 Uticaj krutosti opruge ............................................................................. 113 7.2.B Analiza odziva ‐ model B ....................................................................................... 114 7.2.B.1 Uticaj krutosti opruge ............................................................................. 115 7.3 ODZIV MODELA POKRETNOG OSCILATORA SA KLATNOM .............................. 116 7.3.A Analiza odziva ‐ model A ....................................................................................... 116 7.3.A.1 Odziv sistema usled kretanja kolica konst. brzinom ................. 116 7.3.A.2 Uticaj brzine i ubrzanja kolica ............................................................. 118 7.3.A.3 Uticaj dužine užetnog sistema ............................................................. 121 7.3.A.4 Uticaj krutosti opruge ............................................................................. 122 7.3.B Analiza odziva ‐ model B ....................................................................................... 123 7.3.B.1 Odziv sistema usled kretanja kolica konst. brzinom ................. 123 7.3.B.2 Uticaj brzine i ubrzanja kolica ............................................................. 125 7.3.B.3 Uticaj dužine užetnog sistema ............................................................. 128 7.3.B.4 Uticaj krutosti opruge ............................................................................. 130 8. ZAKLJUČAK ...................................................................................................................... 131 9. LITERATURA ................................................................................................................... 135 DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 1 1 UVOD Portalna dizalica predstavlja jednu od osnovnih vrsta mašina prekidnog transporta. U domaćoj literaturi iz oblasti dizalica još uvek se može naći i naziv ramna dizalica za ovu mašinu [1], pošto njena bazična konfiguracija ima oblik П‐ rama. Osnovna namena portalne dizalice jeste obavljanje transporta tereta na otvorenom prostoru. Na svetskoj mapi proizvođača portalnih dizalica Srbiju predstavlja preduzeće Goša‐Fom, Smederevska Palanka. 1.1. PORTALNE DIZALICE Glavni delovi portalne dizalice su kolica, noseća konstrukcija i pogon kretanja dizalice. Pod pojmom kolica podrazumeva se sistem koji se kreće duž glavnog nosača dizalice, sa mehanizmom za dizanje tereta kao osnovnom komponentom. Noseća konstrukcija portalne dizalice je čelična struktura koja se može sagledati kroz tri konstrukcione celine: glavni nosač, kruta noga i elastična (pendel) noga. Danas, projektanti izbegavaju varijantu portalne dizalice sa istom krutom i elastičnom nogom čak i kod malih raspona glavnog nosača (do 25 m). Osnovni razlog je taj što elastična noga svojom elastičnošću ili zglobom omogućuje da se ceo sistem prilagodi temperaturnim dilatacijama ili nepravilnostima dizalične šine. Na slici 1.1 prikazan je najopštiji konstrukcioni tip portalne dizalice sa glavnim merama koje definišu namenu same dizalice. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 2 Slika 1.1. Skica portalne dizalice sa glavnim merama Glavne konstrukcione mere portalne dizalice (slika 1.1) su:  Raspon između nogu (L)  Visina dizanja od kote terena (H)  Dohvat (prepust) sa strane terena (L1)  Dohvat (prepust) sa strane vode ili terena (L2)  Maksimalna visina dizanja (Hmax) Raspon između nogu (L) i dohvati sa strane terena (L1, L2) definišu širinu, a visina dizanja od kote terena (H) definiše visinu manipulativnog prostora za transport tereta za portalne dizalice opšte namene. U slučajevima upotrebe dizalica na rečnim lukama ili skladištima rasutih materijala, uslovi pretovara definišu dohvat sa strane vode (L2) i maksimalnu visinu dizanja (Hmax). Pretovar rasutih materijala se danas uglavnom vrši pretovarnim mostovima koji predstavljaju posebnu vrstu transportnih mašina, iako su određeni tipovi slični portalnim dizalicama. Takođe, pretovar komadnih materijala u lukama se obavlja kontejnerskim dizalicama koje zbog velikog dohvata (L2 > L) i velike visine dizanja (Hmax > L) predstavljaju posebnu vrstu transportnih mašina. Generalno, portalna dizalica se može koristiti svuda gde se transport materijala odvija na otvorenom prostoru. Ipak, mogu se izdvojiti dva objekta gde se transport materijala ne može zamisliti bez portalne dizalice: brodogradilišta i kontejnerski terminali u lukama. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 3 Portalne dizalice za brodogradilišta su ekstremno velike, pa se često nazivaju i Golijat (Goliath) dizalice. Prva ovakva dizalica je napravljena od strane proizvođača Krupp (Nemačka) i instalirana u brodogradilištu Harland i Wolff (Belfast, V.Britanija) 1969. godine. Ubrzo je dobila ime Golijat, što je kasnije usvojeno kao naziv za posebnu vrstu portalnih dizalica. U istom brodogradilištu je 1974. godine postavljena druga dizalica koja je dobila ime Samson. Obe dizalice su nosivosti po 840 t i imaju rasponi deo dužine 140 m, dok je Golijat visok 96 m, a Samson 106 m. Značaj portalnih dizalice za kontejnerske terminale je u korelaciji sa velikim značajem kontejnerskog transporta u svetskoj ekonomiji. Uz srednje godišnje povećanje kontejnerskog transporta robe od 8%, očekuje se da ukupan broj pretovarenih kontejnera u svetu 2015. godine iznosi preko 700 miliona kontejnera (TEU), [2]. Kako portalne dizalice predstavljaju nezamenljiv sistem u kontejnerskom transportu i koji je uvek prilagođen zahtevima kupca, realno je očekivati dalji razvoj ovih dizalica kod najvećih svetskih proizvođača (Konecranes, Liebherr, Kuenz, Gottwald...). a) b) Slika 1.2. Portalna dizalica za a) brodogradilište, b) kontejnerski terminal Motiv za izradu ove disertacije jeste razvoj postojećih modela za analizu dinamičkog ponašanja portalnih dizalica što predstavlja osnovu za projektovanje lakih a pouzdanih konstrukcija. Takođe, ovaj rad ima motiv da prikaže savremen pristup problemima u dinamici struktura. U užem smislu, ova disertacija je orjentisana na analizu dinamičkog ponašanja portalnih dizalica za kontejnerske terminale jer se razvoj i poboljšanje performansi očekuje kod ove vrste dizalica. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 4 1.2. KLASIFIKACIJA PORTALNIH DIZALICA Na osnovu podataka najvećih svetskih proizvođača moguće je dati predlog klasifikacije portalnih dizalica koji je prikazan na slici 1.3. Slika 1.3. Klasifikacija portalnih dizalica Proizvođači portalnih dizalica za kontejnerske terminale u lukama, uglavnom, prvu podelu portalnih dizalica vrše prema načinu kretanja portalne dizalice na terenu, i to na:  RMG dizalice (Rail Mounted Gantry)  RTG mobilne dizalice (Rubber Tired Gantry) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 5 RMG dizalice (slika 1.4.a) predstavljaju klasičnu varijantu portalne dizalice koja podrazumeva da se cela dizalica kreće po šinama i samim tim je predviđena za transport tereta na unapred određenom lokalitetu. Takođe, one obezbeđuju najveći kapacitet skladištenja i veoma su pogodne za korišćenje u terminalima velikog kapaciteta, [3]. RTG dizalice, slika 1.4.b, su mobilne dizalice gde se kretanje ostvaruje preko pneumatika. Ove dizalice se u evropskim zemljama koriste manje od šinskih dizalica. Veoma su pogodne za terminale gde je potrebno obezbediti veću mobilnost u radu jer omogućuju i poprečno kretanje u prostoru skladišta zaokretanjem točkova do ± 900. One su po pravilu manje od šinskih dizalica na istom objektu. a) b) Slika 1.4. Portalna dizalica a) RMG, b) RTG Prema tipu konstrukcije, uobičajeno je da se portalne dizalice dele na:  Osnovni tip (П ram)  Sa jednim ili dva prepusta Osnovni tip portalne dizalice ima oblik П‐rama. Ovaj tip konstrukcije je jedini oblik Golijat dizalica, tj. ekstremno velikih dizalica za brodogradilišta, slika 1.2.a., kao i RTG portalnih dizalica za kontejnerske terminale u lukama, slika 1.4.b. Varijanta konstrukcije sa prepustima, sa 1 prepustom (slika 1.2.b) ili sa 2 prepusta DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 6 (slika 1.5.b), predstavlja najbolje rešenje sa aspekta odvajanja skladišnog prostora i pretovarnog prostora. Dalje, po načinu izvođenja konstrukcije dizalice se dele na:  Konstrukcije sa punim nosačima  Rešetkaste  Kombinovane konstrukcija Konstrukcija sa punim nosačima je najčešće izvođenje kod portalnih dizalica. Elementi konstrukcije koja se izrađuje od punih nosača su uglavnom kutijaste konstrukcije, mada je negde primetno i izvođenje konstrukcija od cevastih profila i to najčešće za elemente elastične noge. Kompletna rešetkasta konstrukcija (slika 1.5.a) je najređa varijanta, a nešto češća varijanta je kombinovana konstrukcija gde je glavni nosač izveden kao rešetka a noge dizalice od punih nosača, slika 1.5.b, zato što dužine glavnog nosača kod ovih dizalica mogu doprineti uštedi u masi celokupne konstrukcije izvođenjem na ovakav način. a) b) Slika 1.5. Portalna dizalica a) rešetkasta, b) kombinovane konstrukcije, Gottwald Prema broju grednih nosača portalne dizalice se dele na jednogrede i dvogrede, što je slično kao i za mosne dizalice. Jednogrede portalne dizalice mogu biti kombinovane konstrukcije dok su dvogrede portalne dizalice po pravilu izrađene od punih nosača. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 7 Noseće konstrukcije portalnih dizalica su danas fleksibilne (savitljive). Ovo se odnosi na konstrukcije sa punim nosačima kod kojih je smanjenje mase konstrukcije jedan od bitnih faktora za izbor i ocenu performansi. Prema stepenu automatizacije rada portalne dizalice se dele na:  Automatizovane  Poluautomatizovane  Neautomatizovane Namena portalne dizalice diktira stepen automatizacije rada. Automatizovane i poluautomatizovane portalne dizalice predstavljaju neophodnost za kontejnerske terminale velikog kapaciteta, u modernim svetskim lukama. Kontejneri se nakon istovara sa broda skladište u luci u pravougaonim blokovima koje opslužuje portalna dizalica. Pozicioniranje kontejnera u bloku se odvija potpuno automatski, pri čemu se izbor lokacije i upravljanje vrše sistemom koji upravlja procesima skladištenja kontejnera (TOS). Radnje koje se vrše u toku ovog procesa se upravljaju sistemima koji, pomoću informacija od senzora, kamera, lasera i mernih uređaja na kolicima i hvatačima kontejnera, na kraju vrše pozicioniranje tereta sa tolerancijom od ± 5 centimetara (podaci ABB Automation). Razvijeni su softveri koji omogućavaju virtuelnu simulaciju transportnog toka u kontejnerskim terminalima, slika 1.6. Slika 1.6. Prikaz virtualne simulacije transporta kontejnera, TBA DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 8 Pojam poluatomatizacije podrazumeva da se sve radnje na skladištu odvijaju automatski, a da se tek u procesu pretovara kontejnera (do trenutka kada se kontejner nađe do 1 metra iznad vozila) uključuje i rukovaoc dizalice. Neautomatizovane portalne dizalice su sve ostale dizalice gde komandovanje radnjama obavlja rukovaoc iz kabine ili sa poda. Prema vrsti kolica, portalne dizalice je moguće podeliti na dizalice sa kolicima sa užetnim pogonom (RTT) i kolicima sa sopstvenim‐elektro pogonom. Kolica sa sopstvenim pogonom dizanja i kretanja su danas dominantna varijanta, a posebno kod portalnih dizalica koje opslužuju skladišni deo kontejnerskih terminala. Slika. 1.7. Kolica sa elektropogonom mehanizma kretanja i mehanizma dizanja Slika. 1.8. Obrtna platforma u sistemu kolica sa hvatačem kontejnera DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 9 Na slici 1.7 prikazan je klasičan tip ovih kolica koja imaju dva vitla (Dual hoist trolley). Ovaj koncept podizanja tereta je veoma važan jer se na taj način, preko glavnih i pomoćnih užadi i sistema koturova koji se nalaze na konstrukciji hvatača kontejnera omogućuje fino pozicioniranje tereta bez pomeranja kolica ili dizalice. Preko ovih mehaničkih sistema je moguće, u određenoj meri, uticati na sprečavanje njihanja tereta. Opcioni mehanizam kod ovih kolica je obrtni mehanizam kojim se omogućuje obrtanje kontejnera u horizontalnoj ravni (1..2 obrt./min), slika 1.8. 1.3. PERFORMANSE PORTALNIH DIZALICA Pojam performanse dizalice predstavlja tehničke podatke koji određuju njenu namenu. Osnovni podaci su nosivost i glavne mere dizalice (raspon između nogu, visina dizanja, dohvat). Drugi tip podataka koji se obavezno propisuje uz dizalicu su brzina dizanja, brzina kretanja kolica i brzina kretanja dizalice. Ovi parametri su važni za definisanje trajanja ciklusa operacija koji se izvode pomoću dizalice. Sa projektantskog aspekta, ove performanse definišu (ne)potrebnost dinamičke analize pored obavezne kvazistatičke analize noseće konstrukcije portalne dizalice. U tabeli 1.1 prikazane su perfomanse RTG dizalica, a u tabeli 1.2 performanse RMG dizalica, renomiranih evropskih proizvođača portalnih dizalica za kontejnerske terminale. Tabela 1.1. Performanse RTG portalnih dizalica za kontejnerske terminale Konecranes Liebherr Nosivost [t] do 40 do 50 Raspon između nogu [m] do 32 20,8...23,6 Visina dizanja [m] do 21 12,3...21 Brzina dizanja tereta [m/min] 45 28 Brzina kretanja kolica [m/min] 90 70 Brzina kretanja dizalice [m/min] 70 70 DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 10 Tabela 1.2. Performanse RMG portalnih dizalica za kontejnerske terminale Konecranes Liebherr Kuenz Nosivost [t] do 50,8 do 50 t do 46 Raspon između nogu [m] 19...50 22...70 do 45,7 Dohvat [m] do 12 do 20 do 20,9 Visina dizanja [m] 12...18 9,2...28 13,3...15,4 Brzina dizanja tereta [m/min] 30 40 25 Brzina kretanja kolica [m/min] do 150 do 180 do 150 Brzina kretanja dizalice [m/min] do 240 do 240 do 140 Ove vrednosti predstavljaju aktuelni presek stanja nivoa tehničkih rešenja razmatrane vrste dizalica. Brzina kretanja kolica kod RMG portalnih dizalica za kontejnerske terminale (proizvođača Liebherr) dostiže vrednosti od 180 m/min (3 m/s) čime se ova vrednost izjednačava sa brzinama kolica kod obalskih kontejnerskih dizalica Post Panamax klase [3]. Realno je očekivati da brzina kretanja kolica dostigne i veće vrednosti zbog stalne težnje za skraćivanjem ciklusa pretovara kontejnera na terminalima, a povećanjem raspona (dužina kretanja) kod glavnih nosača portalnih dizalica za kontejnerske terminale. Očekivanja se mogu poistovetiti sa tendencijom rasta brzina kretanja kolica kod obalskih kontejnerskih dizalica, pa treba očekivati brzine do 250 m/min, [4], sa vremenima ubzanja i usporenja koja iznose 5 s. U pregledu aktuelnog stanja (2010. god) ponude proizvođača Kocks (www.kockskrane.de) može se izdvojiti podatak (dizalica Boxer 6000) da je već ostvarena brzina kretanja kolica od 240 m/min, sa ubrzanjem pogona kretanja kolica od 0,83 m/s2. Predmet analiza u ovoj disertaciji je ocena mogućnosti povećanja performansi kretanja kolica i posledično njihovim uticajima na konstrukciju portalnih dizalica za kontejnerske terminale. Matematički modeli koji će biti postavljeni u ovom radu se mogu iskoristiti i za analizu dinamičkog ponašanja drugih vrsta dizalica gde brzina kretanja kolica ima tendenciju da dostigne gore pomenute vrednosti. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 11 2 PREDMET I PLAN ISTRAŽIVANJA Predmet istraživanja u ovom radu, u najužem smislu, pripada problematici pokretnog opterećenja. Prikaz relevantnih istraživanja iz ove oblasti je dat u nastavku. 2.1. PREGLED PRETHODNIH ISTRAŽIVANJA IZ OBLASTI POKRETNOG OPTEREĆENJA Problem pokretnog opterećenja predstavlja posebnu oblast dinamike pod kojim se podrazumeva razmatranje dinamičkih odziva elastične strukture usled dejstva opterećenja koje se kreće po strukturi. U literaturi su se izdvojile sledeće osnovne vrste problema pokretnog opterećenja: problem pokretne sile, problem pokretne mase i problem pokretnog oscilatora. Razmatranje oscilacija elastične strukture usled dejstva pokretnog sistema koji u sebi sadrži teret kome je omogućeno klaćenje je predmet savremenih istraživanja problema pokretnog opterećenja kod dizalica. Ovakva konfiguracija celog sistema će ovde biti izdvojena i nazvana problemom pokretnog klatna. 2.1.1. Problem pokretne sile Prinudne transverzalne oscilacije proste grede usled dejstva pokretne sile prvi su razmatrali Krylov [5] i Timoshenko [6], sa željom da daju bolji uvid u dinamičko ponašanje mostova železničkog i putničkog saobraćaja. Uvedene su sledeće pretpostavke: sila je konstantnog intenziteta, sila se kreće po nosaču DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 12 konstantnom brzinom, nosač je Bernuli‐Ojler greda uniformnog poprečnog preseka, a njene oscilacije su u domenu malih oscilacija bez prigušenja. Ovakva postavka problema se u literaturi naziva klasičan problem pokretne sile, slika 5.1, koji je kasnije detaljno opisao Inglis [7]. Važnost ovog koncepta se ogleda u činjenici da je tada analitički dokazano da brzina kretanja sile po gredi ima izuzetan uticaj na poprečne oscilacije grede. Uveden je pojam kritične brzine pod kojom je dobijen maksimalni dinamički ugib grede koji je 1,57 puta veći od najvećeg statičkog ugiba (na sredini grede). Dalje, Timoshenko je ovaj problem poboljšao uvođenjem prinudne pokretne sile čiji se intenzitet menja po harmonijskom (prostom) zakonu [8]. Uticaj kretanja kompozicije vagona po železničkom mostu je razmatran preko modela pokretnog kontinualnog opterećenja koji su postavili Goldenblat i Bolotin, [9,10]. Zastoj u izučavanju problema pokretnog opterećenja je nastao 60‐ih godina prošlog veka jer su analitički (matematički) pristupi dostigli svoju granicu rešivosti. Većina analitičkih formulacija pokretne sile je data u monografiji iz 1972. godine koja se smatra bazičnom knjigom za istraživače problema pokretnog opterećenja [11]. Autor, Ladislav Fryba, u 3. izdanju ove knjige u referencama dodaje radove koji uključuju metodu konačnih elemenata (MKE) u problematici pokretnog opterećenja. Upravo je razvoj nove numeričke metode, MKE, ponovo otvorio vrata ovom problemu. Olsson je u radu [12] verifikovao primenu MKE za rešavanje klasičnog problema pokretne sile. Kako su dobijeni rezultati bili u odličnoj saglasnosti sa analitičkim rešenjem, autor daje prednost novoj metodi i ukazuje na mogućnost korišćenja novih modela pokretnog sistema. Ovo je već bilo pokazano u njegovoj disertaciji gde je uveo nove modele struktura poput višeraspone grede i ramovske strukture po kojima se kreću pokretni sistemi sa elastično ovešenim masama, [13]. Praktično, u njegovim radovima je ukazano na neophodnost primene MKE u savremenim problemima pokretnog opterećenja gde je objekat složena struktura, tj. gde model strukture ne može da se sagleda pomoću modela proste grede. Wu i ostali [14], su slikovito opisali primenu konačnoelementne formulacije pokretne sile koja se konstantnom brzinom kreće po strukturi portalne dizalice. Ovaj rad je važan sa aspekta modeliranja pokretne sile u komercijalnim programima za DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 13 analizu struktura na bazi MKE. Primena ovog koncepta na dinamičko ponašanje mosnih dizalica u softveru SAP 2000 je data u [15]. U savremenim istraživanjima problem pokretne sile egzistira uglavnom kao verifikaciona podloga za dinamičke odzive struktura usled dejstva složenijih modela pokretnog sistema kao što su pokretna masa i pokretni oscilator. I pored toga, ovaj problem ima najveći praktični značaj za većinu inženjera koji se bave projektovanjem struktura izloženih dejstvu pokretnog opterećenja jer je lako primenljiv, za razliku od pomenutih složenijih modela. U tom smislu, još uvek je aktuelno nalaženje približnih, a dovoljno tačnih formula za rešenje klasičnog problema pokretne sile [16]. 2.1.2. Problem pokretne mase Problem pokretne mase koja se konstantnom brzinom kreće po prostoj gredi zanemarljive mase prvi je postavio i rešio Stokes [17]. Ovaj rad je bio od izuzetnog značaja jer je uključivao inercijalne efekte pokretne mase u proračun poprečnih oscilacija grede. Međutim, ovo je ipak aproksimacija generalnog problema pokretne mase koji uključuje inercijalne efekte mase koja se kreće po strukturi čija masa ne može biti zanemarena. Zbog složenosti ovog problema, prva istraživanja su nastala tek u drugoj polovini prošlog veka pri čemu se razmatraju samo jednomaseni modeli pokretnog sistema na prostim modelima struktura. U radovima se često iznosi doprinos ovoj problematici koji je dao čovek našeg porekla Milomir Stanišić. On je u radu [18] rešenje predstavio razvojem sopstvenih funkcija u redove pri čemu one zadovoljavaju konturne uslove proste grede, a koeficijenti su promenljivi u vremenu. Problem pokretne mase na konzolnom nosaču i ostalim konturnim uslovima jednoraspone grede je razmatran od strane Akina i Mofida kombinovanim analitičkim i numeričkim pristupom, [19]. Isti autori su u [20] prešli na rešavanje ovog problema pomoću diskretizovanog modela grede. Devedesetih godina prošlog veka ovaj problem je ponovo dobio na aktuelnosti. Klasičan problem pokretne mase (slika 5.2), koji uključuje razmatranje oscilacija grede usled inercijalnih efekata mase koja se kreće konstantnom brzinom po prostoj Bernuli‐Ojler gredi, je postavljen metodom pretpostavljenih DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 14 metoda pri čemu je rešenje dobijeno metodom Runge‐Kuta (Runge‐Kutta), [21]. Istovremeno, Esmailzadeh i Gorashi su ovaj problem prikazali u [22], pri čemu je rešenje dobijeno metodom konačnih razlika. Lee je uz rešavanje klasičnog problema pokretne mase najavio i razmatranje mogućnosti odvajanja pokretne mase i strukture prilikom kretanja mase veoma velikom brzinom, [23]. Tada, u literaturi je počeo da se usvaja pojam dinamičke interakcije kod problema pokretne mase sa kojim se podrazumevalo da prinudna sila koja deluje na strukturu pored težine mase uključuje, između ostalog, centripetalnu silu i silu od Koriolisovog ubrzanja pokretne mase u diferencijalnu jednačinu oscilacija struktura. Ove sile nastaju kao efekat kretanja mase po deformisanoj strukturi (krivolinijskoj putanji). Efekat ovih sila na oscilacije mostova je prikazao Michaltsos u [24] kojim je poboljšan klasičan problem pokretne mase koji je isti autor dao u [25]. Veći doprinos Michaltsos je dao u radu [26] gde je istraživao uticaj ubrzanja i usporenja u profilu kretanja pokretne mase na prostoj gredi. Ovo istraživanje je važno jer je ukazano da promenljiva brzina kretanja ima izuzetan uticaj u periodu ubrzanja tereta i samim tim predstavlja neophodnost u modeliranju kretanja pokretnog opterećenja. U domaćoj literaturi iz oblasti mašinstva, prvi je problem pokretne mase otvorio Zrnić [27] u svojoj doktorskoj disertaciji. Uticaj kretanja kolica na dinamičko ponašanje strele obalske kontejnerske dizalice je definisan metodom pretpostavljenih modova pri čemu su usvojene odgovarajuće dopustive funkcije prema metodi Rejli‐Rica, [28,29]. Dinamički koeficijenti za maksimalni ugib i moment savijanja na konkretnom primeru mega kontejnerske dizalice su dati u [30]. Na osnovu dobijenih rezultata zaključuje se da je problem pokretne sile dovoljno tačan za praktičnu upotrebu ali da problem pokretne mase predstavlja savremen naučni pristup u problematici pokretnog opterećenja kod dizalica i samim tim mora biti usvojen kao polazni koncept u dinamičkoj analizi. Wu, Whittaker i Cartmel [31] su predstavili koncept kombinacije MKE i analitičkog pristupa za dobijanje dinamičkih odziva strukture portalne dizalice usled dejstva specijalne konstrukcije kolica sa teretom. Prvo je izvršena konačnoelementna formulacija i implementacija interaktivnih sila koje deluju na strukturu pomoću interpolacionih funkcija konačnih elemenata koji su usvojeni za DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 15 diskretizaciju strukture. Clough i Penzien su polaznu ideju za ovaj koncept prikazali u svojoj knjizi [32], koja se u anglosaksonskoj literaturi smatra veoma značajnom knjigom iz oblasti dinamike struktura. Dalje, pomoću jednačina analitičke mehanike, dobijene su interaktivne sile između strukture i specijalnog tipa kolica. Ovaj koncept predstavlja pravi iskorak u naučnom istraživanju problema pokretne mase jer se na taj način omogućuje analiza složenih struktura izloženih dejstvu pokretnog opterećenja. Za razliku od čistog analitičkog pristupa gde je aproksimacija strukture na proste modele obavezna, ovako je moguće implementirati problem pokretne mase na strukture u celoj svojoj specifičnosti. Radovi ovih autora uvek sadrže verifikaciju postavljenih matematičkih modela preko rezultata dobijenih pomoću problema pokretne sile na prostijim modelima struktura, a negde sadrže i eksperimentalnu verifikaciju rezultata dobijenih na umanjenom laboratorijskom modelu. Važnost eksperimentalne verifikacije dinamičkih odziva dobijenih pomoću MKE je prikazan u [33]. Naime, na umanjenom modelu strukture portalne dizalice izvršeno je merenje frekvencija oscilovanja i upoređivanje sa vrednostima dobijenih na KE modelu u programu I‐Deas. Iako se vrednosti razlikuju do 30 %, što se može opravdati sa odstupanjima tačnosti opreme za merenje vibracija i sa izvođenjem spojeva na umenjenom modelu (gde oni postaju veoma kruti), autor naglašava da je potrebno vršiti verifikaciju rezultata dobijenih KE modeliranjem struktura sa eksperimentalnim merenjima gde god je to moguće. Problem pokretne mase na objektu stubne konzolne dizalice je prikazan u [34]. Rezultati pokazuju da je simulacija realnog ciklusa kretanja kolica u dinamičkoj analizi dizalica potreban parametar za upoređenje rezultata sa statičkom analizom. 2.1.3. Problem pokretnog oscilatora Problem oscilacija tereta vezanog pomoću opruge za osovinu točka koji se kotrlja bez klizanja po krutoj i neravnoj podlozi je uveden u istraživanja sa pojavom amortizera u konstrukcijama vozila, početkom 20. veka. Savremeni radovi razmatraju opštiji slučaj koji pretpostavlja da se kretanje odvija po DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 16 elastičnoj strukturi. U tom kontekstu, ovakva postavka se u prevodu sa engleskog može nazvati problem pokretnog oscilatora. Složenost ovog problema zavisi prevashodno od modela pokretnog sistema. Uglavnom su u upotrebi jednoosovinski i dvoosovinski modeli pokretnih sistema, [13], a u sklopu njih prosti (jedna opruga) ili dvostruki (dve opruge) oscilatori. Među prvima je ovu problematiku razmatrao Lin. U radu [35] je data konačnoelementna formulacija interakcije pokretnog opterećenja i jednoraspone grede pri čemu su razmatrani modeli jednoosovinskog i dvoosovinskog prostog oscilatora, ponaosob. Postavljene su jednačine kretanja i dobijeni su dinamički odzivi numeričkom metodom Runge‐Kuta. Ovde je na modelu proste grede opterećene pokretnom silom konstantnog intenziteta pokazano da je konačnoelementna formulacija pokretnog opterećenja u odličnoj korelaciji sa analitičkim rešenjem. Slika 2.1 Pokretni oscilator na obostrano ukleštenoj gredi [36] Nastavak istraživanja ovih autora je dat u [36] gde je razmatran uticaj diskretizacije jednoraspone grede na tačnost rešenja pri konačnoelementnoj formulaciji problema, slika 2.1. Zaključeno je da broj elemenata prilikom diskretizacije mora biti minimum 2 (a poželjno 8) puta veći od broja elemenata koji bi se koristili u statičkoj analizi grede pomoću MKE. Primena ovog koncepta na istraživanje dinamičkog ponašanja mosnih dizalica usled dejstva pokretnog jednoosovinskog dvostrukog oscilatora (kojim je modeliran užetni sistem i opruga u sistemu kolica) je prikazan u [37]. Dalji razvoj ovog problema jeste dinamička analiza mogućnosti odvajanja pokretnog oscilatora i strukture [38]. Zaključeno je da odvajanje ovih dvaju DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 17 sistema može nastati samo pri velikim brzinama kretanja. Takođe, ukoliko je krutost opruge u sistemu pokretnog oscilatora velika, model se može aproksimirati modelom pokretne mase ukoliko je kontakt strukture i pokretnog sistema stalan [39]. Potrebno je pomenuti još jedan parametar koji utiče na dinamičko ponašanje kod problema pokretnog oscilatora, a to je oblik gazeće površine strukture po kojoj se kreće pokretni sistem. Razlikuju se glatka (idealno ravna) i neravna (sa izbočinama ili udubljenjima) površina. Sa akcentom na primeni kod mostova u putničkom saobraćaju, zaključeno je u [40] da neravna površina zahteva razmatranje problema pokretnog oscilatora. 2.1.4. Problem pokretnog klatna Kao što je poznato, matematičko klatno predstavlja osnovni problem analitičke mehanike. Ukoliko se pretpostavi da se tačka vešanja klatna kreće po krutoj podlozi dobija se složeniji dinamički problem, koji može biti iskorišćen za modeliranje podsistema dizalica sa aspekta automatizacije i optimizacije procesa istovara rasutih materijala [41]. Razmatranje ovakve fizičke postavke klaćenja tereta još uvek predstavlja osnovu za formiranje algoritama u sistemu upravljanja kretanja kolica za sprečavanje klaćenja tereta (anti‐sway control). Razmatranje oscilacija elastične strukture usled dejstva pokretnog sistema koji u sebi sadrži teret kome je omogućeno klaćenje ima veliku primenu u radovima iz problematike pokretnog opterećenja kod dizalica. Prvi ovakav model, gde je objekat istraživanja konstrukcija jednogrede mosne dizalice, je predstavio Oguamanam postavljanjem diferencijalnih jednačina oscilacija proste grede pomoću Hamiltonovog principa, [42]. Pretpostavljeno je da se klaćenje tereta vrši samo u vertikalnoj ravni što je kasnije prošireno na razmatranje oscilacija grede u vertikalnoj i horizontalnoj ravni, [43], pri čemu se i klaćenje tereta odvija prostorno, slika 2.2. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 18 Slika 2.2 Model mosne dizalice sa teretom [43] Wu je u radu [44] uveo koncept pokretne matrice mase radi implementacije inercijalnih efekata pokretne mase na dinamičko ponašanje strukture portalne dizalice. Pritom, razmatrani su trapezni profili brzina kretanja pokretnog sistema. Rešenja su dobijena numeričkom metodom direktne integracije. Takođe, razmatran je uticaj klaćenja tereta u ravni upravnoj na osu glavnog nosača sa pretpostavkom poznate promene ugla klaćenja u vremenu po harmonijskom zakonu. Koncept pokretne matrice masa podrazumeva da se u ukupnoj matrici inercije sistema implementira uticaj pokretne mase u zavisnosti od položaja same mase na strukturi. Na taj način ukupna matrica inercije postaje matrica čiji se elementi menjaju u vremenu što predstavlja veoma složen problem sa aspekta određivanja rešenja postavljene diferencijalne jednačine koja opisuje problem. Ovde se daje prednosti metodi Njumarka (Newmark) za direktnu integraciju. Wu je 2008. godine proširio koncept pokretne matrice masa na koncept pokretnog konačnog elementa čime se na adekvatan način simulira dinamička interakcija strukture i pokretne mase čime matrica inercije, matrica prigušenja i matrica krutosti sistema postaju matrice sa promenljivim elementima, [45]. Na objektu strukture portalne dizalice istraživane su podužne i poprečne oscilacije strukture usled dejstva pokretnog sistema koji se idealizuje jednomasenim modelom. Uticaj klaćenja tereta u vertikalnoj ravni je aproksimiran svođenjem na DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 19 ekvivalentnu pokretnu masu sa pretpostavkom poznavanja funkcije ugla klaćenja tereta. Metoda Njumarka je iskorišćena za rešavanje postavljenog matematičkog modela. Koncept pokretnog konačnog elementa, kao inovativan i savremen pristup, je poslužio da se ponovo razmotri problem oscilacija proste grede usled dejstva pokretne mase koja se kreće promenljivom brzinom, [46]. Simulacija ciklusa kretanja kontejnera u softveru ADAMS, sa mogućnošću klaćenja, i uticaj na životni vek konstrukcije portalne dizalice je prikazana u [47]. Yazid je razmatrao oscilacije noseće konstrukcije portalne dizalice usled klaćenja tereta i elastičnosti užetnog sistema. Jednačine kretanja celog sistema su izvedene pomoću metode konačnih elemenata i Langranževih jednačina, [48]. Rešenja su dobijena numerički, kombinacijom metode Njumarka i metode Runge‐ Kuta. Zaključeno je da elastičnost strukture utiče na dinamičko ponašanje klatna i obrnuto. U radu [49], Younesian je postavio uticaj prostornog klaćenja tereta na opterećenje glavnog nosača, pri čemu je uključeno i kretanje sistema kolica sa trapeznim profilom brzina kretanja. 2.2 METODOLOGIJA ISTRAŽIVANJA Polaznu osnovu za rešavanje problema pokretnog opterećenja čine dva osnovna poglavlja teorije oscilacija: 1. Oscilacije elastičnih tela‐poprečne oscilacije prizmatične grede 2. Oscilacije sistema sa konačnim brojem stepeni slobode (MDOF system) Gore navedena problematika u domaćoj literaturi se može naći u knjigama [50,51,52], a u stranoj literaturi u [53]. Jednačine oscilovanja elastičnog tela tipa prizmatične grede se opisuju parcijalnim diferencijalnim jednačinama. Određivanje apsolutno tačnih sopstvenih frekvencija je moguće samo za proste slučajeve oslanjanja grede pa se često pribegava približnim metodama za određivanje frekvencija od kojih se mogu izdvojiti metod Rejlija i metod Rejli‐Rica. Određivanje osnovne frekvencije sistema je veoma važno ali ne i dovoljno sa aspekta oscilacija nekonzervativnih sistema, tj. sistema izloženih dejstvu prinudnih sila. Kod ovih problema, rešavanje DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 20 parcijalne diferencijalne jednačine je izuzetno složeno (nekad i nemoguće) sa matematičkog aspekta. Ovo je posebno primetno kod rešavanja problema pokretnog opterećenja gde je prinudna sila ne samo funkcija od vremena već i funkcija položaja. Uobičajeno je da se prinudna pokretna sila u jednačinama definiše preko specijalne funkcije ‐ Dirakove funkcije. Ovaj matematički pristup u radovima iz oblasti problema pokretnog opterećenja se klasifikuje kao analitički pristup. Dominantna metoda koja se izdvojila za postavljanje jednačina oscilacija usled dejstva pokretnog opterećenja jeste metoda pretpostavljenih modova, čime se sistem prevodi na konačan broj običnih diferencijalnih jednačina. Numeričke metode koje se koriste za rešavanje diferencijalnih jednačina su metoda Runge‐Kuta i metoda Njumarka. Na osnovu svega, potrebno je zaključiti da određeni nivo aproksimacije mora da postoji i u analitičkom pristupu, a rešavanje jednačina prati velike matematičke teškoće. Svođenje bilo kog sistema na konačan broj stepeni slobode predstavlja idealizaciju realnog sistema. Međutim, ovaj pristup je dominantan u teoriji oscilacija pa i dinamici konstrukcija kao njenoj posebnoj disciplini. Jednačine oscilovanja sistema sa konačnim brojem stepeni slobode se opisuju običnim diferencijalnim jednačinama drugog reda. Postavljanje jednačina se vrši pomoću Njutnovih zakona i osnovnih jednačina analitičke mehanike: Dalamberov princip i Langranževe jednačine druge vrste. Razvoj ovog pristupa je nastao sa razvojem računara i posledično sa razvojem posebne matematičke discipline koja se naziva matrični račun. Ovim su oscilatorni sistemi mogli biti razmatrani svođenjem na veliki broj stepeni slobode sa garancijom da dobijena rešenja daju zadovoljavajuću sliku o ponašanju osnovnog modela. Pored idealizacije sistema, u dinamici konstrukcija izuzetnu važnost ima i princip diskretizacije. Diskretizacija masa ili redukcija masa nosača (lumped masses) je pojam koji se vezuje za kinetostatičku metodu koja je prva dala jednačine linearnih oscilacija složenih grednih nosača sa konačnim brojem stepeni slobode i koja se uglavnom koristila u radovima domaćih istraživača. Danas, apsolutno dominantna metoda u dinamici konstrukcija koja se vezuje za pojam diskretizacije jeste metoda konačnih elemenata. Ona je opšteprihvaćena kao najefikasnija i najpraktičnija metoda za analizu statičkih i DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 21 dinamičkih problema u gotovo svim poljima mašinstva. Razvijen je veliki broj softvera koji omogućavaju laku primenu ove metode i sprovođenje odgovarajućih analiza. Ovaj matematički‐numerički pristup u radovima iz oblasti problema pokretnog opterećenja se klasifikuje kao konačnoelementi pristup. Ovde mora biti napomenuto da ovaj pristup, apriori, ne podrazumeva korišćenje gotovih softverskih paketa za analizu dinamičkog ponašanja struktura već formiranje karakterističnih matrica prema usvojenoj diskretizaciji sistema: matrice inercije sistema, matrice prigušenja sistema, matrice krutosti sistema i vektora spoljašnjih sila. Uobičajeno je da se prigušenje razmatra preko teorije Rejlija. Rešenja jednačina se dobijaju metodama direktne integracije (step‐by step methods) pri čemu je metoda Njumarka veoma često u primeni. Može se zaključiti da konačnoelementni pristup iziskuje relativno jednostavan matematički aparat, zasnovan na matričnom računu, ali je obim programiranja izrazito veliki, kako za formiranje jednačina tako i za njihovo rešavanje. Čist konačnoelementni pristup u problematici pokretnog opterećenja je moguć samo za problem pokretne sile i za problem pokretne mase. Svi ostali modeli pokretnih sistema uključuju potrebu za jednačinama analitičke mehanike. Ovo se može nazvati kombinovanim pristupom kao pojmom koji opisuje metodologiju savremenih istraživanja problema pokretnog opterećenja. Odzivi modela u ovom radu su dobijeni na modelima portalne dizalice čije noseće konstrukcije u svakom slučaju predstavljaju složene strukture. Usvojen je kombinovani pristup za formiranje matematičkog modela sistema portalne dizalice, tj. konačnoelementni pristup je iskorišćen za modeliranje noseće konstrukcije portalne dizalice, a jednačine analitičke mehanike su iskorišćene za modeliranje kolica. Razmatraju se tri modela kolica, dva koja su standardna u problemima pokretnog opterećenja i jedan koji predstavlja originalan pristup‐model pokretnog oscilatora sa klatnom. Modeliranje kretanja kolica je izvršeno tako da je moguće razmatrati bilo koji profil brzina kretanja kolica. Postavljene su jednačine kretanja sistema kao sistem diferencijalnih jednačina drugog reda sa promenljivim koeficijentima. Rešenja jednačina su dobijena metodom Njumarka. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 22 2.3 ZAKLJUČNA RAZMATRANJA Na osnovu pregleda relevantne literature koji je dat u poglavlju 2.1. može se zaključiti da je broj radova koji tretiraju problem pokretnog opterećenja na strukturama dizalica veoma mali. Sa druge strane, dinamika dizalica je u savremenim naučnim radovima, uglavnom, upravo prisutna kroz implementaciju problema pokretnog opterećenja za dobijanje odgovarajućih dinamičkih odziva konstrukcija dizalica [42, 43, 44, 45, 48, 49]. U nastavku je dat detaljniji osvrt na modele koji su korišćeni i na analizu podataka u ovim malobrojnim istraživanjima u vezi sa dizalicama. Implementacija problema pokretnog opterećenja kod struktura mosnih dizalica, sa razmatranjem klaćenja tereta, je prvo pokazana u radu Oguamanam‐a, [42, (1998)]. Na slici 2.3 je prikazan model analiziranog sistema. Slika 2.3. Model mosne dizalice [42] Analiza rezultata je izvršena na primeru grednog nosača sa sledećim karakteristikama: ρ=8000 kg/m, E=2,117 1011 Pa, Lb=10 m, A=16 10‐4 m2, I=2,13 10‐7 m4. Kako je u radu rečeno da se radi o kvadratnom poprečnom preseku grednog nosača, može se lako zaključiti da su njegove dimenzije 4x4 cm. Na osnovu ovih podataka mogu se dati sledeće napomene: odnos dužine grede i dužine stranice preseka iznosi 1000/4=250, što višestruko prevazilazi odnos kojim se opisuju gredni nosači, pa je i opisivanje savijanja preko teorije Bernuli‐Ojlera ovde neadekvatno; statičkom analizom najvećeg ugiba samo od sopstvene težine ovog DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 23 grednog nosača (mase 128 kg) dobija se vrednost od 37 cm!; najveći dopušteni ugib dizalice (Lb/400) iznosi 2,5 cm, pa nije zadovoljen čak ni red veličina ovog ugiba. U radu je prikazana složena matematička formulacija spregnutih oscilacija strukture i klaćenja tereta, ali dobijeni rezultati nemaju praktičnog značaja za opisivanje dinamičkog ponašanja dizalica zbog gore navedenih napomena. Takođe, usvojeno je da se kolica kreću konstantnom brzinom čime je zanemaren realan ciklus kretanja. Isti autor je u radu [43, (2001)] razmatrao sličan model koji u odnosu na model iz rada [42] uključuje i poprečne oscilacije glavnog nosača u horizontalnoj ravni, slika 2.2. Dužina grede je popravljena na Lb=6 m, ali ne i ostale karakteristike nosača (momenti inercije su isti kao u [42]). Nema rezultata za pomeranje nosača usled pokretnog klatna, pa se ovaj rad ne može podvesti pod dinamiku dizalica već dinamiku tereta. Dobijeni rezultati za ugao klaćenja tereta, sa maksimalnim vrednostima od 0,03 rad (1,71o), nemaju praktičan značaj jer se kod dizalica dozvoljava njihanje tereta do 5o. Takođe, brzina kolica od 0,133 m/s, sa periodima ubrzanja i usporenja od 15 s, ne zahteva prikazani složeni matematički aparat za određivanje njihanja već je dovoljan model matematičkog klatna. Wu je u radu [44, (2004)] razmatrao odziv strukture portalne dizalice prema uticajima koji nastaju od kretanja kolica i njihanja tereta u vertikalnoj ravni yz. Šematski model je prikazan na slici 2.4. Slika 2.4 Model portalne dizalice; mT ‐masa kolica; msw ‐masa tereta; [44] DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 24 Početna pretpostavka u radu je svođenje problema pokretnog klatna na problem pokretne mase pri čemu je uveden pojam ekvivalentne pokretne mase u obliku ݉௘௤ሺݐሻ ൌ ்݉ ൅݉௦௪ cosଶ ߠሺݐሻ (oznake prema slici 2.4). Ovom aproksimacijom je zanemaren uticaj centrifugalne sile pri klaćenju tereta čime nije zadovoljen uslov da kompletni uticaji klaćenja tereta budu uzeti u obrzir. Čak i uz ovo uprošćenje, ne može se nikako prihvatiti korišćenje ove ekvivalentne mase za obuhvatanje uticaja inercije mase u pravcima ose x i ose z. Posledično, izrazi za sile Peqx(t), Peqy(t), Peqx(t)z su potpuno neadekvatni kao opterećenje strukture. Dodatno, njihanje tereta zavisi i od ubrzanja strukture na mestima kontakta kolica i strukture što je u ovom radu izbegnuto pomenutom aproksimacijom i pretpostavljanjem apsolutnih vrednosti ugla njihanja tereta čime nisu obuhvaćene spregnute oscilacije celog sistema. Ovim su rezultati odziva strukture dovedeni u pitanje, a koji su inače veoma mali (pomeranje središnje tačke nosača iznosi do 3mm) pa nije jasna njihova praktična upotreba. Ipak, potrebno je reći da je prikazan inovativan koncept obuhvatanja uticaja inercije mase definisanjem promenljivih koeficijenata u matrici inercije preko šetajuće matrice masa, pomoću metode konačnih elemenata, čime se model sa konačnim brojem stepeni slobode adekvatno prilagođava za opisivanje problema pokretnog opterećenja. U radu [45 (2008)], Wu je u proračun uključio i podužne oscilacije elemenata strukture, a klaćenje tereta se vrši u ravni yx, slika 2.5a. Ovde je ponovo uvedena ekvivalentna pokretna masa, u istom obliku kao u radu [44], što je neadekvatna aproksimacija zbog datih napomena. Na slici 2.5.b prikazane su karakteristike modela ramovske strukture koji je iskorišćen za analizu odziva. Međutim, razmatrane brzine do 100 m/s predstavljaju potpuno nerealne i nedostižne karakteristike kolica kod dizalica, a i kod bilo koje staze dužine 5,5 m. Na dijagramu koji daje analizu središnjeg pomeranja ramovske strukture preko modela pokretne mase i pokretne sile, [45, Fig. 6], nije jasno prikazana razlika za brzine do 5 m/s, a što bi moglo dati važne zaključke za realan domet brzina dizaličnih kolica u skoroj budućnosti. Takođe, prilikom razmatranja njihanja tereta date su nejasne jedinice za kružnu frekvenciju, tj. frekvenciju oscilovanja klatna, pa su odnosi frekvencija klatna i frekvencija strukture drugačiji i ne mogu poslužiti za donošenje zaključaka bilo koje vrste. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 25 a) b) Slika 2.5 a) Model portalne dizalice, b) Karakteristike KE modela strukture [45] U radu [48 (2011)], Yazid je razmatrao spregnute oscilacije strukture, klaćenja tereta i užetnog sistema preko modela koji je prikazan na slici 2.6. U radu je pretpostavljeno da su kolica stacionirana na sredini raspona, pa se ovaj rad ne može svrstati u oblast pokretnog opterećenja. Ovim je zanemaren realan ciklus kretanja sa ubrzanjem koje utiče na njihanje tereta u vertikalnoj ravni YX. Slika 2.6 Model portalne dizalice sa kolicima na sredini raspona [48] U radu [49 (2010)] se analizira struktura portalne dizalice sa prostornim klaćenjem tereta, slika 2.7. U odnosu na radove Wu‐a [44,45], ovde su uvedeni odgovarajući izrazi koji proizilaze iz dinamike tereta. U radu nije dat način DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 26 implementacije interakcije klatna i strukture u konačnoelementni model strukture, gde je staza modelirana sa samo 4 elementa, kao i podatak o uključivanju mešovitih članova koji proizilaze iz ubrzanja tačke na stazi gde se ostvaruje kontakt sa kolicima. U ovom radu je predstavljen model koji obuhvata veći broj parametara kod dizalica od prethodno pomenutih radova, sa aspekta problema pokretnog opterećenja, jer je objavljen 2010. godine. Slika 2.7 Model portalne dizalice [49] Na osnovu ovih zaključaka može se reći da je problem pokretnog opterećenja nedovoljno, a ponegde i nejasno, uključen u savremena istraživanja dinamičkog ponašanja dizalica. Generalno, sa daljim povećanjem fleksibilnosti struktura dizalica i sa tendencijom rasta brzina kretanja kolica, moguće je očekivati nastavak ovih istraživanja. U užem smislu, uz konstataciju da problem pokretne sile predstavlja dovoljno istraženo naučno poglavlje u dinamici struktura, problemi pokretnog opterećenja imaju utemeljenje u dinamici dizalica iz sledećih razloga:  masa tereta koji se transportuje je često veća od mase same konstrukcije dizalice pa je uticaj inercije mase važan parametar za dinamičku analizu  kolica mogu biti konstruktivno izvedena tako da veza između kućišta vitla i kolica nije kruta već ostvarena pomoću elastičnih opruga ili je veza ostvarena čvrstom vezom ali na mestu gde savojna krutost konstrukcije kolica dozvoljava elastična pomeranja  na svim dizalicama kačenje tereta na vitlo je ostvareno preko pomoćnog užetnog ili lančanog sistema čime je omogućeno klaćenje tereta DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 27 2.4 PLAN ISTRAŽIVANJA Osnovni plan istraživanja je prikazan na slici 2.8. i sprovešće se u saglasnosti sa usvojenim metodama za postavljanje i rešavanje diferencijalnih jednačina kretanja sistema, poglavlje 2.2. Početne pretpostavke za istraživanje naslovnog problema su sledeće: 1) Razmatraju se dinamički uticaji koji nastaju u toku ciklusa kretanja kolica 2) Razmatraju se opterećenja koja deluju u vertikalnoj ravni koja je upravna na pravac kretanja dizalice, kao i odgovarajući dinamički odzivi strukture portalne dizalice 3) Kolica su u stalnom kontaktu sa strukturom dizalice 4) Staza po kojoj se kreću kolica je idealno ravna Slika 2.8. Osnovni plan istraživanja u disertaciji DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 28 2.5 CILJ ISTRAŽIVANJA Osnovni ciljevi ove disertacije imaju praktični i naučni karakter. Afirmacija i razvoj kombinovanog pristupa pri formiranju modela dizalica predstavlja početni cilj ovog rada, čime se omogućava adekvatno istraživanje problema pokretnog opterećenja kod složenih struktura dizalica. Glavni ciljevi ovog rada, u skladu sa planom istraživanja, predstavljaju formiranje tri matematička modela portalne dizalice kroz trojako posmatranje kolica kod dizalica: model pokretne mase, model pokretnog oscilatora i model pokretnog oscilatora sa klatnom. Posledično, daje se komparativan prikaz načina modeliranja kolica kod portalnih dizalica. Zajedničke karakteristike ovih modela su sledeće:  pored poprečnih, uključene su i podužne oscilacije elemenata strukture  uključeni su uticaji inercije masa u dinamičkoj interakciji strukture i kolica  razmatraju se realni ciklusi kretanja kolica, sa ubrzanjem i usporenjem Matematičkim modelima su uključena u razmatranje dva konstrukciona tipa portalnih dizalica: osnovni i tip sa jednim prepustom. Ovo je važno istaći jer konstrukcija sa prepustom nije, kroz modele, zastupljena u literaturi (poglavlje 2.3.). Model pokretne mase, iako prisutan u literaturi iz ove oblasti, je ovde detaljno opisan i formiran jer predstavlja bazični model koji pored razmatranja uticaja inercije mase na dinamičko ponašanje strukture portalne dizalice služi i za verifikaciju rezultata dobijenih preko složenijih modela. U radovima [42,43,44,45] nisu prikazani verodostojni podaci koji mogu poslužiti za razmatranje dinamičkog ponašanja portalnih dizalica za kontejnerske terminale. Model pokretnog oscilatora ima za cilj opisivanje i određivanje uticaja elastičnog oslanjanja, u vertikalnom pravcu, u sistemu kolica kod portalnih dizalica. Sa aspekta problema pokretnog opterećenja, ovim se daje mogućnost koja nije prikazana u dosadašnjim modelima kod dizalica. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 29 Model pokretnog oscilatora sa klatnom, koji je predstavljen u ovoj disertaciji, je originalan pristup modeliranju pokretnog sistema u oblasti problema pokretnog opterećenja. Sa aspekta istraživanja kod dizalica, u odnosu na radove [44,45,49], može se slobodno reći da predstavlja nadgradnju postojećih modela jer je obuhvaćena elastična veza između kolica i vitla. Ovaj model je kombinacija modela pokretnog oscilatora i pokretnog klatna, a od važnijih karakteristika potrebno je pomenuti sledeće:  nije usvojeno oubičajeno ograničenje na male oscilacije klaćenja tereta  uključene su potpune karakteristike klaćenja tereta u vertikalnoj ravni Praktični ciljevi ove disertacije su korišćenje dobijenih rezultata pri projektovanju fleksibilnih struktura portalnih dizalica sa visokim performansama, a posebno kod zahteva za povećanjem brzina kretanja kolica preko 3 m/s. Takođe, postavljeni matematički modeli predstavljaju osnovu za formiranje algoritma upravljanja pogona kretanja kolica. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 30 3 MODELIRANJE NOSEĆIH KONSTRUKCIJA PORTALNIH DIZALICA SISTEMOM ELASTIČNIH TELA Na osnovu pregleda literature, poglavlje 2.1, može se primetiti da je analitički pristup veoma prisutan u problematici pokretnog opterećenja. Međutim, modeli ramovskih struktura uopšte nisu razmatrani što se objašnjava time da tek skorašnje performanse pokretnih sistema (velike brzine i ubrzanja) na ramovskih strukturama daju svrsishodnost upotrebe problema pokretnog opterećenja. Ovde će biti prikazan koncept primene analitičkog pristupa na dinamiku nosećih konstrukcija portalne dizalice. Prvo, geometrija noseće konstrukcije portalne dizalice i dimenzije preseka njenih elemenata apsolutno dozvoljavaju da se model formira kao sistem elastičnih tela tipa prizmatične grede. Idealizacija se zasniva na sledećim pretpostavkama:  materijal elemenata je homogen i izotropan  svi elementi su Bernuli‐Ojler grede, tj. razmatra se samo savijanje elemenata, a time i poprečne oscilacije greda  poprečna pomeranja središta preseka elemenata su upravna na podužnu osu i mala u odnosu na dužinu elementa  poprečni preseci elemenata ostaju ravni i upravni na elastičnu liniju Na slici 3.1 prikazani su tipovi konstrukcija portalnih dizalica koji će biti razmatrani: osnovni tip (slika 3.1.a) i tip sa jednim prepustom (slika 3.1.b). Osnovne geometrijske karakteristike koje su iskorišćenje za definisanje dužina elemenata‐prizmatičnih greda su: L ‐ raspon između nogu, H ‐ visina na strani DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 31 krute noge, h ‐ visina na strani elastične noge i Lp ‐ dužina prepusta glavnog nosača. Iako je veoma redak kod portalnih dizalica, uvođenjem parametra h takvim da je h≠H, obuhvaćen je slučaj kada noge portalnih dizalica nisu postavljene na istom nivou. Glavni elementi noseće konstrukcije su izrađeni od istog materijala (konstrukcionog čelika). Na slici 3.2 prikazani su usvojeni ramovski modeli sistema elastičnih tela koji odgovaraju postavljenim tipovima konstrukcija portalnih dizalica na slici 3.1. Svi konstrukcioni delovi portalne dizalice su modelirani kao elementi konstantnog poprečnog preseka sledećih karakteristika: modul elastičnosti E, gustina materijala ρ, moment inercije Ii i površina preseka Ai (i = 1, 2, 3, 4). Ovde će biti sproveden postupak definisanja poprečnih oscilacija prema tipu konstrukcije sa jednim prepustom, slika 3.2.b, pošto je postupak za osnovni tip noseće konstrukcije portalne dizalice detaljno prikazan u [54]. Slika 3.1. Tipovi noseće konstrukcije portalne dizalice, a) Osnovni, b) Sa jednim prepustom Slika 3.2. Ramovski modeli sistema elastičnih tela prema usvojenim tipovima konstrukcije na slici 3.1.a i 3.1.b. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 32 Parcijalna diferencijalna jednačina slobodnih transverzalnih oscilacija za svaki od elemenata glasi డమఔ೔డ௧మ ൅ ாூ೔ ఘ஺೔ డరఔ೔ డ௭೔ర ൌ 0 (i = 1, 2, 3, 4) (3.1) Opšta rešenja jednačine (3.1), tj. poprečna pomeranja elemenata, se mogu prikazati u sledećem obliku ߥଵ ൌ ߥଵሺݖଵ, ݐሻ ൌ ܼଵሺݖଵሻܶሺݐሻ 0 ൑ ݖଵ ൑ ܮ (3.2) ߥଶ ൌ ߥଶሺݖଶ, ݐሻ ൌ ܼଶሺݖଶሻܶሺݐሻ 0 ൑ ݖଶ ൑ ܪ (3.3) ߥଷ ൌ ߥଷሺݖଷ, ݐሻ ൌ ܼଷሺݖଷሻܶሺݐሻ 0 ൑ ݖଷ ൑ ݄ (3.4) ߥସ ൌ ߥସሺݖସ, ݐሻ ൌ ܼସሺݖସሻܶሺݐሻ 0 ൑ ݖସ ൑ ܮ௣ (3.5) Funkcije oblika oscilovanja su ovde postavljene preko Krilov‐Dankanovih funkcija [55], u sledećem obliku ܼ௜ሺݖ௜ሻ ൌ ܩ௜ܵሺ݇௜ݖ௜ሻ ൅ ܤ௜ܶሺ݇௜ݖ௜ሻ ൅ ܥ௜ܷሺ݇௜ݖ௜ሻ ൅ ܦ௜ܸሺ݇௜ݖ௜ሻ (i = 1,2,3,4) (3.6) gde su ki karakteristične vrednosti, a Gi, Bi, Ci, Di konstante koje se određuju iz konturnih uslova. Vremenska funkcija ima oblik ܶሺݐሻ ൌ ܺ cos߱ݐ ൅ ܻ sin߱ݐ (3.7) pri čemu su X i Y konstante koje se određuju iz početnih uslova kretanja, a kružna frekvencija oscilovanja iznosi ߱ ൌ ݇௜ଶටா ூ೔ఘ ஺೔ (i = 1, 2, 3, 4) (3.8) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 33 S obzirom da se struktura sastoji od 4 gredna elementa potrebno je formulisati 16 konturnih uslova. Na osloncu elementa 2 važe sledeći uslovi ܼଶሺ0ሻ ൌ 0 (3.9a) െܧܫଶܼଶ′′ሺ0ሻ ൌ 0 (3.9b) Na osloncu elementa 3 uslovi su ܼଷሺ݄ሻ ൌ 0 (3.9c) െܧܫଷܼଷ′′ሺ݄ሻ ൌ 0 (3.9d) Na slobodnom kraju elementa 4 uslovi su െܧܫସܼସ′′ሺ0ሻ ൌ 0 (3.9e) െܧܫସܼସ′′′ሺ0ሻ ൌ 0 (3.9f) Na mestima veze elemenata 1 i 3, kao i elemenata 1,2 i 4 dobijaju se ܼଵሺܮሻ ൌ 0 (3.9g) ܼଵ′ ሺܮሻ ൌ ܼଷ′ ሺ0ሻ (3.9h) െܧܫଵܼଵ′′ሺܮሻ ൌ െܧܫଷܼଷ′′ሺ0ሻ (3.9j) ܼଵሺ0ሻ ൌ 0 (3.9k) ܼସ൫ܮ௣൯ ൌ 0 (3.9l) ܼଶ′ ሺܪሻ ൌ ܼଵ′ ሺ0ሻ (3.9m) ܼସ′ ൫ܮ௣൯ ൌ ܼଵ′ ሺ0ሻ (3.9n) െܧܫଵܼଵ′′ሺ0ሻ ൌ െܧܫଶܼଶ′′ሺܪሻ െ ܧܫସܼସ′′൫ܮ௣൯ (3.9o) ܼଶሺܪሻ ൌ െܼଷሺ0ሻ (3.9p) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 34 Poslednji (i ključni) uslov se dobija primenom zakona o kretanju centra mase glavnog nosača (elemenata 1 i 4) pri dejstvu transverzalnih sila na mestima veze sa krutom i elastičnom nogom (elementima 2 i 3), odnosno െ൫ߩܣଵܮ ൅ ߩܣସܮ௣൯ߥሷଷሺ0, ݐሻ ൌ ܧܫଶߥଶ′′′ሺܪ, ݐሻ ൅ ܧܫଷߥଷ′′′ሺ0, ݐሻ (3.9r) koji, preko (3.3, 3.4) i (3.7), postaje ൫ߩܣଵܮ ൅ ߩܣସܮ௣൯ܼଷሺ0ሻ߱ଶ ൌ ܧܫଶܼଶ′′′ሺܪሻ ൅ ܧܫଷܼଷ′′′ሺ0ሻ (3.10) Iz (3.8), karakteristične vrednosti ki (i=2,3,4), se mogu izraziti preko k1, odnosno ݇௜ ൌ ݇ଵට஺೔஺భ ூభ ூ೔ ర (3.11) Na osnovu postavljenih uslova (3.9) i korišćenjem relacija za karakteristične vrednosti (3.11) dobija se karakteristična jednačina koja je veoma složena zbog kombinacije trigonometrijskih i hiperboličnih funkcija i zavisna od velikog broja parametara ܣଵ, ܣଶ, ܣଷ, ܣସ, ܫଵ, ܫଶ, ܫଷ, ܫସ, ܮ, ܮ௣, ܪ, ݄, oblika ݂ሺ݇ଵሻ ൌ 0 (3.12) koji ne dozvoljava eksplicitno određivanje frekvencija u analitičkom obliku. Rešavanje ove karakteristične jednačine, po karakterističnoj vrednosti k1, je omogućeno tek u poslednjih deset godina, sa pojavom naprednih matematičkih programa (npr. Wolfram Mathematica 6). Zbog prirode trigonometrijskih i hiperboličnih funkcija koje su sastavni deo karakteristične jednačine, rešenja je potrebno prvo grafički prikazati, lokalizovati pa tek onda odrediti tačna rešenja. Ovde će biti date prve tri frekvencije oscilovanja za usvojene varijante parametara portalnih dizalica za kontejnerske terminale. Glavne geometrijske mere su definisane preko L i H, a pretpostavljeno je da Lp = 0,25 L i h=H, što inače DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 35 predstavlja veliki broj izvođenja ovih dizalica. S obzirom da je dimenzionisanje glavnog nosača prvi korak u projektovanju portalne dizalice na osnovu usvojenih glavnih mera kojim je definisana njena namena, upravo su statičke karakteristike glavnog nosača usvojene kao polazne i određene prema preporukama za izbor dimenzija poprečnog preseka [56] tako da prezentuju fleksibilnu i krutu varijantu grednog nosača. Sa pretpostavkom da je I4=I1 i A4=A1 (prepust je isti nosač kao i rasponi deo glavnog nosača), varijacija statičkih parametara elemenata krute i pendel noge je izvršena preko sledećih izraza ߙ ൌ ூభூమ , ߚ ൌ ூభ ூయ, ߛ ൌ ஺భ ஺మ, ߜ ൌ ஺భ ஺య, (3.13) U svemu prema prethodnom, prikaz prve tri frekvencije za usvojene varijante izvođenja portalnih dizalica je dat u tabeli 3.1. Sagledavanjem rezultata iz tabele 3.1., u okviru jedne geometrijske varijante, može se potvrditi da fleksibilna konstrukcija generiše niže frekvencije oscilovanja u odnosu na krutu konstrukciju i da osnovna (prva) frekvencija drastično opada sa smanjenjem krutosti krute i pendel noge. Dobijeni rezultati mogu biti iskorišćeni za verifikaciju frekvencija ramovskih struktura dobijenih korišćenjem MKE. Sa ovim konstatacijama može se zaključiti primena analitičkog pristupa zbog složenosti definisanja prinudnih oscilacija ramovskih struktura usled dejstva pokretnog opterećenja. Zbog zaključaka koji su dati u metodologiji istraživanja, poglavlje 2.2, u nastavku će biti dat drugačiji pristup za razmatranje problema pokretnog opterećenja kod dizalica. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 36 Tabela 3.1. Prve tri frekvencije oscilovanja usvojenih varijanti portalnih dizalica L = 30 m H = 18 m I1= 0,024 m4 A1= 0,07 m2 I1= 0,05 m4 A1= 0,09 m2 α γ β δ f1 [Hz] f2 [Hz] f3 [Hz] f1 [Hz] f2 [Hz] f3 [Hz] 1 1 1 1 1,57 7,58 15,69 2,00 9,65 19,97 1 1 10 2 1,12 6,49 10,34 1,43 8,27 13,16 10 2 10 2 0,69 5,25 10,11 0,88 6,68 12,87 10 2 50 4 0,54 5,09 6,64 0,68 6,48 8,45 L = 40 m H = 20 m I1= 0,05 m4 A1= 0,09 m2 I1= 0,13 m4 A1= 0,11 m2 α γ β δ f1 [Hz] f2 [Hz] f3 [Hz] f1[Hz] f2 [Hz] f3 [Hz] 1 1 1 1 1,45 5,76 12,57 2,12 8,41 18,34 1 1 10 2 1,03 4,92 10,26 1,50 7,18 14,96 10 2 10 2 0,65 3,90 9,67 0,95 5,70 14,10 10 2 50 4 0,51 3,77 6,70 0,74 5,51 9,81 L = 45 m H = 18 m I1= 0,07 m4 A1= 0,1 m2 I1= 0,18 m4 A1= 0,14 m2 α γ β δ f1 [Hz] f2 [Hz] f3 [Hz] f1 [Hz] f2 [Hz] f3[Hz] 1 1 1 1 1,74 5,40 12,08 2,36 7,34 16,38 1 1 10 2 1,23 4,61 11,52 1,67 6,25 15,67 10 2 10 2 0,81 3,58 9,86 1,13 4,88 13,66 10 2 50 4 0,63 3,44 9,11 0,85 4,66 12,34 DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 37 4 FORMIRANJE MODELA PORTALNE DIZALICE PRIMENOM KOMBINOVANOG PRISTUPA Glavni delovi portalne dizalice sa aspekta postavljenog problema su noseća čelična konstrukcija i kolica. Postavljanje odgovarajućeg modela portalne dizalice se sastoji iz dva odvojena dela, tj. postavljanje modela za noseću konstrukciju dizalice‐strukturu i postavljanje modela za kolica. 4.1 MODEL NOSEĆE KONSTRUKCIJE PORTALNE DIZALICE Noseća konstrukcija portalne dizalice je prostorna struktura čija složenost zavisi od tipa konstrukcije, načina izvođenja konstrukcije i broja glavnih nosača. Princip diskretizacije pri formiranju ekvivalentnih modela je obavezan u savremenim istraživanjima iz oblasti dinamike struktura portalnih dizalica [31,44,45,48], pa je i ovde izvršeno svođenje na sistem sa konačnim brojem stepeni slobode. Na osnovu početne pretpostavke da se razmatraju samo opterećenja i dinamički odzivi strukture u vertikalnoj ravni koja je upravna na pravac kretanja dizalice, prostorni sistem noseće konstrukcije je preveden u ravanski model. Koncept formiranja uprošćenog (idealizovanog) modela strukture portalne dizalice pomoću principa diskretizacije je prikazan na slikama 4.1 i 4.2. Odnosno, razmatraju se dva konstrukciona tipa portalnih dizalica zbog toga što predstavljaju najčešće izvođenje portalnih dizalica, a uključene su specifičnosti svakog od tipova koje utiču na karakter statičkog/dinamičkog ponašanja strukture. Osnovni konstrukcioni tip portalne dizalice je ovde nazvan model A, slika 4.1. Glavne geometrijske karakteristike ovog P‐ramovskog diskretizovanog modela DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 38 su: Lgn ‐ dužina glavnog nosača između sistemskih linija nogu, H ‐ visina na strani krute noge koja predstavlja dužinu od sistemske linije na glavnom nosaču do mesta kontakta točkova čeonog nosača krute noge sa šinom i h ‐ visina na strani pendel noge koja predstavlja dužinu od sistemske linije na glavnom nosaču do mesta kontakta točkova čeonog nosača pendel noge sa šinom. Sledeći konstrukcioni tip portalne dizalice je ovde nazvan model B, slika 4.2, i razlikuje se od osnovnog modela portalne dizalice po postojanju prepusta glavnog nosača. Glavne geometrijske karakteristike ovog ramovskog diskretizovanog modela su: L ‐ dužina rasponog dela glavnog nosača između sistemskih linija nogu, Lp ‐ dužina prepusta glavnog nosača, a H i h su iste kao i za model A. Slika 4.1 Diskretizovani model strukture osnovnog tipa konstrukcije portalne dizalice‐model A Slika 4.2 Diskretizovani model strukture tipa konstrukcije portalne dizalice sa jednim prepustom‐model B DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 39 Za oba modela strukture portalne dizalice je usvojena podela glavnog nosača na 10 elemenata, krute noge na 2 elementa i pendel noge na 2 elementa. Na ovaj način je u dovoljnoj meri uključena specifičnost konstrukcije sa aspekta dinamičkog proračuna. Kao što je napomenuto u metodologiji (poglavlje 2.2), savremena istraživanja iz oblasti analize dinamičkog ponašanja nosećih konstrukcija portalnih dizalica se uglavnom baziraju na metodi konačnih elemenata (MKE), što će i ovde biti iskorišćeno za definisanje segmenata diskretizovane strukture portalne dizalice prema modelu A i modelu B. 4.1.1 Konačnoelementni model strukture Oba konačnoelementna modela strukture portalne dizalice se sastoje od 15 čvorova i 14 elemenata. Elementi glavnog nosača su označeni brojevima n =1‐10, elementi krute noge sa n =11 i 12, a elementi pendel noge sa n=13 i 14. Zbog svoje sličnosti oba modela su prikazana na slici 4.3. Slika 4.3. Konačnoelementni modeli strukture portalne dizalice Svaki čvor i (i=1‐15) ima 3 stepena slobode, tj. horizontalno pomeranje, vertikalno pomeranje i rotaciju oko ose upravne na ravan modela. Ovo je inicijalna postavka i za formiranje matrice inercije i matrice krutosti modela strukture, tj. uključeno je 45 mogućih pomeranja. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 40 Struktura je oslonjena preko nepokretnih zglobnih oslonaca u čvorovima 13 i 15, pa je ovim čvorovima onemogućeno pomeranje u vertikalnom i horizontalnom pravcu. Kako u oscilovanju strukture učestvuju samo elementi sa slobodnim pomeranjima čvorova model strukture ima efektivno 41 stepen slobode, tj. 41 čvornih pomeranja, zbog nepokretnih oslonaca u čvorovima 13 i 15 (blokirana su po 2 vertikalna i horizontalna pomeranja), što je prikazano na slici 4.4. Slika 4.4. Pomeranja čvorova KE modela; a) čvor i=1‐10, 12 i 14, b) čvorovi 13 i 15 Uvođenjem MKE u dinamičku analizu sistema, pomeranja čvorova strukture postaju generalisane koordinate. Dakle, vektor pomeranja konačnoelementnog modela strukture dizalice‐ vektor generalisanih koordinata strukture se može predstaviti kao ܃ ൌ ሼܷଡ଼ଵ ܷଢ଼ଵ θܷଵ ܷଡ଼ଶܷଢ଼ଶ …ܷଢ଼ଵଶ θܷଵଶ θܷଵଷ … ௑ܷଵସܷଢ଼ଵସ θܷଵସ θܷଵହሽT (4.1) a prvi i drugi izvod po vremenu komponenata ovog vektora čine komponente vektora brzina i ubrzanja, ܃ሶ i ܃ሷ , respektivno. Dužine elemenata odgovaraju ravnomernoj diskretizaciji elemenata strukture i za model A su definisane sledećim izrazima: ln = Lgn/10 (n =1‐10) (4.2a) ln =H/2 (n =11,12) (4.2b) ln = h/2 (n =13,14) (4.2c) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 41 Za model B, dužine elemenata su sledeće: ln = Lp/2 (n =1,2) (4.3a) ln = L/8 (n =3‐10) (4.3b) ln =H/2 (n =11,12) (4.3c) ln = h/2 (n =13,14) (4.3d) Pri formiranju KE modela B je dodatno uvedena pretpostavka da je dužina prepusta glavnog nosača u zavisnosti od dužine rasponog dela glavnog nosača sa odnosom Lp = L/4. Ovo je uvedeno da bi elementi glavnog nosača (n =1‐10) imali istu dužinu, tj. da bi se staza po kojoj se kreće pokretno opterećenje sastojala od 10 elemenata istih dužina. Na ovaj način se lako može KE model B preformulisati u KE model A, sa preformulacijom elemenata 11 i 12 u odgovarajućim matricama. Ova pretpostavka ima utemeljenje u statičkom ponašanju prepusta posmatranog kao konzolnog elementa koji, inženjerski gledano, može da bude opterećen četvrtinom opterećenja grednog nosača koji prezentuje rasponi deo glavnog nosača da bi se održalo slično naponsko stanje u nosaču sa aspekta dominantnog naprezanja konstrukcija dizalica‐savijanja. Ovo se može izbeći podelom staze na veći broj delova istih dužina, ali se na taj način povećava i programerski obim operacija koji su obavezni za postavljanje problema dinamičke interakcije strukture i kolica primenom kombinovanog pristupa, kao i vreme potrebno za dobijanje dinamičkih odziva u programskom paketu. Diskretizacija dela strukture po kojem se kreće pokretno opterećenja na 10 elemenata je uobičajena mera diskretizacije u istraživanjima. Sa praktičnog i inženjerskog aspekta se može zaključiti da su postavljeni KE modeli strukture relativno jednostavni ali oni u svakom slučaju adekvatno oslikavaju osnovne tipove konstrukcija portalnih dizalica i daju mogućnost dovoljno tačnih statičkih i dinamičkih analiza, što će biti pokazano verifikacijom rezultata. Takođe, ispunjene su sve pretpostavke za adekvatnu implementaciju problema pokretnog opterećenja, [36], čime je ispunjen uslov za istraživanje naslovnog problema ove disertacije. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 42 Na ilustracijama koje slede u narednim poglavljima će biti predstavljen samo model B kao model koji predstavlja opštiji konstrukcioni tip portalne dizalice u odnosu na model A. Ovaj model se može lako prilagoditi da odgovara ostalim konstrukcionim tipovima portalnih dizalica (Slika 4.5.a,b,c). Takođe, ovaj model može biti prilagođen za modeliranje osnovnog tipa obalske kontejnerske dizalice (slika 4.5.d). Ovo se može izvesti bilo preformulacijom elemenata postojećeg modela (slika 4.5.a) ili dodavanjem elemenata i čvorova u model (slika 4.5.b,c) ili kombinacijom ova dva pristupa (slika 4.5.d). Primena ovih modifikacija bi podrazumevala da je unapred poznat i definisan tip konstrukcije portalne dizalice što ima smisla samo u konkretnim situacijama‐problemima dizajna portalnih dizalica. Algoritam za opisivanje dinamičke interakcije strukture i kolica portalne dizalice, koji će biti izložen u narednim poglavljima, je apsolutno primenljiv i za ovako modifikovane modele. Slika 4.5. (a), (b), (c) Skice portalnih dizalica prema tipu konstrukcije ili načinu izvođenja konstrukcije, (d) skica osnovnog tipa obalske kontejnerske dizalice DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 43 4.1.2 Tip usvojenog konačnog elementa Struktura portalne dizalice je podeljena na 14 elemenata koji su ovde predstavljeni linijskim KE koji može biti napregnut na aksijalno naprezanje i na savijanje. Osnovne pretpostavke za usvajanje ovog elementa su sledeće:  aksijalne deformacije elementa su saglasne Hukovom zakonu  poprečne deformacije elementa su saglasne teoriji Bernuli‐Ojlera Ovaj element se u domaćoj literaturi, [57], naziva gredni KE element (plane‐frame element), tj. kompletni linijski KE za ravansku strukturnu analizu. Pomenuti KE je kombinacija ravanskog elementa tipa štapa i elementa tipa nosača, tj. čvorovi elementa imaju (u lokalnom koordinatnom sistemu elementa) aksijalno pomeranje, poprečno pomeranje i rotaciju oko ose upravne na ravan elementa, slika 4.6. Vektor pomeranja početnog i krajnjeg čvora n‐tog elementa u lokalnom koordinatnom sistemu (slika 4.6) se može predstaviti kao ܝ ൌ ሾݑ௜ ݒ௜ ߠ௜ ݑ௝ ݒ௝ ߠ௝ሿ் (4.4) Slika 4.6. Pomeranja grednog KE u lokalnom koordinatnom sistemu DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 44 Osnovne mehaničke i statičke karakteristike n‐tog elementa su:  n ‐ gustina materijala,  E ‐ modul elastičnosti,  nA ‐ površina preseka i  nI ‐ aksijalni moment inercije preseka za glavnu osu upravno na ravan. S obzirom da su pomeranja čvorova nepoznate veličine, neophodno je uspostaviti njihovu direktnu vezu sa veličinom pomeranja u bilo kojoj tački polja elementa pomoću interpolacionih funkcija [57]. One su egzaktno određene za gredni linijski KE, a ovde će biti prikazane jer su od izuzetne važnosti za formiranje dinamičkih jednačina postavljenog modela. S obzirom na samu postavku ovog tipa KE jasno je da pomeranje u polju ima svoje dve komponente: wx(x) ‐ pomeranje u pravcu podužne ose elementa (slika 4.7.a) i wy(x) ‐ pomeranje poprečno na osu elementa (slika 4.7.b). Slika 4.7. Pomeranja proizvoljne tačke u polju elementa na rastojanju x od čvora i; (a) pomeranje u pravcu podužne ose elementa, (b) pomeranje poprečno na osu elementa DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 45 Dakle, aksijalno pomeranje bilo koje tačke u polju štapa, na rastojanju x od levog čvora, se može predstaviti kao ݓ௫ሺݔሻ ൌ ଵܰሺݔሻݑ௜ ൅ ସܰሺݔሻݑ௝ (4.5) gde su ଵܰሺݔሻ ൌ 1 െ ௫௟ (4.6) ସܰሺݔሻ ൌ ௫௟ (4.7) Poprečno pomeranje u polju grednog nosača se može predstaviti kao ݓ௬ሺݔሻ ൌ ଶܰሺݔሻݒ௜ ൅ ଷܰሺݔሻߠ௜ ൅ ହܰሺݔሻݒ௝ ൅ ଺ܰሺݔሻߠ௝ (4.8) gde su ଶܰሺݔሻ ൌ 1 െ 3 ቀ௫௟ቁ ଶ ൅ 2 ቀ௫௟ቁ ଷ (4.9) ଷܰሺݔሻ ൌ ൬௫௟ െ 2 ቀ ௫ ௟ቁ ଶ ൅ ቀ௫௟ቁ ଷ൰ ݈ (4.10) ହܰሺݔሻ ൌ 3 ቀ௫௟ቁ ଶ െ 2 ቀ௫௟ቁ ଷ (4.11) ଺ܰሺݔሻ ൌ ൬െቀ௫௟ቁ ଶ ൅ ቀ௫௟ቁ ଷ൰ ݈ (4.12) Korišćenjem izraza (4.9)‐(4.12) (Hermit‐ovih kubnih polinoma) za definisanje poprečnog pomeranja elementa, ostvaren je pun metod, [14]. Prost metod, postavljanjem N2=N1, N5=N4 i N3=N6=0, daje znatno lakši algoritam definisanja problema i zanemarljive razlike u rezultatima, ali ovi zaključci važe samo za problem pokretne sile. S obzirom na važnost prvog izvoda funkcija N1 i N4, kao i drugog izvoda funkcija N2, N3, N5 i N6 za povezivanje MKE sa Teorijom elastičnosti, [58], u tabeli 4.1. je dat prikaz karakterističnih izvoda ovih funkcija. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 46 Tabela 4.1 Prvi i drugi izvodi interpolacionih funkcija ௜ܰ ݀ ௜ܰሺݔሻ݀ݔ ൌ ௜ܰ ′ ݀ ଶ ௜ܰሺݔሻ ݀ݔଶ ൌ ௜ܰ ′′ ଵܰ െ 1݈ 0 ଶܰ 1݈ ൬െ6 ݔ ݈ ൅ 6 ቀ ݔ ݈ ቁ ଶ ൰ 1݈ଶ ቀെ6 ൅ 12 ݔ ݈ ቁ ଷܰ 1 െ 4 ݔ݈ ൅ 3 ቀ ݔ ݈ ቁ ଶ 1݈ ቀെ4 ൅ 6 ݔ ݈ ቁ ସܰ 1݈ 0 ହܰ 1݈ ൬6 ݔ ݈ െ 6 ቀ ݔ ݈ ቁ ଶ ൰ 1݈ଶ ቀ6 െ 12 ݔ ݈ ቁ ଺ܰ െ2 ݔ݈ ൅ 3 ቀ ݔ ݈ ቁ ଶ 1݈ ቀെ2 ൅ 6 ݔ ݈ ቁ Na osnovu postavljenih interpolacionih funkcija, određuje se matrica krutosti ovog tipa KE u lokalnom koordinatnom sistemu [59], u sledećem obliku ܓ௡ ൌ 3 2 3 2 2 2 3 2 3 2 2 2 0 0 0 0 12 6 12 60 0 6 4 6 20 0 0 0 0 0 12 6 12 60 0 6 2 6 40 0                       n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n n EA EA l l EI EI EI EI l l l l EI EI EI EI l l l l EA EA l l EI EI EI EI l l l l EI EI EI EI l l l l (4.13) Za dinamičku analizu, potrebno je odrediti i matricu inercije ovog tipa KE korišćenjem interpolacionih funkcija (čime ona postaje konzistentna matrica) i ima sledeći oblik [60] DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 47 ܕ௡ ൌ ఘ೙஺೙௟೙ସଶ଴ 2 2 2 2 140 0 0 70 0 0 0 156 22 0 54 3 0 22 4 0 13 3 70 0 0 140 0 0 0 54 13 0 156 22 0 13 3 0 22 4              n n n n n n n n n n n n l l l l l l l l l l l l (4.14) Za potrebe transformacije iz lokalnog u globalni koordinatni sistem, slika 4.6, koristi se transformaciona matrica u sledećem obliku ܂௡= 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1                   cos sin ‐sin cos cos sin ‐sin cos (4.15) 4.1.3 Matrica krutosti KE modela strukture Matrice krutosti elemenata u globalnom koordinatnom sistemu su iste kao matrice krutosti u lokalnom sistemu za elemente staze jer se ose ovih sistema poklapaju, a za elemente nogu se dobijaju pomoću transformacione matrice (4.15) za ugao  3π/2, odnosno ۹௡ ൌ ܓ௡, n = 1‐10 (4.16) ۹௡ ൌ ܂௡் ܓ௡܂௡ n =11‐14 (4.17) Proširivanjem i usaglašavanjem ovih matrica [61] dobija se globalna matrica krutosti modela strukture portalne dizalice kao ሾܭ௦௧ሿସହൈସହ ൌ ∑ ۹௡ଵସଵ (4.18) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 48 Eliminacijom vrsta i kolona u globalnoj matrici krutosti koji odgovaraju nepokretnim osloncima, saglasno napomenama datim u poglavlju 4.1.1., dobija se kvadratna matrica krutosti modela čiji je red 41, u svom konačnom obliku ۹ܛܜ ൌ ሾܭ௦௧ሿସଵൈସଵ (4.19) 4.1.4 Matrica inercije KE modela strukture Matrica inercije modela strukture portalne dizalice se dobija analogno principu za dobijanje matrice krutosti iznetim u prethodnom poglavlju, pa imamo ۻ௡ ൌ ܕ௡, n = 1‐10 (4.20) ۻ௡ ൌ ܂௡்ܕ௡܂௡, n =11‐14 (4.21) Globalna matrica inercije strukture dizalice iznosi ሾܯ௦௧ሿସହൈସହ ൌ ∑ ۻ௡ଵସଵ (4.22) U konačnom obliku matrica inercije se može predstaviti kao ۻܛܜ ൌ ሾܯ௦௧ሿସଵൈସଵ (4.23) Matrice inercije i matrice krutosti struktura u ovom radu su formirane u programu Wolfram Mathematica 6, gde se na osnovu glavnih mera dizalice i statičkih karakteristika elemenata generišu elementi matrica prema modelu A ili B. 4.1.5 Slobodne neprigušene oscilacije strukture Saglasno postavkama iznetim u prethodnim poglavljima moguće je razmatranje slobodnih neprigušenih oscilacija KE modela strukture portalne dizalice, tj. konzervativnog sistema sa konačnim brojem stepeni slobode, u svemu prema poznatoj jednačini DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 49 ۻܛܜ܃ሷ ൅ ۹ܛܜ܃ ൌ 0 (4.24) gde su ܃ሷ ‐ vektor ubrzanja generalisanih koordinata strukture ܃ ‐ vektor generalisanih koordinata strukture Na osnovu izraza (4.24) postavlja se frekventna jednačina koja služi za određivanje kružnih sopstvenih frekvencija, u obliku ‖۹ܛܜ െ ߱ଶۻܛܜ‖ ൌ 0 (4.25) Iz izraza (4.25) dobija se spektar kružnih frekvencija ߱௜, i =1‐41, a na osnovu toga i odgovarajući periodi oscilovanja i sopstvene frekvencije strukture, prema izrazima ௜ܶ ൌ ଶగఠ೔ (4.26) ௜݂ ൌ ଵ்೔ (4.27) Najniža frekvencija se naziva i osnovna frekvencija i zajedno sa nekoliko prvih sledećih najnižih frekvencija ima najveći značaj za analizu dinamičkog odziva strukture. Verifikacija programa za postavljanje diskretizovanog KE modela strukture je izvršena preko rezultata dobijenih razmatranjem oscilacija strukture kao sistema elastičnih tela tipa prizmatične grede, glava 3. Prema datim pretpostavkama, za L=30 m, H=18 m, I1=0,024 m4, A1=0,07 m2 i ߙ=1, ߛ=1, ߚ=10, ߜ=2 dobijene su frekvencije oscilovanja f1=1,122 Hz, f2=6,45 Hz, f3=10,386 Hz, f4=15,43 Hz...itd. Upoređenjem sa prve tri frekvencije koje su dobijene za odgovarajuću varijantu, tabela 3.1, može se konstatovati odlično poklapanje rezultata sa odstupanjima manjim od 0,7 %. Na slici 4.8. prikazana su prva dva oblika oscilovanja strukture portalne dizalice prema navedenim parametrima. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 50 Slika 4.8. a) Osnovni oblik oscilovanja, f1=1,12 Hz, b) 2. mod, f2=6,45 Hz 4.1.6. Matrica prigušenja KE modela strukture Određivanje prigušenja u strukturi dizalice je moguće odrediti jedino eksperimentalnim putem na pojedinačnoj konstrukciji u realnim radnim uslovima ili eventualno na umanjenom modelu (scale model) dizalice. Ovaj rad je, kao i većina radova iz ove oblasti, posvećen analizi strukture portalne dizalice kao fazi pre konstruisanja i pre izrade dizalice. Sa aspekta postavke što opštijeg matematičkog modela za analizu dinamičkog ponašanja portalne dizalice, potrebno je uključiti i mogućnost uticaja prigušenja u strukturi. Ovo se izvodi preko teorije lorda Rejlija (Rayleigh), [62], i matrica prigušenja strukture se postavlja u sledećem obliku ۱ܛܜ ൌ ߙۻܛܜ ൅ ߚ۹ܛܜ (4.28) tj. u direktnoj zavisnosti od matrice inercije i matrice krutosti KE modela strukture dizalice. Generalno, konstante ߙ i ߚ se određuju za bilo koje dve frekvencije na sledeći način ߙ ൌ ଶఠ೔ఠೕሺక೔ఠೕିకೕఠ೔ሻఠೕమିఠ೔మ (4.29) ߚ ൌ ଶሺకೕఠೕିక೔ఠ೔ሻఠೕమିఠ೔మ (4.30) gde su ߱௜, ௝߱ kružne frekvencije oscilovanja, a ߦ௜, ߦ௝ odgovarajući koeficijenti prigušenja. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 51 U ovom radu, ove konstante su određene prema prve dve frekvencije, jer su one i najbitnije u dinamici dizalica. Implementacija prigušenja u model podrazumeva da je izvršena modalna analiza neprigušenih oscilacija strukture. Koeficijenti prigušenja se razmatraju u opsegu 0,5‐7 %, [63]. Dodatno, može se predpostaviti da su koeficijenti prigušenja za obe kontrolne frekvencije isti, tj. ߦ௜ ൌ ߦ௝ ൌ ߦ. 4.2 MODELI KOLICA Kolica predstavljaju osnovni funkcionalni sistem na portalnoj dizalici. Ovaj sistem predstavlja jednu kompaktnu mašinu koja se sastoji od velikog broja električnih i mašinskih delova. Teret koji se transportuje se, u svim analizama, posmatra kao sastavni deo ovog sistema, tj. razmatra se slučaj punog opterećenja portalne dizalice. Slika 4.9. Kolica portalne dizalice za transport kontejnera Idealizacija kolica je izvršena do nivoa gde se ovaj sistem sastoji od tri glavna dela (slika 4.9): 1‐kolica (odnosi se na konstrukciju kolica sa točkovima), 2‐ vitla i 3‐tereta. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 52 Sa aspekta dinamike ovaj sistem se može posmatrati kao vezani materijalni sistem, što će ovde biti i pokazano. Veze između ovih delova su ostvarene elastičnim ili krutim vezama, u zavisnosti od načina modeliranja i istraživanja postavljenog problema. U skladu sa terminologijom koja se koristi u problemima pokretnog opterećenja i ciljem istraživanja u ovom radu, kolica će biti posmatrana kao:  Pokretna sila  Pokretna masa  Pokretni oscilator  Pokretni oscilator sa klatnom Model pokretne sile nije eksplicitno obrađen u ovom radu, ali se on može dobiti preko modela pokretne mase ukoliko se zanemare inercijalni efekti pokretne mase. Za postavljanje konkretnih modela biće iskorišćene sledeće pretpostavke: 1. Opterećenje kojim pokretni sistem deluje na glavne nosače je ravnomerno raspodeljeno na glavne nosače, ukoliko se razmatra dvogreda portalna dizalica. Ovo je neophodna pretpostavka jer je usvojen ravanski ramovski model strukture portalne dizalice, poglavlje 4.1. Ova pretpostavka se kod dizalica može smatrati opravdanom jer su razlike u pritiscima točkova veoma male što se postiže odgovarajućim sistemom namotavanja užadi i konstruktivnim parametrima kolica. 2. Razmak točkova na kolicima b, [1], je dovoljno manji od dužine staze, odnosno b ≪ Lgn, i opterećenje kojim pokretni sistem deluje na strukturu se može predstaviti tako da deluje u jednoj tački. 3. Osnovni delovi kolica su predstavljeni kao materijalne tačke, gde je m1 ‐ masa kolica, m2 ‐ masa vitla i m3 ‐ masa tereta. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 53 4.2.1 Pokretna masa Pokretnom masom pretpostavljeno je da je celokupan sistem kolica jedna koncentrisana masa, mss, koja se sastoji od pojedinačnih masa delova kolica ݉௦௦ ൌ ݉ଵ ൅݉ଶ ൅݉ଷ (4.31) i koja se kreće po strukturi, slika 4.10.a. Slika 4.10. Model kolica: (a) ‐ model pokretne mase, (b) ‐ model pokretnog oscilatora, (c) ‐ model pokretnog oscilatora sa klatnom 4.2.2. Pokretni oscilator Ovim modelom predstavljen je celokupni sistem kolica kao sistem od dve koncentrisane mase, m1 i m23, koje su međusobno spojene oprugom koeficijenta krutosti k. Masa m1 je masa kolica, a masa m23 predstavlja masu vitla i tereta, odnosno ݉ଶଷ ൌ ݉ଶ ൅݉ଷ (4.32) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 54 Ovaj model, slika 4.10.b, ima jedan dodatni stepen slobode, u odnosu na pokretnu masu sa slike 4.10.a, i to vertikalno pomeranje mase m23, ovde označeno sa koordinatom y, u odnosu na položaj statičke ravnoteže. Ovaj model ima cilj da opiše uticaj elastičnog ovešenja vitla sa teretom na dinamičko ponašanje strukture portalne dizalice. 4.2.3. Pokretni oscilator sa klatnom Poslednji u nizu modela sistema kolica koji će ovde biti razmatran se sastoji od tri koncentrisane mase m1, m2 i m3, odnosno, mase kolica, vitla i tereta, respektivno. Mase m1 i m2 su spojene elastičnom oprugom krutosti k, čime je simulirano elastično vešanje između ovih dvaju elemenata. Mase m3 i m2 su spojene neistegljivim lakim štapom dužine Lu, čime je predstavljen užetni sistem, a omogućeno klaćenje tereta. Model je nazvan pokretni oscilator sa klatnom jer sadrži osnovne karakteristike pomenutih dinamičkih pojmova. Ovaj model, prikazan na slici 4.10.c, ima dva dodatna stepena slobode, vertikalno pomeranje mase m2 u odnosu na položaj statičke ravnoteže, ovde označeno sa koordinatom y, i ugao klaćenja tereta, označenog kao ߮. Pomeranje mase m1 je definisano profilom brzina kretanja i pomeranjem strukture na mestu položaja kolica. Pokretni oscilator sa klatnom je originalan model kolica u problemima pokretnog opterećenja kod dizalica i ima cilj da opiše uticaj oscilacija ovešene mase u sistemu kolica i klaćenje tereta na oscilacije strukture. 4.3. MODEL KRETANJA KOLICA Ovde je usvojeno da se pokretni sistem‐kolica u toku jednog ciklusa kretanja duž staze kreće sa trapeznim profilom brzina. Ovaj profil se najčešće koristi prilikom proračuna mehanizama za kretanje kod dizalica, [1], slika 4.11.a,b. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 55 Slika 4.11. (a) Profil brzina kretanja, (b) Profil ubrzanja, (c) Profil kretanja Sistem kolica se kreće tako što ubrzava sa konstantnim ubrzanjem au u vremenu tu da bi dostigao nominalnu vrednost brzine kretanja vr, koju održava u vremenu ravnomernog kretanja tr i usporava usporenjem ak u vremenskom periodu tk do zaustavljanja. Na osnovu predstavljenog profila brzina izračunava se put koji pokretni sistem pređe duž staze u proizvoljnom vremenskom trenutku, slika 4.11.c. (a) (b) Slika 4.12. (a) Stepenasti profil, (b) Zarezani profil; [27] DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 56 Pored ovog profila brzina kretanja, postoje i stepenasti i zarezani profil brzine kretanja, slika 4.12, [27]. Trapezni profil je usvojen jer predstavlja oštriji kriterijum za analizu u odnosu na pomenute profile. Jednačine koje će biti prezentovane u sledećim poglavljima dozvoljaju da se koristi bilo koji profil kretanja. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 57 5 DINAMIČKA INTERAKCIJA NOSEĆE KONSTRUKCIJE I KOLICA Pojam dinamičke interakcije strukture i sistema koji se kreće po strukturi podrazumeva razmatranje prinudnih oscilacija strukture usled dejstva pokretnog opterećenja koje pored težine obuhvata i inercijalne efekte od pokretnog sistema. U klasi problema pokretnog opterećenja, model pokretne mase je prvi za koji se vezuje pojam dinamičke interakcije. Prikazom klasičnog problema pokretne sile i klasičnog problema pokretne mase, ovde će biti objašnjena formulacija ovog pojma. Na slici 5.1 prikazan je model proste grede po kojoj se konstantnom brzinom v kreće sila konstantnog intenziteta P. Sila P, u ovom slučaju, predstavlja model pokretnog sistema. Slika 5.1. Prosta greda opterećena vertikalnom silom konstantnog intenziteta P koja se kreće konstantnom brzinom v; Ojler‐Bernulijeva greda, dužine L, površine preseka A, modula elastičnosti E, momenta inercije I, gustine materijala ߩ DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 58 Diferencijalna jednačina transverzalnih oscilacija Ojler‐Bernulijeve grede, w=w(x,t), za 0 ≤ x ≤ L i 0 ≤ t ≤ L/v, ima sledeći oblik ߩܣ డమ௪డ௧మ ൅ ܧܫ డర௪ డ௫ర ൌ ܲߜሺݔ െ ݒݐሻ (5.1) gde se ߜሺݔ െ ߟሻ, Dirak‐ova funkcija (specijalni oblik jediničnog impulsa, [53]), koristi za analitičko rešavanje problema pokretnog opterećenja. Rešenje jednačine (5.1) je u analitičkom obliku dato u [8,12], pri čemu se pretpostavlja da je greda mirovala u početnom trenutku. Rešenje je veoma važno sa aspekta postavke dodatnog parametara kojim brzina kretanja utiče na oscilacije proste grede i koji iznosi ߙ ൌ భ்ଶ ఛ (5.2) gde su T1 osnovni period oscilovanja grede, a ߬ ൌ ܮ/ݒ ukupno vreme koje je potrebno za prelazak opterećenja od levog do desnog oslonca grede. Na slici 5.2 prikazan je isti model grede kao i u prethodnom slučaju, ali je pokretni sistem modeliran kao pokretna masa koja se kreće duž grede konstantnom brzinom v. Diferencijalna jednačina transverzalnih oscilacija Ojler‐Bernulijeve grede, w= w(x,t), usled dejstva pokretne mase, pri čemu je 0 ≤ x ≤ L i 0 ≤ t ≤ L/v, se postavlja u obliku ߩܣ డమ௪డ௧మ ൅ ܧܫ డర௪ డ௫ర ൌ ሾ݉௦௦݃ െ ݉௦௦ ௗమ௪ሺ௩௧,௧ሻ ௗ௧మ ሿߜሺݔ െ ݒݐሻ (5.3) Jednačina (5.3) jasno oslikava dinamičku interakciju između grede i mase jer izraz koji predstavlja kontaktnu‐prinudnu silu koja deluje na gredu ݉௦௦݃ െ ݉௦௦ ௗ మ௪ሺ௩௧,௧ሻ ௗ௧మ ൌ ܲ (5.4) zavisi od težine mase i od inercijalne sile od mase koja nastaje kao posledica transverzalnih oscilacija same grede. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 59 Slika 5.2. (a) Prosta greda opterećena pokretnom masom mss koja se kreće konstantnom brzinom v; Ojler‐Bernulijeva greda, dužine L, površine preseka A, modula elastičnosti E, momenta inercije I, gustine ߩ; (b) protivdejstvo strukture na pokretnu masu, (c) kretanje mase po deformisanoj strukturi Izraz (5.4) predstavlja i jednačinu kretanja pokretne mase u vertikalnom pravcu, pri čemu u ovom slučaju prema trećem Njutnovom zakonu sila P predstavlja uticaj, odnosno dejstvo grede na pokretnu masu (Slika 5.2b) i može se postaviti pomoću teoreme o promeni količine kretanja u vertikalnom pravcu. Član, ݉௦௦ ௗ మ௪ሺ௩௧,௧ሻ ௗ௧మ , koji je dodatno postavljen u (5.3), u odnosu na (5.1), uključuje inerciju same mase u analizu, čime se računa i drugi izvod po vremenu transverzalih oscilacija grede u tački na kojoj se nalazi pokretna masa. Kako se transverzalne oscilacije prizmatične grede određuju u funkciji koordinate x i funkciji od vremena, odnosno w = w(x,t), prvi izvod postaje ௗ௪ሺ௫,௧ሻௗ௧ ൌ డ௪ డ௫ డ௫ డ௧ ൅ డ௪ డ௧ = డ௪ డ௫ ݔሶ ൅ డ௪ డ௧ (5.5) pa drugi izvod po vremenu dobija sledeći oblik ௗమ௪ሺ௫,௧ሻ ௗమ௧ ൌ డమ௪ డ௫మ డ௫ డ௧ ݔሶ ൅ డ௪ డ௫ ݔሷ ൅ 2 డమ௪ డ௫డ௧ డ௫ డ௧ ൅ డమ௪ డ௧మ ൌ ݓ ′′ݔሶ ଶ ൅ ݓ ′ݔሷ ൅ 2ݓ ′ሶ ݔሶ ൅ ݓሷ (5.6) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 60 gde su oznake za diferenciranje po x i za izvod po vremenu date saglasno uobičajenoj praksi u dinamici. Ovde, veličine označene sa ݔሶ i ݔሷ predstavljaju brzinu i ubrzanje pokretne mase duž grede. U slučaju problema definisanog slikom 5.2, imamo da je ݔሶ ൌ ݒ ൌ ܿ݋݊ݏݐ. i ݔሷ ൌ 0, pa je često u radovim kontaktna sila koja deluje na gredu (5.4) prikazana u sledećem obliku ܲ ൌ ݉௦௦݃ െ݉௦௦ሺݓ ′′ݒଶ ൅ 2ݓ ′ሶ ݒ ൅ ݓሷ ሻ ∥௫ୀ௩௧ (5.7) Ovo odražava fizički smisao dinamičke interakcije strukture i pokretne mase jer se masa, u stvari kreće po deformisanoj strukturi, a samim tim javljaju se dodatni članovi u izrazu (5.7) koji se mogu identifikovati, [24] kao centripetalna sila ܨ௖ ൌ ݉௦௦ሺݓ ′′ݒଶሻ ∥௫ୀ௩௧ (5.8) i Koriolisova sila ܨ௖௢௥ ൌ ݉௦௦ሺ2ݓ ′ሶ ݒሻ ∥௫ୀ௩௧ (5.9) Rešenje j‐ne (5.3), prema (5.7) je u drugoj polovini dvadesetog veka dobijeno različitim analitičko‐numeričkim metodama koje su u sebi sadržale glomazan i komplikovan matematički aparat. Osnovne poteškoće dolaze od formulacije prinudne sile koja je, čak i u najjednostavnijoj formulaciji (5.1), funkcija od položaja, tj. od vremena, a dodatno je njen intenzitet funkcija od vremena kada su uključeni i uticaji inercije pokretnog sistema (5.3). I pored toga, korišćenje ovog pristupa ima sledeće nedostatke:  može se primeniti na modele greda i konzola, sa jednim ili više raspona  podužne deformacije grede su zanemarene  uticaj ubrzanja tereta je zanemaren, a njegova implementacija generiše još glomazniji matematički aparat  razmatranje ostalih modela pokretnog opterećenja je moguće samo u aproksimativnoj formi DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 61 Zbog svega navedenog, i u skladu sa ciljem istraživanja ove disertacije, u nastavku će biti izvršena postavka formulacije dinamičke interakcije noseće konstrukcije i kolica portalne dizalice prema kombinovanom pristupu, kao preduslov za definisanje jednačina kretanja sistema. 5.1 KONAČNOELEMENTNO MODELIRANJE POKRETNOG OPTEREĆENJA STRUKTURE Ovde se razmatra najopštiji slučaj dejstva pokretnog opterećenja na strukturu portalne dizalice i saglasno postavkama modela za strukturu, kolica i profil kretanja kolica (poglavlje 4), interaktivni uticaj pokretnog sistema na strukturu na mestu kontakta se predstavlja preko jedno vertikalne sile intenziteta Py(t), usmerene nadole, i jedne horizontalne sile intenziteta Px(t), usmerene udesno. Slika 5.3. Interaktivni uticaj pokretnog sistema na strukturu na mestu kontakta modela strukture i kolica; Py(t) ‐ vertikalna sila, Px(t) ‐ horizontalna sila Sile Py(t) i Px(t) deluju na mestu kontakta modela strukture i pokretnog sistema, odnosno deluju na elemente glavnog nosača dizalice, tj. stazu. Sa aspekta postavljenog KE modela strukture dizalice ovo znači da su samo elementi staze izloženi dejstvu spoljašnjih sila. Kako ove sile stalno menjaju svoj položaj na strukturi staze, primena standardnih softverskih paketa za MKE analizu nije DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 62 moguća čak i u slučaju da su sile Py(t) i Px(t) konstantnih intenziteta (p.a., sa aspekta poznavanja dometa komercijalnih softverskih paketa za statičku/dinamičku analizu konstrukcija, npr. SAP2000, KRASTA, Algor..). U tom smislu, potrebno je proširiti obim primene MKE, sa dodatnim analitičkim pristupom koji se naziva modeliranjem spoljašnjeg pokretnog opterećenja. Modeliranje spoljašnjeg pokretnog opterećenja se vrši formiranjem ekvivalentnih čvornih opterećenja koja su promenljiva u vremenu i zavise od modela staze i profila kretanja opterećenja. Ovde će biti prikazano formiranje ekvivalentnih čvornih opterećenja staze strukture portalne dizalice usled dejstva Py(t) i Px(t), na izdvojenom modelu staze, slika 5.4.a. Kao što je i uobičajeno u sličnim radovima iz ove oblasti, usvojene su sledeće pretpostavke:  pokretno opterećenje počinje svoje kretanje od krajnjeg levog čvora staze  položaj pokretnog sistema na stazi je određen koordinatom xm(t) koja odgovara usvojenom profilu brzina kretanja kolica  sve veličine se razmatraju u vremenskom domenu koji odgovara vremenu koje je potrebno pokretnom sistemu da pređe celu stazu Slika 5.4. (a) KE model staze izložen dejstvu interaktivnih sila, (b) Ekvivalentna čvorna opterećenja elementa s na kojem se nalazi pokretno opterećenje DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 63 Kada su elementi staze podvrgnuti dejstvu sila, sva čvorna opterećenja su po intenzitetu jednaka 0 osim čvornih opterećenja elementa s, na kome se nalazi pokretno opterećenje. Vektor čvornih opterećenja, za element s, se može predstaviti u obliku ሼ݂௦ሺݐሻሽ ൌ ሾ ଵ݂௦ሺݐሻ ଶ݂௦ሺݐሻ ଷ݂௦ሺݐሻ ସ݂௦ሺݐሻ ହ݂௦ሺݐሻ ଺݂௦ሺݐሻሿ் (5.10) Na osnovu [32], imamo sledeće relacije ଵ݂௦ሺݐሻ ൌ ଵܰሺߦሻ ௫ܲሺݐሻ (5.11a) ଶ݂௦ሺݐሻ ൌ െ ଶܰሺߦሻ ௬ܲሺݐሻ (5.11b) ଷ݂௦ሺݐሻ ൌ െ ଷܰሺߦሻ ௬ܲሺݐሻ (5.11c) ସ݂௦ሺݐሻ ൌ ସܰሺߦሻ ௫ܲሺݐሻ (5.11d) ହ݂௦ሺݐሻ ൌ െ ହܰሺߦሻ ௬ܲሺݐሻ (5.11e) ଺݂௦ሺݐሻ ൌ െ ଺ܰሺߦሻ ௬ܲሺݐሻ (5.11f) gde su ௜ܰ (i =1‐6), funkcije oblika za gredni linijski konačni element i definisane su izrazima (4.6), (4.7), (4.9), (4.10), (4.11) i (4.12), gde je iskorišćena sledeća smena ߦ ൌ ௫௟ (5.12) pri čemu je dužina segmenta staze ݈ ൌ ݈௡ (n =1‐10), a x rastojanje duž elementa do tačke dejstva sila. Dalje, potrebno je odrediti ekvivalentna čvorna opterećenja za sve čvorove KE modela strukture portalne dizalice, pri čemu su čvorna opterećenja za sve elemente, osim elemenata staze, po intenzitetu jednaka nuli. Razmatranje DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 64 prinudnih oscilacija sistema sa konačnim brojem stepena slobode, kao u ovom slučaju, se vrši preko jednačine ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൌ ۾ሺݐሻ (5.13) gde su ۻܛܜ, ۱ܛܜ, ۹ܛܜ, matrica inercije, matrica prigušenja i matrica krutosti strukture, respektivno; ܃ሷ , ܃ሶ , ܃, vektori ubrzanja, vektor brzina i vektor generalisanih koordinata strukture, respektivno; a ۾ሺݐሻ vektor spoljašnjih sila. Vektor spoljašnjih sila će biti postavljen u sledećem obliku ۾ሺݐሻ ൌ ሼܨଡ଼ଵ ܨଢ଼ଵ ܯଵ ܨଡ଼ଶ ܨଢ଼ଶ ܯଶ … ܨଡ଼ଵଵ ܨଢ଼ଵଵ ܯଵଵ 0…0ሽT (5.14) gde su ܨଡ଼୩, ܨଢ଼୩,ܯ୩, (k =1‐11), horizontalna sila, vertikalna sila i moment u čvoru k. Kako je poznata zakonitost promene položaja pokretnog opterećenja na gredi, u ovom slučaju (slika 5.4a), moguće je eksplicitno odrediti broj elementa s, na kome deluje pokretno opterećenje na sledeći način ݏ ൌ ܥ݁݋ ܾݎ݋݆ ቂ௫೘ሺ௧ሻ௟ ቃ ൅ 1 (5.15) Čvorovi s‐tog elementa su s i s+1. Odavde, korišćenjem izraza (5.11) možemo odrediti čvorna opterećenja elemenata, kada se pokretno opterećenje nalazi na elementu s, na sledeći način ܨଡ଼ୱሺݐሻ ൌ ଵܰሺߦሻ ௫ܲሺݐሻ (5.16a) ܨଢ଼ୱሺݐሻ ൌ െ ଶܰሺߦሻ ௬ܲሺݐሻ (5.16b) ܯୱሺݐሻ ൌ െ ଷܰሺߦሻ ௬ܲሺݐሻ (5.16c) ܨଡ଼ୱାଵሺݐሻ ൌ ସܰሺߦሻ ௫ܲሺݐሻ (5.16d) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 65 ܨଢ଼ୱାଵሺݐሻ ൌ െ ହܰሺߦሻ ௬ܲሺݐሻ (5.16e) ܯୱାଵሺݐሻ ൌ െ ଺ܰሺߦሻ ௬ܲሺݐሻ (5.16f) ܨଡ଼୨ሺݐሻ ൌ 0, j ൌ 1 െ 11 ⋀ j ് s, s ൅ 1 (5.16g) ܨଢ଼୨ሺݐሻ ൌ 0, j ൌ 1 െ 11 ⋀ j ് s, s ൅ 1 (5.16h) ܯ୨ሺݐሻ ൌ 0, j ൌ 1 െ 11 ⋀ j ് s, s ൅ 1 (5.16i) gde se u izraze za Ni (i=1‐6) postavlja zavisnost od globalne koordinate položaja xm(t) umesto od lokalne koordinate x, u sledećem obliku ߦ ൌ ௫೘ሺ௧ሻିሺ௦ିଵሻ௟௟ (5.17) Da bi se obuhvatio ceo vremenski domen, ukupno vreme koje je potrebno da bi opterećenje prešlo celu gredu, ߬, se mora diskretizovati na p koraka određenog vremenskog intervala ∆t, odnosno ߬ ൌ ݌ Δݐ (5.18) Odavde, sile i momenti sadrže vrednosti za vremenski trenutak ݐ ൌ ሺݎ െ 1ሻ Δݐ, odnosno za svaki vremenski korak r (r =1 ÷ p+1): ܨ௑୩ ൌ ൣܨ௑୩଴ ܨ௑୩Δ௧ ܨ௑୩ଶΔ௧ …ܨ௑୩௥Δ௧ …ܨ௑୩௣Δ௧൧ (5.19a) ܨ௒୩ ൌ ൣܨ௒୩଴ ܨ௒୩Δ௧ ܨ௒୩ଶΔ௧ …ܨ௒୩௥Δ௧ …ܨ௒୩௣Δ௧൧ (5.19b) ܯ୩ ൌ ൣܯ୩଴ ܯ୩Δ௧ ܯ୩ଶΔ௧ …ܯ୩௥Δ௧ …ܯ୩௣Δ௧൧ (5.19c) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 66 Za određivanje ekvivalentnih čvornih opterećenja strukture, usled dejstva pokretnog opterećenja, napravljen je originalni program prema algoritmu prikazanom na slici 5.5. 5.1.1 Prezentacija modeliranja spoljašnjeg opterećenja prema modelu pokretne sile Postavljeni algoritam za određivanje ekvivalentnih čvornih opterećenja je veoma važan za modeliranje problema pokretne sile pri čemu se pretpostavlja da je zakonitost promene intenziteta Px(t), Py(t) poznata. Ovo je posebno pogodno u situaciji kada se struktura modelira u određenom programskom paketu za MKE. Zbog boljeg razumevanja činjenice da se ekvivalentna čvorna opterećenja sračunavaju u zavisnosti od položaja i od vremena, ovde će biti date odgovarajuće ilustracije za opterećenja čvorova za model staze prikazan na slici 5.4, gde se opterećenje kreće konstantnom brzinom od 5 m/s na dužini od 40 m. Pri tome, biće razmatrana dva slučaja intenziteta pokretnog opterećenja: 1. Slučaj ‐ Px =1 N i Py =1 N. 2. Slučaj ‐ Px =1 cos(0.5 π t) N i Py =1 Sin(0.5 π t) N. Prvi slučaj, gde su pokretne sile konstantnog intenziteta, se može primeniti u slučajevima kada je masa od pokretnog opterećenja dovoljno manja od mase strukture pa se inercijalni efekti od mase mogu zanemariti. Na slici 5.6 prikazana su ekvivalentna čvorna opterećenja za ovaj slučaj. Na slici 5.7 prikazana su ekvivalentna čvorna opterećenja usled dejstva pokretnih sila prema 2. slučaju. Ovaj slučaj se može primeniti u situaciji da je pokretnim opterećenjem modelirana pokretna mašina koja u sebi ima, usled neuravnoteženosti sistema, ekcentrično postavljenu obrtnu masu. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 67 Slika 5.5. Algoritam za određivanje ekvivalentnog spoljašnjeg vektora strukture DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 68 a) b) (c) Slika 5.6. 3D prikaz ekvivalentnih čvornih opterećenja za čvorove k (k =1‐11), Lgn=40 m, v=5 m/s, Py=1 N=const., Px=1 N =const., p =200; (a) Horizontalna sila u čvorovima FXk, (b) Vertikalna sila u čvorovima FYk, (c) Moment u čvorovima Mk DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 69 a) b) c) Slika 5.7. 3D prikaz ekvivalentnih čvornih opterećenja za čvorove k (k =1‐11), Lgn=40 m, v=5 m/s, Py =sin(0.5 π t), Px = cos(0.5 π t), p =200; (a) Horizontalna sila u čvorovima FXk, (b) Vertikalna sila u čvorovima FYk, (c) Moment u čvorovima Mk DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 70 Izdvajanjem vrednosti za FX, FY i M za svaki čvor (k=1‐11) i njihovom implementacijom u (5.13) ili na modelu u određenom softveru (koji odgovara slici 5.4), moguće je dobiti odziv usled pokretne sile koji u potpunosti obuhvata promenljivost položaja opterećenja i dinamičke karakteristike prinude. 5.2 MATRIČNA PREZENTACIJA POMERANJA I UBRZANJA TAČKE KONTAKTA STRUKTURE I KOLICA Prinudne sile koje deluju na strukturu dizalice, ili na pokretni sistem, zavise od dinamičke interakcije ovih dvaju sistema i deluju u tački kontakta. Kako se kontakt ostvaruje duž cele staze potrebno je definisati pomeranje svake ove tačke i usaglasiti sa matematičkim modelom sistema pri čemu mora biti zadovoljena matrična postavka diferencijalnih jednačina kretanja. Ovde će biti razmotreno formiranje izraza za matričnu prezentaciju za podužno pomeranje wx i poprečno pomeranje wy ove tačke, kao i odgovarajuća ubrzanja. Na slici 5.8 prikazana su pomeranja tačke na stazi na kojoj se ostvaruje kontakt između strukture i pokretnog sistema. Analogno sa formulacijom iznetom u poglavlju 5.1, ovaj element će biti označen sa s. Slika 5.8. Pomeranje tačke kontakta na strukturi staze Vektor čvornih pomeranja elementa s, kao elementa KE modela staze, se može prikazati kao ܃ܛ ൌ ሾ ௑ܷ௦ ܷ௒௦ ܷఏ௦ ௑ܷ௦ାଵ ܷ௒௦ାଵ ܷఏ௦ାଵሿ் (5.20) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 71 Podužno pomeranje tačke u polju elementa s je zadato samom formulacijom KE koji je iskorišćen za modeliranje strukture, i u lokalnom koordinatnom sistemu, slika 4.7a, je definisano izrazom (4.5). Kako se lokalni i globalni koordinatni sistem elemenata KE modela strukture poklapaju za elemente staze, vektori čvornih pomeranja, izrazi (4.4) i (5.20), odgovaraju jedan drugom bez transformacije. Čvorna pomeranja strukture su generalisane koordinate sistema (i zavise od vremena), pa izraz (4.5) dobija sledeći oblik ݓ௫ሺݔ, ݐሻ ൌ ଵܰሺݔሻ ௑ܷ௦ ൅ ସܰሺݔሻ ௑ܷ௦ାଵ (5.21) gde su: x ‐ položaj tačke u polju na elementu s od levog čvora; N1, N4 funkcije oblika definisane izrazima (4.6) i (4.7); UXs, UXs+1 horizontalna pomeranja čvorova elementa s. Poprečno pomeranje tačke kontakta strukture i modela pokretnog opterećenja je definisano izrazom (4.8), i analogno prethodnom, može se predstaviti u sledećem obliku ݓ௬ሺݔ, ݐሻ ൌ ଶܰሺݔሻܷ௒௦ ൅ ଷܰሺݔሻܷఏ௦ ൅ ହܰሺݔሻܷ௒௦ାଵ ൅ ଺ܰሺݔሻܷఏ௦ାଵ (5.22) gde su: x ‐ položaj tačke u polju na elementu s, od levog čvora; N2, N3, N5, N6 funkcije oblika definisane izrazima (4.9), (4.10), (4.11) i (4.12); UYs, UYs+1 horizontalna pomeranja čvorova elementa s, Uθs, Uθs+1 rotacije čvorova elementa s Proširivanjem i usaglašavanjem sa brojem čvornih pomeranja elementa s, izraz (5.22) se može predstaviti kao proizvod dve matrice, odnosno ݓ௬ ൌ ۼܡ܃ܛ (5.23) gde je ۼܡ ൌ ሾ0 ଶܰ ଷܰ 0 ହܰ ଺ܰሿ (5.24) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 72 Konačnim proširivanjem i usaglašavanjem do ukupnog broja čvornih pomeranja strukture imamo da je ݓ௬ ൌ ۼ܇܃ (5.25) gde je ۼ܇ ൌ ሾ0 0 … ሾ0 ଶܰ ଷܰ 0 ହܰ ଺ܰሿ … 0 0 0ሿ (5.26) matrica vrsta koja u sebi sadrži podmatricu ۼܡ sa elementima koji poziciono odgovaraju čvornim pomeranjima elementa s na kome se nalazi pokretno opterećenje. Svi ostali elementi su jednaki 0. Elementi matrice ۼ܇ se sračunavaju na osnovu trenutnog položaja pokretnog opterećenja na elementu s, u zavisnosti od x. S obzirom da se staza ima 11 elemenata, podmatrica ۼܡ se šeta u okviru elemenata ۼ܇ zaključno sa 33 elementom matrice, a ostali elementi su uvek jednaki 0. Analogno prethodnom, izraz (5.21) se može napisati u obliku ݓ௫ ൌ ۼܠ܃ܛ (5.27) gde je ۼܠ ൌ ሾ ଵܰ 0 0 ସܰ 0 0ሿ (5.28) U svom konačnom obliku izraz (5.27) postaje ݓ௫ ൌ ۼ܆܃ (5.29) gde je ۼ܆ ൌ ሾ0 0 … ሾ ଵܰ 0 0 ସܰ 0 0ሿ … 0 0 0ሿ (5.30) Ovde je interesantno napomenuti da su na ovaj način, izrazima (5.25) i (5.29) aksijalno pomeranje i poprečno pomeranje tačke na stazi dobile karakter DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 73 partikularnog rešenja podužnih oscilacija i transverzalnih oscilacija prizmatičnog tela, [43], koji pretpostavlja rešenje u obliku proizvoda 2 funkcije, jedne koja zavisi od položaja i druge koja zavisi od vremena. Ubrzanje tačke, poprečno na element staze, iznosi ݓሷ ௬ ൌ ௗ మ௪೤ ௗ௧మ (5.31) i zamenom sledećih izraza koji se dobijaju iz (5.25), ݓ௬ᇱᇱ ൌ ۼܡᇱᇱ܃ (5.32a) ݓሶ ௬ᇱ ൌ ۼܡᇱ ܃ሶ (5.32b) ݓ௬ᇱ ൌ ۼܡᇱ ܃ (5.32c) ݓሷ ൌ ۼܡ܃ሷ (5.32d) u izraz (5.31), koji je u analogiji sa (5.6), dobija se ݓሷ ௬ ൌ ۼ܇ᇱᇱ܃ ݔሶ ଶ ൅ ۼ܇ᇱ ܃ ݔሷ ൅ 2ۼ܇ᇱ ܃ሶ ݔሶ ൅ ۼ܇܃ሷ (5.33) U izrazu (5.33), iskorišćene su sledeće oznake ݔሶ ൌ ݔሶ௠, ݔሷ ൌ ݔሷ௠, za brzinu i ubrzanje pokretnog opterećenja, respektivno, i ۼ܇ᇱ ൌ ሾ0 0 … ሾ0 ଶܰᇱ ଷܰᇱ 0 ହܰᇱ ଺ܰᇱሿ … 0 0 0ሿ (5.34) ۼ܇ᇱᇱ ൌ ሾ0 0 … ሾ0 ଶܰᇱᇱ ଷܰᇱᇱ 0 ହܰᇱᇱ ଺ܰᇱᇱሿ … 0 0 0ሿ (5.35) pri čemu se vrednosti za ௜ܰᇱ, ௜ܰᇱᇱ (i =2,3,5,6) računaju prema izrazima iz Tabele 4.1. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 74 Analogno prethodnom, ubrzanje tačke u podužnom pravcu dobija konačni oblik ݓሷ ௫ ൌ ۼ܆ᇱ ܃ ݔሷ ൅ 2ۼ܆ᇱ ܃ሶ ݔሶ ൅ ۼ܆܃ሷ (5.36) gde je ۼ܆ᇱ ൌ ሾ0 0 … ሾ ଵܰᇱ 0 0 ସܰᇱ 0 0ሿ … 0 0 0ሿ (5.37) pri čemu se vrednosti za ௜ܰᇱ (i =1,4) računaju prema izrazima iz Tabele 4.1. Izraz (5.36) ima jedan član manje u odnosu na izraz (5.33) zbog same postavke KE elementa staze, odnosno tabele 4.1, jer ௑ܰᇱᇱ postaje nula matrica. Na sledećoj slici prikazan je algoritam za određivanje matrica ۼ܆ i ۼ܇, kao i matrica ۼ܆ᇱ ۼ܇ᇱ , ۼ܇ᇱᇱ, ukoliko se iskoriste vrednosti u zagradama koje označavaju diferenciranje po x. Algoritam je izveden za izdvojeni model staze, slika 5.4. 5.3 ODREĐIVANJE INTERAKTIVNIH SILA ZA POSTAVLJENE MODELE KOLICA U poglavlju 5.1 dat je prikaz uticaja pokretnog sistema na strukturu portalne dizalice preko jedne horizontalne sile i jedne vertikalne sile, intenziteta Px(t) i Py(t) i smerova prikazanih na slici 5.3. Na usvojenom modelu strukture određena su ekvivalentna čvorna opterećenja, tj. izvršeno je formiranje vektora spoljašnjih sila usled dejstva pokretnog opterećenja u zavisnosti od Px(t) i Py(t). U nastavku će ove sile biti označene samo sa Px i Py. S obzirom na postavke date izrazom (5.16), može se primetiti da vektor spoljašnjeg opterećenja prema (5.14) ima sledeći oblik ۾ሺݐሻ ൌ ሼ0 0 0… ଵ݂ୱሺtሻ ଶ݂ୱሺtሻ ଷ݂ୱሺtሻ ସ݂ୱሺtሻ ହ݂ୱሺtሻ ଺݂ୱሺtሻ … 0 0 0ሽT (5.38) gde su ௜݂௦ (i =1‐6) definisani izrazima (5.11). Odnosno, vektor spoljašnjeg opterećenja se sastoji elemenata koji su jednaki nuli i 6 elemenata koji se računaju na osnovu (5.11), a poziciono odgovaraju čvornim pomeranjima elementa s na kome se nalazi pokretno opterećenje. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 75 Slika 5.9. Algoritam za određivanje matrica NY i NX DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 76 Zamenom izraza (5.11) u izraz (5.10) kojim je predstavljen vektor čvornih opterećenja elementa s, i korišćenjem matrica (5.24) i (5.27) imamo da je ሼ݂௦ሺݐሻሽ ൌ ۼܠ܂ ௫ܲ െ ۼܡ܂ ௬ܲ (5.39) Korišćenjem analogije za postavkama iznetim u poglavlju 5.2, a koje se tiču formiranja matrica vrsta ۼ܇ i ۼ܆ koje su date izrazima (5.26) i (5.30), vektor spoljašnjih sila koji deluje na strukturu dobija pogodan matrični oblik ۾ሺݐሻ ൌ ۼ܆܂ ௫ܲ െ ۼ܇܂ ௬ܲ (5.40) Primetno je da je moguće razmatrati prinudne oscilacije strukture, zamenom (5.40) u (5.13), superpozicijom partikularnih rešenja koja se dobijaju na osnovu zamene pojedinačnih sabiraka iz (5.40). Na osnovu trećeg Njutonovog zakona, (protiv)dejstvo koje struktura vrši na pokretni sistem je dato silama Px i Py, pri čemu one imaju suprotne smerove u odnosu na one prikazane na slici 5.3. Intenziteti sila Px i Py se upravo određuju daljom dinamičkom analizom modela sistema kolica. Ovde će biti postavljeni izrazi ovih sila za sve usvojene modele sistema kolica, tj. model pokretne mase, model pokretnog oscilatora i model oscilatora sa klatnom. 5.3.1 Model pokretne mase Na osnovu postavljenog modela pokretne mase, poglavlje 4.2.1, može se razmatrati kretanje mase mss kao kretanje materijalne tačke koje se sastoji od kretanja xm, i prenosnog kretanja definisanog pomeranjem tačke na strukturi staze na kojoj se nalazi pokretna masa, u poprečnom i podužnom pravcu, označenim sa wy i wx, respektivno, slika 5.10.a. Na materijalnu tačku deluju interaktivne sile Px i Py i težina mase Gss, slika 5.10.b. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 77 Na osnovu teoreme o promeni količine kretanja, i projektovanjem na ose koordinatnog sistema, dobija se sledeće ௗௗ௧ ൫݉௦௦ݓሶ ௬൯ ൌ െܩ௦௦ ൅ ௬ܲ (5.41) ௗௗ௧ ሺ݉௦௦ݓሶ ௫ ൅ ݉௦௦ݔሶ௠ሻ ൌ െ ௫ܲ (5.42) Slika. 5.10. Dinamički model pokretne mase Odnosno, dobija se ௬ܲ ൌ ݉௦௦݃ ൅ ݉௦௦ݓሷ ௬ (5.43) ௫ܲ ൌ െ݉௦௦ݔሷ௠ െ ݉௦௦ݓሷ ௫ (5.44) 5.3.2 Model pokretnog oscilatora Na osnovu postavljenog modela pokretnog oscilatora, poglavlje 4.2.2, može se razmatrati kretanje materijalnog sistema koji se sastoji od dve materijalne tačke masa m1 i m23. Kretanje materijalne tačke m1 se sastoji od kretanja xm, i prenosnog kretanja definisanog pomeranjem tačke na strukturi staze na kojoj se nalazi pokretna masa, u poprečnom i podužnom pravcu, označenim sa wy i wx. Masa m23 se u vertikalnom pravcu kreće prema usvojenoj koordinati y, a u horizontalnom DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 78 pravcu isto kao masa m1, slika 5.11.a. Na materijalni sistem deluju sile Px i Py, i težine G1 i G23, slika 5.11.b. Na osnovu teoreme o promeni količine kretanja, i projektovanjem na ose koordinatnog sistema, dobija se ௗௗ௧ ൫݉ଶଷݕሶ െ ݉ଵݓሶ ௬൯ ൌ ܩଵ ൅ ܩଶଷ െ ௬ܲ (5.45) ௗௗ௧ ሺ݉ଵሺݓሶ ௫ ൅ ݔሶ௠ሻ ൅ ݉ଶଷሺݓሶ ௫ ൅ ݔሶ௠ሻሻ ൌ െ ௫ܲ (5.46) Slika. 5.11. Dinamički model pokretnog oscilatora Odnosno, koristeći da je ukupna masa pokretnog sistema mss, imamo da je ௬ܲ ൌ ݉௦௦݃ ൅ ݉ଵݓሷ ௬ െ ݉ଶଷݕሷ (5.47) ௫ܲ ൌ െ݉௦௦ݔሷ௠ െ ݉௦௦ݓሷ ௫ (5.48) 5.3.3 Model pokretnog oscilatora sa klatnom Na osnovu postavljenog modela pokretnog oscilatora sa klatnom, poglavlje 4.2.3, može se razmatrati kretanje materijalnog sistema koji se sastoji od 3 materijalne tačke mase m1, m2 i m3. Kretanje materijalne tačke m1 se sastoji od DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 79 kretanja xm, i prenosnog kretanja definisanog pomeranjem tačke na strukturi staze na kojoj se nalazi pokretna masa, u poprečnom i podužnom pravcu, označenim sa wy i wx. Masa m2 se u vertikalnom pravcu kreće prema usvojenoj koordinati y, a u horizontalnom pravcu isto kao masa m1. Kretanje mase m3 se sastoji od prenosnog kretanja koje joj saopštava materijalna tačka m2 i od relativnog kretanja koje se može poistovetiti sa matematičkim klatnom i definisano je uglom otklona ߮, slika 5.12.a. Na materijalni sistem deluju sile Px i Py, težine G1, G2 i G3, slika 5.12.b. Slika. 5.12. Model pokretnog oscilatora sa klatnom Na osnovu teoreme o promeni količine kretanja, i projektovanjem na ose koordinatnog sistema, dobija se ௗௗ௧ ሺെ݉ଵݓሶ ௬ ൅ ݉ଶݕሶ ൅ ݉ଷሺݕሶ െ ܮ௨ ሶ߮ sin ߮ሻሻ ൌ ܩଵ ൅ ܩଶ ൅ ܩଷ െ ௬ܲ (5.49) ௗௗ௧ ሺሺ݉ଵ ൅ ݉ଶ ൅݉ଷሻሺݔሶ௠ ൅ ݓሶ ௫ሻ ൅ ݉ଷܮ௨ ሶ߮ cos ߮ሻ ൌ െ ௫ܲ (5.50) Odnosno, koristeći (4.26) dobija se DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 80 ௬ܲ ൌ ݉௦௦݃ ൅ ݉ଵݓሷ ௬ െ ሺ݉ଶ ൅݉ଷሻݕሷ ൅ ݉ଷܮ௨ሺ ሷ߮ sin ߮ ൅ ሺ ሶ߮ ሻଶ cos ߮ሻ (5.51) ௫ܲ ൌ െ݉௦௦ሺݔሷ௠ ൅ ݓሷ ௫ሻ െ ݉ଷܮ௨ሺ ሷ߮ cos߮ െ ሺ ሶ߮ ሻଶ sin߮ሻ (5.52) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 81 6 MATEMATIČKI MODELI PORTALNE DIZALICE Postavljanje matematičkog modela sistema portalne dizalice, sa aspekta naslovnog problema, se vrši kombinovanjem konačnoelementnog modela strukture i dinamičkog modela kolica preko dinamičke interakcije ovih dvaju sistema koja je definisana u prethodnom poglavlju. U opštem obliku, diferencijalna jednačina kretanja sistema sa konačnim brojem stepeni slobode se postavlja kao ۻ ܙሷ ൅ ۱ ܙሶ ൅ ۹ ܙ ൌ ۴ሺݐሻ (6.1) pri čemu su ۻ, ۱, ۹ ‐ matrica inercije, matrica prigušenja i matrica krutosti celog sistema, respektivno; ܙሷ , ܙሶ , ܙ ‐ vektor generalisanih ubrzanja, vektor generalisanih brzina, vektor generalisanih koordinata celog sistema, respektivno; F(t) ‐ vektor spoljašnjeg opterećenja celog sistema. S obzirom da su postavljena tri različita modela kolica, ovde će biti izvršeno formiranje tri različita matematička modela celog sistema. Na slikama 6.1, 6.2 i 6.4 je prikazan model B, kao opštiji KE model strukture portalne dizalice, ali postavljene jednačine podrazumevaju korišćenje oba modela. 6.1 MATEMATIČKI MODEL PORTALNE DIZALICE PREMA MODELU POKRETNE MASE Portalna dizalica prema modelu pokretne mase podrazumeva model prikazan na slici 6.1, a u sebi sadrži sve postavke modela prema poglavljima 4.1 i 4.2.1. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 82 Slika 6.1 Portalna dizalica prema modelu pokretne mase Kako je već predočeno u poglavlju 4.2.1, model pokretne mase ne uključuje dodatne stepene slobode već samo one definisane samom strukturom dizalice, pa je vektor generalisanih koordinata sistema ܙ ൌ ܃ (6.2) Samim tim, jednačina prinudnih oscilacija strukture dizalice je data kao ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൌ ۾ሺݐሻ (6.3) Zamenom (5.40) u (6.3) dobija se sledeće ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൌ ۼ܆܂ ௫ܲ െ ۼ܇܂ ௬ܲ (6.4) Zamenom (5.43, 5.44) u (6.3) dobija se ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൌ െۼ܇܂൫݉௦௦݃ ൅ ݉௦௦ݓሷ ௬൯ ൅ ۼ܆܂ሺെ݉௦௦ݔሷ௠ െ ݉௦௦ݓሷ ௫ሻ (6.5) Zamenom izraza (5.33, 5.36) u (6.5) dobija se ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൅ ݉ୱୱۼ܇܂൫ۼ܇܃ሷ ൅ 2ۼ܇ᇱ ܃ሶ ݔሶ ൅ ۼ܇ᇱᇱ܃ݔሶ ଶ ൅ ۼ܇ᇱ ܃ݔሷ൯ ൅ ൅݉௦௦ۼ܆܂൫ۼ܆ᇱ ܃ ݔሷ ൅ 2ۼ܆ᇱ ܃ሶ ݔሶ ൅ ۼ܆܃ሷ ൯ ൌ െۼ܇܂݉௦௦݃ െ ۼ܆܂݉௦௦ݔሷ௠ (6.6) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 83 U konačnom obliku, prethodni izraz se može predstaviti kao ሺۻܛܜ ൅ ۻ૚ሻ܃ሷ ൅ ሺ۱ܛܜ ൅ ۱૚ሻ܃ሶ ൅ ሺ۹ܛܜ ൅ ۹૚ሻ܃ ൌ െۼ܇܂݉௦௦݃ െ ۼ܆܂݉௦௦ݔሷ௠ (6.7) gde su ۻ૚ ൌ ۼ܇܂ۼ܇݉௦௦൅ۼ܆܂ۼ܆݉௦௦ (6.7a) ۱૚ ൌ 2 ۼ܇܂ۼ܇ᇱ ݉௦௦ ݔሶ ൅ 2 ۼ܆܂ۼ܆ᇱ ݉௦௦ ݔሶ (6.7.b) ۹૚ ൌ ۼ܇܂ۼ܇ᇱᇱ ݉௦௦ ݔሶ ଶ ൅ ۼ܇܂ۼ܇ᇱ ݉௦௦ ݔሷ ൅ ۼ܆܂ۼ܆ᇱ ݉௦௦ ݔሷ (6.7.c) 6.2 MATEMATIČKI MODEL PORTALNE DIZALICE PREMA MODELU POKRETNOG OSCILATORA Portalna dizalica prema modelu pokretnog oscilatora podrazumeva model prikazan na slici 6.2, a u sebi sadrži sve postavke modela prema poglavljima 4.1 i 4.2.2. Slika 6.2. Portalna dizalice prema modelu pokretnog oscilatora DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 84 Kako je već predočeno u poglavlju 4.2.2, model pokretnog oscilatora uključuje 1 dodatni stepen slobode, pa je vektor generalisanih koordinata sistema ܙ ൌ ൜܃ݕൠ (6.8) Diferencijalna jednačina prinudnih oscilacija strukture ima oblik ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൌ ۼ܆܂ ௫ܲ െ ۼ܇܂ ௬ܲ (6.9) Zamenom (5.47) i (5.48) u (6.13) dobija se sledeće ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൌ െۼ܇܂൫݉௦௦݃ ൅ ݉ଵݓሷ ௬ െ ݉ଶଷݕሷ൯ ൅ ۼ܆܂ሺെ݉௦௦ݔሷ௠െ݉௦௦ݓሷ ௫ሻ (6.10) Zamenom (5.33) i (5.36) u (6.14) dobija se ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൅ ݉ଵۼ܇܂൫ۼ܇܃ሷ ൅ 2ۼ܇ᇱ ܃ሶ ݔሶ ൅ ۼ܇ᇱᇱ܃ݔሶ ଶ ൅ ۼ܇ᇱ ܃ݔሷ൯ െ െۼ܇܂݉ଶଷݕሷ ൅ ݉ୱୱۼ܆܂൫ۼ܆܃ሷ ൅ 2ۼ܆ᇱ ܃ሶ ݔሶ ൅ ۼ܆ᇱ ܃ݔሷ൯ ൌ െۼ܇܂݉௦௦݃ െ ۼ܆܂݉௦௦ݔሷ௠ (6.11) Diferencijalna jednačina kretanja mase m23, slika 6.3, dobija sledeći oblik ݉ଶଷݕሷ ൌ ܩଶଷ െ ܨ௢௣ (6.12) pri čemu je intenzitet sile u opruzi ܨ௢௣ ൌ ݇ሺݕ ൅ ݓ௬ ൅ ଵ݂௦௧ሻ (6.13) Kako je izduženje u položaju statičke ravnoteže ଵ݂௦௧ ൌ ீమయ௞ (6.14) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 85 diferencijalna jednačina kretanja (6.12) dobija oblik ݉ଶଷ ݕሷ ൅ ݇൫ݕ ൅ ݓ௬൯ ൌ 0 (6.15) Slika 6.3. Vezani materijalni sistem modela pokretnog oscilatora Zamenom (5.25) u (6.15) dobija se sledeće ݉ ଶଷ ݕሷ ൅ ݇ݕ ൅ ݇ ۼ܇܃ ൌ 0 (6.16) Kombinovanjem izraza (6.11) i (6.16), u oblik prema (6.1) dobija se konačan matematički model oscilacija pokretnog oscilatora u obliku ൤ۻܛܜ ൅ ۻ૛ െ݉ଶଷۼ܇܂0 ݉ଶଷ ൨ ൜ ܃ሷ ݕሷ ൠ ൅ ቂ ۱ܛܜ ൅ ۱૛ 00 0ቃ ൜ ܃ሶ ݕሶ ൠ ൅ ൤ ۹ܛܜ ൅ ۹૛ 0 ݇ۼ܇ ݇൨ ൜ ܃ ݕൠ ൌ ൌ ൜െۼ܇܂݉௦௦݃ െ ۼ܆܂݉௦௦ݔሷ௠0 ൠ (6.17) gde su ۻ૛ ൌ ۼ܇܂ۼ܇ ݉ଵ൅ۼ܆܂ۼ܆ ݉௦௦ (6.17a) ۱૛ ൌ 2 ۼ܇܂ۼ܇ᇱ ݉ଵ ݔሶ ൅ 2 ۼ܆܂ۼ܆ᇱ ݉௦௦ ݔሶ (6.17b) ۹૛ ൌ ݉ଵۼ܇܂ۼ܇ᇱᇱݔሶ ଶ ൅ ݉ଵۼ܇܂ۼ܇ᇱ ݔሷ ൅ ۼ܆܂ۼ܆ᇱ ݉௦௦ ݔሷ (6.17c) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 86 6.3 MATEMATIČKI MODEL PORTALNE DIZALICE PREMA MODELU POKRETNOG OSCILATORA SA KLATNOM Portalna dizalica prema modelu pokretnog oscilatora sa klatnom podrazumeva model prikazan na slici 6.4, a u sebi sadrži sve postavke modela prema poglavljima 4.1 i 4.2.3. Slika 6.4. Portalna dizalica prema modelu pokretnog oscilatora sa klatnom Na osnovu pretpostavki iz poglavlja 4.2.3 vektor generalisanih koordinata sistema ima sledeći oblik ܙ ൌ ൝ ܃ ݕ ߮ ൡ (6.18) Diferencijalna jednačina oscilovanja za strukturu portalne dizalice ima isti oblik kao i u prethodnim slučajevima, tj. ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൌ ۼ܆܂ ௫ܲ െ ۼ܇܂ ௬ܲ (6.19) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 87 Zamenom (5.51) i (5.52) u (6.19), i zamenom (5.33) i (5.36) u (6.19) dobija se sledeće ۻܛܜ܃ሷ ൅ ۱ܛܜ܃ሶ ൅ ۹ܛܜ܃ ൅ ݉௦௦ۼ܆܂ሺۼ܆܃ሷ ൅ 2ۼ܆ᇱ ܃ሶ ݔሶ ൅ ۼ܆ᇱ ܃ݔሷሻ ൅ ݉ଵۼ܇܂൫ۼ܇܃ሷ ൅ 2ۼ܇ᇱ ܃ሶ ݔሶ ൅ ۼ܇ᇱᇱ܃ݔሶ ଶ ൅ ۼ܇ᇱ ܃ݔሷሻ െ ۼ܇܂ሺ݉ଶ ൅ ݉ଷሻݕሷ ൅ ۼ܆܂݉ଷܮ௨ሺ ሷ߮ cos ߮ െ ሺ ሶ߮ ሻଶ sin߮ሻ ൅ ൅ۼ܇܂݉ଷܮ௨ሺ ሷ߮ sin߮ ൅ ሺ ሶ߮ ሻଶ cos߮ሻ ൌ െۼ܆܂݉௦௦ ݔሷ௠ െ ۼ܇܂݉௦௦݃ (6.20) Jednačine kretanja elemenata pokretnog sistema se postavljaju na osnovu ekvivalentnog dinamičkog modela prikazanog na slici 6.5. Slika 6.5. Dinamički model pokretnog oscilatora sa klatnom Na masu m2 deluju sopstvena težina G2, sila u opruzi Fop i reakcija užeta S. Diferencijalna jednačina kretanja mase m2 je oblika ݉ଶݕሷ ൌ ܩଶ െ ܨ௢௣ ൅ ܵ cos߮ (6.21) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 88 Intenzitet sile u opruzi iznosi ܨ௢௣ ൌ ݇ሺݕ ൅ ݓ௬ ൅ ଶ݂௦௧ሻ (6.22) gde je f2st izduženje opruge u položaju statičke ravnoteže . Dalje, posmatra se kretanje mase m3 izdvojene iz sistema (slika 6.5b). Primenom jednačina relativnog kretanja materijalne tačke, diferencijalna jednačina kretanja mase m3 je definisana sa ݉ଷ Ԧܽ௥ ൌ ܨԦ௣௑௜௡ ൅ ܨԦ௣௒௜௡ ൅ ܩଷሬሬሬሬԦ ൅ Ԧܵ (6.23) gde su intenziteti inercijalnih sila ܨ௣௑௜௡ ൌ ݉ଷሺݔሷ௠ ൅ ݓሷ ௫ሻ (6.24) ܨ௣௒௜௡ ൌ ݉ଷݕሷ (6.25) Projektovanjem jednačine (6.23) na ose prirodnog triedra dobija se ݉ଷܮ௨ ሷ߮ ൌ ݉ଷݕሷ ݏ݅݊ ߮ െ ܩଷ ݏ݅݊ ߮ െ݉ଷሺݔሷ௠ ൅ ݓሷ ௫ሻ ܿ݋ݏ ߮ (6.26) ݉ଷܮ௨ሺ ሶ߮ ሻଶ ൌ ݉ଷݕሷ cos ߮ െ ܩଷ cos߮ ൅݉ଷሺݔሷ௠ ൅ ݓሷ ௫ሻ sin߮ ൅ ܵ (6.27) Odavde, preko izraza (6.26) i (6.27), može se prikazati sila u užetu u sledećem obliku ܵ cos ߮ ൌ െ݉ଷݕሷ ൅ ܩଷ ൅݉ଷܮ௨ሺ ሷ߮ sin߮ ൅ ሺ ሶ߮ ሻଶ cos ߮ሻ (6.28) Kako su ܩଶ ൌ ݉ଶ݃, ܩଷ ൌ ݉ଷ݃ i statičko izduženje opruge DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 89 ଶ݂௦௧ ൌ ሺ௠మା௠యሻ௚௞ (6.29) diferencijalna jednačina (6.21) dobija sledeći oblik ሺ݉ଶ ൅ ݉ଷሻݕሷ ൌ െ݇ሺݕ ൅ ݓ௬ሻ ൅ ݉ଷܮ௨ሺ ሷ߮ sin ߮ ൅ ሺ ሶ߮ ሻଶ cos ߮ሻ (6.30) Zamenom (5.25) u (6.30) dobija se ሺ݉ଶ ൅ ݉ଷሻݕሷ ൅ ݇ݕ ൅ ݇ۼ܇܃ െ݉ଷܮ௨ሺ ሷ߮ sin ߮ ൅ ሺ ሶ߮ ሻଶ cos ߮ሻ ൌ 0 (6.31) Izraz (6.31) postaje ܮ௨ ሷ߮ ൌ ݕሷ sin߮ െ ݃ sin߮ െ ݔሷ௠ cosφ െ ݓሷ ௫ cos ߮ (6.32) Zamenom (5.29) u (6.37) dobija se ۼ܆܃ሷ cos ߮ ൅ 2ۼ܆ᇱ ܃ሶ ݔሶ cos ߮ ൅ ۼ܆ᇱ ܃ݔሷ cos߮ ൅ ܮ௨ ሷ߮ െ ݕሷ sin ߮ ൅ ݃ sin߮ ൌ െݔሷ௠ cos߮ (6.33) Kombinovanjem izraza (6.20), (6.31) i (6.33), u oblik prema (6.1) dobija se konačan matematički model sistema u obliku ቎ ۻܛܜ ൅ ۻ૜ 0 ݉ଷܮ௨ۼ܆܂ cosφ 0 ݉ଶ ൅݉ଷ െ݉ଷܮ௨ sin߮ۼ܆ cosφ െ sin߮ ܮ௨ ቏ ቐ ܃ሷ ݕሷ ሷ߮ ቑ ൅ ൅ ቎ ۱ܛܜ ൅ ۱૜ 0 െ݉ଷܮ௨ۼ܆܂φሶ sinφ0 0 െ݉ଷܮ௨ ሶ߮ cos ߮ 2ݔሶۼ܆ᇱ cosφ 0 0 ቏ ቐ ܃ሶ ݕሶ ሶ߮ ቑ ൅ ൅ ቎ ۹ܛܜ ൅ ۹૜ ൅݇ۼ܇܂ 0 ൅݇ۼ܇ ݇ 0 ݔሷۼ܆ᇱ cosφ 0 0 ቏ ൝ ܃ ݕ ߮ ൡ ൌ ቐ െۼ܇܂݉௦௦݃ െ ۼ܆܂݉௦௦ ݔሷ௠ 0 െݔሷ௠ cosφ െ݃ sinφ ቑ (6.34) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 90 gde su ۻ૜ ൌ ۼ܇܂ۼ܇ ݉ଵ൅ۼ܆܂ۼ܆ ݉௦௦ (6.34a) ۱૜ ൌ 2 ۼ܇܂ۼ܇ᇱ ݉ଵ ݔሶ ൅ 2 ۼ܆܂ۼ܆ᇱ ݉௦௦ ݔሶ (6.34b) ۹૜ ൌ ݉ଵۼ܇܂ۼ܇ᇱᇱݔሶ ଶ ൅ ݉ଵۼ܇܂ۼ܇ᇱ ݔሷ ൅ ۼ܆܂ۼ܆ᇱ ݉௦௦ ݔሷ ൅ ݇ۼ܇܂ۼ܇ (6.34c) 6.4 POSTUPAK REŠAVANJA DIFERENCIJALNIH JEDNAČINA Na osnovu matematičkih modela koji su prikazani u prethodnim poglavljima, mogu se odrediti generalisana pomeranja, brzine i ubrzanja sistema. Zbog složenosti rešavanja postavljenih diferencijalnih jednačina drugog reda sa promenljivim koeficijentima, ovde je prikazan postupak za određivanje rešenja prema izabranom modelu portalne dizalice. Generalno, potrebno je: 1. Odrediti matricu krutosti i matricu inercije KE modela izabrane strukture, Kst i Mst (poglavlje 4.1.3 i 4.1.4) 2. Odrediti sopstvene frekvencije izabranog modela strukture pomoću izraza (4.25), u određenom programu 3. Izdvojiti prve dve najniže kružne frekvencije ω1 i ω2, pretpostaviti koeficijent prigušenja ξ i odrediti matricu prigušenja modela strukture, Cst, preko (4.23) 4. Podeliti ukupni vremenski domen τ na p koraka tako da se dobije što manji vremenski interval ∆t 5. Odrediti vrednosti za ݔሷ௠, ݔሶ௠, ݔ௠ i sve komponente matrica Nଡ଼, Nଡ଼ᇱ , kao i matrica Nଢ଼, Nଢ଼ᇱ , Nଢ଼ᇱᇱ za svaki vremenski trenutak t=(r‐1) ∆t, gde je r =1 ÷ p+1 6. Rešiti jednačine (6.7) ili (6.17) ili (6.34) pomoću Njumarkove metode pri čemu je ovo potrebno uraditi za svaki vremenski korak r, gde je r =1 ÷ p+1 DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 91 7 IDENTIFIKACIJA I ANALIZA ODZIVA MODELA Određivanje odziva postavljenih modela portalne dizalice je izvršeno u originalnom programu MovMass (u kodu Wolfram Mathematica 6, [64]) kojim su obuhvaćene sve postavke date u poglavljima 4, 5 i 6. U poglavlju 4.1 prikazani su opšti modeli A i B kojim su predstavljena dva osnovna konstrukciona tipa strukture portalne dizalice. Dodatno, ovde su tim modelima pridružene glavne mere portalne dizalice, statičke karakteristike elemenata strukture i nazivna nosivost (mss). U tabelama 7.1 i 7.2 prikazani su podaci na osnovu kojih će biti izvršena analiza odziva matematičkih modela portalne dizalice (poglavlje 6). Usvojeni materijal strukture (čelik) ima gustinu ߩ = 7850 kg/m3, a modul elastičnosti iznosi E = 2,1 1011 N/m2. Tabela 7.1. Karakteristike usvojenog modela portalne dizalice (A) Model A n 1‐10 11 12 13,14 mss = 60 t An [m 2] 0,090 0,085 0,070 0,048 In [m4] 0,041 0,036 0,024 0,010 DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 92 Tabela 7.2. Karakteristike usvojenog modela portalne dizalice (B) Model B n 1‐10 11 12 13,14 mss = 40 t An [m 2] 0,0584 0,055 0,051 0,0432 In [m4] 0,0255 0,024 0,0154 0,01 Vremenski interval algoritma za određivanje rešenja iznosi ∆t=0,01 s, osim u slučajevima gde je eksplicitno navedeno drugačije. Inicijalni modeli uključuju prigušenje u strukturi dizalice sa koeficijentim prigušenja koji iznosi ߦଵ ൌ ߦଶ ൌ 0,005, koji podrazumeva prigušenje materijala i aerodinamičko prigušenje, [63]. U svim proračunima je usvojeno ubrzanje zemljine teže g = 9,81 m/s2. Na slici 7.1 prikazani su dva izabrana profila brzina kretanja pokretnog sistema, vm1 i vm2, koji su formirani saglasno dužinama staza na modelima A i B i prema sledećim karakteristikama. Profil vm1 sa nominalnom brzinom od 3 m/s i ubrzanjem/usporenjem od 0,6 m/s2 predstavlja aktuelni domet brzina kretanja kolica (www.liebherr.de, Tabela 1.2), sa nominalnim vremenima ubrzanja i usporenja od 5 s [3]. Profil vm2 sa nominalnom brzinom od 5 m/s i ubrzanjem/usporenjem od 1,25 m/s2 pretpostavlja ekstremne performanse koje kolica portalnih dizalica mogu imati u skoroj budućnosti, iako je brzina od 4 m/s već ostvarena kod kolica kontejnerskih dizalica (www.kockskrane.de), a ubrzanje od 1,2 m/s2 (iako naglašeno kao ektremnom vrednošću) iskorišćeno u radu [27]. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 93 A) B) Slika 7.1. Profili brzina kretanja kolica vm1 i vm2, A) model A, B) model B 7.1 ODZIV MODELA POKRETNE MASE 7.1.A. ANALIZA ODZIVA ‐ MODEL A Prvo je izvršena analiza odziva portalne dizalice prema modelu A, sa karakteristikama datim u tabeli 7.1, rešavanjem j‐ne (6.7). Određene su sve sopstvene frekvencije modela strukture, a u tabeli 7.3 dati su prvih 7 odgovarajućih perioda oscilovanja. Na slici 7.2 prikazana su prva četiri oblika oscilovanja, sa odgovarajućim frekvencijama. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 94 Verifikacija postavljenog modela strukture izvršena je tako što je formiran linijski KE model u programskom paketu SAP 2000 [65], prema podacima iz tabele 7.1, gde su nakon izvršene modalne analize dobijene sopstvene frekvencije strukture. Odgovarajući periodi oscilovanja su prikazani u tabeli 7.3. Upoređivanjem vrednosti dobijenih preko matematičkog modela korišćenjem konačnoelementnog pristupa (program MovMass) i rezultata dobijenih čistim numeričkim putem preko MKE softvera (SAP 2000), može se uočiti poklapanje rezultata sa relativnim odstupanjem do 1,15%. Tabela 7.3. Uporedni prikaz rezultata prvih 7 perioda oscilovanja i 1 2 3 4 5 6 7 MovMass, Ti [s] 0.5788 0.2025 0.0677 0.0430 0.0348 0.0286 0.019 SAP2000, Ti [s] 0.5794 0.2026 0.0678 0.0433 0.0352 0.0288 0.019 1. mod; ߱ଵ= 10,438 s‐1, ଵܶ = 0,579 s, ଵ݂ = 1,72 Hz 2. mod; ߱ଶ= 30,998 s‐1, ଶܶ = 0,202 s, ଶ݂ = 4,93 Hz 3. mod; ߱ଷ= 92,77 s‐1, ଷܶ = 0,0677 s, ଷ݂ = 14,76 Hz 4. mod; ߱ସ= 146,1 s‐1, ସܶ = 0,043 s, ସ݂ = 23,725 Hz Slika 7.2. Prva 4 oblika oscilovanja strukture dizalice DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 95 Na slici 7.3 je dat 3D prikaz vertikalnih pomeranja čvorova staze strukture u vremenu (dijagram je okrenut zbog bolje prezentacije) pri čemu se masa mss=60 t kreće profilom brzina vm1 (slika 7.1.A). Slika 7.3. Vertikalna pomeranja čvorova staze, UYi (i=1‐11); vm1 Evidentno je da najveća vertikalna pomeranja ima čvor 6 (sredina rasponog dela nosača) i da ona nastaju kada se pokretna masa nalazi na sredini raspona glavnog nosača. Na slici 7.4 prikazana su pomeranja čvorova staze u horizontalnom pravcu. Prikaz je dat na 2D dijagramu jer su pomeranja svih čvorova na stazi gotovo identična što se objašnjava velikom aksijalnom krutošću elemenata staze, tj. neznatnim aksijalnim deformacijama samih elemenata. Slika 7.4. Horizontalna pomeranja čvorova staze, UXi (i=1‐11) U i=1-11Xi 0 5 10 15 -0.02 0.00 0.02 0.04 0.06 Vreme t s Ho riz on ta ln a po m er an ja cv or ov as ta ze m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 96 Na osnovu prethodnog, može se izdvojiti čvor 6 kao karakteristično mesto na strukturi preko čijih odziva će biti izvršena analiza rezultata ovog modela. Izdvojena pomeranja ovog čvora, u vertikalnom i horizontalnom pravcu, data su na slici 7.5. Slika 7.5. Vertikalno i horizontalno pomeranje čvora 6, UY6 i UX6; vm1 Pomeranja u horizontalnom pravcu čvora 12 (na sredini krute noge) i čvora 14 (na sredini elastične noge) su prikazana na slici 7.6. Vertikalna pomeranja ovih čvorova nisu prikazana jer imaju veoma zanemarljive vrednosti zbog velike aksijalne krutosti elemenata. Središnje pomeranje krute noge, UX12, ima niže vrednosti od pomeranja UX14, zbog veće savojne krutosti krute noge u odnosu na elastičnu nogu. Slika 7.6. Horizontalno pomeranje čvora 12 i čvora 14, UX12 i UX14 U U X6 Y6 0 5 10 15 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 Vreme t s Po m er an ja -c vo r6 m  U U X14 X12 0 5 10 15 -0.02 -0.01 0.00 0.01 0.02 0.03 0.04 Vreme t s Po m er an ja -c vo r1 2 ic vo r1 4 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 97 7.1.A.1 Uticaj brzine i ubrzanja kolica Uticaj brzine i ubrzanja pokretne mase se sagledava kroz razmatranje odziva modela usled implementacije profila brzine kretanja vm2 (slika 7.1.A) u proračun. Na slici 7.7. dato je vertikalno pomeranje čvora 6, a na slici 7.8 je dato horizontalno pomeranje čvora 6, usled kretanja mase mss=60 t prema postavljenim profilima brzina kretanja sa slike 7.1.A, u funkciji odnosa vremena i ukupnog vremena koje je potrebno da masa pređe ceo nosač (t/τ). Slika 7.7. Ukupno vertikalno pomeranja čvora 6, UY6 ; vm1 i vm2 Slika 7.8. Horizontalno pomeranja čvora 6, UX6 ; vm1 i vm2 vm1 vm2 0.0 0.2 0.4 0.6 0.8 1.0 -0.05 -0.04 -0.03 -0.02 -0.01 0.00 Bezdimenziono vreme tt Ve rti ka ln op om er an je -c vo r6 m  vm1 vm2 Usporenje 0.0 0.2 0.4 0.6 0.8 1.0 -0.04 -0.02 0.00 0.02 0.04 0.06 Bezdimenziono vreme tt Ho riz on ta ln o po m er an je -c vo r6 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 98 Amplitude oscilovanja su veće pri većoj brzini pokretne mase (vm2), za oba pomeranja čvora 6. Povećanje amplituda za UY6 je najveće u periodu usporenja što predstavlja uticaj horizontalne sile od ubrzanja mase. Takođe, najveće amplitude za UX6 nastaju u periodu usporenja, gde su dostignute najveće vrednosti pomeranja čvora 6 u horizontalnom pravcu. Najveće vrednosti ovih pomeranja su date u tabeli 7.4. Analiza dobijenih rezultata prikazanih na slikama 7.7 i 7.8 će biti izvršena uporedno sa vrednostima horizontalnog i vertikalnog pomeranja čvora 6 koji se dobijaju statičkim proračunom za sledeća dva slučaja opterećenja: 1‐ vertikalna sila intenziteta mss g (=588,6 kN) deluje na sredini raspona glavnog nosača (na čvoru 6); 2‐ horizontalna sila intenziteta ±mss a2 (=±75 kN), gde je a2=1,25 m/s2, deluje u čvoru 1. Rezultati su dati u tabeli 7.4. Tabela 7.4. Prikaz rezultata za pomeranja čvora 6 Statički proračun Pokretna masa 1: mss g 2: ± mss a2 vm1 vm2 UX6 2,22 cm ± 1,8 cm Maks. UX6 4 cm 5,44 cm UY6 ‐5,04 cm ± 0,28 cm Maks. UY6 ‐5,25 cm ‐5,41 cm Najveće pomeranje UY6 za vm2 iznosi 5,41 cm i veće je od vrednosti dobijene za vm1 (5,25 cm). Vrednost od 5,41 cm predstavlja povećanje od 7,3 % u odnosu na vrednost od 5,04 cm koja se dobija statičkim proračunom za slučaj 1. S obzirom da najveća vrednost vertikalnog središnjeg pomeranja nastaje kada se teret nalazi kod sredine raspona u periodu ravnomernog kretanja kolica, ova vrednost ne predstavlja bitno povećanje sa inženjerske tačke gledišta dokaza krutosti. Najveće horizontalno pomeranje čvora 6 za vm2 nastaje u periodu kočenja (usporenja) kolica i iznosi 5,44 cm, a za vm1 najveća vrednost za UX6 od 4 cm nastaje kada se kolica kreću ravnomerno. Vrednost od 5,44 cm, u ovom slučaju, predstavlja bitno povećanje i dostiže red veličina najvećeg vertikalnog pomeranja čvora 6. Takođe, poređenjem sa vrednošću od 1,8 cm dobijenom statičkim proračunom, slučaj 2 (Tabela 7.4), vrednost od 5,44 cm predstavlja bitno povećanje. Potrebno je DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 99 naglasiti da vrednost horizontalnog pomeranja zavisi i od uticaja težine tereta, što je obuhvaćeno dinamičkih modelom pokretne mase. Ovim se ostvaruje bolja simulacija ponašanja strukture u nestacionarnim periodima kretanja kolica, što je ovde prezentovano uvođenjem ubrzanja/usporenja od 1,25 m/s2, iako ova vrednost predstavlja ekstremnu performansu kolica portalnih dizalica. Sa ovako određenim pomeranjima, brzinama i ubrzanjima čvorova strukture, moguće je odrediti promenu intenziteta interaktivnih sila koja deluju na strukturu staze. Razmatra se profil brzina kretanja vm2. Na slici 7.9 je prikazan odnost intenziteta vertikalne sile Py i težine mase (=588,6 kN). Slika 7.9. Odnos vertikalne sile Py i težine mss g; m=mss = 60 t; vm2 Slika 7.10. Horizontalna sila, Px; mss = 60 t; vm2 0.0 0.2 0.4 0.6 0.8 1.0 0.97 0.98 0.99 1.00 1.01 1.02 1.03 Bezdimenziono vreme tt Py m g P - m xx 0 2 4 6 8 10 12 -100 -50 0 50 100 150 Vreme t s Ho riz on ta ln as ila Px kN  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 100 Najveće povećanje vertikalne sile Py koja deluje na strukturu iznosi 3,0 % u odnosu na težinu mase. Na slici 7.10 prikazana je ukupna horizontalna sila Px koja deluje na strukturu u ovom slučaju. Ona osciluje oko vrednosti koja sadrži samo uticaj ubrzanja kolica (koje može imati vrednosti do ±75 kN), a najveće amplitude nastaju u periodu usporenja. Odzivi strukture usled kretanja kolica konstantnom brzinom Odziv strukture pri kretanju kolica konstantnom brzinom se određuju rešavanjem j‐ne (6.7) pri čemu se u algoritmu koristi da je ubrzanje jednako nuli, tj. ݔሷ௠ ൌ ݔሷ ൌ 0, a položaj mase na strukturi je određen sa ݔ௠ሺݐሻ ൌ ݒ ݐ. Vremenski korak, za ovaj slučaj, iznosi ∆t=0,0032 s. Za vrednost brzine od ݒ ൌ 5 ݉/ݏ ൌ ܿ݋݊ݏݐ., koja predstavlja nominalnu brzinu kolica u profilu brzina kretanja vm2, dobijena su pomeranja čvora 6 u horizontalom i vertiklanom pravcu, slika 7.11. Slika 7.11. Pomeranja čvora 6, UY6 , UX6 ; vm=5 m/s=const. Najveća vrednost za UY6 iznosi 5,05 cm. Poređenjem ove vrednosti sa pomeranjem čvora 6 dobijenim statičkim proračunom za slučaj 1 (5,04 cm, tabela 7.4) može se konstatovati zanemarljiva razlika. UY6 UX6 0 2 4 6 8 -0.04 -0.02 0.00 0.02 Vreme t s Po m er an ja -c vo r6 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 101 Najveće horizontalno pomeranje čvora 6 od 3,41 cm (slika 7.11) nastaje 2,6 s po početku kretanja kolica, odnosno kada se kolica nalaze na 4 elementu KE modela staze. Statičkim proračunom kada vertikalna sila mss g (=588,6 kN) deluje u čvoru 4, dobija se horizontalno pomeranje čvora 6 od 3 cm. Dakle, najveće pomeranje UX6 dobijeno preko matematičkog modela ima povećanje od 13,6% u odnosu na vrednost dobijenu statičkim proračunom. 7.1.A.2 Uticaj prigušenja u strukturi Pored osnovnog modela koji uključuje prigušenje u strukturi dizalice sa koeficijentim prigušenja ߦଵ ൌ ߦଶ ൌ 0,005, ovde se dodatno razmatra slučaj sa mogućim koeficijentom prigušenja u zavarenoj strukturi ߦଵ ൌ ߦଶ ൌ 0,02 i maksimalno mogućim prigušenjem u strukturi ߦଵ ൌ ߦଶ ൌ 0,07, koji se koriste u dinamičkim analizama primenom MKE, [66]. Pokretna masa od 60 t se kreće profilom brzina vm2. Vremenski interval algoritma za određivanje rešenja iznosi ∆t=0,0048 s. Na slici 7.12 prikazana su vertikalna pomeranja čvora 6, a na slici 7.13 horizontalna pomeranja čvora 6 sa detaljem vrednosti pomeranja u periodu usporenja kolica. Slika 7.12. Vertikalno pomeranja čvora 6, UY6 ; ߦ = 0,005; 0,02 i 0,07 0 2 4 6 8 10 12 -0.05 -0.04 -0.03 -0.02 -0.01 0.00 Vreme t s Ve rti ka ln op om er an je -c vo r6 m  -0.055 -0.050 -0.045 = 7 % =2 % =0,5 %    DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 102 Slika 7.13. Horizontalno pomeranja čvora 6, UX6 ; ߦ = 0,005, 0,02 i 0,07 U tabeli 7.5 prikazana su maksimalna vertikalna pomeranja i maksimalna horizontalna pomeranja u periodu ubrzanja, za čvor 6, sa procentualnim smanjenjem ovih vrednosti u odnosu na inicijalni model sa prigušenjem ߦ ൌ 0,5%. Tabela 7.5. Prikaz rezultata za pomeranje čvora 6 ߦ ൌ 0,005 ߦ ൌ 0,02 Smanj. ߦ ൌ 0,07 Smanj. Maks. UX6 5,44 cm 4,9 cm 9,9 % 3,9 cm 28,3 % Maks. UY6 5,41 cm 5,32cm 1,67 % 5,16 cm 4,8 % Smanjenje maksimalnog vertikalnog pomeranja za ߦ ൌ 2% iznosi 1,67 % što predstavlja zanemarljiv podatak i tek sa prigušenjem od ߦ ൌ 7% dostiže vrednost od 4,8 %. Horizontalno pomeranje čvora 6, slika 7.13, prikazuje da se amplitude oscilovanja smanjuju sa povećanjem prigušenja u strukturi. Primetno je značajno smanjenje maksimalne vrednosti pomeranja za koeficijent prigušenja ߦ ൌ 2%, od 9,9 %. Daljim povećanjem prigušenja do maksimalne vrednosti od ߦ ൌ 7% dolazi se do smanjenja od 28,3 % u odnosu na neprigušene oscilacije (ߦ ൌ 0,5%). Rezultati sa slike 7.13 su bitni jer pokazuju da je moguće, u određenoj meri, uticati na smanjenje amplituda pomeranja usled povećanja ubrzanja kolica (poglavlje 7.1.A.1) povećanjem prigušenja u strukturi modifikacijom materijala ili veza između strukturalnih elemenata. Međutim, potrebno je ovde napomenuti da je tehnologija i način povećanja prigušenja u strukturi nepoznat (sa aspekta 0 2 4 6 8 10 12 -0.02 0.00 0.02 0.04 Vreme t s Ho riz on ta ln o po m er an je -c vo r6 m  0.030 0.035 0.040 0.045 0.050 0.055 0.060 = 0,5 % = 2 % = 7 %   DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 103 autora) u oblasti projektovanja dizalica. Generalno, razmatranje oscilacija bez prigušenja dovodi do sigurnijih konstrukcija u smislu dokaza napona i deformacija, jer dinamički odzivi (pomeranja i opterećenja elemenata) imaju veće amplitude oscilovanja oko statičkih vrednosti. Međutim, savremena istraživanja koja se bave pouzdanošću mašinskih delova i struktura se upravo baziraju na smanjenju područja sigurnosti, [67], u smislu projektovanja tzv. racionalnih delova i struktura. 7.1.A.3 Uticaj inercijalnih efekata pokretne mase Ovde se razmatraju odzivi strukture dobijeni preko modela pokretne mase i modela pokretne sile. Masa od 60 t se kreće brzinom vm2. Odzivi modela pokretne sile su dobijeni u istom programu MovMass, koji je postavljen za model pokretne mase, pri čemu su zanemareni članovi koji predstavljaju inerciju masu u izrazu (6.7). Na slici 7.14 uporedno su prikazana vertikalna pomeranja središnjeg čvora za pomenute modele. Slika 7.14. Vertikalno pomeranje čvora 6, UY6 ; ݉ =60 t; vm2 Maksimalna vrednost pomeranja UY6 prema modelu pokretne sile u ovom slučaju iznosi 5,2 cm i niža je od maksimalne vrednosti dobijene za model pokretne mase (5,41 cm). Pokretna sila Pokretna masa 0 2 4 6 8 10 12 -0.05 -0.04 -0.03 -0.02 -0.01 0.00 Vreme t s Ve rti ka ln op om er an je -c vo r6 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 104 Slika 7.15. Ukupno horizontalno pomeranja čvora 6, UX6 ; ݉ =60 t; vm2 Na slici 7.18 prikazana su horizontalna pomeranja čvora 6 prema modelu pokretne mase i pokretne sile. Maksimalna vrednost pomeranja u periodu usporenja, za UX6, iznosi 4,65 cm za model pokretne sile. Ova vrednost je niža od odgovarajuće dobijene preko modela pokretne masa koja iznosi 5,44 cm. Oblik grafika pokazuje da model pokretne sile predstavlja dobru aproksimacuju posmatranog modela portalne dizalice. Ovo je posebno važno sa praktičnog aspekta jer model pokretne mase daje tačnije rezultate ali je obim programiranja i složenost rešavanja izuzetno veliki. Sa aspekta posmatranja uticaja kolica kroz jednu masu, ovaj model prevazilazi okvire potrebnog matematičkog aparata za opisivanje dinamičkog ponašanja portalnih dizalica pri normalim radnim uslovima. Osnovna prednost modela pokretne mase je dobijanje spektra frekvencija celog sistema, dok su u modelu pokretne sile uključene samo sopstvene frekvencije strukture. Ovo može imati odgovarajuću primenu u analizama dinamike dizalice usled slučajne pobude. Pokretna sila Pokretna masa 0 2 4 6 8 10 12 -0.02 0.00 0.02 0.04 Vreme t s Ho riz on ta ln op om er an je -c vo r6 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 105 7.1.B. ANALIZA ODZIVA ‐ MODEL B Ovde će biti izvršena analiza odziva portalne dizalice prema modelu pokretne mase na modelu B, tabela 7.2, rešavanjem j‐ne (6.7). Određene su sve sopstvene frekvencije strukture pri čemu su periodi oscilovanja prvih sedam frekvencija dati u tabeli 7.6. Na sledećoj slici prikazana su prva 2 oblika oscilovanja sa odgovarajućim frekvencijama. Tabela 7.6. Prikaz rezultata prvih 7 perioda oscilovanja i 1 2 3 4 5 6 7 Ti [s] 0.8950 0.2173 0.0953 0.0739 0.0556 0.0506 0.0311 1. mod; ߱ଵ= 7,02 s‐1, ଵܶ = 0,895 s, ଵ݂ = 1,117 Hz 2. mod; ߱ଶ= 28,91 s‐1, ଶܶ = 0,217 s, ଶ݂ = 4,6 Hz Slika 7.16. Prva 2 oblika oscilovanja strukture dizalice Slika 7.17. Vertikalna pomeranja čvorova staze, UYi (i=1‐11) DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 106 Na slici 7.17 je dat 3D prikaz vertikalnih pomeranja čvorova staze strukture (dijagram je okrenut zbog bolje prezentacije) pri čemu se pokretna masa mss=40 t kreće profilom brzina vm1, slika 7.1.B. Na slici 7.18 prikazana su pomeranja čvorova staze u horizontalnom pravcu. Slika 7.18. Horizontalna pomeranja čvorova staze, UXi (i=1‐11); vm1 Na ovom modelu, a saglasno graficima sa slike 7.17 i slike 7.18 mogu se izdvojiti čvor 1 (slobodni kraj prepusta dizalice) i čvor 7 (sredina rasponog dela dizalice) kao karakteristična mesta na strukturi dizalice za analizu odziva jer se ekstremi vertikalnog pomeranja javljaju za ta dva čvora, a horizontalna pomeranja su ista za sve čvorove staze strukture zbog velike aksijalne krutosti elemenata. Na sledećoj slici prikazana su vertikalna pomeranja čvorova 1 i 7. Slika 7.19. Vertikalna pomeranja čvora 1, UY1 i čvora 7, UY7 ; vm1 U i=1-11Xi 0 5 10 15 20 -0.10 -0.05 0.00 0.05 Vreme t s Ho riz on ta ln ap om er an ja cv or ov as ta ze m  U U Y1 Y7 0 5 10 15 20 -0.08 -0.06 -0.04 -0.02 0.00 0.02 0.04 Vreme t s Ve rti ka ln ap om er an ja -c vo r1 ic vo r7 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 107 Najveća vertikalna pomeranja čvora 1 nastaju kada se teret nalazi na početku svog kretanja, tj. na slobodnom kraju prepusta, a najveće vertikalno pomeranje čvora 7 nastaje kada se teret nalazi na sredini rasponog dela nosača. Na sledećoj slici prikazana su samo pomeranja u horizontalnom pravcu čvora 12 i čvora 14 jer vertikalna pomeranja imaju niske vrednosti (do 0,045 cm) zbog velike aksijalne krutosti elemenata. Slika 7.20 Horizontalno pomeranje čvora 12 i čvora 14, UX12 i UX14 7.1.B.1 Uticaj brzine i ubrzanja kolica Ovde su dobijeni odzivi strukture, pored korišćenja profila brzina vm1, implemeniranjem profila brzina kretanja kolica vm2, slika 7.1.B. Slika 7.21 Vertikalno pomeranje čvora 1; UY1; vm1 i vm2 U U X14 X12 0 5 10 15 20 -0.08 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 Vreme t s Po m er an ja -c vo r1 2 ic vo r1 4 m  v v m1 m2 0.0 0.2 0.4 0.6 0.8 1.0 -0.08 -0.06 -0.04 -0.02 0.00 0.02 0.04 Bezdimenziono vreme tt Ve rti ka ln op om er an je -c vo r1 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 108 Za pomenute profile brzina kretanja na slici 7.21 je dato vertikalno pomeranje čvora 1, na slici 7.22 vertikalno pomeranje čvora 7 i na slici 7.23 horizontalno pomeranje čvora 7. Slika 7.22 Vertikalno pomeranje čvora 7; UY7; vm1 i vm2 Slika 7.23 Horizontalno pomeranje čvora 7; UX7; vm1 i vm2 Analiza rezultata sa slika 7.21, 7.21 i 7.23 će biti izvršena uporedno sa rezultatima dobijenim statičkim proračunom za tri slučaja opterećenja: 1‐ vertikalna sila intenziteta 392,4 kN (=mss g) deluje na slobodnom kraju prepusta (na čvoru 1); 2 ‐ vertikalna sila intenziteta mss g deluje na sredini rasponog dela glavnog nosača (na čvoru 7); 3‐ horizontalna sila intenziteta ± 50 kN (= ±mss a2, pri čemu je a2=1,25 m/s2) deluje u čvoru 1. v v m1 m2 0.0 0.2 0.4 0.6 0.8 1.0 -0.04 -0.02 0.00 0.02 Bezdimenziono vreme tt Ve rti ka ln op om er an je -c vo r7 m  Ubrzanje Usporenje v v m1 m2 0.0 0.2 0.4 0.6 0.8 1.0 -0.20 -0.15 -0.10 -0.05 0.00 0.05 0.10 Bezdimenziono vreme tt Ho riz on ta ln op om er an je -c vo r7 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 109 Tabela 7.7. Prikaz rezultata za pomeranja čvora 1 i čvora 7 Statički proračun Pokretna masa 1: mss g 2: mss g 3: ± mss a2 vm1 vm2 UX1 ‐7 cm +2,15 cm ± 3,42 cm Maks. UX1 ‐10,3 cm ‐13,78 cm UY1 ‐7,8 cm +2,75 cm ± 0,9 cm Maks. UY1 ‐8,6 cm ‐9,47 cm UY7 +2,75 cm ‐5,27 cm ∓ 0,27 cm Maks. UY7 ‐5,53 cm ‐5,38 cm Maksimalno pomeranje staze u horizontalnom pravcu iznosi 13,78 cm i nastaje pri kretanju mase profilom brzina vm2. S obzirom da se ova vrednost javlja pri početku kretanja može se izvršiti njeno poređenje sa rezultatima statičkog proračuna prema slučajevima 1 i 3. Kako je sumarna vrednost prema ovim slučajevima opterećenja 10,42 cm (tabela 7.7), vrednost od 13,78 cm predstavlja povećanje od 32,2 %. Takođe, najveća vrednost za UX1 pri kretanju kolica profilom vm2 (13,78 cm) je veća za 33,7 % od vrednosti 10,3 cm koja se dobija pri vm1. Grafici vertikalnog pomeranja čvora 1 (UY1, slika 7.21) pokazuju mala odstupanja u amplitudama u periodima ubrzanja i usporenja, a najveća vrednost ovog pomeranja od 9,47 cm, pri kretanju profilom brzina vm2, predstavlja povećanje od 10,1 % u odnosu na odgovarajuću vrednost (8,6 cm) pri vm1. Najveće vertikalno pomeranje čvora 1 (9,47 cm) pri vm2 predstavlja povećanje od 8,85 % u odnosu na vrednost sumarnog vertikalnog pomeranja pri statičkom proračunu prema slučajevima 1 i 3, a koja iznosi od 8,7 cm. Pomeranja sredine rasponog dela glavnog nosača prema profilima brzina kretanja vm1 i vm2, UY7 (slika 7.22), pokazuju veoma male razlike. Najveća vrednost, pri vm1, iznosi 5,53 cm i predstavlja povećanje od 4,9 % u odnosu na vrednost dobijenu statičkim proračunom (5,27 cm) prema slučaju opterećenja 2. Odzivi strukture usled kretanja kolica konstantnom brzinom Rešavanjem j‐ne (6.7) pri čemu se koristi da je ubrzanje jednako nuli, tj. ݔሷ௠ ൌ ݔሷ ൌ 0, a položajem mase na strukturi određenim sa ݔ௠ሺݐሻ ൌ ݒ ݐ, dobija se odziv strukture pri kretanju kolica konstantnom brzinom. Vremenski korak, za ovaj slučaj, iznosi ∆t=0,005 s. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 110 Za vrednost brzine od ݒ ൌ 5 ݉/ݏ ൌ ܿ݋݊ݏݐ., prikazani su horizontalno pomeranje čvora 1, vertikalno pomeranje čvora 1 i vertikalno pomeranje čvora 7 na slikama 7.24, 7.25 i 7.26, respektivno. Slika 7.24. Horizontalno pomeranje čvora 1‐UX1, v=5 m/s=const. Slika 7.25. Vertikalno pomeranje čvora 1‐UY1, v=5 m/s=const. Slika 7.26. Vertikalno pomeranje čvora 7‐UY7, v=5 m/s=const. 0 2 4 6 8 10 -0.06 -0.04 -0.02 0.00 0.02 0.04 Vreme t s Ho riz on ta ln op om er an je cv or -1 m  0 2 4 6 8 10 -0.08 -0.06 -0.04 -0.02 0.00 0.02 Vreme t s Ve rti ka ln op om er an je cv or -1 m  0 2 4 6 8 10 -0.04 -0.02 0.00 0.02 Vreme t s Ve rti ka ln op om er an je cv or -7 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 111 Najveće pomeranje za UX1 na početku kretanja iznosi 7 cm i odgovara vrednosti dobijenoj statičkim proračunom za slučaj 1 (tabela 7.7), što se može zaključiti i za najveće pomeranje UY1 koje u ovom slučaju iznosi 7,8 cm. Pomeranje UY7 dostiže svoj maksimum od 5,4 cm kada se kolica nalaze na sredini rasponog dela nosača što predstavlja povećanje od 2,4 % u odnosu na vrednost dobijenu statičkim proračunom za slučaj opterećenja 2 (tabela 7.7). 7.1.B.2 Uticaj prigušenja u strukturi Na slici 7.27 prikazano je horizontalno pomeranje čvora 1‐UX1, pri kretanju kolica profilom brzina kretanja vm2, a prigušenje pored inicijalnog (0,5 %) iznosi 7 %. Pri ovom prigušenju, najveća vrednost ovog pomeranja iznosi 13,36 cm što predstavlja smanjenje od 3,1 % u odnosu na vrednost za prigušenje od 0,5 % (tabela 7.6). Slika 7.27. Horizontalno pomeranje čvora 1, vm2; ߦ = 0,005 i ߦ = 0,07 7.2 ODZIV MODELA POKRETNOG OSCILATORA Odzivi modela portalne dizalice prema modelu pokretnog oscilatora se dobijaju rešavanjem j‐ne (6.17). 0 2 4 6 8 10 12 14 -0.10 -0.05 0.00 0.05 Vreme t s Ho riz on ta ln o po m er an je -c vo r1 m  -0.14 -0.12 -0.10 -0.08 -0.06 =0,5 % = 7%  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 112 7.2.A. ANALIZA ODZIVA ‐ MODEL A Prvo se razmatra odziv modela A, tabela 7.1. Ukupna masa mss=60 t je podeljena na masu kolica m1=3 t i masu vitla i tereta m23=57 t. Inicijalna krutost opruge u modelu iznosi k=109 N/m. Na slici 7.28 su prikazana pomeranja čvora 6 pri profilu brzina kretanja vm1, a na slici 7.29 za vm2. Slika 7.28. Pomeranja čvora 6, vm1; k=109 N/m Slika 7.29. Pomeranja čvora 6, vm2; k=109 N/m U tabeli 7.8 su prikazani rezultati za najveća pomeranja čvora 6, u horizontalnom i vertikalnom pravcu. U U X6 Y6 0 5 10 15 -0.04 -0.02 0.00 0.02 0.04 Vreme t s Po m er an ja -c vo r6 m  UX6 UY6 0 2 4 6 8 10 12 -0.04 -0.02 0.00 0.02 0.04 Vreme t s Po m er an ja -c vo r6 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 113 Tabela 7.8. Prikaz rezultata za pomeranja čvora 6 Pokretni oscilator vm1 vm2 Maks. UX6 4 cm 5,44 cm Maks. UY6 ‐5,25 cm ‐5,41 cm Poređenjem sa vrednostima iz tabele 7.4 može se primetiti poklapanje rezultata. Kako i grafici sa slika 7.28 i 7.29 odgovaraju graficima na slikama 7.7 i 7.8 može se zaključiti, sa aspekta odziva strukture, da nema razlike između modela pokretnog oscilatora i modela pokretne mase. 7.2.A.1. Uticaj krutosti opruge Određivanjem svih pomeranja čvorova strukture, moguće je odrediti vertikalno pomeranje duž staze u funkciji vremena, tj. odrediti wy preko j‐ne (5.25). Na taj način dobija se vertikalno pomeranje tačke na strukturi koja se nalazi iznad tereta m23 koji se pomera u vertikalnom pravcu sa koordinatom y. Na slici 7.30 su prikazani wy i y, kada se razmatra profil brzina kretanja vm1, a koeficijent krutosti opruge iznosi k=107 N/m. Primetna su veoma mala odstupanja ovih dvaju pomeranja. Slika 7.30. Pomeranje y i wy; vm2; k=107 N/m y wy 0 2 4 6 8 10 12 -0.05 -0.04 -0.03 -0.02 -0.01 0.00 Vreme t s wy ,y m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 114 Smanjivanjem vrednosti krutosti opruge sa k=107 N/m na vrednost od k=105 N/m, dobija se uporedni prikaz pomeranja y (slika 7.31). Primetno je odstupanje oblika ovih grafika. Slika 7.31. Pomeranje y; vm2; k=105, k=107 N/m 7.2.B ANALIZA ODZIVA ‐ MODEL B Ovde se razmatra model B, tabela 7.2. Masa kolica iznosi m1=2 t, a masa vitla i tereta m23 = 38 t. Inicijalno, krutost opruge u modelu je k = 109 N/m. Na slici 7.32 prikazana su pomeranja čvora 1 i čvora 7, pri kretanju kolica profilom brzina vm1, a na slici 7.33 su data ova pomeranja pri kretanju kolica profilom brzina vm2. a) b) Slika 7.32. Pomeranja a) čvora 1, b) čvora 7; vm1; k=109 N/m k=10 k=10 7 5 0.0 0.2 0.4 0.6 0.8 1.0 -0.06 -0.04 -0.02 0.00 Bezdimenziono vreme tt Po m er an je y m  UX1 UY1 0 5 10 15 20 -0.10 -0.05 0.00 0.05 Vreme t s Po m er an ja -c vo r1 m  UY7 UX7 0 5 10 15 20 -0.10 -0.05 0.00 0.05 Vreme t s Po m er an ja -c vo r7 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 115 a) b) Slika 7.33 Pomeranja a) čvora 1, b) čvora 7; vm2; k=109 N/m Pomeranja sa slike 7.32 su identična odgovarajućim pomeranjima datim na slikama 7.18 i 7.19, a pomeranja sa slike 7.33 su identična odgovarajućim na slikama 7.26 i 7.27, na osnovu čega se može zaključiti da nema razlika u odzivima prema modelu pokretnog oscilatora i modelu pokretne mase. 7.2.B.1 Uticaj krutosti opruge Ukoliko se razmatra profil brzina kretanja vm2, i koeficijent krutosti k=107 N/m, uporedni prikaz pomeranja wy i y je dat na slici 7.34. Veoma male razlike između ovih dijagrama potvrđuju da se ovešena masa ponaša kao i tačka na stazi strukture ispod koje se masa nalazi, pri ovoj krutosti opruge. Slika 7.34. Pomeranje y i wy, vm2; k=107 N/m U X1 U Y1 0 2 4 6 8 10 12 14 -0.10 -0.05 0.00 0.05 Vreme t s Po m er an ja -c vo r1 m  UY7 UX7 0 2 4 6 8 10 12 14 -0.10 -0.05 0.00 0.05 Vreme t s Po m er an ja -c vo r7 m  wy y 0 2 4 6 8 10 12 14 -0.08 -0.06 -0.04 -0.02 0.00 Vreme t s wy ,y m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 116 Daljim smanjivanjem krutosti opruge do vrednosti k=106 N/m dolazi do povećanja pomeranja y, a primetna su velika odstupanja od dijagrama za krutost od k=107 N/m, slika 7.35. Slika 7.35 Uticaj krutosti opruge na pomeranje y; vm2; k=107 N/m, k=106 N/m 7.3 ODZIV MODELA POKRETNOG OSCILATORA SA KLATNOM Odzivi modela portalne dizalice prema modelu pokretnog oscilatora sa klatnom se dobijaju rešavanjem j‐ne (6.34). 7.3.A. ANALIZA ODZIVA ‐ MODEL A Određivanje odziva modela A, tabela 7.2, je izvršeno sa podelom ukupne mase mss na masu kolica m1=3 t, masu vitla m2=5 t i masu tereta m3= 52 t. 7.3.A.1 Odziv sistema usled kretanja kolica konstantnom brzinom Ovde su dobijeni odzivi sistema usled kretanja kolica konstantnom brzinom, za vrednosti brzina v=3 m/s i v=5 m/s. Usvojena je krutost opružnog elementa u sistemu kolica, k =109 N/m, a dužina užetnog sistema iznosi Lu=12 m. Na slici 7.36. prikazana su pomeranja čvora 6 i promena ugla klaćenja φ. k=10 k=106 7 0 2 4 6 8 10 12 14 -0.08 -0.06 -0.04 -0.02 0.00 0.02 Vreme t s ym  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 117 a) b) c) Slika 7.36. a) Horizont. pomeranje čvora 6, UX6, b) Vertik. pomeranje čvora 6, UY6, c)Ugao klaćenja φ; v=const.=3 m/s i v=const.=5 m/s v=3 ms=const. v=5 ms=const. 0.0 0.2 0.4 0.6 0.8 1.0 -0.02 -0.01 0.00 0.01 0.02 0.03 0.04 Bezdimenziono vreme tt Ho riz on ta ln op om er an je -c vo r6 m  v=3 ms=const. v=5 ms=const. 0.0 0.2 0.4 0.6 0.8 1.0 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0.00 0.01 Bezdimenziono vreme tt Ve rti ka ln o po m er an ja -c vo r6 m  v=3 ms =const. v=5 ms=const. 0.0 0.2 0.4 0.6 0.8 1.0 -0.004 -0.002 0.000 0.002 0.004 0.006 Bezdimenziono vreme tt Ug ao j ra d DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 118 Vertikalno pomeranje središnjeg čvora glavnog nosača, slika 7.36.b, pokazuje neznatna odstupanja za navedene konstantne brzine, sa vrednostima koje su veoma bliske vrednostima sa slike 7.11. Vrednosti ugla klaćenja tereta pokazuju rast amplitude za veću brzinu kolica, slika 7.36.c, a maksimalna vrednost koja se javlja na kraju kretanja iznosi 0,004 rad (0,230) predstavlja neznatan otklon što nema uticaja na dinamiku tereta i njegovo pozicioniranje. Minimalna odstupanja nastaju kod horizontalnog pomeranja središneg čvora, tj. horizontalnih pomeranja čvorova staze, slika 7.36.a. Primetan je rast amplitude za veću brzinu kretanja kolica, ali se i ovo može zanemariti u odnosu na vrednost statičkih pomeranja. 7.3.A.2. Uticaj brzine i ubrzanja kolica Ovde se određuju odzivi pri kretanju kolica profilima brzina vm1 i vm2, slika 7.1.a. Krutost opružnog elementa u sistemu kolica iznosi k =109 N/m, a dužina užetnog sistema iznosi Lu=12 m. Na slici 7.37. prikazana je promena ugla klaćenja tereta, a na slikama 7.38 i 7.39 prikazani su vertikalno i horizontalno pomeranje čvora 6, respektivno. Slika 7.37. Promena ugla klaćenja φ; vm1, vm2 v vm2 m1 0.0 0.2 0.4 0.6 0.8 1.0 -0.4 -0.2 0.0 0.2 0.4 Bezdimenziono vreme tt Ug ao j ra d DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 119 Slika 7.38. Vertikalno pomeranje čvora 6, UY6; vm1, vm2 Slika 7.39. Horizontalno pomeranje čvora 6, UX6; vm1, vm2 U tabeli 7.9 prikazane su maksimalne vrednosti koordinata ߮ , UY6, UX6, za profile brzina kretanja vm1 i vm2. Ugao φ dostiže vrednost od 0,25 rad (ൎ 14,30) što predstavlja povećanje od 56 % u odnosu na njegovu najveću vrednost za vm1 (9,10). Najveća vrednost za UY6 iznosi 5,51 cm što predstavlja povećanje od 6,7 % u odnosu na vrednost za vm1. vm2 vm1 0.0 0.2 0.4 0.6 0.8 1.0 -0.06 -0.04 -0.02 0.00 0.02 Bezdimenziono vreme tt Ve rti ka ln o po m er an je -c vo r6 m  v vm1 m2 0.0 0.2 0.4 0.6 0.8 1.0 -0.02 0.00 0.02 0.04 0.06 Bezdimenziono vreme tt Ho riz on ta ln op om er an ja -c vo r6 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 120 Tabela 7.9. Maksimalne vrednosti pomeranja koordinata vm1 vm2 Ugao ߮ 0,16 rad 0,25 rad Maks. UY6 ‐5,16 cm ‐5,51 cm Maks. UX6 4,1 cm 4,4 cm Maksimalno pomeranje za UX6 iznosi 4,4 cm, tabela 7.5, i predstavlja malo povećanje u odnosu na vrednost za vm1 (4,1 cm). Međutim, poređenjem oblika grafika za UX6 i UY6 za model pokretnog oscilatora sa klatnom i za model pokretne mase, slike 7.7 i 7.8 može se primetiti da vertikalno pomeranje UY6 ima slične grafike, a horizontalno pomeranje UX6 veoma različite grafike. Postavkom horizontalne sile koja deluje na strukturu, izraz (5.52), može se izdvojiti njen član koji predstavlja uticaj klaćenja tereta ௫ܲఝ ൌ ݉ଷܮ௨ሺ ሷ߮ cos߮ െ ሺ ሶ߮ ሻଶ sin߮ሻ (7.1) Izostavljanjem ovog uticaja, tj. zanemarivanjem ovog izraza u jednačini (6.34), za vm2, dobijaju se odzivi strukture UX6 i UY6 prikazani na slici 7.40. Slika 7.40. Pomeranja čvora 6, UX6 i UY6; vm2 ; zanemareno ௫ܲఝ Ovako dobijeni grafici, slika 7.41, odgovaraju graficim za UX6 i UY6 na slikama 7.7 i 7.8. UX6 UY6 0 2 4 6 8 10 12 -0.04 -0.02 0.00 0.02 0.04 Vreme s Po m er an ja -c vo r6 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 121 7.3.A.3. Uticaj dužine užetnog sistema Ovde se razmatraju odzivi sistema za različite dužine užetnog sistema (Lu), s obzirom da ova veličina ima uticaja na promenu perioda oscilovanja klatna. Grafici ugla klaćenja za dužine užetnog sistema za Lu =6 m i Lu=9 m su prikazani na slici 7.41. Kolica se kreću profilom vm2. Za matematičko klatno, prema gore navedenim dužinama i izrazu ௜ܶ ൌ 2ߨ/ඥ݃/ܮ௨௜, periodi oscilovanja bi iznosili T1=4,94 s, i T2=6,03 s. Slika 7.41. Promena ugla klaćenja φ; vm2; Lu= 6 m, Lu=9 m Ugao klaćenja dostiže vrednost od 0,4 rad (ൎ 230) što predstavlja veliku vrednost kod dizalica. S obzirom na zaključke da klaćenje tereta utiče najviše na horizontalna pomeranja (poglavlje 7.3.A.2), na sledećoj slici prikazano je pomeranje UX6 za navedene dužine užetnog sistema. Lu=6 m Lu=9 m 0.0 0.2 0.4 0.6 0.8 1.0 -0.4 -0.2 0.0 0.2 0.4 0.6 Bezdimenziono vreme tt Ug ao kla ce nj a j ra d DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 122 Slika 7.42. Horizontalno pomeranje čvora 6 ‐ UX6; vm2; Lu= 6 m, Lu=9 m 7.3.A.4. Uticaj krutosti opruge Ovde se razmatra slučaj kada je krutost opruge k=107 N/m, a kolica se kreću profilom vm2, za minimalnu dužinu užeta od Lu=6 m. Na slici 7.43 prikazana su pomeranja wy i y, u ovom slučaju. Slika 7.43. Pomeranje wy i y, k=107 N/m; vm2; Lu=6 m Lu= 6 mLu=9 m 0.0 0.2 0.4 0.6 0.8 1.0 -0.02 0.00 0.02 0.04 0.06 Bezdimenziono vreme tt Ho riz on ta ln op om er an je -c vo r6 m  w y y 0 2 4 6 8 10 12 -0.05 0.00 0.05 Vreme t s wy , ym  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 123 Koordinata y pokazuje velika odstupanja u odnosu na vertikalno pomeranje duž staze wy, za razliku od karaktera dijagrama prikazanog na slici 7.31, i dostiže vrednost od 9 cm pri kraju ciklusa kretanja. Ukoliko se izraz ௌܲఝ ൌ ݉ଷܮ௨ሺ ሷ߮ sin ߮ ൅ ሺ ሶ߮ ሻଶ cos ߮ሻ zanemari u jednačinama (6.20) i (6.34), na sledećoj slici je prikazan grafik za y i wy dobijen preko (6.39). S obzirom da se grafici neznatno razlikuju, može se zaključiti da odziv za koordinatu y na slici 7.43 predstavlja uticaj klaćenja tereta. Slika 7.44 Pomeranje wy i y, k=107 N/m; vm2; Lu=6 m; zanemareno PSφ 7.3.B ANALIZA ODZIVA ‐ MODEL B Masa kolica je podeljena na masu kolica m1=2 t, masu vitla m2=4 t i masu tereta m3= 34 t. 7.3.B.1 Odziv sistema usled kretanja kolica konstantnom brzinom Ovde su dobijeni odzivi sistema usled kretanja kolica konstantom brzinom, za vrednosti brzina v=3 m/s i v=5 m/s. Za ovaj slučaj usvojena je krutost opružnog elementa u sistemu kolica od k =109 N/m, a dužina užetnog sistema iznosi Lu=10. Na slici 7.45 su prikazana pomeranja čvora 1, čvora 7 i promena ugla klaćenja φ. wy y 0 2 4 6 8 10 12 -0.05 -0.04 -0.03 -0.02 -0.01 0.00 Vreme t s wy , ym  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 124 Vertikalno pomeranje čvora 1, slika 7.45.b, pokazuje neznatna odstupanja za navedene konstantne brzine, sa vrednostima koje su veoma bliske vrednostima sa slike 7.25. Vertikalno pomeranja središnjeg čvora rasponog dela glavnog nosača, slika 7.45.c, takođe pokazuje mala odstupanja. Ugao klaćenja tereta, slika 7.45.d, pokazuje odstupanja za navedene brzine, ali maksimalna vrednost od 0,0073 rad (0,420) predstavlja neznatan otklon i može se reći da nema uticaja na dinamiku tereta i njegovo pozicioniranje. Sva pomeranja sa slike 7.45 pokazuju povećanje amplitude sa povećanjem brzine na v=5 m/s, a najveće povećanje ima horizontalno pomeranje čvora 1 (cele staze). Maksimalna vrednost ovog pomeranja iznosi 7,21 cm kada se teret nalazi na prepustu i 3,96 cm kada se tereta nalazi na rasponom delu nosača. a) b) v=3 ms=const. v=5 ms=const. 0.0 0.2 0.4 0.6 0.8 1.0 -0.10 -0.05 0.00 0.05 Bezdimenziono vreme tt Ho riz on ta ln op om er an je -c vo r1 m  v=3 ms =const. v=5 ms=const. 0.0 0.2 0.4 0.6 0.8 1.0 -0.10 -0.05 0.00 0.05 Bezdimenziono vreme tt Ve rti ka ln o po m er an je -c vo r1 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 125 c) d) Slika 7.45. Pomeranja sistema, a) Horiz. pomeranje čvora 1‐UX1, b) Vertikalno pomeranje čvora 1 ‐ UY1, c) Vertikalno pomeranje središnjeg čvora rasponog dela glavnog nosača ‐ UY7, d) Ugao klaćenja tereta φ ; v =const.=3 m/s, v =const.=5 m/s 7.3.B.2. Uticaj brzine i ubrzanja kolica Razmatraju se odzivi sistema usled kretanja kolica profilima brzina vm1 i vm2, slika 7.1.a. Na slici 7.46 je prikazana promena ugla klaćenja, pri čemu se vrednosti povećavaju sa porastom brzina i ubrzanja. v=3 ms=const. v=5 ms=const. 0.0 0.2 0.4 0.6 0.8 1.0 -0.06 -0.04 -0.02 0.00 0.02 0.04 Bezdimenziono vreme tt Ve rti ka ln o po m er an ja -c vo r7 m  v=5 ms=const. v=3 ms=const. 0.0 0.2 0.4 0.6 0.8 1.0 -0.010 -0.005 0.000 0.005 0.010 Bezdimenziono vreme tt Ug ao j ra d DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 126 Slika 7.46. Promena ugla klaćenja, φ; vm1, vm2 Na slikama 7.47 i 7.48 prikazani su vertikalno i horizontalno pomeranje čvora 1, respektivno. Slika 7.47. Horizontalno pomeranje čvora 1; vm2, vm1 Slika 7.48. Vertikalno pomeranje čvora 1; vm2, vm1 vm2 vm1 0.0 0.2 0.4 0.6 0.8 1.0 -0.4 -0.2 0.0 0.2 0.4 Bezdimenziono vreme tt Ug ao j ra d vm2vm1 0.0 0.2 0.4 0.6 0.8 1.0 -0.10 -0.05 0.00 0.05 0.10 Bezdimenziono vreme tt Ho riz on ta ln op om er an je -c vo r1 m  vm2 vm1 0.0 0.2 0.4 0.6 0.8 1.0 -0.10 -0.05 0.00 0.05 Bezdimenziono vreme tt Ve rti ka ln o po m er an je -c vo r1 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 127 Na slici 7.49 je dato vertikalo pomeranje središnjeg čvora rasponog dela glavnog nosača, tj. čvora 7. Slika 7.49. Vertikalno pomeranje čvora 7; vm2, vm1 U tabeli 7.10 prikazane su maksimalne vrednosti koordinata φ, UX1, UY1, UY7, za profile brzina vm1 i vm2. Tabela 7.10. Maksimalne vrednosti pomeranja koordinata vm1 vm2 Ugao ߮ 0,128 rad 0,26 rad Maks. UX1 ‐8,11 cm ‐9,25 cm Maks. UY1 ‐7,93 cm ‐8,07 cm Maks. UY7 ‐5,04 cm ‐5,46 cm Ugao ߮ dostiže vrednost od 0,26 rad (14,90) što predstavlja povećanje od 103 % u odnosu na vrednost za vm1 (7,330). Maksimalno pomeranje za UX1 iznosi 9,25 cm i prestavlja povećanje od 14 % u odnosu na vrednost za vm1. Takođe, pomeranje UX1, kada se kolica nalaze na rasponom delu nosača, dostiže vrednost od 7,82 cm (slika 7.47) za vm2, a 6,17 cm za vm1. Maksimalna pomeranja UY1 nastaju u trenutku kada se kolica nalaze na slobodnom kraju prepusta, a kada se kolica nalaze na rasponom delu nosača ova pomeranja iznose 4,34 cm za vm2 i 3,68 cm za vm1. vm2 vm1 0.0 0.2 0.4 0.6 0.8 1.0 -0.06 -0.04 -0.02 0.00 0.02 0.04 Bezdimenziono vreme tt Ve rti ka ln o po m er an je -c vo r7 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 128 Na slici 7.50 su prikazani odzivi karakterističnih čvorova strukture pri čemu su rezultati dobijeni izostavljanjem izraza (7.1) u jednačini (6.34), za vm2. Slika 7.50. Pomeranja UX1, UY1 i UY7; vm2; zanemareno Pxφ U ovom slučaju, grafici za pomeranja UX1, UY1, UY7 odgovaraju graficima prikazanim na slikama 7.21, 7.22 i 7.23. 7.3.B.3. Uticaj dužine užetnog sistema Razmatranje odziva sistema za različite dužine užetnog sistema je izvršeno za Lu=14 m i Lu=11 m. Kolica se kreću profilom brzina vm2, a krutost opruge u sistemu kolica iznosi k=109 N/m. Za matematičko klatno, prema dužinama Lu1=14m i Lu2=11 m i izrazu ௜ܶ ൌ 2ߨ/ඥ݃/ܮ௨௜, periodi oscilovanja bi iznosili T1=7,5s i T2=6,64 s. Grafici ugla klaćenja su prikazani na sledećoj slici. UX1 U UY1 Y7 0 2 4 6 8 10 12 14 -0.15 -0.10 -0.05 0.00 0.05 Vreme , t s Po m er an ja cv or -1 iv er tik al no po m .c vo r-7 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 129 Slika 7.51. Promena ugla klaćenja φ; vm2; k=109 N/m; Lu=14 m i Lu=11 m Ugao klaćenja, pri dužini užetnog sistema od 11 m, dostiže vrednost od 0,5 rad (ൎ 28,60) što predstavlja veliki otklon tereta. Slika 7.52. Horizontalno pomeranje čvora 1; Lu=14 m i Lu=11 m Najveća vrednost horizontalnog pomeranja čvora 1, za Lu=11m, iznosi 11,88 cm i nastaje u trenutku kada i klaćenje tereta dostiže svoju najveću vrednost, slika 7.45. Niže vrednosti dužine užetnog sistema bi generisale još veće vrednosti ovih pomeranja, ali ovo nema praktičnog značaja zbog postojanja mehanizama za sprečavanje klaćenja tereta koji svoju funkciju obavljaju upravo na malim dužinama užetnog sistema. Lu=11 m Lu=14 m 0.0 0.2 0.4 0.6 0.8 1.0 -0.4 -0.2 0.0 0.2 0.4 0.6 Bezdimenziono vreme tt Ug ao kla ce nj a j ra d Lu=11 m Lu=14 m 0.0 0.2 0.4 0.6 0.8 1.0 -0.10 -0.05 0.00 0.05 0.10 0.15 Bezdimenziono vreme tt Ho riz on ta ln op om er an je -c vo r1 m  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 130 7.3.B.4 Uticaj krutosti opruge Za parametre sistema koji su postavljeni u prethodnom poglavlju, Lu=11 m i profil brzina kretanja vm2, dodatno se razmatra odziv sistema kada je krutost opruge k=106 N/m. Slika 7.53 Pomeranje wy i y, k=106 N/m; vm2; Lu=11 m Koordinata y pokazuje velika odstupanja u odnosu na vertikalno pomeranje wy i dostiže vrednost od 6,8 cm pri kraju ciklusa kretanja. w y y 0 2 4 6 8 10 12 14 -0.08 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 Vreme t s wy , ym  DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 131 8 ZAKLJUČAK U disertaciji se analizira dinamičko ponašanje konstrukcije portalne dizalice usled dejstva kolica kao pokretnog opterećenja. U užem smislu, razmatrane su portalne dizalice za kontejnerske terminale zbog visokih performansi ovih mašina koje se ogledaju u velikim rasponima, dohvatima, visinama dizanja, a posebno zbog izuzetnih brzina kretanja kolica. Oscilacije ramovskih struktura su nedovoljno opisane u literaturi iz oblasti dinamike konstrukcija. Zbog toga je u ovom radu izvršeno određivanje karakteristične jednačine poprečnih oscilacija ravanske ramovske strukture sa jednim prepustom preko sistema elastičnih tela tipa prizmatične grede. Određene su prve tri frekvencije za pretpostavljene geometrijske varijante, tabela 3.1, koje mogu poslužiti u budućim istraživanjima iz ove oblasti kao verifikaciona podloga za određivanje frekvencija modela. U savremenoj literaturi koja je usko povezana sa predmetom istraživanja u ovoj disertaciji, prezentovan je kombinovani pristup za postavljanje jednačina kretanja sistema. U ovom radu je kombinovani pristup takođe usvojen za formiranje modela portalnih dizalica i objašnjen u potpunosti. Za razliku od koncepta svođenja složene strukture na jednostavne modele uz veliki broj aproksimacija koje dovode u pitanje tačnost rešenja, pristup u ovoj disertaciji daje mogućnost razmatranja strukture u celoj svojoj složenosti. Pored toga, eksplicitno su prikazani algoritmi za određivanje karakterističnih matrica sa elementima promenljivim u vremenu, glava 5, čime jedan od doprinosa ovog rada predstavlja i razvoj kombinovanog pristupa problemima dinamike dizalica. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 132 Prema usvojenoj koncepciji formirana su dva konačnoelementna modela strukture portalne dizalice koja predstavljaju dva osnovna konstrukciona tipa ovih dizalica. Kolica su razmatrana preko tri modela na osnovu čega su i formirani modeli sistema kao model pokretne mase, model pokretnog oscilatora i model pokretnog oscilatora sa klatnom. Dinamička interakcija između kolica i strukture je definisana za svaki od modela. Konačnoelementi modeli struktura su definisani preko odgovarajućih matrica, u analitičkom obliku. Kombinovanjem sa dinamičkim modelima kolica preko dinamičke interakcije ovih dvaju modela postavljene su diferencijalne jednačine kretanja sistema. Za sva tri modela postavljene jednačine su drugog reda sa promenljivim koeficijentima u vremenu, a za model pokretnog oscilatora sa klatnom dobijena je nelinearna diferencijalna jednačina. Matematički model sistema prezentovan u ovom radu kao pokretni oscilator sa klatnom predstavlja nadgradnju postojećih savremenih modela sistema u problemima pokretnog opterećenja jer je pored transverzalnih i podužnih oscilacija elemenata strukture i inercijalnih efekata masa pokretnog sistema uzeto u obzir i sledeće: elastičnost u konstrukciji kolica ili prisustvo opruge u sistemu kolica i uticaj klaćenja tereta. U ovom radu je kretanje kolica uključeno kroz trapezni profil brzina kretanja, slika 4.11, kojim je obuhvaćeno i ubrzanje i usporenje kolica. Razmotren je, takođe, i slučaj kretanja kolica konstantnom brzinom na modelima pokretne mase i pokretnog oscilatora sa klatnom. Rezultati i napomene iz poglavlja 7.1.A.1 i 7.1.B.1, ukazuju da je najveće vertikalno pomeranje strukture identično a najveće horizontalno pomeranje približno odgovarajućim pomeranjima dobijenim statičkim proračunom. Dodatno, u poglavljima 7.3.A.1 i 7.3.B.1 su dobijeni odzivi usled kretanja kolica konstantnom brzinom v=3 m/s i v=5 m/s, gde se pored navedenih pomeranja dobijaju zanemarljive vrednosti za klaćenje tereta. S obzirom da većina radova iz oblasti pokretnog opterećenja pretpostavlja kretanje kolica konstantnom brzinom, ovde je potrebno zaključiti sledeće:  Sa pretpostavkom kretanja kolica konstantnom brzinom do 5 m/s, što predstavlja inače ekstremnu performansu kod dizalica, nije potrebno uvoditi problem pokretnog opterećenja u analizu dinamičkog ponašanja dizalica DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 133 Rezultati iz poglavlja 7.1.A.2 i 7.1.B.2 ukazuju da povećanje strukturalnog prigušenja u KE modelu struktura dovodi do smanjenja amplituda horizontalnih pomeranja. Iako se ovo može smatrati povoljnim uticajem na konstrukciju, potrebno je napomenuti da bi slične analize trabalo da budu praćene eksperimentalnom verifikacijom, a posebno razmatranjem novih materijala ili dizajnom veza u strukturi kojim bi se povećalo strukturalno prigušenje što može predstavljati osnovu za dalja istraživanja iz oblasti dizalica. Generalno, za uobičajene vrednosti strukturalnog prigušenja (do 2 %), može se zaključiti da ovaj uticaj nije potrebno uključiti u analizu dinamičkog ponašanja dizalica. Uticaj krutosti opruge u sistemu kolica ili elastičnost konstrukcije kolica, poglavlja 7.2.A.1 i 7.2.B.1, na pomeranja staze je zanemarljiv, a bitan samo za vertikalno pomeranje tereta. Kako ova poslednja napomena važi samo za vrednosti koeficijenta krutosti koji predstavlja nerealan podatak kod dizalica, može se zaključiti da nije potrebno razmatrati ovaj uticaj u modelima portalnih dizalica, a samim tim i model pokretnog oscilatora je od malog značaja. Sa uvođenjem ekstremnih performansi kolica, nominalnom brzinom od vnom=5m/s i ubrzanjem/usporenjem a=1,25 m/s2, odzivi strukture dobijeni pomoću modela pokretne mase pokazuju:  Povećanje performanski kretanja kolica ima slab uticaj na vertikalna pomeranja strukture, i samim tim nije potrebno razmatrati ovaj slučaj pri uobičajenom dokazu deformacija strukture.  Ovaj uticaj je najveći i veoma bitan za horizontalna pomeranja strukture, pri čemu maksimalne vrednosti ovih pomeranja mogu biti veće od vertikalnih pomeranja. Ovo je pokazano i analizom rezultata pomoću modela pokretnog oscilatora sa klatnom, poglavlja 7.3.A.2 i 7.3.B.2, gde se pored ovih zaključaka može primetiti i povećanje klaćenja tereta sa povećanjem performansi kretanja kolica. Ovaj matematički model, (6.34), obezbeđuje adekvatan odgovor sistema, jer pokazuje dvojaki uticaj povećanja ubrzanja, tj. povećanje ubrzanja utiče na povećanje dinamičkih veličina klatna koji se dalje odražavaju na DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 134 horizontalno pomeranje strukture menjajući pritom karakter pomeranja koji se dobija ako se zanemari ovaj uticaj. Zbog toga, ovaj model ima najveći značaj od svih postavljenih modela, kako u ovom radu, tako i odnosu na modele iz literature [44,45]. Odzivi su dobijeni na modelima struktura koje predstavlju realne konstrukcije, pa za razliku od modela prikazanim u [42,43], ovo se može izdvojiti kao realan pristup u problemima pokretnog opterećenja. Vrednosti ugla klaćenja tereta dobijene, ovim, tzv. numeričkim eksperimentom potvrđuju da je primena mehanizama za sprečavanje njihanja tereta obavezna kod portalnih dizalica za kontejnerske terminale. Matematički model pokretnog oscilatora sa klatnom se može primeniti i za druge tipove konstrukcije portalnog tipa (obalske kontejnerske dizalice, pretovarni mostovi) kao osnova (ne samo polazna) za formiranje upravljačkih algoritama za kretanje kolica. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 135 9 LITERATURA [1] Ostrić, D., Tošić, S.: Dizalice, Mašinski fakultet, Beograd, 2005. [2] Deutsche Bank Research, Container shipping: Successsful turnaround. http://www.dbresearch.com/PROD/DBR_INTERNET_EN‐ PROD/PROD0000000000271589.pdf [3] Georgijević, M: Pretovar kontejnera, skripte, FTN Novi Sad, 2011. [4] Bhimani, A.K., Hsieh, J.K.: Cranes to Serve Ship in the Slip Ceres Paragon Terminal, Amsterdam, Proc. of the Conference, PORTS '01, ASCE, section 30, chapter 3, Norforlk, USA, 2001. [5] Krylov A.N.: Mathematical collection of Papers of the Academy of Sciences, Vol. 61, Peterburg, 1905. [6] Timoshenko S.P.: Forced vibration of Prismatic Bars (in Russian), Izvestiya Kievskogo politekhnicheskogo instituta, 1908. [7] Inglis C.E.: A Mathematical Treatise on Vibration in Railway Bridges, Cambridge University Press, 1934 [8] Timoshenko S.P.: On the forced Vibrations of Bridges, Philosoph.Magazine, Ser. 6, 43, 1922. [9] Goldenblat I.I.: Some New Problems of Structural Dynamics (in Russian), Izvestiya AN SSSR, Otd.techn.nauk, No.6, 819‐833,1950. [10] Bolotin V.V.: On the effect of a Moving Load on Bridges (in Russian), Trudy Moskovskogo instituta inzhenerov zheleznodorozhnogo transporta, Vol. 74, 269‐296, 1950. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 136 [11] Fryba, L.: Vibration of solids and structures under moving loads, 3rd edition, Thomas Telford, 1999. [12] Olsson, M.: On the fundamental moving load problem, Journal of Sound and Vibration, 14 5 (2), 299‐307, 1991. [13] Olsson, M.: Analysis of structures subjected to moving loads, Doctoral Thesis, Lund Institute of technology, Sweden, 1986. [14] Wu, J.J., Whittaker A.R., Cartmell M.P.: The use of finite element techniques for calculating the dynamic response of structures to moving loads, Computers and Structures, 78, 789‐799, 2000. [15] Gašić, V., Zrnić N., Obradović A., Milovančević M.: Revisiting the use of finite element packages for moving load problem at bridge cranes, XIX International Conference MHCL 09, 71‐74, 2009. [16] Pesterev A.V., Yang B., Bergman L.A., Tan C.A.: Revisiting the moving force problem, Journal of Sound and Vibration 261, 75‐91, 2003. [17] Stokes G.G.: Discussion of a Differential Equation Relating to the Breaking of Railway Bridges. Trans. Cambridge Philosoph. Soc.,8, Part 5, 1849. [18] Stanišić M.M., Harding J.C.: On the Response of Beams to an Arbitrary Number of Concentrated Moving Masses, Journal of Franklin Institute, Vol 287, No.2, 115‐123, 1969. [19] Akin, J.E., Mofid, M.: Numerical solution for response of beams with moving mass, Journal of Structural Engineering. (1989), pp. 120‐131. [20] Mofid, M., Akin J.E.: Discrete element response of beams with traveling mass, Advances in Engineerign Software, 25, 321‐331, 1996. [21] Lee, H.P.: Dynamic Response of a Beam with a moving Mass, Journal of Sound and Vibration, 191 (2), 289‐294, 1996. [22] Esmailzadeh, E., Ghorashi, M.: Vibration analysis of beams traversed by a moving mass, Journal of Engineering 8, 213‐220, 1995. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 137 [23] Lee, U.: Separation Between the Flexible Structure and the Moving Mass Sliding on It, Journal of Sound and Vibration, 209 (5), 867‐877, 1998. [24] Michaltsos, G.T., Kounadis, A.N.: The Effects of Centripetal and Coriolis Forces on the Dynamic Response of Light Bridges Under Moving Loads, Journal of Vibration and Control, 7, 315‐326, 2001. [25] Michaltos, G., Sophianopoulos, D., Kounadis, A.N.: The effect of a moving mass and other parameters on the dynamic response of a simply supported beam, Journal of Sound and Vibration. (1996); 191:357‐362. [26] Michaltsos, G.T.: Dynamic Behaviour of a Single‐Span Beam Subjected to Loads Moving with Variable Speeds, Journal of Sound and Vibration, 258 (2), 359‐371, 2002. [27] Zrnić, N.: Uticaj kretanja kolica na dinamičko ponašanje obalskih kontejnerskih dizalica, Doktorska disertacija, Mašinski fakultet Beograd, 2005. [28] Zrnić, N., Hoffmann, K, Bošnjak, S.: Modelling of dynamic interaction between structure and trolley for mega container cranes, Math. Comput. Model. Dyn. Syst. 15(2009), pp. 295‐311. [29] Zrnić, N., Bošnjak, S., Hoffmann, K.: Parameter sensitivity analysis of non‐ dimensional models of quayside container cranes, Math. Comput. Model. Dyn. Syst., Vol. 16, No. 2, 2010, pp. 145‐160. [30] Zrnić, N., Gašić, V., Obradović, A., Bošnjak, S.: Appropriate modeling of dynamic behaviour of quayside container cranes boom under a moving trolley. // Springer Proceedings in Physics 139, Vibration problems ICOVP 2011, pp.81‐86. [31] Wu, J.J., Whittaker, A.R., Cartmell, M.P.: Dynamic responses of structures to moving bodies using combined finite element and analytical methods, International Journal of Mechanical Sciences 43, 2555‐2579, 2001. [32] Clough, RW., Penzien, J.: Dynamics of structures, McGraw‐Hill, New York, 1993. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 138 [33] Wu, J.J.: Finite element analysis and vibration testing of a three‐ dimensional crane structure, Measurements 39, 740‐749, 2006. [34] Gašić V., Znić N., Rakin M.: Consideration of a Moving Mass Effect on Dynamic Behaviour of a Jib Crane Structure, Tehnički Vjesnik‐Technical Gazette, 19 (1), 115‐121, 2012. [35] Lin, Y.‐H., Tretheway, M. W.: Finite element analysis of elastic beams subjected to moving dynamic loads, Journal of Sound and Vibration, 136 (2), 323‐342, 1990. [36] Rieker, J.R, Lin, Y‐H., Tretheway, M.W.: Discretization considerations in moving load finite element beam models, Finite Element in Analysis and Design 21, 129‐144, 1996. [37] Gašić, V., Zrnić, N., Obradović, A., Bošnjak, S.: Consideration of Moving Oscillator Problem in Dynamic Responses of Bridge Cranes, FME Transactions, 39(1), 17‐24, 2011. [38] Stancioiu, D., Ouyang H., Mottershead, J.E.: Vibration of a beam excited by a moving oscillator considering separatio and reattachment, Journal of Sound and Vibration 2310, 1128‐1140, 2008. [39] Pesterev, A.V., Bergman, L.A., Tan, C.A, Tsao, T.‐C., Yang, B.: On asymptotics of the solution of the moving oscillator problem, Journal of Sound and Vibration 260, 519‐536, 2003. [40] Pesterev, A.V., Bergman, L.A., Tan, C.A, Tsao, T.‐C., Yang, B.: Some recent results in moving load problems with application to highway bridges, Proc Int Conf Motion Vib Control, 6, Vol.1, K20‐K27, 2002. [41] Bugaruć, U.: Prilog optimizaciji procesa istovara rasutih tereta u rečnim lukama, Magistarska teza, Mašinski fakultet u Beogradu, 1996. [42] Oguamanam, D. C. D., Hansen, J. S.: Dynamic response of an overhead crane system, Journal of Sound and Vibration, 213(5) (1998), pp. 889‐906. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 139 [43] Oguamanam, D.C.D, Hansen, J.S., Heppler, G.R.: Dynamics of a three‐ dimensional overhead crane system, Journal of Sound and Vibration, 242(3), 411‐426, 2001. [44] Wu, J.J.: Dynamic responses of a three‐dimensional framework due to a moving carriage hoisting a swinging object, International journal for numerical methods in engineering, 59, 1679‐1702, 2004. [45] Wu, J.J.: Transverse and longitudinal vibrations of a frame structure due to a moving trolley and the hoisted object using moving finite element, International Journal of Mechanical Sciences 50, 613‐625, 2008. [46] Esen, I.: Dynamic response of a beam due to an accelerating moving mass using moving finite element approximation, Mathematical and Computational Applications, Vol.16, No.1, pp.171‐182, 2011. [47] Georgijević, M., Bojanić, V., Bojanić, G., Novković M.: Container cranes for river ports, control and calculation of life time as the optimization base, International conference MHCL 09, Belgrade, pp. 53‐60, 2009. [48] Yazid, E., Parman, S., Fuad, K.: Vibration Analysis of flexible Gantry Crane System Subjected Swinging Motion of Payload, Journal of Applied Sciences, 10, 1707‐1715, 2011. [49] Younesian, D., Ghafoori, E., Sadeghpour, M.: Nonlinear vibration of a three‐ dimensional moving gantry crane subjected to a travelling trolley hoisting a swinging object, Transactions of the Canadiana Society for Mechanical Engineering, Vol. 34, No. 3‐4, pp. 333‐350, 2010. [50] Vuković, J., Obradović, A.: Linearne oscilacije mehaničkih sistema, Mašinski fakultet u Beogradu, 2007. [51] Radosavljević, LJ.: Teorija oscilacija, Mašinski fakultet, Beograd, 1981. [52] Rašković, D.: Teorija oscilacija, Naučna knjiga, Beograd, 1957. [53] Meirovich, L.: Elements of vibration analysis, McGraw‐Hill, 1986. DINAMIČKA INTERAKCIJA NOSEĆE STRUKTURE I KOLICA PORTALNIH DIZALICA VISOKIH PERFORMANSI 140 [54] Gašić, V., Obradović, A., Petković, Z.: Mathematical modelling of the in‐ plane vibrations of portal cranes with FEM verification, Machine Design 2009, Novi Sad, Serbia, 2009., pp. 121‐126 [55] Karnovsky, I., Lebed, O.: Formulas for structural dynamics, McGraw‐Hill, 2004. [56] Petković, Z., Ostrić, D.: Metalne konstrukcije u mašinogradnji, Mašinski fakultet u Beogradu, 1996. [57] Petković, Z.: Metalne konstrukcije u mašinogradnji 2, Mašinski fakultet u Beogradu, 2005. [58] Paz, M., Leigh, W.: Structural dynamics, theory and computation, Springer, 2004. [59] Spyrakos, C, Raftoyiannis, J.: Linear and nonlinear finite element analysis in engineering practice, Algor Inc. Pittsburgh, 1997. [60] Liu, G.R., Quek, S.S.: The finite element method, a practical course, Butterworth‐Heinemann, 2003. [61] Przemieniecki, JS.: Theory of Matrix Structural Analysis, McGraw‐Hill, New York, 1985. [62] Bath, KJ: Finite element procedures in engineering analysis, Prentice‐Hall, 1982. [63] Georgijević, M.: Dinamika dizalica, eksperimentalna i modelska analiza, Zadužbina Andrejević, 1996. [64] Wolfram Mathematica 6 Documentation. http://www. wolfram.com [65] http://www.csiberkeley.com/sap2000 [66] Spyrakos, C.: Finite Element Modeling in Engineering Practice, Algor, Inc., Pittsburgh, USA, 1994. [67] Ognjanović, M.: Mašinski elementi, Mašinski fakultet, Beograd, 2003.               BIOGRAFIJA AUTORA        Vlada M. Gašić  je rođen 8.10.1975. godine u Beogradu. Osnovnu  i srednju  školu (gimnazija‐matematički smer) je završio u Kraljevu sa odličnim uspehom.    Mašinski  fakultet u Beogradu  je upisao 1994. godine. Diplomirao  je 1999.  godine na Katedri za mehanizaciju sa prosečnom ocenom 9,26. Diplomski  rad  iz  oblasti  transpotnih mašina  je  odbranio  sa  ocenom  10.  Kao  student  generacije  1994/95 dobitnik je plakete prof. dr Vojislav K. Stojanović.     Magistarske  studije na Mašinskom  fakultetu na  smeru  za mehanizaciju  je  upisao  1999.  godine,  a  2004.  godine  stekao  je  titulu magistra  tehničkih  nauka,  odbranivši  magistarsku  tezu  pod  nazivom  "Analiza  dinamičkog  ponašanja  pretovarnih mostova  za ugalj u  termoelektranama", urađenu pod mentorstvom  prof. dr Zorana Petkovića.    U  Institutu  za  mehanizaciju  angažovan  je  2000.  godine  kao  istraživač‐ talenat.  Na  Katedri  za  mehanizaciju  je  izabran  za  asistenta‐pripravnika  2001.  godine. Na  istoj  katedri  je  izabran  za  asistenta  2005.  godine  i  reizabran  2009.  godine.     Tokom rada na Mašinskom fakultetu je bio angažovan u izvođenju nastave  (vežbi)  na  sledećim  predmetima:  Metalne  konstrukcije,  Transportne  mašine,  Osnove metalnih konstrukcija u mašinogradnji, Računarsko projektovanje mašina  za  transport  i mehanizaciju, Transportni uređaji  i Mašinski elementi. Kao autor  i  saradnik učestvovao je u realizaciji više naučnih i stručnih radova, kao i projekata  za  privredu.  Od  2006.  godine  poseduje  licencu  odgovornog  projektanta  transportnih sredstava, skladišta i mašinskih konstrukcija i tehnologije.    Vojni rok je odslužio 2002/03. godine. Aktivno govori engleski jezik i služi se  francuskim jezikom. Oženjen je i otac je dve kćerke.                    Прилог 1. Изјава о ауторству Потписани ВЛАДА М. ГАШИЋ број индекса Изјављујем да је докторска дисертација под насловом ДИНАМИЧКА ИНТЕРАКЦИЈА НОСЕЋЕ СТРУКТУРЕ И КОЛИЦА ПОРТАЛНИХ ДИЗАЛИЦА ВИСОКИХ ПЕРФОРМАНСИ  резултат сопственог истраживачког рада,  да предложена дисертација у целини ни у деловима није била предложена за добијање било које дипломе према студијским програмима других високошколских установа,  да су резултати коректно наведени и  да нисам кршио ауторска права и користио интелектуалну својину других лица. Потпис докторанда У Београду, 01.07.2012.год. мр Влада М. Гашић Прилог 2. Изјава o истоветности штампане и електронске верзије докторског рада Име и презиме аутора Влада М. Гашић Број индекса _________________________________________________________ Студијски програм ____________________________________________________ Наслов рада Динамичка интеракција носеће структуре и колица порталних дизалица високих перформанси Ментор проф. др Ненад Зрнић Потписани мр Влада Гашић Изјављујем да је штампана верзија мог докторског рада истоветна електронској верзији коју сам предао за објављивање на порталу Дигиталног репозиторијума Универзитета у Београду. Дозвољавам да се објаве моји лични подаци везани за добијање академског звања доктора наука, као што су име и презиме, година и место рођења и датум одбране рада. Ови лични подаци могу се објавити на мрежним страницама дигиталне библиотеке, у електронском каталогу и у публикацијама Универзитета у Београду. Потпис докторанда У Београду, 01.07.2012. год. мр Влада М. Гашић Прилог 3. Изјава о коришћењу Овлашћујем Универзитетску библиотеку „Светозар Марковић“ да у Дигитални репозиторијум Универзитета у Београду унесе моју докторску дисертацију под насловом: Динамичка интеракција носеће структуре и колица порталних дизалица високих перформанси која је моје ауторско дело. Дисертацију са свим прилозима предао сам у електронском формату погодном за трајно архивирање. Моју докторску дисертацију похрањену у Дигитални репозиторијум Универзитета у Београду могу да користе сви који поштују одредбе садржане у одабраном типу лиценце Креативне заједнице (Creative Commons) за коју сам се одлучио. 1. Ауторство 2. Ауторство - некомерцијално 3. Ауторство – некомерцијално – без прераде 4. Ауторство – некомерцијално – делити под истим условима 5. Ауторство – без прераде 6. Ауторство – делити под истим условима (Молимо да заокружите само једну од шест понуђених лиценци, кратак опис лиценци дат је на полеђини листа). Потпис докторанда У Београду, 01.07.2012.год. мр Влада М. Гашић 1. Ауторство - Дозвољавате умножавање, дистрибуцију и јавно саопштавање дела, и прераде, ако се наведе име аутора на начин одређен од стране аутора или даваоца лиценце, чак и у комерцијалне сврхе. Ово је најслободнија од свих лиценци. 2. Ауторство – некомерцијално. Дозвољавате умножавање, дистрибуцију и јавно саопштавање дела, и прераде, ако се наведе име аутора на начин одређен од стране аутора или даваоца лиценце. Ова лиценца не дозвољава комерцијалну употребу дела. 3. Ауторство - некомерцијално – без прераде. Дозвољавате умножавање, дистрибуцију и јавно саопштавање дела, без промена, преобликовања или употребе дела у свом делу, ако се наведе име аутора на начин одређен од стране аутора или даваоца лиценце. Ова лиценца не дозвољава комерцијалну употребу дела. У односу на све остале лиценце, овом лиценцом се ограничава највећи обим права коришћења дела. 4. Ауторство - некомерцијално – делити под истим условима. Дозвољавате умножавање, дистрибуцију и јавно саопштавање дела, и прераде, ако се наведе име аутора на начин одређен од стране аутора или даваоца лиценце и ако се прерада дистрибуира под истом или сличном лиценцом. Ова лиценца не дозвољава комерцијалну употребу дела и прерада. 5. Ауторство – без прераде. Дозвољавате умножавање, дистрибуцију и јавно саопштавање дела, без промена, преобликовања или употребе дела у свом делу, ако се наведе име аутора на начин одређен од стране аутора или даваоца лиценце. Ова лиценца дозвољава комерцијалну употребу дела. 6. Ауторство - делити под истим условима. Дозвољавате умножавање, дистрибуцију и јавно саопштавање дела, и прераде, ако се наведе име аутора на начин одређен од стране аутора или даваоца лиценце и ако се прерада дистрибуира под истом или сличном лиценцом. Ова лиценца дозвољава комерцијалну употребу дела и прерада. Слична је софтверским лиценцама, односно лиценцама отвореног кода.