Fracture mechanics of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in an infinite piezoelectric or on the interface of piezoelectric bimaterials. For homogeneous materials it is found that the normal electric displacement D-2, induced by the crack, is constant along the crack faces which depends only on the remote applied stress fields. Within the crack slit, the perturbed electric fields induced by the crack are also constant and not affected by the applied electric displacement fields. For bimaterials, generally speaking, an interface crack exhibits oscillatory behavior and the normal electric displacement D-2 is a complex function along the crack faces. However, for bimaterials, having certain symmetry, in which an interface crack displays no oscillatory behavior, it is observed that the normal electric displacement D-2 is also constant along the crack faces and the electric field E-2 has the singularity ahead of the crack tip and has a jump across the interface. Energy release rates are established for homogeneous materials and bimaterials having certain symmetry. Both the crack front parallel to the poling axis and perpendicular to the poling axis are discussed. It is revealed that the energy release rates are always positive for stable materials and the applied electric displacements have no contribution to the energy release rates.

Crack problems of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in piezoelectrics or on the interfaces of piezoelectric bimaterials. A class of boundary problems involving many cracks is also solved. For homogeneous materials it is found that the normal electric displacement D-2 induced by the crack is constant along the crack faces which depends only on the applied remote stress field. Within the crack slit, the electric fields induced by the crack are also constant and not affected by the applied electric field. For the bimaterials with real H, the normal electric displacement D-2 is constant along the crack faces and electric field E-2 has the singularity ahead of the crack tip and a jump across the interface.

Analysis of electromagnetic propagation in birefringent thin film

In general, the propagating behavior of extraordinary wave in anisotropic materials is different from that in isotropic materials. With the tangential continuity of Maxwell's equations, the electromagnetic propagating behaviors have been investigated at the incident and exit interfaces of the uniaxial anisotropic thin film. The emphasis was placed on two interesting optical phenomena such as homolateral refraction behavior and wide-angle Brewster's phenomenon, which occurred at the interfaces of uniaxial anisotropic thin film.

Y-chromosome evidence for no independent origin of modern human in China

East Asia is one of the few regions in the world where a large number of human fossils have been unearthed. The continuity of hominid fossils in East Asia, particularly in China has been presented as strong evidence supporting an independent origin of mod

Genetic Stability of Progeny from an Artificial Allotetraploid Carp Using Sperm from Five Fish Species

Microsatellite markers and D-loop sequences of mtDNA from a female allotetraploid parent carp and her progenies of generations 1 and 2 induced by sperm of five distant fish species were analyzed. Eleven microsatellite markers were used to identify 48 alleles from the allotetraploid female. The same number of alleles (48) appeared in the first and second generations of the gynogenetic offspring, regardless of the source of the sperm used as an activator. The mtDNA D-loop analysis was performed on the female tetraploid parent, 25 gynogenetic offspring, and 5 sperm-donor species. Fourteen variable sites from the 1,018 bp sequences were observed in the offspring as compared to the female tetraploid parent. Results from D-loop sequence and microsatellite marker analysis showed exclusive maternal transmission, and no genetic information was derived from the father. Our study suggests that progenies of artificial tetraploid carp are genetically stable, which is important for genetic breeding of this tetraploid fish.

Constructing Thin Polythiophene Film Composed of Aligned Lamellae via Controlled Solvent Vapor Treatment

Thin poly(3-butylthiophene) (P3BT) film composed of aligned lamellae attached to the edge of the original film has been achieved via a controlled solvent vapor treatment (C-SVT) method. The polarized optical microscopy operated at both single-polarization and cross-polarization modes has been used to investigate the alignment of the fiber-like lamellae. A numerical simulation method is used to quantitatively calculate angle distributions of the lamellae deviated from the film growth direction. Prepatterned P3BT film edge acts as nuclei which densely initialize subsequent crystal growth by exhausting the materials transported from the partially dissolved film. The growth of new film upon crystallization is actually a self-healing process where the two-dimensional geometric confinement is mainly responsible for this parallel alignment of P3BT crystals. The solvent vapor pressure should be carefully chosen so as to induce crystal growth but avoid liquid instability which will destroy the continuity of the film.

The effects of shell characteristics on the current-voltage behaviors of dye-sensitized solar cells based on ZnO/TiO2 core/shell arrays

