The prediction of conveyor trajectories

Size: px
Start display at page:

Download "The prediction of conveyor trajectories"

Transcription

1 University of Wollongong Researh Online Faulty of Engineering - Papers (Arhive) Faulty of Engineering and Information Sienes 007 The predition of onveyor trajetories David B. Hastie University of Wollongong, dhastie@uow.edu.au Peter W. Wypyh University of Wollongong, peter_wypyh@uow.edu.au Puliation Details Hastie, DB & Wypyh, PW, The predition of onveyor trajetories, International Materials Handling Conferene: Belton 14, 007, p 1-16, South Afria: International Materials Handling. Researh Online is the open aess institutional repository for the University of Wollongong. For further information ontat the UOW Lirary: researh-pus@uow.edu.au

2 The Predition of Conveyor Trajetories DB Hastie and PW Wypyh Centre for Bulk Solids and Partiulate Tehnologies Faulty of Engineering, University of Wollongong Northfields Avenue, Wollongong, NSW, 5 1 INTRODUCTION From transportation of ulk materials in underground and open-ut mines to appliations in power stations and proessing plants to name ut a few, elt onveyors are widely used as an eonomial ontinuous onveying method. Most of these industries rely strongly on onveyor transfers to divert material from one onveyor to another at some point in their system. The unfortunate reality is that some ompanies will e driven y ost minimisation for the short-term rather than planning for the long-term. This will result in ad deisions eing made and the inorret onveyor transfer eing installed with flow-on effets suh as downtime for modifiations or even replaement. There are many faets to the design of a suessful transfer hute inluding minimising produt impat, degradation, hute wear, noise, dust and spillage while maximising the material veloity to allow the produt to leave the hute at or near the speed of the reeiving onveyor. Fully understanding the ehaviour of a material is paramount to designing a suessful transfer hute. Aurately determining the material disharge and trajetory from the head pulley is the first step in this design proess. This paper will fous on the determination of the material trajetory as it leaves the elt onveyor head pulley. This proess will also determine the point at whih the material leaves the elt, referred to as the disharge angle. There are numerous methods availale in the literature fousing on the modelling of material disharge and trajetory, inluding C.E.M.A. (1966; 1979; 1994; 1997; 005), M.H.E.A. (1977; 1986), Booth (1934), Golka (199; 1993), Korzen (1989), Goodyear (1975) and Dunlop (198). These methods have een evaluated for oth horizontal and inlined elts for a range of elt speeds and pulley diameters to evaluate oth low and high-speed onditions. The result of these evaluations have een ompared to allow omment on their ease of use, ompleteness and potential auray. TRAJECTORY METHODS EXAMINED.1 Disharge Angles Within eah of the methods presented, the disharge angle of material leaving the head pulley is determined analytially or graphially whih then ditates the starting position of the material trajetory. As will e evident in setion 3, it was found that there are numerous approahes to determine this disharge angle with some methods omparale to eah other while some eing vastly different. The resulting trajetories generated will vary signifiantly, also detailed in setion 3. Of ourse there should only e one disharge angle for a given material and set of onditions.. Horizontal and Inlined Conveyor Geometries The Conveyor Equipment Manufaturers Assoiation has released six editions of the C.E.M.A. guide, Belt Conveyors for Bulk Materials sine The first five editions follow the same proedure for determination of material trajetory with only slight adjustments to various values in referene tales. In the 6 th Edition of the C.E.M.A. guide (005), there is a hange to how the time interval is alulated for high-speed elts for determining the trajetory profile, now eing alulated ased on the elt speed rather than the tangential veloity. For all other onditions the time intervals are alulated as per the previous editions. The C.E.M.A. method addresses seven onveyor onditions; - low-speed elts where material wraps around the pulley efore disharge (horizontal, inlined and delined elts), 1

3 - high-speed elts where material disharges at the point of tangeny etween the elt and head pulley (horizontal, inlined and delined elts), and - inlined elts where material disharges at the top most vertial point of ontat etween the elt and the pulley (inlined elts). The equations presented in Tale 1 represent the onditions whih need to e satisfied to determine whih of the seven ases aove is to e used. If low-speed onditions exist, equation (7) is used to determine the disharge angle, α d. Horizontal Inlined gr s s gr Low-speed Medium High-speed < 1 (1) N/A s 1 gr () < osα (3) s 1 gr = (4) s 1 gr > (5) s osα < < 1 (6) gr Tale 1 Belt speed onditions for the C.E.M.A. (1966; 1979; 1994; 1997; 005) and M.H.E.A. (1986) methods s gr = osα (7) d The trajetory produed y the C.E.M.A. method is ased on the entroidal path. To plot the trajetory, a tangent line is drawn from the point at whih the material leaves the pulley. At regular time intervals along this tangent line, vertial lines are projeted down to fall distanes supplied in the C.E.M.A. guide. A urve is then drawn through these points to produe the entroid trajetory. The upper and lower trajetory limits an also e plotted y offsetting from the entroidal urve the distane to the elt and to the load height. It is evident from this proedure that a onstant width trajetory path results. C.E.M.A. notes that for light fluffy materials, a high elt speed will alter the upper and lower limits with oth vertial and lateral spread due to air resistane. The Mehanial Handling Engineer s Assoiation guide, Reommended Pratie for Troughed Belt Conveyors (M.H.E.A. 1977) addresses oth low-speed and high-speed elts via entripetal aeleration. The speed onditions of Tale 1 are used, this time using pulley diameter rather than material entroid, resulting in the determination of the lower trajetory. This method (M.H.E.A. 1977) also provides an approximation for the outer trajetory of the material y first determining the angle, α d, at whih the upper surfae of material starts its trajetory. The M.H.E.A. released a seond edition of the guide (M.H.E.A. 1986) whih etter explains the method of determining the material trajetory. On inspetion, it is idential to the C.E.M.A. method up to and inluding the 5 th Edition. However it uses metri rather than imperial units and the onversion results in minor differenes, ultimately ausing minor variations to the trajetory urves whih are determined exatly the same as for the C.E.M.A. method. Booth (1934) found that while using availale theory, a large disrepany was present etween the theory and that of the atual trajetory. After areful investigation and onfirmation of these errors, Booth onluded that, for one, the effets of the material slip were not eing addressed as material disharged over the head pulley. This led to an analytial analysis to develop a more representative theory.

4 Booth s method egins y determining the angle, α, at whih the partile will leave the head pulley, again ased on the low-speed and high-speed onditions desried y the C.E.M.A. method (1966; 1979; 1994; 1997; 005). For low-speed onditions, an initial estimate of the disharge angle is found using equation (8), followed y the angle at whih material slip first ours on the elt, α r, equation (9). The analytial analysis y Booth produes equation (10) and y setting ( ψ ) = and ψ = αr the onstant of integration, C, an e determined. One C has een defined, equations (8) and (10) are solved simultaneously using ( ψ ) = d and ψ = αd to determine the disharge angle and disharge veloity. osα d = (8) gr gr 1 = osαr sinα r (9) μ ( ψ ) ( ) = ( μ + ) μ 1 osψ 3 μ sinψ gr 4 1 μψ + (10) The trajetory is then plotted alulating a time interval from the tangent line and projeting down the freefall distanes provided y Booth. However there is no mention of how the upper trajetory is determined or how the material height at the disharge point is alulated. Booth aknowledged that this method was tedious and ompliated and as an alternative, developed a hart to minimise the time required to analyse a partiular elt onveyor geometry still with a reasonale auray. Golka s method (199; 1993) for determining material trajetory is ased on the Cartesian oordinate system and is for materials without ohesion or adhesion. For low-speed elts, this method one again follows that of C.E.M.A. in determining the disharge angle of the lower trajetory using equation (). It also alulates a separate disharge angle for the upper trajetory y sustituting the radius of the upper trajetory into equation (). The material height on the elt efore disharge is taken diretly from C.E.M.A. however an adjusted material height, h d, is also alulated for the point at whih the upper trajetory disharges from the pulley, equation (11). Ce h d 0.5 h = R p 1+ 1 (11) R p Two divergent oeffiients have een introdued y Golka (199; 1993), ε 1 for the lower and ε for the upper trajetory, whih takes into aount variales suh as air resistane, size distriution, permeaility and partile segregation. Unfortunately no explanation is given as to how they are determined. A set of urves has een produed using the parameters utilised in setion 3 and varying the divergent oeffiients, as an e seen in Figure 1. Not knowing with any ertainty what divergent oeffiients should e used for a given produt/ondition an result in a sustantial variation of the predited trajetory urves. 3

