Seismic Response Reduction of a Building Using. Top-Story Isolation System with MR Damper

Similar documents
Hyun-Su Kim Division of Architecture, Sunmoon Univ., Asan, , Korea,

A Comparison between a New Semi-Active Tuned Mass Damper and an Active Tuned Mass Damper

Development of Experimental Benchmark Problems for International Collaboration in Structural Response Control

PERFORMANCE OF TUNED MASS DAMPERS FOR RESPONSE REDUCTION OF STRUCTURES UNDER NEAR-FIELD AND FAR-FIELD SEISMIC EXCITATIONS

Application of Isolation Technology in High-voltage Electrical Equipments

A Strategy for the Optimization of Damper Configurations based on Building Performance Indices

Seismic Control of Continuous Bridges Using Variable Radius Friction Pendulum Systems and Viscous Fluid Dampers

Hybrid Shake Table for the Testing of Midlevel Seismic Isolation Systems

Smart Outriggers for Seismic Protection of High-rise Buildings

Intelligent Base Isolation Systems

Influence of Passive Isolator Parameters on Seismic Response Control of Multi-Span Bridges

Experimental Study on the Outrigger Damping System for High-Rise Building

Optimal Performance-based

Dynamic Analysis of High Rise Seismically Isolated Buildings

Nonlinear Dynamic Analysis of Base Isolated Reinforced Concrete Building

BEHAVIOUR OF SOIL-STRUCTURE SYSTEM WITH TUNED-MASS DAMPERS DURING NEAR-SOURCE EARTHQUAKES

The Structural Performance of RC Coupling Beam. Using Cold-Formed Steel Channel Sections for. Alternative Detail

Kiran K. Shetty & Krishnamoorthy Department of Civil Engineering, Manipal Institute of Technology, Manipal , Karnataka, India

SEISMIC DESIGN OF STRUCTURE

Multi-storey structures with compound seismic isolation

Seismic Analysis of Structure by Introducing Weak Story on Top as Tuned Mass Damper

Study of TMD-controlled Reinforced Concrete Framed Structure under Code-based Artificial Earthquake

DYNAMIC TEST OF MULTIPLE TUNED MASS DAMPERS FOR VIBRATION CONTROL OF HIGH-RISE BUILDINGS

Effectiveness of Base Isolation Technique and Influence of Isolator Characteristics on Response of a Base Isolated Building

Investigation of PTMD System Affected by Parkfield Near-Field Earthquake

EQUIVALENT UNIFORM DAMPING RATIOS for Irregular in Height Concrete / Steel Structural systems

DESIGN OF OPTIMUM PARAMETERS OF TUNED MASS DAMPER FOR A G+8 STORY RESIDENTIAL BUILDING

SEISMIC ANALYSIS OF BUILDING USING TWO TYPES OF PASSIVE ENERGY DISSIPATION DEVICES

ADVANCED PASSIVE CONTROL TECHNIQUES FOR RETROFIT OF EXISTING BUILDINGS IN SEISMIC ZONE

R. J. McNamara, S.E.

IJSRD - International Journal for Scientific Research & Development Vol. 4, Issue 05, 2016 ISSN (online):

OPTIMUM DAMPING IN ISOLATION SYSTEM

Earthquake Response Analysis of Spherical Tanks with Seismic Isolation

EXPERIMENTAL STUDY OF AN ADAPTIVE BASE ISOLATION SYSTEM FOR BUILDINGS

SEISMIC RESPONSE OF MULTISTOREYED STEEL FRAME WITH VISCOUS FLUID SCISSOR JACK DAMPERS

Effective Viscous Damping Ratio in Seismic Response of Reinforced Concrete Structures

1554. Influence of friction variability on isolation performance of a rolling-damper isolation system

Frequency Control of Smart Base Isolation System Employing A Novel Adaptive Magneto-rheological Elastomer Base Isolator

Performance Improvement of Tall Reinforced Concrete Structures with Multiple Tuned Mass Dampers (MTMD)

A NOVEL ELASTOMERIC BASE ISOLATION SYSTEM FOR SEISMIC MITIGATION OF LOW-RISE BUILDINGS

Osamu Takahashi, Hiromasa Aida 1 Junji Suhara, Ryoichiro Matsumoto 2 Yasuo Tsuyuki 3 Takafumi Fujita 4. Kozo Keikaku Engineering 2

