Lower-Bound Axisymmetric Formulation for Geomechanics Problems Using Nonlinear Optimization

A lower-bound limit analysis formulation, by using two-dimensional finite elements, the three-dimensional Mohr-Coulomb yield criterion, and nonlinear optimization, has been given to deal with an axisymmetric geomechanics stability problem. The optimization was performed using an interior point method based on the logarithmic barrier function. The yield surface was smoothened (1) by removing the tip singularity at the apex of the pyramid in the meridian plane and (2) by eliminating the stress discontinuities at the corners of the yield hexagon in the pi-plane. The circumferential stress (sigma(theta)) need not be assumed. With the proposed methodology, for a circular footing, the bearing-capacity factors N-c, N-q, and N-gamma for different values of phi have been computed. For phi = 0, the variation of N-c with changes in the factor m, which accounts for a linear increase of cohesion with depth, has been evaluated. Failure patterns for a few cases have also been drawn. The results from the formulation provide a good match with the solutions available from the literature. (C) 2014 American Society of Civil Engineers.

Electronic Structure Evolution across the Peierls Metal-Insulator Transition in a Correlated Ferromagnet

Transition metal compounds often undergo spin-charge-orbital ordering due to strong electron-electron correlations. In contrast, low-dimensional materials can exhibit a Peierls transition arising from low-energy electron-phonon-coupling-induced structural instabilities. We study the electronic structure of the tunnel framework compound K2Cr8O16, which exhibits a temperature-dependent (T-dependent) paramagnetic-to-ferromagnetic- metal transition at T-C = 180 K and transforms into a ferromagnetic insulator below T-MI = 95 K. We observe clear T-dependent dynamic valence (charge) fluctuations from above T-C to T-MI, which effectively get pinned to an average nominal valence of Cr+3.75 (Cr4+:Cr3+ states in a 3:1 ratio) in the ferromagnetic-insulating phase. High-resolution laser photoemission shows a T-dependent BCS-type energy gap, with 2G(0) similar to 3.5(k(B)T(MI)) similar to 35 meV. First-principles band-structure calculations, using the experimentally estimated on-site Coulomb energy of U similar to 4 eV, establish the necessity of strong correlations and finite structural distortions for driving the metal-insulator transition. In spite of the strong correlations, the nonintegral occupancy (2.25 d-electrons/Cr) and the half-metallic ferromagnetism in the t(2g) up-spin band favor a low-energy Peierls metal-insulator transition.

Observation of correlated spin-orbit order in a strongly anisotropic quantum wire system

Quantum wires with spin-orbit coupling provide a unique opportunity to simultaneously control the coupling strength and the screened Coulomb interactions where new exotic phases of matter can be explored. Here we report on the observation of an exotic spin-orbit density wave in Pb-atomic wires on Si(557) surfaces by mapping out the evolution of the modulated spin-texture at various conditions with spin-and angle-resolved photoelectron spectroscopy. The results are independently quantified by surface transport measurements. The spin polarization, coherence length, spin dephasing rate and the associated quasiparticle gap decrease simultaneously as the screened Coulomb interaction decreases with increasing excess coverage, providing a new mechanism for generating and manipulating a spin-orbit entanglement effect via electronic interaction. Despite clear evidence of spontaneous spin-rotation symmetry breaking and modulation of spin-momentum structure as a function of excess coverage, the average spin polarization over the Brillouin zone vanishes, indicating that time-reversal symmetry is intact as theoretically predicted.

Imprint of modified Einstein's gravity on white dwarfs: Unifying Type Ia supernovae

We establish the importance of modified Einstein's gravity (MG) in white dwarfs (WDs) for the first time in the literature. We show that MG leads to significantly sub- and super-Chandrasekhar limiting mass WDs, depending on a single model parameter. However, conventional WDs on approaching Chandrasekhar's limit are expected to trigger Type Ia supernovae (SNeIa), a key to unravel the evolutionary history of the universe. Nevertheless, observations of several peculiar, under-and over-luminous SNeIa argue for the limiting mass widely different from Chandrasekhar's limit. Explosions of MG induced sub-and super-Chandrasekhar limiting mass WDs explain under-and over-luminous SNeIa respectively, thus unifying these two apparently disjoint sub-classes. Our discovery questions both the global validity of Einstein's gravity and the uniqueness of Chandrasekhar's limit.

