Unmanned Aerial Systems (UASs) are increasingly starting to dominate the lower airspace. This increases the chance that a UAS will hit the means of transport of people e.g. aircraft, helicopters. For air traffic, the European Union Aviation Safety Agency (EASA) and the Federal Aviation Administration (FAA) are in the process of determining the threat of a UAS impact on aircraft. This will result in new rules that may impose new, or additional, requirements on the ‘wetted zones’ of aircraft and the like. Current research suggests that aircraft wetted-areas (e.g. wing leading edge) that are certified for so called ‘bird-impact’ may not sustain a ‘UAS-impact’, such an UAS-impact may even damage the primary load carrying structure. But what would happen if multiple UASs are flying close to one other? To the authors’ knowledge, the damage caused by multiple UAS impacts on a wetted surface zone has not yet been established.

A finite element modelling approach is chosen for the UAS; specifically, a Lagrangian approach is applied using material nonlinearity and damage. A comparison is made between the damage caused by a bird impactor and a UAS impactor model. To establish the resulting damage of multiple UAS impacts on a wing leading edge, a multiple-UAS impact scenario is executed. The results show that a wing leading edge capable of sustaining a bird impact may not be capable of sustaining a UAS impact, which supports previous findings. Furthermore, for all simulated cases the front spar was not penetrated by components that managed to enter the leading edge. However, for the heavier drone some deformation of the front spar was observed. The multiple UAS impact scenario causes additional damage to the leading edge with respect to the single UAS impact, with greater deformation of the front spar being observed.

In this paper, the influence of strain rate on the mechanical behavior of high-strength low-alloy (HC420LA) steel were studied. Quasi-static and dynamic tensile experiments were performed with strain rates ranging from 0.001 to 500 s−1 at room temperature. The digital image correlation technique was used to obtain the full-field strain. The experimental results showed that HC420LA steel exhibited positive strain rate sensitivity. Based on experimental results, the modified Johnson–Cook (J–C) model was used to model the constitutive behavior of HC420LA steel. Predictions of the standard and modified J–C models were compared using standard statistical parameters. The modified J–C model showed better agreement with the experimental data. Then, numerical simulation of the representative tensile test at a strain rate of 100 s−1 was performed using the finite element code LS-DYNA. Good correlation between the experimental and numerical simulation results was achieved.

Full-scale simulations of a (Magnetorheological) MR damper are carried out for revealing its hysteretic behaviors associated with implementation of semi-active control using the routine of computational fluid dynamics. By virtue of the structural symmetry of the MR damper, a two-dimensional configuration for finite element simulation is built up. Herschel-Bulkley model is employed to represent the property of the MR fluid, of which the control parameters and their relevances to the input current are addressed. Typical cases involving sinusoidal and irregular displacements, steady and transient currents loaded upon the MR damper are investigated. Numerical investigations reveal that the damper force has a positive correlation with input current, excitation amplitude and excitation frequency. The full-scale simulation is proved to exhibit a sound accuracy through the validation of experimental data. It provides a logical manner revealing the true performance of MR dampers under desirable operating modes in practice, and can be readily integrated with the gain design of the associated semi-actively controlled structure. This progress bypasses the technical challenge inherent in the traditional tests with low-frequency cyclic loadings due to the limitation of experimental setup. Besides, comparative study between two-dimensional and three-dimensional configuration simulations of the MR damper shows that former has a better applicability, which can be carried out on a low-cost platform.

A numerical study was undertaken to investigate the mechanical properties of metal-ceramic multilayer micro-pillars under compression. The model system contains alternating layers of aluminum (Al) and silicon carbide (SiC) above a thick silicon (Si) substrate. The finite element analyses show that the pillar deformed in a non-uniform manner under compression, especially when a tapered side wall was included in the model. The heterogeneous feature of the structure plays a significant role in influencing the apparent stress-strain response. The elastic property of the pillar compares favorably with the true composite modulus, but a large discrepancy was found in the plastic response. The difference in deformation behavior between homogeneous pillars and multilayer pillars, along with other geometric effects, were examined in detail.

The study aims to estimate the stress-strain constitutive behavior of intermetallic compounds (IMCs) observed in a solder interconnect from experimental nanoindentation responses through a modified analysis procedure for improved solution robustness based on Cheng and Cheng's and Dao et al.'s models. On the basis of parametric finite element nanoindentation simulation and dimensional analysis together with the concept of representative strain, a set of universal dimensionless functions are established, by which a forward and reverse analysis algorithm are created to predict nanoindentation responses from given elastoplastic properties and vice versa, respectively. The proposed analysis procedure is validated through comparison with the experimental nanoindentation responses and limited literature data. The results show that the proposed analysis procedure is an effective means for plastic property characterization of micro/nanoscale IMCs. The representative strain is found to be 0.056, which differs from the Dao et al.'s and Giannakopoulos and Suresh's estimate. Besides, though generally brittle and hard in nature, the IMCs in a micro/nanoscale thickness show high plasticity, and comprise a yield strength surpassing most typical engineering metals.

