首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
This paper summarizes the results of an extensive study on the inelastic seismic response of X‐braced steel buildings. More than 100 regular multi‐storey tension‐compression X‐braced steel frames are subjected to an ensemble of 30 ordinary (i.e. without near fault effects) ground motions. The records are scaled to different intensities in order to drive the structures to different levels of inelastic deformation. The statistical analysis of the created response databank indicates that the number of stories, period of vibration, brace slenderness ratio and column stiffness strongly influence the amplitude and heightwise distribution of inelastic deformation. Nonlinear regression analysis is employed in order to derive simple formulae which reflect the aforementioned influences and offer a direct estimation of drift and ductility demands. The uncertainty of this estimation due to the record‐to‐record variability is discussed in detail. More specifically, given the strength (or behaviour) reduction factor, the proposed formulae provide reliable estimates of the maximum roof displacement, the maximum interstorey drift ratio and the maximum cyclic ductility of the diagonals along the height of the structure. The strength reduction factor refers to the point of the first buckling of the diagonals in the building and thus, pushover analysis and estimation of the overstrength factor are not required. This design‐oriented feature enables both the rapid seismic assessment of existing structures and the direct deformation‐controlled seismic design of new ones. A comparison of the proposed method with the procedures adopted in current seismic design codes reveals the accuracy and efficiency of the former. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

2.
This paper presents the main results of the evaluation of residual inter‐story drift demands in typical moment‐resisting steel buildings designed accordingly to the Mexican design practice when subjected to narrow‐band earthquake ground motions. Analytical 2D‐framed models representative of the study‐case buildings were subjected to a set of 30 narrow‐band earthquake ground motions recorded on stations placed in soft‐soil sites of Mexico City, where most significant structural damage was found in buildings as a consequence of the 1985 Michoacan earthquake, and scaled to reach several levels of intensity to perform incremental dynamic analyses. Thus, results were statistically processed to obtain hazard curves of peak (maximum) and residual drift demands for each frame model. It is shown that the study‐case frames might exhibit maximum residual inter‐story drift demands in excess of 0.5%, which is perceptible for building's occupants and could cause human discomfort, for a mean annual rate of exceedance associated to peak inter‐story drift demands of about 3%, which is the limiting drift to avoid collapse prescribed in the 2004 Mexico City Seismic Design Provisions. The influence of a member's post‐yield stiffness ratio and material overstrength in the evaluation of maximum residual inter‐story drift demands is also discussed. Finally, this study introduces response transformation factors, Tp, that allow establishing residual drift limits compatible with the same mean annual rate of exceedance of peak inter‐story drift limits for future seismic design/evaluation criteria that take into account both drift demands for assessing a building's seismic performance. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

3.
The effect of peak ground velocity (PGV) on single‐degree‐of‐freedom (SDOF) deformation demands and for certain ground‐motion features is described by using a total of 60 soil site records with source‐to‐site distances less than 23 km and moment magnitudes between 5.5 and 7.6. The observations based on these records indicate that PGV correlates well with the earthquake magnitude and provides useful information about the ground‐motion frequency content and strong‐motion duration that can play a role on the seismic demand of structures. The statistical results computed from non‐linear response history analyses of different hysteretic models highlight that PGV correlates better with the deformation demands with respect to other ground motion intensity measures. The choice of PGV as ground motion intensity decreases the dispersion due to record‐to‐record variability of SDOF deformation demands, particularly in the short period range. The central tendencies of deformation demands are sensitive to PGV and they may vary considerably as a function of the hysteretic model and structural period. The results provided in this study suggest a consideration of PGV as a stable candidate for ground motion intensity measure in simplified seismic assessment methods that are used to estimate structural performance for earthquake hazard analysis. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

