Finite element analysis of interfacial friction and slip in composites with fully unbonded filler

A finite element model is developed to predict the stress-strain behaviour of particulate composites with fully unbonded filler particles. This condition can occur because of the lack of adhesion property of the filler surface. Whilst part of the filler particle is separated from the matrix, another section of filler keeps in contact with the matrix because of the lateral compressive displacement of the matrix. The slip boundary condition is imposed on the section of the interface that remains closed. The states of stress and displacement fields are obtained. The location of any further deformation through crazing or shear band formation is identified. A completely unbonded inclusion with partial slip at a section of the interface reduces the concentration of the stress at the interface significantly. Whereas this might lead to slightly higher strength, it decreases the load transfer efficiency and stiffness of this type of composite.

3D coupled thermomechanical finite element analysis of ultrasonic consolidation

A 3-D coupled temperature-displacement finite element analysis is performed to study an ultrasonic consolidation process. Results show that ultrasonic wave is effective in causing deformation in aluminum foils. Ultrasonic vibration leads to an oscillating stress field. The oscillation of stress in substrate lags behind the ultrasonic vibration by about 0.1 cycle of ultrasonic wave. The upper foil, which is in contact with the substrate, has the most severe deformation. The substrate undergoes little deformation. Apparent material softening by ultrasonic wave, which is of great concern for decades, is successfully simulated. The higher the friction coefficient, the more obvious the apparent material softening effect.

Finite element analysis of drilling of titanium alloy

Drilling is a highly demanding machining process due to complex tool geometry and the progressive material failure on the work piece. In this study, a 3D model is developed using commercial finite element software ABAQUS/Explicit. The proposed model simulates the drilling process by taking into account of the damage initiation and evolution of the work piece material, a contact model at the interface between drill bit and work piece and the process parameters. The results of the simulations demonstrate the effects of machining parameters on drilling. The results also confirm the capability and advantage of FE simulation of the drilling process. © 2011 Published by Elsevier Ltd.

Finite element analysis of ultrasonic insertion of SiC fibre in aluminium alloy 6061

In this work, a computational framework has been proposed to successfully simulate the fibre embedding using ultrasonic consolidation process. The main components of the proposed computational approach are a developed constitutive model and a friction model which are especially suitable for the condition of ultrasonic process. The effect of different process parameters, such as velocity of sonotrode, displacement amplitude of ultrasonic vibration and applied loads are studied. The presented work especially focuses on the quality of the developed weld and the fibre coverage due to the plastic flow around the fibre. The areas of maximum plastic flow predicted by the simulation are confirmed by the EBSD microstructural studies. © 2011 Inderscience Enterprises Ltd.

Cold-formed steel portal frames in fire: Full-scale testing and finite element analysis

A full-scale, non-uniform natural fire test on a cold-formed steel portal frame building is described. The results of the test are used to validate a non-linear, elasto-plastic, finite element shell idealisation, for the purposes of later forming the basis of a performance-based design approach for cold-formed steel portal frames at elevated temperatures.

Interface management for automating finite element analysis workflows

This paper outlines the importance of robust interface management for facilitating finite element analysis workflows. Topological equivalences between analysis model representations are identified and maintained in an editable and accessible manner. The model and its interfaces are automatically represented using an analysis-specific cellular decomposition of the design space. Rework of boundary conditions following changes to the design geometry or the analysis idealization can be minimized by tracking interface dependencies. Utilizing this information with the Simulation Intent specified by an analyst, automated decisions can be made to process the interface information required to rebuild analysis models. Through this work automated boundary condition application is realized within multi-component, multi-resolution and multi-fidelity analysis workflows.

Micromechanical analysis of cohesive granular materials using the discrete element method with an adhesive elasto-plastic contact model

An adhesive elasto-plastic contact model for the discrete element method with three dimensional non-spherical particles is proposed and investigated to achieve quantitative prediction of cohesive powder flowability. Simulations have been performed for uniaxial consolidation followed by unconfined compression to failure using this model. The model has been shown to be capable of predicting the experimental flow function (unconfined compressive strength vs. the prior consolidation stress) for a limestone powder which has been selected as a reference solid in the Europe wide PARDEM research network. Contact plasticity in the model is shown to affect the flowability significantly and is thus essential for producing satisfactory computations of the behaviour of a cohesive granular material. The model predicts a linear relationship between a normalized unconfined compressive strength and the product of coordination number and solid fraction. This linear relationship is in line with the Rumpf model for the tensile strength of particulate agglomerate. Even when the contact adhesion is forced to remain constant, the increasing unconfined strength arising from stress consolidation is still predicted, which has its origin in the contact plasticity leading to microstructural evolution of the coordination number. The filled porosity is predicted to increase as the contact adhesion increases. Under confined compression, the porosity reduces more gradually for the load-dependent adhesion compared to constant adhesion. It was found that the contribution of adhesive force to the limiting friction has a significant effect on the bulk unconfined strength. The results provide new insights and propose a micromechanical based measure for characterising the strength and flowability of cohesive granular materials. 

