UTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.
|
|
- Edwin Jennings
- 6 years ago
- Views:
Transcription
1 UTC Engineering 329 Proportional Controller Design for Speed System By John Beverly Green Team John Beverly Keith Skiles John Barker 24 Mar 2006
2 Introdution This experiment is intended test the variable speed mor that powers the aerar at the waste-water treatment plant in Chattanooga TN and design a proportional ontroller. The aerar mixes the water enable steady growth of the bateria that onsume the unwanted waste in the water. The management of the faility require the output of the mixer be variable. Previous tests developed the Steady State Operating Curve (SSOC) from whih the system gain was determined. The Step Response test was also previously performed and baseline values for the system gain, dead time and time onstant were determined. This test will use the data from the previous tests as a basis determine the ontroller gain and the effets it has on the system. This data an then be used tune the system ahieve the desired response for any given input. This report inludes Bakground and Theory, the Proedure used ollet the empirial data, the Results of the tests and the mathematial modeling, a Disussion of the test findings and a Conlusion with Reommendations. An Appendix ontaining illustrations and referenes is loated in the bak. The Bakground and Theory setion is broken down in explanations of the previous SSOC and Step Response tests along with the results of the experiments and the tests performed in this experiment. (Frequeny Response, Root Lous and Proportional Control and Offset) Likewise, the Proedure setion lists eah test that was performed, inluding the previous tests, and is broken in eah setion.
3 Bakground and Theory The following desription refers the diagram below. (Figure 1) Power input the pump drive mor, (M-247) is ontrolled by the ontroller (SRC 247). Power input is adjusted manually provide a perentage of tal power the mor drive. The power input is provided by a 3-phase drive. (SCZ 247) As the input is adjusted, the mor (M- 247) hanges speed aordingly. The mor speed is monired by an optial sensor (ST 247) whih reports the speed in revolutions per minute. (rpm) As the mor speed hanges the mixer/aerar speed also hanges eventually affeting the oxygenation of the waste water. (Note: The gear-box ontains a gear set that redues the speed of the mixer 1/100 of the mor speed. Normal operating range is between 2 and 17 rpm of the mixer and rpm for the mor.) Figure 1 Shemati of the Aerar/Mixer
4 Figure 2 is a blok diagram of the Aerar/Mixer System. The manipulated (manual) input is represented by m(t). The Aerar Mixer in the red box is the transfer funtion. The output in rpm is represented by (t). Figure 2 Blok Diagram of the Aerar/Mixer SSOC The graph below illustrates the Steady State Operating Curve (SSOC). It is the response of the mor in revolutions per minute (rpm) a given input power (%). Several tests were performed through a range of input power develop the data neessary form the urve. It allows the output be predited for any given input. The error bars indiate a onfidene range of 95%. The slope (system gain (K)) is onstant throughout the range and is 17.3 rpm/%. (See Table 1 in the Appendix setion of this report for a omplete list of the data.)
5 Steady State Operating Curve 1800 High Speed Output (RPM) Low Mid Mor Input (%) Figure 3 Steady State Operating Curve Step Response The effet of a stepped input was measured and modeled mathematially. (Step Response Test) The First Order Plus Dead Time (FOPDT) method was used develop the model. Figure 4 is an example of the step response and the derived mathemati model. The model was determined aurately represent the output at any input in the operating range ( rpm). The equation used model the response urve is shown below (Equation 1) where A is the amplitude of the step input, t is time, t d is the dead time (lag time between input and output), t o is the lag of the model step, τ is the time onstant (length of time reah steady state), K is the system gain at steady state. The
6 variables A, K and t o were adjusted manually until the urves mathed the shape of the output response urve. (The blue model urve lays on p of the magenta urve and mathes the shape Figure 4) t t o t d Equation 1: ( ) (1 ) τ Y t = AK e The results of the Step Response are as follows: (also Table 2 in Appendix) K = 17.1 rpm/% t o =.14 se τ =.19 se Step Up 50-70% Speed (RPM) Output Model Output Input Model Input Input (%) Tim e(s) Figure 4 Step Response
7 Frequeny Response The aerar mor was subjeted a sine wave input at various frequenies and input power levels throughout the operating range. The response of the mor was reorded and plotted. Figure 5 is a sample of one of the plots. Two important piees of information were gleaned from the sine wave plots; phase angle and amplitude ratio. The phase angle is the degree of lag between the output and input. It was determined using Equation 2, where t is the distane between the peaks of the input and output and T is the distane between the peaks of the input. (See Figure 5) The amplitude ratio is the height of the output sine wave ( Δ ) divided by the height of the input sine wave ( Δ m). (Equation 3) t T 0 Equation 2: PhaseAngle ( PA) = (360 ) Equation 3: AmplitudeRatio( AR) = Δ Δm
8 Frequeny Response Speed (RPM) t Δ T Δm Input (%) Time (se) Figure 5 Sine Response Graph The amplitude ratio and phase angle are then graphed against the frequeny in a Bode Plot. (See Figure 6) The plot for the amplitude ratio is log-log and the phase angle is semi-log. Four piees of information an be taken from the Bode Plots; f u, Ku, Order and K. (See Figure 6) The value for fu an be determined by loating the frequeny at the intersetion of a line drawn horizontally through the phase angle plot at 180 degrees and the data line. The Ku is the inverse of the intersetion point of the f u line (transferred the amplitude ratio plot) and the data line. (Inverse of the point value 1/Ku) The differential order of the plot an be found by the maximum slope of the urved portion. (This is a guideline develop the equation neessary model the urve.) The gain (k) an be found with a horizontal line drawn aross the p of the urve where it no longer inreases. (See Table 3-Frequeny Test Data in Results)
9 Bode Plot degree K Amplitude Ration (rpm/%) 10 1/Ku Frequeny (Hz) 1 Bode Plot fu Phase Angle (deg) -180 d Frequeny (Hz) -250 Figure 6 Bode Plots Amplitude Ratio and Phase Angle
10 Root Lous *NOTE: A detailed derivation of the following equations, along with more explanation is listed in the Appendix at the end of this report. The following diagram (Figure 7) depits the Open Loop System. The value of G (s) is K (ontroller gain) and the G(s) is ke. (System Transfer Funtion) The Open Loop τs + 1 Transfer Funtion (OLTF) is a produt of the transfer funtions of the system. (See Equation 4) Figure 7 Open Loop Blok Diagram Equation 4: OLTF = G( s) G( s) = K Ke τs + 1 The following diagram depits the Closed Loop System. The negative feedbak loop enables the system respond the output. The Closed Loop Transfer Funtion (CLTF) is listed below in Equation 5 in terms of the OLTF and the system values. The CLTF will permit further evaluation of the system mathematially.
11 Figure 8 Closed Loop Blok Diagram Equation 5: CLTF K Ke OLTF = τs + 1 K Ke = 1 + OLTF K Ke 1+ ( τs + 1) + K Ke ( τs + 1) τs + 1 = The following equation (Equation 6) is the Charateristi Equation (CE) used develop the Root Lous urve desribed in the next paragraph. The roots of the equation were determined for several values of K and the value of K=17.3 rpm/%. Table 4 in the Results setion lists the data olleted for eah K seleted. ***NOTE: This Equation inorporates Pade s Approximation in plae of e -. See Equation 6a in the appendix. 2 Equation 6: CE 2 τ s + ( τ + K K ) s + 1+ K K = The roots of Equation 4 at eah value of K are plotted on a graph in whih the Y-axis represents the imaginary portion of the root and the x-axis represents the real portion of the root. The figure below (Figure 9) shows the method of finding the value of K for a speified Deay Ratio. (Deay Ratio is the initial overshoot divided by the seond
12 overshoot. Overshoot is the height of the output overage above the set point. See Figure 10) A line drawn through the origin at 78 o from the real axis will ross the root lous urve at the K value for ¼ deay ratio. 1/500 and 1/10 deay ratios are likewise found at 45 o and 70 o respetively. Critial damping, where there is no overshoot, an be found at the point where the root lous urve would ross the real axis. Ku (Ultimate K ) (not shown in the illustration) is the point where the root lous urve rosses the imaginary (vertial) axis. 1/500 (45 deg) K =.023 ROOT LOCUS PLOT 1/10 (70 deg) K =.075 1/4 (78 deg) K = IMAGINARY AXIS Critial Damping REAL AXIS -16 Figure 9 Root Lous Plot
13 Figure 10 Deay Ratio **Note: Smaller deay ratios result in longer times for the system reah steady state while larger ratios result in higher overshoot. Proportional Controller and Offset Offset is the steady state error that ours due a hange in set point. Equation 7 is used alulate the offset for the hosen value of K. (System gain K is fixed at this point a value determined with the development of the SSOC and Step Response tests) Higher values of K result in less offset, but may ause higher overshoot. Also, the K should not be higher than the K u or the system will be unstable and never ahieve a steady state output. Equation 7: Δr offset = 1 + K K Equation 8 desribes the proportional ontroller. The variable m(t) is the output and is onsidered the manipulated variable with respet the ontroller.(src 247 Figure 1) K
14 is the ontroller gain, m is the bias (input power %) The variable e(t) equals r(t) (t), where r(t) is the hange in rpm from the starting point r ( Δ r ) and (t) refers the transmitter output (ST 247 Figure 1) that detets the mor speed. Equation 8: m( t) = m + k e( t) = m + K ( r( t) ( t)) If the set point is hanged ( Δ r) the proportional ontroller equation will hange aordingly. This means the (t) value from the transmitter will affet the steady state output (m(t)) proportionally aording the new values of the proportional ontroller equation and there will be an offset value per Equation 7. (See Appendix for a sample of the effets of hanging set point) Figure 11 - Offset
15 **NOTE: Seletion of the K will have signifiant effet on the offset and the deay ratio and must not exeed the value of Ku. Table 5 in the Results setion ontains a list offset values for eah K at the major divisions of deay ratio and eah hange in set point referred as delta r. ( Δr )
16 Proedure **NOTE: The SSOC and Step Response tests were previously performed but are reviewed here. 1. SSOC: The first set of tests performed on the system was the SSOC development. A set of data was olleted from several tests through the entire operating range of the system at inrements of 10%. The results were plotted. (see Figure 3) Table 1 in the Results setion ontains the average values and error information for eah input. 2. Step Response: The next step of the system evaluation was the Step Response test. FOPDT modeling was performed using an Exel spreadsheet. (See Figure 4) The modeling resulted in the data listed in Table 2 of the Results setion. 3. Frequeny Response: The frequeny test was performed at three baseline input levels. (35%, 65% and 85%) The range of frequenies between.01 Hz 3Hz was tested. The frequeny was varied in inrements of approximately.03 Hz.3 Hz with a minimum of ten tests per baseline. The data for eah test was graphed. (See Figure 5) The Amplitude Ratio and Phase Angle were determined for eah test and plotted in a Bode Plot. (See Figure 6) The data found during these tests is listed in Table 3 of the Results setiont. 4. Root Lous Modeling: The root lous graph was generated using an Exel spreadsheet and the values of K were determined for ¼, 1/10, and 1/500 deay ratios, Ku and ritial damping. The values of K ranged from Table 4 in the Results lists the data.
17 5. Proportional Controller and Offset: The values of offset were alulated at eah of the K values that orrespond Ku, Critial Deay and ¼, 1/10 and 1/500 deay ratio. Table 5 of the Results setion lists the data.
18 Results Table 3 list the results of the Frequeny Response Experiment. The values were ahieved using a First Order Plus Dead Time (FOPDT) mathematial model. The urve mathing shown in Figure 5 is representative of all the urves generated for eah of the three baselines tested. (35%, 65% and 85% input power) Table 3 Frequeny Response Data K 17.3 rpm/% t o.2 se.17 se τ The Root Lous testing resulted in the values of K for the important deay ratios. (Critial, Ku, ¼, 1/10 and 1/500) The urve generated an be seen in Figure 9. The table below lists all values of K that were used graph the urve by hand (for verifiation), but many more points were used generate the Exel graph in Figure 9. Table 4 Root Lous Data K Real Root Imaginary Root (+) Deay Ratio i i i ¼ deay i i i i i 1/500 deay i 1/10 deay & -8 0 Critial Deay Ku The table below lists the values of offset for a hange in set point at eah value of K for the orresponding deay ratio. These an be used determine whih value of K needs be used ahieve the needed output with onsideration of the offset. The offset for Ku
19 provides the least offset but will reate signifiant overshoot and a larger deay time before settling in a steady state. Table 5 Offset Values in rpm for eah deay ratio and Δr Δr Offset Values for hange in set point (r) (rpm) Δr deay K ¼ / / Ku Critial Δr Δr Δr Δr
20 Disussion This experiment reveals the offset values for eah K. The use of the Ku puts the system at the edge of stability but provided very little offset. This means the reliability of the output is good, but the system may osillate and experiene signifiant overshoot before it ahieves a steady state operation. Less stability problems exist at lower deay ratios but the system experienes offset values that may be o high. That is, the atual output may vary o far from the required output for signifiant hanges in the set point. Critial damping is very stable but produes very high values of offset and may take a o long reah steady state output.
21 Conlusions and Reommendations The system models throughout these experiments provide an aurate representation of the atual output. The proportional ontroller model allows reasonable preditions of the system response onsidering the offset. The table of offsets an be used adjust the proportional ontroller equations determine the system response. If the output needs be very stable and the deay time is unimportant then the ¼ deay value for K an be used. *NOTE: The K value of Ku is unstable and will not settle in a steady state response. (ontinuos osillation) Some value less than Ku ould be used but speifi values that are lose need be disreetly tested before reliable results an be determined that will not osillate ontinuously. Overshoot at values approahing Ku will be large and settling times will inrease. In general ¼ deay and 1/10 deay should work well for this system sine a small differene in the speed of the aerar will our. For example, at 1/10 deay the speed of the mor offset for hange in set point of 100 is 42 rpm while the aerar is a mere.42 rpm. Likewise for ¼ deay, the aerar hange in speed is very small. (1/100) These are stable values of ontroller gain and should provide aeptable results for all hanges in set point.
22 Appendies Table 1 lists the data used develop the SSOC. The Mor Input olumn lists the input test values. The Output Speed lists the atual mor speed in revolutions per minute (rpm). The Standard Deviation olumn lists the deviation of the olleted data for eah of the tests. Table 1. Experimental values used reate the SSOC. Mor Input (%) Output Speed (RPM) Standard Deviation The step response experiments led the development of the baseline values listed below: Table 2 Step Response Data K 17.3 rpm/% t o.2 se.17 se τ Root Lous Equation Derivations The following steps will develop the mathematial model method used in this experiment. The Laplae Domain equivalents will be used. The final values will be expressed in terms of time.
23 The illustration below depits an open loop system. The Open Loop Transfer Funtion (OLTF) in Equation 1, is the produt of the system transfer funtions. Equation 2 is further development of the equation using the values of the transfer funtions. Equation 1a: OLTF = G ( s) G( s) Equation 2a: OLTF = K Ke τs +1 The following illustration depits the Closed Loop system. The blok diagram above an represented algebraially as follows: (Equation 3) Equation 3a: C( s) R( s) G ( s) G( s) = 1+ G ( s) G( s) = CLTF This equation is the Closed Loop Transfer Funtion (CLTF). The equation an be represented in terms of the OLTF as follows in Equation 4: Equation 4a: OLTF CLTF = 1+ OLTF
24 The values for the OLTF an be substituted develop Equation 5. Equation 5a: CLTF K Ke τs + 1 K Ke 1+ τs + 1 K Ke = ( τs + 1) + K Ke = ( τs + 1) The denominar of the CLTF (Equation 5) will be used with Pade s Approximation (Equation 6) develop an equation that will be used for several possible values of K. The roots of the polynomial equation set equal zero (Equation 8) will be determined using the quadrati equation (Equation 9) and plotted for eah value of K. Pade s Approximation is a substitution for the exponential funtion (e - ) allow the CLTF denominar be expanded in a meaningful polynomial where roots an be found. Equation 6a: Pade s Approximation e 1 s = 2 1+ s 2 Equation 7a: CLTF Denominar with Pade s Approximation 1 s ( τs + 1) + K K 2 ( τs + 1) = s 2
25 ,whih expands the polynomial in Equation 8: Equation 8a: τ s + ( τ + K K ) s + 1+ K K = 0,whih is the Charateristi equation (CE) in terms of the Laplae Domain that is used develop the Root Lous Curve for several values of K at the predetermined value of K. (17.3 rpm/%) 2 Equation 9a: Quadrati Equation 2 b 4a b + = roots 2a, where the oeffiients a,b and are: a = τ, b = ( τ + K K ), = 1+ K K Refer Figure 10 Root Lous urve see a plot of the roots found for eah frequeny. Changes in Set Point Proportional Controller m( t) = m + k e( t) = m + K ( r( t) ( t)) r = initial set point and r = the new set point Example: r = 680, m = 40%, r = 400, Δr = 280, K =.07, K = 17 r ( t) = Δr = r r, = 0 initial setting m( t) = (680 ( t))
26 m( t) = ( t) where (t) is the speed reported by ST 247 of Figure 1 the ontroller SRC 247. So, when the speed hanges the ontroller hanges the drives output at SCZ 247 aording the equation. If the setpoint hanges r ( t) = Δr = r r, = where the new set point is Offset Δr 280 = = 128rpm 1 + KK 1+.07(17) The new proportional equation is m ( t) = ( 280 ( t)) m ( t) = ( t) so aording the new proportional equation values the response will our with an offset of 128 rpm and the speed reported by ST 247 will ause the ontroller adjust the response aordingly with a 128 rpm offset. (differene in atual response) See Figure 11 for a graphial illustration.
27 List of Figures 1- System Shemati 2- System Blok Diagram 3- SSOC 4- Step Response 5-5Frequeny Response 6- Bode Diagrams 7- Open Loop 8- Closed Loop 9- Root Lous Plot 10- Deay Ratio 11- Offset List of Tables 1 - SSOC Data 2 Step Response Data 3 Frequeny Response Data 4 Root Lous Data 5 Offset Data List of Equations 1 FOPDT 2- Phase Angle 3 - Amplitude Ratio 4 - OLTF 5 - CLTF 6 - Charateristi Equation 7 - Offset 8 - Proportional Controller Referenes Smith and Corripio Priniples and Praties of Aumati Proess Control, Edition 3E, Wiley Publishing 2006 Dr. Henry 328 Website, UTC Blak Board
Yellow Team Pressure Control System Proportional Control: Model and Experiment
Yellow Team Pressure Control System Proportional Control: Model and Experiment Team Members: Jason Hixson Laura Amini Mike Bradley UTC ENGR 8 Nov. th, Outline Bakground System Diagram Operating Range SSOC
More informationController Design Based on Transient Response Criteria. Chapter 12 1
Controller Design Based on Transient Response Criteria Chapter 12 1 Desirable Controller Features 0. Stable 1. Quik responding 2. Adequate disturbane rejetion 3. Insensitive to model, measurement errors
More information23.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 informationMillennium 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 informationEvaluation 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 informationProportional Controller Performance for Aerator Mixer System
1 Proportional Controller Performance for Aerator Mixer System By Nicholas University of Tennessee at Chattanooga ENGR 329-1 Green Team (Monty Veal, TJ Hurless) April 2th, 21 2 Introduction- The experiment
More informationCHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS. Professor Dae Ryook Yang
CHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS Professor Dae Ryook Yang Spring 208 Dept. of Chemial and Biologial Engineering 0- Road Map of the Leture X Stability of losed-loop ontrol system
More informationThe 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 informationMeasuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach
Measuring & Induing Neural Ativity Using Extraellular Fields I: Inverse systems approah Keith Dillon Department of Eletrial and Computer Engineering University of California San Diego 9500 Gilman Dr. La
More informationIMPEDANCE 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 informationWavetech, 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 informationAcoustic Waves in a Duct
Aousti Waves in a Dut 1 One-Dimensional Waves The one-dimensional wave approximation is valid when the wavelength λ is muh larger than the diameter of the dut D, λ D. The aousti pressure disturbane p is
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationDIGITAL 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 informationECE-320 Linear Control Systems. Winter 2013, Exam 1. No calculators or computers allowed, you may leave your answers as fractions.
ECE-320 Linear Control Systems Winter 2013, Exam 1 No alulators or omputers allowed, you may leave your answers as frations. All problems are worth 3 points unless noted otherwise. Total /100 1 Problems
More informationTheory. Coupled Rooms
Theory of Coupled Rooms For: nternal only Report No.: R/50/TCR Prepared by:. N. taey B.., MO Otober 00 .00 Objet.. The objet of this doument is present the theory alulations to estimate the reverberant
More informationGeneral 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 informationA 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 informationFiber 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 informationThe simulation analysis of the bridge rectifier continuous operation in AC circuit
Computer Appliations in Eletrial Engineering Vol. 4 6 DOI 8/j.8-448.6. The simulation analysis of the bridge retifier ontinuous operation in AC iruit Mirosław Wiślik, Paweł Strząbała Kiele University of
More informationStabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer
Proeedings of the 4th IIAE International Conferene on Industrial Appliation Engineering 216 Stabilization of the Preision Positioning Stage Working in the Vauum Environment by Using the Disturbane Observer
More informationDynamic Simulation and Composition Control in A 10 L Mixing Tank
roeeding of International Conferene on Chemial and Material Engineering ISBN : 978--97-8-7 SE.4 - Dynami Simulation and Composition Control in A L Mixing Tank Yulius Deddy Hermawan Chemial Engineering
More informationMODELLING 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 informationKINETICS 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 informationRESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED
The Seventh Asia-Paifi Conferene on Wind Engineering, November 8-1, 9, Taipei, Taiwan RESEARCH ON RANDOM FORIER WAVE-NMBER SPECTRM OF FLCTATING WIND SPEED Qi Yan 1, Jie Li 1 Ph D. andidate, Department
More informationSpray Boot Pressure Station UTC -ENGR 3280-L Week 10 March 20, 2013 Blue Team. Ethan Tummins Jeff Clowdus Jerry Basham
Spray Boot Pressure Station UTC -ENGR 3280-L Week 10 March 20, 2013 Blue Team Ethan Tummins Jeff Clowdus Jerry Basham Presentation Overview Spray Booth Pressure Station Overview Steady State Operating
More informationWhat 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 informationINTERNATIONAL 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 informationImprovements in the Modeling of the Self-ignition of Tetrafluoroethylene
Exerpt from the Proeedings of the OMSOL onferene 010 Paris Improvements in the Modeling of the Self-ignition of Tetrafluoroethylene M. Bekmann-Kluge 1 *,. errero 1, V. Shröder 1, A. Aikalin and J. Steinbah
More informationAn Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems
An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,
More informationCanimals. borrowed, with thanks, from Malaspina University College/Kwantlen University College
Canimals borrowed, with thanks, from Malaspina University College/Kwantlen University College http://ommons.wikimedia.org/wiki/file:ursus_maritimus_steve_amstrup.jpg Purpose Investigate the rate of heat
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationThe 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 informationAn Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances
An aptive Optimization Approah to Ative Canellation of Repeated Transient Vibration Disturbanes David L. Bowen RH Lyon Corp / Aenteh, 33 Moulton St., Cambridge, MA 138, U.S.A., owen@lyonorp.om J. Gregory
More informationSupplementary Materials
Supplementary Materials Neural population partitioning and a onurrent brain-mahine interfae for sequential motor funtion Maryam M. Shanehi, Rollin C. Hu, Marissa Powers, Gregory W. Wornell, Emery N. Brown
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationUPPER-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 information2. 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 informationSOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS
SOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS Prepared by S. Broverman e-mail 2brove@rogers.om website http://members.rogers.om/2brove 1. We identify the following events:. - wathed gymnastis, ) - wathed baseball,
More informationUNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS. V. N. Matveev and O. V. Matvejev
UNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS V. N. Matveev and O. V. Matvejev Joint-Stok Company Sinerta Savanoriu pr., 159, Vilnius, LT-315, Lithuania E-mail: matwad@mail.ru Abstrat
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationWave 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 informationSupplementary information for: All-optical signal processing using dynamic Brillouin gratings
Supplementary information for: All-optial signal proessing using dynami Brillouin gratings Maro Santagiustina, Sanghoon Chin 2, Niolay Primerov 2, Leonora Ursini, Lu Thévena 2 Department of Information
More informationf 2 f n where m is the total mass of the object. Expression (6a) is plotted in Figure 8 for several values of damping ( ).
F o F o / k A = = 6 k 1 + 1 + n r n n n RESONANCE It is seen in Figure 7 that displaement and stress levels tend to build up greatly when the oring requeny oinides with the natural requeny, the buildup
More information9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry
9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived
More informationControl Systems. Control Systems Design Lead-Lag Compensator.
Design Lead-Lag Compensator hibum@seoulteh.a.kr Outline Lead ompensator design in frequeny domain Lead ompensator design steps. Example on lead ompensator design. Frequeny Domain Design Frequeny response
More informationDeveloping 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 informationA Novel Process for the Study of Breakage Energy versus Particle Size
Geomaterials, 2013, 3, 102-110 http://dx.doi.org/10.4236/gm.2013.33013 Published Online July 2013 (http://www.sirp.org/journal/gm) A Novel Proess for the Study of Breakage Energy versus Partile Size Elias
More informationEE 321 Project Spring 2018
EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,
More informationResearch Report Frank McCullough, and Ned H. Burns
Te:hni:al Report Do:umentotion Page 1, Report No. 2. Government Aession No. FHWATX-91+556-3 4. Title and Subtitle ANALYSIS OF URLING MOVEMENTS AND ALIBRATION OF PP PROGRAM 3. Reipient's atalog No. 5. Report
More informationSimplified 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 informationPart G-4: Sample Exams
Part G-4: Sample Exams 1 Cairo University M.S.: Eletronis Cooling Faulty of Engineering Final Exam (Sample 1) Mehanial Power Engineering Dept. Time allowed 2 Hours Solve as muh as you an. 1. A heat sink
More informationThe 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 informationCalculation of Desorption Parameters for Mg/Si(111) System
e-journal of Surfae Siene and Nanotehnology 29 August 2009 e-j. Surf. Si. Nanoteh. Vol. 7 (2009) 816-820 Conferene - JSSS-8 - Calulation of Desorption Parameters for Mg/Si(111) System S. A. Dotsenko, N.
More informationMotor 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 informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More information10.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 informationSTUDY OF INTERFACIAL BEHAVIOR OF CNT/POLYMER COMPOSITE BY CFE METHOD
THE 19TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS STUDY OF INTERFACIAL BEHAVIOR OF CNT/POLYMER COMPOSITE BY CFE METHOD Q. S. Yang*, X. Liu, L. D. Su Department of Engineering Mehanis, Beijing University
More informationAyan Kumar Bandyopadhyay
Charaterization of radiating apertures using Multiple Multipole Method And Modeling and Optimization of a Spiral Antenna for Ground Penetrating Radar Appliations Ayan Kumar Bandyopadhyay FET-IESK, Otto-von-Guerike-University,
More informationRobust Flight Control Design for a Turn Coordination System with Parameter Uncertainties
Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo
More informationTHERMAL 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 informationCombined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2
Sensor and Simulation Notes Note 53 3 May 8 Combined Eletri and Magneti Dipoles for Mesoband Radiation, Part Carl E. Baum University of New Mexio Department of Eletrial and Computer Engineering Albuquerque
More informationAC : A GRAPHICAL USER INTERFACE (GUI) FOR A UNIFIED APPROACH FOR CONTINUOUS-TIME COMPENSATOR DESIGN
AC 28-1986: A GRAPHICAL USER INTERFACE (GUI) FOR A UNIFIED APPROACH FOR CONTINUOUS-TIME COMPENSATOR DESIGN Minh Cao, Wihita State University Minh Cao ompleted his Bahelor s of Siene degree at Wihita State
More informationA NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS
Vietnam Journal of Mehanis, VAST, Vol. 4, No. (), pp. A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Le Thanh Tung Hanoi University of Siene and Tehnology, Vietnam Abstrat. Conventional ship autopilots are
More informationFailure 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 informationFour-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 informationExperiment 3: Basic Electronic Circuits II (tbc 1/7/2007)
Experiment 3: Basi Eletroni iruits II (tb /7/007) Objetive: a) To study the first-order dynamis of a apaitive iruits with the appliation of Kirhoff s law, Ohm s law and apaitane formula. b) To learn how
More informationSpeed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion
Speed-feedbak Diret-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Y. Yamamoto, T. Nakamura 2, Y. Takada, T. Koseki, Y. Aoyama 3, and Y. Iwaji 3
More informationUniversities of Leeds, Sheffield and York
promoting aess to White Rose researh papers Universities of Leeds, Sheffield and York http://eprints.whiterose.a.uk/ This is an author produed version of a paper published in Journal of Composites for
More informationDuct Acoustics. Chap.4 Duct Acoustics. Plane wave
Chap.4 Dut Aoustis Dut Aoustis Plane wave A sound propagation in pipes with different ross-setional area f the wavelength of sound is large in omparison with the diameter of the pipe the sound propagates
More informationCHAPTER 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 informationSlenderness Effects for Concrete Columns in Sway Frame - Moment Magnification Method
Slenderness Effets for Conrete Columns in Sway Frame - Moment Magnifiation Method Slender Conrete Column Design in Sway Frame Buildings Evaluate slenderness effet for olumns in a sway frame multistory
More informationConveyor 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 informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationEXPERIMENTAL INVESTIGATION OF CRITICAL STRAIN FOR VISCOPLASTIC MATERIALS WITH DIFFERENT GEOMETRIES AND RHEOMETERS
EXPERIMENTAL INVESTIGATION OF CRITICAL STRAIN FOR VISCOPLASTIC MATERIALS WITH DIFFERENT GEOMETRIES AND RHEOMETERS Guilherme A. S. Balvedi, 2 Diogo E. V. Andrade, 3 Admilson T. Frano e 3 Cezar O. R. Negrão
More informationSingular 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 informationQCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines
DOI.56/sensoren6/P3. QLAS Sensor for Purity Monitoring in Medial Gas Supply Lines Henrik Zimmermann, Mathias Wiese, Alessandro Ragnoni neoplas ontrol GmbH, Walther-Rathenau-Str. 49a, 7489 Greifswald, Germany
More informationAnalysis 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 information7 Max-Flow Problems. Business Computing and Operations Research 608
7 Max-Flow Problems Business Computing and Operations Researh 68 7. Max-Flow Problems In what follows, we onsider a somewhat modified problem onstellation Instead of osts of transmission, vetor now indiates
More information12.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 informationChapter 2 Lecture 9 Longitudinal stick fixed static stability and control 6 Topics
hapter Leture 9 Longitudinal stik fied stati stability and ontrol 6 Topis Eample.4 Eample.4 Referene.4 desribes the stability and ontrol data for ten airplanes. This inludes a general aviation airplane
More informationON-LINE ESTIMATION OF THE VENTILATION RATE OF GREENHOUSES
ON-LINE ESTIMATION OF THE VENTILATION RATE OF GREENHOUSES J. Bontsema *, E.J. van Henten *, J.G. Kornet *, J. Budding ** and T. Rieswijk ** * Agrotehnology and Food Innovations B.V., Greenhouse Tehnology
More informationExperiment 03: Work and Energy
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physis Department Physis 8.01 Purpose of the Experiment: Experiment 03: Work and Energy In this experiment you allow a art to roll down an inlined ramp and run into
More informationINFLUENCE 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 information6.4 Dividing Polynomials: Long Division and Synthetic Division
6 CHAPTER 6 Rational Epressions 6. Whih of the following are equivalent to? y a., b. # y. y, y 6. Whih of the following are equivalent to 5? a a. 5, b. a 5, 5. # a a 6. In your own words, eplain one method
More informationAircraft CAS Design with Input Saturation Using Dynamic Model Inversion
International Journal of Control, Automation, and Systems Vol., No. 3, September 003 35 Airraft CAS Design with Input Saturation sing Dynami Model Inversion Sangsoo Lim and Byoung Soo Kim Abstrat: This
More informationDevelopment 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 informationCritical 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 information3 Tidal systems modelling: ASMITA model
3 Tidal systems modelling: ASMITA model 3.1 Introdution For many pratial appliations, simulation and predition of oastal behaviour (morphologial development of shorefae, beahes and dunes) at a ertain level
More informationANALYSIS OF THE SEISMIC BEHAVIOR OF KEDDARA DAM USING STRONG MOTION RECORDS
ANALYSIS OF THE SEISMIC BEHAVIOR OF KEDDARA DAM USING STRONG MOTION RECORDS S. Louadj, R. Bahar, E. Vinens, N. Laouami 4 Ph.D Student, Geomaterials and Environment Laboratory LGEA, University of Tizi-Ouzou,
More informationA 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 informationPhysical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena
Page 1 of 10 Physial Laws, Absolutes, Relative Absolutes and Relativisti Time Phenomena Antonio Ruggeri modexp@iafria.om Sine in the field of knowledge we deal with absolutes, there are absolute laws that
More informationEinstein 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 informationCALCULATION 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 informationThe coefficients a and b are expressed in terms of three other parameters. b = exp
T73S04 Session 34: elaxation & Elasti Follow-Up Last Update: 5/4/2015 elates to Knowledge & Skills items 1.22, 1.28, 1.29, 1.30, 1.31 Evaluation of relaxation: integration of forward reep and limitations
More informationSubject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (
16.50 Leture 0 Subjet: Introdution to Component Mathing and Off-Design Operation At this point it is well to reflet on whih of the many parameters we have introdued (like M, τ, τ t, ϑ t, f, et.) are free
More informationExamining 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 informationMath 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 informationMoment Curvature Characteristics for Structural Elements of RC Building
Moment Curvature Charateristis for Strutural Elements of RC Building Ravi Kumar C M 1,*, Vimal Choudhary 2, K S Babu Narayan 3 and D. Venkat Reddy 3 1 Researh Sholar, 2 PG Student, 3 Professors, Department
More informationADHESION MEASURES OF ELASTO-PLASTIC THIN FILM VIA BUCKLE-DRIVEN DELAMINATION
ADHESION MEASURES OF ELASTO-PLASTIC THIN FILM VIA BUCKLE-DRIVEN DELAMINATION Yu Shouwen and Li Qunyang Department of Engineering Mehanis, Tsinghua University, Beijing 184, China Yusw@mail.tsinghua.edu.n
More information