Modified precise time step integration method of structural dynamic analysis   总被引：2，自引：1，他引：1  

汪梦甫  周锡元 《地震工程与工程振动(英文版)》2005,4(2):287-293

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.  相似文献   

Probabilistic seismic response and reliability assessment of isolated bridges     

GiuseppeCarloMarano 《地震工程与工程振动(英文版)》2005,4(1):95-106

Bridge seismic isolation strategy is based on the reduction of shear forces transmitted from the superstructure to the piers by two means: shifting natural period and earthquake input energy reduction by dissipation concentrated in protection devices. In this paper, a stochastic analysis of a simple isolated bridge model for different bridge and device parameters is conducted to assess the efficiency of this seismic protection strategy. To achieve this aim, a simple nonlinear softening constitutive law is adopted to model a wide range of isolation devices, characterized by only three essential mechanical parameters. As a consequence of the random nature of seismic motion, a probabilistic analysis is carried out and the time modulated Kanai-Tajimi stochastic process is adopted to represent the seismic action. The response covariance in the state space is obtained by solving the Lyapunov equation for a stochastic linearized system. After a sensitivity analysis, the failure probability referred to extreme displacement and the mean value of dissipated energy are assessed by using the introduced stochastic indices of seismic bridge protection efficiency. A parametric analysis for protective devices with different mechanical parameters is developed for a proper selection of parameters of isolation devices under different situations.  相似文献   

Centrifuge modeling of PGD response of buried pipe   总被引：1，自引：1，他引：0  

MichaelO＇Rourke VikramGadicherla TarekAbdoun 《地震工程与工程振动(英文版)》2005,4(1):69-73

A new centrifuge based method for determining the response of continuous buried pipe to PGD is presented. The physical characteristics of the RPI‘s 100 g-ton geotechnical centrifuge and the current lifeline experiment split-box are described: The split-box contains the model pipeline and surrounding soil and is manufactured such that half can be offset, in flight, simulating PGD. In addition, governing similitude relations which allow one to determine the physical characteristics, (diameter, wall thickness and material modulus of elasticity) of the model pipeline are presented. Finally, recorded strains induced in two buried pipes with prototype diameters of 0.63 m and 0.95 m (24 and 36 inch) subject to 0.6 and 2.0 meters (2 and 6 feet) of full scale fault offsets and presented and compared to corresponding FE results.  相似文献   

Measurement of six degree-of-freedom ground motion by using eight accelerometers     

杨振宇沈毅  刘志言 《地震工程与工程振动(英文版)》2005,4(2):229-232

A new integrated measuring system with eight force-balance accelerometers is proposed to obtain a direct measurement of six degree-of freedom （DOF） ground motions, including three rotational and three actual translational acceleration components without gyroscopes. In the proposed measuring system, the relationship between the output from eight force-balance accelerometer and the six DOF motion of the measuring system under an earthquake are described by differential equations. These equations are derived from the positions and directions of the eight force-balance accelerometers in the measuring system. The third-order Runge-Kutta algorithm is used to guarantee the accuracy of the numerical calculation. All the algorithms used to compute the six DOF components of the ground motion are implemented in a real-time in Digital Signal Processor （DSP）. The distortion of the measured results caused by position and direction errors of the accelerometers in the measuring system are reduced by multiplying a compensation coefficient C to the output and subtracting static zero drift from the measured results, respectively.  相似文献   

Geochemical characteristics of crustal anatexis during the formation of migmatite at the Southern Sierra Nevada,California   总被引：14，自引：0，他引：14  

Lingsen?ZengEmail author  Jason?B.?Saleeby  Mihai?Ducea 《Contributions to Mineralogy and Petrology》2005,150(4):386-402