5 ε1 = 0, ε = 0 ε1 = 0.1, ε = 0.1 ε1 = 0., ε = 0. ε1 = 0.3, ε = 0.3 ε1 = 0.4, ε = DIMENSIONS IN MILLIMETERS Figure 1 ariation in trajetories ased on different divergent oeffiients Golka (199; 1993) uses three ases, Tale, to predit the material trajetory from the pulley depending on predefined onditions. The ritial veloity, r, needs to e determined from equation (1) where the radius represents either the lower or upper trajetory. This determines the veloity at whih the transition from low-speed to high-speed ours. 1 = and is determined from equation (13). Condition α d1 α d CASE 1 1 < r1 and < r Use equation (14) Use equation (15) CASE 1 > r 1 point of tangeny point of tangeny CASE 3 1 < r1 and > r Use equation (14) point of tangeny Tale Disharge angle determination for the Golka method (199; 1993) = grosα (1) r 0.5 h = (13) R p α = (14) os 1 d 1 gr p osα d = g Rp hd ( + ) (15) 4

6 The trajetory is then determined from a series of Cartesian oordinate ased equations for pre-seleted time steps. Golka (199; 1993) does not inlude elt thikness when determining the radii of the lower and upper stream. Material height is inluded, though ased on values taken from the early C.E.M.A. guide (1966; 1979). Of all the methods reviewed, Korzen (1989) is the most omplex in its approah, addressing the issues of adhesive materials, inertia, slip and air drag in its alulations. There is also a distintion etween stati frition, μ s, and kinemati frition, μ k, used in the determination of the disharge veloity, d, and disharge angle, α d. For high-speed elts, α d = α and d =. For low-speed elts the angle at whih material egins to slip on the elt efore disharge, α r, is determined from equation (16). To evaluate the onstant of integration, C, the onditions ( ψ ) = and ψ = α r are used in equation (17). To determine the disharge angle, the relationships ( ψ ) = R g osψ and ψ = α d are used. α μ μ r = tan s ± sin sin ( tan s) R g ( k ) ( 1+ 16μk ) k 4 k ( ψ ) = gr Ce μ ψ σ a γ h (16) 4μ 1 osψ 5μ sinψ + (17) The detailed numerial analysis developed y Korzen is ahieved y a series of suessive approximations whih inorporate orreted air drag oeffiients ased on partile shape and a proportionality fator for air drag. Using the X-displaement for any given position of the trajetory, the Y-displaement, trajetory angle and resultant veloity an e otained. The first approximation is for a free falling partile in a vauum, whih is used as the initial trajetory estimate, where air drag is negleted, for all other approximations air drag is applied. The analysis is ontinued until the differential error etween suessive approximations have deviations no greater than 1% or %. One the analysis has een ompleted for a suitale range of X-displaements, the X and Y oordinates are plotted to produe the trajetory for the entral path. The disharge angle and disharge veloity alulated aove are for the entre height of the material stream. Korzen (1989) also allows for alulation of the upper and lower trajetory limit disharge veloities. For high-speed elts the upper and lower disharge veloities are the same as that for the entre height trajetory ut for a low-speed elt, the disharge veloities are alulated ased on ratios of the radius of the lower and upper trajetories. For partiles over 1g in mass, the effets of air drag an e dismissed (Korzen 1989). This eing the ase, only the first approximation is performed and the trajetory is plotted from those values. This will also result in the trajetory urves having a loser math to those generated y other methods. Although the disharge veloity for the lower and upper trajetories an e determined, there is no method desried for the determination of the orresponding trajetories. The assumption has een made that the same approximation method is used, sustituting the alulated disharge veloity for the appropriate trajetory. Korzen (1989) does not inlude the elt thikness when determining the trajetory urves. Although as a perentage of the head pulley diameter the elt thikness is quite small, this omission will result in a minor vertial offset of the plotted trajetories ompared to methods whih inlude elt thikness. There also appears to e several errors in the worked examples presented y Korzen (1989). As a result of one group of errors, additional approximations were required to ahieve an adequate error level. Also, in several quoted equations there also 5

7 appears to e some typographial errors, whih an raise questions as to whether the method is eing applied orretly. The Handook of Conveyor and Elevator Belting (Goodyear 1975) is quite simplisti in its approah to determining material trajetory. Using the priniples of projetile motion, equation (18) and equation (19) are used to determine the X and Y oordinates of the trajetory. The disharge angle is determined for low-speed elt onditions y satisfying one of the ases in Tale 1, this time using the radius of the entral height of the material. X = t (18) gt Y = (19) Goodyear (1975) provides two ases for horizontal onveyor geometries whih are idential to those of equations (1) and (). For low-speed inlined onveyors equation (3) is again used, however for high-speed inlined onveyors equation (0) is used. The Goodyear method determines the disharge angle for the entral material stream, (i.e. h/), however there is no referene to the determination of the disharge angles or trajetories for the lower or upper trajetories. The assumption has een made that all disharge angles are of equal value, thus resulting in parallel trajetory streams. One the X and Y oordinates and the disharge angle have een determined, produing the trajetory urve is straight forward. gr s > osα (0) Belt thikness is not used in the Goodyear method (1975) and as this method determines only the entre trajetory path, a material height must e known. The worked examples y Goodyear (1975) inorporate half the material height into the radius used to determine the entral trajetory path without impliitly providing details of the origin of the material height. Goodyear states that the atual trajetory may e different to that of the one alulated due to other fores ating on the partile stream whih haven t een used in these alulations. The Dunlop Conveyor Manual (198) uses a graphial method to determine the material trajetory leaving the onveyor for low-speed elts and analytial method for high-speed elts. For high-speed elts, material again leaves the elt at the point where the elt is at a tangent to the pulley. To alulate the X oordinate, the distane travelled tangentially from the elt, equation (18) is used and the Y oordinate, the distane material falls elow the tangent line, is alulated using equation (19). Dunlop states that the mathematis ehind the trajetory for low-speed onveyors is omplex and so developed a graphial method. Knowing the elt speed and pulley diameter, it is a simple proess to determine the disharge angle and X, the inremental distane along the disharge tangent line. At eah of these inremental distanes, Y is projeted vertially down with the same values as for high-speed onveyors. If it is found that the desired elt speed does not interset with the required pulley diameter then the method for high-speed elts should e used. For the low-speed ondition, pulley diameters etween 31mm and 1600mm are presented. Outside of this range there is no way to estimate the disharge angle or X. The resulting trajetory is a predition of the lower oundary. There is referene to material depth in the worked examples in the Dunlop Conveyor Manual ut no explanation of how it is determined and there is an indiation that the flow pattern is onvergent..3 Delined Conveyor Belts Only three methods, C.E.M.A. (1966; 1979; 1994; 1997; 005), M.H.E.A. (1986) and Goodyear (1975) allow for the determination of the disharge angle for delined elt onveyors and while Dunlop (198) makes no speifi mention of delined onveyor elts in its guide, it alludes to a delined onveyor elt having the same disharge angle as a 6

8 horizontal or inlined onveyor via a worked example. For this reason delined onveyors have not een inorporated into the omparisons of setion 3, however will riefly e disussed. Up to and inluding the fifth edition of the C.E.M.A. guide (1966; 1979; 1994; 1997) and M.H.E.A (1986) the low-speed or high-speed onditions are determined from equation (1). High-speed disharge is again at the point of tangeny and for low-speed onditions equation (3) is used, followed y equation (7). gr s osα (1) The sixth edition of C.E.M.A. (005) has a variation to the low-speed ondition, now adding the disharge angle to the elt delination angle. However, in the worked example for this speifi ase, the disharge angle alulated is idential to the disharge angle found in the previous five editions of the C.E.M.A. guide whih have the disharge angle determined from the vertial. So this would in atual fat indiate that the graphial representation of this ase is inorret. The Goodyear (1975) method of determining the disharge angle for delined onveyor elts uses equation (0) for high-speed onditions resulting in disharge from the tangeny point and for low-speed uses equation (3) and equation (7) and is plotted from the vertial as with the earlier C.E.M.A. methods and M.H.E.A..4 Critial Belt Speed The ritial elt speed, touhed on previously, refers to the point of transition from low-speed to high-speed onditions and Tale 3 summarises the five unique equations used. As previously explained, Golka (199; 1993) determines two distint disharge angles and following from this there are also two ritial elt speeds. The method y Korzen (1989) inorporates an adhesive stress omponent, however when the adhesive stress equals zero, the equation is idential to that of Goodyear (1975). For the Dunlop method (198), the ritial elt speed is graphially determined for the lower and upper disharge angles. There are however limitations to it s determination due to the maximum pulley diameter harted eing 1600mm. C.E.M.A. and M.H.E.A. r, gr osα = () R Booth and Golka (low), = gr osα (3) Golka (high), r h = gr 1+ osα (4) R r σ a Korzen r, = gr osα + γ h (5) Goodyear, = gr osα (6) r Tale 3 Critial elt speeds for the trajetory methods 7

