Nonlinear finite element analysis of reinforced concrete corbel

Reinforced concrete corbels have been analysed using the nonlinear finite element method. An elasto-plastic-cracking constitutive formulation using Huber-Hencky-Mises yield surface augmented with a tension cut-off is employed. Smeared-fixed cracking with mesh-dependent strain softening is employed to obtain objective results. Multiple non-orthogonal cracking and opening and closing of cracks are permitted. The model and the formulation are verified with respect to available numerical solution for an RC corbel. Results of analyses of nine reinforced concrete corbels are presented and compared with experimental results. Nonlinear finite element analysis of reinforced concrete structures is shown to be a complement and also a feasible alternative to laboratory testing.

Nonlinear hydrodynamics of a hard-sphere fluid near the glass transition

We conduct a numerical study of the dynamic behavior of a dense hard-sphere fluid by deriving and integrating a set of Langevin equations. The statics of the system is described by a free-energy functional of the Ramakrishnan-Yussouff form. We find that the system exhibits glassy behavior as evidenced through a stretched exponential decay and a two-stage relaxation of the density correlation function. The characteristic times grow with increasing density according to the Vogel-Fulcher law. The wave-number dependence of the kinetics is extensively explored. The connection of our results with experiment, mode-coupling theory, and molecular-dynamics results is discussed.

Path integral description of polymers using fractional Brownian walks

The statistical properties of fractional Brownian walks are used to construct a path integral representation of the conformations of polymers with different degrees of bond correlation. We specifically derive an expression for the distribution function of the chains’ end‐to‐end distance, and evaluate it by several independent methods, including direct evaluation of the discrete limit of the path integral, decomposition into normal modes, and solution of a partial differential equation. The distribution function is found to be Gaussian in the spatial coordinates of the monomer positions, as in the random walk description of the chain, but the contour variables, which specify the location of the monomer along the chain backbone, now depend on an index h, the degree of correlation of the fractional Brownian walk. The special case of h=1/2 corresponds to the random walk. In constructing the normal mode picture of the chain, we conjecture the existence of a theorem regarding the zeros of the Bessel function.

Spectral approach for the soliton and periodic solutions of the nonlinear wave equation

A spectral method that obtains the soliton and periodic solutions to the nonlinear wave equation is presented. The results show that the nonlinear group velocity is a function of the frequency shift as well as of the soliton power. When the frequency shift is a function of time, a solution in terms of the Jacobian elliptic function is obtained. This solution is periodic in nature, and, to generate such an optical pulse train, one must simultaneously amplitude- and frequency-modulate the optical carrier. Finally, we extend the method to include the effect of self-steepening.

A formula for the solution of general Abel integral equation

A new formula for the solution of the general Abel Integral equation is derived, and an important special case is checked with the known result.

A general procedure for modified crack closure integral in 3D problems with cracks

In linear elastic fracture mechanics (LEFM), Irwin's crack closure integral (CCI) is one of the signficant concepts for the estimation of strain energy release rates (SERR) G, in individual as well as mixed-mode configurations. For effective utilization of this concept in conjunction with the finite element method (FEM), Rybicki and Kanninen [Engng Fracture Mech. 9, 931 938 (1977)] have proposed simple and direct estimations of the CCI in terms of nodal forces and displacements in the elements forming the crack tip from a single finite element analysis instead of the conventional two configuration analyses. These modified CCI (MCCI) expressions are basically element dependent. A systematic derivation of these expressions using element stress and displacement distributions is required. In the present work, a general procedure is given for the derivation of MCCI expressions in 3D problems with cracks. Further, a concept of sub-area integration is proposed which facilitates evaluation of SERR at a large number of points along the crack front without refining the finite element mesh. Numerical data are presented for two standard problems, a thick centre-cracked tension specimen and a semi-elliptical surface crack in a thick slab. Estimates for the stress intensity factor based on MCCI expressions corresponding to eight-noded brick elements are obtained and compared with available results in the literature.

Growth and characterization of some novel crystals for nonlinear optical applications

The advent of high intensity lasers coupled with the recent advances in crystal technology has led to rapid progress in the field of nonlinear optics. This article traces the history of materials development that has taken place over the past forty odd years and dwells on the current status in this important area. The materials aspect is discussed under three classes viz. inorganic, organic and semiorganic crystals. In the end, some of the crystal growth work that has been carried out in author's laboratory is presented.

Modified Crack Closure Integral (MCCI) For 3-D Problems Using 20-Noded Brick Elements

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

Modified crack closure integral (MCCI) for 3-d problems using 20 noded brick elememnts

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

Nonlinear effects on rotating disk flow in a cylindrical casing

The flow due to a finite disk rotating in an incompressible viscous fluid has been studied. A modified Newton-gradient finite difference scheme is used to obtain the solution of full Navier-Stokes equations numerically for different disk and cylinder sizes for a wide range of Reynolds numbers. The introduction of the aspect ratio and the disk-shroud gap, significantly alters the flow characteristics in the region under consideration, The frictional torque calculated from the flow data reveals that the contribution due to nonlinear terms is not negligible even at a low Reynolds number. For large Reynolds numbers, the flow structure reveals a strong boundary layer character.

Performance of a pipelined ring algorithm for Fast Fourier Transform on transputer arrays

In this paper a pipelined ring algorithm is presented for efficient computation of one and two dimensional Fast Fourier Transform (FFT) on a message passing multiprocessor. The algorithm has been implemented on a transputer based system and experiments reveal that the algorithm is very efficient. A model for analysing the performance of the algorithm is developed from its computation-communication characteristics. Expressions for execution time, speedup and efficiency are obtained and these expressions are validated with experimental results obtained on a four transputer system. The analytical model is then used to estimate the performance of the algorithm for different number of processors, and for different sizes of the input data.

Physical characterization of the organic nonlinear optical crystal � 3-methoxy 4-hydroxy benzaldehyde

We highlight our recent experimental work on an efficient molecular nonlinear optical crystal, 3-methoxy 4-hydroxy benzaldehyde (MHBA). Optical quality single crystals of MHBA were grown from mixtures of solvents and from melt. The overall absorption and transparency window were improved by growing them in a mixture of chloroform and acetone. The grown crystals were characterized for their optical transmission, mechanical hardness and laser damage. We have observed a strong correlation between mechanical properties and laser induced damage.

A single network adaptive critic (SNAC) architecture for optimal control synthesis for a class of nonlinear systems

Even though dynamic programming offers an optimal control solution in a state feedback form, the method is overwhelmed by computational and storage requirements. Approximate dynamic programming implemented with an Adaptive Critic (AC) neural network structure has evolved as a powerful alternative technique that obviates the need for excessive computations and storage requirements in solving optimal control problems. In this paper, an improvement to the AC architecture, called the �Single Network Adaptive Critic (SNAC)� is presented. This approach is applicable to a wide class of nonlinear systems where the optimal control (stationary) equation can be explicitly expressed in terms of the state and costate variables. The selection of this terminology is guided by the fact that it eliminates the use of one neural network (namely the action network) that is part of a typical dual network AC setup. As a consequence, the SNAC architecture offers three potential advantages: a simpler architecture, lesser computational load and elimination of the approximation error associated with the eliminated network. In order to demonstrate these benefits and the control synthesis technique using SNAC, two problems have been solved with the AC and SNAC approaches and their computational performances are compared. One of these problems is a real-life Micro-Electro-Mechanical-system (MEMS) problem, which demonstrates that the SNAC technique is applicable to complex engineering systems.

A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.