Effectiveness of Tuned Mass Dampers in Mitigating Earthquake Ground Motions in Low and Medium Rise Buildings

Development and Application of Seismic Isolation and Response Control of Buildings in Japan

Seismic Vibration Control of a MDOF Building Connected with Viscous and SAVFD Dampers

RESPONSE CONTROL BY DAMPER DEVICES OF HIGH-RISE BUILDING UNDER LONG-PERIOD GROUND MOTION

Investigating the Effect of Pounding on Seismic Response of Isolated Structures

IMPROVING SEISMIC PERFORMANCE: ADD STIFFNESS OR DAMPING? Trevor E. Kelly 1

EARTHQUAKE RESPONSE ANALYSIS OF A MID-STORY SEISMIC ISOLATED BUILDING

FULL-SCALE VERIFICATION TEST OF DYNAMIC RESPONSE CONTOROL TECHNIQUES FOR STRONG EARTHQUAKES

PSEUDO NEGATIVE STIFFNESS DAMPERS FOR EARTHQUAKE RESPONSE CONTROL OF CABLE-STAYED BRIDGES

Chapter 3 Passive Protective Schemes for Rail- Counterweight System

STUDY ON SHAKING TABLE TESTS OF ISOLATED BRIDGE MODEL WITH LRB

Seismic Response of Low-Rise Base Isolated Structures

Dynamics of structures

COMPARISON BETWEEN TWO OPTIMIZATION PROCEDURES FOR DAMPER LOCATION IN STEEL FRAME STRUCTURES

ROBUST DESIGN OF TMD SYSTEMS FOR THE SEISMIC RESPONSE CONTROL OF ASYMMETRIC SOFT STOREY BUILDING

Study Response of Fixed Base and Isolation Base

Development of a novel multi-layer MRE isolator for suppression of building vibrations under seismic events

Investigation of Seismic Response of Steel MRF Building Structures with Nonlinear Viscous Dampers Using Real- Time Hybrid Simulation

Seismic Behaviour Of A Soft Storey Building With & Without Viscous Dampers

Received 28 August 2017; received in revised form 16 October 2017; accepted 27 October 2017 DOI

SEISMIC RESPONSE OF UNSYMMETRIC BUILDING WITH OPTIMALLY PLACED FRICTION DAMPERS

Choice of in-structure Damping Model: Do we have an answer? Answer?

Chinmoy Kolay Research Engineer

Integrated Optimal Life Cycle Design of Super Tall Buildings with Viscous Dampers

EFFECT OF SOFT STOREY LEVEL VARIATION ON TORSIONAL ROTATION OF LEAD RUBBER BASE ISOLATED STRUCTURE

Assistant Professor, Applied Mechanics Department, Government College of Engineering, Amravati, India 2

K. Ikago 1, N. Inoue 2 1. INTRODUCTION 2. THE CONCEPT OF DISPLACEMENT CONTORL DESIGN

Smart Base Isolation Strategies Employing Magnetorheological Dampers

Effect of Building Height on Torsional Response of Lead Rubber Bearing Base Isolated Structures-A Study

Performance of Tuned Mass Dampers in Vibration Response Control of Base-Excited Structures

COMPARITIVE STUDY OF BASE ISOLATORS AND VISCOUS FLUID DAMPERS ON SEISMIC RESPONSE OF RC STRUCTURES

MULTI-LEVEL FORTIFICATION INTENSITIES SEISMIC PERFORMANCE ASSESSMENT FOR REINFORCED CONCRETE FRAME- SHEAR WALL STRUCTURE WITH VISCOUS DAMPERS

Floor Slab Isolation for Mitigating the Seismic Response of Building Systems

DYNAMIC ANALYSIS OF 25 STOREY RCC BUILDING WITH AND WITHOUT VISCOUS DAMPERS. Naziya Ghanchi1, Shilpa Kewate2

STFT Based Real-time Control with Variable Stiffness

A MULTI-FUNCTION UNSEATING PREVENTION DEVICE FOR ISOLATED BRIDGES

MID-STORY SEISMIC ISOLATION FOR STRENGTHENING OF A MULTI-STORY BUILDING