Probing Photoexcited Carriers in a Few-Layer MoS2 Laminate by Time-Resolved Optical Pump-Terahertz Probe Spectroscopy

We report the dynamics of photoinduced carriers in a free-standing MoS2 laminate consisting of a few layers (1-6 layers) using time-resolved optical pump-terahertz probe spectroscopy. Upon photoexcitation with the 800 nm pump pulse, the terahertz conductivity increases due to absorption by the photoinduced charge carriers. The relaxation of the non-equilibrium carriers shows fast as well as slow decay channels, analyzed using a rate equation model incorporating defect-assisted Auger scattering of photoexcited electrons, holes, and excitons. The fast relaxation time occurs due to the capture of electrons and holes by defects via Auger processes, resulting in nonradiative recombination. The slower relaxation arises since the excitons are bound to the defects, preventing the defect-assisted Auger recombination of the electrons and the holes. Our results provide a comprehensive understanding of the non-equilibrium carrier kinetics in a system of unscreened Coulomb interactions, where defect-assisted Auger processes dominate and should be applicable to other 2D systems.

Behavior of Brick-Mortar Interfaces in FRP-Strengthened Masonry Assemblages under Normal Loading and Shear Loading

Determination of shear strength of brick-mortar bed joint is critical to overcome the sliding-shear or joint-shear failure in masonry. In the recent past, researchers have attempted to enhance the shear strength and deformation capacity of brick-mortar bed joints by gluing fiber-reinforced polymer (FRP) composite across the bed joint. FRP composites offer several advantages like high strength-to-weight ratio, and ease of application in terms of labor, time, and reduced curing period. Furthermore, FRP composites are desirable for strengthening old masonry buildings having heritage value because of its minimal interference with the existing architecture. A majority of earlier studies on shear strengthening of masonry available in the literature adopted masonry having the ratio of modulus of elasticity of masonry unit (Emu) to modulus of elasticity of mortar (Em) greater than one. Information related to shear behavior of FRP glued masonry composed of masonry units having Young's modulus lower than mortar is limited. Hence the present study is focused on characterizing the interfacial behavior of brick-mortar bed joint of masonry assemblages composed of solid burnt clay bricks and cement-sand mortar (E-mu/E-m ratio less than one), strengthened with FRP composites. Masonry triplets and prisms with bed joint inclined to loading axis (0 degrees, 30 degrees, 45 degrees, 60 degrees and 90 degrees) are employed in this study. Glass and carbon FRP composites composed of bidirectional FRP fabric with equal density in both directions are used for strengthening masonry. Masonry triplets are glued with glass and carbon FRP composites in two configurations: (1) both faces of the triplet specimens are fully glued with GFRP composites; and (2) both faces of the triplet specimens are glued with GFRP and CFRP composites in strip form. The performance of masonry assemblages strengthened with FRP composites is assessed in terms of gain in shear strength, shear displacement, and postpeak behavior for various configurations and types of FRP composites considered. A semianalytical model is proposed for the prediction of shear strength of masonry bed joints glued with FRP composites. A composite failure envelope consisting of a Coulomb friction model and a compression cap is obtained for unreinforced masonry and GFRP-strengthened masonry based on the test results of masonry triplets and masonry prisms with bed joints having various inclinations to the loading (C) 2015 American Society of Civil Engineers.

Uplift Resistance of Interfering Pipelines Buried in Sand

The vertical uplift resistance of interfering pipelines buried in sands has been computed using the lower-bound limit analysis in conjunction with finite elements and nonlinear optimization. The soil mass is assumed to follow the Mohr-Coulomb failure criterion and an associated flow rule. It is specified that all the pipes fail simultaneously at the same magnitude of the failure load. For different clear spacing (S) between the pipes, the magnitude of the efficiency factor (xi(gamma)) is determined. Because of pipes' interference, with a reduction in the spacing between the pipelines, the magnitude of xi(gamma) is found to decrease continuously. The results were found to compare quite well with the available data from literature for horizontal strip anchors. (C) 2015 American Society of Civil Engineers.

