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.

In-plane elastic responses of circular cell honeycombs and bundled circular tubes in a diamond array structure

Assemblages of circular tubes and circular honeycombs in close packed arrangement are presently both competing and complementing regular honeycomb structures (HCS). The intrinsic isotropy of bundled tubes/rings in hexagonal arrays restricts their use to applications with isotopic need. With the aim of extending the utility of tubes/rings assemblages to anisotropic needs, this paper explores the prospects of bundled tubes and circular honeycombs in a general diamond array structure (DAS) to cater these needs. To this end, effective transverse Young's moduli and Poisson's ratio for thick/thin DAS are obtained theoretically. Analysis frameworks including thin ring theory (TRT), curved beam theory (CBT) and elasticity formulations are tested and corroborated by FEA employing contact elements. Results indicate that TRT and CBT are reasonable for thin tubes and honeycombs. Nevertheless, TRT yields compact formulae to study the anisotropy ratio, moduli spectrum and sensitivity of the assemblage as a function of thicknesses and array structure. These formulae supplement designers as a guide to tailor the structures. On the other hand, elasticity formulation can estimate over a larger range including very thick tubes/rings. In addition, this formulation offers to estimate refined transverse strengths of assemblages. (C) 2015 Elsevier Ltd. All rights reserved.

In-plane effective shear modulus of generalized circular honeycomb structures and bundled tubes in a diamond array structure

Use of circular hexagonal honeycomb structures and tube assemblies in energy absorption systems has attracted a large number of literature on their characterization under crushing and impact loads. Notwithstanding these, effective shear moduli (G*) required for complete transverse elastic characterization and in analyses of hierarchical structures have received scant attention. In an attempt to fill this void, the present study undertakes to evaluate G* of a generalized circular honeycomb structures and tube assemblies in a diamond array structure (DAS) with no restriction on their thickness. These structures present a potential to realize a spectrum of moduli with minimal modifications, a point of relevance for manufactures and designers. To evaluate G* in this paper, models based on technical theories - thin ring theory and curved beam theory - and rigorous theory of elasticity are investigated and corroborated with FEA employing contact elements. Technical theories which give a good match for thin HCS offer compact expressions for moduli which can be harvested to study sensitivity of moduli on topology. On the other hand, elasticity model offers a very good match over a large range of thickness along with exact analysis of stresses by employing computationally efficient expressions. (C) 2015 Elsevier Ltd. All rights reserved.

Behavior of Brick-Mortar Interfaces in FRP-Strengthened Masonry Assemblages under Normal Loading and Shear Loading

Determination of shear strength of brick-mortar bed joint is critical to overcome the sliding-shear or joint-shear failure in masonry. In the recent past, researchers have attempted to enhance the shear strength and deformation capacity of brick-mortar bed joints by gluing fiber-reinforced polymer (FRP) composite across the bed joint. FRP composites offer several advantages like high strength-to-weight ratio, and ease of application in terms of labor, time, and reduced curing period. Furthermore, FRP composites are desirable for strengthening old masonry buildings having heritage value because of its minimal interference with the existing architecture. A majority of earlier studies on shear strengthening of masonry available in the literature adopted masonry having the ratio of modulus of elasticity of masonry unit (Emu) to modulus of elasticity of mortar (Em) greater than one. Information related to shear behavior of FRP glued masonry composed of masonry units having Young's modulus lower than mortar is limited. Hence the present study is focused on characterizing the interfacial behavior of brick-mortar bed joint of masonry assemblages composed of solid burnt clay bricks and cement-sand mortar (E-mu/E-m ratio less than one), strengthened with FRP composites. Masonry triplets and prisms with bed joint inclined to loading axis (0 degrees, 30 degrees, 45 degrees, 60 degrees and 90 degrees) are employed in this study. Glass and carbon FRP composites composed of bidirectional FRP fabric with equal density in both directions are used for strengthening masonry. Masonry triplets are glued with glass and carbon FRP composites in two configurations: (1) both faces of the triplet specimens are fully glued with GFRP composites; and (2) both faces of the triplet specimens are glued with GFRP and CFRP composites in strip form. The performance of masonry assemblages strengthened with FRP composites is assessed in terms of gain in shear strength, shear displacement, and postpeak behavior for various configurations and types of FRP composites considered. A semianalytical model is proposed for the prediction of shear strength of masonry bed joints glued with FRP composites. A composite failure envelope consisting of a Coulomb friction model and a compression cap is obtained for unreinforced masonry and GFRP-strengthened masonry based on the test results of masonry triplets and masonry prisms with bed joints having various inclinations to the loading (C) 2015 American Society of Civil Engineers.

Three-dimensional warping in strip-based composite four-bar mechanisms

This work intends to demonstrate the effect of geometrically non-linear cross-sectional analysis of certain composite beam-based four-bar mechanisms in predicting the three-dimensional warping of the cross-section. The only restriction in the present analysis is that the strains within each elastic body remain small (i.e., this work does not deal with materials exhibiting non-linear constitutive laws at the 3-D level). Here, all component bars of the mechanism are made of fiber-reinforced laminates. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. Each component of the mechanism is modeled as a beam based on geometrically non-linear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and non-linear 1-D analyses along the three beam reference curves. The splitting of the three-dimensional beam problem into two- and one-dimensional parts, called dimensional reduction, results in a tremendous savings of computational effort relative to the cost of three-dimensional finite element analysis, the only alternative for realistic beams. The analysis of beam-like structures made of laminated composite materials requires a much more complicated methodology. Hence, the analysis procedure based on Variational Asymptotic Method (VAM), a tool to carry out the dimensional reduction, is used here. The representative cross-sections of all component bars are analyzed using two different approaches: (1) Numerical Model and (2) Analytical Model. Four-bar mechanisms are analyzed using the above two approaches for Omega = 20 rad/s and Omega = pi rad/s and observed the same behavior in both cases. The noticeable snap-shots of the deformation shapes of the mechanism about 1000 frames are also reported using commercial software (I-DEAS + NASTRAN + ADAMS). The maximum out-of-plane warping of the cross-section is observed at the mid-span of bar-1, bar-2 and bar-3 are 1.5 mm, 250 mm and 1.0 mm, respectively, for t = 0:5 s. (C) 2015 Elsevier Ltd. All rights reserved.