This paper presents a 3D damage model utilized for studying the failure characteristics of GLARE. The current damage model adopted the 3D forms of Hashin’s and Puck’s failure criteria for predicting the onset of failure of fibres and matrix in the composite plies. Whilst the damage evolution is modelled based on the dissipation of fracture energy. In addition, a ductile damage model was employed to study the failure of metal layers and the delamination was assessed via a cohesive interface model. The current damage model was adopted to predict the failure modes and the blunt notch strength of GLARE; where various failure modes were observed, such as Fibre breakage, matrix cracking, delamination and plastic damage of aluminium layers. The model showed strong agreement with experimental results.
Abstract This paper presents a 3D damage model utilized for studying the failure characteristics of GLARE. The current damage model adopted the 3D forms of Hashin’s and Puck’s [...]
Investigating cure shrinkage-induced stress in thick composite beams by virtual manufacturing is the focus of this study. The research aims to understand the behaviour of thick-walled composite structures, particularly in relation to curing shrinkage-induced damages. The curing process of resin is simulated thermally and mechanically to investigate the residual cure-induced stress. The study utilizes a finite element model in Abaqus, considering material properties, mesh, boundary conditions, and user subroutines. Ten different cure cycles are investigated, showing improvements in reducing internal stresses after curing compared to the manufacturer's cycle of about 20%. However, during curing, the investigated cycles provide marginal improvements. This study demonstrates the potential for optimizing cure cycles to reduce internal stresses in thick-walled applications. It is important to note that the proposed method is not experimentally validated and requires accurate measurements for validation.
Abstract Investigating cure shrinkage-induced stress in thick composite beams by virtual manufacturing is the focus of this study. The research aims to understand the behaviour of [...]
Free vibration analysis is an essential requirement to capture the behaviour of composite structures subject to dynamic loading environment. To enhance the vibratory behaviour of composite structures, variable stiffness (VS) concept offers increased design flexibilities to tailor the structural response to meet a wide range of applications. Mechanically, the increased design space created by VS techniques leads to complexities of non-classical stiffness couplings which necessitate robust computational frameworks with enriched kinematics to predict the dynamic response accurately and efficiently. In this regard, this study proposes an enhanced differential quadrature based Strong Unified Formulation (SUF) to investigate the free vibration behaviour of thermally prestressed constant and variable stiffness composite beams. The proposed SUF model exploits the flexible kinematical description of the Theory of Unified Formulation to combine a hierarchical serendipity Lagrange-based 2D finite element (FE) with 1D differential quadrature method beam element for efficient free vibration characterisation of composite beams induced with prestress at different temperatures. The proposed SUF free vibration solutions of constant stiffness and VS beams demonstrate satisfactory accuracy and achieved improved efficiency with up to 99.9% computational savings when benchmarked against ABAQUS 3D FE solutions. Finally, a numerical study reveals that the effects of thermal prestress significantly contribute to the free vibration response of constant stiffness and VS laminated beams underscoring the importance of the study.
Abstract Free vibration analysis is an essential requirement to capture the behaviour of composite structures subject to dynamic loading environment. To enhance the vibratory behaviour [...]
Composite materials, used in primary aircraft structures, produce weight reduction and improved fuel efficiency over legacy metal airframes but are more susceptible to lightning strike damage. Therefore, research into lightning strike damage and protection systems, through experiments and simulations, is an important research topic. For any FE simulation appropriate representation of the material behaviour, the loading and boundary conditions are key to accurate predictions. In addition, an aspect which has been under reported in many studies is the meshing strategy. Fibre direction mesh alignment has been reported to yield more accurate results in the modelling of mechanical damage (intralaminar damage initiation and propagation) in unidirectional fibre reinforced composite structures. However, this model meshing strategy has not found wide application and has not been used for the modelling of thermal damage events, e.g. lightning strike direct effect simulation. Instead, authors have typically refined the mesh around the arc attachment area. This paper, for the first time, examines the influence of fibre direction mesh alignment for artificial lightning strike simulations and the prediction of thermal damage. Initially, the mesh alignment is introduced partially in the central region of the specimen. The paper uses a mature modelling approach with a transient, fully coupled, thermal-electric step in ABAQUS with a lightning test Waveform A (40 kA, 4/20 µs) applied to the specimen. Specimen boundary conditions match those typically used in experiments and a mesh convergence study is undertaken to ensure no element size influence on the results. The use of this meshing strategy has been shown to significantly improve the prediction of both moderate and severe thermal damage profiles, when compared with the standard meshes used in previous research. The predicted moderate (2659 mm2 vs 2833 mm2 ) and severe (1059 mm2 vs 1061 mm2 ) damage areas were improved to within 4% and 1% of experimental results, respectively, using this meshing strategy.
Abstract Composite materials, used in primary aircraft structures, produce weight reduction and improved fuel efficiency over legacy metal airframes but are more susceptible to lightning [...]
In this note, we are concerned with the solvability of multifield coupled problems with different, often conflictual types of non-linearities. We bring into focus the challenges of getting EFM numerical solutions. As for instance, we share our investigations of the solvability of thixoviscoplastic flow problems in FEM settings. On one hand, nonlinear multifield coupled problems are often lacking unified FEM analysis due to the presence of different nonlinearities. Thus, the importance of treating auxiliary subproblems with different analysis tools to guarantee existence of solutions. Moreover, the nonlinear multifield problems are extremely sensitive to the coupling. On other hand, monolithic Newton-multigrid FEM solver shows a great success in getting numerical solutions for multifield coupled problems. Thixoviscoplastic flow problem is a perfect example in this regard. It is a two field coupled problem, by means of microstructure dependent plastic-viscosity as well as microstructure dependent yield stress, and microstructure and shear rate dependent buildup and breakdown functions. We adapt different numerical techniques to show the solvability of the problem, and expose the accuracy of FEM numerical solutions via the simulations of thixoviscoplastic flow problems in channel configuration.
Abstract In this note, we are concerned with the solvability of multifield coupled problems with different, often conflictual types of non-linearities. We bring into focus the challenges [...]
We consider the Kirsch problem, taking into account the surface stresses at the boundary of the circular hole and on the front surfaces of the plate, in the framework of the original Gurtin–Murdoch model. The boundary conditions on the cylindrical surface of a circular hole in a nanoplate are derived in terms of a complex variable in the case of the plane stress state. The solution of the two-dimensional problem for an infinite plane with a circular hole under remote loading is explicitly obtained. Based on the analytical solution, we investigated the dependence of the elastic stress field on the nanosised plate thickness and dimension of the hole. Numerical examples are given in the paper to illustrate quantitatively the effect of the plate thickness at the nanoscale on the stress field at and near the cylindrical surface. The results are presented graphically as the dependence of the components of the stress tensor on the polar angle.
Abstract We consider the Kirsch problem, taking into account the surface stresses at the boundary of the circular hole and on the front surfaces of the plate, in the framework of the [...]
This paper focus on mathematical modelling and numerical simulation of human phonation process. The mathematical FSI model is presented consisting of the description of the structural model, the flow model and the coupling conditions. In order to treat the VFs contact, the problem of the glottis closure is addressed. To this end several ingredients are used including the use of suitable boundary conditions, modification of the flow model and robust mesh deformation algorithm. The FSI model is extended to FSAI problem by inclusion of the Lighthill model of aeroacoustics. The numerical approximation of the problem is presented and several numerical results are shown.
Abstract This paper focus on mathematical modelling and numerical simulation of human phonation process. The mathematical FSI model is presented consisting of the description of the [...]
This research aims to develop an advanced numerical model to accurately predict and optimize the acoustic insulation performance of roller shutter boxes, which are important for thermal and acoustic insulation in building facades. Traditional laboratory tests for evaluating sound transmission can be expensive and lack repeatability, particularly at low frequencies. To overcome these limitations, the proposed numerical approach utilizes the finite element method to model solid and fluid domains within the roller shutter box structure. Poroelastic layers are accounted for using a mixed displacement-pressure formulation of the Biot poroelasticity equations. Excitation and sound radiation are simulated using a diffuse field of plane waves with random phases and directions, employing the infinite elements method. The numerical model is validated by comparing its results with laboratory tests, which are described in detail. The practical application of this numerical method includes investigating factors such as assembly conditions, positioning of poroelastic layers, and the inclusion of heavy masses on the acoustic behavior of roller shutter boxes.
Abstract This research aims to develop an advanced numerical model to accurately predict and optimize the acoustic insulation performance of roller shutter boxes, which are important [...]