共查询到20条相似文献，搜索用时 15 毫秒
1.
《Soil Dynamics and Earthquake Engineering》2003,23(4):249-253
A parallel soil–structure interaction (SSI) model is presented for applications on distributed computer systems. Substructring method is applied to the SSI system and a coupled finite–infinite element based parallel computer program is developed. In the SSI system, infinite elements are used to represent the soil which extends to infinity. In this case, a large finite element mesh is required to define the near field for reliable predictions. The resulting large-scale problems are solved on distributed computer systems in this study. The domain is represented by separated substructures and an interface. The number of substructures are determined by the available processors in the parallel platform. To avoid the formation of large interface equations, smaller interface equations are distributed to processors while substructure contributions are performed. This saves a lot of memory storage and computational effort. Direct solution techniques are used for the solution of interface and substructure equation systems. The program is investigated through some example problems. The example problems exposed the need for solving large-scale problems in order to reach better results. The results of the example problems demonstrated the benefits of the parallel SSI algorithm. 相似文献
2.
In this paper, we show that supercritical fluids have a greater significance in the generation of pegmatites,and for ore-forming processes related to granites than is usually assumed. We show that the supercritical melt or fluid is a silicate phase in which volatiles; principally H_2O are completely miscible in all proportions at magmatic temperatures and pressures. This phase evolves from felsic melts and changes into hydrothermal fluids, and its unique properties are particularly important in sequestering and concentrating low abundance elements, such as metals. In our past research, we have focused on processes observed at upper crustal levels, however extensive work by us and other researchers have demonstrated that supercritical melt/fluids should be abundant in melting zones at deep-crustal levels too. We propose that these fluids may provide a connecting link between lower and upper crustal magmas,and a highly efficient transport mechanism for usually melt incompatible elements. In this paper, we explore the unique features of this fluid which allow the partitioning of variouselements and compounds, potentially up to extreme levels,and may explain various features both of mineralization and the magmas that produced them. 相似文献
3.
《Advances in water resources》1999,23(3):229-237
It can be very time consuming to use the conventional numerical methods, such as the finite element method, to solve convection–dispersion equations, especially for solutions of large-scale, long-term solute transport in porous media. In addition, the conventional methods are subject to artificial diffusion and oscillation when used to solve convection-dominant solute transport problems. In this paper, a hybrid method of Laplace transform and finite element method is developed to solve one- and two-dimensional convection–dispersion equations. The method is semi-analytical in time through Laplace transform. Then the transformed partial differential equations are solved numerically in the Laplace domain using the finite element method. Finally the nodal concentration values are obtained through a numerical inversion of the finite element solution, using a highly accurate inversion algorithm. The proposed method eliminates time steps in the computation and allows using relatively large grid sizes, which increases computation efficiency dramatically. Numerical results of several examples show that the hybrid method is of high efficiency and accuracy, and capable of eliminating numerical diffusion and oscillation effectively. 相似文献
4.
Partitioning behavior between amphibole and silicate glass of thirty-three minor and trace elements(Sc,Ti, V, Cr, Co, Rb, Sr, P, Y, Zr, Nb, Cs, Ba, K, La, Ce, Pr,Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta, Pb,Th, and U) have been determined experimentally. Products of crystallization of hydrous basalt melts from 0.6 GPa/860 °C up to 2.6 GPa/970 °C were obtained in a multianvil apparatus. Major and trace element compositions of amphibole and glass were determined with a combination of electron microprobe and laser ablation inductively coupled plasma mass spectrometry. The main mineral phase is calcic amphibole, and the coexisting glass compositions are tonalite, granodiorite, and granite. The compatibility of rare earth elements increase at 915 °C and then decrease at 970 °C, but the compatibility of most of these elements shows a continued, significant increase with increasing pressure. For high-field strength elements, large ion lithophile elements, actinide compatibility decrease with increasing temperature or pressure, but transition metals show a continued increase in compatibility within the temperature–pressure conditions. From mathematical and graphical fitting, we determined best-fit values for the ideal ionic radius(r_0, 1.01–1.04 ?), the strain-free partitioncoefficient(D_0, 1.18–1.58), and apparent Young's modulus(E, 142–370 GPa) for the M4 site in amphibole according to the lattice strain model. The D_0M4 for rare earth elements rises at 915 °C and then drops at 970 °C at 0.6 GPa.However, the D_0M4 values are positively proportional to the pressure for rare earth elements in the amphibole-glass pairs at 0.6–2.6 GPa and 970 °C. Furthermore, the derived best-fit values for r_0M4 and EM4 are almost constant and trend to increase with rising temperature and pressure,respectively. The partition coefficient is distinctly different for different melt compositions. The rare earth elements become more enriched in amphibole if the quenched glass is granodiorite or granite compared to the tonalitic glasses. 相似文献
5.
The dynamic response of a seismic soil–pile–structure interaction (SSPSI) system is investigated in this paper by conducting nonlinear 3D finite element numerical simulations. Nonlinear behaviors such as non-reflecting boundary condition and soil–pile–structure interaction modeled by the penalty method have been taken into account. An equivalent linear model developed from the ground response analysis and the modified Drucker–Prager model are separately used for soil ground. A comparison of the two models shows that the equivalent linear soil model results in an underestimated acceleration response of the structure under this ground shaking and the soil behavior should be considered as a fully-nonlinear constitutive model in the design process of the SSPSI system. It was also observed that the dynamic response of the system is greatly affected by the nonlinearity of soil–pile interface and is not sensitive to the dilation angle of the soil. Furthermore, the effect of the presence of pile foundations on SSPSI response is also analyzed and discussed. 相似文献
6.
Vir&#; Axelle Xiang Jiansheng Milthaler Frank Farrell Patrick Emmet Piggott Matthew David Latham John-Paul Pavlidis Dimitrios Pain Christopher Charles 《Ocean Dynamics》2012,62(10):1487-1501
Fluid–structure interactions are modelled by coupling the finite element fluid/ocean model ‘Fluidity-ICOM’ with a combined finite–discrete element solid model ‘Y3D’. Because separate meshes are used for the fluids and solids, the present method is flexible in terms of discretisation schemes used for each material. Also, it can tackle multiple solids impacting on one another, without having ill-posed problems in the resolution of the fluid’s equations. Importantly, the proposed approach ensures that Newton’s third law is satisfied at the discrete level. This is done by first computing the action–reaction force on a supermesh, i.e. a function superspace of the fluid and solid meshes, and then projecting it to both meshes to use it as a source term in the fluid and solid equations. This paper demonstrates the properties of spatial conservation and accuracy of the method for a sphere immersed in a fluid, with prescribed fluid and solid velocities. While spatial conservation is shown to be independent of the mesh resolutions, accuracy requires fine resolutions in both fluid and solid meshes. It is further highlighted that unstructured meshes adapted to the solid concentration field reduce the numerical errors, in comparison with uniformly structured meshes with the same number of elements. The method is verified on flow past a falling sphere. Its potential for ocean applications is further shown through the simulation of vortex-induced vibrations of two cylinders and the flow past two flexible fibres.
相似文献7.
Axelle Viré Jiansheng Xiang Frank Milthaler Patrick Emmet Farrell Matthew David Piggott John-Paul Latham Dimitrios Pavlidis Christopher Charles Pain 《Ocean Dynamics》2012,62(10-12):1487-1501
Fluid–structure interactions are modelled by coupling the finite element fluid/ocean model ‘Fluidity-ICOM’ with a combined finite–discrete element solid model ‘Y3D’. Because separate meshes are used for the fluids and solids, the present method is flexible in terms of discretisation schemes used for each material. Also, it can tackle multiple solids impacting on one another, without having ill-posed problems in the resolution of the fluid’s equations. Importantly, the proposed approach ensures that Newton’s third law is satisfied at the discrete level. This is done by first computing the action–reaction force on a supermesh, i.e. a function superspace of the fluid and solid meshes, and then projecting it to both meshes to use it as a source term in the fluid and solid equations. This paper demonstrates the properties of spatial conservation and accuracy of the method for a sphere immersed in a fluid, with prescribed fluid and solid velocities. While spatial conservation is shown to be independent of the mesh resolutions, accuracy requires fine resolutions in both fluid and solid meshes. It is further highlighted that unstructured meshes adapted to the solid concentration field reduce the numerical errors, in comparison with uniformly structured meshes with the same number of elements. The method is verified on flow past a falling sphere. Its potential for ocean applications is further shown through the simulation of vortex-induced vibrations of two cylinders and the flow past two flexible fibres. 相似文献
8.
An unstructured mesh model of the west coast of Britain, covering the same domain and using topography and open boundary forcing
that are identical to a previous validated uniform grid finite difference model of the region, is used to compare the performance
of a finite volume (FV) and a finite element (FE) model of the area in determining tide–surge interaction in the region. Initial
calculations show that although qualitatively both models give comparable tidal solutions in the region, comparison with observations
shows that the FV model tends to under-estimate tidal amplitudes and hence background tidal friction in the eastern Irish
Sea. Storm surge elevations in the eastern Irish Sea due to westerly, northerly and southerly uniform wind stresses computed
with the FV model tend to be slightly higher than those computed with the FE model, due to differences in background tidal
friction. However, both models showed comparable non-linear tide–surge interaction effects for all wind directions, suggesting
that they can reproduce the extensive tide–surge interaction processes that occur in the eastern Irish Sea. Following on from
this model comparison study, the physical processes contributing to surge generation and tide–surge interaction in the region
are examined. Calculations are performed with uniform wind stresses from a range of directions, and the balance of various
terms in the hydrodynamic equations is examined. A detailed comparison of the spatial variability of time series of non-linear
bottom friction and non-linear momentum advection terms at six adjacent nodes at two locations in water depths of 20 and 6 m
showed some spatial variability from one node to another. This suggests that even in the near coastal region, where water
depths are of the order of 6 m and the mesh is fine (of order 0.5 km), there is significant spatial variability in the non-linear
terms. In addition, distributions of maximum bed stress due to tides and wind forcing in nearshore regions show appreciable
spatial variability. This suggests that intensive measurement campaigns and very high-resolution mesh models are required
to validate and reproduce the non-linear processes that occur in these regions and to predict extreme bed stresses that can
give rise to sediment movement. High-resolution meshes will also be required in pollution transport problems. 相似文献
9.
The Nanling Range in South China is well known for its rich granite-related W–Sn deposits.To elucidate the controls of different granite-related W–Sn metallogenesis in the region,we chose five representative orerelated granites(Yanbei,Mikengshan,Tieshanlong,Qianlishan,and Yaogangxian intrusions)in the Hunan–Jiangxi region,and studied their magmatic zircon ages and trace element geochemistry.Our new zircon data showed the differences in ages,temperatures and oxygen fugacity of the ore-forming magmas.Zircon U–Pb ages of the Yanbei and Mikengshan intrusions are characterized by 142.4±2.4 and 143.0±2.3 Ma,respectively,whereas the Tieshanlong and Qianlishan intrusions are 159.5±2.3and 153.2±3.3 Ma,respectively.The Sn-related intrusions were younger than the W-related intrusions.The Tiin-zircon thermometry showed that there was no systematic difference between the Sn-related Yanbei(680–744℃)and Mikengshan(697–763℃)intrusions and the W-related Tieshanlong(730–800℃),Qianlishan(690–755℃)and Yaogangxian(686–751℃)intrusions.However,the zircon Ce^4+/Ce^3+ratios of the Yanbei(averaged at 18.3)and Mikengshan(averaged at 18.8)intrusions are lower than those of the Tieshanlong(averaged at 36.9),Qianlishan(averaged at 38.4)and Yaogangxian(averaged at 37)intrusions,indicating that the Sn-related granitic magmas might have lower oxygen fugacities than those of the W-related.This can be explained by that,in more reduced magmas,Sn is more soluble than W and thus is more enriched in the residual melt to form Sn mineralization.The difference in source materials between the Sn-related and the W-related granites seems to have contributed to the different redox conditions of the melts. 相似文献
10.
In this paper, a fiber beam-column element considering flexure–shear interaction and bond-slip effect is developed for cyclic analysis of reinforced concrete (RC) structures. The element is based on conventional displacement-based Timoshenko beam theory, where the transverse shear deformation is included, and adopts the fiber model to describe the section force–deformation behavior. In the fiber model, shear deformation is assumed to be uniformly distributed along the section and is only resisted by concrete, thus the multi-dimensional concrete damage model is used for concrete fibers and therefore flexure–shear interaction is reflected naturally at the material level. Meanwhile, to account for the significant bond-slip effect at critical regions, the anchorage slip of bars at these regions is analytically derived. Then it is used to modify the uniaxial stress–strain model for steel fibers by assuming that the total strain can be treated as the sum of the bar deformation and anchorage slip, therefore the bond-slip effect is implicitly but simply represented. To validate the proposed element, a series of RC member and structure tests under cyclic loading are simulated. The results indicate that the proposed element can predict cyclic responses of RC structures, and can be used as a reliable tool for analysis of RC structures. 相似文献
11.
A comprehensive dynamic three dimensional finite element model, which includes the effect of lots of important parameters on the micropiles seismic performance, has been presented. The validation of the built model has been carried out using remodeling a single degree of freedom shaking table test done by Mc Manus at the University of Canterbury. The gained results proved the accuracy of the constructed model. Then, using the parametric analysis, effects of all the earthquake characteristics, soil properties, superstructure and micropiles' cap and micropiles structure on the seismic performance of micropiles have been investigated by means of presenting internal forces and displacements which occurred as the main result of earthquake. Furthermore, using the data analysis, the most and the least influential parameters on internal forces are obtained based on the Cosine Amplitude Method (CAM). 相似文献
12.
Quantitative textural analyses including crystal size distributions (CSDs) provide insights into crystallisation kinetics of magmatic systems. Investigations of volcanic crystal textures often rely on greyscale variations on backscattered electron images to identify crystal phases, which must then be thresholded and/or traced manually, a laborious task, and investigations are typically restricted to a single crystal phase. A method is presented that uses energy-dispersive X-ray element maps to generate textural data. Each pixel is identified as a crystal phase, glass or vesicle according to relative chemical composition enabling concurrent acquisition of multiphase CSD, crystallinity and mineral mode data. Data processing is less time intensive for the operator but considerable instrument time is required to generate element maps. The method is applied to 17 dacite samples from the 1980–1986 and 3 from the 2004–2005 eruptive periods of Mount St. Helens volcano (USA) to provide quantitative insights into multiphase textural evolution. All of the CSDs are curved and concave-up in the standard CSD plot with curvature increasing with plagioclase content. To facilitate comparisons with previous studies, CSDs for microlites (<50?μm length crystals) are approximated as straight lines. The line intercepts and slopes provide information on n 0 (nucleation density) and characteristic length or Gτ (the product of growth rate (G) and residence time (τ)), respectively. These parameters, as well as the total groundmass crystallinity, show distinct differences between explosive deposits from summer 1980 and post-summer 1980 domes. Post-summer 1980 microlite n 0 values are mostly at the lower end of the range of those measured for summer 1980 samples. Total groundmass crystallinities during summer 1980 are between 10 and 30?vol.%, whereas post-summer 1980 crystallinity increases to between 39 and 51?vol.%. The range of n 0 values is similar to those previously published for Mount St. Helens but Gτ is consistently higher. Gτ of a May 1985 sample analysed in this study is approximately 2?μm higher compared with previously published data for the same sample when processed using similar methodologies. Groundmass crystallinity data show similar trends to those previously published for the 1980 to 1986 eruption, increasing sharply after summer 1980 then increasing more gradually during the dome-building phase of the eruption. The effects of varying L, the apparent crystal size, and crystal aspect ratio on resultant CSDs are also investigated. Whilst relative temporal variations in CSDs for a given dataset can be investigated, absolute values from different studies cannot be compared unless methods of data acquisition and processing are exactly the same. 相似文献
13.
ＦｉｎｉｔｅｅｌｅｍｅｎｔｓｉｍｕｌａｔｉｏｎｏｆｓｔｅａｄｙｓｔａｔｅＳＨｗａｖｅｍｏｔｉｏｎＺｈｅｎ－ＰｅｎｇＬＩＡＯ（廖振鹏）ａｎｄＧｕａｎｇＹＡＮＧ（杨光）（ＩｎｓｔｉｔｕｔｅｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈａｎｉｃｓ，ＳｔａｔｅＳｅｉｓｍ... 相似文献
14.
《Journal of Volcanology and Geothermal Research》1999,88(4):325-342
This work reports a study on the rare-earth element (REE) behaviour in the high-sulfidation-type superficial alteration and in the phyllic and propylitic alteration in the subsurface of the island of Vulcano. The results of this study demonstrate notable differences in the REE behaviour in the different alteration facies. In silicic alteration, all REE are equally strongly depleted; in advanced argillic alteration, LREE are quite immobile whereas HREE are strongly depleted; in intermediate argillic, phyllic and propylitic alteration, REE patterns are almost unchanged compared with their fresh rock equivalents. These features indicate that the behaviour of REE in hydrothermal alteration facies at La Fossa is mainly controlled by pH, availability of complexing ions in the fluid (mainly SO42− and, in the silicic facies, F−) and the presence of secondary minerals able to host REE in their structures. The origin of the acid fluids is related to the hypogenic introduction of prevalently magmatic acid gases produced by the degassing of the shallow magmatic system of La Fossa volcano. In the active high-temperature fumarolic field of La Fossa, REE behaviour illustrates the transition from a dry vapor system, corresponding to the high-temperature fumaroles, to a wet condensing system in the surrounding area. The results of this research, carried out on an active high-sulfidation system, provide a useful tool for the study of high-sulfidation epithermal ore deposits and particularly help to understand the alteration processes during the main stage of hypogene wall rock alteration. 相似文献
15.
Based on presumed active fault and corresponding model, this paper predicted the near-fault ground motion filed of a scenario earthquake （Mw=6 3/4 ） in an active fault by the explicit finite element method in combination with the source time function with improved transmitting artificial boundary and with high-frequency vibration contained. The results indicate that the improved artificial boundary is stable in numerical computation and the predicted strong ground motion has a consistent characteristic with the observed motion. 相似文献
16.
Dynamic contact theory is applied to simulate the sliding of surface fault.Finite element method is used to analyze the effect of surface fault to site ground motions.Calculated results indicate that amplification effect is obvious in the area near surface fault,especially on the site that is in the downside fault.The results show that the effect of surface fault should be considered when important structure is constructed in the site with surface fault. 相似文献
17.
In this paper, progress in strain study of blocks and faults by GPS data are discussed, and the concept that active structures between blocks are the main body of crustal strain is clarified. By energy transfer principle of elastic mechanics, the relation between strain around faults and tectonic force on fault surfaces is set up and main body element model of crustal strain is constructed. Finally, the relation between mechanical evolution of model and seismogenic process of Kunlun earthquake （Ms=8.1） is discussed by continuous GPS data of datum stations. The result suggests that the relatively relaxed change under background of strong compressing and shearing may help to trigger moderate-strong earthquakes. 相似文献
18.
The concept of structure–soil–structure dynamic interaction was introduced, and the research methods were discussed. Based on several documents, a systematic summary of the history and status of the structure–soil–structure dynamic interaction research that considers adjacent structures was proposed as a reference for researchers. This study is in the initial stage, given its complexity and excessive simplification of the model for soil and structures, and should be carried forward for its significance. An attempt was made to summarize the common major computer programs in this area of study. Furthermore, the advantages, disadvantages, and applicability of such programs were discussed. The existing problems and the future research trend in this field were also examined. 相似文献
19.
Liyan Tian Paterno R. Castillo James W. Hawkins David R. Hilton Barry B. Hanan Aaron J. Pietruszka 《Journal of Volcanology and Geothermal Research》2008
New major and trace element and Sr–Nd isotope data are presented for basaltic glasses from active spreading centers (Central Lau Spreading Center (CLSC), Relay Zone (RZ) and Eastern Lau Spreading Center (ELSC)) in the Central Lau Basin, SW Pacific. Basaltic lavas from the Central Lau Basin are mainly tholeiitic and are broadly similar in composition to mid-ocean ridge basalts (MORB). Their generally high 87Sr/86Sr ratios, combined with relatively low 143Nd/144Nd ratios are more akin to MORB from the Indian rather than Pacific Ocean. In detail, the CLSC, RZ and ELSC lavas are generally more enriched in large ion lithophile elements (Rb, Ba, Sr, and K) than average normal-MORB, which suggests that the mantle beneath the Central Lau Basin was modified by subducted slab-derived components. Fluid mobile/immobile trace element and Sr – Nd isotope ratios suggest that the subduction components were essentially transferred into the mantle via hydrous fluids derived from the subducted oceanic crust; contributions coming from the subducted sediments are minor. Compared to CLSC lavas, ELSC and RZ lavas show greater enrichment in fluid mobile elements and depletion in high field strength elements, especially Nb. Thus, with increasing distance away from the arc, the influence of subduction components in the mantle source of Lau Basin lavas diminishes. The amount of hydrous fluids also influences the degree of partial melting of the mantle beneath the Central Lau Basin, and hence the degree of melting also decreases with increasing distance from the arc. 相似文献
20.
High precision trace element data obtained by inductively coupled plasma mass spectrometry and Sr–Nd isotope analyses are presented for mafic volcanic rocks from Gough Island, South Atlantic. The new data reveal negative Ce anomalies, with Ce/Ce? values in Gough lavas extending down to ～ 0.92. Ce is only fractionated from other rare earth elements (REE) due to formation of Ce4+ under oxidizing conditions of near-surface environments while other REE remain trivalent. Ce anomalies in convergent margin magmas have been shown to indicate a contribution of a subducted sediment component. In contrast, Ce anomalies in intra-plate basalts have been attributed to weathering processes, but can be excluded here based on element–element systematics indicating magmatic trends rather than weathering-induced element mobility. Shallow-level contamination by local marine sediments with negative Ce anomaly inherited from seawater can be excluded because Gough lavas with increasingly negative Ce anomalies do not trend towards low Ce/Pb ratios characterizing such sediments. Rather, it is argued that the negative Ce anomalies in Gough Island lavas are consistent with variable amounts of a sediment component in the mantle plume source. Mixtures between estimates of subducting sediment columns with negative Ce anomaly and mantle capable of giving rise to Gough Island magmas without Ce anomalies reproduce the Gough compositional array with the exception of highly fluid-mobile elements. The calculated trace element composition of the deeply recycled sediment in the Gough plume source is depleted in fluid-mobile elements such as Ba and Pb relative to the composition of some present-day subducting sediments. This loss is attributed to the dehydration or flushing of sediment in the subduction factory, consistent with constraints from arc magmas. 相似文献