Extending lagrange interpolation to develop a 4-node quadrilateral element

Today finite element method is a well established tool in engineering analysis and design. Though there axe many two and three dimensional finite elements available, it is rare that a single element performs satisfactorily in majority of practical problems. The present work deals with the development of 4-node quadrilateral element using extended Lagrange interpolation functions. The classical univariate Lagrange interpolation is well developed for 1-D and is used for obtaining shape functions. We propose a new approach to extend the Lagrange interpolation to several variables. When variables axe more than one the method also gives the set of feasible bubble functions. We use the two to generate shape function for the 4-node arbitrary quadrilateral. It will require the incorporation of the condition of rigid body motion, constant strain and Navier equation by imposing necessary constraints. The procedure obviates the need for isoparametric transformation since interpolation functions are generated for arbitrary quadrilateral shapes. While generating the element stiffness matrix, integration can be carried out to the accuracy desired by dividing the quadrilateral into triangles. To validate the performance of the element which we call EXLQUAD4, we conduct several pathological tests available in the literature. EXLQUAD4 predicts both stresses and displacements accurately at every point in the element in all the constant stress fields. In tests involving higher order stress fields the element is assured to converge in the limit of discretisation. A method thus becomes available to generate shape functions directly for arbitrary quadrilateral. The method is applicable also for hexahedra. The approach should find use for development of finite elements for use with other field equations also.

NMR relaxation study of disorder in condensed matter: solid solutions of betaine phosphate and phosphite

Proton NMR relaxation measurements have been carried out in anti-ferroelectric Betaine phosphate (BP), ferroelectric Betaine phosphite (BPI) and the mixed system BPI(1-x)BPx, at 11.4MHz and 23.3MHz from 300K to 80K for x=0.0, 0.25, 0.45, 0.85, and 1.0. The temperature dependence of spin lattice relaxation time T, exhibits two minima as expected from the BPP model in BP and BPI. The Larmor frequency dependence of T, in the mixed system is rather unusual and exhibits different slopes for the low temperature wings at the two frequencies, which is a clear experimental evidence of the presence of different methyl groups with different activation energies (E-a) indicating disorder.

Stability and transition in the flow of polymer solutions

The linear stability analysis of a plane Couette flow of viscoelastic fluid have been studied with the emphasis on two dimensional disturbances with wave number k similar to Re-1/2, where Re is Reynolds number based on maximum velocity and channel width. We employ three models to represent the dilute polymer solution: the classical Oldroyd-B model, the Oldroyd-B model with artificial diffusivity and the non-homogeneous polymer model. The result of the linear stability analysis is found to be sensitive to the polymer model used. While the plane Couette flow is found to be stable to infinitesimal disturbances for the first two models, the last one exhibits a linear instability.

A finite element analysis of quasistatic crack growth in a pressure sensitive constrained ductile layer

Polymeric adhesive layers are employed for bonding two components in a wide variety of technological applications, It has been observed that, unlike in metals, the yield behavior of polymers is affected by the state of hydrostatic stress. In this work, the effect of pressure sensitivity of yielding and layer thickness on quasistatic interfacial crack growth in a ductile adhesive layer is investigated. To this end, finite deformation, finite element analyses of a cracked sandwiched layer are carried out under plane strain, small-scale yielding conditions for a wide range of mode mixities. The Drucker-Prager constitutive equations are employed to represent the behavior of the layer. Crack propagation is simulated through a cohesive zone model, in which the interface is assumed to follow a prescribed traction-separation law. The results show that for a given mode mixity, the steady state Fracture toughness [K](ss) is enhanced as the degree of pressure sensitivity increases. Further, for a given level of pressure sensitivity, [K](ss) increases steeply as mode Il loading is approached. (C) 2000 Elsevier Science Ltd. All rights reserved.

Adhesion-Induced Instability in Asperities

