Automatic finite element formulation and assembly of hyperelastic higher order structural models

The formulation of higher order structural models and their discretization using the finite element method is difficult owing to their complexity, especially in the presence of non-linearities. In this work a new algorithm for automating the formulation and assembly of hyperelastic higher-order structural finite elements is developed. A hierarchic series of kinematic models is proposed for modeling structures with special geometries and the algorithm is formulated to automate the study of this class of higher order structural models. The algorithm developed in this work sidesteps the need for an explicit derivation of the governing equations for the individual kinematic modes. Using a novel procedure involving a nodal degree-of-freedom based automatic assembly algorithm, automatic differentiation and higher dimensional quadrature, the relevant finite element matrices are directly computed from the variational statement of elasticity and the higher order kinematic model. Another significant feature of the proposed algorithm is that natural boundary conditions are implicitly handled for arbitrary higher order kinematic models. The validity algorithm is illustrated with examples involving linear elasticity and hyperelasticity. (C) 2013 Elsevier Inc. All rights reserved.

Analysis of cracked and un-cracked semicircular rings under symmetric loading

Edge cracked specimens have been widely utilized for fracture testing. Edge cracked semicircular disk (ECSD) specimen has now been well characterized with regard to its form factor and weight function. This paper presents a modified semicircular ring version of this specimen to enhance the form factor in general while retaining other desirable features. The efficacy of the modified design is proved by combining theory of elasticity solutions with finite element results to arrive at the optimum design geometry. New insights emerging from this work are used to theoretically re-examine the arch-tension and the four-point bend specimens. (C) 2014 Elsevier Ltd. All rights reserved.

An unusual demixing behavior in PS-PVME blends in the presence of nanoparticles

The effect of silver nanoparticles (sNP) on the demixing and the evolution of morphology in off-critical blends of 90/10 and 10/90 (wt/wt) PS/PVME polystyrene/poly(vinyl methyl ether)] was probed here using shear rheology and optical microscopy. The faster component (PVME) has a higher molecular weight (80 kDa) than the slower component (PS, 35 kDa), which makes this system quite interesting to study with respect to the evolving morphology, as the blends transit through the binodal and the spinodal envelopes. An unusual demixing behavior was observed in both PVME rich and PS rich blends. Temperature modulated differential scanning calorimetry measurements showed that the T-g value for the blends with sNP was slightly lower than that of the neat blends. A decreased volume of cooperativity at T-g suggests confined segmental dynamics in the presence of sNP. Although, the addition of sNP had no influence on the thermodynamic demixing temperature, it significantly altered the elasticity of the minor component during the transition of the blend from the homogeneous to the heterogeneous state. This is manifested from energetically driven localization of the sNP in the PVME phase during demixing. As a direct consequence of this, the formation of the microstructures upon demixing was observed to be delayed in the presence of sNP. Interestingly, in the intermediate quench depth, the higher viscoelastic phase evolved as an interconnected network, which subsequently coarsened into discrete droplets in the late stages for the 90/10 PS/PVME blends. Similar observations were made for 10/90 PS/PVME blends where threads of PVME appeared at deeper quench depths in the presence of sNP. The interconnected network formation of the minor phase (here PVME), which is also the faster component in the blend, was different from the usual demixing behavior.

Designing Elastic Organic Crystals: Highly Flexible Polyhalogenated N-Benzylideneanilines

The intermolecular interactions and structural features in crystals of seven halogenated N-benzylideneanilines (Schiff bases), all of which exhibit remarkable flexibility, were examined to identify the common packing features that are the raison d'etre for the observed elasticity. The following two features, in part related, were identified as essential to obtain elastic organic crystals: 1)A multitude of weak and dispersive interactions, including halogen bonds, which may act as structural buffers for deformation through easy rupture and reformation during bending; and 2)corrugated packing patterns that would get interlocked and, in the process, prevent long-range sliding of molecular planes.

A Smooth Discretization Bridging Finite Element and Mesh-free Methods Using Polynomial Reproducing Simplex Splines

