Damage Mechanism of Ultra-thin Asphalt Overlay(UTAO) based on Discrete Element Method

Aiming to analyze the damage mechanism of UTAO from the perspective of meso-mechanical mechanism using discrete element method(DEM),we conducted study of diseases problems of UTAO in several provinces in China,and found that aggregate spalling was one of the main disease types of UTAO.A discrete element model of UTAO pavement structure was constructed to explore the meso-mechanical mechanism of UTAO damage under the influence of layer thickness,gradation,and bonding modulus.The experimental results show that,as the thickness of UTAO decreasing,the maximum value and the mean value of the contact force between all aggregate particles gradually increase,which leads to aggregates more prone to spalling.Compared with OGFC-5 UTAO,AC-5 UTAO presents smaller maximum and average values of all contact forces,and the loading pressure in AC-5 UTAO is fully diffused in the lateral direction.In addition,the increment of pavement modulus strengthens the overall force of aggregate particles inside UTAO,resulting in aggregate particles peeling off more easily.The increase of bonding modulus changes the position where the maximum value of the tangential force appears,whereas has no effect on the normal force.

Polarization characteristics and controllability mechanism of passive scattering elements

Polarization feature is one of the important features of radar targets,which has been used in many fields.In this paper,the grid models of some typical foreign moving targets are constructed on the simulation platform,such as glider,cruiser,fixed wing aircraft,and rotorcraft.The electromagnetic scattering characteristics of the moving platforms under the incidence of circular polarization waves are calculated.The typical polarization characteristics which the orthogonal and in-phase components have in the echoes are analyzed and proved.Based on the polarization scattering matrix(PSM)theory,from the point of view of the physical reproduction,the technical status quo that the existing technical approaches are difficult to realize the passive simulation of polarization characteristic of the target is summarized.To solve this problem,combined with the vector synthesis law,the realization mechanism of controllable polarization characteristic of target echoes is proposed,the analytical expressions of polarization control matrix and polarization ratio are deduced,and the controllability of polarization ratio feature in the case of circular polarization is verified by simulation calculation.

An 8-Node Plane Hybrid Element for StructuralMechanics Problems Based on the Hellinger-Reissner Variational Principle

The finite element method (FEM) plays a valuable role in computer modeling and is beneficial to the mechanicaldesign of various structural parts. However, the elements produced by conventional FEM are easily inaccurate andunstable when applied. Therefore, developing new elements within the framework of the generalized variationalprinciple is of great significance. In this paper, an 8-node plane hybrid finite element with 15 parameters (PHQ8-15β) is developed for structural mechanics problems based on the Hellinger-Reissner variational principle.According to the design principle of Pian, 15 unknown parameters are adopted in the selection of stress modes toavoid the zero energy modes.Meanwhile, the stress functions within each element satisfy both the equilibrium andthe compatibility relations of plane stress problems. Subsequently, numerical examples are presented to illustrate theeffectiveness and robustness of the proposed finite element. Numerical results show that various common lockingbehaviors of plane elements can be overcome. The PH-Q8-15β element has excellent performance in all benchmarkproblems, especially for structures with varying cross sections. Furthermore, in bending problems, the reasonablemesh shape of the new element for curved edge structures is analyzed in detail, which can be a useful means toimprove numerical accuracy.

Clay minerals and elemental composition of sediments on different sedimentary units in the northern East China Sea shelf:provenance tracing and genetic mechanism analysis