9 .5 Effet of Transition Geometry As the onveyor elt moves through the transition zone, stresses are generated at the elt edges and an ause permanent deformation if the elasti limit is surpassed. To minimise or even overome this, the head pulley an e raised to ½ the trough depth thus reduing the stresses on the elt (C.E.M.A. 005). Although this is stated y C.E.M.A. (005), there is no indiation whether this is applied in the determination of onveyor trajetories. Roerts, Wihe et al. (004) state that the transition angle, ε, does effet elt disharge trajetories and as an e seen in Figure, the point of tangeny of the elt to the head pulley eomes α + ε, where ε is determined from equation (7). 1 z ε = tan L (7) t For low-speed onditions, the angle θ 0 is determined, defining the angle that disharge ommenes. For high-speed onditions disharge is at the point of tangeny (i.e. θ 0 = 0º). Figure Transition Geometry (Roerts, Wihe et al. 004) The trajetory oordinates an then e determined using equation (8). X ( α ε θ ) = Y tan os gy ( α + ε θ ) 0 (8) The method presented y Roerts et al. (004) is the only approah whih diretly addresses the transition arrangement when onsidering the determination of the trajetory and even though a half trough transition is reommended y C.E.M.A. (005), there is no indiation it is applied to the trajetory method provided. 3 TRAJECTORY METHODS COMPARED To allow for a diret omparison etween all the methods desried, a onstant set of parameters was seleted as shown in Tale 4. The only variations allowed for are elt speed ( = 1.5, 3.0 and 6.0ms -1 ), pulley diameter (D p = 0.5, 1.0 and 1.5m) and elt inlination angle (0 30 ), all others are fixed. As adhesive stress is only aounted for y the Korzen (1989) method, it has een set to zero to allow omparison of all methods. Although not quantified here, if adhesive stress were to e applied, the result would e an inreased disharge angle whih would result in a trajetory stream falling loser to the head pulley ompared with the other methods. 8

10 PARAMETERS ALUES Belt Width, w 0.76 m Belt Thikness, 0.01 m Belt Speed, 1.5 ms -1 Pulley Diameter, D p 0.5 m Surharge Angle 0 Troughing Angle 0 Belt Inlination Angle, α 0 Divergent Coeffiient, Lower, ε Divergent Coeffiient, Upper, ε Coeffiient of Frition, Stati, μ s Coeffiient of Frition, Kinemati, μ k Adhesion Stress, σ a 0 kpa Equivalent Spherial Partile Diameter m Atmospheri Temperature, T atm 0 C Atmospheri Pressure, P atm 101 kpa Air isosity, η f 1.80E-05 Nsm - Produt Density, ρ 000 kgm -3 Tale 4 Conveyor parameters used for omparisons For the omparisons made, a numer of graphs have een seleted to represent the oserved differenes for oth low-speed and high-speed elt onditions, all for D p = 1.0m. Where trajetory urves for different models vary y no more than % for oth the lower and upper trajetories, at a vertial drop height of 4m, they have een grouped together. As mentioned in setion, some methods do not inlude elt thikness when determining the radius for the lower, entre or upper trajetory radii and as suh minor modifiations to those methods has een implemented y adding the elt thikness, allowing for diret omparisons to e ahieved. Of further note is the radius for the upper trajetory streams varies for some methods regardless of the inlusion of elt thikness due to the method used to determine the height of the material stream. 3.1 Disharge Angles A set of aritrary parameters has een used to ompare the disharge angles for the trajetory methods desried and a seleted group is shown in Figures 3 to 8. Using a onstant pulley diameter and varying the elt veloity only, the following oservations have een made; a) Figure 3 shows two distint groupings of disharge angles with slight variations evident in oth groups ut as the elt veloity inreases, there is a notieale spread, moving away from the initial two groupings, see Figures 3 to 7, ) Figure 7 shows that some methods, C.E.M.A. (1966; 1979; 1994; 1997; 005), M.H.E.A. (1986) and Golka (upper) (199; 1993) have already reahed high-speed onditions, ) Golka (upper) (199; 1993) reahes high-speed onditions ased on the alulated tangential veloity of the upper stream, whereas Golka (lower) (199; 1993) is still under low-speed onditions, refer to Tale 3, d) Figure 8 shows that all methods are now under high-speed onditions, noting that all ritial elt speeds from Tale 4 are elow the given elt speed, ( =1.75ms -1 ), 9

11 CEMA 1,,4,5 BOOTH CEMA 6 GOLKA MHEA86 KORZEN DIMENSIONS IN METRES DUNLOP GOODYEAR Figure 3 Figure 4 Figure 5 D p =0.5m, =0.50ms -1 D p =0.5m, =0.75ms -1 D p =0.5m, =1.00ms Figure 6 Figure 7 Figure 8 D p =0.5m, =1.5ms -1 D p =0.5m, =1.50ms -1 D p =0.5m, =1.75ms -1 e) referring to the two disharge angles for Dunlop (198) in Figure 7, there is an indiation that there is a high onvergene of the lower and upper paths ased on the authors assumption that two distint disharge angles should e determined. There will in all likelihood e a rossing of the streams whih in reality would not our. Low-speed onditions are ahievale over a wider range of elt speeds as the pulley diameter inreases, for example a pulley diameter of m will operate under low-speed onditions up to approximately 3.1ms Low-Speed Trajetories A omparison graph for the low-speed trajetory preditions is presented in Figure 9. Although having a larger pulley diameter than the example of Figure 3, Figure 9 learly shows two main groupings of methods for the produed trajetories. It is also lear that these two groupings have an approximate horizontal variation of 500mm at a drop height of 4000mm 10

12 whih will have a marked effet on the design of a transfer hute dependant on whih trajetory method is applied. Although the figures are not presented in this paper, a further trend found that as pulley diameter inreases for a given elt speed, the trajetories for the models fall loser to the head pulley. The methods of Booth (1934) and Dunlop (198) produe a near idential urve regardless of the fat that one uses a highly analytial approah and the other uses a straightforward graphial approah respetively CEMA1,,4,5,6,MHEA1986 BOOTH,DUNLOP GOODYEAR GOLKA,GOLKA (divergent) KORZEN KORZEN (air drag) DIMENSIONS IN MILLIMETERS Figure 9 Low-speed ondition, horizontal onveyor, lower and upper trajetory path Pulley diameter, D p = 1.0m, elt veloity, = 1.5ms High-Speed Trajetories In order to investigate the hange in trajetory profile due to elt speed, two high-speed elt onditions have een seleted, = 3.0 ms -1 and = 6.0 ms -1. In Figure 10 and Figure 11 the inlusion of air drag y Korzen (1989) results in a trajetory predition learly muh shallower than all other methods. If air drag is negleted, the resulting trajetory predition is loated amongst the other methods. In setion it was explained that Golka (199; 1993) uses divergent oeffiients to otain a etter approximation of the trajetory paths. Without explanation of how these divergent oeffiients have een determined, it is hard to justify using this method with any level of auray. If on the other hand the divergent oeffiients are negleted in the alulations, the resulting trajetories are idential to the Korzen (1989) method when air drag is negleted, as the equations are idential. In oth ases shown, the early C.E.M.A. (1966; 1979; 1994; 1997) and M.H.E.A. (1986) methods generate the highest trajetory urve and as the disharge veloity inreases, the variation from the other urves eomes more defined. As with the low-speed onditions, the Booth (1934) and Dunlop (198) methods again produe near idential trajetory urves. The inorporation of adhesive stress y Korzen (1989) has een previously explained ut it is lear in Figure 1 and Figure 11 that if adhesive 11

