基于D-S证据理论的岩爆预测方法研究

为了有效预测岩爆,提出基于D-S证据理论的岩爆预测方法.首先,选取与岩爆发生相关的6个指标因素作为证据体,并通过模糊物元框架和正态型隶属度函数构建证据体的基本概率分配.然后,利用K均值将证据体分类,并提出簇内证据用传统方式融合而簇间证据用权重方式融合的组合融合规则,以减轻高冲突证据融合的不利影响.最后,将模型应用在秦岭终南山公路隧道2号竖井工程,且与经验方法对比.为了分析预测过程的不确定性和估计岩爆发生概率,采用蒙特卡洛模拟进行抽样仿真,并通过Spearman秩相关系数衡量输入指标的全局敏感性.研究结果表明:输入指标在不同的岩爆案例的影响程度差异较大且方向不同;5个岩爆案例的发生概率在40.8%~70.1%之间.该模型表现出优异的预测分类性能,可为深埋地下工程岩爆预测提供参考.

Analysis of ground vibrations due to underground trains by 2.5D finite/infinite element approach

The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, are modeled by the finite elements, and the wave propagation properties of the far field extending to infinity are modeled by the infinite elements. One particular feature of the 2.5D approach is that it enables the computation of the three-dimensional response of the half-space, taking into account the load-moving effect, using only a two-dimensional profile. Although the 2.5D finite/infinite element approach shows a great advantage in studying the wave propagation caused by moving trains, attention should be given to the calculation aspects, such as the rules for mesh establishment, in order to avoid producing inaccurate or erroneous results. In this paper, some essential points for consideration in analysis are highlighted, along with techniques to enhance the speed of the calculations. All these observations should prove useful in making the 2.5D finite/infinite element approach an effective one.

Trace-element geochemistry and S–O isotopes in the fluorite-barite mineralization of Merguechoum,Moroccan eastern Meseta:insights into ore genesis to the Pangea rifting

The Merguechoum fluorite-barite mineralization,located in the Eastern Meseta of Morocco,is hosted in the Late Hercynian granite.The ore consists of fine crystals of fluorite 1,massive barite 1,euhedral crystals of fluorite 2,and barite 2 with calcite and minor quartz and sulfides.The Merguechoum ore deposits have never been investigated.This study was the first contribution that studied the genesis of fluorite and barite.The ore occurs as dissemination within granite intrusion and also fills the NE-SWtrending meter-sized fractures and faults.The values of the total Rare Earth Elements and Yttrium(REY)and the ratios of LREY/HREY,Y/Ho,Tb/Ca,and Tb/La indicate that the Merguechoum fluorite precipitated from hydrothermal fluids,likely basinal brines,which interacted with the Hercynian granite.The REY data indicate that the ore-forming fluids of the early stage have intensely interacted with the Hercynian granite compared to those of the late ore stage.The gradual decrease in the europium(Eu/Eu^(*)),yttrium(Y/Y^(*)),and cerium(Ce/Ce^(*))anomalies and a low concentration ofΣREY observed in the second ore stage compared to the first ore stage suggest an increase in p H and fO_(2)and by inference a decrease in temperature during the evolution of the hydrothermal system.This evolution could be explained by fluid mixing between the ascending basinal hydrothermal fluids and the diluted sulfate-rich meteoric water barite separates from selected samples reveal that the dissolved sulfates(SO_(4)^(2-))were derived from Permian–Triassic sulfates and/or coeval poreseawater sulfates.The proposed fluid mixing triggered the precipitation of an early-stage F-Ba assemblage followed by the second-stage F-Ba mineralization.Geologic fieldwork,REY inventories,and isotope data point to the ore genesis during the Permian–Triassic extensional tectonic activity concerning the Pangea rifting.This extensional tectonic environment is likely the driving force that mobilized a large amount of the ore-forming basinal brines along the available faults and fractures to the loci of ore deposition.