The composition,provenance,and genetic mechanism of sediment on different sedimentary units of the East China Sea(ECS)shelf are essential for understanding the depositional dynamics environment in the ECS.The sediments in the northern ECS shelf are distributed in a ring-shaped distribution centered on the southwestern Cheju Island Mud.From the inside to the outside,the grain size goes from fine to coarse.Aside from the“grain size effect”,hydrodynamic sorting and mineral composition are important restrictions on the content of rare earth elements(REEs).Based on the grain size,REEs,and clay mineral composition of 300 surface sediments,as well as the sedimentary genesis,the northern ECS shelf is divided into three geochemical zones:southwestern Cheju Island Mud Area(ZoneⅠ),Changjiang Shoal Sand Ridges(ZoneⅡ-1),Sand Ridges of the East China Sea shelf(ZoneⅡ-2).The northern ECS shelf is mostly impacted by Chinese mainland rivers(the Changjiang River and Huanghe River),and the provenance and transport mechanism of sediments of different grain sizes is diverse.The bulk sediments come primarily from the Changjiang River,with some material from the Huanghe River carried by the Yellow Sea Coastal Current and the North Jiangsu Coastal Current,and less from Korean rivers.Among them,surface sediments in the southwestern Cheju Island Mud Area(ZoneⅠ)come mostly from the Changjiang River and partly from the Huanghe River.It was formed by the counterclockwise rotating cold eddies in the northern ECS shelf,which caused the sedimentation and accumulation of the fine-grained sediments of the Changjiang River and the Huanghe River.The Changjiang Shoal Sand Ridges(ZoneⅡ-1)were developed during the early-middle Holocene sea-level highstand.It is the modern tidal sand ridge sediment formed by intense hydrodynamic action under the influence of the Yellow Sea Coastal Current,North Jiangsu Coastal Current,and Changjiang Diluted Water.The surface sediments mainly originate from the Changjiang River and Huanghe River,with the Changjiang River dominating,and the Korean River(Hanjiang River)influencing just a few stations.Sand Ridges of the East China Sea shelf(ZoneⅡ-2)are the relict sediments of the paleo-Changjiang River created by sea invasion at the end of the Last Deglaciation in the Epipleistocene.The clay mineral composition of the surface sediments in the study area is just dominated by the Changjiang River,with the North Jiangsu Coastal Current and the Changjiang Diluted Water as the main transporting currents.

Mechanism of interaction relation between the rare-earth element Ce and impurity elements Pb and Bi in Ag-based filler metal

The mechanism of interaction relation between the rare-earth element Ce and elements Pb and Bi in Ag-based filler metal has been studied. The results show that the compounds CePb and CeBi with high melting point can be easily produced between these three elements in the filler metal, which greatly limited the formation of the isolated phase Pb or Bi and also eliminated the bad effect of impurity elements Pb and Bi on the spreading property of Ag-based filler metal. The metallurgical and quantum-mechanical bond formation analysis show that a strong chemical affinity was existed between the rare-earth element Ce and impurity elements Pb and Bi, which was proved by the XRD analysis results.

Rare earth elements-rich phase and enriching mechanism in sediments from CC area, the Pacific Ocean

Compared to North American shale composition (NASC), REE contents of sediments from the CC area in the Pacific Ocean are obviously high except that cerium has equal content to that of NASC. Three-valence rare earth elements were completely enriched in phosphate-phase and cerium in iron-phase. Rare earth elements in the sediments were originally derived from seawater. During lithi- genic and minerogenic processes of metalliferous nodules, three-valence rare earth elements in sediments mobilized and incorporated into sediments as authigenous biogenic-apatite, while cerium had change from Ce3+ to Ce4+ and directly precipitated from seawater and entered metalliferous nodules and caused Ce anomalies in REE pattern in sediments.

STUDY ON THE SEPARATION AND EXTRACTION MECHANISM OF RARE EARTH ELEMENTS BY MEANS OF REVERSED-PHASE PAPER CHROMATOGRAPHY

Reversed-phase paper chromatography technique is used for study on the extraction mechanism and sep- aration of rare earth elements.As the stationary phase,chromatographic paper strips are impregnated with a solution of monomyristyl phosphoric acid (MPA) in chloroform.Mineral acids are used as developers. The effect of concentration of acids and/or salts upon R_f has been investigated.According to the re- sults of R_f values for a given rare earth element in various acids,the order of extraction ability is HCl>HNO_3>H_2SO_4.A tetrad effect is clearly observed.for the R_f value of rare earth elements.The effects of other parameters on the R_f value,such as the quantities of extractant retained by the paper and the temperature are also examined.Based on the determination of the molar ratio of MPA to rare earth elements and the number of H^+ ions released in extraction reaction,a reasonable mechanism is proposed.The mutual separation of heavy rare earth elements will be better than that of the light rare earth group because of the larger separation coefficient of the former.A mixture of Ho-Er-Tm-Lu is successfully separated by the present method.