13 stress were inluded in the alulations, the resulting trajetory urves would deviate even further from the trajetory urves of the other methods CEMA1,,4,5,MHEA1986 BOOTH,DUNLOP CEMA6,GOODYEAR GOLKA,KORZEN GOLKA(divergent) KORZEN (air drag) DIMENSIONS IN MILLIMETERS Figure 10 High-speed ondition, horizontal onveyor, lower and upper trajetory path Pulley diameter, D p = 1.0m, elt veloity, = 3.0ms DIMENSIONS IN MILLIMETERS CEMA1,,4,5,MHEA1986 CEMA6,BOOTH,DUNLOP,GOODYEAR,GOLKA,KORZEN GOLKA(divergent) KORZEN (air drag) Figure 11 High-speed ondition, horizontal onveyor, lower and upper trajetory path Pulley diameter, D p = 1.0m, elt veloity, = 6.0ms -1 4 ALIDATION 4.1 Experimental A unique onveyor transfer researh faility, Figure 1, will soon e ommissioned at the University of Wollongong whih will allow oth trajetories and variale geometry transfer 1

14 hutes to e experimentally investigated. It is the intention to generate a sustantial dataase of trajetory information from this faility addressing as many of the parameters inorporated into the trajetory methods presented. This, it is hoped, will allow diret validation of the trajetory methods and ultimately determine whih of the methods is most aurate. Figure 1 Experimental onveyor transfer researh faility at University of Wollongong 4. Simulation The distint element method (DEM) is an ideal simulation tool for approahing this type of granular flow prolem. There are a range of a ommerial pakages availale whih will allow trajetory profiles to e simulated. This will produe another method of omparing the trajetory methods presented to the dataase of experimental results. It should e highlighted however, that there is proaly no option for seleting how a trajetory is determined, solely eing ased on the physis of the ommerial DEM odes (whih proaly do not allow for differenes in moisture, ohesion, et). The simulations may or may not orrespond to the trajetory methods presented and will e evaluated as required. 5 CONCLUSION This paper has presented a numer of the more widely used and/or readily availale trajetory predition methods pulished in the literature. They range in omplexity from the asi, Dunlop (198) and Goodyear (1975) to the omplex, Booth (1934), Golka (199; 1993) and Korzen (1989). Some methods inlude a multitude of parameters suh as C.E.M.A. (1966; 1979; 1994; 1997; 005) and M.H.E.A. (1986) while others inorporate additional parameters suh as divergent oeffiients (Golka 199; 1993) and air drag (Korzen 1989). For low-speed elt onditions there appear to e two distint groupings of trajetory preditions, see Figure 8, when the elt speed is relatively low. However as the elt speed inreases, so does the satter of disharge angles as seen in Figure 3 to Figure 8. 13

15 For high-speed elt onditions an oservation made was that some of the asi methods approximated a trajetory whih was also predited y the more omplex methods, see Figure 10 and Figure 11. Inorporating air drag into the trajetory preditions (Korzen 1989) auses a dramati variation from the other trajetory methods. When onsidering air drag, it makes sense that material will drop away muh quiker than the other methods ut whether it is truly representative of an atual trajetory needs to e explored further. It is fair to say that the differing approahes to determining the material trajetory result in some signifiant differenes whih annot all e orret. Based solely on the omparisons made in this paper, it is impossile to say with any ertainty whih of the methods will produe the most aurate predition of an atual trajetory. It may even e the ase that one method predits aurately the low-speed onditions whereas another is etter suited for high-speed onditions. Further researh is definitely needed and with the aid of a unique onveyor transfer researh faility and simulation via DEM at the University of Wollongong, this should provide a means to apply improvements to the availale trajetory methods 6 ACKNOWLEDGEMENTS The authors wish to aknowledge the support of the Australian Researh Counil, Rio Tinto OTX and Rio Tinto Iron Ore Expansion Projets for their finanial and in-kind ontriutions to the Linkage Projet whih allows this researh to e pursued. 7 NOMENCLATURE elt thikness m C onstant of integration - D p pulley diameter m g gravity ms - h material depth m h d material depth at disharge m Lt transition length m P atm atmospheri pressure kpa R aritrary radius m R radius to outer elt surfae m R radius of material entroid/entre m R p pulley radius m t inrement time for trajetory path s T atm atmospheri temperature C 1 disharge veloity of lower oundary ms -1 disharge veloity of upper oundary ms -1 elt veloity ms -1 r ritial veloity ms -1 d veloity of material at disharge point ms -1 s tangential veloity of material at disharge point ms -1 w elt width m z height of the transition m α initial material disharge angle measured from the vertial α elt inlination angle α d material disharge angle measured from the vertial trajetory α d1 material disharge angle measured from the vertial for lower trajetory α d material disharge angle measured from the vertial for upper trajetory α r angle at whih slip egins to our γ speifi gravity of ulk solid knm -3 ε transition angle ε 1 divergent oeffiient - ε divergent oeffiient - η f air visosity Nsm - θ 0 angle at whih disharge ommenes μ oeffiient of frition - 14

16 μ k oeffiient of kinemati frition - μ s oeffiient of stati frition - σ a adhesion stress kpa 8 REFERENCES Booth, E. P. O. (1934). "Trajetories from Conveyors - Method of Calulating Them Correted". Engineering and Mining Journal, ol. 135, No. 1, Deemer, pp C.E.M.A. (1966). Belt Conveyors for Bulk Materials. 1st Ed, Conveyor Equipment Manufaturers Assoiation. C.E.M.A. (1979). Belt Conveyors for Bulk Materials. nd Ed, Conveyor Equipment Manufaturers Assoiation. C.E.M.A. (1994). Belt Conveyors for Bulk Materials. 4th Ed, Conveyor Equipment Manufaturers Assoiation. C.E.M.A. (1997). Belt Conveyors for Bulk Materials. 5th Ed, Conveyor Equipment Manufaturers Assoiation. C.E.M.A. (005). Belt Conveyors for Bulk Materials. 6th Ed, Conveyor Equipment Manufaturers Assoiation. Dunlop (198). Dunlop Industrial Conveyor Manual. Golka, K. (199). "Disharge Trajetories of Bulk Solids". 4th International Conferene on Bulk Materials Storage, Handling and Transportation, Wollongong, NSW, Australia, 6th - 8th July, pp Golka, K. (1993). "Predition of the Disharge Trajetories of Bulk Materials". Bulk Solids Handling, ol. 13, No. 4, Novemer, pp Goodyear (1975). Goodyear Handook of Conveyor & Elevator Belting: Setion 11. Korzen, Z. (1989). "Mehanis of Belt Conveyor Disharge Proess as Affeted y Air Drag". Bulk Solids Handling, ol. 9, No. 3, August, pp M.H.E.A. (1977). Reommended Pratie for Troughed Belt Conveyors, Mehanial Handling Engineer s Assoiation. M.H.E.A. (1986). Reommended Pratie for Troughed Belt Conveyors, Mehanial Handling Engineer s Assoiation. Roerts, A. W., Wihe, S. J., Ili, D. D. and Plint, S. R. (004). "Flow Dynamis and Wear Considerations in Transfer Chute Design". 8th International Conferene on Bulk Materials Storage, Handling and Transportation, Wollongong, NSW, Australia, 5th - 8th July, Institution of Engineers, Australia, pp

17 9 AUTHORS C Mr. David Hastie B.E. (Mehanial), M.E. (Honours) has een employed with the Centre for Bulk Solids and Partiulate Tehnologies at the University of Wollongong, Australia sine He is a Memer of the Institution of Engineers Australia and a memer of the Australian Soiety for Bulk Solids Handling. In July of 005 he took up an Aademi Researh Fellow position within the Faulty of Engineering. Key areas of interest inlude pneumati onveying, rotary valve air leakage, partile segregation and onveyor transfers and trajetories and has extensive experiene in experimental investigations, instrumentation, data aquisition and analysis, omputer programming and digital video imaging and proessing. He was the main researher on the suessful IFPRI, In. funded projet, on the predition of transport oundaries in fluidised dense-phase and low-veloity pneumati onveying and urrently he is Researh Assoiate on a three year projet titled Quantifiation And Modelling Of Partile Flow Mehanisms In Conveyor Transfers. He is also a part-time PhD andidate working on transfer hute quantifiation and modelling. Peter Wypyh B.E. (Mehanial, Honours 1), PhD is the Diretor of the ARC endorsed Key Centre for Bulk Solids and Partiulate Tehnologies at the University of Wollongong. He has een involved with the researh and development of solids handling and proessing tehnology sine Peter Wypyh has pulished over 80 artiles. He is urrently Chair of the Australian Soiety for Bulk Solids Handling. Peter Wypyh is also the General Manager of Bulk Materials Engineering Australia and has ompleted over 450 industrial projets, involving R&D of new tehnologies, feasiility studies, trouleshooting, general/onept design, optimisation, deottleneking, safety/hazard audits and/or rationalisation of plants and proesses for ompanies all around Australia and in the USA, Hong Kong, New Zealand, China, Singapore and Korea. 10 AUTHOR CONTACTS Mr David Hastie Centre for Bulk Solids and Partiulate Tehnologies Faulty of Engineering University of Wollongong Northfields Avenue, Wollongong, NSW, 5, Australia dhastie@uow.edu.au A/Prof Peter Wypyh Centre for Bulk Solids and Partiulate Tehnologies Faulty of Engineering University of Wollongong Northfields Avenue, Wollongong, NSW, 5, Australia wypyh@uow.edu.au 16