Computing of the Anchor by the Method of Three-Dimension Point-Radiate Infinite Elements

On the basis of the one-dimension infinite element theory, the coordinate translation and shape function of 3D point-radiate 8-node and 4-node infinite elements are derived. They are coupled with 20-node and 8-node finite elements to compute the compression distortion of the prestressed anchorage segment. The results indicate that when the prestressed force acts on the anchorage head and segment, the stresses and the displacements in the rock around the anchorage head and segment concentrate on the zone center with the anchor axis, and they decrease with exponential forms. Therefore,the stresses and the displacement spindles are formed. The calculating results of the infinite element are close to the theoretical results. This indicates the method is right. This article introduces a new way to study the mechanism of prestressed anchors. The obtained results have an important role in the research of the anchor mechanism and engineering application.

INFINITE ELEMENT METHOD FOR PROBLEMS ON UNBOUNDED AND MULTIPLY CONNECTED DOMAINS

The author studies the infinite element method for the boundary value problems of second order elliptic equations on unbounded and multiply connected domains. The author makes a partition of the domain into infinite number of elements. Without dividing the domain, as usual, into a bounded one and an exterior one, he derives an initial value problem of an ordinary differential equation for the combined stiffness matrix, then obtains the approximate solution with a small amount of computer work. Numerical examples are given.

ADAPTIVE ELLIPSOIDAL ACOUSTIC INFINITE ELEMENT

It is shown that the basis of the ellipsoidal acoustic infinite elementBurnett method, the multipole expansion, cannot represent real ellipsoidal acoustic field exactly.To solve the problem, a weight of angular direction is added to the multipole expansion. Thecomparison of the modified method and the prime method shows that the modified method can describeand solve the ellipsoidal acoustic field more accurately than ever. A dilating sphere is used totest the new method further. Unlike other infinite element methods, varied ratio of the ellipsoidalartificial boundary instead of sphere is used. The pressure value of the artificial boundary isutilized as the initial value of the new method. Then the radiating phenomena of the ellipsoidalacoustic field can be researched using the new method. These examples show the feasibility of theadaptive method.

四川康定偏岩子造山型金矿床成因:石英微量元素和硫化物S同位素证据

大渡河金矿带是中国重要的金矿化富集区,发育一系列造山型金矿床。偏岩子金矿床位于大渡河南段的康定地区,矿体赋存于震旦系灯影组白云岩中。金成矿过程分为金-黄铁矿-石英阶段、金-氟镁石阶段、金-硫化物-氟镁石-石英阶段和黄铁矿-白云石-白云母阶段4个成矿阶段。通过石英中微量元素和硫化物S同位素的原位分析,旨在约束成矿流体性质、成矿物质来源和矿床成因。石英中Li,B,Na,Al,P和K相对富集,主要以Al^(3+)^(+)(Li^(+),Na^(+),K^(+),H^(+))Si^(4+)和Al^(3+)^(+)P 5^(+)2Si^(4+)的替换机制进入石英晶格。石英中Al与Li呈明显正相关关系,指示其形成温度较低。围岩中硫化物的δ^(34) S值介于-22.3‰~-13.1‰之间,表明S来源于海相硫酸盐的细菌还原作用;矿石中硫化物的δ^(34) S值范围为5.4‰~7.0‰,与交代岩石圈地幔的硫同位素值一致。综合矿床地质特征、石英微量元素和硫化物硫同位素组成,研究认为偏岩子属于造山型金矿床。成矿流体在上升过程中与赋矿白云岩接触,发生水-岩反应,导致自然金、含金硫化物、氟镁石和萤石的沉淀结晶。

Simultaneous nitrification and autotrophic denitrification in fluidized bed reactors using pyrite and elemental sulfur as electron donors