4.
Results of an analytical study aimed at evaluating residual displacement ratios, Cr, which allow the estimation of residual displacement demands from maximum elastic displacement demands is presented. Residual displacement ratios were computed using response time‐history analyses of single‐degree‐of‐freedom systems having 6 levels of relative lateral strength when subjected to an ensemble of 240 earthquake ground motions recorded in stations placed on firm sites. The results were statistically organized to evaluate the influence of the following parameters: period of vibration, level of relative lateral strength, site conditions, earthquake magnitude, and distance to the source. In addition, the influence of post‐yield stiffness ratio in bilinear systems and of the unloading stiffness in stiffness‐degrading systems was also investigated. A special emphasis is given to the uncertainty of these ratios. From this study, it is concluded that mean residual displacement ratios are more sensitive to changes in local site conditions, earthquake magnitude, distance to the source range and hysteretic behaviour than mean inelastic displacement ratios. In particular, residual displacement ratios exhibit large levels of record‐to‐record variability and, therefore, this dispersion should be taken into account when estimating residual displacements. A simplified expression is presented to estimate mean residual displacements ratios for elastoplastic systems during the evaluation of existing structures built on firm soil sites. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

5.
The variance of collapse capacity is an important constituent of probabilistic methodologies used to evaluate the probability of collapse of structures subjected to earthquake ground motions. This study evaluates the effect of ground motion randomness (i.e. record‐to‐record (RTR) variability) and uncertainty in the deterioration parameters of single‐degree‐of‐freedom (SDOF) systems on the variance of collapse capacity. Collapse capacity is evaluated in terms of a relative intensity defined as the ratio of ground motion intensity to a structure strength parameter. The effect of RTR variability on the variance of collapse capacity is directly obtained by performing dynamic analyses of deteriorating hysteretic models for a set of representative ground motions. The first‐order second‐moment (FOSM) method is used to quantify the effect of deterioration parameter uncertainty. In addition to RTR variability, the results indicate that uncertainty in the displacement at the peak (cap) strength and the post‐capping stiffness significantly contribute to the variance of collapse capacity. If large dispersion of these parameters exists, the effect of uncertainty in deterioration parameters on the variance of collapse capacity may be comparable to that caused by RTR variability. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

6.
A process is outlined and evaluated for the estimation of seismic roof and storey drift demands for frame structures from the spectral displacement demand at the first mode period of the structure. The spectral displacement demand is related to the roof drift demand for the multi‐degree‐of‐freedom (MDOF) structure using three modification factors, accounting for MDOF effects, inelasticity effects, and P‐delta effects. Median values and measures of dispersion for the factors are obtained from elastic and inelastic time history analyses of nine steel moment resisting frame structures subjected to sets of ground motions representative of different hazard levels. The roof drift demand is related to the storey drift demands, with the results being strongly dependent on the number of stories and the ground motion characteristics. The relationships proposed in this paper should prove useful in the conceptual design phase, in estimating deformation demands for performance assessment, and in improving basic understanding of seismic behaviour. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

7.
The seismic performance of conventional wood‐frame structures in south‐western British Columbia is analytically investigated through incremental dynamic analysis by utilizing available UBC‐SAWS models, which were calibrated based on experimental test results. To define an adequate target response spectrum that is consistent with information from national seismic hazard maps, record selection/scaling based on the conditional mean spectrum (CMS) is implemented. Furthermore, to reflect complex seismic hazard contributions from different earthquake sources (i.e. crustal events, interface events, and inslab events), we construct CMS for three earthquake types, and use them to select and scale an adequate set of ground motion records for the seismic performance evaluation. We focus on the impacts of adopting different record selection criteria and of using different shear‐wall types (Houses 1–4; in terms of seismic resistance, House 1>House 2>House 3>House 4) on the nonlinear structural response. The results indicate that the record selection procedures have significant influence on the probabilistic relationship between spectral acceleration at the fundamental vibration period and maximum inter‐story drift ratio, highlighting the importance of taking into account response spectral shapes in selecting and scaling ground motion records. Subjected to ground motions corresponding to the return period of 2500 years, House 1 is expected to experience very limited extent of damage; Houses 2 and 3 may be disturbed by minor damage; whereas House 4 may suffer from major damage occasionally. Finally, we develop statistical models of the maximum inter‐story drift ratio conditioned on a seismic intensity level for wood‐frame houses, which is useful for seismic vulnerability assessment. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