Numerical analysis of complex failure mechanisms in composite panels

A three-dimensional continuum damage mechanics-based material model has been implemented in an implicit Finite Element code to simulate the progressive degradation of advanced composite materials. The damage model uses seven damage variables assigned to tensile, compressive and non-linear shear damage at a laminae level. The objectivity of the numerical discretization is assured using a smeared formulation. The material model was benchmarked against experimental uniaxial coupon tests and it is shown to reproduce key aspects observable during failure, such as the inclined fracture plane in matrix compression and the shear band in a ±45° tension specimen.

Non-linear finite element analysis of thin composite structures

A non-conforming three-node triangular finite element with 18 degree of freedom, is used in conjugation with the Kirchhoff theory for the non-linear analysis of thin composite plate-shell structure. The formulation of the geometrically non-linear analysis is based on an updated Lagrangian formulation associated with the Newton-Raphson iterative technique, which incorporates an automatic arc-length control procedure.

Optimizing the embedded heaters position in the pultrusion process using finite element analysis

This study is based on a previous experimental work in which embedded cylindrical heaters were applied to a pultrusion machine die, and resultant energetic performance compared with that achieved with the former heating system based on planar resistances. The previous work allowed to conclude that the use of embedded resistances enhances significantly the energetic performance of pultrusion process, leading to 57% decrease of energy consumption. However, the aforementioned study was developed with basis on an existing pultrusion die, which only allowed a single relative position for the heaters. In the present work, new relative positions for the heaters were investigated in order to optimize heat distribution process and energy consumption. Finite Elements Analysis was applied as an efficient tool to identify the best relative position of the heaters into the die, taking into account the usual parameters involved in the process and the control system already tested in the previous study. The analysis was firstly developed with basis on eight cylindrical heaters located in four different location plans. In a second phase, in order to refine the results, a new approach was adopted using sixteen heaters with the same total power. Final results allow to conclude that the correct positioning of the heaters can contribute to about 10% of energy consumption reduction, decreasing the production costs and leading to a better eco-efficiency of pultrusion process.

Flexible Multibody Simulation Approach in the Dynamic Analysis of Bone Strains Activity

The objective of this study is to show that bone strains due to dynamic mechanical loading during physical activity can be analysed using the flexible multibody simulation approach. Strains within the bone tissue play a major role in bone (re)modeling. Based on previous studies, it has been shown that dynamic loading seems to be more important for bone (re)modeling than static loading. The finite element method has been used previously to assess bone strains. However, the finite element method may be limited to static analysis of bone strains due to the expensive computation required for dynamic analysis, especially for a biomechanical system consisting of several bodies. Further, in vivo implementation of strain gauges on the surfaces of bone has been used previously in order to quantify the mechanical loading environment of the skeleton. However, in vivo strain measurement requires invasive methodology, which is challenging and limited to certain regions of superficial bones only, such as the anterior surface of the tibia. In this study, an alternative numerical approach to analyzing in vivo strains, based on the flexible multibody simulation approach, is proposed. In order to investigate the reliability of the proposed approach, three 3-dimensional musculoskeletal models where the right tibia is assumed to be flexible, are used as demonstration examples. The models are employed in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model is developed using the actual geometry of the subject’s tibia, which is obtained from 3 dimensional reconstruction of Magnetic Resonance Images. Inverse dynamics simulation based on motion capture data obtained from walking at a constant velocity is used to calculate the desired contraction trajectory for each muscle. In the forward dynamics simulation, a proportional derivative servo controller is used to calculate each muscle force required to reproduce the motion, based on the desired muscle contraction trajectory obtained from the inverse dynamics simulation. Experimental measurements are used to verify the models and check the accuracy of the models in replicating the realistic mechanical loading environment measured from the walking test. The predicted strain results by the models show consistency with literature-based in vivo strain measurements. In conclusion, the non-invasive flexible multibody simulation approach may be used as a surrogate for experimental bone strain measurement, and thus be of use in detailed strain estimation of bones in different applications. Consequently, the information obtained from the present approach might be useful in clinical applications, including optimizing implant design and devising exercises to prevent bone fragility, accelerate fracture healing and reduce osteoporotic bone loss.

Development and Analysis of a Novel Isosceles Trapezoidal Dielectric Resonator Antenna For Wireless Communication

