首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 500 毫秒
1.
白冰  李春峰 《岩土力学》2007,28(Z1):715-718
通过数值计算方法研究了地铁列车振动荷载作用下,由左右平行隧道过渡到上下平行隧道情形下的三维弹性动力响应。建议了三维情形下列车动力荷载的施加方式,对隧道结构变形特征以及振动加速度响应特征进行研究。计算表明,在拱底中央位置处,模型中部断面的位移峰值比起始段的位移峰值有一定减小,而其加速度峰值则变化不大;双列列车交会动载作用下的位移峰值比单行动载作用下的位移峰值要大,而相应的峰值加速度则出现得要早。  相似文献   

2.
为研究列车振动荷载作用下盾构隧道结构及周围土体的动力响应特性,采用模型试验方法,通过布置在盾构隧道底部的激振器施加扫频激振荷载和列车振动荷载,采用频率响应函数FRF与最大加速度分析了盾构隧道衬砌结构与周围土体不同位置处的动力响应及其衰减规律。研究结果表明:FRF是隧道衬砌结构和周围土体自身的动力响应特性的体现,与激振力的大小、扫频方向及扫频时间无关;在隧道管片衬砌结构的底部和顶部均体现出高频响应大于低频响应的特性,隧道顶部加速度响应沿隧道纵向衰减幅度明显小于隧道底部;隧道周围土体的动力响应沿深度有明显变化,但均表现出沿隧道轴向衰减的规律。隧道结构上部第1层测点土体的动力响应在全频域内随频率的增加逐渐增大。但在第2层和地表的第3层测点,土体的动力响应在30~90 Hz区段线性增长,在90~300 Hz区段出现波动变化,并无明显增大趋势;与扫频激振荷载引起的动力响应规律一致,由列车运行振动荷载引起的隧道管片衬砌结构和周围土体的振动沿隧道轴向逐渐衰减,隧道底部的加速度响应大于顶部,随着列车车速的增大,在隧道内部引起的加速度响应将显著增大。同时,在列车振动荷载作用下发现地表存在加速度放大效应,地表第3层测点的加速度响应均大于隧道结构上部第1层测点。  相似文献   

3.
地铁列车振动作用下近距离平行隧道的弹塑性动力响应   总被引:13,自引:1,他引:12  
白冰  李春峰 《岩土力学》2009,30(1):123-128
利用有限元数值方法,对地铁列车振动荷载作用下近距离平行隧道结构的二维弹塑性动力响应进行模拟,分析了单列列车荷载以及双列列车荷载作用下土体-隧道结构体系的加速度和内力响应时程曲线。研究表明,列车动力荷载作用下的内力值要比无列车荷载作用下的内力值有较大程度的增加;在单列列车荷载作用下,另一相邻平行隧道的加速度响应、弯矩、轴力峰值要比列车运行所在隧道相应点处的响应值要小;在双列列车荷载作用下隧道结构体系动力响应值要比单列列车荷载作用下的动力响应值有显著增加。  相似文献   

4.
以西安地铁工程为背景,设计了穿越地裂缝隧道——地层动力响应试验模型,开展了地铁振动作用下穿越地裂缝隧道——地层相互作用的动力模型试验,揭示了地裂缝活动及地铁列车振动时对位于地裂缝处地层的动力响应规律。试验结果表明:在地裂缝未活动时,隧道拱顶位置的加速度响应与拱底相比,其值相对较小,表明地铁列车振动引发拱底部的振动加速度,通过衬砌传递至上覆岩土体时,加速度发生了显著的衰减;当地裂缝上盘下降时,隧道拱底及拱顶测点产生的振动响应比地裂缝未活动时明显更为强烈。表明地裂缝的活动对隧道结构振动特性具有显著影响;受地铁列车运行位置变化和地裂缝上盘下降的双重影响,地裂缝两侧土体振动加速度幅值有明显的差异,这会对隧道结构的振动特性造成不利影响,在设计中应采取适当措施,防止造成隧道衬砌的局部损伤或破坏。  相似文献   