Correlations between mineral composition and mechanical properties of granite using digital image processing and discrete element method

This study investigated the correlations between mechanical properties and mineralogy of granite using the digital image processing(DIP) and discrete element method(DEM). The results showed that the X-ray diffraction(XRD)-based DIP method effectively analyzed the mineral composition contents and spatial distributions of granite. During the particle flow code(PFC2D) model calibration phase, the numerical simulation exhibited that the uniaxial compressive strength(UCS) value, elastic modulus(E), and failure pattern of the granite specimen in the UCS test were comparable to the experiment. By establishing 351 sets of numerical models and exploring the impacts of mineral composition on the mechanical properties of granite, it indicated that there was no negative correlation between quartz and feldspar for UCS, tensile strength(σ_(t)), and E. In contrast, mica had a significant negative correlation for UCS, σ_(t), and E. The presence of quartz increased the brittleness of granite, whereas the presence of mica and feldspar increased its ductility in UCS and direct tensile strength(DTS) tests. Varying contents of major mineral compositions in granite showed minor influence on the number of cracks in both UCS and DTS tests.

Effects of Alloying Elements on the Microstructures and Mechanical Properties of Heavy Section Ductile Cast Iron

The effects of alloying elements on the as-cast microstructures and mechanical properties of heavy section ductile cast iron were investigated to develop press die material having high strength and high ductility. Measurements of ultimate tensile strength, 0.2% proof strength, elongation and unnotched Charpy impact energy are presented as a function of alloy amounts within 0.25 to 0.75 wt pct range. Hardness is measured on the broken tensile specimens. The small additions of Mo, Cu, Ni and Cr changed the as-cast mechanical properties owing to the different as-cast matrix microstructures. The ferrite matrix of Mo and Ni alloyed cast iron exhibits low strength and hardness as well as high elongation and impact energy. The increase in Mo and Ni contents developed some fractions of pearlite structures near the austenite eutectic cell boundaries, which caused the elongation and impact energy to drop in a small range. Adding Cu and Cr elements rapidly changed the ferrite matrix into pearlite matrix, so strength and hardness were significantly increased. As more Mo and Cr were added, the size and fraction of primary carbides in the eutectic cell boundaries increased through the segregation of these elements into the intercellular boundaries.

Numerical investigations on mechanical characteristics and failure mechanism of outwash deposits based on random meso-structures using discrete element method

Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.

Effect of RE elements on the microstructural and mechanical properties of as-cast and age hardening processed Mg-4Al-2Sn alloy

In the present article,the effect of Rare Earth elements on the microstructural development and mechanical properties of as cast and age treated Mg-4Al-2Sn(AT42)alloy is studied.Investigation has been conducted by optical and scanning electron microscope,XRD and tensile tests.Analysis of the data showed that alloy’s dendrites turn into larger dendritic structure with sharp and narrow arms from equiaxed rosette type by addition of RE elements.In contrast to the base alloy,aging treatment shows a positive effect on the mechanical properties(yield strength,tensile strength and elongation)of AT42+1RE alloy mainly because of retention of the thermally stable RE containing intermetallics as strong barriers to grain growth.Also,increase of solute aluminum due to the decomposition of Mg 17 Al 12 along with saturated RE elements led to formation of blocky shape Al 2 RE in the microstructure during aging which enhanced the mechanical properties.It was found that the best result(yield of 70 MPa,tensile strength of 168 MPa and elongation of 14%)could be achieved by aging the AT42+1RE alloy at 443 K(170℃)for 8 h.However,mechanical properties of AT42+1RE alloy starts to decrease after exceeding its optimum aging conditions due to the coarsening of intermetallics.

Investigation of the block toppling evolution of a layered model slope by centrifuge test and discrete element modeling

Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.