Using friction dampers for improving earthquake response of self-variable stiffness RC framed buildings

EARTHQUAKE RESPONSE ANALYSIS OF MID-STORY BUILDINGS ISOLATED WITH VARIOUS SEISMIC ISOLATION TECHNIQUES

sixteen seismic design Earthquake Design Earthquake Design Earthquake Design dynamic vs. static loading hazard types hazard types: ground shaking

SEISMIC RISK AND COST-EFFECTIVENESS OF BASE ISOLATED BUILDINGS IN PERU

SIGNIFICANCE OF ACCELERATION IN THE DESIGN OF DAMPERS FOR LIFELINE BUILDINGS IN SEISMIC GROUND MOTIONS

Vibration Control of Framed Structure Using Tuned Mass Damper

Field Dynamic Test for A Base-Isolated 15-Story Steel Structure

Wataru Goto. NTT FACILITIES, INC. in Japan. 15th World Conference on Earthquake Engineering 24 to 28 September 2012 LISBON PORTUGAL

Sensitivity Analysis of LRB Isolator for RC Shear Frame

Effectiveness of Tuned Liquid Column Dampers on Vibration Behavior of One-story Steel Structure under Harmonic Force and an Earthquake Record

Title. Author(s)M. OHSAKI; S. TSUDA; N. SUGIYAMA. Issue Date Doc URL. Type. Note. File Information DIRECTIONAL SEISMIC VIBRATION CONTROL

Application of Viscous Dampers in Seismic Rehabilitation of Steel Moment Resisting Frames

Seismic Response Control of the Building using Passive Devices

Large-Scale Real-Time Hybrid Simulations

Vibration control of a building model with base isolation

SEISMIC EXCITATION OF LOW TO HIGH RISE RCC STRUCTURE WITH LEAD RUBBER BEARING BASE ISOLATION

THE APPLICATION OF RESPONSE CONTROL DESIGN USING MIDDLE-STORY ISOLATION SYSTEM TO HIGH-RISE BUILDING

EFFECT OF FLUID VISCOUS DAMPERS IN MULTI-STOREYED BUILDINGS

Proceedings of the 3rd International Conference on Environmental and Geological Science and Engineering

Parametric Analysis of the Combination of Rubber and Sliding Bearings for Seismic Isolation of Steel Structures

Transcription:

Contemporary Engineering Sciences, Vol. 7, 04, no., 979-986 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ces.04.493 Seismic Response Reduction of a Building Using Top-Story Isolation System with MR Damper Hyun-Su Kim Division of Architecture, Sunmoon University, Asan-si, Korea Copyright 04 Hyun-Su Kim. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract In this study, the adaptability of a semi-active top-story isolation system with an MR damper for reduction of seismic responses of tall buildings has been investigated. To this end, a 0-story building structure was selected as an example. Low damping elastomeric bearings and an MR damper were used to compose a smart top-story isolation system. Artificial earthquakes generated based on design spectrum were used for structural analyses. A fuzzy logic controller (FLC) was used to control an MR damper and the FLC was optimized by multi-objective genetic algorithms (MOGA). Numerical simulations show that a smart top-story isolation system can effectively reduce both structural responses and isolation story drifts of the building structure compared to a passive system. Keywords: Vibration control, Semi-active control system, Top-story isolation, Seismic response reduction, Genetic algorithms, Multi-objective optimization Introduction The use of a base isolation system is a well-known approach in the seismic design of low-rise buildings. However, for technical reasons, this base isolation system cannot be effectively applied to high-rise buildings [, ]. The natural periods of tall buildings are relatively long, and changing the first mode of vibration requires a large flexibility in the isolation mechanism. This is considered a major technical problem, especially if a large misalignment of vertical force on the isolator is expected. The vertical stiffness of the isolator is also a technical concern. Based on these problems, a top- or mid-story seismic isolation system has been proposed [3, 4, 5]. The isolated upper stories of the tall building structure

