COMPLAS 2021 is the 16th conference of the COMPLAS Series.
The COMPLAS conferences started in 1987 and since then have become established events in the field of computational plasticity and related topics. The first fifteen conferences in the COMPLAS series were all held in the city of Barcelona (Spain) and were very successful from the scientific, engineering and social points of view. We intend to make the 16th edition of the conferenceanother successful edition of the COMPLAS meetings.
The objectives of COMPLAS 2021 are to address both the theoretical bases for the solution of nonlinear solid mechanics problems, involving plasticity and other material nonlinearities, and the numerical algorithms necessary for efficient and robust computer implementation. COMPLAS 2021 aims to act as a forum for practitioners in the nonlinear structural mechanics field to discuss recent advances and identify future research directions.
Scope
COMPLAS 2021 is the 16th conference of the COMPLAS Series.
Reinforced concrete (RC) shear walls are often used as the main lateral-resisting component in the seismic design of buildings. They provide a large percentage of the lateral stiffness of the structure, and therefore, they may experience large shear stresses at some point under earthquake loading. Consequentially, to accurately predict their behavior, it is recommended to use detailed finite element (FE) modeling with appropriate non-linear constitutive models for concrete and steel. However, such types of simulations are challenging and could significantly increase the computational time required to obtain the analysis results. In this paper, we study the viability of creating an artificial neural-network-based surrogate model of the RC shear wall that is able to capture its nonlinear behavior and predict the results obtained with a detailed FE model, offering a much lower computational effort. For this purpose, we develop a detailed parametric non-linear FE model based on well-established practices and validated studies using the OpenSees finite element software framework. The FE model consists of multi-layer shell elements with the vertical and transverse reinforcement included as smeared rebar layers. The concrete layers implement a damage mechanism with smeared crack constitutive model, whereas the rebar layers consider a uniaxial plasticity material law. The parametric FE model is used to build a large database of RC walls of different sizes and characteristics, with their corresponding lateral load capacity that is obtained through the detailed non-linear pushover analysis. Finally, the obtained database is used to train and validate the ANN-surrogate model. The developed model is able to accurately predict the lateral load capacity of RC shear walls without the need of detailed FE modeling, thus drastically reducing the complexity and the computational time required for the numerical solution and providing a reliable and robust analysis alternative, with only small compromise of accuracy.
Abstract Reinforced concrete (RC) shear walls are often used as the main lateral-resisting component in the seismic design of buildings. They provide a large percentage of the lateral [...]
In accidents involving cars with pedestrians, the impact of the head on structural parts of the vehicle takes a significant risk of injury. If the head hits the windscreen, the injury is highly influenced by glass fracture. In pedestrian protection tests, a head impactor is shot on the windscreen while the resultant acceleration at the COG of the head is measured. To assess the risk of fatal or serious injury, a head injury criterion (HIC) as an explicit function of the measured acceleration can be determined. The braking strength of glass which has a major impact on the head acceleration, however, is not deterministic but depends on production-related micro-cracks on the glass surface as well as on the loading rate. The aim of the present paper is to show a pragmatic method, how to include the stochastic failure of glass in crash and impact simulations. The methodology includes a fracture mechanical model for the strain rate-dependent failure of glass, the experimental determination of the glass strength for the different areas of a windscreen (surface, edge, and screen-printing area), the statistical evaluation of the experimental data and the computation of a HIC probability distribution by stochastic simulation.
Abstract In accidents involving cars with pedestrians, the impact of the head on structural parts of the vehicle takes a significant risk of injury. If the head hits the windscreen, [...]
The Discrete Material Optimization (DMO) and the Shape Function with Penalization (SFP) constitute the state-of-the-art material interpolation techniques for identifying from a list of pre-defined candidate materials the most suitable one(s) for the structural domain. The candidate materials are represented on this list through their mechanical properties, and are interpolated within the domain of interest (DOI), whether that is the finite element (FE) domain or groups of FEs, so-called patches. Depending on the technique preferred to interpolate the mechanical properties within the DOI, a different type of weights is selected. Goal of the discrete material optimization problem (MOP) is to solve for these weights and determine for each FE/patch a unique material from the list. The current work extends the concept of the SFP technique by employing as weights the shape functions of the hyper-tetrahedral FE, the dimension of which is dynamically adapted depending on the number of candidate materials considered for the structural domain. This generalized hyper-tetrahedral FE constitutes what is defined as a simplex, and similar to the SFP technique each of its nodes is tied to a specific candidate material. In the context of discrete optimization and utilizing the shape functions of an abstract high-dimensional FE as weights for the candidate materials, the proposed interpolation technique secures the continuity between the number of candidate materials that can be considered for the structure, a feature lacking in the SFP technique. Additionally, given that the number of nodes forming the simplex FE is always one unit greater than the dimension of the space it is defined within, the dimension of the resulting MOP drops by one per DOI. The developed material interpolation technique is combined with the topology optimization problem (TOP) to formulate the concurrent material and topology optimization problem for compliance minimization of the structure. Finally, the latter is examined on the academic case study of the 3D Messerchmitt-B¨olkow-Blohm (MBB) beam for the case of the concurrent topology and discrete fiber orientation optimization problem.
Abstract The Discrete Material Optimization (DMO) and the Shape Function with Penalization (SFP) constitute the state-of-the-art material interpolation techniques for identifying from [...]
In the present work, analytical and numerical models have been developed to calculate shock wave stresses caused by laser plasma in a material interface. The input pressure causing the shock wave has been derived from a laser-matter interaction model. The material velocity, stresses, and strains versus time in the two materials and the interface have been computed by solving the jump equations for conservation of mass and balance of momentum. The material interface has been considered as an immovable boundary while the back free surface as an unrestrained boundary. Spall fracture strength of the interface was evaluated and was used as stripping criterion. The model has been used for the fast computation of interfacial stresses causing paint stripping on aluminium substrates and the subsequent fast assessment of stripping initiation. An explicit finite element model combined with the cohesive zone modelling method and a spall fracture model have been developed. These models have been compared to each other in terms of time, accuracy and input properties demand and to the experimental results.
Abstract In the present work, analytical and numerical models have been developed to calculate shock wave stresses caused by laser plasma in a material interface. The input pressure [...]
Spoolable thermoplastic composite pipe (TCP) is an ideal alternative to traditional, heavier metallic counterparts for deepwater riser applications. During operation the pipe is subjected to mechanical loads simultaneously with through-wall thermal gradients arising from the mismatch between temperatures of hot pipe contents and cool surrounding ocean. In this work, structural analysis of TCP under coupled thermomechanical loads is performed using the finite element method (FEM). Temperature-dependent material properties are considered. Material safety factors for different laminate stacking sequences are compared and multi-angle stacking is shown to be effective for both pressure- and tensiondominated scenarios. Safety factors are also generated for TCP bent at reduced and elevated temperatures illustrative of spooling in different environments. It is clear that optimising the laminate for operation will adversely affect spooling capacity and vice-versa, i.e. TCP intended for extreme in-service conditions will require large spools.
Abstract Spoolable thermoplastic composite pipe (TCP) is an ideal alternative to traditional, heavier metallic counterparts for deepwater riser applications. During operation the pipe [...]
Two chemo-mechanical coupled models for electrode particles of lithium-ion batteries are compared. On the one hand a CahnHilliard-type phase-field approach models lithium intercalation, phase separation and large deformations in phase transforming cathode materials like lithium iron phosphate. On the other hand a chemo-mechanical particle model for lithium intercalation and large deformations for an anode material such as silicon is studied. The comparison of two different ways to define the deformation gradient for the large deformation approach and the two different material properties lead to differences in the resulting quantities and equations for the coupling of the chemo-mechanical model. The usage of an adaptive solution algorithm as well as the parallelization of the finite element solver via the message passing interface concept results in a more reasonable computation time to perform two-dimensional simulations. Both materials are numerically investigated and the results are compared from a physical point of view. When fast charging a battery, higher stress values are reached, which can cause a shorter cycle life. A strong scalability analysis shows good performance for the assembling, however a saturation occurs in the performance of the solver used.
Abstract Two chemo-mechanical coupled models for electrode particles of lithium-ion batteries are compared. On the one hand a CahnHilliard-type phase-field approach models lithium [...]
The present paper aims at deriving closed-form expressions for the optimum winding angles of fibres in laminated cylindrical pressure vessels subjected to internal pressure, axial load and torque. To achieve this goal, the state of the stress on the surface of the vessel is represented by a composite layered element subjected to general in-plane loading. The symbolic software Mathematica is used to formulate, solve the optimization problem and to generate the data of optimum fibre angles for different loading combinations. The generated data is fitted with simple polynomial expressions capable of accurately predicting the optimum winding angles. The optimization is performed for three candidate composite materials: carbon/epoxy, glass/epoxy, and aramid/epoxy laminates.
Abstract The present paper aims at deriving closed-form expressions for the optimum winding angles of fibres in laminated cylindrical pressure vessels subjected to internal pressure, [...]
V. Giglioni, I. Venanzi, V. Poggioni, A. Baia, A. Milani, F. Ubertini
eccomas2022.
Abstract
During their life-cycle, civil infrastructures are continuously prone to significant functionality losses, primarily due to material's degradation and exposure to several natural hazards. Following these concerns, many researchers have attempted to develop reliable monitoring strategies, as integration to visual inspections, to efficiently ensure bridge maintenance and early-stage damage detection. In this framework, recent improvements in sensor technologies and data science have stimulated the use of Machine Learning (ML) algorithms for Structural Health Monitoring (SHM). Among unsupervised learning techniques, the potential of autoencoder networks has been attracting notable interest in the context of anomaly detection. In this light, the present paper proposes two different autoencoder-based damage detection techniques, focused on the Multi-Layer Perceptron (MLP) and the Convolutional Autoencoder (CAE) networks, respectively. During the training, the selected ML models learn how reconstructing raw acceleration sequences acquired from sound conditions. Unknown data, including both healthy and damaged bridge responses, are afterwards used to test the implemented networks and to detect damage occurrence. To this aim, a specific index of reconstruction loss is selected as a damage sensitive feature with the aim to quantify the errors between the original and reconstructed sequences. The performance exhibited by the two approaches is compared and evaluated by application to the Z24 benchmark bridge. Results demonstrate the effectiveness of the proposed methodology to perform feature classification and real time damage detection at the level of macro-sequences as new sensor data is collected, resulting suitable for continuous assessment of full-scale monitored bridges.
Abstract During their life-cycle, civil infrastructures are continuously prone to significant functionality losses, primarily due to material's degradation and exposure to several [...]
We advance this field by systematically exploring nonlinear interactions of A0 Lamb mode signals with delamination defects of various sizes and interlaminar locations, to facilitate the characterization of delaminations. The interrogation signal, in this regard, is a modulated sinusoid whose frequency varies in steps of 20 kHz between 40 kHz and 100 kHz. Commercial FEM software is used for modelling the contact at delamination interfaces and for simulating the Lamb wave propagation through the waveguide with delamination defect. It is demonstrated that the intermittently acting contact pressure between the two surfaces of delamination acts as a source of nonlinearity, resulting in generation of higher order harmonics of interrogation frequency. A metric for measuring nonlinearity, the nonlinearity index (NI), is used to quantify the strength of wave-damage interactions over a range of interrogation frequencies. The NI index is observed to vary with both the interlaminar location as well as the width of the delamination. The maximum value of NI is further influenced by the frequency of excitation signal. To infer the effect of delamination parameters on the NI, a concept of a concept of contact energy intensity is introduced, which is largely dependent on the size and the interlaminar position of the delamination. The nonlinearity index patterns are explained by combining the intensity of the contact energy with the phase difference between waves traveling through the sub-laminations and the flexural rigidity of the two sub-laminates at the delamination location. The inferences provided can potentially be used for determining the interlaminar location and width of delamination employing higher harmonic Lamb wave signals generated by the breathing delamination.
Abstract We advance this field by systematically exploring nonlinear interactions of A0 Lamb mode signals with delamination defects of various sizes and interlaminar locations, to [...]
In this paper a multi-stage numerical analysis is presented with the aim to investigate effects of forming on the high-cycle fatigue performance of a deep-drawn aluminum sheet structure for use in a floating photovoltaic system. Forming simulations of a subsection of the full geometry are performed in a realistic two-step drawing-springback cycle. A simplified global analysis of service load response is performed to obtain displacements at the submodel boundary, that are used to generate boundary conditions for a local service load analysis. The local analysis is then performed on three different models of the subsection: (I) excluding forming effects, (II) including effects of thinning, and (III) including effects of thinning and residual stresses. The critical location with respect to the fatigue limit criterion in the Ottosen-Stenström-Ristinmaa high-cycle fatigue model was identified for case (I), and this location was used to compare the different models to assess effects of forming on high-cycle fatigue performance. Furthermore, the dynamic friction coefficient , and the isotropickinematic mixing coefficient were varied in order to investigate their respective effects.
Abstract In this paper a multi-stage numerical analysis is presented with the aim to investigate effects of forming on the high-cycle fatigue performance of a deep-drawn aluminum sheet [...]