Adhesive forces between two approaching asperities will deform the asperities, and under certain conditions this will result in a sudden runaway deformations leading to a jump-to-contact instability. We present finite element-based numerical studies on adhesion-induced deformation and instability in asperities. We consider the adhesive force acting on an asperity, when it is brought near a rigid half-space, due to van der Waals interaction between the asperity and the half-space. The adhesive force is considered to be distributed over the volume of the asperity (body force), thus resulting in more realistic simulations for the length scales considered. Iteration scheme based on a ``residual stress update'' algorithm is used to capture the effect of deformation on the adhesion force, and thereby the equilibrium configuration and the corresponding force. The numerical results are compared with the previous approximate analytical solutions for adhesion force, deformation of the asperity and adhesion-induced mechanical instability (jump-to-contact). It is observed that the instability can occur at separations much higher,and could possibly explain the higher value of instability separation observed in experiments. The stresses in asperities, particularly in case of small ones, are found to be high enough to cause yielding before jump -to-contact. The effect of roughness is considered by modeling a spherical protrusion on the hemispherical asperity.This small-scale roughness at the tip of the asperities is found to control the deformation behavior at small separations, and hence are important in determining the friction and wear due to the jump-to-contact instability.

A simple technique for approximate solutions of the Falkner-Skan equation

A simple, sufficiently accurate and efficient method for approximate solutions of the Falkner-Skan equation is proposed here for a wide range of the pressure gradient parameter. The proposed approximate solutions are obtained utilising a known solution of another differential equation.

Semi-similar solutions of the unsteady compressible second-order boundary layer flow at the stagnation point

Semi-similar solutions of the unsteady compressible laminar boundary layer flow over two-dimensional and axisymmetric bodies at the stagnation point with mass transfer are studied for all the second-order boundary layer effects when the free stream velocity varies arbitrarily with time. The set of partial differential equations governing the unsteady compressible second-order boundary layers representing all the effects are derived for the first time. These partial differential equations are solved numerically using an implicit finite-difference scheme. The results are obtained for two particular unsteady free stream velocity distributions: (a) an accelerating stream and (b) a fluctuating stream. It is observed that the total skin friction and heat transfer are strongly affected by the surface mass transfer and wall temperature. However, their variation with time is significant only for large times. The second-order boundary layer effects are found to be more pronounced in the case of no mass transfer or injection as compared to that for suction. Résumé Des solutions semi-similaires d'écoulement variable compressible de couche limite sur des corps bi-dimensionnels thermique, sont étudiées pour tous les effets de couche limite du second ordre, lorsque la vitesse de l'écoulement libre varie arbitrairement avec le temps. Le systéme d'équations aux dérivées partielles représentant tous les effets est écrit pour la premiére fois. On le résout numériquement á l'aide d'un schéma implicite aux différences finies. Les résultats sont obtenus pour deux cas de vitesse variable d'écoulement libre: (a) un écoulement accéléré et (b) un écoulement fluctuant. On observe que le frottement pariétal total et le transfert de chaleur sont fortement affectés par le transfert de masse et la température pariétaux. Néanmoins, leur variation avec le temps est sensible seulement pour des grandes durées. Les effets sont trouvés plus prononcés dans le cas de l'absence du transfert de masse ou de l'injection par rapport au cas de l'aspiration.

Extension of Darken's Quadratic Formalism to Dilute Multicomponent Solutions

Darken's quadratic formalism is extended to multicomponent solutions. Equations are developed for the representation of the integral and partial excess free energies, entropies and enthalpies in dilute multicomponent solutions. Quadratic formalism applied to multicomponent solutions is thermodynamically consistent. The formalism is compared with the conventional second order Maclaurin series or interaction parameter representation and the relations between them are derived. Advantages of the quadratic formalism are discussed.

Frequency-Dependent Viscosity of Membrane Solutions

It is shown that dilute suspensions of membranes have strongly frequency-dependent viscosities. This behaviour should be seen in a variety of measurements such as capillary flow, mechanical impedance and ultrasound damping.

Biosorption of iron(III) from aqueous solutions using the husk of Cicer arientinum

Iron is a major pollutant released as a by-product during several industrial operations especially during acid mining of metal ores. In this paper, the use of Bengal gram husk (husk of channa dal, Cicer arientinum) in the biosorption of Fe(III) from aqueous solutions is discussed. Parameters like agitation time, adsorbent dosage and pH were studied at different Fe(Ill) concentrations. The adsorption data fit well with Langmuir and Freundlich isotherm models. The adsorption capacity (q(max)) calculated from the Langmuir isotherm was 72.16 mg of Fe(III)/g of the biosorbent at an initial pH of 2.5. Desorption Studies were performed at different concentrations of hydrochloric acid showing that quantitative recovery of the metal ion is possible. The infrared spectra of the biomass before and after treatment with Fe(III), revealed that hydroxyl, carboxyl and amide bonds are involved in the uptake of Fe(III) ions.

A novel fiber-coated strong base-type anion exchanger with superfast kinetics. Removal and recovery of silver thiosulfate from aqueous solutions