In this study, simultaneous nitrification and autotrophic denitrification (SNAD) with either elemental sulfur or pyrite were investigated in fluidized bed reactors in mesophilic conditions. The reactor performance was evaluated at different ammonium (12-40 mg/L of NH4+-N), nitrate (35-45 mg/L of NO3--N), and dissolved oxygen (DO) (0.1-1.5 mg/L) concentrations, with a hydraulic retention time of 12 h. The pyrite reactor supported the SNAD process with a maximum nitrogen removal efficiency of 139.5 mg/(L·d) when the DO concentration was in the range of 0.8-1.5 mg/L. This range, however, limited the denitrification efficiency of the reactor, which decreased from 90.0% ± 5.3% in phases II-V to 67.9% ± 7.2% in phases VI and VII. Sulfate precipitated as iron sulfate (FeSO4/Fe2(SO4)3) and sodium sulfate (Na2SO4) minerals during the experiment. The sulfur reactor did not respond well to nitrification with a low and unstable ammonium removal efficiency, while denitrification occurred with a nitrate removal efficiency of 97.8%. In the pyrite system, the nitrifying bacterium Nitrosomonas sp. was present, and its relative abundance increased from 0.1% to 1.1%, while the autotrophic denitrifying genera Terrimonas, Ferruginibacter, and Denitratimonas dominated the community. Thiobacillus, Sulfurovum, and Trichlorobacter were the most abundant genera in the sulfur reactor during the entire experiment.

基于S-R和分解定理的二维几何非线性问题的虚单元法求解

应变-旋转(Strain-Rotation,S-R)和分解定理为分析几何非线性问题提供了合理可靠的理论基础,但用有限元求解时会遇到大变形发生后的网格畸变问题。近年提出的虚单元法(Virtual element method,VEM)适用于一般的多边形网格,因此,该文尝试使用一阶虚单元求解基于S-R和分解定理的二维几何非线性问题,以克服网格畸变的影响。基于重新定义的多项式位移空间基函数,推演获得一阶虚单元分析线弹性力学问题时允许位移空间向多项式位移空间的投影表达式;按照虚单元法双线性格式的计算规则,分析处理基于更新拖带坐标法和势能率原理的增量变分方程;进而建立离散系统方程及其矩阵表达形式,并编制MATLAB求解程序;采用常规多边形网格和畸变网格,应用该文算法分析均布荷载下的悬臂梁和均匀内压下的厚壁圆筒变形。结果与已有文献和ANSYS软件的对比表明:该文算法在两种网格中均可有效执行且具备足够数值精度。总体该文算法为基于S-R和分解定理的二维几何非线性问题求解提供了一种鲁棒方法。

小红石砬子铅锌(银)矿床成因:来自S-Pb同位素和元素特征的证据

小红石砬子铅锌(银)矿床是吉中-延边铜镍多金属成矿带中具有代表性的中型矿床,前人关于该矿床类型划分的观点有火山成因块状硫化物型(Volcanogenic Massive Sulfide,VMS)、中-低温热液脉型、喷流沉积-热液叠加改造型,这三种观点的争论严重制约了对该矿区及区域成矿和找矿的认识。本文选取该矿床中金属硫化物开展原位微量和S-Pb同位素分析,赋矿围岩开展全岩S同位素分析。闪锌矿Fe-Zn含量(4.48%~10.33%、53.67%~62.12%)与中温热液矿床一致(3.00%~10.00%、50.00%~60.00%,成矿温度约200℃),Zn/Cd比值(53.68~76.75)同样指示中低温成矿作用(高温:Zn/Cd>500;低温:Zn/Cd<250),Ga/Ge比值在lg(Ga/Ge)-T图解中,对应温度约为100~180℃,上述温度范围与前人测定的流体包裹体均一温度基本吻合(90~250℃)。闪锌矿Zn/Cd的比值(53.68~76.75)与Cd含量(7577×10-6~11266×10-6)特征与密西西比河谷型(Mississippi Valley-Type,MVT)矿床一致;黄铁矿Fe/S比值(0.76~1.29,平均值为0.89)接近典型MVT矿床的理论值(0.88)。金属硫化物δ34 S V-CDT值介于-9.6‰~-0.5‰之间,总体偏负;赋矿围岩δ34 S V-CDT值介于-8.5‰~4.2‰之间,两者重叠范围较大(主要集中于-9.0‰~-2.0‰),说明围岩提供了主要的还原硫,而古老基底可能提供了部分较轻的硫。金属硫化物Pb同位素比值为206 Pb/204 Pb=18.30~18.39、207 Pb/204 Pb=15.61~15.63、208 Pb/204 Pb=38.38~38.44,指示铅为混合来源,成矿物质可能主要来自上地壳和造山带,模式年龄计算显示上太古界乌拉山群变质岩可能为铅的初始来源。综上所述,小红石砬子铅锌(银)矿床地球化学特征与MVT型矿床基本一致,地质特征略有差别,因此推测其为类MVT型矿床。