Effects of rare-earth elements on the glass-forming ability and mechanical properties of Cu_(46)Zr_(47-x)Al_7M_x(M = Ce,Pr,Tb,and Gd) bulk metallic glasses

Cu46Zr47-xA17Mx （M = Ce, Pr, Tb, and Gd） bulk metallic glassy （BMG） alloys were prepared by copper-mold vacuum suction casting. The effects of rare-earth elements on the glass-forming ability （GFA）, thermal stability, and mechanical properties of Cu46Zr47-xA17Mx were investigated. The GFA of Cu46Zr47-xA17Mx （M = Ce, Pr） alloys is dependent on the content of Ce and Pr, and the optimal content is 4 at.%. Cu46Zr47-xA17Thx（X = 2, 4, and 5） amorphous alloys with a diameter of 5 mm can be prepared. The GFA of Cu46Zr47-xA17Gdx（x = 2, 4, and 5） increases with increasing Gd. Tx and Tp of all decrease. Tg is dependent on the rare-earth element and its content. ATx for most of these alloys decreases except the Cu46Zra2Al7Gd5 alloy. The activation energies △Eg, △Ex, and △Ep for the Cu46Zr42A17Gd5 BMG alloy with Kissinger equations are 340.7, 211.3, and 211.3 kJ/mol, respectively. These values with Ozawa equations are 334.8, 210.3, and 210.3 kJ/mol, respec- tively. The Cu46Zr45Al7Tb2 alloy presents the highest microhardness, Hv 590, while the Cu46Zr43A17Pr4 alloy presents the least, Hv 479. The compressive strength （at.f.） of the Cu46Zra3A17Gd4 BMG alloy is higher than that of the Cu46Zr43Al7Tb4 BMG alloy.

Analysis of mesoscopic mechanical dynamic characteristics of ballast bed with under sleeper pads

The meso-dynamical behaviour of a high-speed rail ballast bed with under sleeper pads(USPs)was studied.The geometrically irregular refined discrete element model of the ballast particles was constructed using 3D scanning techniques,and the 3D dynamic model of the rail-sleeper-ballast bed was constructed using the coupled discrete element method-multiflexible-body dynamics(DEM-MFBD)approach.We analyse the meso-mechanical dynamics of the ballast bed with USPs under dynamic load on a train and verify the correctness of the model in laboratory tests.It is shown that the deformation of the USPs increases the contact area between the sleeper and the ballast particles,and subsequently the number of contacts between them.As the depth of the granular ballast bed increases,the contact area becomes larger,and the contact force between the ballast particles gradually decreases.Under the action of the elastic USPs,the contact forces between ballast particles are reduced and the overall vibration level of the ballast bed can be reduced.The settlement of the granular ballast bed occurs mainly at the shallow position of the sleeper bottom,and the installation of the elastic USPs can be effective in reducing the stress on the ballast particles and the settlement of the ballast bed.

Discrete element analysis on press and fracture mechanism of propellant grain

To analyze fracture mechanism of propellant grain and study the mechanical properties of propellant grain, the press and fracture processes of propellant grain with and without initial defects are modeled using the discrete element method. On the basis of the appropriate constitutive relationships, the discrete element model of the propellant grain was established. Compared with experimental measurements, the micro-parameters of the bonded-particle model of the propellant grain under unconfined uniaxial compression tests were calibrated. The propellant grains without initial defects, with initial surface defects, and with initial internal defects were studied numerically through a series of unconfined uniaxial compression tests. Results show that the established discrete element model is an efficient tool to study the press and fracture processes of the propellant grain. The fracture process of the propellant grain without initial defects can be divided into the elastic deformation phase, crack initiation phase, crack stable propagation phase, and crack unstable propagation phase. The fracture mechanism of this grain is the global shear failure along the direction of the maximum shear stress. Initial defects have significant effects on both the fracture mechanism and peak strength of the propellant grain. The major fracture mechanism of the propellant grain with initial surface defects is local shear failure, whereas that of the propellant grain with initial internal defects is global tensile failure. Both defects weaken the peak strengths of the propellant grain. Therefore, the carrying and filling process of the propellant grain needs to minimize initial defects as far as possible.

