Evolutionary Inverse Material Identification: Bespoke Characterization of Soft Materials Using a Metaheuristic Algorithm

The growing interest in soft robotics has resulted in an increased demand for accurate and reliable material modelling. As soft robots experience high deformations, highly nonlinear behavior is possible. Several analytical models that are able to capture this nonlinear behavior have been proposed, however, accurately calibrating them for specific materials and applications can be challenging. Multiple experimental testbeds may be required for material characterization which can be expensive and cumbersome. In this work, we propose an alternative framework for parameter fitting established hyperelastic material models, with the aim of improving their utility in the modelling of soft continuum robots. We define a minimization problem to reduce fitting errors between a soft continuum robot deformed experimentally and its equivalent finite element simulation. The soft material is characterized using four commonly employed hyperelastic material models (Neo Hookean; Mooney–Rivlin; Yeoh; and Ogden). To meet the complexity of the defined problem, we use an evolutionary algorithm to navigate the search space and determine optimal parameters for a selected material model and a specific actuation method, naming this approach as Evolutionary Inverse Material Identification (EIMI). We test the proposed approach with a magnetically actuated soft robot by characterizing two polymers often employed in the field: Dragon Skin™ 10 MEDIUM and Ecoflex™ 00-50. To determine the goodness of the FEM simulation for a specific set of model parameters, we define a function that measures the distance between the mesh of the FEM simulation and the experimental data. Our characterization framework showed an improvement greater than 6% compared to conventional model fitting approaches at different strain ranges based on the benchmark defined. Furthermore, the low variability across the different models obtained using our approach demonstrates reduced dependence on model and strain-range selection, making it well suited to application-specific soft robot modelling.

Quasi-Analytical Calculation of Frequency-Dependent Resistance of Rectangular Conductors Considering the Edge Effect

This paper presents an accurate quasi-analytical approximation of frequency-dependent ac resistance of single rectangular conductors. In this work, first, a two-dimensional analytical ac resistance of rectangular conductors is derived. Unlike circular conductors, where current density distributes evenly in each layer of the conductor’s cross-section, the edge effect is involved for rectangular conductors. Due to the edge effect, one cannot define an accurate boundary condition for solving the two-dimensional partial differential equation of magnetic field or current density of rectangular conductors. Hence, the calculated two-dimensional analytical current density result is not accurate and is modified and fitted on FEM simulation, taking the conductor’s thickness into account using the least-square problem to improve its accuracy. Unlike numerical approaches, the proposed method yields an easy-to-use formula applicable to industrial applications in different fields. Contrary to the one-dimensional approach, which is only valid for very thin rectangular conductors, this method takes edge effect into account and can be used for any thickness (from square to very thin rectangular conductors). The proposed method can be used in applications where an accurate ac resistance of rectangular conductors over a wide frequency range is required, such as white-box modeling of power transformers and interpreting its frequency response analysis (FRA), and calculating the resistance of electric machine winding, busbars, and printed circuit board traces.

Development of fully superconducting 5 MW aviation generator with liquid hydrogen cooling

The use of liquid hydrogen as a fuel will be inevitable in the aviation of the future. This statement means that manufacturers will also implement liquid hydrogen for cooling all superconducting aviation equipment of an electric propulsion system. The development of fully electric aircraft is the most promising solution in this case. Scientists from the Department of electrical machines and power electronics of Moscow aviation institute have conducted calculations and theoretical researches of critical specific mass-dimensional parameters (MW/ton and MW/m3 at 21 K) of fully superconducting aviation synchronous generator of the electric propulsion system. The results are in this article. The article discusses the results 3D finite element modeling (FEM) simulation of a 5 MW fully superconducting synchronous generator with combined excitation. Superconducting armature and axial excitation windings based on second generation high temperature superconductors (HTS-2G) are located on the stator, which makes it possible to contactlessness and the absence of sliding seals. A dry gap will reduce gas-dynamic losses and increase the nominal peripheral speed of the rotor. The use of liquid hydrogen as a coolant makes it possible to significantly increase the linear load of the generator, and high current densities to reduce the cross-sectional area of the coils, which will make it possible to place them in individual cryostats in the future. Individual cryostats will allow to remove the heat release of magnetic losses from the cryogenic zone and reduce the consumption of refrigerant. For the purpose of internal redundancy of the HTS coils, the machine has a complete set of reserve winding made of ultrapure aluminum, also cooled by liquid hydrogen. If the superconducting coils get out of the stand, the generator will provide 15 % power on standby

Coupled CFD-FEM Simulation of Steel Box Bridge Exposed to Fire

Increasing fire-induced bridge failures are demanding more precise behavior prediction for the bridges subjected to fires. However, current numerical methods are limited to temperature curves prescribed for building structures, which can misestimate the fire impact significantly. This paper developed a framework coupling the computational dynamics (CFD) method and finite element method (FEM) to predict the performance of fire-exposed bridges. The fire combustion was simulated in CFD software, Fire Dynamic Simulator, to calculate the thermal boundary required by the thermomechanical simulation. Then, the adiabatic surface temperatures and heat transfer coefficient were applied to the FEM model of the entire bridge girder. A sequential coupled thermomechanical FEM simulation was then carried out to evaluate the performance of the fire-exposed bridge, thermally and structurally. The methodology was then validated through a real fire experiment on a steel beam. The fire performance of a simply supported steel box bridge was simulated using the proposed coupled CFD-FEM methodology. Numerical results show that the presented method was able to replicate the inhomogeneous thermomechanical response of box bridges exposed to real fires. The girder failed due to the buckling of a central diaphragm after the ignition of the investigated tanker fire in no more than 10 min. The framework presented in this study is programmatic and friendly to researchers and can be applied for the estimation of bridges in different fire conditions.

