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 Legendre spectral element model for sloshing and acoustic analysis in nearly incompressible fluids

A new spectral finite element formulation is presented for modeling the sloshing and the acoustic waves in nearly incompressible fluids. The formulation makes use of the Legendre polynomials in deriving the finite element interpolation shape functions in the Lagrangian frame of reference. The formulated element uses Gauss-Lobatto-Legendre quadrature scheme for integrating the volumetric stiffness and the mass matrices while the conventional Gauss-Legendre quadrature scheme is used on the rotational stiffness matrix to completely eliminate the zero energy modes, which are normally associated with the Lagrangian FE formulation. The numerical performance of the spectral element formulated here is examined by doing the inf-sup test oil a standard rectangular rigid tank partially filled with liquid The eigenvalues obtained from the formulated spectral element are compared with the conventional equally spaced node locations of the h-type Lagrangian finite element and the predicted results show that these spectral elements are more accurate and give superior convergence The efficiency and robustness of the formulated elements are demonstrated by solving few standard problems involving free vibration and dynamic response analysis with undistorted and distorted spectral elements. and the obtained results are compared with available results in the published literature (C) 2009 Elsevier Inc All rights reserved

Fracture analysis of stiffened panels under combined tensile, bending, and shear loads

The objective is to present the formulation of numerically integrated modified virtual crack closure integral technique for concentrically and eccentrically stiffened panels for computation of strain-energy release rate and stress intensity factor based on linear elastic fracture mechanics principles. Fracture analysis of cracked stiffened panels under combined tensile, bending, and shear loads has been conducted by employing the stiffened plate/shell finite element model, MQL9S2. This model can be used to analyze plates with arbitrarily located concentric/eccentric stiffeners, without increasing the total number of degrees of freedom, of the plate element. Parametric studies on fracture analysis of stiffened plates under combined tensile and moment loads have been conducted. Based on the results of parametric,studies, polynomial curve fitting has been carried out to get best-fit equations corresponding to each of the stiffener positions. These equations can be used for computation of stress intensity factor for cracked stiffened plates subjected to tensile and moment loads for a given plate size, stiffener configuration, and stiffener position without conducting finite element analysis.

Structural stability of slender aerospace vehicles: Part I Mathematical modeling

A detailed mechanics based model is developed to analyze the problem of structural instability in slender aerospace vehicles. Coupling among the rigid-body modes, the longitudinal vibrational modes and the transverse vibrational modes due to asymmetric lifting-body cross-section are considered. The model also incorporates the effects of aerodynamic pressure and the propulsive thrust of the vehicle. The model is one-dimensional, and it can be employed to idealized slender vehicles with complex shapes. Condition under which a flexible body with internal stress waves behaves like a perfect rigid body is derived. Two methods are developed for finite element discretization of the system: (1) A time-frequency Fourier spectral finite element method and (2) h-p finite element method. Numerical results using the above methods are presented in Part II of this paper. (C) 2010 Elsevier Ltd. All rights reserved.

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.

Wave propagation analysis in carbon nanotube embedded composite using wavelet based spectral finite elements

In this paper, elastic wave propagation is studied in a nanocomposite reinforced with multiwall carbon nanotubes (CNTs). Analysis is performed on a representative volume element of square cross section. The frequency content of the exciting signal is at the terahertz level. Here, the composite is modeled as a higher order shear deformable beam using layerwise theory, to account for partial shear stress transfer between the CNTs and the matrix. The walls of the multiwall CNTs are considered to be connected throughout their length by distributed springs, whose stiffness is governed by the van der Waals force acting between the walls of nanotubes. The analyses in both the frequency and time domains are done using the wavelet-based spectral finite element method (WSFEM). The method uses the Daubechies wavelet basis approximation in time to reduce the governing PDE to a set of ODEs. These transformed ODEs are solved using a finite element (FE) technique by deriving an exact interpolating function in the transformed domain to obtain the exact dynamic stiffness matrix. Numerical analyses are performed to study the spectrum and dispersion relations for different matrix materials and also for different beam models. The effects of partial shear stress transfer between CNTs and matrix on the frequency response function (FRF) and the time response due to broadband impulse loading are investigated for different matrix materials. The simultaneous existence of four coupled propagating modes in a double-walled CNT-composite is also captured using modulated sinusoidal excitation.

Elasto-plastic Cracking Analysis of Reinforced Concrete

A new elasto-plastic cracking constitutive model for reinforced concrete is presented. The nonlinear effects considered cover almost all the nonlinearities exhibited by reinforced concrete under short term monotonic loading. They include concrete cracking in tension, plasticity in compression, aggregate interlock, tension softening, elasto-plastic behavior of steel, bond-slip between concrete, and steel reinforcement and tension stiffening. A new procedure for incorporating bondslip in smeared steel elements is described. A modified Huber-Hencky-Mises failure criterion for plastic deformation of concrete, which fits the experimental results under biaxial stresses better, is proposed. Multiple cracking at Gauss points and their opening and closing are considered. Matrix expressions are developed and are incorporated in a nonlinear finite element program. After the objectivity of the model is demonstrated, the model is used to analyze two different types of problems: one, a set of four shear panels, and the other, a reinforced concrete beam without shear reinforcement. The results of the analysis agree favorably with the experimental results.

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.

Thermomechanical generalization of pin joints in finite plates

The understanding of thermoelastic behaviour of joints is significant in order to ensure the integrity of large and complex structures exposed to a thermal environment, particularly in fields such as aerospace and nuclear engineering. Thermomechanical generalization of partial contact behaviour of a pin joint under combined in-plane mechanical loading and on-axis unidirectional heat flow has already been established by the authors for the analytically simpler domains of large plates. This paper successfully extends the on-going investigation to a single pin in a finite rectangular isotropic plate as a two-dimensional abstraction from a practical situation of a multipin fastener joint. The finite element method is used to analyse the joint problem under on-axis thermomechanical loading and unified load-contact relationships are established for a class of loading conditions.