8.
This paper presents the results of a multi‐level pseudo‐dynamic seismic test program that was performed to assess the performance of a full‐scale three‐bay, two‐storey steel–concrete composite moment‐resisting frame built with partially encased composite columns and partial‐strength beam‐to‐column joints. The system was designed to develop a ductile response in the joint components of beam‐to‐column joints including flexural yielding of beam end plates and shear yielding of the column web panel zone. The ground motion producing the damageability limit state interstorey drift caused minor damage while the ultimate limit state ground motion level entailed column web panel yielding, connection yielding and plastic hinging at the column base connections. The earthquake level chosen to approach the collapse limit state induced more damage and was accompanied by further column web panel yielding, connection yielding and inelastic phenomena at column base connections without local buckling. During the final quasi‐static cyclic test with stepwise increasing displacement–amplitudes up to an interstorey drift angle of 4.6%, the behaviour was ductile although cracking of beam‐to‐end‐plate welds was observed. Correlations with numerical simulations taking into account the inelastic cyclic response of beam‐to‐column and column base joints are also presented in the paper together. Inelastic static pushover and time history analysis procedures are used to estimate the structural behaviour and overstrength factors of the structural system under study. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

9.
This paper presents an analytical study evaluating the influence of ground motion duration on structural damage of 3‐story, 9‐story, and 20‐story SAC steel moment resisting frame buildings designed for downtown Seattle, WA, USA, using pre‐Northridge codes. Two‐dimensional nonlinear finite element models of the buildings are used to estimate the damage induced by the ground motions. A set of 44 ground motions is used to study the combined effect of spectral acceleration and ground motion significant duration on drift and damage measures. In addition, 10 spectrally equivalent short‐duration shallow crustal ground motions and long‐duration subduction zone records are selected to isolate duration effect and assess its effect on the response. For each ground motion pair, incremental dynamic analyses are performed at at least 20 intensity levels and response measures such as peak interstory drift ratio and energy dissipated are tracked. These response measures are combined into two damage metrics that account for the ductility and energy dissipation. Results indicate that the duration of the ground motion influences, above all, the combined damage measures, although some effect on drift‐based response measures is also observed for larger levels of drift. These results indicate that because the current assessment methodologies do not capture the effects of ground motion duration, both performance‐based and code‐based assessment methodologies should be revised to consider damage measures that are sensitive to duration. Copyright © 2016 John Wiley & Sons, Ltd  相似文献   

10.
The assessment of earthquake loss often requires the definition of a relation between a measure of damage and a quantity of loss, usually achieved through the employment of a damage‐to‐loss model. These models are frequently characterized by a large variability, which inevitably increases the uncertainty in the vulnerability assessment and earthquake loss estimation. This study provides an insight on the development of damage‐to‐loss functions for moment‐frame reinforced concrete buildings through an analytical methodology. Tri‐dimensional finite element models of existing reinforced concrete buildings were subjected to a number of ground motion records compatible with the seismicity in the region of interest, through nonlinear dynamic analysis. These results were used to assess, for a number of damage states, the probability distribution of loss ratio, taking into consideration member damage and different repair techniques, as well as to derive sets of fragility functions. Then, a vulnerability model (in terms of the ratio of cost of repair to cost of replacement, conditional on the level of ground shaking intensity) was derived and compared with the vulnerability functions obtained through the combination of various damage‐to‐loss models with the set of fragility functions developed herein. In order to provide realistic estimates of economic losses due to seismic action, a comprehensive study on repair costs using current Portuguese market values was also carried out. The results of this study highlight important issues in the derivation of vulnerability functions, which are a fundamental component for an adequate seismic risk assessment. © 2015 The Authors. Earthquake Engineering & Structural Dynamics published by John Wiley & Sons Ltd.  相似文献   