980 Hyun-Su Kim with a top- or mid-story seismic isolation system are rolled to act as a large tuned mass damper (TMD). Recently, smart base isolation strategies have been widely studied as novel mitigation techniques to reduce structural damage caused by severe loads. Current literature shows that a smart base isolation system can reduce base drifts without accompanying increases in acceleration that are seen with passive strategies [6, 7]. Several researchers have investigated the use of semi-active smart dampers for seismic response mitigation as a component of a hybrid control system [7, 8]. Control methodologies for these semi-active smart base isolation systems include fuzzy control [8], sliding mode control [9], clipped-optimal control [7], and others. It has been shown that smart base isolation can protect low-rise building structures in extreme earthquakes without sacrificing performance during more frequent, moderate seismic events. However, a smart top-story isolation system for tall buildings has not been fully investigated to date. Therefore, the control performance of a smart top-story isolation system for a building structure subjected to seismic excitation was investigated in this paper. A smart top-story isolation system consists of MR dampers and low damping elastomeric bearings and a fuzzy logic controller (FLC) was used to control an MR damper [0]. A multi-objective genetic algorithm (MOGA) was selected to optimize the fuzzy rules and membership functions of FLC. Analytical models and artificial earthquake load In this study, a typical 0-story building, as shown in Figure (a), and an example building with a smart top-story isolation system, as shown in Figure (b), are employed as example structures for evaluating the control performance of a smart top-story isolation system. Each building is represented by a shear building model with 0 DOFs. As shown in Figure (b), a smart top-story isolation system consists of an MR damper and elastic bearings. In order to compare the control performance of a smart top-story isolation system, a passive top-story isolation system is used. In the passive system, an MR damper, as shown in Figure (b) is replaced by a proper viscous damper. m 0 MR Damper k 0 c 0 m k b 9 k 9 c 9 k 9 c 9 m 0 m 9 Elastic bearing m 3 m 3 k 3 c 3 m k 3 c 3 m k c k c m m k c k c (a) General building (b) Smart top-story isolation Figure. Analytical model of 0-story example structure

Seismic response reduction of a building 98 The equations of motion of the general building and the top-story isolated building are given by equations () and (). M x () Cx Κx ΜΕx g M x () F Cx Κ x ΜΕx g E FC where E is a unit vector with a size of [0x] for applying ground acceleration to each story mass, F is the MR damper force in the smart top-story isolation C F system and E is a location vector of the isolation system. This location vector contains element for the top story and 0 elements for the other stories. In the equations, x, x, x and x, x, x are the displacement, velocity and acceleration of the general building and the top-story isolated building, respectively, and x represents the ground acceleration. M is the 0 by 0 g dimensional mass matrix for both the general building and top-story isolated building. It is represented by a typical diagonal matrix with the diagonal element, 5 m i, that is the mass of the ith floor (i=,,.,0) and is 8 0 kg. K and K are the stiffness matrices of the general building and the top-story isolated building, respectively. In this example, the stiffness coefficient 9 of each floor ( k i ) is.35 0 N/m. In general, a top-story isolation system provides the best control performance when its dynamic behavior is similar to a tuned mass damper (TMD). Generally, the frequency of a TMD is tuned to the first modal frequency of the structure because the motion of a long period structure is generally governed by the first modal response. Accordingly, the stiffness coefficient ( k b ) of the elastic bearing shown in Figure (b) is selected to tune the frequency of the isolated top-story structure to the first modal frequency of the substructure (i.e., 9-story building structure without a top story). The first five natural periods of the 0-story building are.00, 0.67, 0.40, 0.9 and 0.3 sec. C and C are the damping matrices of the general building and the top-story isolated building, respectively, and are defined using Rayleigh damping. Classical Rayleigh damping is assumed for the example building with a % critical damping ratio in the first two modes. Optimal parameters for the design of a passive TMD proposed by Warburton [] are employed as shown in equation 3 for passive top-story isolation system. In this equation, the parameter refers to the mass ratio of TMD with respect to the main structure, and is 5.6%. The optimal damping ratio for the passive top-story isolation system is calculated to be.5% using equation (3). ( - /4) ( 3 ) 4( )( - /) As shown in Figure, the Bouc-Wen model [] was used to describe how the

