Analysis of adhesively bonded CFRE composite scarf joints modified with MWCNTs

The main objective of the present work is to improve the performance of bonded joints in carbon fiber composite structures through introducing Multi-Walled Carbon Nanotubes (MWCNTs) into Epocast 50-A1/946 epoxy, which was primarily developed for joining and repairing of composite aircraft structures. Results from tension characterizations of structural adhesive joints (SAJs) with different scarf angles (5–45°) showed improvement up to 40% compared to neat epoxy (NE)–SAJs. Special attention was considered to investigate the performance of SAJs with 5° scarf angle under different environments. The tensile strength and stiffness of both NE-SAJs and MWCNT/E-SAJs were dramatically decreased at elevated temperature. Water absorption showed a marginal drop of about 2.0% in the tensile strength of the moist SAJs compared to the dry one. Cracks initiation and propagation were detected effectively using instrumented-SAJs with eight strain gauges. The experimental results agree well with the predicted using three-dimensional finite element analysis model.     

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.     

Design criterion for fatigue strengthening of riveted beams in a 120-year-old railway metallic bridge using pre-stressed CFRP plates

This study presents a design criterion developed for fatigue strengthening of a 120-year-old metallic railway bridge in Switzerland and presents a pre-stressed un-bonded reinforcement (PUR) system developed to apply the strengthening. The PUR system uses carbon fiber reinforced polymer (CFRP) plates; however, unlike conventional pre-stressed CFRP reinforcement methods, preparation of the existing metallic bridge surface is not required. This decreases the time required for on-site strengthening procedures. The principle of the constant life diagram (CLD) and two fatigue failure criteria (Johnson and Goodman) are described. Analytical formulations are developed based on the CLD method to determine the minimum CFRP pre-stress level required to prevent fatigue crack initiation. The PUR system uses an applied pre-stress force to reduce the mean stress level (and stress ratio) to shift an existing fatigue-susceptible metallic detail from the ‘at risk’ finite life regime to the ‘safe’ infinite life regime. The applied CLD method is particularly valuable when the stress history of the detail is not known and it is difficult to assess the remaining fatigue life. Moreover, it is shown that the currently adopted approach in many structural codes which emphasizes stress range as the dominant parameter influencing fatigue life are non-conservative for tension–tension stress patterns (i.e., stress ratios of 0 < R < 1). Analyses show that the modified Johnson formula accurately reflects the combined effect of stress range, mean stress level, and material properties, and offers a relatively easy design procedure. Details of a retrofit field application on members of a riveted wrought iron railway bridge are given. A wireless sensor network (WSN) system is used for long-term monitoring of the on-site CFRP stress levels and temperature of the retrofitted details. WSN measurements indicate that increases in ambient temperature result in increased CFRP pre-stress levels.     

Nonlinear response in glass fibre non-crimp fabric reinforced vinylester composites

This paper addresses the nonlinear stress-strain response in glass fibre non-crimp fabric reinforced vinylester composite laminates subjected to in-plane tensile loading. The nonlinearity is shown to be a combination of brittle and plastic failure. It is argued that the shift from plastic to brittle behaviour in the vinylester is caused by the state of stress triaxiality caused by the interaction between fibre and vinylester. A model combining damage and plasticity is calibrated and evaluated using data from extensive experimental testing. The onset of damage is predicted using the Puck failure criterion, and the evolution of damage is calibrated from the observed softening in plies loaded in transverse tension. Shear loading beyond linear elastic response is observed to result in irreversible strains. A yield criterion is implemented for shear deformation. A strain hardening law is fitted to the stress-strain response observed in shear loaded plies. Experimental results from a selection of laminates with different layups are used to verify the numerical models. A complete set of model parameters for predicting elastic behaviour, strength and post failure softening is presented for glass fibre non-crimped fabric reinforced vinylester. The predicted behaviour from using these model parameters are shown to be in good agreement with experimental results.     

A high performance crashworthiness simulation system based on GPU

