A finite element method for the Saint-Venant torsion and bending problems for prismatic beams

In this work, we present a finite element formulation for the Saint-Venant torsion and bending problems for prismatic beams. The torsion problem formulation is based on the warping function, and can handle multiply-connected regions (including thin-walled structures), compound and anisotropic bars. Similarly, the bending formulation, which is based on linearized elasticity theory, can handle multiply-connected domains including thin-walled sections. The torsional rigidity and shear centers can be found as special cases of these formulations. Numerical results are presented to show the good coarse-mesh accuracy of both the formulations for both the displacement and stress fields. The stiffness matrices and load vectors (which are similar to those for a variable body force in a conventional structural mechanics problem) in both formulations involve only domain integrals, which makes them simple to implement and computationally efficient. (C) 2014 Elsevier Ltd. All rights reserved.

Analytical solutions for transient temperature distribution in a geothermal reservoir due to cold water injection

An analytical solution to describe the transient temperature distribution in a geothermal reservoir in response to injection of cold water is presented. The reservoir is composed of a confined aquifer, sandwiched between rocks of different thermo-geological properties. The heat transport processes considered are advection, longitudinal conduction in the geothermal aquifer, and the conductive heat transfer to the underlying and overlying rocks of different geological properties. The one-dimensional heat transfer equation has been solved using the Laplace transform with the assumption of constant density and thermal properties of both rock and fluid. Two simple solutions are derived afterwards, first neglecting the longitudinal conductive heat transport and then heat transport to confining rocks. Results show that heat loss to the confining rock layers plays a vital role in slowing down the cooling of the reservoir. The influence of some parameters, e.g. the volumetric injection rate, the longitudinal thermal conductivity and the porosity of the porous media, on the transient heat transport phenomenon is judged by observing the variation of the transient temperature distribution with different values of the parameters. The effects of injection rate and thermal conductivity have been found to be profound on the results.

Role of different factors affecting interdiffusion in Cu(Ga) and Cu(Si) solid solutions

A detailed diffusion study was carried out on Cu(Ga) and Cu(Si) solid solutions in order to assess the role of different factors in the behaviour of the diffusing components. The faster diffusing species in the two systems, interdiffusion, intrinsic and impurity diffusion coefficients, are determined to facilitate the discussion. It was found that Cu was more mobile in the Cu-Si system, whereas Ga was the faster diffusing species in the Cu-Ga system. In both systems, the interdiffusion coefficients increased with increasing amount of solute (e.g. Si or Ga) in the matrix (Cu). Impurity diffusion coefficients for Si and Ga in Cu, found out by extrapolating interdiffusion coefficient data to zero composition of the solute, were both higher than the Cu tracer diffusion coefficient. These observed trends in diffusion behaviour could be rationalized by considering: (i) formation energies and concentration of vacancies, (ii) elastic moduli (indicating bond strengths) of the elements and (iii) the interaction parameters and the related thermodynamic factors. In summary, we have shown here that all the factors introduced in this paper should be considered simultaneously to understand interdiffusion in solid solutions. Otherwise, some of the aspects may look unusual or even impossible to explain.

Comparative Study of Protein Unfolding in Aqueous Urea and Dimethyl Sulfoxide Solutions: Surface Polarity, Solvent Specificity, and Sequence of Secondary Structure Melting

Elucidation of possible pathways between folded (native) and unfolded states of a protein is a challenging task, as the intermediates are often hard to detect. Here, we alter the solvent environment in a controlled manner by choosing two different cosolvents of water, urea, and dimethyl sulfoxide (DMSO) and study unfolding of four different proteins to understand the respective sequence of melting by computer simulation methods. We indeed find interesting differences in the sequence of melting of alpha helices and beta sheets in these two solvents. For example, in 8 M urea solution, beta-sheet parts of a protein are found to unfold preferentially, followed by the unfolding of alpha helices. In contrast, 8 M DMSO solution unfolds alpha helices first, followed by the separation of beta sheets for the majority of proteins. Sequence of unfolding events in four different alpha/beta proteins and also in chicken villin head piece (HP-36) both in urea and DMSO solutions demonstrate that the unfolding pathways are determined jointly by relative exposure of polar and nonpolar residues of a protein and the mode of molecular action of a solvent on that protein.