98 Hyun-Su Kim damping force is related to the velocity and applied command voltage. The MR damper model has a maximum generated force of about 500 kn, depending on the relative velocity across the MR damper with a saturation voltage of 5 V. 600 passive-on (5V) MR Damper Force (kn) 400 00 0-00 -400 passive-off (0V) -600-40 -0 0 0 40 Velocity (cm/sec) Figure. Force-velocity relationship of MR damper In numerical analysis, the model of the example structure is subjected to an artificial ground motion. An artificial earthquake record was created to fit the design response spectrum, representing the seismic characteristics in the regions of low-to-moderate seismicity using the SIMQKE software program [3]. Figure 3 shows the acceleration time history graph of the artificial earthquake. The peak ground acceleration (PGA) of the artificially-generated earthquake was 0.44g, which suggests that the artificial earthquake represents the ground motions for regions of low-to-moderate seismic activity. The artificial ground motion was generated with a time step of 0.0 seconds and its duration was 30 seconds. 0. Acceleration (g) 0. 0-0. -0. 0 0 0 30 Time (sec) Figure 3. Acceleration time history of an artificial earthquake 3 Development of control algorithm for smart top-story isolation As explained previously, the FLC is used to control the smart top-story isolation system. In the FLC, the 9th floor acceleration and the isolator drift are selected for two input variables and the output variables are the command voltages sent to the MR damper. Fuzzy rules and membership functions are opti-

Seismic response reduction of a building 983 mized using a multi-objective genetic algorithm (MOGA). The reduction of the peak drift of the isolator and the RMS acceleration of the 9th floor are selected as two objectives in a multi-objective optimization process. Each response controlled by the MOGA-optimized FLC is normalized by the corresponding response of the building with a passive top-story isolation system in each objective function. Therefore, if the calculated objective function values are less than, the control performance of the smart top-story isolation system is better than that of the passive top-story isolation system. The optimization process of MOGA is presented in Figure 4. Set initial parameters (population size, final generation, number of variables, range of variables, rate of mutation and crossover, types of selection, mutation, crossover operator) Generate initial population Select one chromosome in current population Decoding chromosome & Make analytical model of top-story isolation system Run time history analysis with earthquake load Evaluate multi-objective functions Evaluated all chromosome? NO Figure 4. MOGA optimization process 4 Control performance evaluation of top-story isolation system A numerical model of the 0-story example building structure with a smart top-story isolation system is implemented in SIMULINK and MATLAB. The MOGA based optimization is performed with the population size of 00 individuals. The upper limit on the number of generations is specified to be 000. As the number of generations increases, the control performance of the elite (i.e. non-dominated) individuals is improved. After the optimization run, the Pareto-optimal front (a set of Pareto-optimal solutions) is obtained as shown in YES Apply selection, mutation, crossover operators Create the next generation

984 Hyun-Su Kim Figure 5. Figure 5 shows that the J and J of each solution are less than ; the smart top-story isolation system can therefore provide better control performance in reducing both the isolation story drift and the 9 th story acceleration compared to the passive system. Consequently, one controller, that can appropriately control both the isolator drift and the 9th floor acceleration responses, has been selected among the Pareto optimal FLCs. J (Normalized 9 th story acceleration) 0.93 0.9 0.9 0.9 J=0.38 J=0.90 0.8 0.9 0.3 0.3 0.3 J (Nomalized isolation story drift) Figure 5. Pareto optimal solutions The values of J and J of the selected FLC (marked as a solid circle in Figure 5) are 0.38 and 0.90, respectively. Accordingly, the smart top-story isolation system controlled by the selected FLC can reduce the peak drift of the isolator and RMS 9th story acceleration by approximately 68% and 0%, respectively, compared to the passive top-story isolation system. The peak displacement of the 0th story with a passive isolation system is. cm as shown in Figure 6. This can be reduced by more than 50% (i.e. 0.9 cm) when a smart isolation system is used. Because the RMS accelerations of the isolated story with the smart system are also less than those of the passive system, it is expected that the smart control system increases the safety of the building structure with a top-story isolation system. The smart top-story isolation system can effectively reduce both the peak displacements and the RMS accelerations of the main structure below the isolation story as shown in Figure 7. Passive Isolation. 7.8 General Building 9.69 66.7 Passive Isolation 8.99 63.6 Smart Isolation 0.9 55.8 Smart Isolation 7.94 57.3 0 0 0 30 0 0 40 60 80 Peak displacement (cm)rms acceleration (cm/s ) Figure 6. 0th story responses 0 4 6 8 0 50 55 60 65 70 Peak displacement (cm) RMS Acceleration (cm/s ) Figure 7. 9th story responses