11.
This paper presents the seismic behaviour of a concentrically braced frame system with self‐centring capability, in which a special type of bracing element termed reusable hysteretic damping brace (RHDB) is used. The RHDB is a passive energy dissipation device with its core energy‐dissipating component made of superelastic Nitinol wires. Compared with conventional bracing in steel structures, RHDB has a few prominent performance characteristics: damage free under frequent and design basis earthquakes in earthquake prone areas; minimal residual drifts due to the self‐centring capability of RHDB frame; and ability to survive several strong earthquakes without the need for repair or replacement. This paper also includes a brief discussion of the RHDB's mechanical configuration and analytical model for RHDB. The seismic performance study of RHDB frame was carried out through a non‐linear time history analysis of 3‐storey and 6‐storey RHDB frame buildings subjected to two suites of 20 earthquake ground motions. The analysis results were compared with buckling‐restrained brace (BRB) frames. This study shows that RHDB frame has a potential to outperform BRB frames by eliminating the residual drift problem. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

12.
A probabilistic representation of the entire ground‐motion time history can be constructed based on a stochastic model that depends on seismic source parameters. An advanced stochastic simulation scheme known as Subset Simulation can then be used to efficiently compute the small failure probabilities corresponding to structural limit states. Alternatively, the uncertainty in the ground motion can be represented by adopting a parameter (or a vector of parameters) known as the intensity measure (IM) that captures the dominant features of the ground shaking. Structural performance assessment based on this representation can be broken down into two parts, namely, the structure‐specific part requiring performance assessment for a given value of the IM, and the site‐specific part requiring estimation of the likelihood that ground shaking with a given value of the IM takes place. The effect of these two alternative representations of ground‐motion uncertainty on probabilistic structural response is investigated for two hazard cases. In the first case, these two approaches are compared for a scenario earthquake event with a given magnitude and distance. In the second case, they are compared using a probabilistic seismic hazard analysis to take into account the potential of the surrounding faults to produce events with a range of possible magnitudes and distances. The two approaches are compared on the basis of the probabilistic response of an existing reinforced‐concrete frame structure, which is known to have suffered shear failure in its columns during the 1994 Northridge Earthquake in Los Angeles, California. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

13.
In a companion paper, an overview and problem definition was presented for ground motion selection on the basis of the conditional spectrum (CS), to perform risk‐based assessments (which estimate the annual rate of exceeding a specified structural response amplitude) for a 20‐story reinforced concrete frame structure. Here, the methodology is repeated for intensity‐based assessments (which estimate structural response for ground motions with a specified intensity level) to determine the effect of conditioning period. Additionally, intensity‐based and risk‐based assessments are evaluated for two other possible target spectra, specifically the uniform hazard spectrum (UHS) and the conditional mean spectrum (CMS, without variability).It is demonstrated for the structure considered that the choice of conditioning period in the CS can substantially impact structural response estimates in an intensity‐based assessment. When used for intensity‐based assessments, the UHS typically results in equal or higher median estimates of structural response than the CS; the CMS results in similar median estimates of structural response compared with the CS but exhibits lower dispersion because of the omission of variability. The choice of target spectrum is then evaluated for risk‐based assessments, showing that the UHS results in overestimation of structural response hazard, whereas the CMS results in underestimation. Additional analyses are completed for other structures to confirm the generality of the conclusions here. These findings have potentially important implications both for the intensity‐based seismic assessments using the CS in future building codes and the risk‐based seismic assessments typically used in performance‐based earthquake engineering applications. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
The scarcity of strong ground motion records presents a challenge for making reliable performance assessments of tall buildings whose seismic design is controlled by large‐magnitude and close‐distance earthquakes. This challenge can be addressed using broadband ground‐motion simulation methods to generate records with site‐specific characteristics of large‐magnitude events. In this paper, simulated site‐specific earthquake seismograms, developed through a related project that was organized through the Southern California Earthquake Center (SCEC) Ground Motion Simulation Validation (GMSV) Technical Activity Group, are used for nonlinear response history analyses of two archetype tall buildings for sites in San Francisco, Los Angeles, and San Bernardino. The SCEC GMSV team created the seismograms using the Broadband Platform (BBP) simulations for five site‐specific earthquake scenarios. The two buildings are evaluated using nonlinear dynamic analyses under comparable record suites selected from the simulated BBP catalog and recorded motions from the NGA‐West database. The collapse risks and structural response demands (maximum story drift ratio, peak floor acceleration, and maximum story shear) under the BBP and NGA suites are compared. In general, this study finds that use of the BBP simulations resolves concerns about estimation biases in structural response analysis which are caused by ground motion scaling, unrealistic spectral shapes, and overconservative spectral variations. While there are remaining concerns that strong coherence in some kinematic fault rupture models may lead to an overestimation of velocity pulse effects in the BBP simulations, the simulations are shown to generally yield realistic pulse‐like features of near‐fault ground motion records.  相似文献   