Modeling of concrete cracking-A hybrid technique of using displacement discontinuity element method and direct boundary element method

This paper presents a new approach by making use of a hybrid method of using the displacement discontinuity element method and direct boundary element method to model concrete cracking by incorporating fictitious crack model. Fracture mechanics approach is followed using the Hillerborg's fictitious crack model. A boundary element based substructure method and a hybrid technique of using displacement discontinuity element method and direct boundary element method are compared in this paper. In order to represent the process zone ahead of the crack, closing forces are assumed to act in such a way that they obey a linear normal stress-crack opening displacement law. Plain concrete beams with and without initial crack under three-point loading were analyzed by both the methods. The numerical results obtained were shown to agree well with the results from existing finite element method. The model is capable of reproducing the whole range of load-deflection response including strain-softening and snap-back behavior as illustrated in the numerical examples. (C) 2011 Elsevier Ltd. All rights reserved.

Effect of Dilatation on the Elasto-Plastic Response of Bulk Metallic Glasses Under Indentation

Elasto-plastic response of bulk metallic glasses (BMGs) follows closely the response of granular materials through pressure dependent (or normal stress) yield locus and shear stress induced material dilatation. On a micro-structural level, material dilatation is responsible for stress softening and formation of localized shear band, however its influence on the macro-scale flow and deformation is largely unknown. In this work, we systematically analyze the effect of material dilatation on the gross indentation response of Zr-based BMG via finite element simulation. The strengthening/softening effect on the load-depth response and corresponding stress-strain profiles are presented in light of differences in elastic-plastic regimes under common indenters. Through comparison with existing experimental results, we draw conclusions regarding selection of suitable dilatation parameters for accurately predicting the gross response of BMGs

Design enhancement of the underwater noise control linings through a parsimonious parametric study

Control of sound transmission through the structure and reflection from the structure immersed in fluid media impose highly conflicting requirements on the design of the carpeted noise control linings. These requirements become even more stringent if the structure is expected to be moving with considerable speed particularly under intense hydrostatic pressure. Numerous configurations are possible for designing these linings. Therefore, in this paper, a few lining configurations are identified from the literature for parametric study so that the designer is provided with an environment to analyze and design the lining. A scheme of finite element analysis is used to analyze these linings for their acoustic performance. Commercial finite element software, NISA®, is used as a platform to develop a customized environment wherein design parameters of different configurations can be varied with consistency checks and generate the finite element meshes using the 8-noded hexahedral element. Four types of designs proposed and analysed here address the parameters of interest such as the echo reduction and the transmission loss. Study of the effect of different surface distributions of the cavities is carried out. Effect of static pressure on different designs is reported.

Influence of Spatial Scales on the Performance of Saturated Subsurface Flow and Transport Models

Numerical modeling of saturated subsurface flow and transport has been widely used in the past using different numerical schemes such as finite difference and finite element methods. Such modeling often involves discretization of the problem in spatial and temporal scales. The choice of the spatial and temporal scales for a modeling scenario is often not straightforward. For example, a basin-scale saturated flow and transport analysis demands larger spatial and temporal scales than a meso-scale study, which in turn has larger scales compared to a pore-scale study. The choice of spatial-scale is often dictated by the computational capabilities of the modeler as well as the availability of fine-scale data. In this study, we analyze the impact of different spatial scales and scaling procedures on saturated subsurface flow and transport simulations.

On an Integrated Transfer Matrix method for multiply connected mufflers

The commercial automotive mufflers are generally of a complicated shape with multiply connected parts and complex acoustic elements. The analysis of such complex mufflers has always been a great challenge. In this paper, an Integrated Transfer Matrix method has been developed to analyze complex mufflers. Integrated transfer matrix relates the state variables across the entire cross-section of the muffler shell, as one moves along the axis of the muffler, and can be partitioned appropriately in order to relate the state variables of different tubes constituting the cross-section. The paper presents a generalized one-dimensional (1-D) approach, using the transfer matrices of simple acoustic elements, which are available from the literature. The present approach is robust and flexible owing to its capability to construct an overall matrix of the muffler with the transfer matrices of individual acoustic elements and boundary conditions, which can then be used to evaluate the transmission loss, insertion loss, etc. Results from the present approach have been validated through comparisons with the available experimental and three-dimensional finite element method (FEM) based results. The results show good agreement with both measurements and FEM analysis up to the cut-off frequency. (C) 2011 Elsevier Ltd. All rights reserved.

Spectral Element Approach to Wave Propagation in Uncertain Composite Beam Structures

This paper presents a study of the wave propagation responses in composite structures in an uncertain environment. Here, the main aim of the work is to quantify the effect of uncertainty in the wave propagation responses at high frequencies. The material properties are considered uncertain and the analysis is performed using Neumann expansion blended with Monte Carlo simulation under the environment of spectral finite element method. The material randomness is included in the conventional wave propagation analysis by different distributions (namely, the normal and the Weibul distribution) and their effect on wave propagation in a composite beam is analyzed. The numerical results presented investigates the effect of material uncertainties on different parameters, namely, wavenumber and group speed, which are relevant in the wave propagation analysis. The effect of the parameters, such as fiber orientation, lay-up sequence, number of layers, and the layer thickness on the uncertain responses due to dynamic impulse load, is thoroughly analyzed. Significant changes are observed in the high frequency responses with the variation in the above parameters, even for a small coefficient of variation. High frequency impact loads are applied and a number of interesting results are presented, which brings out the true effects of uncertainty in the high frequency responses. [DOI: 10.1115/1.4003945]