A C-0 Interior Penalty Method for a Fourth-order Variational Inequality of the Second Kind

In this article, we propose a C-0 interior penalty ((CIP)-I-0) method for the frictional plate contact problem and derive both a priori and a posteriori error estimates. We derive an abstract error estimate in the energy norm without additional regularity assumption on the exact solution. The a priori error estimate is of optimal order whenever the solution is regular. Further, we derive a reliable and efficient a posteriori error estimator. Numerical experiments are presented to illustrate the theoretical results. (c) 2015Wiley Periodicals, Inc.

Yielding of sensitive clays: micromechanical considerations

Test results reported on several natural sensitive soils show significant anisotropy of the yield curves, which are generally oriented along the coefficient of earth pressure at rest (K-0) axis. An attempt is made in this paper to explain the anisotropy in yielding from microstructural considerations. An elliptic pore, with particle domains aligned along the periphery of the pore, and with the major axis of the pore being oriented along the direction of the in situ major principal stress, is chosen as the unit of microstructure. An analysis of forces at the interdomain contacts around the ellipse is carried out with reference to experimentally determined yield stress conditions of one soil, and a yield criteria is defined. The analysis, with the proposed yield criteria, enables one to define the complete yield curve for any other soil from the results of only two tests (one constant eta compression test with eta close to eta(K?0), where eta is the stress ratio (= q/p) and eta(K?0) is the stress ratio corresponding to anisotropic K-0 compression, and another undrained shear test). Predicted yield curves are compared with experimental yield curves of several soils reported in the literature.

Approximate solution of u'=-A2u using a relation between exp(-tA2) and exp(itA)

Let A be a positive definite operator in a Hilbert space and consider the initial value problem for u(t) = -A(2)u. Using a representation of the semigroup exp(-A(2)t) in terms of the group exp(iAt) we express u in terms of the solution of the standard heat equation w(t) = W-yy, with initial values v solving the initial value problem for v(y) = iAv. This representation is used to construct a method for approximating u in terms of approximations of v. In the case that A is a 2(nd) order elliptic operator the method is combined with finite elements in the spatial variable and then reduces the solution of the 4(th) order equation for u to that of the 2(nd) order equation for v, followed by the solution of the heat equation in one space variable.

A comparative study of the energetics, structures, and mechanisms of the HCN H HNC and LiCN <-> LiNC isomerizations

The potential energy surfaces of the HCN<->HNC and LiCN<->LiNC isomerization processes were determined by ab initio theory using fully optimized triple-zeta double polarization types of basis sets. Both the MP2 corrections and the QCISD level of calculations were performed to correct for the electron correlation. Results show that electron correlation has a considerable influence on the energetics and structures. Analysis of the intramolecular bond rearrangement processes reveals that, in both cases, H (or Li+) migrates in an almost elliptic path in the plane of the molecule. In HCN<->HNC, the migrating hydrogen interacts with the in-plane pi,pi* orbitals of CN, leading to a decrease in the C-N bond order. In LiCN<->LiNC, Li+ does not interact with the corresponding pi,pi* orbitals of CN.

Energy-momentum conserving algorithm for nonlinear transient analysis within the framework of hybrid elements

This work deals with the formulation and implementation of an energy-momentum conserving algorithm for conducting the nonlinear transient analysis of structures, within the framework of stress-based hybrid elements. Hybrid elements, which are based on a two-field variational formulation, are much less susceptible to locking than conventional displacement-based elements within the static framework. We show that this advantage carries over to the transient case, so that not only are the solutions obtained more accurate, but they are obtained in fewer iterations. We demonstrate the efficacy of the algorithm on a wide range of problems such as ones involving dynamic buckling, complicated three-dimensional motions, et cetera.

Differential Games of Mixed Type with Control and Stopping Times

We study a zero sum differential game of mixed type where each player uses both control and stopping times. Under certain conditions we show that the value function for this problem exists and is the unique viscosity solution of the corresponding variational inequalities. We also show the existence of saddle point equilibrium for a special case of differential game.