The ZnO/TiO2 core/shell structure was formed through deposition of a TiO2 coating layer on the hydrothermally fabricated ZnO nanorod arrays through radio frequency magnetron sputtering. The effects of the TiO2 shell's characteristics on the current-voltage behaviors of the core/shell-based dye-sensitized solar cells (CS-DSSC) were investigated. As the rates of injection, transfer, and recombination of electrons of such CS-DSSC were affected significantly by the crystallization, morphology, and continuity of the TiO2 shells, the photovoltaic efficiency was accordingly varied remarkably. In addition, the efficiency was further improved by enhancing the surface area in the core/shell electrode.

A Constrained Particle Dynamics for Continuum-Particle Hybrid Method in Micro- and Nano-Fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier-Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.

Co和Cu薄膜电学特性和连续性特征的判据

采用离子束溅射沉积了不同厚度的Co膜和Cu膜，利用四电极法测量了薄膜的电阻率，从而得到了Co膜和Cu膜的电导率随薄膜厚度的变化关系。实验结果表明，Co膜和Cu膜的电学特性都具有明显的尺寸效应。比较了同时考虑表面散射和晶界散射的电导理论得到的电导率公式与实验结果，不同薄膜厚度电导率的理论结果与实验结果符合较好。提出了厚度作为金属薄膜生长从不连续膜进入连续膜的一个特征判据，并利用原子力显微镜（AFM）观测了膜厚在特征厚度附近的Co膜和Cu膜的表面形貌。

Rashba electron transport in one-dimensional quantum waveguides

The properties of Rashba wave function in the planar one-dimensional waveguide are studied, and the following results are obtained. Due to the Rashba effect, the plane waves of electron with the energy E divide into two kinds of waves with the wave vectors k(1)=k(0)+k(delta) and k(2)=k(0)-k(delta), where k(delta) is proportional to the Rashba coefficient, and their spin orientations are +pi/2 (spin up) and -pi/2 (spin down) with respect to the circuit, respectively. If there is gate or ferromagnetic contact in the circuit, the Rashba wave function becomes standing wave form exp(+/- ik(delta)l)sin[k(0)(l-L)], where L is the position coordinate of the gate or contact. Unlike the electron without considering the spin, the phase of the Rashba plane or standing wave function depends on the direction angle theta of the circuit. The travel velocity of the Rashba waves with the wave vector k(1) or k(2) are the same hk(0)/m*. The boundary conditions of the Rashba wave functions at the intersection of circuits are given from the continuity of wave functions and the conservation of current density. Using the boundary conditions of Rashba wave functions we study the transmission and reflection probabilities of Rashba electron moving in several structures, and find the interference effects of the two Rashba waves with different wave vectors caused by ferromagnetic contact or the gate. Lastly we derive the general theory of multiple branches structure. The theory can be used to design various spin polarized devices.

紫菜壳孢子萌发过程及其世代差异研究

紫菜是一类海洋藻类研究的模式系统，具有重要的经济价值和理论研究意义。本研究基于紫菜的世代交替生活史，对其壳孢子萌发过程进行了研究，并对丝状孢子体与叶状配子体的世代差异进行了分析，主要包括以下三部分： 1） 对紫菜丝状孢子体与叶状配子体的连接枢纽——壳孢子的发育过程及其光影响进行了研究与讨论，发现壳孢子只有在附着后才能形成细胞壁并发育，说明壳孢子的附着是触发了壳孢子细胞壁的形成以及后期发育的开关，亦即附着是触发紫菜世代交替过程中配子体转录组表达的开关；纤维素酶和果胶酶均能抑制壳孢子附着，但影响机制各不相同，推测果胶质主要介导壳孢子的初始附着，而纤维素则与永久附着相关；波长≥580 nm的高强度（200 μmol•m-2•s-1）可见光有利于壳孢子早期发育。 2） 结合现有藻类数据，对坛紫菜丝状孢子体阶段11000 EST数据进行了大规模的生物信息学分析，结果首次发现坛紫菜丝状孢子体中可能存在PCK型C4光合固碳途径，并筛选出44条在紫菜孢子体中表达上调的代表基因。 3） 结合坛紫菜、条斑紫菜、海带和红毛菜，对红藻和褐藻等大型海藻孢子体与配子体阶段代表基因Rubisco的表达与羧化酶活性差异进行了研究分析，结果表明Rubisco的表达量和初始羧化酶活在其配子体中均显著高于其孢子体世代，即与藻体不同世代的相对复杂度无关，而与染色体倍性相关，说明Rubisco的世代差异极有可能与染色体倍性相连锁，因而可能是海藻世代交替过程中的重要功能基因。