Conveyor trajectory discharge angles

Conveyor trajectory discharge angles University of Wollongong Researh Online Faulty of Engineering - Papers (Arhive) Faulty of Engineering and Information Sienes 007 Conveyor trajetory disharge angles David B. Hastie University of Wollongong,

More information

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six

More information

Buckling loads of columns of regular polygon cross-section with constant volume and clamped ends

Buckling loads of columns of regular polygon cross-section with constant volume and clamped ends 76 Bukling loads of olumns of regular polygon ross-setion with onstant volume and lamped ends Byoung Koo Lee Dept. of Civil Engineering, Wonkwang University, Iksan, Junuk, 7-79, Korea Email: kleest@wonkwang.a.kr

More information

UTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.

UTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker. UTC Engineering 329 Proportional Controller Design for Speed System By John Beverly Green Team John Beverly Keith Skiles John Barker 24 Mar 2006 Introdution This experiment is intended test the variable

More information

Fiber Optic Cable Transmission Losses with Perturbation Effects

Fiber Optic Cable Transmission Losses with Perturbation Effects Fiber Opti Cable Transmission Losses with Perturbation Effets Kampanat Namngam 1*, Preeha Yupapin 2 and Pakkinee Chitsakul 1 1 Department of Mathematis and Computer Siene, Faulty of Siene, King Mongkut

More information

Verka Prolović Chair of Civil Engineering Geotechnics, Faculty of Civil Engineering and Architecture, Niš, R. Serbia

Verka Prolović Chair of Civil Engineering Geotechnics, Faculty of Civil Engineering and Architecture, Niš, R. Serbia 3 r d International Conferene on New Developments in Soil Mehanis and Geotehnial Engineering, 8-30 June 01, Near East University, Niosia, North Cyprus Values of of partial fators for for EC EC 7 7 slope

More information

EFFECTS OF COUPLE STRESSES ON PURE SQUEEZE EHL MOTION OF CIRCULAR CONTACTS

EFFECTS OF COUPLE STRESSES ON PURE SQUEEZE EHL MOTION OF CIRCULAR CONTACTS -Tehnial Note- EFFECTS OF COUPLE STRESSES ON PURE SQUEEZE EHL MOTION OF CIRCULAR CONTACTS H.-M. Chu * W.-L. Li ** Department of Mehanial Engineering Yung-Ta Institute of Tehnology & Commere Ping-Tung,

More information

Simplified Buckling Analysis of Skeletal Structures

Simplified Buckling Analysis of Skeletal Structures Simplified Bukling Analysis of Skeletal Strutures B.A. Izzuddin 1 ABSRAC A simplified approah is proposed for bukling analysis of skeletal strutures, whih employs a rotational spring analogy for the formulation

More information

Calculation of Magnetic Field of Two Coaxially Located Circular Permanent Magnets Using a Method of Physical Analogies

Calculation of Magnetic Field of Two Coaxially Located Circular Permanent Magnets Using a Method of Physical Analogies World Applied Sienes Journal 3 (): 345-35, 3 ISSN 88-495 IDOSI Puliations, 3 DOI: 589/idosiwasj3335 Calulation of Magneti Field of Two Coaxially Loated Cirular Permanent Magnets Using a Method of Physial

More information

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install

More information

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field Four-dimensional equation of motion for visous ompressible substane with regard to the aeleration field, pressure field and dissipation field Sergey G. Fedosin PO box 6488, Sviazeva str. -79, Perm, Russia

More information

MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP OF CONCRETE IN UNIAXIAL COMPRESSION

MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP OF CONCRETE IN UNIAXIAL COMPRESSION VIII International Conferene on Frature Mehanis of Conrete and Conrete Strutures FraMCoS-8 J.G.M. Van Mier, G. Ruiz, C. Andrade, R.C. Yu and X.X. Zhang Eds) MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP

More information

Ordering Generalized Trapezoidal Fuzzy Numbers. Using Orthocentre of Centroids

Ordering Generalized Trapezoidal Fuzzy Numbers. Using Orthocentre of Centroids International Journal of lgera Vol. 6 no. 69-85 Ordering Generalized Trapezoidal Fuzzy Numers Using Orthoentre of Centroids Y. L. P. Thorani P. Phani Bushan ao and N. avi Shanar Dept. of pplied Mathematis

More information

Compression Members Local Buckling and Section Classification

