Stochastic free vibration analyses of composite shallow doubly curved shells – A Kriging model approach

This paper presents the Kriging model approach for stochastic free vibration analysis of composite shallow doubly curved shells. The finite element formulation is carried out considering rotary inertia and transverse shear deformation based on Mindlin’s theory. The stochastic natural frequencies are expressed in terms of Kriging surrogate models. The influence of random variation of different input parameters on the output natural frequencies is addressed. The sampling size and computational cost is reduced by employing the present method compared to direct Monte Carlo simulation. The convergence studies and error analysis are carried out to ensure the accuracy of present approach. The stochastic mode shapes and frequency response function are also depicted for a typical laminate configuration. Statistical analysis is presented to illustrate the results using Kriging model and its performance.     

Microwave assisted-in situ synthesis of porous titanium/calcium phosphate composites and their in vitro apatite-forming capability

Microwave irradiation has been proven to be an effective heating source in synthetic chemistry, and can accelerate the reaction rate, provide more uniform heating and help in developing better synthetic routes for the fabrication of bone-grafting implant materials. In this study, a new technique, which comprises microwave heating and powder metallurgy for in situ synthesis of Ti/CaP composites by using Ti powders, calcium carbonate (CaCO₃) powders and dicalcium phosphate dihydrate (CaHPO₄·2H₂O) powders, has been developed. Three different compositions of Ti:CaCO₃:CaHPO₄·2H₂O powdered mixture were employed to investigate the effect of the starting atomic ratio of the CaCO₃ to CaHPO₄·2H₂O on the phase, microstructural formation and compressive properties of the microwave synthesized composites. When the starting atomic ratio reaches 1.67, composites containing mainly alpha-titanium (α-Ti), hydroxyapatite (HA), beta-tricalcium phosphate (β-TCP) and calcium titanate (CaTiO₃) with porosity of 26%, pore size up to 152 μm, compressive strength of 212 MPa and compressive modulus of 12 GPa were formed. The in vitro apatite-forming capability of the composite was evaluated by immersing the composite into a simulated body fluid (SBF) for up to 14 days. The results showed that biodissolution occurred, followed by apatite precipitation after immersion in the SBF, suggesting that the composites are suitable for bone implant applications as apatite is an essential intermediate layer for bone cells attachment. The quantity and size of the apatite globules increased over the immersion time. After 14 days of immersion, the composite surface was fully covered by an apatite layer with a Ca/P atomic ratio approximately of 1.68, which is similar to the bone-like apatite appearing in human hard tissue. The results suggested that the microwave assisted-in situ synthesis technique can be used as an alternative to traditional powder metallurgy for the fabrication of Ti/CaP biocomposites.     

Questioning numerical integration methods for microsphere (and microplane) constitutive equations

In the last few years, more and more complex microsphere models have been proposed to predict the mechanical response of various polymers. Similarly than for microplane models, they consist in deriving a one-dimensional force vs. stretch equation and to integrate it over the unit sphere to obtain a three-dimensional constitutive equation. In this context, the focus of authors is laid on the physics of the one-dimensional relationship, but in most of the case the influence of the integration method on the prediction is not investigated.Here we compare three numerical integration schemes: a classical Gaussian scheme, a method based on a regular geometric meshing of the sphere, and an approach based on spherical harmonics. Depending on the method, the number of integration points may vary from 4 to 983,040! Considering simple quantities, i.e. principal (large) strain invariants, it is shown that the integration method must be carefully chosen. Depending on the quantities retained to described the one-dimensional equation and the required error, the performances of the three methods are discussed. Consequences on stress–strain prediction are illustrated with a directional version of the classical Mooney–Rivlin hyperelastic model. Finally, the paper closes with some advices for the development of new microsphere constitutive equations.     

Exactly satisfying initial conditions neural network models for numerical treatment of first Painlevé equation

In this paper, novel computing approach using three different models of feed-forward artificial neural networks (ANNs) are presented for the solution of initial value problem (IVP) based on first Painlevé equation. These mathematical models of ANNs are developed in an unsupervised manner with capability to satisfy the initial conditions exactly using log-sigmoid, radial basis and tan-sigmoid transfer functions in hidden layers to approximate the solution of the problem. The training of design parameters in each model is performed with sequential quadratic programming technique. The accuracy, convergence and effectiveness of the proposed schemes are evaluated on the basis of the results of statistical analyses through sufficient large number of independent runs with different number of neurons in each model as well. The comparisons of these results of proposed schemes with standard numerical and analytical solutions validate the correctness of the design models.     

Characteristics and formation mechanism of porosities in CFRP during laser joining of CFRP and steel