砂泥岩薄互层储层预测方法与应用研究

With the deeply development of exploration and development in petroleum in China, new increasing reserves are found in old oil fields and the verge of the old ones through re-study of geological property. It is more and more important to discovery and develop thin layer or thin inter-bedded layers reservoirs. All of the targets are thin sand-shale inter-bedded reservoirs and the core technology is reservoir predictions between wells in thin sand-shale inter-bedded layers. The continuity of the thin sand-shale inter-bedded layers in space or separating and heterogeneity is the key of reservoir geology research. The seismic reflection, high resolution analysis method and inversion method to thin sand-shale inter-bedded layers are thorough discussed and deeply studied in this paper to try to find the methods and resolutions of reservoir geology research. The below is followed. 1. Based on the pre-research of other people, five models are created: the sand sphenoid body, interlay sandstone and interlay shale of the equal thickness, interlay sandstone of the equal thickness and interlay shale of the unequal thickness, interlay sandstone of the unequal thickness and interlay shale of the unequal thickness, interlay sandstone of the changing thickness in sequence and interlay shale of the changing thickness in sequence. Then the study of the forward modeling are conducted on the thin layer and thin inter-bedded layers geological characters and seismic reflections including amplitude, frequency, phase, wave shape and time-frequency responding in the domains of time and frequency. The affect of petro-physics difference of layers, single thin layer thickness, thickness of inter-bedded, layer number of inter-bedded, incident wavelet domain frequency and types, sample interval to seismic reflection characters, frequency spectrum and time-frequency respond of reflectivity is theoretically discussed. 2. Qualitatively analyzing the sedimentary rhythm of the thin inter-bedded layers in vertical orientation and computing the single layer thickness or the average thickness with the method of generalized S transform. Identifying the reflecting interface or lithology interface using the amplitude value of amplitude spectrum domain frequency. 3. Based on the seismic respond of thin sand-shale inter-bedded layers, bring out the high resolution analysis method of seismic data in thin sand-shale inter-bedded layers using wavelet analysis and the idea of affecting low and high frequency with middle frequency. Then analyzing the effect to the method and testing some wavelets in the method. This method is applied to the theoretical models and the field data. 4. Bring forward one improved very fast simulated annealing method (IVFSA) to resolve the problem nonlinearity and multi-parameters of the inversion in thin inter-bedded layers. And IVFSA is more productive and higher precision than general ways. 5. New target constrained function is used in the inversion based on the property of the inversion in thin inter-bedded layers. 6. Making the full use of geological and logging information, IVFSA and the new function are applied in the non-linear inversion to improve reservoir prediction and evaluation in thin inter-bedded formations combined with the idea of logging and seismic inversion. This method was applied to the field data and got good results.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

急倾斜金属矿体开采岩体移动规律与变形趋势预测研究