5.
针对地铁建设中典型的马蹄形断面隧道建立比例尺为1∶20的分析模型,并采用数值分析方法研究马蹄形隧道处于单一土层及工程所处区域典型的成层土体中时的动力响应。分析结果表明,在单一土层中由于土体的约束作用,结构产生的位移以整体沉降为主,在成层土体中除产生一定的整体变形外还伴随一定的扭转变形。在两种地层情况下马蹄形地铁隧道在地震动力作用下的动力加速度响应、竖向位移均在拱顶处产生最大值,其中在单一土层中的加速度响应最大值为结构中部加速度的2.29倍。结构在顶部和侧板处所产生的动应力响应值也较大。研究表明,地震动力荷载作用下顶板、侧板均为受力较大部位,在设计和施工中应予以充分重视。  相似文献   

6.
列车动载作用下交叠隧道动力响应数值模拟   总被引:24,自引:3,他引:24  
陈卫军  张璞 《岩土力学》2002,23(6):770-774
采用上海地铁某区间隧道现场实测振动加速度数据,对南浦大桥近距离交叠隧道在列车振动荷载作用下的动力响应进行了有限元数值模拟。模拟结果表明:上下重叠隧道在动载作用下相互影响显著;随着隧道间净距的增加,相互影响逐渐减弱;对隧道结构而言,土压静载是结构设计的控制荷载,列车振动引起的附加内力增幅较小。  相似文献   

7.
通过进行地裂缝与斜交地铁隧道的物理模型试验,研究地铁列车荷载作用下地裂缝与斜交马蹄形地铁隧道的动力相互作用特性。试验结果表明:地铁行驶产生的振动在土层中各个方向传播时会有不同程度的衰减,地裂缝对地铁振动具有阻隔作用;地裂缝附近隧道下方土层的振动要比上部土层强烈;地铁隧道的拱底部位相比拱腰和拱顶部位振动响应更强烈。地裂缝未活动时,隧道底部与土体的接触附加压力较大;地裂缝上盘下降时,位于地裂缝附近的下盘隧道底部和上盘隧道顶部与土体的接触附加压力较大。地裂缝未活动时,激振作用产生的隧道顶部和底部的附加应变均较小;地裂缝上盘下降后,位于上盘的隧道顶部和位于下盘的隧道底部产生负的附加应变,位于下盘的隧道顶部和位于上盘的隧道底部产生正的附加应变,且随上盘下降量的增大,附加应变逐渐变大。  相似文献   

8.
列车移动荷载对下穿公路隧道稳定性影响研究   总被引:1,自引:0,他引:1  
将列车荷载模拟为单一的移动荷载,考虑围岩为黏弹性成层地基,利用分层法研究了移动荷载作用下围岩的动力响应,得到了围岩在动力荷载下的受影响范围,从而可评价列车荷载对下穿公路隧道稳定性的影响程度。最后将理论分析结果与现场量测结果(主要是上方列车荷载通过时量测下方公路隧道的振动响应)进行对比分析,二者吻合较好。  相似文献   

9.
为探明地铁列车荷载作用下粉细砂地层的动力响应特征及长期累积变形发展规律,首先基于南京地铁2号线沉降监测数据,分析了位于粉细砂地层中某区段隧道结构的沉降发展规律;然后采用2.5维数值方法,模拟分析了地铁列车荷载作用下粉细砂地层的振动响应特征;最后,基于振动荷载作用下粉细砂的累积变形经验公式,采用分层总和法研究了粉细砂地层由列车荷载诱发长期累积变形的发展规律。研究结果表明:(1)列车运行引发地铁隧道地基粉细砂地层的动偏应力随着深度的增加呈先增大后减小的变化趋势,其最大值通常出现在隧道底部以下1 m处,最大值约为3.5 kPa;(2)列车振动荷载诱发粉细砂地层的长期变形在运营初期发展较快,并在运营1年左右达到稳定,随后缓慢增长,其稳定后的地基总变形量约为15.1 mm;(3)由列车振动荷载长期作用诱发的隧道粉细砂地基沉降量约占运营7年后总沉降量值的26%。  相似文献   

10.
霍元坤  吴志坚  陈拓  马巍 《冰川冻土》2011,33(4):908-913
为研究多年冻土区机车振动荷载作用下桥梁桩基的动力响应特征,选取了多年冻土区清水河特大桥的两个桥梁桩基进行现场强震动测试.对测试数据进行加速度响应及傅里叶谱响应分析,得出桥梁桩基的振动响应衰减规律及振动能量的分布规律.在此基础上,结合有限元结构分析软件,运用瞬态动力分析法,选取一个测试桩为实例进行二维动力响应的数值计算,...  相似文献   

