Alternative positional FEM applied to thermomechanical impact of truss structures

This paper presents a formulation to deal with dynamic thermomechanical problems by the finite element method. The proposed methodology is based on the minimum potential energy theorem written regarding nodal positions, not displacements, to solve the mechanical problem. The thermal problem is solved by a regular finite element method. Such formulation has the advantage of being simple and accurate. As a solution strategy, it has been used as a natural split of the thermomechanical problem, usually called isothermal split or isothermal staggered algorithm. Usual internal variables and the additive decomposition of the strain tensor have been adopted to model the plastic behavior. Four examples are presented to show the applicability of the technique. The results are compared with other authors` numerical solutions and experimental results. (C) 2010 Elsevier B.V. All rights reserved.

A three dimensional finite volume approach to the thermo-mechanical modelling of the shape casting of metals

This paper presents a three dimensional, thermos-mechanical modelling approach to the cooling and solidification phases associated with the shape casting of metals ei. Die, sand and investment casting. Novel vortex-based Finite Volume (FV) methods are described and employed with regard to the small strain, non-linear Computational Solid Mechanics (CSM) capabilities required to model shape casting. The CSM capabilities include the non-linear material phenomena of creep and thermo-elasto-visco-plasticity at high temperatures and thermo-elasto-visco-plasticity at low temperatures and also multi body deformable contact with which can occur between the metal casting of the mould. The vortex-based FV methods, which can be readily applied to unstructured meshes, are included within a comprehensive FV modelling framework, PHYSICA. The additional heat transfer, by conduction and convection, filling, porosity and solidification algorithms existing within PHYSICA for the complete modelling of all shape casting process employ cell-centred FV methods. The termo-mechanical coupling is performed in a staggered incremental fashion, which addresses the possible gap formation between the component and the mould, and is ultimately validated against a variety of shape casting benchmarks.

Chaos and Unpredictability in the Vibration of an Elasto-Plastic Beam

This contribution discusses the nonlinear dynamics of a pin-ended elasto-plastic beam with both kinematic and isotropic hardening. An iterative numerical procedure based on the operator split technique is developed in order to deal with the nonlinearities in the equations of motion. Free and forced responses for harmonic sinusoidal and square wave excitations are investigated. Numerical simulations present many interesting behaviors such as jump phenomena, sensitivity to initial conditions, chaos and transient chaos. These results indicate that there are practical problems in predicting the response of the beam even when periodic steady state response is expected.

Hygro-elasto-plastic behavior of planar orthotropic material

The mechanical and hygroscopic properties of paper and board are factors affecting the whole lifecycle of a product, including paper/board quality, production, converting, and material and energy savings. The progress of shrinkage profiles, loose edges of web, baggy web causing wrinkling and misregistration in printing are examples of factors affecting runnability and end product quality in the drying section and converting processes, where paper or board is treated as a moving web. The structural properties and internal stresses or plastic strain differences built up during production also cause the end-product defects related to distortion of the shape of the product such as sheet or box. The objective of this work was to construct a model capable of capturing the characteristic behavior of hygroscopic orthotropic material under moisture change, during different external in-plane stretch or stress conditions. Two independent experimental models were constructed: the elasto-plastic material model and the hygroexpansivity-shrinkage model. Both describe the structural properties of the sheet with a fiber orientation probability distribution, and both are functions of the dry solids content and fiber orientation anisotropy index. The anisotropy index, introduced in this work, simplifies the procedure of determining the constitutive parameters of the material model and the hygroexpansion coefficients in different in-plane directions of the orthotropic sheet. The mathematically consistent elasto-plastic material model and the dry solids content dependent hygroexpansivity have been constructed over the entire range from wet to dry. The presented elastoplastic and hygroexpansivity-shrinkage models can be used in an analytical approach to estimate the plastic strain and shrinkage in simple one-dimensional cases. For studies of the combined and more complicated effects of hygro-elasto-plastic behavior, both models were implemented in a finite element program for a numerical solution. The finite element approach also offered possibilities for studying different structural variations of orthotropic planar material, as well as local buckling behavior and internal stress situations of the sheet or web generated by local strain differences. A comparison of the simulation examples presented in this work to results published earlier confirms that the hygro-elasto-plastic model provides at least qualitatively reasonable estimates. The application potential of the hygro-elasto-plastic model is versatile, including several phenomena and defects appearing in the drying, converting and end-use conditions of the paper or board webs and products, or in other corresponding complex planar materials.

Plastic-damage smeared crack model to simulate the behaviour of structures made by cement based materials

This work proposes a constitutive model to simulate nonlinear behaviour of cement based materials subjected to different loading paths. The model incorporates a multidirectional fixed smeared crack approach to simulate crack initiation and propagation, whereas the inelastic behaviour of material between cracks is treated by a numerical strategy that combines plasticity and damage theories. For capturing more realistically the shear stress transfer between the crack surfaces, a softening diagram is assumed for modelling the crack shear stress versus crack shear strain. The plastic damage model is based on the yield function, flow rule and evolution law for hardening variable, and includes an explicit isotropic damage law to simulate the stiffness degradation and the softening behaviour of cement based materials in compression. This model was implemented into the FEMIX computer program, and experimental tests at material scale were simulated to appraise the predictive performance of this constitutive model. The applicability of the model for simulating the behaviour of reinforced concrete shear wall panels submitted to biaxial loading conditions, and RC beams failing in shear is investigated.

Application of the finite element method in cold forging processes

The demand for more efficient manufacturing processes has been increasing in the last few years. The cold forging process is presented as a possible solution, because it allows the production of parts with a good surface finish and with good mechanical properties. Nevertheless, the cold forming sequence design is very empirical and it is based on the designer experience. The computational modeling of each forming process stage by the finite element method can make the sequence design faster and more efficient, decreasing the use of conventional "trial and error" methods. In this study, the application of a commercial general finite element software - ANSYS - has been applied to model a forming operation. Models have been developed to simulate the ring compression test and to simulate a basic forming operation (upsetting) that is applied in most of the cold forging parts sequences. The simulated upsetting operation is one stage of the automotive starter parts manufacturing process. Experiments have been done to obtain the stress-strain material curve, the material flow during the simulated stage, and the required forming force. These experiments provided results used as numerical model input data and as validation of model results. The comparison between experiments and numerical results confirms the developed methodology potential on die filling prediction.

A generalized anisotropic quadric yield criterion and its application to bone tissue at multiple length scales

Nonlinear computational analysis of materials showing elasto-plasticity or damage relies on knowledge of their yield behavior and strengths under complex stress states. In this work, a generalized anisotropic quadric yield criterion is proposed that is homogeneous of degree one and takes a convex quadric shape with a smooth transition from ellipsoidal to cylindrical or conical surfaces. If in the case of material identification, the shape of the yield function is not known a priori, a minimization using the quadric criterion will result in the optimal shape among the convex quadrics. The convexity limits of the criterion and the transition points between the different shapes are identified. Several special cases of the criterion for distinct material symmetries such as isotropy, cubic symmetry, fabric-based orthotropy and general orthotropy are presented and discussed. The generality of the formulation is demonstrated by showing its degeneration to several classical yield surfaces like the von Mises, Drucker–Prager, Tsai–Wu, Liu, generalized Hill and classical Hill criteria under appropriate conditions. Applicability of the formulation for micromechanical analyses was shown by transformation of a criterion for porous cohesive-frictional materials by Maghous et al. In order to demonstrate the advantages of the generalized formulation, bone is chosen as an example material, since it features yield envelopes with different shapes depending on the considered length scale. A fabric- and density-based quadric criterion for the description of homogenized material behavior of trabecular bone is identified from uniaxial, multiaxial and torsional experimental data. Also, a fabric- and density-based Tsai–Wu yield criterion for homogenized trabecular bone from in silico data is converted to an equivalent quadric criterion by introduction of a transformation of the interaction parameters. Finally, a quadric yield criterion for lamellar bone at the microscale is identified from a nanoindentation study reported in the literature, thus demonstrating the applicability of the generalized formulation to the description of the yield envelope of bone at multiple length scales.

Elasto-plastic structural stability of arches

The elastic-plastic structural stability behaviour of arches is analysed in the present work.The application of the developed mathematical model, allows to determine the elastic-plastic equilibrium paths, looking for critical points, bifurcation or limit, along those paths, associated to the critical load, in case it comes to happen.The equilibrium paths in the elastic-plastic behaviour when compared with the paths in the linear elastic behaviour, may show that, due to influence of the material plasticity, modifications paths appear and consequently alterations in the values of its critical loads.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

Impact Of Pregabalin Treatment On Synaptic Plasticity And Glial Reactivity During The Course Of Experimental Autoimmune Encephalomyelitis.

Multiple sclerosis (MS) is an autoimmune and neurodegenerative disease that affects young adults. It is characterized by generating a chronic demyelinating autoimmune inflammation in the central nervous system. An experimental model for studying MS is the experimental autoimmune encephalomyelitis (EAE), induced by immunization with antigenic proteins from myelin. The present study investigated the evolution of EAE in pregabalin treated animals up to the remission phase. The results demonstrated a delay in the onset of the disease with statistical differences at the 10th and the 16th day after immunization. Additionally, the walking track test (CatWalk) was used to evaluate different parameters related to motor function. Although no difference between groups was obtained for the foot print pressure, the regularity index was improved post treatment, indicating a better motor coordination. The immunohistochemical analysis of putative synapse preservation and glial reactivity revealed that pregabalin treatment improved the overall morphology of the spinal cord. A preservation of circuits was depicted and the glial reaction was downregulated during the course of the disease. qRT-PCR data did not show immunomodulatory effects of pregabalin, indicating that the positive effects were restricted to the CNS environment. Overall, the present data indicate that pregabalin is efficient for reducing the seriousness of EAE, delaying its course as well as reducing synaptic loss and astroglial reaction.

Comparing Near-infrared Conventional Diffuse Reflectance Spectroscopy And Hyperspectral Imaging For Determination Of The Bulk Properties Of Solid Samples By Multivariate Regression: Determination Of Mooney Viscosity And Plasticity Indices Of Natural Rubber.

Conventional reflectance spectroscopy (NIRS) and hyperspectral imaging (HI) in the near-infrared region (1000-2500 nm) are evaluated and compared, using, as the case study, the determination of relevant properties related to the quality of natural rubber. Mooney viscosity (MV) and plasticity indices (PI) (PI0 - original plasticity, PI30 - plasticity after accelerated aging, and PRI - the plasticity retention index after accelerated aging) of rubber were determined using multivariate regression models. Two hundred and eighty six samples of rubber were measured using conventional and hyperspectral near-infrared imaging reflectance instruments in the range of 1000-2500 nm. The sample set was split into regression (n = 191) and external validation (n = 95) sub-sets. Three instruments were employed for data acquisition: a line scanning hyperspectral camera and two conventional FT-NIR spectrometers. Sample heterogeneity was evaluated using hyperspectral images obtained with a resolution of 150 × 150 μm and principal component analysis. The probed sample area (5 cm(2); 24,000 pixels) to achieve representativeness was found to be equivalent to the average of 6 spectra for a 1 cm diameter probing circular window of one FT-NIR instrument. The other spectrophotometer can probe the whole sample in only one measurement. The results show that the rubber properties can be determined with very similar accuracy and precision by Partial Least Square (PLS) regression models regardless of whether HI-NIR or conventional FT-NIR produce the spectral datasets. The best Root Mean Square Errors of Prediction (RMSEPs) of external validation for MV, PI0, PI30, and PRI were 4.3, 1.8, 3.4, and 5.3%, respectively. Though the quantitative results provided by the three instruments can be considered equivalent, the hyperspectral imaging instrument presents a number of advantages, being about 6 times faster than conventional bulk spectrometers, producing robust spectral data by ensuring sample representativeness, and minimizing the effect of the presence of contaminants.

Arachnoid cyst: adversity ans plasticity

We report four cases of surgically treated intracranial arachnoid cysts, one with cyst-peritoneal shunt and three with craniotomy and arachnoid membrane resection. Their classification and etiopathogeny are discussed, and especially the different methods of treatment comparing the drastic complications (adversities) with the favorable solutions in severe clinical cases (plasticity) treated at our institution.

Thermo-optical characteristics and concentration quenching effects in Nd(3+)doped yttrium calcium borate glasses

In this work we present a comprehensive study of the spectroscopic and thermo-optical properties of a set of samples with composition xNd(2)O(3)-(5-x)Y(2)O(3-)40CaO-55B(2)O(3) (0 <= x <= 1.0 mol%). Their fluorescence quantum efficiency (eta) values were determined using the thermal lens technique and the dependence on the ionic concentration was analyzed in terms of energy transfer processes, based on the Forster-Dexter model of multipolar ion-ion interactions. A maximum eta = 0.54 was found to be substantially higher than for yttrium aluminoborate crystals and glasses with comparable Nd(3+) content. As for the thermo-optical properties of yttrium calcium borate, they are comparable to other well-known laser glasses. The obtained energy transfer microparameters and the weak dependence of. on the Nd(3+) concentration with a high optimum Nd(3+) concentration put this system as a strong candidate for photonics applications. (C) 2011 American Institute of Physics. [doi:10.1063/1.3567091]

Luminescent and thermo-optical properties of Nd(3+)-doped yttrium aluminoborate laser glasses

In this work we performed a thorough spectroscopic and thermo-optical investigation of yttrium aluminoborate glasses doped with neodymium ions. A set of samples, prepared by the conventional melt-quenching technique and with Nd(2)O(3) concentrations varying from 0.1 to 0.75 mol %, were characterized by ground state absorption, photoluminescence, excited state lifetime measurements, and thermal lens technique. For the neodymium emission at 1064 nm ((4)F(3/2) -> (4)I(11/2) transition), no significant luminescence concentration quenching was observed and the experimental lifetime values ranged around 70 mu s. The obtained values of thermal conductivity and diffusivity of approximately 10.3 x 10(-3) W / cm K and 4.0 x 10(-3) cm(2) / s, respectively, are comparable to those of commercial laser glasses. Moreover, the fluorescence quantum efficiency of the glasses, calculated using the Judd-Ofelt formalism and luminescence decay, lies in the range from 0.28 to 0.32, larger than the typical values obtained for Nd(3+) doped YAl(3)(BO(3))(4) crystals. (c) 2009 American Institute of Physics. [DOI: 10.1063/1.3176503]

Thermal lens and interferometric method for glass transition and thermo physical properties measurements in Nd(2)O(3) doped sodium zincborate glass

In this work the time resolved thermal lens method is combined with interferometric technique, the thermal relaxation calorimetry, photoluminescence and lifetime measurements to determine the thermo physical properties of Nd(2)O(3) doped sodium zincborate glass as a function of temperature up to the glass transition region. Thermal diffusivity, thermal conductivity, fluorescence quantum efficiency, linear thermal expansion coefficient and thermal coefficient of electronic polarizability were determined. In conclusion, the results showed the ability of thermal lens and interferometric methods to perform measurements very close to the phase transition region. These techniques provide absolute values for the measured physical quantities and are advantageous when low scan rates are required. (c) 2008 Optical Society of America