Compression Members Local Buckling and Section Classification Compression Memers Loal Bukling and Setion Classifiation Summary: Strutural setions may e onsidered as an assemly of individual plate elements. Plate elements may e internal (e.g. the wes of open eams

More information

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple

More information

Analysis of Leakage Paths Induced by Longitudinal Differential Settlement of the Shield-driven Tunneling

Analysis of Leakage Paths Induced by Longitudinal Differential Settlement of the Shield-driven Tunneling 2016 rd International Conferene on Engineering Tehnology and Appliation (ICETA 2016) ISBN: 978-1-60595-8-0 Analysis of Leakage Paths Indued by Longitudinal Differential Settlement of the Shield-driven

More information

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 4, 2012

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 4, 2012 INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume, No 4, 01 Copyright 010 All rights reserved Integrated Publishing servies Researh artile ISSN 0976 4399 Strutural Modelling of Stability

More information

The Laws of Acceleration

The Laws of Acceleration The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the

More information

On the Buckling of Axially Restrained Steel Columns in Fire. University of Bath, Bath BA2 7AY, UK. University of Sheffield, Sheffield S1 3JD, UK

On the Buckling of Axially Restrained Steel Columns in Fire. University of Bath, Bath BA2 7AY, UK. University of Sheffield, Sheffield S1 3JD, UK Title On the Bukling of Axially Restrained Steel Columns in Fire Authors Affiliations P.G. Shepherd a * and I.W. Burgess a Department of Arhiteture & Civil Engineering University of Bath, Bath BA2 7AY,

More information

10.2 The Occurrence of Critical Flow; Controls

10.2 The Occurrence of Critical Flow; Controls 10. The Ourrene of Critial Flow; Controls In addition to the type of problem in whih both q and E are initially presribed; there is a problem whih is of pratial interest: Given a value of q, what fators

More information

The Complete Energy Translations in the Detailed. Decay Process of Baryonic Sub-Atomic Particles. P.G.Bass.

The Complete Energy Translations in the Detailed. Decay Process of Baryonic Sub-Atomic Particles. P.G.Bass. The Complete Energy Translations in the Detailed Deay Proess of Baryoni Su-Atomi Partiles. [4] P.G.Bass. PGBass P12 Version 1..3 www.relativitydomains.om August 218 Astrat. This is the final paper on the

More information

A Spatiotemporal Approach to Passive Sound Source Localization

A Spatiotemporal Approach to Passive Sound Source Localization A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,

More information

Examining Applied Rational Functions

Examining Applied Rational Functions HiMAP Pull-Out Setion: Summer 1990 Eamining Applied Rational Funtions Flod Vest Referenes Environmental Protetion Agen. Gas Mileage Guide. (Copies an usuall e otained from a loal new ar dealer.) Information

More information

Process engineers are often faced with the task of

Process engineers are often faced with the task of Fluids and Solids Handling Eliminate Iteration from Flow Problems John D. Barry Middough, In. This artile introdues a novel approah to solving flow and pipe-sizing problems based on two new dimensionless

More information

Analysis of discretization in the direct simulation Monte Carlo

Analysis of discretization in the direct simulation Monte Carlo PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,

More information

To investigate the relationship between the work done to accelerate a trolley and the energy stored in the moving trolley.

To investigate the relationship between the work done to accelerate a trolley and the energy stored in the moving trolley. SP2h.1 Aelerating trolleys Your teaher may wath to see if you an follow instrutions safely take areful measurements. Introdution The work done y a fore is a measure of the energy transferred when a fore

More information

Relativity in Classical Physics

Relativity in Classical Physics Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of

More information

WRAP-AROUND GUSSET PLATES

WRAP-AROUND GUSSET PLATES WRAP-AROUND GUSSET PLATES Where a horizontal brae is loated at a beam-to-olumn intersetion, the gusset plate must be ut out around the olumn as shown in Figure. These are alled wrap-around gusset plates.

More information

Experimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading

Experimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading 2th International LS-DYNA Users Conferene Constitutive Modeling(2) xperimental Investigation and F Analysis of Fiber Woven Layered Composites under Dynami Loading Pavel A. Mossakovsky, Fedor K. Antonov,

More information

The Second Postulate of Euclid and the Hyperbolic Geometry

The Second Postulate of Euclid and the Hyperbolic Geometry 1 The Seond Postulate of Eulid and the Hyperboli Geometry Yuriy N. Zayko Department of Applied Informatis, Faulty of Publi Administration, Russian Presidential Aademy of National Eonomy and Publi Administration,

More information

UNIVERSAL RELATIONSHIP BETWEEN COLLECTION EFFICIENCY AND THE CORONA POWER OF THE ELECTROSTATIC PRECIPITATOR

UNIVERSAL RELATIONSHIP BETWEEN COLLECTION EFFICIENCY AND THE CORONA POWER OF THE ELECTROSTATIC PRECIPITATOR Australia 006 Paper 5B UNIVERSAL RELATIONSHIP BETWEEN COLLECTION EFFICIENCY AND THE CORONA POWER OF THE ELECTROSTATIC PRECIPITATOR YAKOV S. KHODORKOVSKY & MICHAEL R. BELTRAN Beltran, In., U.S.A. ABSTRACT

More information

Effect of Rotation, Magnetic Field and Initial Stresses on Propagation of Plane Waves in Transversely Isotropic Dissipative Half Space

Effect of Rotation, Magnetic Field and Initial Stresses on Propagation of Plane Waves in Transversely Isotropic Dissipative Half Space Applied Mathematis 4 7- http://dx.doi.org/.436/am..48a5 Published Online August (http://www.sirp.org/journal/am) Effet of otation Magneti Field and Initial Stresses on Propagation of Plane Waves in Transversely

More information

2. The Energy Principle in Open Channel Flows

2. The Energy Principle in Open Channel Flows . The Energy Priniple in Open Channel Flows. Basi Energy Equation In the one-dimensional analysis of steady open-hannel flow, the energy equation in the form of Bernoulli equation is used. Aording to this

More information

Chapter 3 Lecture 7. Drag polar 2. Topics. Chapter-3

Chapter 3 Lecture 7. Drag polar 2. Topics. Chapter-3 hapter 3 eture 7 Drag polar Topis 3..3 Summary of lift oeffiient, drag oeffiient, pithing moment oeffiient, entre of pressure and aerodynami entre of an airfoil 3..4 Examples of pressure oeffiient distributions

More information

Natural Convection Experiment Measurements from a Vertical Surface

Natural Convection Experiment Measurements from a Vertical Surface OBJECTIVE Natural Convetion Experiment Measurements from a Vertial Surfae 1. To demonstrate te basi priniples of natural onvetion eat transfer inluding determination of te onvetive eat transfer oeffiient.

More information

Click here to order this book in one of two formats: softcover ISBN: $50.00 ISBN: $50.00

Click here to order this book in one of two formats: softcover ISBN: $50.00 ISBN: $50.00 ere is a sample hapter from Letures on Radiation Dosimetry Physis: A Deeper Look into the Foundations of Clinial Protools his sample hapter is opyrighted and made availale for personal use only No part

More information

A Differential Equation for Specific Catchment Area

A Differential Equation for Specific Catchment Area Proeedings of Geomorphometry 2009. Zurih, Sitzerland, 3 ugust - 2 September, 2009 Differential Equation for Speifi Cathment rea J. C. Gallant, M. F. Huthinson 2 CSIRO Land and Water, GPO Box 666, Canberra

More information

STUDY OF THE CAPABILITIES OF AVO-METHODS FOR THE DETECTION OF HYDROCARBON ACCUMULATIONS

STUDY OF THE CAPABILITIES OF AVO-METHODS FOR THE DETECTION OF HYDROCARBON ACCUMULATIONS STUDY OF THE CAPABILITIES OF AVO-METHODS FOR THE DETECTION OF HYDROCARBON ACCUMULATIONS Martin Toshev University of Mining and Geology St. Ivan Rilski, 1700 Sofia, E-mail: martin.toshev86@gmail.om ABSTARCT.

More information

COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION

COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION 4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and

More information

1. INTRODUCTION. l t t r. h t h w. t f t w. h p h s. d b D F. b b d c. L D s

1. INTRODUCTION. l t t r. h t h w. t f t w. h p h s. d b D F. b b d c. L D s Rapid Assessment of Seismi Safety of Elevated ater Tanks with FRAME Staging 1. NTRODUCTON 1.1 ntrodution ater tanks are lifeline items in the aftermath of earthquakes. The urrent pratie of designing elevated

More information

Stability Analysis of Orbital Motions around Uniformly Rotating Irregular Asteroids

Stability Analysis of Orbital Motions around Uniformly Rotating Irregular Asteroids Stability Analysis of Orbital Motions around Uniformly Rotating Irregular Asteroids By Xiyun HoU, ) Daniel J. SCHEERES, ) Xiaosheng XIN, ) Jinglang FENG, ) Jingshi TANG, ) Lin LIU, ) ) Shool of Astronomy

More information

The gravitational phenomena without the curved spacetime

The gravitational phenomena without the curved spacetime The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,

More information

IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION

IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION 09-1289 Citation: Brilon, W. (2009): Impedane Effets of Left Turners from the Major Street at A TWSC Intersetion. Transportation Researh Reord Nr. 2130, pp. 2-8 IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE

More information

CHAPTER 26 The Special Theory of Relativity

CHAPTER 26 The Special Theory of Relativity CHAPTER 6 The Speial Theory of Relativity Units Galilean-Newtonian Relativity Postulates of the Speial Theory of Relativity Simultaneity Time Dilation and the Twin Paradox Length Contration Four-Dimensional

More information

Theory of Thermodynamic Variables of Rubber Band Heat Engine

Theory of Thermodynamic Variables of Rubber Band Heat Engine Journal of Physis: onferene Series PAPER OPEN AESS Theory of Thermodynami Variales of Ruer Band Heat Engine To ite this artile: Nurhidayah Muharayu et al 016 J. Phys.: onf. Ser. 739 01091 Related ontent

More information

Motor Sizing Application Note

Motor Sizing Application Note PAE-TILOGY Linear Motors 70 Mill orest d. Webster, TX 77598 (8) 6-7750 ax (8) 6-7760 www.trilogysystems.om E-mail emn_support_trilogy@parker.om Motor Sizing Appliation Note By Jak Marsh Introdution Linear

More information

A Heuristic Approach for Design and Calculation of Pressure Distribution over Naca 4 Digit Airfoil

A Heuristic Approach for Design and Calculation of Pressure Distribution over Naca 4 Digit Airfoil IOSR Journal of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 PP 11-15 www.iosrjen.org A Heuristi Approah for Design and Calulation of Pressure Distribution over Naa 4 Digit Airfoil G.

More information

University of Bristol - Explore Bristol Research

University of Bristol - Explore Bristol Research Chen, J., Ravey, E., Hallett, S. R., Wisnom, M. R., & Grassi, M. (2009). Predition of delamination in braided omposite T-piee speimens. Composites Siene and Tehnology, 69(14), 2363-2367. 10.1016/j.ompsiteh.2009.01.027

More information

Probabilistic analysis of the face stability of circular tunnels

Probabilistic analysis of the face stability of circular tunnels Probabilisti analysis of the fae stability of irular tunnels Guilhem Mollon 1, Daniel Dias 2 and Abdul-Hamid Soubra 3, M. ASCE 1 PhD Student, INSA Lyon, LGCIE Site Coulomb 3, Géotehnique, Bât. J.C.A. Coulomb,

More information

Developing Excel Macros for Solving Heat Diffusion Problems

Developing Excel Macros for Solving Heat Diffusion Problems Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper

More information

THE REFRACTION OF LIGHT IN STATIONARY AND MOVING REFRACTIVE MEDIA

THE REFRACTION OF LIGHT IN STATIONARY AND MOVING REFRACTIVE MEDIA HDRONIC JOURNL 24, 113-129 (2001) THE REFRCTION OF LIGHT IN STTIONRY ND MOVING REFRCTIVE MEDI C. K. Thornhill 39 Crofton Road Orpington, Kent, BR6 8E United Kingdom Reeived Deember 10, 2000 Revised: Marh