Freeform Skeletal Shape Optimization of Compliant Mechanisms

Compliant mechanisms are elastic continua used to transmit or transform force and motion mechanically. The topology optimization methods developed for compliant mechanisms also give the shape for a chosen parameterization of the design domain with a fixed mesh. However, in these methods, the shapes of the flexible segments in the resulting optimal solutions are restricted either by the type or the resolution of the design parameterization. This limitation is overcome in this paper by focusing on optimizing the skeletal shape of the compliant segments in a given topology. It is accomplished by identifying such segments in the topology and representing them using Bezier curves. The vertices of the Bezier control polygon are used to parameterize the shape-design space. Uniform parameter steps of the Bezier curves naturally enable adaptive finite element discretization of the segments as their shapes change. Practical constraints such as avoiding intersections with other segments, self-intersections, and restrictions on the available space and material, are incorporated into the formulation. A multi-criteria function from our prior work is used as the objective. Analytical sensitivity analysis for the objective and constraints is presented and is used in the numerical optimization. Examples are included to illustrate the shape optimization method.

Phase transformation, improved ferroelectric and magnetic properties of (1-x) BiFeO3-xPb(Zr0.52Ti0.48)O-3 solid solutions

The authors prepared (1 - x) BiFeO3 - (x)Pb(Zr0.52Ti0.48)O-3 for x <= 0.30 by sol-gel method and investigated the material's structures, magnetic and electrical properties. Detailed Rietveld analysis of X-ray diffraction data revealed that the system retains distorted rhombohedral R3c structure for x <= 0.10 but transforms to monoclinic (Cc) structure for x > 0.10. Disappearance of some Raman modes corresponding to A1 modes and the decrease in the intensities of the remaining A1 modes with increasing x in the Raman spectra, which is a clear indication of structural modification and symmetry changes brought about by PZT doping. Enhanced magnetization with PZT doping content may be attributed to the gradual change and destruction in the spin cycloid structure of BiFeO3. The leakage current density at 3.5 kV/cm was reduced by approximately three orders of magnitude by doping PZT (x = 0.30), compared with BFO ceramics. (C) 2014 AIP Publishing LLC.

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.

Wave propagation analysis in laminated composite plates with transverse cracks using the wavelet spectral finite element method

This paper presents a newly developed wavelet spectral finite element (WFSE) model to analyze wave propagation in anisotropic composite laminate with a transverse surface crack penetrating part-through the thickness. The WSFE formulation of the composite laminate, which is based on the first-order shear deformation theory, produces accurate and computationally efficient results for high frequency wave motion. Transverse crack is modeled in wavenumber-frequency domain by introducing bending flexibility of the plate along crack edge. Results for tone burst and impulse excitations show excellent agreement with conventional finite element analysis in Abaqus (R). Problems with multiple cracks are modeled by assembling a number of spectral elements with cracks in frequency-wavenumber domain. Results show partial reflection of the excited wave due to crack at time instances consistent with crack locations. (C) 2014 Elsevier B.V. All rights reserved.

Maximum mass of stable magnetized highly super-Chandrasekhar white dwarfs: stable solutions with varying magnetic fields

We address the issue of stability of recently proposed significantly super-Chandrasekhar white dwarfs. We present stable solutions of magnetostatic equilibrium models for super-Chandrasekhar white dwarfs pertaining to various magnetic field profiles. This has been obtained by self-consistently including the effects of the magnetic pressure gradient and total magnetic density in a general relativistic framework. We estimate that the maximum stable mass of magnetized white dwarfs could be more than 3 solar mass. This is very useful to explain peculiar, overluminous type Ia supernovae which do not conform to the traditional Chandrasekhar mass-limit.

Evaluating FuncSION: A software for automated synthesis of design solutions for stimulating ideation during mechanical conceptual design

The goal of the work reported in this paper is to use automated, combinatorial synthesis to generate alternative solutions to be used as stimuli by designers for ideation. FuncSION, a computational synthesis tool that can automatically synthesize solution concepts for mechanical devices by combining building blocks from a library, is used for this purpose. The objectives of FuncSION are to help generate a variety of functional requirements for a given problem and a variety of concepts to fulfill these functions. A distinctive feature of FuncSION is its focus on automated generation of spatial configurations, an aspect rarely addressed by other computational synthesis programs. This paper provides an overview of FuncSION in terms of representation of design problems, representation of building blocks, and rules with which building blocks are combined to generate concepts at three levels of abstraction: topological, spatial, and physical. The paper then provides a detailed account of evaluating FuncSION for its effectiveness in providing stimuli for enhanced ideation.