A NOVEL ELLIPSOIDAL ACOUSTIC INFINITE ELEMENT

A novel ellipsoidal acoustic infinite element is proposed. It is based a new pressure representation, which can describe and solve the ellipsoidal acoustic field more exactly. The shape functions of this novel acoustic infinite element are similar to the (Burnett's) method, while the weight functions are defined as the product of the complex conjugates of the shaped functions and an additional weighting factor. The code of this method is cheap to generate as for 1-D element because only 1-D integral needs to be numerical. Coupling with the standard finite element, this method provides a capability for very efficiently modeling acoustic fields surrounding structures of virtually any practical shape. This novel method was deduced in brief and the conclusion was kept in detail. To test the feasibility of this novel method efficiently,in the examples the infinite elements were considered,excluding the finite elements relative. This novel ellipsoidal acoustic infinite element can deduce the analytic solution of an oscillating sphere. The example of a prolate spheroid shows that the novel infinite element is superior to the boundary element and other acoustic infinite elements. Analytical and numerical results of these examples show that this novel method is feasible.

THE h-p VERSION OF THE INFINITE ELEMENT METHOD AND ITS APPLICATION

The paper investigates the h-p version of the infinite element method. An exponential conver gence can be abtained by a small amount of computing work.

PLANE INFINITE ANALYTICAL ELEMENT AND HAMILTONIAN SYSTEM

It is not convenient to solve those engineering problems defined in an infinite field by using FEM. An infinite area can be divided into a regular infinite external area and a finite internal area. The finite internal area was dealt with by the FEM and the regular infinite external area was settled in a polar coordinate. All governing equations were transformed into the Hamiltonian system. The methods of variable separation and eigenfunction expansion were used to derive the stiffness matrix of a new infinite analytical element.This new element, like a super finite element, can be combined with commonly used finite elements. The proposed method was verified by numerical case studies. The results show that the preparation work is very simple, the infinite analytical element has a high precision, and it can be used conveniently. The method can also be easily extended to a three-dimensional problem.

S对铸造高温合金K417G凝固行为的影响

在1340~1230℃采用等温凝固实验,研究了w[S]分别为0.0001%、0.0018%和0.0037%对K417G合金凝固特性的影响。结果表明,S降低合金的固相线温度,扩大固液两相区,当w[S]从0.0001%增加到0.0037%时,合金的固相线温度由1255℃降低至1245℃。在凝固过程中S元素偏析于液相,并在终凝区析出Y相,且S促进Al、Cr元素在固液界面处的富集。在相同等温温度下,S元素对MC碳化物和γ+γ’共晶的析出温度与形貌无明显影响,但S含量的增加促使MC碳化物、γ+γ’共晶和Y相数量的增加,且Y相的析出温度从1230℃提升至1240℃。

INFINITE SPLINE BOUNDARY ELEMENT METHOD AND ITS APPLICATION TO STRUCTURE-FOUNDATION INTERACTION

In view of the infinity behaviors of 3-D Kelvin solution, we constructed an infinite spline boundary element which has fine precision in the analysis of the half space foundation subjected to uniform pressure on the circular domain. We also analysed a square plate resting on elastic half space foundation. The results indicate that this model not only fits for the coupled analysis of foundation and structures but also has the advantage of fewer degrees of freedom and fine precision.

Eigenvalue Analysis of Thin Plate with Complicated Shapes By a Novel Infinite Element Method

A novel infinite element method(IEM)is presented for solving plate vibration problems in this paper.In the proposed IEM,the substructure domain is partitioned into multiple layers of geometrically similar finite elements which use only the data of the boundary nodes.A convergence criterion based on the trace of the mass matrix is used to determine the number of layers in the IE model partitioning process.Furthermore,in implementing the Craig-Bampton(CB)reduction method,the inversion of the global stiffness matrix is calculated using only the stiffness matrix of the first element layer.The validity and performance of the proposed method are investigated by means of four illustrative problems.The first example considers the case of a simple clamped rectangular plate.It is observed that the IEM results are consistent with the theoretical results for first six natural frequencies.The second example considers the frequency response of a clamped rectangular plate with a crack.The main feature of IEM is that a very fine and good quality virtual mesh can be created around the crack tip.The third and fourth examples consider the natural frequency of a multiple point supported plate and a perforated plate,respectively.The results are obtained just need to adjust the reference point or boundary nodes.The parametric analyses for various geometric profiles are easy to be conducted using these numerical techniques.In general,the results presented in this study have shown that the proposed method provides a direct,convenient and accurate tool for eigenvalue analysis of thin plate structure with complicated shapes.

Finite Element Analysis of the Base

This paper describes the structure of the base,on which two working arms are installed simultaneously.To ensure structura safety,the fatigue failure analysis and statics analysis are finished using the finite element method.The calculation can make sure that the structure of the base meets the design standard,and the material can be reduced one grade.

On the Total Dynamic Response of Soil-Structure Interaction System in Time Domain Using Elastodynamic Infinite Elements with Scaling Modified Bessel Shape Functions

This paper is devoted to a new approach—the dynamic response of Soil-Structure System (SSS), the far field of which is discretized by decay or mapped elastodynamic infinite elements, based on scaling modified Bessel shape functions are to be calculated. These elements are appropriate for Soil-Structure Interaction problems, solved in time or frequency domain and can be treated as a new form of the recently proposed elastodynamic infinite elements with united shape functions (EIEUSF) infinite elements. Here the time domain form of the equations of motion is demonstrated and used in the numerical example. In the paper only the formulation of 2D horizontal type infinite elements (HIE) is used, but by similar techniques 2D vertical (VIE) and 2D corner (CIE) infinite elements can also be added. Continuity along the artificial boundary (the line between finite and infinite elements) is discussed as well and the application of the proposed elastodynamical infinite elements in the Finite element method is explained in brief. A numerical example shows the computational efficiency and accuracy of the proposed infinite elements, based on scaling Bessel shape functions.

ESD保护电路在HDMI板级信号完整性中的影响分析及其布局优化研究

为了解决HDMI接入ESD保护电路后信号完整性受破坏等问题,从器件的空间布局对HDMI信号完整性影响进行研究分析,同时考虑了瞬态电压抑制(TVS)高频寄生参数的影响,搭建了ESD放电模型和TVS高频等效电路模型,并对其可靠性进行了验证。从差分器件接入信号线旋转角度和彼此间错开距离研究其对信号完整性的影响,设计了25套不同夹角和4套不同错开距离的板级模型,在不同特性的传输频率下进行S参数仿真,并从中选取出垂直型、水平型、错开型三种具有代表性的空间布局模型,利用有限元仿真得到差分信号线的S参数和眼图。仿真结果表明,垂直型排布相比于其他两种典型空间布局,回波损耗平均降低了248.1%,插入损耗平均降低了20.6%,眼图的眼宽和眼高在三种空间布局中最大。研究成果为PCB静电放电保护电路分析与设计提供了布局优化指导。