Seismic response reduction of a building 985 The isolator drift and the 9 th story acceleration time histories of the passive and smart isolation systems are presented in Figures 8 and 9, respectively. The isolation story drift, which is closely related to the safety of a top-story isolation system, can be considerably reduced by the smart isolation system. If additional damping devices are used in the passive isolation system, the isolation story drift can be reduced. However, the addition of damping may deteriorate control performance of the top-story isolation system. Consequently, a smart top-story isolation system can effectively reduce the seismic responses of both the main structure and the isolated top story. Isolation story dift (m) 0. 0. 0-0. -0. 0 0 0 30 Time (sec) Figure 8. Isolator drift Passive Isolation System Smart Isolation System 9 th Story Acceleration (m/s ) 3 0 - - -3 Passive Isolation System Smart Isolation System 0 0 0 30 Time (sec) Figure 9. 9th story acceleration 5 Conclusions In this study, a smart top-story isolation system consisting of an MR damper and low damping elastomeric bearings was proposed to reduce the seismic responses of a building structure. A fuzzy logic control algorithm was developed to effectively control the smart top-story isolation system using a multi-objective genetic algorithm. Both the isolator drift and acceleration response of the example structure were used as objective functions for the multi-objective optimization problem. The optimally-designed passive system was used as a reference for a comparative study. Based on numerical simulations, it can be seen that the smart top-story isolation system can effectively reduce both the main structure and the isolated top story responses compared to the passive top-story isolation system. The MOGA-optimized FLC shows good performance for controlling smart top-story isolation. After a single optimization run using MOGA, an engineer can simply select another FLC that satisfies the desired performance requirements from among a number of Pareto optimal solutions. This may be the most important characteristic of the MOGA based optimization compared to other optimization methods. It is expected that a top-story isolation system is an effective means for mitigating the seismic-induced responses of a tall building. Acknowledgements. This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korea government (MEST) (NRF-03RAA0583).

986 Hyun-Su Kim References [] F. Naeim and J.M. Kelly, Design of Seismic Isolated Structures: from Theory to Practice. Wiley: New York. (999). [] J.M. Kelly, The current state of base isolation in The United States, Proc. nd World Conference on Structural Control, Kyoto, Japan, (999), 043-05. [3] M. Ziyaeifar and H. Noguchi, Partial mass isolation in tall buildings, Earthquake Engineering and Structural Dynamics, 7 (998), 49-65. [4] X. Huang, F. Zhou, S. Wang and X. Luo, Theoretical and experimental investigation on mid-story seismic isolation structures, The 4th World Conference on Earthquake Engineering, Beijing, China, (008). [5] R. Villaverde, Implementation study of aseismic roof isolation system in 3-story building, Journal of Seismology and Earthquake Engineering, (000), 7-7. [6] B.F. Spencer, Jr. and S. Nagarajaiah, State of the art of structural control, Journal of Structural Engineering, ASCE, 9 (003), 845-856. [7] B.F. Spencer, Jr., E.A. Johnson and J.C. Ramallo, Smart isolation for seismic control, JSME International Journal Series C, 43, (000), 704-7. [8] H.S. Kim and P.N. Roschke, GA-fuzzy control of smart base isolated benchmark building using supervisory control technique, Advances in Engineering Software, 38 (007), 453-465. [9] J. N. Yang, J. C. Wu, A. M. Reinhorn and M. Riley, Control of sliding-isolated buildings using sliding-mode control, Journal of Structural Engineering, (996), 79-86. [0] L.M. Jansen and S.J. Dyke, Semi-active control strategies for MR dampers: a comparative study", Journal of Engineering Mechanics, ASCE, 6 (000), 795-803. [] G.B. Warburton, Optimum absorber parameters for various combinations of response and excitation parameters, Earthquake Engineering and Structural Dynamics, 0 (98), 38-40. [] R.H. Sues, S.T. Mau and Y.K. Wen, System Identification of Degrading Hysteretic Restoring Forces, Journal of Engineering Mechanics, ASCE, 4, (988), 833-846. [3] E.H. Vanmarcke and D.A. Gasparini, SIMQKE: A program for artificial motion generation, Dept. of Civil Engineering, Massachusettes Institute of Technology, (976). Received: September, 04