More information

Critical Reflections on the Hafele and Keating Experiment

Critical Reflections on the Hafele and Keating Experiment Critial Refletions on the Hafele and Keating Experiment W.Nawrot In 1971 Hafele and Keating performed their famous experiment whih onfirmed the time dilation predited by SRT by use of marosopi loks. As

More information

Development of a user element in ABAQUS for modelling of cohesive laws in composite structures

Development of a user element in ABAQUS for modelling of cohesive laws in composite structures Downloaded from orbit.dtu.dk on: Jan 19, 2019 Development of a user element in ABAQUS for modelling of ohesive laws in omposite strutures Feih, Stefanie Publiation date: 2006 Doument Version Publisher's

More information

THE EFFECT OF CONSOLIDATION RATIOS ON DYNAMIC SHEAR MODULUS OF SOIL

THE EFFECT OF CONSOLIDATION RATIOS ON DYNAMIC SHEAR MODULUS OF SOIL Otober 12-17, 28, Beijing, China THE EFFECT OF CONSOLIDATION RATIOS ON DYNAMIC SHEAR MODULUS OF SOIL J. Sun 1 and X.M. Yuan 2 1 Assoiate Professor, Institute of Civil Engineering, Heilongjiang University,

More information

Improving the worthiness of the Elder problem as a benchmark for buoyancy driven convection models

Improving the worthiness of the Elder problem as a benchmark for buoyancy driven convection models Improving the worthiness of the Elder problem as a benhmark for buoyany driven onvetion models Maarten van Reeuwijk, Simon A. Mathias The Department of Civil and Environmental Engineering Imperial College

More information

Optimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach

Optimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)

More information

Relativistic Dynamics

Relativistic Dynamics Chapter 7 Relativisti Dynamis 7.1 General Priniples of Dynamis 7.2 Relativisti Ation As stated in Setion A.2, all of dynamis is derived from the priniple of least ation. Thus it is our hore to find a suitable

More information

CALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS

CALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS International Journal of Modern Physis A Vol. 24, No. 5 (2009) 974 986 World Sientifi Publishing Company CALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS PAVEL SNOPOK, MARTIN

More information

Lesson 23: The Defining Equation of a Line

Lesson 23: The Defining Equation of a Line Student Outomes Students know that two equations in the form of ax + y = and a x + y = graph as the same line when a = = and at least one of a or is nonzero. a Students know that the graph of a linear

More information

Edexcel GCSE Maths Foundation Skills Book Ratio, proportion and rates of change 1

Edexcel GCSE Maths Foundation Skills Book Ratio, proportion and rates of change 1 Guidane on the use of odes for this mark sheme ethod mark A C P ao oe ft Auray mark ark awarded independent of method Communiation mark Proof, proess or justifiation mark Corret answer only Or equivalent

More information

23.1 Tuning controllers, in the large view Quoting from Section 16.7:

23.1 Tuning controllers, in the large view Quoting from Section 16.7: Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output

More information

THEORETICAL PROBLEM No. 3 WHY ARE STARS SO LARGE?

THEORETICAL PROBLEM No. 3 WHY ARE STARS SO LARGE? THEORETICAL PROBLEM No. 3 WHY ARE STARS SO LARGE? The stars are spheres of hot gas. Most of them shine beause they are fusing hydrogen into helium in their entral parts. In this problem we use onepts of

More information

IN-PLANE VIBRATIONS OF CURVED BEAMS WITH VARIABLE CROSS-SECTIONS CARRYING ADDITIONAL MASS

IN-PLANE VIBRATIONS OF CURVED BEAMS WITH VARIABLE CROSS-SECTIONS CARRYING ADDITIONAL MASS 11 th International Conferene on Vibration Problems Z. Dimitrovová et al. (eds.) Lisbon, Portugal, 9-1 September 013 IN-PLANE VIBRATIONS OF CURVED BEAMS WITH VARIABLE CROSS-SECTIONS CARRYING ADDITIONAL

More information

Failure Assessment Diagram Analysis of Creep Crack Initiation in 316H Stainless Steel

Failure Assessment Diagram Analysis of Creep Crack Initiation in 316H Stainless Steel Failure Assessment Diagram Analysis of Creep Crak Initiation in 316H Stainless Steel C. M. Davies *, N. P. O Dowd, D. W. Dean, K. M. Nikbin, R. A. Ainsworth Department of Mehanial Engineering, Imperial

More information

Lecture 3 - Lorentz Transformations

Lecture 3 - Lorentz Transformations Leture - Lorentz Transformations A Puzzle... Example A ruler is positioned perpendiular to a wall. A stik of length L flies by at speed v. It travels in front of the ruler, so that it obsures part of the

More information

arxiv:gr-qc/ v7 14 Dec 2003

arxiv:gr-qc/ v7 14 Dec 2003 Propagation of light in non-inertial referene frames Vesselin Petkov Siene College, Conordia University 1455 De Maisonneuve Boulevard West Montreal, Quebe, Canada H3G 1M8 vpetkov@alor.onordia.a arxiv:gr-q/9909081v7

More information

12.1 Events at the same proper distance from some event

12.1 Events at the same proper distance from some event Chapter 1 Uniform Aeleration 1.1 Events at the same proper distane from some event Consider the set of events that are at a fixed proper distane from some event. Loating the origin of spae-time at this

More information

Einstein s Three Mistakes in Special Relativity Revealed. Copyright Joseph A. Rybczyk

Einstein s Three Mistakes in Special Relativity Revealed. Copyright Joseph A. Rybczyk Einstein s Three Mistakes in Speial Relativity Revealed Copyright Joseph A. Rybzyk Abstrat When the evidene supported priniples of eletromagneti propagation are properly applied, the derived theory is

More information

Fig Review of Granta-gravel

Fig Review of Granta-gravel 0 Conlusion 0. Sope We have introdued the new ritial state onept among older onepts of lassial soil mehanis, but it would be wrong to leave any impression at the end of this book that the new onept merely

More information

5.2 THE EFFECT OF TURBULENCE ON CLOUD MICROSTRUCTURE, PRECIPITATION FORMATION AND THE ORGANISATION OF STRATOCUMULUS AND SHALLOW CUMULUS CONVECTION

5.2 THE EFFECT OF TURBULENCE ON CLOUD MICROSTRUCTURE, PRECIPITATION FORMATION AND THE ORGANISATION OF STRATOCUMULUS AND SHALLOW CUMULUS CONVECTION 5.2 THE EFFECT OF TURBULENCE ON CLOUD MICROSTRUCTURE, PRECIPITATION FORMATION AND THE ORGANISATION OF STRATOCUMULUS AND SHALLOW CUMULUS CONVECTION Charmaine N. Franklin * Centre for Australian Weather

More information

A population of 50 flies is expected to double every week, leading to a function of the x

A population of 50 flies is expected to double every week, leading to a function of the x 4 Setion 4.3 Logarithmi Funtions population of 50 flies is epeted to doule every week, leading to a funtion of the form f ( ) 50(), where represents the numer of weeks that have passed. When will this

More information

IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED ON THE KELVIN- VOIGHT PAIR

IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED ON THE KELVIN- VOIGHT PAIR Bulletin of the Transilvania University of Braşov Series II: Forestry Wood Industry Agriultural Food Engineering Vol. 9 (58) No. - 06 IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED

More information

A NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL

A NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL Journal of Geongineering, Vol. Yi-Chuan 4, No. 1, Chou pp. 1-7, and April Yun-Mei 009 Hsiung: A Normalized quation of Axially Loaded Piles in lasto-plasti Soil 1 A NORMALIZD QUATION OF AXIALLY LOADD PILS

More information

INFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE BEHAVIOR OF HYDRAULIC CYLINDER SUBJECTED TO JERKY MOTION

INFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE BEHAVIOR OF HYDRAULIC CYLINDER SUBJECTED TO JERKY MOTION Proeedings of ICFDP 8: 8 th International Congress of Fluid Dynamis & Propulsion Deember 14-17, 006, Sharm El-Shiekh, Sinai, Egypt ICFDP8-EG-154 INFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE

More information

UPPER-TRUNCATED POWER LAW DISTRIBUTIONS

UPPER-TRUNCATED POWER LAW DISTRIBUTIONS Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.

More information

Singular Event Detection

Singular Event Detection Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate

More information

KINETICS OF IRON OXIDE DIRECT REDUCTION BY COAL E.R. ABRIL 1