This work sets forth a `hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (Cp-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of mesh-free (MF) schemes. Numerical experiments show the proposed method requires lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form. Numerical demonstrations of optimal convergence rates for a few test cases are given and the method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model.

Hybrid Elements for Modelling Squeeze Film Effects Coupled with Structural Interactions in Vibratory MEMS Devices

We present a hybrid finite element based methodology to solve the coupled fluid structure problem of squeeze film effects in vibratory MEMS devices, such as gyroscopes, RF switches, and 2D resonators. The aforementioned devices often have a thin plate like structure vibrating normally to a fixed substrate, and are generally not perfectly vacuum packed. This results in a thin air film being trapped between the vibrating plate and the fixed substrate which behaves like a squeeze film offering both stiffness and damping. For accurate modelling of such devices the squeeze film effects must be incorporated. Extensive literature is available on squeeze film modelling, however only a few studies address the coupled fluid elasticity problem. The majority of the studies that account for the plate elasticity coupled with the fluid equation, either use approximate mode shapes for the plate or use iterative solution strategies. In an earlier work we presented a single step coupled methodology using only one type of displacement based element to solve the coupled problem. The displacement based finite element models suffer from locking issues when it comes to modelling very thin structures with the lateral dimensions much larger than the plate thickness as is typical in MEMS devices with squeeze film effects. In this work we present another coupled formulation where we have used hybrid elements to model the structural domain. The numerical results show a huge improvement in convergence and accuracy with coarse hybrid mesh as compared to displacement based formulations. We further compare our numerical results with experimental data from literature and find them to be in good accordance.

CROSS-SECTIONAL ANALYSIS OF PRE-TWISTED THICK BEAMS USING VARIATIONAL ASYMPTOTIC METHOD

An asymptotically-exact methodology is presented for obtaining the cross-sectional stiffness matrix of a pre-twisted moderately-thick beam having rectangular cross sections and made of transversely isotropic materials. The anisotropic beam is modeled from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy of the beam is computed making use of the constitutive law and the kinematical relations derived with the inclusion of geometrical nonlinearities and initial twist. Large displacements and rotations are allowed, but small strain is assumed. The Variational Asymptotic Method is used to minimize the energy functional, thereby reducing the cross section to a point on the reference line with appropriate properties, yielding a 1-D constitutive law. In this method as applied herein, the 2-D cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged as orders of the small parameters. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that renders the 1-D strain measures well-defined. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

COMPUTATIONAL MODELING OF A MULTI-LAYERED PIEZO-COMPOSITE BEAM MADE UP OF MFC

The Variational Asymptotic Method (VAM) is used for modeling a coupled non-linear electromechanical problem finding applications in aircrafts and Micro Aerial Vehicle (MAV) development. VAM coupled with geometrically exact kinematics forms a powerful tool for analyzing a complex nonlinear phenomena as shown previously by many in the literature 3 - 7] for various challenging problems like modeling of an initially twisted helicopter rotor blades, matrix crack propagation in a composite, modeling of hyper elastic plates and various multi-physics problems. The problem consists of design and analysis of a piezocomposite laminate applied with electrical voltage(s) which can induce direct and planar distributed shear stresses and strains in the structure. The deformations are large and conventional beam theories are inappropriate for the analysis. The behavior of an elastic body is completely understood by its energy. This energy must be integrated over the cross-sectional area to obtain the 1-D behavior as is typical in a beam analysis. VAM can be used efficiently to approximate 3-D strain energy as closely as possible. To perform this simplification, VAM makes use of thickness to width, width to length, width multiplied by initial twist and strain as small parameters embedded in the problem definition and provides a way to approach the exact solution asymptotically. In this work, above mentioned electromechanical problem is modeled using VAM which breaks down the 3-D elasticity problem into two parts, namely a 2-D non-linear cross-sectional analysis and a 1-D non-linear analysis, along the reference curve. The recovery relations obtained as a by-product in the cross-sectional analysis earlier are used to obtain 3-D stresses, displacements and velocity contours. The piezo-composite laminate which is chosen for an initial phase of computational modeling is made up of commercially available Macro Fiber Composites (MFCs) stacked together in an arbitrary lay-up and applied with electrical voltages for actuation. The expressions of sectional forces and moments as obtained from cross-sectional analysis in closed-form show the electro-mechanical coupling and relative contribution of electric field in individual layers of the piezo-composite laminate. The spatial and temporal constitutive law as obtained from the cross-sectional analysis are substituted into 1-D fully intrinsic, geometrically exact equilibrium equations of motion and 1-D intrinsic kinematical equations to solve for all 1-D generalized variables as function of time and an along the reference curve co-ordinate, x(1).

Analysis of Transverse Shear Strains in Pre-Twisted Thick Beams using Variational Asymptotic Method

The cross-sectional stiffness matrix is derived for a pre-twisted, moderately thick beam made of transversely isotropic materials and having rectangular cross sections. An asymptotically-exact methodology is used to model the anisotropic beam from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy is computed making use of the beam constitutive law and kinematical relations derived with the inclusion of geometrical nonlinearities and an initial twist. The energy functional is minimized making use of the Variational Asymptotic Method (VAM), thereby reducing the cross section to a point on the beam reference line with appropriate properties, forming a 1-D constitutive law. VAM is a mathematical technique employed in the current problem to rigorously split the 3-D analysis of beams into two: a 2-D analysis over the beam cross-sectional domain, which provides a compact semi-analytical form of the properties of the cross sections, and a nonlinear 1-D analysis of the beam reference curve. In this method, as applied herein, the cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged in orders of the small parameters. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that render the 1-D strain measures well-defined. The zeroth-order 3-D warping field thus yielded is then used to integrate the 3-D strain energy density over the cross section, resulting in the 1-D strain energy density, which in turn helps identify the corresponding cross-sectional stiffness matrix. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

VAM Applied to Dimensional Reduction of Nonlinear Multifunctional Film-fabric Laminates

This work aims at asymptotically accurate dimensional reduction of non-linear multi-functional film-fabric laminates having specific application in design of envelopes for High Altitude Airships (HAA). The film-fabric laminate for airship envelope consists of a woven fabric core coated with thin films on each face. These films provide UV protection and Helium leakage prevention, while the core provides required structural strength. This problem is both geometrically and materially non-linear. To incorporate the geometric non-linearity, generalized warping functions are used and finite deformations are allowed. The material non-linearity is handled by using hyper-elastic material models for each layer. The development begins with three-dimensional (3-D) nonlinear elasticity and mathematically splits the analysis into a one-dimensional through-the-thickness analysis and a two-dimensional (2-D) plate analysis. The through-the-thickness analysis provides the 2-D constitutive law which is then given as an input to the 2-D reference surface analysis. The dimensional reduction is carried out using Variational Asymptotic Method (VAM) for moderate strains and very small thickness-to-wavelength ratio. It features the identification and utilization of additional small parameters such as ratio of thicknesses and stiffness coefficients of core and films. Closed form analytical expressions for warping functions and 2-D constitutive law of the film-fabric laminate are obtained.

Unified Continuum Damage Model for Matrix Cracking in Composite Rotor Blades

This paper deals with modeling of the first damage mode, matrix micro-cracking, in helicopter rotor/wind turbine blades and how this effects the overall cross-sectional stiffness. The helicopter/wind turbine rotor system operates in a highly dynamic and unsteady environment leading to severe vibratory loads present in the system. Repeated exposure to this loading condition can induce damage in the composite rotor blades. These rotor/turbine blades are generally made of fiber-reinforced laminated composites and exhibit various competing modes of damage such as matrix micro-cracking, delamination, and fiber breakage. There is a need to study the behavior of the composite rotor system under various key damage modes in composite materials for developing Structural Health Monitoring (SHM) system. Each blade is modeled as a beam based on geometrically non-linear 3-D elasticity theory. Each blade thus splits into 2-D analyzes of cross-sections and non-linear 1-D analyzes along the beam reference curves. Two different tools are used here for complete 3-D analysis: VABS for 2-D cross-sectional analysis and GEBT for 1-D beam analysis. The physically-based failure models for matrix in compression and tension loading are used in the present work. Matrix cracking is detected using two failure criterion: Matrix Failure in Compression and Matrix Failure in Tension which are based on the recovered field. A strain variable is set which drives the damage variable for matrix cracking and this damage variable is used to estimate the reduced cross-sectional stiffness. The matrix micro-cracking is performed in two different approaches: (i) Element-wise, and (ii) Node-wise. The procedure presented in this paper is implemented in VABS as matrix micro-cracking modeling module. Three examples are presented to investigate the matrix failure model which illustrate the effect of matrix cracking on cross-sectional stiffness by varying the applied cyclic

FRP-Confined Clay Brick Masonry Assemblages under Axial Compression: Experimental and Analytical Investigations

Unreinforced masonry (URM) structures that are in need of repair and rehabilitation constitute a significant portion of building stock worldwide. The successful application of fiber-reinforced polymers (FRP) for repair and retrofitting of reinforced-concrete (RC) structures has opened new avenues for strengthening URM structures with FRP materials. The present study analyzes the behavior of FRP-confined masonry prisms under monotonic axial compression. Masonry comprising of burnt clay bricks and cement-sand mortar (generally adopted in the Indian subcontinent) having E-b/E-m ratio less than one is employed in the study. The parameters considered in the study are, (1) masonry bonding pattern, (2) inclination of loading axis to the bed joint, (3) type of FRP (carbon FRP or glass FRP), and (4) grade of FRP fabric. The performance of FRP-confined masonry prisms is compared with unconfined masonry prisms in terms of compressive strength, modulus of elasticity and stress-strain response. The results showed an enhancement in compressive strength, modulus of elasticity, strain at peak stress, and ultimate strain for FRP-confined masonry prisms. The FRP confinement of masonry resulted in reducing the influence of the inclination of the loading axis to the bed joint on the compressive strength and failure pattern. Various analytical models available in the literature for the prediction of compressive strength of FRP-confined masonry are assessed. New coefficients are generated for the analytical model by appending experimental results of the current study with data available in the literature. (C) 2014 American Society of Civil Engineers.

Dual Stress and Thermally Driven Mechanical Properties of the Same Organic Crystal: 2,6-Dichlorobenzylidene-4-fluoro-3-nitroaniline

An elastic organic crystal, 2,6-dichlorobenzylidine-4-fluoro-3-nitroaniline (DFNA), which also shows thermosalient behavior, is studied. The presence of these two distinct properties in the same crystal is unusual and unprecedented because they follow respectively from isotropy and anisotropy in the crystal packing. Therefore, while both properties lead from the crystal structure, the mechanisms for bending and thermosalience are quite independent of one another. Crystals of the low-temperature (a) form of the title compound are bent easily without any signs of fracture with the application of deforming stress, and this bending is within the elastic limit. The crystal structure of the a-form was determined (P2(1)/c, Z = 4, a = 3.927(7) angstrom, b = 21.98(4) angstrom, c = 15.32(3) angstrom). There is an irreversible phase transition at 138 degrees C of this form to the high-temperature beta-form followed by melting at 140 degrees C. Variable-temperature X-ray powder diffraction was used to investigate the structural changes across the phase transition and, along with an FTIR study, establishes the structure of the beta-form. A possible rationale for strain build-up is given. Thermosalient behavior arises from anisotropic changes in the three unit cell parameters across the phase transition, notably an increase in the b axis parameter from 21.98 to 22.30 angstrom. A rationale is provided for the existence of both elasticity and thermosalience in the same crystal. FTIR studies across the phase transition reveal important mechanistic insights: (i) increased pi...pi repulsions along 100] lead to expansion along the a axis; (ii) change in alignment of C-Cl and NO2 groups result from density changes; and (iii) competition between short-range repulsive (pi...pi) interactions and long-range attractive dipolar interactions (C-Cl and NO2) could lie at the origin of the existence of two distinctive properties.

Micromechanical analysis of self compacting concrete

The micro-level properties of different self compacting concrete (SCC) mixes with and without mineral admixures are studied. The study considers SCC as a two phase material consisting of matrix and aggregate. Micro indentation technique is employed to obtain the hardness of individual phases and to compute the micro-property (modulus of elasticity). Using a self consistent homogenization procedure, the micro-property is scaled-up to obtain the macro-property which is shown to agree with the experimentally obtained macro values. It is seen that there exists a smaller interfacial transition zone at different ages of curing across all the mixes due to the presence of more fines in SCC. Also, there is no significant change in the property of the SCC having no fly ash or silica fume beyond 28 days whereas a substantial change in the micro and macro properties are seen in the SCC having fly ash and silica fume.

Interplay of Substrate Conductivity, Cellular Microenvironment, and Pulsatile Electrical Stimulation toward Osteogenesis of Human Mesenchymal Stem Cells in Vitro

The influences of physical stimuli such as surface elasticity, topography, and chemistry over mesenchymal stem cell proliferation and differentiation are well investigated. In this context, a fundamentally different approach was adopted, and we have demonstrated the interplay of inherent substrate conductivity, defined chemical composition of cellular microenvironment, and intermittent delivery of electric pulses to drive mesenchymal stem cell differentiation toward osteogenesis. For this, conducting polyaniline (PANI) substrates were coated with collagen type 1 (Coll) alone or in association with sulfated hyaluronan (sHya) to form artificial extracellular matrix (aECM), which mimics the native microenvironment of bone tissue. Further, bone marrow derived human mesenchymal stem cells (hMSCs) were cultured on these moderately conductive (10(-4)10(-3) S/cm) aECM coated PANI substrates and exposed intermittently to pulsed electric field (PEF) generated through transformer-like coupling (TLC) approach over 28 days. On the basis of critical analysis over an array of end points, it was inferred that Coll/sHya coated PANI (PANI/Coll/sHya) substrates had enhanced proliferative capacity of hMSCs up to 28 days in culture, even in the absence of PEF stimulation. On the contrary, the adopted PEF stimulation protocol (7 ms rectangular pulses, 3.6 mV/cm, 10 Hz) is shown to enhance osteogenic differentiation potential of hMSCs. Additionally, PEF stimulated hMSCs had also displayed different morphological characteristics as their nonstimulated counterparts. Concomitantly, earlier onset of ALP activity was also observed on PANI/Coll/sHya substrates and resulted in more calcium deposition. Moreover, real-time polymerase chain reaction results indicated higher mRNA levels of alkaline phosphatase and osteocalcin, whereas the expression of other osteogenic markers such as Runt-related transcription factor 2, Col1A, and osteopontin exhibited a dynamic pattern similar to control cells that are cultured in osteogenic medium. Taken together, our experimental results illustrate the interplay of multiple parameters such as substrate conductivity, electric field stimulation, and aECM coating on the modulation of hMSC proliferation and differentiation in vitro.