11.
In order to avoid the pollution of trace metals in marine environment, it is necessary to establish the data and understand the mechanisms influencing the distribution of trace metals in marine environment. The concentration of heavy metals (Fe, Mn, Cr, Cu, Ni, Pb, Zn, Co and Cd) were studied in sediments of Ennore shelf, to understand the metal contamination due to heavily industrialized area of Ennore, south-east coast of India. Concentration of metals shows significant variability and range from 1.7 to 3.7% for Fe, 284–460 μg g−1 for Mn, 148.6–243.2 μg g−1 for Cr, 385–657 μg g−1 for Cu, 19.8–53.4 μg g−1 for Ni, 5.8–11.8 μg g−1 for Co, 24.9–40 μg g−1 for Pb, 71.3–201 μg g−1 for Zn and 4.6–7.5 μg g−1 for Cd. For various metals the contamination factor (CF) and geoaccumulation index (I geo) has been calculated to assess the degree of pollution in sediments. The geoaccumulation index shows that Cd, Cr and Cu moderately to extremely pollute the sediments. This study shows that the major sources of metal contamination in the Ennore shelf are land-based anthropogenic ones, such as discharge of industrial wastewater, municipal sewage and run-off through the Ennore estuary. The intermetallic relationship revealed the identical behavior of metals during its transport in the marine environment.  相似文献   

12.
Isothermal thermogravimetric experiments have been carried out to determine the reaction kinetics of the dehydration processes in fuller's earth, a natural Ca-montmorillonite. Dehydration in swelling clays is a complex reaction, and analysis of the thermogravimetric data using empirical rate equations and time-transformation analysis reveals that the nature of the rate controlling mechanism is dependent upon both the temperature regime of the sample as well as the extent of reaction. For fuller's earth, we find that the dehydration kinetics are dominated by a nucleation and growth mechanism at low temperatures and fractions transformed (stage I), but above 90 °C the last stages of the reaction are diffusion controlled (stage II). The activation energy for dehydration during stage I is around 35 kJ · mol−1, whereas the removal of water during stage II requires an activation energy of around 50 kJ · mol−1. These two stages of dehydration are associated with primary collapse of the interlayer (stage I) and movement of water that is hydrated to cations within the interlayer (stage II). Received: 28 August 1998 / Revised, accepted: 27 January 1999  相似文献   

13.
A pressure-volume-temperature data set has been obtained for lawsonite [CaAl2Si2O7(OH)2.H2O], using synchrotron X-ray diffraction and an externally heated diamond anvil cell. Unit-cell volumes were measured to 9.4 GPa and 767 K by angle dispersive X-ray diffraction using imaging plates. Phase changes were not observed within this pressure-temperature range, and lawsonite compressed almost isotropically at constant temperature. The P-V-T data have been analyzed using a Birch- Murnaghan equation of state and a linear equation of state expressed as β=–1/V0 (∂V/∂P) T . At room temperature, the derived equation of state parameters are: K 0=124.1 (18) GPa K'0 set to 4) and β–1=142.0(24) GPa, respectively. Our results are intermediate between previously reported measurements. The high-temperature data show that the incompressibility of lawsonite decreases with increasing temperature to ∼500 K and then increases above. Hence, the second order temperature derivative of the bulk modulus is taken into account in the equation of state; a fit of the volume data yields K 0=123.9(18) GPa, (∂K/∂T)P=–0.111(3) GPa K–1, (∂2 K/∂T 2)P=0.28(6) 10–3 GPa K–2, α0=3.1(2) 10–5 K–1, assuming K'0=4. Received: 2 June 1998 / Revised, accepted: 12 Ocotber 1998  相似文献   