KINETICS OF IRON OXIDE DIRECT REDUCTION BY COAL E.R. ABRIL 1 KINETICS OF IRON OXIDE DIRECT REDUCTION BY COAL E.R. ABRIL 1 CIMM- Av.Velez Sarsfield 1561 C.P.5000 Córdoba, Argentina. aabril@intiemor.gov.ar Abstrat - A new interpretation to the kinetis of iron oxide

More information

DEVELOPMENT OF A MULTI-FEED P-T WELLBORE MODEL FOR GEOTHERMAL WELLS

DEVELOPMENT OF A MULTI-FEED P-T WELLBORE MODEL FOR GEOTHERMAL WELLS PROCEEDINGS, Thirty-First Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 3-February 1, 6 SGP-TR-179 DEVELOPMENT OF MULTI-FEED P-T WELLBORE MODEL FOR GEOTHERML

More information

A Time-Dependent Model For Predicting The Response Of A Horizontally Loaded Pile Embedded In A Layered Transversely Isotropic Saturated Soil

A Time-Dependent Model For Predicting The Response Of A Horizontally Loaded Pile Embedded In A Layered Transversely Isotropic Saturated Soil IOSR Journal of Mehanial and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 232-334X, Volume 16, Issue 2 Ser. I (Mar. - Apr. 219), PP 48-53 www.iosrjournals.org A Time-Dependent Model For Prediting

More information

UNIT 1 OPEN CHANNEL FLOW 2 MARK QUESTIONS AND ANSWERS

UNIT 1 OPEN CHANNEL FLOW 2 MARK QUESTIONS AND ANSWERS DEPARTMENT: CIVIL ENGINEERING SEMESTER: IV- SEMESTER SUBJECT CODE / Name: CE53 / Applied Hydrauli Engineering 1. Define open hannel flow with examples. Examples: UNIT 1 OPEN CHANNEL FLOW MARK QUESTIONS

More information

STUDY OF INHERENT FREQUENCY OF HELMHOLTZ RESONATOR

STUDY OF INHERENT FREQUENCY OF HELMHOLTZ RESONATOR 005 WJTA Amerian Waterjet Conferene August -3, 005! Houston, Texas Paper 6B-4 STUDY OF INHERENT FREQUENCY OF HELMHOLT RESONATOR Gong Weili An Liqian Cui Longlian Xie Guixin Shool of Mehanis, Arhiteture

More information

Chapter 2 Linear Elastic Fracture Mechanics

Chapter 2 Linear Elastic Fracture Mechanics Chapter 2 Linear Elasti Frature Mehanis 2.1 Introdution Beginning with the fabriation of stone-age axes, instint and experiene about the strength of various materials (as well as appearane, ost, availability

More information

Wavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013

Wavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013 Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it

More information

THERMAL MODELING OF PACKAGES FOR NORMAL CONDITIONS OF TRANSPORT WITH INSOLATION t

THERMAL MODELING OF PACKAGES FOR NORMAL CONDITIONS OF TRANSPORT WITH INSOLATION t THERMAL MODELING OF PACKAGES FOR NORMAL CONDITIONS OF TRANSPORT WITH INSOLATION THERMAL MODELING OF PACKAGES FOR NORMAL CONDITIONS OF TRANSPORT WITH INSOLATION t Tehnial Programs and Servies/Engineering

More information

General Equilibrium. What happens to cause a reaction to come to equilibrium?

General Equilibrium. What happens to cause a reaction to come to equilibrium? General Equilibrium Chemial Equilibrium Most hemial reations that are enountered are reversible. In other words, they go fairly easily in either the forward or reverse diretions. The thing to remember

More information

Math 151 Introduction to Eigenvectors

Math 151 Introduction to Eigenvectors Math 151 Introdution to Eigenvetors The motivating example we used to desrie matrixes was landsape hange and vegetation suession. We hose the simple example of Bare Soil (B), eing replaed y Grasses (G)

More information

Wave Propagation through Random Media

Wave Propagation through Random Media Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene

More information

What are the locations of excess energy in open channels?

What are the locations of excess energy in open channels? Leture 26 Energy Dissipation Strutures I. Introdution Exess energy should usually be dissipated in suh a way as to avoid erosion in unlined open hannels In this ontext, exess energy means exess water veloity

More information

SEDIMENT TRANSPORT CALCULATION CONSIDERING COHESIVE EFFECTS AND ITS APPLICATION TO WAVE-INDUCED TOPOGRAPHIC CHANGE

SEDIMENT TRANSPORT CALCULATION CONSIDERING COHESIVE EFFECTS AND ITS APPLICATION TO WAVE-INDUCED TOPOGRAPHIC CHANGE Proeedings of the 7 th International Conferene on Asian and Paifi Coasts (APAC 03) Bali, Indonesia, September 4-6, 03 SEDIMENT TRANSPORT CALCULATION CONSIDERING COHESIVE EFFECTS AND ITS APPLICATION TO

More information

ON DYNAMICALLY EQUIVALENT FORCE SYSTEMS AND THEIR APPLICATION TO THE BALANCING OF A BROOM OR THE STABILITY OF A SHOE BOX

ON DYNAMICALLY EQUIVALENT FORCE SYSTEMS AND THEIR APPLICATION TO THE BALANCING OF A BROOM OR THE STABILITY OF A SHOE BOX Proeedings of DEC 04 ASME 004 Design Engineering ehnial Conferenes and Computers and Information in Engineering Conferene September 8-Otober, 004, Salt Lake City, Utah, USA DE C0 04-5 7 188 ON DYNAMICALLY

More information

The Hanging Chain. John McCuan. January 19, 2006

The Hanging Chain. John McCuan. January 19, 2006 The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a

More information

Method of removing Fe particulate matters from subway environments

Method of removing Fe particulate matters from subway environments International Journal of Environmental onitoring and Analysis 5; 3(): -6 Published online Deember 3, 4 (http://www.sienepublishinggroup.om/j/ijema) doi:.648/j.ijema.53. ISSN: 38-7659 (Print); ISSN: 38-7667

More information

13.Prandtl-Meyer Expansion Flow

13.Prandtl-Meyer Expansion Flow 3.Prandtl-eyer Expansion Flow This hapter will treat flow over a expansive orner, i.e., one that turns the flow outward. But before we onsider expansion flow, we will return to onsider the details of the

More information

Testing of the External Turning Tool Holder Equipped With Alternate Bars in Its Construction

Testing of the External Turning Tool Holder Equipped With Alternate Bars in Its Construction ISSN 1392 1320 MATERIALS SCIENCE (MEDŽIAGOTYRA). Vol. 21, No. 3. 2015 Testing of the External Turning Tool Holder Equipped With Alternate Bars in Its Constrution Marten MADISSOO 1, Argo ROSSNER 1, Jüri

More information

Time and Energy, Inertia and Gravity

Time and Energy, Inertia and Gravity Time and Energy, Inertia and Gravity The Relationship between Time, Aeleration, and Veloity and its Affet on Energy, and the Relationship between Inertia and Gravity Copyright 00 Joseph A. Rybzyk Abstrat

More information

Evaluation of a Dual-Load Nondestructive Testing System To Better Discriminate Near-Surface Layer Moduli

Evaluation of a Dual-Load Nondestructive Testing System To Better Discriminate Near-Surface Layer Moduli 52 TRANSPORTATION RESEARCH RECORD 1355 Evaluation of a Dual-Load Nondestrutive Testing System To Better Disriminate Near-Surfae Layer Moduli REYNALDO ROQUE, PEDRO ROMERO, AND BYRON E. RUTH Theoretial analyses

More information

Chapter 26 Lecture Notes

Chapter 26 Lecture Notes Chapter 26 Leture Notes Physis 2424 - Strauss Formulas: t = t0 1 v L = L0 1 v m = m0 1 v E = m 0 2 + KE = m 2 KE = m 2 -m 0 2 mv 0 p= mv = 1 v E 2 = p 2 2 + m 2 0 4 v + u u = 2 1 + vu There were two revolutions

More information

A simple expression for radial distribution functions of pure fluids and mixtures

A simple expression for radial distribution functions of pure fluids and mixtures A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.

More information

Modeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers

Modeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers A. E. Romanov et al.: Threading Disloation Density Redution in Layers (II) 33 phys. stat. sol. (b) 99, 33 (997) Subjet lassifiation: 6.72.C; 68.55.Ln; S5.; S5.2; S7.; S7.2 Modeling of Threading Disloation

More information

arxiv:gr-qc/ v2 6 Feb 2004

arxiv:gr-qc/ v2 6 Feb 2004 Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this

More information