Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis

The ultimate bearing capacity of a circular footing, placed over rock mass, is evaluated by using the lower bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The generalized Hoek-Brown (HB) failure criterion, but by keeping a constant value of the exponent, alpha = 0.5, was used. The failure criterion was smoothened both in the meridian and pi planes. The nonlinear optimization was carried out by employing an interior point method based on the logarithmic barrier function. The results for the obtained bearing capacity were presented in a non-dimensional form for different values of GSI, m(i), sigma(ci)/(gamma b) and q/sigma(ci). Failure patterns were also examined for a few cases. For validating the results, computations were also performed for a strip footing as well. The results obtained from the analysis compare well with the data reported in literature. Since the equilibrium conditions are precisely satisfied only at the centroids of the elements, not everywhere in the domain, the obtained lower bound solution will be approximate not true. (C) 2015 Elsevier Ltd. All rights reserved.

A FRAMEWORK FOR THE ERROR ANALYSIS OF DISCONTINUOUS FINITE ELEMENT METHODS FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS AND APPLICATIONS TO C-0 IP METHODS

In this article, an abstract framework for the error analysis of discontinuous Galerkin methods for control constrained optimal control problems is developed. The analysis establishes the best approximation result from a priori analysis point of view and delivers a reliable and efficient a posteriori error estimator. The results are applicable to a variety of problems just under the minimal regularity possessed by the well-posedness of the problem. Subsequently, the applications of C-0 interior penalty methods for a boundary control problem as well as a distributed control problem governed by the biharmonic equation subject to simply supported boundary conditions are discussed through the abstract analysis. Numerical experiments illustrate the theoretical findings.

Bearing capacity factors for a conical footing using lower- and upper-bound finite elements limit analysis

Bearing capacity factors, N-c, N-q, and N-gamma, for a conical footing are determined by using the lower and upper bound axisymmetric formulation of the limit analysis in combination with finite elements and optimization. These factors are obtained in a bound form for a wide range of the values of cone apex angle (beta) and phi with delta = 0, 0.5 phi, and phi. The bearing capacity factors for a perfectly rough (delta = phi) conical footing generally increase with a decrease in beta. On the contrary, for delta = 0 degrees, the factors N-c and N-q reduce gradually with a decrease in beta. For delta = 0 degrees, the factor N-gamma for phi >= 35 degrees becomes a minimum for beta approximate to 90 degrees. For delta = 0 degrees, N-gamma for phi <= 30 degrees, as in the case of delta = phi, generally reduces with an increase in beta. The failure and nodal velocity patterns are also examined. The results compare well with different numerical solutions and centrifuge tests' data available from the literature.

Power Allocation in MIMO Wiretap Channel with Statistical CSI and Finite-Alphabet Input

In this paper, we consider the problem of power allocation in MIMO wiretap channel for secrecy in the presence of multiple eavesdroppers. Perfect knowledge of the destination channel state information (CSI) and only the statistical knowledge of the eavesdroppers CSI are assumed. We first consider the MIMO wiretap channel with Gaussian input. Using Jensen's inequality, we transform the secrecy rate max-min optimization problem to a single maximization problem. We use generalized singular value decomposition and transform the problem to a concave maximization problem which maximizes the sum secrecy rate of scalar wiretap channels subject to linear constraints on the transmit covariance matrix. We then consider the MIMO wiretap channel with finite-alphabet input. We show that the transmit covariance matrix obtained for the case of Gaussian input, when used in the MIMO wiretap channel with finite-alphabet input, can lead to zero secrecy rate at high transmit powers. We then propose a power allocation scheme with an additional power constraint which alleviates this secrecy rate loss problem, and gives non-zero secrecy rates at high transmit powers.

Intrinsic large gap quantum anomalous Hall insulators in LaX (X = Br,Cl,I)

We report a theoretical prediction of a new class of bulk and intrinsic quantum anomalous Hall (QAH) insulators LaX (X=Br, Cl, and I) via relativistic first-principles calculations. We find that these systems are innate long-ranged ferromagnets which, with the help of intrinsic spin-orbit coupling, become QAH insulators. A low-energy multiband tight-binding model is developed to understand the origin of the QAH effect. Finally, integer Chern number is obtained via Berry phase computation for each two-dimensional plane. These materials have the added benefit of a sizable band gap of as large as similar to 25 meV, with the flexibility of enhancing it to above 75 meV via strain engineering. The synthesis of LaX materials will provide the impurity-free single crystals and thin-film QAH insulators for versatile experiments and functionalities.

A Study of Pressure-Dependent Squeeze Film Stiffness as a Resonance Modulator Using Static and Dynamic Measurements