We provide data on the geochemical and isotopic consequences of nonmodal partial melting of a thick Jurassic pelite unit at mid-crustal levels that produced a migmatite complex in conjunction with the intrusion of part of the southern Sierra Nevada batholith at ca. 100 Ma. Field relations suggest that this pelitic migmatite formed and then abruptly solidified prior to substantial mobilization and escape of its melt products. Hence, this area yields insights into potential mid-crustal level contributions of crustal components into Cordilleran-type batholiths. Major and trace-element analyses in addition to field and petrographic data demonstrate that leucosomes are products of partial melting of the pelitic protolith host. Compared with the metapelites, leucosomes have higher Sr and lower Sm concentrations and lower Rb/Sr ratios. The initial ⁸⁷Sr/⁸⁶Sr ratios of leucosomes range from 0.7124 to 0.7247, similar to those of the metapelite protoliths (0.7125–0.7221). However, the leucosomes have a much wider range of initial ε_Nd values, which range from −6.0 to −11.0, as compared to −8.7 to −11.3 for the metapelites. Sr and Nd isotopic compositions of the leucosomes, migmatites, and metapelites suggest disequilibrium partial melting of the metapelite protolith. Based on their Sr, Nd, and other trace-element characteristics, two groups of leucosomes have been identified. Group A leucosomes have relatively high Rb, Pb, Ba, and K₂O contents, Rb/Sr ratios (0.15<Rb/Sr<1.0), and initial ε_Nd values. Group B leucosomes have relatively low Rb, Pb, Ba, and K₂O contents, Rb/Sr ratios (<0.15), and initial ε_Nd values. The low Rb concentrations and Rb/Sr ratios of the group B leucosomes together suggest that partial melting was dominated by water-saturated or H₂O-fluxed melting of quartz + feldspar assemblage with minor involvement of muscovite. Breakdown of quartz and plagioclase with minor contributions from muscovite resulted in low Rb/Sr ratios characterizing both group A and group B leucosomes. In contrast, group A leucosomes have greater contributions from K-feldspar, which is suggested by: (1) their relatively high K concentrations, (2) positive or slightly negative Eu anomalies, and (3) correlation of their Pb and Ba concentrations with K₂O contents. It is also shown that accessory minerals have played a critical role in regulating the partitioning of key trace elements such as Sm, Nd, Nb, and V between melt products and residues during migmatization. The various degrees of parent/daughter fractionations in the Rb–Sr and Sm–Nd isotopic systems as a consequence of nonmodal crustal anatexis would render melt products with distinct isotopic signatures, which could profoundly influence the products of subsequent mixing events. This is not only important for geochemical patterns of intracrustal differentiation, but also a potentially important process in generating crustal-scale as well as individual pluton-scale isotopic heterogeneities.  相似文献   

Chemical composition of rock-forming minerals in gold-related granitoid intrusions,southwestern New Brunswick,Canada: implications for crystallization conditions,volatile exsolution,and fluorine-chlorine activity     

Xue-Ming?YangEmail author  David?R.?Lentz 《Contributions to Mineralogy and Petrology》2005,150(3):287-305

Chemical composition of rock-forming minerals in Appalachian Siluro-Devonian granitoid intrusions, southwestern New Brunswick, was systematically determined by electron microprobe. The mineral chemical data together with petrographic examination was used to test magmatic equilibration and to constrain crystallization conditions, volatile exsolution, and fluorine-chlorine activity of fluids associated with these intrusions. Mineralogical distinction between Late Silurian to Early Devonian granodioritic to monzogranitic series (GMS) and Late Devonian granitic series (GS) rocks is evident, although both are subsolvus I-type to evolved I-type granitoids. Oxidized to reduced GMS rocks consist of quartz, plagioclase (An>10), K-feldspar, biotite, apatite, titanite, zircon, monazite, ± hornblende, ± pyroxene, ± magnetite, ± ilmenite, and ± sulfide. GS rocks comprise quartz, K-feldspar, plagioclase (An<10), mica group minerals, zircon, monazite, apatite, sulfide, ± ilmente, ± magnetite, ± topaz, ± columbite, and ± xenotime. Inter-intrusion and intra-intrusion variations in mineral chemistry are interpreted to reflect petrogenetic processes (e.g., assimilation and fractional crystallization) during granitoid evolution. Although magmatic equilibration among rock-forming minerals are disturbed by subsolidus hydrothermal processes, GMS rocks appear to have higher magmatic temperatures, variable levels of emplacement, a range of (i.e., reduced intrusions 10^−16.7∼10^−13.4 and oxidized intrusions 10^−14.0∼10^−10.5 bars), and relatively low f _HF/f _HCl ratios (10^−3.0∼10^−1.0) in exsolved fluids, compared to GS rocks. Reduced GMS intrusions bear higher gold potential and thus may be prospective targets for intrusion-related gold systems. Electronic Supplementary Material Supplementary material is available for this article at  相似文献   

Refining the timing of eclogite metamorphism: a geochemical,petrological, Sm-Nd and U-Pb case study from the Pohorje Mountains,Slovenia (Eastern Alps)   总被引：4，自引：2，他引：2  

Christine?MillerEmail author  Roland?Mundil  Martin?Th?ni  Jürgen?Konzett 《Contributions to Mineralogy and Petrology》2005,150(1):70-84

