Innovative Technologies for Reducing Structural Vibrations due to Natural Events and Human Activities
View this Special IssueResearch Article  Open Access
Xiangxiu Li, Ping Tan, Xiaojun Li, Aiwen Liu, "Mechanism Analysis and Parameter Optimization of MegaSubIsolation System", Shock and Vibration, vol. 2016, Article ID 2605839, 12 pages, 2016. https://doi.org/10.1155/2016/2605839
Mechanism Analysis and Parameter Optimization of MegaSubIsolation System
Abstract
The equation of motion of megasubisolation system is established. The working mechanism of the megasubisolation system is obtained by systematically investigating its dynamic characteristics corresponding to various structural parameters. Considering the number and location of the isolated substructures, a procedure to optimally design the isolator parameters of the megasubisolation system is put forward based on the genetic algorithm with base shear as the optimization objective. The influence of the number and locations of isolated substructures on the control performance of megasubisolation system has also been investigated from the perspective of energy. Results show that, with increase in substructure mass, the working mechanism of the megasubisolation system is changed from tuned vibration absorber and energy dissipation to seismic isolation. The locations of the isolated substructures have little influence on the optimal frequency ratio but have great influence on the optimal damping ratio, while the number of isolated substructures shows great impact on both the optimal frequency ratio and damping ratio. When the number of the isolated substructures is determined, the higher the isolated substructures, the more the energy that will be consumed by the isolation devices, and with the increase of the number of isolated substructures, the better control performance can be achieved.
1. Introduction
With the development of economics and the advancement of science and technology, more and more tall buildings and super tall buildings have been built in recent years. The most significant engineering concerns in construction of tall and super tall buildings are the safety of the building structures and the comfort of occupants under external forces such as winds and earthquakes. Isolation technology is known as an effective way to improve the structural seismic response. The objective of seismic isolation is to decouple the structure from the ground motion, preventing the structure from absorbing the earthquake energy. Recently some significant research has been done on the isolation structures. Khoshnoudian and Mehrparvar [1] proposed a new control algorithm to protect nonlinear baseisolated structures against earthquakes. Intensity measures for the seismic response prediction of baseisolated buildings were investigated by Mollaioli et al. [2]. Leopa et al. [3] have done some research studies on damage identification in passive vibroisolation devices, which dealt with the theoretical aspects combined with experimental analysis regarding early damage identification in passive vibroisolation devices. Briman and Ribakov [4] have developed a method for building retrofits, that is, replacing weak conventional columns with low static and dynamic loadcarrying capacities by seismic isolation columns. The effectiveness of various control strategies in hybrid base isolation systems including isolators and semiactive variable friction dampers has been investigated by Ribakov [5]. The disorder and damage of baseisolated medical facilities when subjected to nearfault and longperiod ground motions have been investigated by Shi et al. [6]. Yan and Chen [7] have studied the seismic performance of midstory isolated structures under nearfield pulselike ground motion, in which a mechanical model has been put forward for this protective system based on the Kelvin pounding model, and a new method has been proposed that synthesizes artificial nearfield pulselike ground motion by combining the real nearfield nonpulse ground motion with simple equivalent pulses. Ozdemir and Akyuz [8] have analyzed the dynamic responses of isolated structures under bidirectional excitations of nearfield ground motions, in which nonlinear response history analyses of a 3story isolated reinforced concrete building have been carried out under both uni and bidirectional earthquake excitations of nearfield records. A comparative analytical study of several control strategies for semiactive devices installed in baseisolated buildings aiming at reducing earthquake induced vibrations is presented by Oliveira et al. [9]. Castaldo et al. [10] have studied the seismic reliability of a baseisolated structure with friction pendulum isolators considering both isolator properties and earthquake main characteristics as random variables.
For buildings of modest height, implementation of passive control devices offers a potential improvement in structural safety and human comfort. But the structural characteristics common to most tall and super tall buildings, such as high shear rigidity, tend to prevent the application of the traditional control devices. It is still a serious challenge for structural engineers to effectively reduce the structural responses of tall and super tall buildings to further improve the structural safety. A new method for controlling the response of tall and super tall buildings under severe external loads was first introduced by Feng and Mita [11]. Feng and Mita first proposed releasing the connections between the megastructure and the substructures in a megasubstructure, but without installing dampers between the megastructure and the substructures. Chai and Feng [12] subsequently improved this configuration and presented a megasubcontrolled system based on a conventional megasubframe and undertook a study of its dynamic response to random wind load excitations. Recently, some studies on the optimal parameters between the substructure and megastructure have been done in order to achieve the best performance by Tian [13]. Lan et al. proposed a multifunction megasubcontrolled structure; this structure has the function of the mass dampers and base isolation as well as damping energy dissipation. However, there is a lack of theoretical certification [14]. A new connection form between the substructures and megastructure was put forward by Pei and Wang [15], in which the top substructure was connected with the megastructure by dampers. And the studies showed that this new connection form can achieve better damping effect and also can prevent collisions between the top substructures and the megastructure. Qin et al. [16, 17] proposed a new control method, which employs active control and passive dampers together to form a huge control system based on the particular conformation characteristic of megasubcontrolled system, and also presented a new kind of structural configuration, named megasubcontrolled structure, which was constructed by applying the structural control principle to structural configuration itself, to form a new structure with obvious response selfcontrol ability, instead of employing the conventional method. Xun’an et al. [18, 19] have investigated the control performances of the megasubcontrolled structure with different control strategies. Lian et al. [20, 21] have analyzed the seismic responses of the megasubcontrolled frame with friction damper, in which friction damped structure is first designed as the substructure in the megasubcontrolled frame to further reduce structural dynamical responses. A more reasonable and realistic scaled model has been designed to investigate the dynamical characteristics and controlling performances of the megasubcontrolled system when subjected to strong earthquake motion by Limazie et al. [22], and the control parameters of the structure system, such as the modulated substructures disposition, the damping coefficient ratio, the stiffness ratio, and the mass ratio of the megastructure and substructure, have been investigated. A parametric study of the relative stiffness ratio and relative mass ratio between the megaframe and the substructures, as well as the additional column stiffness ratio that influences the response control effectiveness of the megasubcontrolled structure, was discussed in the literature [23], but the number and locations of the isolated substructures were not considered.
The research of the control mechanism and parameter optimization of the megasubisolation system are insufficient in the past few years. In this paper, the working mechanism of the megasubisolation system is investigated by systematically studying its dynamic characteristics corresponding to various structural parameters. By considering the number and locations of the isolated substructures, the parameter optimization of megasubisolation system has been studied based on the genetic algorithm theory where base shear was used as the optimization objective function. The influence of the number and locations of isolated substructures on the control performance of megasubisolation system has also been investigated from the perspective of energy.
2. Control Mechanism Study of MegaSubIsolation System
2.1. Equations of Motion of MegaSubIsolation System
In a megasubisolation system, substructures are connected with the megastructure by isolation devices. For the analysis purposes, the isolation devices are simulated with Kelvin model, and the megastructure is simplified as a series of pointed models, while each substructure is simplified as a lumped mass. Figure 1 shows the simplified analysis model of the megasubisolation system.
According to the D’Alembert principle, motion equation of the analysis model shown in Figure 1 can be established as follows.
For substructure where is the equivalent mass of each substructure; and are the equivalent damping and horizontal stiffness of the isolation device, respectively; , , and are the acceleration, velocity, and displacement of the substructure relative to the ground; , , and are the acceleration, velocity, and displacement of the megastructure relative to the ground; is the ground motion acceleration.
For megastructure in which represents the unit vector; is the force vector (substructures apply to the megastructures); , , and are the equivalent mass, equivalent damping, and equivalent stiffness of each megastructure, respectively. There is no substructure installed at the first floor of the megastructure because the bottom story of substructure is supported on the ground, and its movement will not directly influence the structural responses. So the number of the substructures is fewer than that of the megastructures by one. The mass matrix, stiffness matrix of the megastructure, and the force vector can be expressed as follows, and the damping matrix of the megastructure is determined by the Rayleigh damping:
According to (1) and (2), the motion equation of the analysis model shown in Figure 1 is given aswhere , and the mass matrix, stiffness matrix, and damping matrix of the megasubisolation system can be obtained as follows:
2.2. Control Mechanism Analysis of MegaSubIsolation System
A typical project is selected as an example [14], which is composed of five megastories and the substructures are attached to megastructure from the second floor to the fifth floor. The mass and shear stiffness of each megastructure are kg and N/m, respectively, and the mass of top megastructure is kg. The mass of each substructure is determined by the mass ratio , the frequency ratio is defined as the ratio of the substructure’s frequency to the firstorder frequency of the megastructure, and each substructure has the same parameter values. When the megastructure is simplified as a series of particlebased models, its first period is 2 s, and the fundamental period of the structure is 2.8 s when the substructures are rigidly connected with the megastructures.
The modal analysis results of the structure when the mass ratio and frequency ratio are shown in Table 1 and Figure 2. It can be observed from Table 1 that the first mode and the fifth mode are the main vibration modes of the structure, and these two vibration modes are the global motion of the structure, while the others are local modes of the substructure or megastructure, which can also be seen from Figure 2. Therefore we believe that the control mechanism of the megasubisolation system is related to these two vibration modes.