14.
Although, the Eastern Desert of Egypt forms about 22% of the surface area of the country, the area is undeveloped due to the limited availability of water. The morphologic units of the Eastern Desert consist of a number of drainage basins covering about 147,820 km2 (66.5% of the total surface area of the Eastern Desert). The basins drain the occasional rainwater, either towards the Nile Valley or to the Red Sea, causing flood hazards. The availability of water from the hydrologic systems of these basins could be improved by constructing runoff controlling systems in these areas (e.g. dykes and partially effective dams), which could save and make use of a considerable amount of water. The groundwater resources in the Eastern Desert can be divided into four main water-bearing units: the fractured crystalline Pre-Cambrian aquifer, the Nubian sandstone aquifer, the fractured limestone and sandstone aquifer and the Quaternary aquifer. The most productive aquifer is the Nubian sandstone while the fractured limestone and sandstone (Miocene) are only productive along the eastern part of the desert. The Quaternary aquifer occurs along the major dry washes (wadis) and is considered of limited potential as it is recharged mainly from the occasional rainfall. Detailed assessment of these aquifers should be carried out locally for further development of the area.
Resumen Aunque el Desierto Oriental de Egipto constituye alrededor del 22% de la superficie del país, esta área no se encuentra desarrollada por causa de la escasez de agua. Las unidades morfológicas del Desierto Oriental consisten de un número de cuencas de drenaje que cubren alrededor de 147,820 km2 (66.5% del área superficial total del Desierto Oriental). Estas cuencas drenan las aguas lluvias ocasionales, bien hacia el Valle del Nilo o hacia el Mar Rojo, causando amenazas de inundación. La disponibilidad de agua a partir de los sistemas hidrológicos de estas cuencas, puede ser mejorada al construir sistemas de control de escorrentía allí (Ej. Diques y presas de efectividad parcial), las cuales podrían almacenar y permitir el uso de una cantidad considerable de agua. Los recursos de agua subterránea en el Desierto Oriental, pueden dividirse en cuatro unidades principales portadoras de agua: El acuífero fracturado cristalino Pre – Cámbrico, el acuífero de la Arenisca de Nubia, el acuífero de arenisca y caliza fracturadas y el acuífero Cuaternario. El acuífero más productivo es la arenisca de Nubia, mientras que la arenisca y caliza fracturadas (Mioceno), son productivas únicamente a lo largo de la parte oriental del desierto. El acuífero Cuaternario se encuentra a lo largo de las corrientes intermitentes mayores (Wadis) y se le considera de potencial limitado por ser recargado principalmente a partir de lluvia ocasional. Una evaluación detallada de estos acuíferos debe llevarse a cabo localmente, para el desarrollo adicional de esta área.

Résumé Bien que le Désert Est occupe 22% de la surface de lEgypte, la région est sousdeveloppée à cause de la ressource en eau limitée. Les unités morphologiques consistent en des bassins drainant qui couvrent une surface de 147,820 km2 ,représentant 66.6% de la surface du Désert Est. Les bassins drainent les plues intermittentes vers la vallée du Nil ou vers la Mer Rouge en provocant des inondations. On peut augmenter la ressource en eau dans cette région en réalisant des systèmes de contrôle de ruissellements (des digues...) qui peuvent sauver et utiliser un volume considérable deau. La ressource en eau souterraine du Désert Est est cantonnée dans quatre grandes structures: l› aquifère fracturé cristallin Précambrien, l› aquifère gréseux Nubien, l› aquifère calcaire et gréseux et l› aquifère quaternaire. Le plus productif est l› aquifère gréseux Nubien tendis que les calcaire et les grèses dage miocène sont productives seulement dans la partie est du désert. L› aquifère quaternaire se trouve au long des oueds sèches et on considéré quil a un potentiel limité, étant rechargé seulement par les plues intermittentes. Pour le futur développement de la région il est nécessaire une évaluation détaillée, à léchelle locale de ces aquifères.
  相似文献   

15.
The lithostratigraphy, depositional environment and age of the Marwar Supergroup have been reviewed in the light of report of δ13C depletion recorded in the carbonates of the Bilara Group (middle part of Marwar Supergroup) and discovery of trilobite-like trace fossils from the ·Red bedsŽ of Nagaur Group (upper part of Marwar Supergroup). The δ13C depletion observed in Bilara carbonates is not a result of glaciation rather due to rapid burial and poor water circulation in the low energy water of the protected basin. Secondly, the trace fossils are, in fact, traces of notostracan crustaceans found in shallow fluvial and shallow lacustrine environment. The present paper also records a spiral, burrowing trace-fossil, possibly Gyrolithes, from a cross-bedded sandstone of the Jodhpur Group.  相似文献   