15.
Near‐fault ground motions with forward directivity are characterized by a large pulse. This pulse‐like motion may cause a highly non‐uniform distribution of story ductility demands for code‐compliant frame structures, with maximum demands that may considerably exceed the level of code expectations. Strengthening techniques for multi‐story frame structures are explored with the objective of reducing maximum drift demands. One option is to modify the code‐based SRSS distribution of story shear strength over the height by strengthening of the lower stories of the frame. The modified distribution reduces the maximum story ductility demand, particularly for weak and flexible structures. However, this strengthening technique is less effective for stiff structures, and is almost ineffective in cases in which the maximum demand occurs in the upper stories, i.e. strong and flexible structures. As an alternative, the benefits of strengthening frames with elastic and inelastic walls are evaluated. The effects of adding walls that are either fixed or hinged at the base are investigated. It is demonstrated that strengthening with hinged walls is very effective in reducing drift demands for structures with a wide range of periods and at various performance levels. Wall inelastic behavior only slightly reduces the benefits of strengthening with hinged walls.Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

16.
The seismic response of single‐degree‐of‐freedom (SDOF) systems incorporating flag‐shaped hysteretic structural behaviour, with self‐centring capability, is investigated numerically. For a SDOF system with a given initial period and strength level, the flag‐shaped hysteretic behaviour is fully defined by a post‐yielding stiffness parameter and an energy‐dissipation parameter. A comprehensive parametric study was conducted to determine the influence of these parameters on SDOF structural response, in terms of displacement ductility, absolute acceleration and absorbed energy. This parametric study was conducted using an ensemble of 20 historical earthquake records corresponding to ordinary ground motions having a probability of exceedence of 10% in 50 years, in California. The responses of the flag‐shaped hysteretic SDOF systems are compared against the responses of similar bilinear elasto‐plastic hysteretic SDOF systems. In this study the elasto‐plastic hysteretic SDOF systems are assigned parameters representative of steel moment resisting frames (MRFs) with post‐Northridge welded beam‐to‐column connections. In turn, the flag‐shaped hysteretic SDOF systems are representative of steel MRFs with newly proposed post‐tensioned energy‐dissipating connections. Building structures with initial periods ranging from 0.1 to 2.0s and having various strength levels are considered. It is shown that a flag‐shaped hysteretic SDOF system of equal or lesser strength can always be found to match or better the response of an elasto‐plastic hysteretic SDOF system in terms of displacement ductility and without incurring any residual drift from the seismic event. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

17.
Over the past 20 years, significant socio-economic losses have been encountered in Turkey due to several moderate to large earthquakes. The studies published after the earthquakes concurringly emphasized that multistory reinforced concrete (RC) buildings, mostly 3–7 story ones, collapsed or were heavily damaged as a result of inadequate seismic performance. Global drift ratio demands are mostly used as a representative quantity for determining the behavior of structures when subjected to earthquakes. In this study, three representative mid-rise RC buildings are analyzed by nonlinear time history analysis using code-compatible real ground motion record sets and the calculated global drift ratio demands of these buildings are statistically evaluated. Ground motion record sets compatible with the design spectrum defined for local soil classes in the Turkish Earthquake Code (TEC-2007) are used for the analyses. In order to evaluate the effect of the number of ground motions on drift ratio demands, five different ground motion record sets with 7, 11 and 15 ground motion records are used separately for each local soil class. Results of this study indicate that (1) the dispersion of global drift ratio demands calculated for individual ground motion records in record sets is high, (2) local soil class has no significant effect on dispersion. However, dispersion increases in a direct proportion to the number of ground motion records in a record set, (3) the mean of global drift ratio demands calculated for different ground motion record sets may differ although they are compatible with the same design spectrum, (4) the mean of the drift demands obtained from different ground motion record sets compatible with a particular design spectrum can be accepted as simply random samples of the same population at 95% confidence level.  相似文献   

