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.

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.

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.

A rectangular laminated anisotropic shallow thin shell finite element

This report contains the details of the development of the stiffness matrix for a rectangular laminated anisotropic shallow thin shell finite element. The derivation is done under linear thin shell assumptions. Expressing the assumed displacement state over the middle surface of the shell as products of one-dimensional first-order Hermite interpolation polynomials, it is possible to insure that the displacement state for the assembled set of such elements, to be geometrically admissible. Monotonic convergence of the total potential energy is therefore possible as the modelling is successively refined. The element is systematically evaluated for its performance considering various examples for which analytical or other solutions are available

A novel approach to design of cis-acting DNA structural element for regulation of gene expression in vivo

Taking advantage of the degeneracy of the genetic code we have developed a novel approach to introduce, within a gene, DNA sequences capable of adopting unusual structures and to investigate the role of such sequences in regulation of gene expression in vivo. We used a computer program that generates alternative codon sequences for the same amino-acid sequence to convert a stretch of nucleotides into an inverted-repeat sequence with the potential to adopt cruciform structure. This approach was used to replace a 51-base-pair EcoRI-HindIII segment in the N-terminal region of the beta-galactosidase gene in plasmid pUC19 with a 51-bp synthetic oligonucleotide sequence with the potential to adopt a cruciform structure with 18 bp in the stem region. In selecting the 51-bp sequence, care was taken to include those codons that are preferred in E. coli. E. coli DH5-alpha cells harbouring the plasmid containing the redesigned sequence showed drastic reduction in expression of the beta-galactosidase gene compared to cells harbouring the plasmid with the native sequence. This approach demonstrates the possibility of introducing DNA secondary-structure elements to alter regulation of gene expression in vivo.

New twelve node serendipity quadrilateral plate bending element based on Mindlin-Reissner theory using Integrated Force Method

The Integrated Force Method (IFM) is a novel matrix formulation developed for analyzing the civil, mechanical and aerospace engineering structures. In this method all independent/internal forces are treated as unknown variables which are calculated by simultaneously imposing equations of equilibrium and compatibility conditions. This paper presents a new 12-node serendipity quadrilateral plate bending element MQP12 for the analysis of thin and thick plate problems using IFM. The Mindlin-Reissner plate theory has been employed in the formulation which accounts the effect of shear deformation. The performance of this new element with respect to accuracy and convergence is studied by analyzing many standard benchmark plate bending problems. The results of the new element MQP12 are compared with those of displacement-based 12-node plate bending elements available in the literature. The results are also compared with exact solutions. The new element MQP12 is free from shear locking and performs excellent for both thin and moderately thick plate bending situations.

Spectral element approach to wave propagation in uncertain beam structures

This paper presents a study on the uncertainty in material parameters of wave propagation responses in metallic beam structures. Special effort is made to quantify the effect of uncertainty in the wave propagation responses at high frequencies. Both the modulus of elasticity and the density are considered uncertain. The analysis is performed using a Monte Carlo simulation (MCS) under the spectral finite element method (SEM). The randomness in the material properties is characterized by three different distributions, the normal, Weibull and extreme value distributions. Their effect on wave propagation in beams is investigated. The numerical study shows that the CPU time taken for MCS under SEM is about 48 times less than for MCS under a conventional one-dimensional finite element environment for 50 kHz loading. The numerical results presented investigate effects of material uncertainties on high frequency modes. A study is performed on the usage of different beam theories and their uncertain responses due to dynamic impulse load. These studies show that even for a small coefficient of variation, significant changes in the above parameters are noticed. A number of interesting results are presented, showing the true effects of uncertainty response due to dynamic impulse load.

A Smooth Finite Element Method Based on Reproducing Kernel DMS-Splines

The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials There is thus a case for combining these advantages in a so-called hybrid scheme or a `smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C-p (p >= 1) continuity One such recent attempt, a NURBS based parametric bridging method (Shaw et al 2008b), uses polynomial reproducing, tensor-product non-uniform rational B-splines (NURBS) over a typical FE mesh and relies upon a (possibly piecewise) bijective geometric map between the physical domain and a rectangular (cuboidal) parametric domain The present work aims at a significant extension and improvement of this concept by replacing NURBS with DMS-splines (say, of degree n > 0) that are defined over triangles and provide Cn-1 continuity across the triangle edges This relieves the need for a geometric map that could precipitate ill-conditioning of the discretized equations Delaunay triangulation is used to discretize the physical domain and shape functions are constructed via the polynomial reproduction condition, which quite remarkably relieves the solution of its sensitive dependence on the selected knotsets Derivatives of shape functions are also constructed based on the principle of reproduction of derivatives of polynomials (Shaw and Roy 2008a) Within the present scheme, the triangles also serve as background integration cells in weak formulations thereby overcoming non-conformability issues Numerical examples involving the evaluation of derivatives of targeted functions up to the fourth order and applications of the method to a few boundary value problems of general interest in solid mechanics over (non-simply connected) bounded domains in 2D are presented towards the end of the paper

Development of a New Finite Element for the Analysis of Sandwich Beams with Soft Core

A new super convergent sandwich beam finite element formulation is presented in this article. This element is a two-nodded, six degrees of freedom (dof) per node (3 dof u(0), w, phi for top and bottom face sheets each), which assumes that all the axial and flexural loads are taken by face sheets, while the core takes only the shear loads. The beam element is formulated based on first-order shear deformation theory for the face sheets and the core displacements are assumed to vary linearly across the thickness. A number of numerical experiments involving static, free vibration, and wave propagation analysis examples are solved with an aim to show the super convergent property of the formulated element. The examples presented in this article consider both metallic and composite face sheets. The formulated element is verified in most cases with the results available in the published literature.

A field-consistent and variationally correct representation of transverse shear strains in the nine-noded plate element

We present a biquadratic Lagrangian plate bending element with consistent fields for the constrained transverse shear strain functions. A technique involving expansion of the strain interpolations in terms of Legendre polynomials is used to redistribute the kinematically derived shear strain fields so that the field-consistent forms (i.e. avoiding locking) are also variationally correct (i.e. do not violate the variational norms). Also, a rational method of isoparametric Jacobian transformation is incorporated so that the constrained covariant shear strain fields are always consistent in whatever general quadrilateral form the element may take. Finally the element is compared with another formulation which was recently published. The element is subjected to several robust bench mark tests and is found to pass all the tests efficiently.

A finite element analysis of dynamic fracture initiation by ductile failure mechanisms in a 4340 steel

In some recent dropweight impact experiments [5] with pre-notched bend specimens of 4340 steel, it was observed that considerable crack tunneling occurred in the interior of the specimen prior to gross fracture initiation on the free surfaces. The final failure of the side ligaments happened because of shear lip formation. The tunneled region is characterized by a flat, fibrous fracture surface. In this paper, the experiments of [5] (corresponding to 5 m/s impact speed) are analyzed using a plane strain, dynamic finite element procedure. The Gurson constitutive model that accounts for the ductile failure mechanisms of micro-void nucleation, growth and coalescence is employed. The time at which incipient failure was observed near the notch tip in this computation, and the value of the dynamic J-integral, J d, at this time, compare reasonably well with experiments. This investigation shows that J-controlled stress and deformation fields are established near the notch tip whenever J d , increases with time. Also, it is found that the evolution of micro-mechanical quantities near the notch root can be correlated with the time variation of J d .The strain rate and the adiabatic temperature rise experienced at the notch root are examined. Finally, spatial variations of stresses and deformations are analyzed in detail.