16.
In situ X-ray diffraction measurements of Fe- and Al-bearing MgSiO3-rich perovskite (FeAl-Pv), which was synthesized from a natural orthopyroxene, were performed at pressures of 19–32 GPa and temperatures of 300–1,500 K using a combination of a Kawai-type apparatus with eight sintered-diamond anvils and synchrotron radiation. Two runs were performed using a high-pressure cell with two sample chambers, and both MgSiO3 perovskite (Mg-Pv) and FeAl-Pv were synthesized simultaneously in the same cell. Thus we were able to measure specific volumes (V/V 0) of Mg-Pv and FeAl-Pv at the same P−T conditions. At all the measurement conditions, values of the specific volume of FeAl-Pv are consistent with those of Mg-Pv within 2 Standard Deviation, strongly suggesting that effect of incorporation of iron and aluminum on the thermoelastic properties of magnesium silicate perovskite is undetectable in this composition, pressure, and temperature range. Two additional runs were performed using a high-pressure cell that has one sample chamber and unit-cell volumes of FeAl-Pv were measured at pressures and temperatures up to 32 GPa and 1,500 K, respectively. All the unit-cell volume data of FeAl-Pv perovskite were fitted to the high temperature Birch–Murnaghan equation of state and a complete set of thermoelastic parameters of this perovskite was determined with an assumption of K′ 300,0 = 4. The determined parameters are K 300,0 = 243(3) GPa, (∂K T,0/∂T) P = −0.030(8) GPa/K, a 0 = 2.78(18) × 10−5 K−1, and b 0 = 0.88(28) × 10−8 K−2, where a 0 and b 0 are the coefficients of the following expression describing the zero-pressure thermal expansion: α T,0 = a 0 + b 0 T. The equation-of-state parameters of FeAl-Pv are in good agreement with those of MgSiO3 perovskite at the conditions corresponding to the uppermost part of the lower mantle.  相似文献   

17.
On the south-eastern edge of Russia, the chemical composition of rainwater is controlled by sea salts, terrestrial material, as well as volcanic (Kuril islands volcanic area) and anthropogenic emissions, mostly in the southern part of the area. The predominant major ions of the Primorye, Sakhalin and the Kuril Islands rainwaters were respectively HCO3–SO42−, Ca–Na, and of Cl–Na. Concentration of trace elements changes within 1–2 orders of magnitude but some difference in the distribution of the elements between continental and island rainwater is found. The concentration of the chemical elements in the particulate fraction varies from < 10% to 90% of the total concentration (dissolved + particulate) with the following distribution: Tl, Na, Ca, Sr, Zn, Cd (< 10%)–Be, Th, Bi, Rb, U, K, Sc (10–20%)–Cu, Mn, Mg, Mo, Se, Ba, Ni, As, Ag, Cs, Co, Y, Ga, V (20–50%)–Sb, Pb, Ge, Cr, Fe, Al (50–90%).The concentration of elements of the particulate fraction of the rainwater usually is significantly different from concentrations in the crust, including both higher and lower concentrations. The terrestrial contribution to dissolved elements was evaluated and follows the decreasing order: Fe > K, Mg, Ca > Ba, Sr > Na (65–1%). Close order was found for total (dissolved and solid) concentrations. Sea salt contribution to dissolved element concentration in the rainwater decrease in the following order: Cl, Mg > K, SO4 > Ca > HCO3, Ba, Fe (78–0.1%). Calculation of anthropogenic and volcanic inputs for two ions (Cl and SO42−) shows that anthropogenic inputs for the Vladivostok and Yuzno-Sakhalinsk cities can be evaluated as 15–20% of Cl and up to 80–90% of SO42−. Volcanic components in the Kuril Islands, where anthropogenic inputs are absent, can reach up to 76% of SO42− and 36% of Cl.  相似文献   