18.
The efficacy of various ground motion intensity measures (IMs) in the prediction of spatially distributed seismic demands (engineering demand parameters, (EDPs)) within a structure is investigated. This has direct implications to building‐specific seismic loss estimation, where the seismic demand on different components is dependent on the location of the component in the structure. Several common IMs are investigated in terms of their ability to predict the spatially distributed demands in a 10‐storey office building, which is measured in terms of maximum interstorey drift ratios and maximum floor accelerations. It is found that the ability of an IM to efficiently predict a specific EDP depends on the similarity between the frequency range of the ground motion that controls the IM and that of the EDP. An IMs predictability has a direct effect on the median response demands for ground motions scaled to a specified probability of exceedance from a ground motion hazard curve. All of the IMs investigated were found to be insufficient with respect to at least one of magnitude, source‐to‐site distance, or epsilon when predicting all peak interstorey drifts and peak floor accelerations in a 10‐storey reinforced concrete frame structure. Careful ground motion selection and/or seismic demand modification is therefore required to predict such a spatially distributed demands without significant bias. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

19.
The self‐centering rocking steel frame is a seismic force resisting system in which a gap is allowed to form between a concentrically braced steel frame and the foundation. Downward vertical force applied to the rocking frame by post‐tensioning acts to close the uplifting gap and thus produces a restoring force. A key feature of the system is replaceable energy‐dissipating devices that act as structural fuses by producing high initial system stiffness and then yielding to dissipate energy from the input loading and protect the remaining portions of the structure from damage. In this research, a series of large‐scale hybrid simulation tests were performed to investigate the seismic performance of the self‐centering rocking steel frame and in particular, the ability of the controlled rocking system to self‐center the entire building. The hybrid simulation experiments were conducted in conjunction with computational modules, one that simulated the destabilizing P‐Δ effect and another module that simulated the hysteretic behavior of the rest of the building including simple composite steel/concrete shear beam‐to‐column connections and partition walls. These tests complement a series of quasi‐static cyclic and dynamic shake table tests that have been conducted on this system in prior work. The hybrid simulation tests validated the expected seismic performance as the system was subjected to ground motions in excess of the maximum considered earthquake, produced virtually no residual drift after every ground motion, did not produce inelasticity in the steel frame or post‐tensioning, and concentrated the inelasticity in fuse elements that were easily replaced. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

20.
A two‐story buckling‐restrained brace (BRB) frame was tested under bidirectional in‐plane and out‐of‐plane loading to evaluate the BRB stability and gusset plate design. The test comprised pseudo‐dynamic loadings using the 1999 Chi‐Chi earthquake scaled to the 50%, 10%, and 2% probability of exceedance in 50 years and a cyclic regime of increasing amplitudes of up to 3.0% story drift ratio (SDR). The specimen had a unique configuration where the beams were connected to the columns through shear tabs welded to the column flanges and bolted to the beam webs. Stable hysteretic behavior with only minor cracking at the gusset‐to‐column welds was observed under the pseudo‐dynamic tests, with maximum in‐plane and out‐of‐plane SDRs of 2.24% and 1.47% respectively. Stable behavior continued into the cyclic test where fracture of the gusset‐to‐column welds occurred in the first cycle to simultaneous bidirectional SDR of 3.0%. The observed BRB stability is consistent with a methodology developed for BRB frames under simultaneous in‐plane and out‐of‐plane drifts. The specimen behavior was studied using a finite element model. It was shown that gusset plates are subjected to a combination of BRB force and frame action demands, with the latter increasing the gusset‐to‐beam and gusset‐to‐column interface demands by an average of 69% and 83% respectively. Consistent with the test results, failure at the gusset‐to‐column interfaces is computed when frame action demands are included, thus confirming that not considering frame action demands may results in unconservative gusset plate designs. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号