Figure 3 shows the relationship between the modal mass participation factors and mass ratio and the relationship between the modal mass participation factors and frequency ratio. From Figure 3, it can be confirmed that the first vibration mode and the fifth vibration mode are the main modes of the structure and it is found that when the frequency ratio increases, the modal mass participation factor of the first vibration mode will increase, while that of the fifth vibration mode will decrease.
In order to study the control mechanism of the megasubisolation system, a 2DOF equivalent model is adopted based on the modal analysis results. The mass of the bottomstructure in the equivalent model is the total mass of the megastructures, the firstorder equivalent stiffness of the megastructure is adopted as the stiffness of the bottomstructure, and the mass of the superstructure is determined by the mass ratio . The connection between the superstructure and bottomstructure is simulated by a Kelvin model.
From Figure 4, it can be seen that the relationship between the modal mass participation factors and frequency ratio of the megasubisolation system is consistent with that of the 2DOF equivalent model when the mass ratio is taken as , , and . Then, the 2DOF equivalent model can be used to explain the control mechanism of the megasubisolation system; that is, the working mechanism of the megasubisolation system is changed from tuned vibration absorber and energy dissipation to seismic isolation with increase in substructure mass, according to the literature [24], in which the working mechanism of the 2DOF equivalent model has been systematically investigated.
3. Parameters Optimization
Base shear is very important in earthquake resistant design of building structures. As for isolation system, changing the stiffness and damping of isolators will shift the dynamic characteristic of the system and the distribution of the story shear in the system. In this paper, a procedure to optimally design the isolator parameters of the megasubisolation system is put forward using the genetic algorithm based on the minimum base shear variance criterion.
The base shear of the megasubisolation system is given bywhere and are the damping and stiffness of the bottom megastructure, respectively. Assuming earthquake excitation to be a stationary stochastic process with zero mean, the responses can also be considered as two statistically stationary processes with zero mean. According to superposition principle of a stationary stochastic process, the following formulas can be achieved:where “−” represents a mean value of a variable and means taking the real part of complex number. Let represent the power spectral density function of seismic excitation and let , represent the transfer functions of the displacement response and velocity response of the bottom megastructure, respectively, which can be obtained by complex modal analysis methods [25]. The mean square value of the velocity response, the displacement response and the mean value of product can be evaluated as follows by stochastic vibration method:Noting that is the complex conjugate function, then only the real part of contributes to the integration in (10). The power spectral function of the base shear of the megasubisolation system can be obtained by the formula . The base shear variance is given by In order to minimize the base shear of the megasubisolation system, the optimal isolator parameters such as frequency ratio and isolator damping ratio can be solved by a nonlinear mathematical programming method: where and represent the lower and upper limit of the frequency ratio which are taken as 0.01 and 1.2, respectively, and , represent the lower and upper limit of the damping ratio which are taken as 0.01 and 0.3, respectively. The damping ratio of the megastructure is supposed as 0.05. A modified KanaiTajimi model suggested by Hu Yuxian is employed as the stochastic ground motion model, and the parameters of , , , and are taken as 15.74 cm^{2}/(rads^{3}), 19.07 rad/s, 0.784, and 2.015 rad/s, respectively. Figure 5 shows the optimization process by genetic algorithm method and Figure 6 shows the optimization results when the mass ratio is taken as 1.
(a) Case of parameter aggregation
(b) Target value evolution
As can be seen from Figure 6(b), the convergence rate of the genetic algorithm is fast, and an optimal design can be obtained in less than 50 generations. From Figure 6(a), it can be seen that the optimal values of the frequency ratio and damping ratio are 0.44 and 0.3, respectively, when the mass ratio is taken as 1.
The optimal results with the changes of the number and location of the isolated substructures can also be obtained by this method, which are shown in Table 2.