We report on the resonant frequency modulation of inertial microelectromechanical systems (MEMS) structures due to squeeze film stiffness over a range of working pressures. Squeeze film effects have been studied extensively, but mostly in the context of damping and Q-factor determination of dynamic MEMS structures, typically suspended over a fixed substrate with a very thin air gap. Here, we show with experimental measurements and analytical calculations how the pressure-dependent air springs (squeeze film stiffness) change the resonant frequency of an inertial MEMS structure by as much as five times. For capturing the isolated effect of the squeeze film stiffness, we first determine the static stiffness of our structure with atomic force microscope probing and then study the effect of the air spring by measuring the dynamic response of the structure, thus finding the resonant frequencies while varying the air pressure from 1 to 905 mbar. We also verify our results by analytical and Finite Element Method calculations. Our findings show that the pressure-dependent squeeze film stiffness can affect a rather huge range of frequency modulation (>400%) and, therefore, can be used as a design parameter for exploiting this effect in MEMS devices. 2014-0310]

Analytic bootstrap at large spin

We use analytic conformal bootstrap methods to determine the anomalous dimensions and OPE coefficients for large spin operators in general conformal field theories in four dimensions containing a scalar operator of conformal dimension Delta(phi). It is known that such theories will contain an in finite sequence of large spin operators with twists approaching 2 Delta(phi) + 2n for each integer n. By considering the case where such operators are separated by a twist gap from other operators at large spin, we analytically determine the n, Delta(phi) dependence of the anomalous dimensions. We find that for all n, the anomalous dimensions are negative for Delta(phi) satisfying the unitarity bound. We further compute the first subleading correction at large spin and show that it becomes universal for large twist. In the limit when n is large, we find exact agreement with the AdS/CFT prediction corresponding to the Eikonal limit of a 2-2 scattering with dominant graviton exchange.

Efficient density matrix renormalization group algorithm to study Y junctions with integer and half-integer spin

An efficient density matrix renormalization group (DMRG) algorithm is presented and applied to Y junctions, systems with three arms of n sites that meet at a central site. The accuracy is comparable to DMRG of chains. As in chains, new sites are always bonded to the most recently added sites and the superblock Hamiltonian contains only new or once renormalized operators. Junctions of up to N = 3n + 1 approximate to 500 sites are studied with antiferromagnetic (AF) Heisenberg exchange J between nearest-neighbor spins S or electron transfer t between nearest neighbors in half-filled Hubbard models. Exchange or electron transfer is exclusively between sites in two sublattices with N-A not equal N-B. The ground state (GS) and spin densities rho(r) = < S-r(z)> at site r are quite different for junctions with S = 1/2, 1, 3/2, and 2. The GS has finite total spin S-G = 2S(S) for even (odd) N and for M-G = S-G in the S-G spin manifold, rho(r) > 0(< 0) at sites of the larger (smaller) sublattice. S = 1/2 junctions have delocalized states and decreasing spin densities with increasing N. S = 1 junctions have four localized S-z = 1/2 states at the end of each arm and centered on the junction, consistent with localized states in S = 1 chains with finite Haldane gap. The GS of S = 3/2 or 2 junctions of up to 500 spins is a spin density wave with increased amplitude at the ends of arms or near the junction. Quantum fluctuations completely suppress AF order in S = 1/2 or 1 junctions, as well as in half-filled Hubbard junctions, but reduce rather than suppress AF order in S = 3/2 or 2 junctions.

SYNGAP1: Mind the Gap

A cardinal feature of early stages of human brain development centers on the sensory, cognitive, and emotional experiences that shape neuronal-circuit formation and refinement. Consequently, alterations in these processes account for many psychiatric and neurodevelopmental disorders. Neurodevelopment disorders affect 3-4% of the world population. The impact of these disorders presents a major challenge to clinicians, geneticists, and neuroscientists. Mutations that cause neurodevelopmental disorders are commonly found in genes encoding proteins that regulate synaptic function. Investigation of the underlying mechanisms using gain or loss of function approaches has revealed alterations in dendritic spine structure, function, and plasticity, consequently modulating the neuronal circuit formation and thereby raising the possibility of neurodevelopmental disorders resulting from synaptopathies. One such gene, SYNGAP1 (Synaptic Ras-GTPase-activating protein) has been shown to cause Intellectual Disability (ID) with comorbid Autism Spectrum Disorder (ASD) and epilepsy in children. SYNGAP1 is a negative regulator of Ras, Rap and of AMPA receptor trafficking to the postsynaptic membrane, thereby regulating not only synaptic plasticity, but also neuronal homeostasis. Recent studies on the neurophysiology of SYNGAP1, using Syngapl mouse models, have provided deeper insights into how downstream signaling proteins and synaptic plasticity are regulated by SYNGAP1. This knowledge has led to a better understanding of the function of SYNGAP1 and suggests a potential target during critical period of development when the brain is more susceptible to therapeutic intervention.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