An experimental investigation on the mechanism of porosity formation during the laser joining of carbon fiber reinforced polymer (CFRP) and steel is presented. The porosity morphology and distribution were characterized by optical and scanning electron microscopy, and the thermal pyrolysis behaviors were investigated by thermal analysis and designed back-side cooling experiments. The results show that there are two types of porosities in CFRP. Porosity I only appears when the heat input is more than 77.8 J/mm. It has a smooth inner wall and distributes near the bonding interface between CFRP and steel at the central area of melted zone, which is caused by gaseous products such as CO₂, NH₃, H₂O, and hydrocarbons produced by the pyrolysis of CFRP. Porosity II can be seen under all joining conditions. It has a rough inner wall and distributes far away from the bonding interface, concentrating at the final solidification locations. Porosity II is caused by the shrinkage of melted CFRP during solidification stage.     

Impact behavior and fractographic study of carbon nanotubes grafted carbon fiber-reinforced epoxy matrix multi-scale hybrid composites

Carbon nanotubes are the most promising reinforcement for high performance composites. Multiwall carbon nanotubes were directly grown onto the carbon fiber surface by catalytic thermal chemical vapor deposition technique. Multi-scale hybrid composites were fabricated using the carbon nanotubes grown fibers with epoxy matrix. Morphology of the grown carbon nanotubes was investigated using field emission scanning electron microscopy and transmission electron microscopy. The fabricated composites were subjected to impact tests which showed 48.7% and 42.2% higher energy absorption in Charpy and Izod impact tests respectively. Fractographic analysis of the impact tested specimens revealed the presence of carbon nanotubes both at the fiber surface and within the matrix which explained the reason for improved energy absorption capability of these composites. Carbon nanotubes presence at various cracks formed during loading provided a direct evidence of micro crack bridging. Thus the enhanced fracture strength of these composites is attributed to stronger fiber–matrix interfacial bonding and simultaneous matrix strengthening due to the grown carbon nanotubes.     

Mechanical damage characterization of glass fiber-reinforced polymer laminates by ultrasonic maps

A method for simultaneous measurement of the thickness and density for Glass Fiber-Reinforced Polymer (GFRP) laminate plates with ultrasonic waves in C-Scan mode is presented in the form of maps. The method uses three different signals in immersion pulse-echo C-Scan mode. The maps obtained based on the density show the heterogeneity of the material at high resolution at the pixel level (1 × 1 mm²) and therefore they represent an efficient tool to assess and evaluate the damage of the composite structures after manufacturing and after an applied mechanical loading.     

Homogenization with uncertainty in Poisson ratio for polymers with rubber particles

The main aim of this work is dual computer analysis of probabilistic coefficients for the homogenized tensor of the polymer filled with the rubber particles having randomized Poisson ratios of both constituents. The major issue is to verify an influence of a randomness in rubber Poisson ratio close to the compressibility limit on the uncertainty of the effective tensor probabilistic characteristics. Probabilistic analysis presented here is carried out using mainly the stochastic perturbation technique provided by the common application of the traditional FEM commercial code ABAQUS and the symbolic computations package MAPLE. This FEM-based technique employs polynomial response function of the optimum order recovered from the weighted least squares method and following a set of deterministic solutions obtained for various values of the randomized input parameter. Optimization procedure is released entirely into a symbolic environment, where maximization of the correlation factor together with minimization of the fitting variance and approximation error are applied. Homogenization technique consists in equating of deformation energies for the real composite and the artificial one characterized by the effective elasticity tensor with uncertainty.     

Covariance matrix adaptation evolution strategy based design of fixed structure robust H∞ loop shaping controller

This paper proposes the application of Covariance Matrix Adaptation Evolution Strategy (CMA-ES) in fixed structure H_∞ loop shaping controller design. Integral Time Absolute Error (ITAE) performance requirement is incorporated as a constraint with an objective of maximization of stability margin in the fixed structure H_∞ loop shaping controller design problem. Pneumatic servo system, separating tower process and F18 fighter aircraft system are considered as test systems. The CMA-ES designed fixed structure H_∞ loop-shaping controller is compared with the traditional H_∞ loop shaping controller, non-smooth optimization and Heuristic Kalman Algorithm (HKA) based fixed structure H_∞ loop shaping controllers in terms of stability margin. 20% perturbation in the nominal plant is used to validate the robustness of the CMA-ES designed H_∞ loop shaping controller. The effect of Finite Word Length (FWL) is considered to show the implementation difficulties of controller in digital processors. Simulation results demonstrated that CMA-ES based fixed structure H_∞ loop shaping controller is suitable for real time implementation with good robust stability and performance.     

Systematic convergence of nonlinear stochastic estimators on the Special Orthogonal Group SO(3)

This paper introduces two novel nonlinear stochastic attitude estimators developed on the Special Orthogonal Group with the tracking error of the normalized Euclidean distance meeting predefined transient and steady‐state characteristics. The tracking error is confined to initially start within a predetermined large set such that the transient performance is guaranteed to obey dynamically reducing boundaries and decrease smoothly and asymptotically to the origin in probability from almost any initial condition. The proposed estimators produce accurate attitude estimates with remarkable convergence properties using measurements obtained from low‐cost inertial measurement units. The estimators proposed in continuous form are complemented by their discrete versions for the implementation purposes. The simulation results illustrate the effectiveness and robustness of the proposed estimators against uncertain measurements and large initialization error, whether in continuous or discrete form.