On an analytical-numerical procedure for the analysis of cylindrical shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.

An absolute method for the determination of surface tension of liquids using pendent drop profiles

Drop formation at conical tips which is of relevance to metallurgists is investigated based on the principle of minimization of free energy using the variational approach. The dimensionless governing equations for drop profiles are computer solved using the fourth order Runge-Kutta method. For different cone angles, the theoretical plots of XT and ZT vs their ratio, are statistically analyzed, where XT and ZT are the dimensionless x and z coordinates of the drop profile at a plane at the conical tip, perpendicular to the axis of symmetry. Based on the mathematical description of these curves, an absolute method has been proposed for the determination of surface tension of liquids, which is shown to be preferable in comparison with the earlier pendent-drop profile methods.

The hole spin model for oxides

A model of mobile 0-holes hybrized with Cu-spins on a square lattice is examined. A variational groundstate wavefunction which interpolates smoothly between n.n. RVB and Néel limits gives a Néellike minimum. A hole in an AF lattice polarizes it locally and becomes quite mobile. Two n.n. holes attract. Finally we speculate how holes can stabilize a spin liquid state.

An investigation of bond formation in the weakly bound first excited 1Σ and lowest 3Σ states of HeH+

The role of the electronic kinetic energy and its Cartesian components is examined during the formation of the first excited 1�£ and the lowest 3�£ states of HeH+ employing wavefunctions of multi-configuration type with basis orbitals in elliptic coordinates. Results show that the bond formation in these states is preceded primarily by a charge transfer from H to He+ rather than by polarisation of the H-orbital by He+

On an analytical-numerical procedure for the analysis of cylindrical-shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.

Stress and strain-driven algorithmic formulations for finite strain viscoplasticity for hybrid and standard finite elements

This work deals with the formulation and implementation of finite deformation viscoplasticity within the framework of stress-based hybrid finite element methods. Hybrid elements, which are based on a two-field variational formulation, are much less susceptible to locking than conventional displacement-based elements. The conventional return-mapping scheme cannot be used in the context of hybrid stress methods since the stress is known, and the strain and the internal plastic variables have to be recovered using this known stress field.We discuss the formulation and implementation of the consistent tangent tensor, and the return-mapping algorithm within the context of the hybrid method. We demonstrate the efficacy of the algorithm on a wide range of problems.

A refined analysis of composite laminates

The purpose of this paper is to develop a sufficiently accurate analysis, which is much simpler than exact three-dimensional analysis, for statics and dynamics of composite laminates. The governing differential equations and boundary conditions are derived by following a variational approach. The displacements are assumed piecewise linear across the thickness and the effects of transverse shear deformations and rotary inertia are included. A procedure for obtaining the general solution of the above governing differential equations in the form of hyperbolic-trigonometric series is given. The accuracy of the present theory is assessed by obtaining results for free vibrations and flexure of simply supported rectangular laminates and comparing them with results from exact three-dimensional analysis.

Vibration of skew plates

The vibration problems of skew plates with different edge conditions involving simple support and clamping have been considered by using the variational method of Ritz, a double series of beam characteristic functions being employed appropriate to the combination of the edge conditions. Natural frequencies and modes of vibration have been obtained for different combinations of side ratio and skew angle. These detailed studies reveal several interesting features concerning the frequency curves and nodal patterns. The results presented should, in addition, be of considerable value and practical significance in design applications.

On the definition of a small orifice

An attempt is made to draw a line of demarcation between small orifices and large orifices. It is proposed that an orifice can be considered 'small' if the discharge through it calculated on the small-orifice assumption differs from the exact discharge by less than half of one per cent. Using this criterion, it is shown that a circular or elliptic orifice can be deemed 'small' as long as the ratio of the depth of the orifice to the head causing the flow (measured from the center of the orifice to the liquid surface) is less than 0.8; a rectangular orifice can be deemed 'small' if the ratio is less than 0.7. A correction factor is suggested for the coefficient of discharge to account for the deviation from the exact discharge.