The numbers shown in Table 2 represent the isolated floors of the substructures; for example, number 2 indicates that the second floor of the substructure is connected with the megastructure by isolation device, and the other substructures are rigidly connected with the megastructures; numbers 4 and 5 indicate that the fourth and the fifth floors of the substructures are connected with the megastructure by isolation device, while the other substructures are rigidly connected with the megastructures. From Table 2, it can be seen that the locations of isolated substructures have little influence on the optimal frequency ratio but have great influence on the optimal damping ratio; when the location of the isolated substructures is lower, the optimal damping ratio shows a trend of decreasing, while the number of isolated substructures shows great impact on both the optimal frequency ratio and damping ratio. With the increasing number of the isolated substructures, the optimal frequency ratio decreases, and the optimal damping ratio increases.
Meanwhile, the mass ratio shows great impact on both the optimal frequency ratio and optimal damping ratio; for the same arrangement scheme of the isolated substructures, the optimal frequency ratio decreases, and the optimal damping ratio increases with the increase of mass ratio.
After the parameters optimization, the power spectral density for displacement of top megastructure and the power spectral density for acceleration of top substructure are compared to those before optimization, as shown in Figures 7 and 8. From Figures 7 and 8, it can be seen that the structural responses after parameter optimization can be reduced obviously compared with that before parameter optimization, so these optimization results can give useful information to designers.
4. Energy Equation
In accordance with the energy balance principle, the energy of each part of structure can be obtained by taking integration for the relative displacement of (4) on both sides at the same time as follows:in which the total kinetic energy of the structure is defined by , the damping energy is , and the elastic strain energy of the structure is defined by ; at the same time, is defined as the total input energy, and is defined as the transfer energy.
5. Influence of the Number and Locations of Isolated Substructures on Control Performance
The number and locations of isolated substructures are important factors influencing the damping effect of the megasubisolation system. Ten layout programs of the isolated substructures are adopted in this paper and, based on the optimization results obtained from Section 2, the maximum displacement of top megastructure, the maximum acceleration of top substructure, the total input energy of the system, and the transfer energy of the substructure have been calculated when the mass ratio is taken as 1, which are shown in Tables 3 and 4. Both the El Centro ground motion and Taft ground motion are employed for the seismic excitation and the PGA of them are taken as 0.4 g.