High-pressure metamorphism in the Pohorje Mountains of Slovenia (Austroalpine unit, Eastern Alps) affected N-MORB type metabasic and metapelitic lithologies. Thermodynamic calculations and equilibrium phase diagrams of kyanite–phengite-bearing eclogites reveal PT conditions of >2.1 GPa at T<750°C, but within the stability field of quartz. Metapelitic eclogite country rocks contain the assemblage garnet + phengite + kyanite + quartz, for which calculated peak pressure conditions are in good agreement with results obtained from eclogite samples. The eclogites contain a single population of spherical zircon with a low Th/U ratio. Combined constraints on the age of metamorphism come from U/Pb zircon as well as garnet–whole rock and mineral–mineral Sm-Nd analyses from eclogites. A coherent cluster of single zircon analyses yields a ²⁰⁶Pb/²³⁸U age of 90.7±1.0 Ma that is in good agreement with results from Sm-Nd garnet–whole rock regression of 90.7±3.9 and 90.1±2.0 Ma (εNd: +8) for two eclogite samples. The agreement between U-Pb and Sm-Nd age data strongly suggests an age of approximately 90 Ma for the pressure peak of the eclogites in the Pohorje Mountains. The presence of garnet, omphacite and quartz inclusions in unfractured zircon indicates high-pressure rather than ultrahigh pressure conditions. The analysed metapelite sample yields a Sm-Nd garnet–whole rock scatterchron age of 97±15 Ma. These data probably support a single P-T loop for mafic and pelitic lithologies of the Pohorje area and a late Cretaceous high-pressure event that affected the entire easternmost Austroalpine basement including the Koralpe and Saualpe eclogite type locality in the course of the complex collision of the Apulian microplate and Europe.  相似文献   

A Mössbauer study of pseudotachylytes: redox conditions during seismogenic faulting     

Kieran?D.?O’HaraEmail author  Frank?E.?Huggins 《Contributions to Mineralogy and Petrology》2005,148(5):602-614

The redox conditions during frictional melting provide information on the physical and chemical conditions during seismic slip in the crust. Here we examine frictional melts from five localities by analyzing host rocks and corresponding pseudotachylytes using Mössbauer spectroscopy. The faults examined are located at South Mountain, Arizona; Fort Foster, Maine (two localities); Long Ridge fault, North Carolina; and the Homestake shear zone, Colorado. The main iron-bearing phases in the pseudotachylytes are phyllosilicates (biotite, muscovite and clays) and iron oxides (magnetite and hematite) and minor pyrite. The ferrous/ferric ratios of the phyllosilicates in the host rocks are the same as those in the pseudotachylytes, except for the hematite-bearing pseudotachylyte from the Long Ridge fault, which is more oxidized. The magnetites in the host rocks and the corresponding pseudotachylytes have different ferric and ferrous iron distributions, which is attributed to different cation chemistry, rather than redox conditions. With the exception of the South Mountain locality, the ferric/ferrous ratios of the micas are interpreted to record the primary redox state of the pseudotachylyte melt as the calculated oxygen fugacities are consistent with magnetite and hematite equilibria. Pyrite-bearing pseudotachylytes plot ~0–1 log₁₀ units above the fayalite-magnetite-quartz (FMQ) buffer. Magnetite-bearing pseudotachylytes plot ~2–4 log₁₀ units above the FMQ buffer, and hematite-bearing pseudotachylytes plot 3.5 log₁₀ units above the hematite-magnetite (HM) buffer. Hematite-bearing pseudotachylytes, together with previous oxygen isotope data, are inferred to represent melting in the presence of externally derived pressurized water. Other localities are inferred to represent melting under rock-buffered, closed system, conditions. If the localities studied are representative of seismogenic faulting, the calculated oxygen fugacities indicate that, in the system C–O–H–S, H₂O and CO₂ should be the dominant fluid species. This is the first detailed study of the redox state of pseudotachylytes.  相似文献   

Contrasting interactions of sodium and potassium with H<Subscript>2</Subscript>O in haplogranitic liquids and glasses at 200 MPa from hydration–diffusion experiments     

Antonio?Acosta-Vigil  David?LondonEmail author  George?B.?Morgan VI 《Contributions to Mineralogy and Petrology》2005,149(3):276-287