High-pressure die-casting (HPDC) AM50 tensile specimens were used to investigate characteristics of gas pores and its effect on mechanical properties of HPDC AM50 magnesium alloy. Combining microstructure morphology gained from optical microscopy, scanning electron microscopy (SEM), and three-dimensional (3D) reconstruction with the experimental data from uniaxial tensile testing, we pursued the relationship between gas pores and the mechanical properties of HPDC AM50 Mg alloy. Results indicate that comparing with 3D reconstruction models, 2D images like optical metallography images and SEM images have one-sidedness. Furthermore, the size and maximum areal fraction of gas pores have negative effects on the mechanical properties of HPDC AM50 Mg alloy. With increase of the maximum size of gas pores in the specimen, the ultimate tensile strength (UTS) and elongation decrease. In addition, with the maximum areal fraction becoming larger, both the UTS and elongation decrease linearly.

Mesh generation is a bottleneck for finite element simulations of biomolecules. A robust and efficient approach, based on the immersed boundary method proposed in [8], has been developed and implemented to generate large-scale mesh body-fitted to molecular shape for general parallel finite element simulations. The molecular Gaussian surface is adopted to represent the molecular surface, and is finally approximated by piecewise planes via the tool phgSurfaceCut in PHG [43], which is improved and can reliably handle complicated molecular surfaces, through mesh refinement steps. A coarse background mesh is imported first and then is distributed into each process using a mesh partitioning algorithm such as space filling curve [5] or METIS [22]. A bisection method is used for the mesh refinements according to the molecular PDB or PQR file which describes the biomolecular region. After mesh refinements, the mesh is optionally repartitioned and redistributed for load balancing. For finite element simulations, the modification of region mark and boundary types is done in parallel. Our parallel mesh generation method has been successfully applied to a sphere cavity model, a DNA fragment, a gramicidin A channel and a huge Dengue virus system. The results of numerical experiments show good parallel efficiency. Computations of electrostatic potential and solvation energy also validate the method. Moreover, the meshing process and adaptive finite element computation can be integrated as one PHG project to avoid the mesh importing and exporting costs, and improve the convenience of application as well.

Cet article présente un modèle de connaissance pour mieux comprendre et expliquer les mécanismes à l'origine de la formation de la soudure dans le procédé de soudage par point en nous limitant ici au cas d'assemblages symétriques de tôles nues en acier extra-doux. Dans une première approche, on considère un modèle simplifié prenant en compte les couplages entre les phénomènes électrique, thermique et métallurgique dans l'assemblage mais en supposant les surfaces de contact constantes aux interfaces électrode-tôle et tôle-tôle. Dans une seconde approche le modèle est complété par un calcul mécanique qui permet de calculer les évolutions des surfaces de contact. Les résistances de contact électrique et thermique sont mesurées en fonction de la température et de la pression sur un dispositif expérimental conçu à cet effet. La variation des profils des électrodes en fonction du nombre de soudures réalisées et ses conséquences sur les surfaces de contact sont de même présentées. Le modèle numérique est validé en comparant, d'une part, la taille et la forme des points soudés relevées expérimentalement et calculées numériquement et, d'autre part, les cinétiques thermiques mesurées et calculées dans les électrodes et les tôles. Les techniques d'instrumentation mises au point visent à limiter les perturbations inévitables générées par l'intrusion des capteurs. Les hypothèses et les jeux de données considérés dans le modèle pour définir les conditions interfaciales sont enfin discutés.

An energy-based damage model is proposed and applied to predict the fracture initiation of the sheet metal forming process. The fracture mechanism is investigated through the plastic energy dissipation. The concepts of the damaging work and the fracture energy are proposed for the quantitative description of damage evolution and crack initiation. The developed formulations are implemented into the finite element program ABAQUS to simulate the biaxial stretching of sheet metals and to predict the fracture strains. The material parameter needed in the damage model for fracture prediction is determined by the stress-strain history of the uniaxial tensile test. The forming limits for aluminum alloy sheets under various strain paths are predicted by the present approach and then compared to the measured data quoted from the literatures [1,2]. Good agreements are found between this study and the previous results.