It can be seen from Tables 3 and 4 that when the number of the isolated substructures is given, locations of isolated substructures have little effect on the maximum displacement of the top megastructure, the maximum acceleration of the top substructure, and the total input energy, but they do have great impact on the transfer energy of substructure, and the transfer energy of substructure is larger when the substructures are located in the upper part of the system. The larger transfer energy of substructure indicates that the more energy consumed by the isolators, the less energy consumed by the megastructure itself when the amount of the total input energy is determined, and thus the safety of the megastructure can be enhanced.
The number of the isolated substructures shows great impact on the control performance of the megasubisolation system. With the increase of the number of isolated substructures, the maximum displacement of the top megastructure, the maximum acceleration of the top substructure, and the total input energy decrease, while the transfer energy increases. Figure 9 shows the average values of maximum displacement of top megastructure, maximum acceleration of top substructure, the total input energy of the system, and the transfer energy of the substructure with the changes of the number of isolated substructures. The number in axis represents the number of the isolated substructures. From Figure 9, it can be confirmed that the number of isolated substructures can influence the control performance of the megasubisolation system greatly and, with the increasing number of isolated substructures, the average values of maximum displacement of top megastructure and the total input energy show a decreasing trend, while the average value of the transfer energy presents the opposite trend. That is, when the number of isolated substructures increases, the total input energy can be reduced and larger energy can be transferred to the substructure which is consumed by the isolation devices, and then the energy consumed by the megastructure itself becomes less. Therefore, in practical engineering applications, all the substructures should be connected with the megastructures by isolation devices in order to achieve the best control performance.
(a) Comparison of the average displacement
(b) Comparison of the average acceleration
(c) Comparison of the average input energy
(d) Comparison of the average transfer energy
6. Conclusions
The control mechanism of the megasubisolation system has been studied from its dynamic characteristics in this paper, and the parameters optimization and the control performance of a megasubisolation system have also been studied when the number and locations of the isolated substructures are different. The following conclusions can be made.
With increasing substructure mass, the working mechanism of the megasubisolation system changes from tuned vibration absorber and energy dissipation to seismic isolation, which is confirmed by theoretical analysis.
Locations of isolated substructures have little influence on the optimal frequency ratio and have great influence on the optimal damping ratio, while the number of isolated substructures and the mass ratio show great impact on both the optimal frequency ratio and damping ratio. The results of parameter optimization can give some useful design suggestions to this novel system.
When the number of isolated substructures is given, locations of the isolated substructures have little effect on the maximum displacement of the top megastructure, the maximum acceleration of the top substructure, and the total input energy, but they have great impact on the transfer energy of substructure, and the transfer energy of substructure is larger when the isolated substructures are located in the upper part of the system. That is, when the isolated substructures are located in the upper part of the structure, more energy will be consumed by the isolation devices. With the increase of the number of isolated substructures, the maximum displacement of top main structure and the total input energy show a decreasing trend, but the transfer energy presents the opposite trend. That is, the more the isolated substructures are, the better the control performance can be achieved. Hence, in practical engineering applications, when the number of isolated substructures is determined, the isolation devices should be installed in the upper substructures and, in order to achieve the best control performance, all the substructures should be connected to the megastructures with isolation devices.
Competing Interests
The authors declare that there are no competing interests regarding the publication of this paper.
Acknowledgments
The work was primarily supported by the National Basic Research Program of China (2011CB013601) and the Major Research Plan of the National Natural Science Foundation of China under Grant no. 9131530107, partially by the National Science Foundation under Grants nos. 51408560 and 51578514, and also in part by the Basic Research and Business Projects of Public Welfare Institutes under Central Level by Grant no. DQJB15B11.
References
 F. Khoshnoudian and B. Mehrparvar, “A new control algorithm to protect nonlinear baseisolated structures against earthquakes,” The Structural Design of Tall and Special Buildings, vol. 22, no. 18, pp. 1376–1389, 2013. View at: Publisher Site  Google Scholar
 F. Mollaioli, A. Lucchini, Y. Cheng, and G. Monti, “Intensity measures for the seismic response prediction of baseisolated buildings,” Bulletin of Earthquake Engineering, vol. 11, no. 5, pp. 1841–1866, 2013. View at: Publisher Site  Google Scholar
 A. Leopa, S. Nastac, and C. Debeleac, “Researches on damage identification in passive vibroisolation devices,” Shock and Vibration, vol. 19, no. 5, pp. 803–809, 2012. View at: Publisher Site  Google Scholar
 V. Briman and Y. Ribakov, “Using seismic isolation columns for retrofitting buildings with soft stories,” The Structural Design of Tall and Special Buildings, vol. 18, no. 5, pp. 507–523, 2009. View at: Publisher Site  Google Scholar
 Y. Ribakov, “Baseisolated structures with selective controlled semiactive friction dampers,” The Structural Design of Tall and Special Buildings, vol. 20, no. 7, pp. 757–766, 2011. View at: Publisher Site  Google Scholar
 Y. D. Shi, M. Kurata, and M. Nakashima, “Disorder and damage of baseisolated medical facilities when subjected to nearfault and longperiod ground motions,” Earthquake Engineering and Structural Dynamics, vol. 43, no. 11, pp. 1683–1701, 2014. View at: Publisher Site  Google Scholar
 G. Y. Yan and F. Q. Chen, “Seismic performance of midstory isolated structures under nearfield pulselike ground motion and limiting deformation of isolation layers,” Shock and Vibration, vol. 2015, Article ID 730612, 17 pages, 2015. View at: Publisher Site  Google Scholar
 G. Ozdemir and U. Akyuz, “Dynamic analyses of isolated structures under bidirectional excitations of nearfield ground motions,” Shock and Vibration, vol. 19, no. 4, pp. 505–513, 2012. View at: Publisher Site  Google Scholar
 F. Oliveira, P. Morais, and A. Suleman, “A comparative study of semiactive control strategies for base isolated buildings,” Earthquake Engineering and Engineering Vibration, vol. 14, no. 3, pp. 487–502, 2015. View at: Publisher Site  Google Scholar
 P. Castaldo, B. Palazzo, and P. Della Vecchia, “Seismic reliability of baseisolated structures with friction pendulum bearings,” Engineering Structures, vol. 95, pp. 80–93, 2015. View at: Publisher Site  Google Scholar
 M. Q. Feng and A. Mita, “Vibration control of tall buildings using mega subconfiguration,” Journal of Engineering Mechanics, vol. 121, no. 10, pp. 1082–1088, 1995. View at: Publisher Site  Google Scholar
 W. Chai and M. Q. Feng, “Vibration control of super tall buildings subjected to wind loads,” International Journal of NonLinear Mechanics, vol. 32, no. 4, pp. 657–668, 1997. View at: Publisher Site  Google Scholar
 Y. J. Tian, Vibration Damping Theory and Optimization Design Method for Multifunctional VibrationAbsorption of RC Megaframe Structures, College of Civil Engineering, Southeast University, Nanjing, China, 2002.
 Z. J. Lan, Y. J. Tian, S. Y. Cao, and H. H. Wang, “The mechanism and efficiency of multifunctional vibration absorption for megaframe structures,” China Civil Engineering Journal, vol. 35, no. 6, pp. 1–5, 2002. View at: Google Scholar
 X.Z. Pei and L. Wang, “Earthquake response analysis of megasuspended braced frame attached with damper,” Journal of Vibration and Shock, vol. 30, no. 11, pp. 190–197, 2011. View at: Google Scholar
 X. J. Qin, X. A. Zhang, L. Zou, Q. Q. Wang, and P. Jiao, “Research on active control of megasub controlled structure with dampers subjected to seismic loads,” in International MultiConference of Engineers and Computer Scientists (IMECS '10), pp. 1055–1060, March 2010. View at: Google Scholar
 X. J. Qin and X. A. Zhang, “Vibration controlling characteristics of a new improved megasub controlled structure subjected to random wind loads,” in Proceedings of the International MultiConference of Engineers and Computer Scientists (IMECS '09), pp. 1342–1347, Hong Kong, March 2009. View at: Google Scholar
 Z. Xun'an, W. Dong, and J. Jiesheng, “The controlling mechanism and the controlling effectiveness of passive megasubcontrolled frame subjected to random wind loads,” Journal of Sound and Vibration, vol. 283, no. 3–5, pp. 543–560, 2005. View at: Publisher Site  Google Scholar
 X. A. Zhang, X. J. Qin, S. Cherry, Y. D. Lian, J. L. Zhang, and J. S. Jiang, “A new proposed passive megasub controlled structure and response control,” Journal of Earthquake Engineering, vol. 13, no. 2, pp. 252–274, 2009. View at: Publisher Site  Google Scholar
 Y. D. Lian, X. A. Zhang, and C. Sheldon, “Damping characteristics of friction damped braced frame and its effectiveness in the megasub controlled structure system,” Earthquake Engineering and Engineering Vibration, vol. 6, no. 2, pp. 171–181, 2007. View at: Publisher Site  Google Scholar
 Y. D. Lian, X. J. Wang, and T. Li, “Analysis of seismic response for the megasub controlled frame with friction damper,” in Chapter 6: Landscape Planning and Design, vol. 450451 of Advanced Materials Research, pp. 1269–1273, 2012. View at: Publisher Site  Google Scholar
 T. Limazie, X. A. Zhang, and X. J. Wang, “Vibration control parameters investigation of the megasub controlled structure system (MSCSS),” Earthquake and Structures, vol. 5, no. 2, pp. 225–237, 2013. View at: Publisher Site  Google Scholar
 X. J. Qin, X. A. Zhang, and S. Cherry, “Study on semiactive control of megasub controlled structure by MR damper subject to random wind loads,” Earthquake Engineering and Engineering Vibration, vol. 7, no. 3, pp. 285–294, 2008. View at: Publisher Site  Google Scholar
 F. L. Zhou, Y. Zhang, and P. Tan, “Theoretical study on story isolation system,” China Civil Engineering Journal, vol. 42, no. 8, pp. 1–8, 2009. View at: Google Scholar
 J. P. Ou and G. Y. Wang, Structural Random Vibration, High Education Press, Beijing, China, 1998.
Copyright
Copyright © 2016 Xiangxiu Li et al. 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.