Finite element analysis of free vibration of the delaminated composite plate with variable kinematic multilayered plate elements

Composite laminates are prone to delamination. Implementation of delamination in the Carrera Unified Formulation frame work using nine noded quadrilateral MITC9 element is discussed in this article. MITC9 element is devoid of shear locking and membrane locking. Delaminated as well as healthy structure is analyzed for free mode vibration. The results from the present work are compared with the available experimental or/and research article or/and the three dimensional finite element simulations. The effect of different kinds and different percentages of area of delamination on the first three natural frequencies of the structure is discussed. The presence of open-mode delamination mode shape for large delaminations within the first three natural frequencies is discussed. Also, the switching of places between the second bending mode, with that of the first torsional mode frequency is discussed. Results obtained from different ordered theories are compared in the presence of delamination. Advantage of layerwise theories as compared to equivalent single layer theories for very large delaminations is stated. The effect of different kinds of delamination and their effect on the second bending and first torsional mode shape is discussed. (C) 2014 Elsevier Ltd. 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.

Mixed finite elements for electromagnetic analysis

The occurrence of spurious solutions is a well-known limitation of the standard nodal finite element method when applied to electromagnetic problems. The two commonly used remedies that are used to address this problem are (i) The addition of a penalty term with the penalty factor based on the local dielectric constant, and which reduces to a Helmholtz form on homogeneous domains (regularized formulation); (ii) A formulation based on a vector and a scalar potential. Both these strategies have some shortcomings. The penalty method does not completely get rid of the spurious modes, and both methods are incapable of predicting singular eigenvalues in non-convex domains. Some non-zero spurious eigenvalues are also predicted by these methods on non-convex domains. In this work, we develop mixed finite element formulations which predict the eigenfrequencies (including their multiplicities) accurately, even for nonconvex domains. The main feature of the proposed mixed finite element formulation is that no ad-hoc terms are added to the formulation as in the penalty formulation, and the improvement is achieved purely by an appropriate choice of finite element spaces for the different variables. We show that the formulation works even for inhomogeneous domains where `double noding' is used to enforce the appropriate continuity requirements at an interface. For two-dimensional problems, the shape of the domain can be arbitrary, while for the three-dimensional ones, with our current formulation, only regular domains (which can be nonconvex) can be modeled. Since eigenfrequencies are modeled accurately, these elements also yield accurate results for driven problems. (C) 2014 Elsevier Ltd. All rights reserved.

Finite element simulations of notch tip fields in magnesium single crystals

Recent experiments using three point bend specimens of Mg single crystals have revealed that tensile twins of {10 (1) over bar2}-type form profusely near a notch tip and enhance the fracture toughness through large plastic dissipation. In this work, 3D finite element simulations of these experiments are carried out using a crystal plasticity framework which includes slip and twinning to gain insights on the mechanics of fracture. The predicted load-displacement curves, slip and tensile twinning activities from finite element analysis corroborate well with the experimental observations. The numerical results are used to explore the 3D nature of the crack tip stress, plastic slip and twin volume fraction distributions near the notch root. The occurrence of tensile twinning is rationalized from the variation of normal stress ahead of the notch tip. Further, deflection of the crack path at twin-twin intersections observed in the experiments is examined from an energy standpoint by modeling discrete twins close to the notch root.

Large deformation dynamic finite element analysis of delaminated composite plates using contact-impact conditions

A new C-0 composite plate finite element based on Reddy's third order theory is used for large deformation dynamic analysis of delaminated composite plates. The inter-laminar contact is modeled with an augmented Lagrangian approach. Numerical results show that the widely used ``unconditionally stable'' beta-Newmark method presents instability problems in the transient simulation of delaminated composite plate structures with large deformation. To overcome this instability issue, an energy and momentum conserving composite implicit time integration scheme presented by Bathe and Baig is used. It is found that a proper selection of the penalty parameter is very crucial in the contact simulation. (C) 2014 Elsevier Ltd. All rights reserved.