FEM Simulation of the Riveting Process and Structural Analysis of Low-Carbon Steel Tubular Rivets Fracture

Riveted joints are a common way to connect elements and subassemblies in the automotive industry. In the assembly process, tubular rivets are loaded axially with ca. 3 kN forces, and these loads can cause cracks and delamination in the rivet material. Such effects at the quality control stage disqualify the product in further assembly process. The article presents an analysis of the fracture mechanism of E215 low-carbon steel tubular rivets used to join modules of driver and passenger safety systems (airbags) in vehicles. Finite element method (FEM) simulation and material testing were used to verify the stresses and analysis of the rivet fracture. Numerical tests determined the state of stress during rivet forming using the FEM-EA method based on the explicit integration of central differences. Light microscopy (LM), scanning electron microscopy (SEM) and chemical composition analysis (SEM-EDS) were performed to investigate the microstructure of the rivet material and to analyze the cracks. Results showed that the cause of rivet cracking is the accumulation and exceeding of critical tensile stresses in the rivet flange during the tube processing and the final riveting (forming) process. Moreover, it was discovered that rivet fracture is largely caused by structural defects (tertiary cementite Fe,Mn3CIII along the boundaries of prior austenite grains) in the material resulting from the incorrectly selected parameters of the final heat treatment of the prefabricate (tube) from which the rivet was produced. The FEM simulation of the riveting and structural characterization results correlated well, so the rivet forming process and fracture mechanism could be fully investigated.

Effect of Process Parameters on the Microstructure and Mechanical Properties of Bars from Al-Cu-Mg Alloy Processed by Multipass Radial-Shear Rolling

The study of microstructure and mechanical properties formation of A2024 alloy obtained by the multipass radial-shear rolling (RSR) method is discussed in this article. FEM simulation was carried out that made it possible to evaluate the influence degree of rolling temperature-velocity parameters on the strain state of material. It has been found the increase in rotary velocity of rolls significantly influences on the deformation heating of bar after RSR (predominantly in its surface layer). The combination of rolling temperature-velocity conditions at selection of deformation regime has complex effect on structure and properties formation. The analysis of sizes and distribution of phase particles has shown that the rolling at lower temperatures allowed to increase the mechanical strength due to the more intensive refinement of undissolved Fe-containing phase. The gradual decrease in the rolling temperature in each pass makes possible to achieve the high strength (UTS~430 MPa and YS~255 MPa) while maintaining the ductility level ~15%, that are comparable to ones obtained at some severe plastic deformation (SPD) methods.

Influence of 3D Printing Topology by DMLS Method on Crack Propagation

The presented text deals with research into the influence of the printing layers’ orientation on crack propagation in an AlSi10Mg material specimen, produced by additive technology, using the Direct Metal Laser Sintering (DMLS) method. It is a method based on sintering and melting layers of powder material using a laser beam. The material specimen is presented as a Compact Tension test specimen and is printed in four different defined orientations (topology) of the printing layers—0°, 45°, 90°, and twice 90°. The normalized specimen is loaded cyclically, where the crack length is measured and recorded, and at the same time, the crack growth rate is determined. The evaluation of the experiment shows an apparent influence of the topology, which is essential especially for possible use in the design and technical preparation of the production of real machine parts in industrial practice. Simultaneously with the measurement results, other influencing factors are listed, especially product postprocessing and the measurement method used. The hypothesis of crack propagation using Computer Aided Engineering/Finite Element Method (CAE/FEM) simulation is also stated here based on the achieved results.

Numerical Investigation of Plastic Strain Homogeneity during Equal-Channel Angular Pressing of a Cu-Zr Alloy

A three-dimensional finite element method (3D FEM) simulation was carried out using ABAQUS/Explicit software to simulate multi-pass processing by equal-channel angular pressing (ECAP) of a circular cross-sectional workpiece of a Cu-Zr alloy. The effective plastic strain distribution, the strain homogeneity and the occurrence of a steady-state zone in the workpiece were investigated during ECAP processing for up to eight passes. The simulation results show that a strain inhomogeneity was developed in ECAP after one pass due to the formation of a corner gap in the outer corner of the die. The calculations show that the average effective plastic strain and the degree of homogeneity both increase with the number of ECAP passes. Based on the coefficient of variance, a steady-state zone was identified in the middle section of the ECAP workpiece, and this was numerically evaluated as extending over a length of approximately 40 mm along the longitudinal axis for the Cu-Zr alloy.

Analysis of Interference-Fit Joints

Interference fit joints have been widely used in many engineering constructions, in particular in electric motors. It is of particular importance to calculate the load capacity of press-fit joints, especially in the overload ranges of construction to estimate the safety factor. The article presents a FEM numerical simulation of pressing the shaft into the hub, taking into account various types of fits. The results of numerical simulations presented in the article were positively verified with the MTS measuring device, which confirmed the correctness of the numerical model. So far, the load-bearing capacity of press-fit joints has been calculated from Lame’s formulas. The results of the load capacity of the joints obtained by the FEM simulation were compared with the results obtained from Lame’s formula. The comparison shows that when designing interference fit joints, attention should be paid to the fact that the press-in process, depending on the type of fit, may be elastic-plastic. Plastic deformations in the contact zone of the joint affect its load-bearing capacity. Therefore, the design of press-fit joints should not be based on Lame’s formulas, which do not take into account the range of plastic work of the material.