This study examines hydration–diffusion in the metaluminous haplogranite system at 200 MPa H₂O and 800–300°C. At 800°C hydration is accompanied by melting and uphill diffusion of sodium from anhydrous glass toward the region of hydration and melting, whereas potassium diffuses away from the hydration front and into anhydrous glass. Silicon and aluminum are simply diluted upon hydration. There is no change in molecular Al/(Na + K) throughout the entire hydration-diffusion aureole and, therefore, (1) there is no loss of alkalis to the vapor, and (2) K migrates to replace Na in order to maintain local charge balance required by ^IVAl. Alkali diffusion occurs over a viscosity contrast from 10^4.1 Pa s in hydrous liquid to 10^11.8–10^13.5 Pa s in anhydrous glass. From these results, we interpret that: (1) Na is structurally or energetically favored over K as a charge-balancing cation for ^IVAl in hydrous granitic liquids, whereas the opposite behavior has been observed for anhydrous melts, and (2) the diffusion of alkalis through silicate melts is largely independent of viscosity. Results from 600°C are similar to those at 800°C, but hydration at 300°C involves a loss of Na and concomitant increase in molar Al/(Na + K) in the hydration zone due to hydrogen-alkali exchange between fluid and glass. Hydration behavior at 400°C is transitional between those at 300°C and 600°C, suggesting that the change in hydration mechanism occurs near the glass transition.  相似文献   

Glassy orthopyroxene granodiorites of the Pannonian Basin: tracers of ultra-high-temperature deep-crustal anatexis triggered by Tertiary basaltic volcanism     

Monika?HuraiováEmail author  Jean?Dubessy  Patrik?Kone?ny  Klaus?Simon  Ján?Král’  Grzegorz?Zielinski  Jozef?Lipka  Vratislav?Hurai 《Contributions to Mineralogy and Petrology》2005,148(5):615-633

Glassy orthopyroxene granodiorite-tonalite (named pincinite after type locality) was described from basaltic lapilli tuffs of the Pliocene maar near Pinciná village in the Slovakian part of the Pannonian Basin. Two pincinite types exhibit a qualitatively similar mineral composition (quartz, An_20–55 plagioclase, intergranular silicic glass with orthopyroxene and ilmenite, ±K-feldspar), but strongly different redox potential and formation PT conditions. Peraluminous pincinite is reduced (6–7% of total iron as Fe³⁺ in corundum-normative intergranular dacitic glass) and contains ilmenite with 8–10 mol% Fe₂O₃ and orthopyroxene dominated by ferrosilite. High-density (up to 0.85 g/cm³) primary CO₂ inclusions with minor H₂, CH₄, H₂S, CO and N₂ (<2 mol% total) are present in Qtz and Plg. Equilibrium PT conditions inferred from the intergranular Opx–Ilm–Glass assemblage and fluid density correspond to 1,170±50°C, 5.6±0.4 kbar, respectively. Metaluminous pincinite is more oxidised (25–27% of total iron as Fe³⁺ in diopside-normative intergranular glass of rhyolite–trachyte–dacite composition) and contains Fe₂O₃-rich ilmenite (17–29 mol%) associated with enstatite. Fluid inclusions are composed of CO₂–H₂O mixtures with up to 38 mol% H₂O. Raman spectroscopy revealed H₂S along with dominant CO₂ in the carbonic phase. Equilibrium PT parameters for the intergranular Opx–Ilm–Glass assemblage correspond to 740±15°C, 2.8±0.1 kbar, respectively. Reducing gas species (<2 mol% total) in the CO₂-inclusions of the peraluminous pincinite resulted from hydrogen diffusion due to fH₂ gradient imposed during decrease of redox potential from the log fO₂ values near QFM during Qtz + Plg growth, to QFM-2 incidental to the superimposed Opx + Ilm assemblage in the intergranular melt. The decrease in oxygen fugacity was recorded also in the metaluminous pincinite, where log fO₂ values changed from ~QFM + 2.6 to QFM + 0.4, but hydrogen diffusion did not occur. Absence of OH-bearing minerals, major and trace element abundances (e.g. REE 300–320, Nb 55–57, Th 4–31, Zr 240–300 ppm, FeO_tot/MgO up to 11), and Sr–O isotope ratios in the pincinites are diagnostic of high-temperature anorogenic magmas originated by dehydration melting of biotite in quartz-feldspathoid crust (⁸⁷Sr/⁸⁶Sr>0.705–0.706, ¹⁸O>9 V-SMOW) around alkali basalt reservoir in depths between 17 and 20 km, and around late stage derivatives of the basalt fractionation, intruding the crust up to depths of 10–11 km. Low water activity in the pincinite parental melt was caused by CO₂-flux from the Tertiary basaltic reservoirs and intrusions. The anatexis leads to generation of a melt-depleted granulitic crust beneath the Pannonian Basin, and the pincinites are interpreted as equivalents of igneous charnockites and enderbites quenched at temperatures above solidus and unaffected by sub-solidus re-equilibration and metamorphic overprint.  相似文献