This thesis describes the development and analysis of an Isosceles Trapezoidal Dielectric Resonator Antenna (ITDRA) by realizing different DR orientations with suitable feed configurations enabling it to be used as multiband, dual band dual polarized and wideband applications. The motivation for this work has been inspired by the need for compact, high efficient, low cost antenna suitable for multi band application, dual band dual polarized operation and broadband operation with the possibility of using with MICs, and to ensure less expensive, more efficient and quality wireless communication systems. To satisfy these challenging demands a novel shaped Dielectric Resonator (DR) is fabricated and investigated for the possibility of above required properties by trying out different orientations of the DR on a simple microstrip feed and with slotted ground plane as well. The thesis initially discusses and evaluates recent and past developments taken place within the microwave industry on this topic through a concise review of literature. Then the theoretical aspects of DRA and different feeding techniques are described. Following this, fabrication and characterization of DRA is explained. To achieve the desired requirements as above both simulations and experimental measurements were undertaken. A 3-D finite element method (FEM) electromagnetic simulation tool, HFSSTM by Agilent, is used to determine the optimum geometry of the dielectric resonator. It was found to be useful in producing approximate results although it had some limitations. A numerical analysis technique, finite difference time domain (FDTD) is used for validating the results of wide band design at the end. MATLAB is used for modeling the ITDR and implementing FDTD analysis. In conclusion this work offers a new, efficient and relatively simple alternative for antennas to be used for multiple requirements in the wireless communication system.

Design and Finite Element Analysis of Hybrid Stepper Motor for Spacecraft Applications

This paper presents the design and analysis of a 400-step hybrid stepper motor for spacecraft applications. The design of the hybrid stepper motor for achieving a specific performance requires the choice of appropriate tooth geometry. In this paper, a detailed account of the results of two-dimensional finite-element (FE) analysis conducted with different tooth shapes such as square and trapezoidal, is presented. The use of % more corresponding increase in detent torque and distorted static torque profile. For the requirements of maximum torque density, less-detent torque, and better positional accuracy and smooth static torque profile, different pitch slotting with equal tooth width has to be provided. From the various FE models subjected to analysis trapezoidal teeth configuration with unequal tooth pitch on the stator and rotor is found to be the best configuration and is selected for fabrication. The designed motor is fabricated and the experimental results is compared with the FE results

Generic Finite Element Model for Mechanically Consistent Scaling of Composite Beam-to-column Joint with Welded Connection

To study the behaviour of beam-to-column composite connection more sophisticated finite element models is required, since component model has some severe limitations. In this research a generic finite element model for composite beam-to-column joint with welded connections is developed using current state of the art local modelling. Applying mechanically consistent scaling method, it can provide the constitutive relationship for a plane rectangular macro element with beam-type boundaries. Then, this defined macro element, which preserves local behaviour and allows for the transfer of five independent states between local and global models, can be implemented in high-accuracy frame analysis with the possibility of limit state checks. In order that macro element for scaling method can be used in practical manner, a generic geometry program as a new idea proposed in this study is also developed for this finite element model. With generic programming a set of global geometric variables can be input to generate a specific instance of the connection without much effort. The proposed finite element model generated by this generic programming is validated against testing results from University of Kaiserslautern. Finally, two illustrative examples for applying this macro element approach are presented. In the first example how to obtain the constitutive relationships of macro element is demonstrated. With certain assumptions for typical composite frame the constitutive relationships can be represented by bilinear laws for the macro bending and shear states that are then coupled by a two-dimensional surface law with yield and failure surfaces. In second example a scaling concept that combines sophisticated local models with a frame analysis using a macro element approach is presented as a practical application of this numerical model.

A Galerkin boundary element method for high frequency scattering by convex polygons

In this paper we consider the problem of time-harmonic acoustic scattering in two dimensions by convex polygons. Standard boundary or finite element methods for acoustic scattering problems have a computational cost that grows at least linearly as a function of the frequency of the incident wave. Here we present a novel Galerkin boundary element method, which uses an approximation space consisting of the products of plane waves with piecewise polynomials supported on a graded mesh, with smaller elements closer to the corners of the polygon. We prove that the best approximation from the approximation space requires a number of degrees of freedom to achieve a prescribed level of accuracy that grows only logarithmically as a function of the frequency. Numerical results demonstrate the same logarithmic dependence on the frequency for the Galerkin method solution. Our boundary element method is a discretization of a well-known second kind combined-layer-potential integral equation. We provide a proof that this equation and its adjoint are well-posed and equivalent to the boundary value problem in a Sobolev space setting for general Lipschitz domains.