Formation Mechanism of NiAl/TiB_2 Nanocomposite by Mechanical Alloying Elemental Powders

A NiAl/TiB2 nanocomposite is synthesized by mechanical alloying elemental powders. Upon milling for a certain time, an abrupt exothermic reaction occurs and a large amount of NiAl and TiB2 compounds form simultaneously. It is suggested that two separate chemical reactions,i.e. Ni+Al →NiAl and Ti+2B→TiB2, are involved during the exothermic reaction. Additionof Ti and B to Ni-Al system impedes the structural evolution of Ni and Al powders and delays the abrupt reaction. The final products are equilibrium phases without any metastable phases formed. This type of reaction is suggested to be suitable for alloy systems with two large heatrelease reactions.

Effect of Size Change on Mechanical Properties of Monolayer Arsenene

The Young's modulus, shear modulus and Poisson's ratio of monolayer arsenene with different sizes were calculated by finite element method, so as to explore the influence of dimension and orientation on the mechanical properties of monolayer arsenene. The calculation results show that the small size has a significant effect on the mechanical properties of the monolayer arsenene. The smaller the size, the larger the Young's modulus and Poisson's ratio of the monolayer arsenene. The size change has a great influence on the Young's modulus of the arsenene handrail direction, and the Young's modulus of the zigzag direction is not sensitive to the size change. Similarly, the size change has a significant effect on the shear modulus of arsenene in the handrail direction, while the shear modulus in the zigzag direction has no significant effect on its size change. For the Poisson's ratio, the situation is just the opposite, and the effect of the size change on the Poisson's ratio of the arsenene zigzag direction is greater than that of the handrail direction.

Sources, enrichment mechanisms, and resource effects of rare metal elements-enriched geothermal springs in Xizang, China

Rare metals such as lithium(Li), rubidium(Rb), and cesium(Cs) are strategically crucial mineral resources for the development of emerging industries in China. Ensuring a stable long-term supply of these resources is essential. The geothermal systems in Xizang, China are well-developed, with a wide distribution of various types. Most high-temperature geothermal systems in Xizang are exceptionally enriched in rare metal elements(RMEs) and have the potential to become a new source of rare metals to secure China's strategic mineral resource supply in the future. A close relationship also exists between the geothermal system and the special salt lake resources on the Tibetan Plateau. Geothermal springs thus play a key role in the migration and enrichment of RMEs from deep to shallow parts of the crust, in the transition between endogenous and exogenous mineralization, and source-to-sink processes. However, the mechanisms of element enrichment and evolution in these springs have not been systematically discussed, and many theoretical issues remain to be investigated. Based on summarizing and analyzing previous research, this study employs hydrochemical and isotopic geochemistry methods to investigate typical geothermal springs across Xizang and explore the anomalous enrichment mechanism of RMEs, and the resource effects of geothermal springs. Comprehensive analysis shows that the total dissolved solids(TDS) and hydrochemical types of geothermal springs are similar to those of major geothermal fields worldwide, but the Tibetan springs are abnormally rich in Li(averaging5.48 mg/L), Rb(averaging 0.75 mg/L), and Cs(averaging 3.58 mg/L), which are hundreds to thousands of times more concentrated than natural waters. The distribution of these enriched geothermal springs is controlled by the Yarlung Zangbo suture zone and the extended N-S trending rifts, especially in the intersection zone of the two, where the geothermal springs are the most enriched. Based on the spatial distribution, isotopic, and elemental geochemistry, the RMEs enriched in Tibetan geothermal springs are mainly derived from the magmatic-hydrothermal fluids generated by the partial melting of the subducted Indian plate under the Eurasian continent. These fluids not only maintain geothermal activities as a heat source but also participate in the material cycle of the geothermal spring as a material source. Against the background of regional crustal enrichment in RMEs,incompatible elements such as Li, Rb, and Cs are gradually enriched in magmatic-hydrothermal processes including partial melting in the source, magmatic differentiation, and hydrothermal fluid exsolution, and some ore-forming elements are further extracted from surrounding rocks through deep high-temperature water-rock interactions. Eventually, an eruption occurs, and these fluids move to the surface to form a geothermal spring rich in RMEs. With the drainage of geothermal springs, the RMEs are continuously transported to the lake basin by surface runoff and continue to concentrate and evolve into salt lake brines under an extremely arid climate environment, constituting an endogenous source and exogenous accumulation salt lake metallogenic model. This comprehensive explanation of the sources, migration, enrichment mechanisms, and resource effects of geothermal springs will deepen the understanding of rare metal mineralization processes, and aid in the advancement of theoretical models for key rare metal mineral resources in various geological bodies of the Tibetan Plateau, significantly expanding exploration scopes and accurately assessing the resource potential of RMEs.