A commercial acrylic fiber with 92% (w/w) acrylonitrile content was partially hydrolyzed converting a fraction of the nitrile (-CN) groups to carboxylic acid (-COOH) groups, to coat the fiber with polyethylenimine (PEI) resin, which was then crosslinked with glutaraldehyde and further quaternized with ethyl chloroacetate to produce a novel strong-base anionic exchanger in the form of fiber. Designated as PAN(QPEI.XG)(Cl-), the fibrous sorbent was compared with a commercial bead-form resin Amberlite IRA-458(Cl-) in respect of sorption capacity, selectivity, and kinetics for removal of silver thiosulfate complexes from aqueous solutions. Though the saturation level of [Ag(S2O3)(2)](3-) on PAN(QPEI.XG)(Cl-) is considerably less than that on IRA-458(Cl-), the gel-coated fibrous sorbent exhibits, as compared to the bead-form sorbent, a significantly higher sorption selectivity for the silver thiosulfate complex in the presence of excess of other anions Such as S2O32-, SO42-, and Cl-, and a remarkably faster rate of both sorption and stripping. The initial uptake of the sorbate by the fibrous sorbent is nearly instantaneous, reaching up to similar to 80% of the saturation capacity within 10 s, as compared to only similar to 12% on the bead-form sorbent. The high initial rate of uptake fits a shell-core kinetic model for sorption on fiber of cylindrical geometry. With 4M HCl, the stripping of the sorbed silver complex from the fibrous sorbent is clean and nearly instantaneous, while, in contrast, a much slower rate of stripping on the bead-form sorbent leads to its fouling due to a slow decomposition of the silver thiosulfate complex in the acidic medium.

Finite element analysis of shear critical prestressed SFRC beams

This study reports the details of the finite element analysis of eleven shear critical partially prestressed concrete T-beams having steel fibers over partial or full depth. Prestressed concrete T-beams having a shear span to depth ratio of 2.65 and 1.59 and failing in the shear have been analyzed Using 'ANSYS'. The 'ANSYS' model accounts for the nonlinear phenomenon, such as, bond-slip of longitudinal reinforcements, post-cracking tensile stiffness of the concrete, stress transfer across the cracked blocks of the concrete and load sustenance through the bridging of steel fibers at crack interlace. The concrete is modeled using 'SOLID65'-eight-node brick element, which is capable Of simulating the cracking and crushing behavior of brittle materials. The reinforcements such as deformed bars, prestressing wires and steel fibers have been modeled discretely Using 'LINK8' - 3D spar element. The slip between the reinforcement (rebar, fibers) and the concrete has been modeled using a 'COMBIN39'-non-linear spring element connecting the nodes of the 'LINK8' element representing the reinforcement and nodes of the 'SOLID65' elements representing the concrete. The 'ANSYS' model correctly predicted the diagonal tension failure and shear compression failure of prestressed concrete beams observed in the experiment. I-lie capability of the model to capture the critical crack regions, loads and deflections for various types Of shear failures ill prestressed concrete beam has been illustrated.

Using Water as a Design Element in Crystal Engineering. Host-Guest Compounds of Hydrated 3,5-Dihydroxybenzoic Acid

Thirteen host guest compounds of 3,5-dihydroxybenzoic acid (DHBA) have been structurally characterized. Water molecules occupy the peripheries of a hexagonal void, created with DHBA molecules, and act as ``hooks'' to connect the guest molecules with the host-framework via hydrogen bonding. The ``water hook'' is an OH group acting as a donor. Consequently, the guest molecules were chosen so that they contain good hydrogen bond acceptor functionalities. A number of multicomponent hydrates were isolated with stoichiometries (DHBA)(x)(H2O). (guest),. Of these, compounds with the following as guests were obtained as crystals that were good enough for single crystal work: ethyl acetate (EtOAc), diethyl oxalate, dimethyl oxalate, di(n-propyl) oxalate, diethyl malonate, diethyl succinate, chloroacetonitrile, N,N-dimethyl formamide (DMF), acetone, dimethyl sulfoxide (DMSO), 1-propanol, and 2-butanol. From 2-butanol, a hemihydrate, (DHBA)(2)(H2O), was also obtained concomitantly. Further to guest stabilization, water acts as a good mediator of effective crystal packing and also determines the topology of the host framework. En the present series of compounds, the role of water is wide ranging, and it is not easy to classify it specifically as a host or as a guest.