Crashworthiness simulation system is one of the key computer-aided engineering (CAE) tools for the automobile industry and implies two potential conflicting requirements: accuracy and efficiency. A parallel crashworthiness simulation system based on graphics processing unit (GPU) architecture and the explicit finite element (FE) method is developed in this work. Implementation details with compute unified device architecture (CUDA) are considered. The entire parallel simulation system involves a parallel hierarchy-territory contact-searching algorithm (HITA) and a parallel penalty contact force calculation algorithm. Three basic GPU-based parallel strategies are suggested to meet the natural parallelism of the explicit FE algorithm. Two free GPU-based numerical calculation libraries, cuBLAS and Thrust, are introduced to decrease the difficulty of programming. Furthermore, a mixed array and a thread map to element strategy are proposed to improve the performance of the test pairs searching. The outer loop of the nested loop through the mixed array is unrolled to realize parallel searching. An efficient storage strategy based on data sorting is presented to realize data transfer between different hierarchies with coalesced access during the contact pairs searching. A thread map to element pattern is implemented to calculate the penetrations and the penetration forces; a double float atomic operation is used to scatter contact forces. The simulation results of the three different models based on the Intel Core i7-930 and the NVIDIA GeForce GTX 580 demonstrate the precision and efficiency of this developed parallel crashworthiness simulation system.     

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.     

Pull out of FRP reinforcement from masonry pillars: Experimental and numerical results

In this paper, debonding phenomena between carbon fiber reinforced polymer (CFRP) strips and masonry support were investigated on the basis of single-lap shear tests, considering different dimensions of the bond length. To capture the post-peak response of the CFRP–masonry joint, the slip between the support and the reinforcement strip was controlled using a clip gauge positioned at the end of the reinforcement. The tests were simulated by means of a finite element model able to capture the post-peak snap-back behavior due to the failure process. The numerical model is based on zero-thickness interface elements and on a proper non-linear cohesive law. The comparison between experimental and numerical results was performed in terms of overall response, measured by both the machine stroke and the clip gauge positioned at the free end of the reinforcement. The cases of effective bond length greater and lesser than the minimum anchorage length, suggested by the CNR Italian recommendation, were considered.     

A comparison of transformation methods for evaluating two‐dimensional weakly singular integrals

Accurate numerical evaluation of integrals arising in the boundary element method is fundamental to achieving useful results via this solution technique. In this paper, a number of techniques are considered to evaluate the weakly singular integrals which arise in the solution of Laplace's equation in three dimensions and Poisson's equation in two dimensions. Both are two‐dimensional weakly singular integrals and are evaluated using (in a product fashion) methods which have recently been used for evaluating one‐dimensional weakly singular integrals arising in the boundary element method. The methods used are based on various polynomial transformations of conventional Gaussian quadrature points where the transformation polynomial has zero Jacobian at the singular point. Methods which split the region of integration into sub‐regions are considered as well as non‐splitting methods. In particular, the newly introduced and highly accurate generalized composite subtraction of singularity and non‐linear transformation approach (GSSNT) is applied to various two‐dimensional weakly singular integrals. A study of the different methods reveals complex relationships between transformation orders, position of the singular point, integration kernel and basis function. It is concluded that the GSSNT method gives the best overall results for the two‐dimensional weakly singular integrals studied. Copyright © 2002 John Wiley & Sons, Ltd.     

Prediction of brittle-to-ductile transitions in polystyrene

In this study it is attempted to predict brittle-to-ductile transitions (BDTs) in polystyrene blends, induced either by an increase in temperature or by a decrease in inter-particle distance. A representative, two-dimensional volume element (RVE) of a polystyrene matrix with 20% circular voids, is deformed in tension. During deformation a hydrostatic-stress based craze-nucleation criterion [1] is evaluated. The simulations demonstrate that crazes initiate at low temperatures while a transition from crazing to shear yielding (BDT) is found around 75 °C. The numerical results correlate well with tensile tests on similar heterogeneous polystyrene. The presence of an absolute length, as experimentally found, is more difficult to explain. Near a free surface a T_g-depression is measured for polystyrene and also the resistance to indentation in polystyrene is lower than expected from bulk properties. Both observations are rationalised by an enhanced segmental mobility of chains near a free surface. As a consequence of these findings, an absolute length-scale could be incorporated in the numerical simulations. For simplicity, the length-scale is modelled by taking a temperature gradient over a thin layer near the internal free surfaces of the RVE. Deformation of the RVE with different absolute length-scales shows that indeed also the experimentally found brittle-to-ductile transition can be predicted if the ligament thickness between the inclusions (‘voids’) in polystyrene is below a critical value of ca. 15 nm.