Study on Dynamic Mechanical Behavior of Al-Mg-Si Alloy

The dynamic mechanical behavior of Al-Mg-Si alloy was investigated under different strain rates by mechanical property and microstructure characterization,constitutive behavior analysis and numerical simulation in the present study.As the strain rate increases,the yield strength,ultimate tensile strength and elongation increase first,then remain almost constant,and finally increase.The alloy always exhibits a typical ductile fracture mode,not depending on the strain rate.However,as the strain rate increases,the number of dimples gradually increases.Tensile deformation can refine grains,however,the grain structure is slightly affected by the strain rate.An optimized Johnson-Cook constitutive equation was used to describe the mechanical behavior and obtained by fitting the true stress-strain curves.The parameter C was described by a function related to the strain rate.The fitting true stress-strain curves by the JC model agree very well with the experimental true stress-strain curves.The true stress-strain curves calculated by the finite element numerical simulation agree well with the experimental true stress-strain curves.

Trace element compositions of pyrite and stibnite:implications for the genesis of antimony mineralization in the Yangla Cu skarn deposit,Northwestern Yunnan,China

The Yangla Cu skarn deposit is located in the central part of the Jinshajiang Suture Zone,southwest China,with a total reserve of 150 Mt Cu@1.03%.The newly discovered antimony orebodies at the depth of Yangla are strictly controlled by the stratum,structure,and lithology,which are lenticular and vein-like within the marble fracture zone,which can provide a window into multistage miner-alization and ore genesis at Yangla.Mineralization can be divided into three types,Cu–Pb–Zn(skarn)pyrite,galena,and sphalerite,Cu(porphyry)chalcopyrite and pyrite,and Sb(hydrothermal)stibnite and pyrite.The mineral assem-blages were stibnite+pyrite+calcite+quartz±minor scheelite in antimony ores.This study presents quantitative measurements of the trace element compositions of pyrite and stibnite from the Yangla antimony ores.Analysis of pyrite with electron probe microanalysis(EPMA)showed enrichment in Co,Ni,Sb,As,and Mo,and deficit in its S and Fe contents when compared to the stoichiometric con-centrations of S and Fe in pyrite.The Sb-related pyrite may belong to sedimentary-reworked genesis and may be modi-fied by hydrothermalfluids,thereby presenting a certain dif-ference(i.e.,crystal morphology,texture,and chemical com-position)compared to the skarn and porphyry Cu-related pyrite in the Yangla Cu skarn deposit.Analysis of stibnite with EPMA and inductively coupled plasma-mass spectrom-etry showed enrichment in As,Pb,Sn,Pb,Cu,and Zn,and presented much higher Sb contents and slightly lower S con-tents when compared to the stoichiometric concentrations of Sb and S in stibnite.Statistical analysis of the stibnite trace elements showed correlations for the elemental pairs Cu–Pb,As–Sb,and Sn–Pb,and the coupled substitution equations Sb^(3+)↔Cu^(+)+Pb^(2+),Sb^(3+)↔As^(3+),and Sn^(2+)↔Pb^(2+)may be the major factors governed the incorporating Cu,Pb,As and Sn within the stibnite.Moreover,this study preliminary shows that the antimony mineralization may belong to a car-bonate replacement hydrothermal genesis at Yangla.