Rockmass movement due to mining steep metallic ore body is a considerable question in the surface movement and deformation issue caused by underground mining. Research on coal mining induced rockmass movement and its prediction problem have been performed for a long-term, and have achieved great progress at home and abroad. However, the rockmass movement caused by mining steep metal mine is distinctivly different from coal seam mining.. Existing surface movement laws and deformation prediction methods are not applicable to the rockmass movement caused by mining steep metal mine. So far the home and abroad research to this theory is presently at an early stage, and there isn’t mature theory or practical prediction method, which made a great impact on production. In this paper, the research object—Jinchuan nickel mine, which is typical steep metal mine, characterized by complex geological conditions, developed faults, cracked rockmass, high geostress, and prominent engineering stability problems. In addition, backfill mining method is used in the mine, the features of rockmass movement caused by this mining method are also different from other mining methods. In this paper, the laws of rock mass movement, deformation and destroy mechanism, and its prediction were analyzed based on the collection of data, detailed in-sit engineering geology survey, ground movement monitoring by GPS, theoretical analysis and numerical simulation. According to the GPS monitoring of ground surface movement, ground subsidence basin with apparent asymmetry is developing, the influence scope is larger in the upper faulted block than in the lower faulted block, and the center of ground movement is moving along the upper faulted block direction with increasing depth of mining. During the past half and seven years, the largest settlement has amounted to 1287.5mm, and corresponding horizontal displacement has amounted to 664.6mm. On the ground surface, two fissure belts show a fast-growing trend of closure. To sum up, mining steep metal mine with backfill method also exist the same serious problem of rockmass movement hazards. Fault, as a low intensity zone in rockmass, when it located within the region of mining influence, the change of potential energy mainly consumed in fault deformation associated with rockmass structure surface friction, which is the essence of displacement and stress barrier effects characterized by fault rupture zone. when steep fault located in the tensile deformation region incurred by underground excavation, no matter excavation in hangingwall or in footwall of the fault, there will be additional tensile stress on the vertical fault plane and decrease in the shear strength, and always showing characteristics of normal fault slip, which is the main reason of fault escarpment appeared on the ground surface. The No.14 shaft deformation and failure is triggered by fault activation, which showed with sidewall move, rupture, and break down features as the main form of a concentrated expression of fault effects. The size and orientation of principal stress in surrounding rock changed regularly with mining; therefore, roadway deformation and damage at different stages have different characteristics and distribution models. During the process of mining, low-intensity weak structures surface always showed the most obvious reaction, accompany with surface normal stress decrease and shear strength bring down, to some extent, occurred with relative slide and deformation. Meanwhile, the impact of mining is a relatively long process, making the structure surface effect of roadway deformation and damage more prominent than others under the influence of mining. Roadway surrounding rockmass deformation caused by the change of strain energy density field after excavation mainly belongs to elastic deformation, and the correspondented damage mainly belongs to brittle rupture, in this circumstance, surrounding rockmass will not appear large deformation. The large deformation of surrounding rockmass can only be the deformation associated with structure surface friction or the plastic deformation of itself, which mainly caused by the permanent self-weigh volume force，and long-term effect of mining led to the durability of this deformation Good pitting fill effect and supporting effect of backfill, as well as the friction of rockmass structure surface lead to obvious macro-rockmass movement with long-lag characteristics. In addition, the loss of original intensity and new structure surface arisen increased flexibility in rockmass and fill deformation in structure surface, which made the time required for rockmass potential energy translate into deformation work associated with plastic deformation and structure surface friction consumed much, and to a large extent, eliminated the time needed to do those plastic work during repeated mining, all of which are the fundamental reason of rockmass movement aftereffect more significant than before. Mining steep deposits in high tectonic stress area and in gravity stress area have different movement laws and deformation mechanism. The steep deposit, when the vertical size of the mining areas is smaller than the horizontal size of the orebody, no matter mining in gravity stress area or in high tectonic stress area, they have similar features of ground movement with mining horizontal orebody; contrarily, there will appear double settlement centers on the ground surface under the condition of mining in high tectonic stress area, while there will always be a single center under the other condition. Meanwhile the ground movement lever, scale of mining influence area and macro features of ground movement, deformation and fracture are also different from mining in gravity stress area, and the fundamental reason lies in the impact of orientation of the maximum principal stress on rock movement features in in-site rock stress field. When mining thick and steep deposit, the ground surface movement and deformation characteristic curves are significantly different from excavating the horizontal ore bed and thin steep deposit. According to the features of rockmass movement rate, the development process of mining-induced rockmass movement is divided into three stages: raising stage, steadily stage and gradually decay stage. Considering the actual exploitation situation, GPS monitoring results and macro-characteristics of surface movement, the current subsidence pattern of Jinchuan No.2 mine is in the early stage of development. Based on analysis of surface movement rate, surface subsidence rate increase rapidly when mining in double lever at the same time, and reach its peak until the exploitation model ended. When double lever mining translate into single, production decreased, surface subsidence rate suddenly start to reduce and maintain a relatively low value, and the largest subsidence center will slowly move along with the hangingwall ore body direction with increasing depth of mining, at the same time, the scope and extent of subsidence in footwall ore body will begin magnify, and a sub-settlement center will appear on ground surface, accompanied with the development and closure trend of ground fissure, the surrounding rockmass of shaft and roadway will be confronted to more frequent and severe deformation and failure, and which will have a negative impact on the overall stability of No.2 mine mining. On the premise of continuity of rockmass movement, gray system model can be used in ground rockmass movement prediction for good results. Under the condition of backfill mining step by step, the loose effect of compact status of the hard, broken rockmass led to lower energy release rate, although surrounding rockmass has high elastic energy, loose and damage occurred in the horizontal ore body, which made the mining process safety without any large geological hazards. During the period of mining the horizontal ore body to end, in view of its special “residual support role”, there will be no large scale rockmass movement hazards. Since ground surface movement mainly related to the intensity of mining speed and backfill effect, on the premise of constant mining speed, during the period of mining the horizontal ore body to end, the rate of ground surface rockmass movement and deformation won’t have sudden change.