18.
 Ca-chloride waters are defined as those in which Q=rCa/r(SO4+HCO3)>1, rNa/rCl<0.80, rMg/rCa<0.5 and wCl/wBr<286 (r=meq l–1 and w=mg l–1). Throughout the last 50 years, different models for the formation of such waters have been suggested. These models include: (1) filtration through semipermeable membranes under conditions of highly compacted argillaceous sediments, (2) deaquation of seawater by evaporation and/or by freezing followed by dolomitization, (3) hydrolysis of plagioclase and biotites in igneous metamorphic rock masses, (4) radiolytic modification of residual metamorphic fluids, and (5) dissolution of chalks followed by ion exchange on smectites. The better understanding of processes and of natural environments leading to the evolution and natural occurrence of such brines, is imperative for the prospection and further sustained exploitation of such waters. Received: 11 October 1996 · Accepted: 24 February 1997  相似文献   

19.
This study presents a numerical investigation of the effect of volatilization and the supercritical state of pore water on maturation of organic matter in host rocks based on the heat flow models assuming the instantaneous and finite-time intrusion mechanisms of magma. A 15 m thick, well-dated basic sill in the DSDP 41-368 hole near Cape Verde Rise, eastern Atlantic is selected as an example due to the sufficient thermophysical parameters of rocks and the definite burial and thermal history of the shale host rocks. Results indicate: (1) The effect of the temperature-dependent thermal properties of pore water at a hydrostatic pressure of 414 bar on the predicted vitrinite reflectance (Rr) is less than 0.1% no matter which intrusion mechanism of magma is assumed and can hence be ignored reasonably; (2) The consideration of volatilization of pore water can reduce the predicted Rr of host rocks significantly. In case of the instantaneous intrusion mechanism, the maximum deviation of the predicted Rr caused by pore-water volatilization reaches 1.3% at the location of half the sill thickness away from the contact (i.e. X/D = 0.5), and the deviation above 0.5% can occur in the region from 0.3 to 1.0 in the form of X/D. In case of the finite-time intrusion mechanism, the maximum deviation of the predicted Rr due to pore-water volatilization attains 1.15% at X/D = 0.25, and the region where the deviation is larger than 0.5% lies between 0.15 and 0.6 in the form of X/D; (3) If hydrothermal convection in the host rocks is allowed for, the predicted Rr of the overlying host rocks is less than that of the underlying host rocks at the same X/D in the inner region of the contact aureole of igneous intrusions, whereas the phenomenon is converse in the outer region. In contrast, the measured Rr profile shows that at the same X/D, Rr of the overlying host rocks is totally higher than that of the underlying host rocks. Thus, it is not the hydrothermal convection in the overlying host rocks that resulted in the asymmetry of the current Rr profiles below and above the sill; (4) The predicted Rr based on the heat conduction model assuming the finite-time intrusion mechanism and pore-water volatilization matches well with the measured one out of the region where the Rr geothermometer is unreliable due to the effect of volatilization of pore water. This demonstrates that the finite-time intrusion mechanism of magma, together with pore-water volatilization, possibly represents natural conditions.  相似文献   

20.
激光焦平面变化对LA-ICPMS锆石U-Pb定年准确度的影响   总被引:3,自引:0,他引:3  
元素分馏是影响LA-ICPMS锆石U-Pb定年准确度的重要因素之一,通常利用标样进行校正。激光聚焦位置变化会引起剥蚀坑形貌及U-Pb分馏的改变,标样和样品聚焦条件不一致将导致标样难以准确校正样品,并最终影响定年结果的准确性,但影响的程度、机制及可容忍范围目前尚不清楚。为此,文章以91500为标样、GJ-1为样品,详细研究了聚焦偏离30μm范围内剥蚀标样与样品锆石的剥蚀坑形貌变化以及由此导致的U-Pb定年误差。实验表明,在距离锆石表面30μm范围内,标样和样品焦平面同步变化时,二者U-Pb分馏形式及程度基本一致,激光焦平面偏离所引起的样品年龄与TIMS推荐值的偏差小于1%;当二者聚焦不同步时,标样与样品的U-Pb分馏差别显著,年龄偏差最大可超过3%。激光聚焦不同步导致的标样与样品剥蚀坑纵横比差异是引起年龄误差的根本原因,激光焦平面偏离锆石表面超过15μm,剥蚀坑坑口明显变大,纵横比减小,U-Pb分馏形式及程度发生改变。通过预剥蚀锆石,观察剥蚀坑轮廓,使激光焦平面在距离锆石样品表面15μm范围内,可确保标样与样品剥蚀坑形貌及U-Pb分馏状态一致,提高LA-ICPMS定年的准确度。  相似文献   

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

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