Increasing Critical Sensitivity of the Load/Unload Response Ratio before Large Earthquakes with Identified Stress Acclumulation Pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.

Comparison of the Dynamic Response of a Saturated Soil Between Two Models

By comparing the dynamic responses of saturated soil to Biot's and Yamamoto's models, the properties of the two models have be pointed out. First of all, an analysis has been made for energy loss of each model from the basic equations. Then the damping of elastic waves in coarse sand and fine sand with loading frequency and soil's parameters have been calculated and the representation of viscous friction and Coulomb friction in the two models has been concluded. Finally, the variations of loading wave damping and stress phase angles with water depth and soil's parameters have been obtained as loading waves range in ocean waves.

Comparison of existing poroelastic models for wave damping in a porous seabed

Mechanism of wave-seabed interaction has been extensively studied by coastal geotechnical engineers in recent years. Numerous poro-elastic models have been proposed to investigate the mechanism of wave propagation on a seabed in the past. The existing poro-elastic models include drained model, consolidation model, Coulomb-damping model, and full dynamic model. However, to date, the difference between the existing models is unclear. In this paper, the fully dynamic poro-elastic model for the wave-seabed interaction will be derived first. Then, the existing models will be reduced from the proposed fully dynamic model. Based on the numerical comparisons, the applicable range of each model is also clarified for the engineering practice.

Modification of Bertrand's theorem and extended Runge-Lenz vector

It is shown that for a particle with suitable angular moments in the screened Coulomb potential or isotropic harmonic potential, there still exist closed orbits rather than ellipse, characterized by the conserved aphelion and perihelion vectors, i.e. extended Runge-Lenz vector, which implies a higher dynamical symmetry than the geometrical symmetry O-3. The closeness of a planar orbit implies the radial and angular motional frequencies are commensurable.

经典轨道的封闭性和径向Schrodinger方程的因式分解

研究表明，保证经典轨道具有封闭性的Bertrand定理可以进一步推广，在适当的角动量睛，仍存在着非椭圆的闭合轨道。对于屏蔽Coulomb场，可获得广义Runge-Lenz矢量。这种轨道封闭性与径向Schrodinger方程因式分解相对应。

一种修正的莫尔-库仑本构模型及其应用

莫尔-库仑本构模型是理想弹塑性的，只考虑了强度随围压的变化，但没有考虑弹性模量随围压和应力水平的变化情况，这不符合地质材料的实际特性，借助于邓肯-张非线性弹性模型对切线模量的处理方法，对M-C模型中的弹性模量进行了修正。修正以后，弹性阶段的应力-应变曲线将表现为双曲线。为了将该修正模型写进FLAC^3D程序，特编制了前处理程序，并用来对刚性桩的横向承载力进行了数值分析。因为该模型没有对卸载模量重新定义，所以该模型较适合于持续加载的情况。

波浪、海洋土参数对海床稳定性影响

基于Yamamoto的多孔弹性介质模型，研究了波生底床的稳定性。通过给出的有限深底床下土响应分析解，针对三种土质底床，讨论了主要波参数和土参数对这些底床稳定性的影响。与其他土模型计算结果进行了比较，分析了海洋土内部Coulomb摩擦因素的影响。

爆炸排淤填石法中岩土介质的本构模型

以爆炸排淤填石法为背景，对相关岩土介质的本构模型进行了探讨。分析认为，在以LS.DYNA 动态有限元分析程序对爆炸排淤填石法的数值模拟中，岩土介质不适于采用Mohr.Coulomb 模型，堆石体宜采用Drucker-Prager 模型，淤泥宜采用Prandtl-Reuss 模型；利用LS-DYNA 程序对淤泥的本构模型进行了验证和确认，计算表明：在形成爆炸空腔的高应变率阶段，淤泥表现为理想不可压缩流体的性质；在小药量小抵抗线条件下，在淤泥自重作用下的低应变率变形阶段，其粘性效应可以忽略不计。