Reply to Discussion on ``Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis'' Comput. Geotech. 70 (2015) 138-149]

A discussion has been provided for the comments raised by the discusser (Clausen, 2015)1] on the article recently published by the authors (Chakraborty and Kumar, 2015). The effect of exponent alpha for values of GSI approximately smaller than 30 becomes more critical. On the other hand, for greater values of GSI, the results obtained by the authors earlier remain primarily independent of alpha and can be easily used. (C) 2015 Elsevier Ltd. All rights reserved.

SYNGAP1: Mind the Gap

A cardinal feature of early stages of human brain development centers on the sensory, cognitive, and emotional experiences that shape neuronal-circuit formation and refinement. Consequently, alterations in these processes account for many psychiatric and neurodevelopmental disorders. Neurodevelopment disorders affect 3-4% of the world population. The impact of these disorders presents a major challenge to clinicians, geneticists, and neuroscientists. Mutations that cause neurodevelopmental disorders are commonly found in genes encoding proteins that regulate synaptic function. Investigation of the underlying mechanisms using gain or loss of function approaches has revealed alterations in dendritic spine structure, function, and plasticity, consequently modulating the neuronal circuit formation and thereby raising the possibility of neurodevelopmental disorders resulting from synaptopathies. One such gene, SYNGAP1 (Synaptic Ras-GTPase-activating protein) has been shown to cause Intellectual Disability (ID) with comorbid Autism Spectrum Disorder (ASD) and epilepsy in children. SYNGAP1 is a negative regulator of Ras, Rap and of AMPA receptor trafficking to the postsynaptic membrane, thereby regulating not only synaptic plasticity, but also neuronal homeostasis. Recent studies on the neurophysiology of SYNGAP1, using Syngapl mouse models, have provided deeper insights into how downstream signaling proteins and synaptic plasticity are regulated by SYNGAP1. This knowledge has led to a better understanding of the function of SYNGAP1 and suggests a potential target during critical period of development when the brain is more susceptible to therapeutic intervention.

Ellipse Fitting Using the Finite Rate of Innovation Sampling Principle

Standard approaches for ellipse fitting are based on the minimization of algebraic or geometric distance between the given data and a template ellipse. When the data are noisy and come from a partial ellipse, the state-of-the-art methods tend to produce biased ellipses. We rely on the sampling structure of the underlying signal and show that the x- and y-coordinate functions of an ellipse are finite-rate-of-innovation (FRI) signals, and that their parameters are estimable from partial data. We consider both uniform and nonuniform sampling scenarios in the presence of noise and show that the data can be modeled as a sum of random amplitude-modulated complex exponentials. A low-pass filter is used to suppress noise and approximate the data as a sum of weighted complex exponentials. The annihilating filter used in FRI approaches is applied to estimate the sampling interval in the closed form. We perform experiments on simulated and real data, and assess both objective and subjective performances in comparison with the state-of-the-art ellipse fitting methods. The proposed method produces ellipses with lesser bias. Furthermore, the mean-squared error is lesser by about 2 to 10 dB. We show the applications of ellipse fitting in iris images starting from partial edge contours, and to free-hand ellipses drawn on a touch-screen tablet.

A hybrid finite element formulation for flexible multibody dynamics

This work deals with the transient analysis of flexible multibody systems within a hybrid finite element framework. Hybrid finite elements are based on a two-field variational formulation in which the displacements and stresses are interpolated separately yielding very good coarse mesh accuracy. Most of the literature on flexible multibody systems uses beam-theory-based formulations. In contrast, the use of hybrid finite elements uses continuum-based elements, thus avoiding the problems associated with rotational degrees of freedom. In particular, any given three-dimensional constitutive relations can be directly used within the framework of this formulation. Since the coarse mesh accuracy as compared to a conventional displacement-based formulation is very high, the scheme is cost effective as well. A general formulation is developed for the constrained motion of a given point on a line manifold, using a total Lagrangian method. The multipoint constraint equations are implemented using Lagrange multipliers. Various kinds of joints such as cylindrical, prismatic, and screw joints are implemented within this general framework. Hinge joints such as spherical, universal, and revolute joints are obtained simply by using shared nodes between the bodies. In addition to joints, the formulation and implementation details for a DC motor actuator and for prescribed relative rotation are also presented. Several example problems illustrate the efficacy of the developed formulation.