This paper presents an algorithm for the fully automatic mesh generation for the finite element analysis of ships and offshore structures. The quality requirements on the mesh generator are imposed by the acceptance criteria of the classification societies as well as the need to avoid shear locking when using low degree shell elements. The meshing algorithm will be generating quadrilateral dominated meshes (consisting of quads and triangles) and the mesh quality requirements mandate that quadrilaterals with internal angles close to $90\degree$ are to be preferred. The geometry is described by a dictionary containing points, rods, surfaces and openings. The first part of the proposed method consist of an algorithm to automatically clean the geometry. The corrected geometry is then meshed by the frontal Delaunay mesh generator as implemented in the gmsh package. We present a heuristic method to precondition the cross field of the frontal quadrilateral mesher. Also the influence of the order in which the plates are meshed will be explored as a preconditioning step.

Rail impedance directly affects the transmission performance of track circuit . Considering the condition of earth stratification, for the difficult to calculate the rail impedance due to the semi-infinite integration interval and the oscillation of the integrand by using the Carson formula, The truncation method is proposed to divide the impedance formula is divided into definite integral and tail integral. The integral is approximated by the spline function, and the tail integral is calculated by using the exponential integral and Euler formula. Based on it, the rail impedance calculation formula of track circuit is obtained. The electromagnetic field model of track circuit with earth stratification is simulated by finite element method, and the correctness of the method is verified. Based on the formula, the influence of current frequency, soil depth and conductivity on rail impedance is studied. The relative error between the calculated results of rail impedance and the simulation results of finite element is within 5%. It can be seen that the formula has high accuracy and correctly reflects the law of rail impedance variation with current frequency, soil depth and resistivity. It provides a reliable reference for the theoretical calculation of rail impedance of track circuit.

The analytical method for the calculation of the properties of a bolted connection established by the structural Eurocodes proposes the T-stub as a component for the characterization of the tension and compression zones in moment joints. In this article, a review of the state of the art on the T-stub component is developed, where the works developed since it was initially defined, and from the perspectives of formulation, experimentation and numerical simulation are summarized and discussed. Additionally, possible future lines of work are proposed.

The non-linear analysis of hollow-core concrete slabs requires the use of advanced numerical techniques, proper constitutive models both for concrete and steel as well as particular computational skills. If prestressing, cracking, crack opening, material softening, etc. are also to be taken into account, then the computational task can far exceed the capabilities of an ordinary engineer. In order for the calculations to be carried out in a traditional design office, simplified calculation methods are needed. Preferably based on the linear finite element (FE) method with a simple approach that takes into account material nonlinearities. In this paper the simplified analysis of hollow-core slabs based on the generalized nonlinear constitutive law is presented. In the proposed method a simple decomposition of the traditional iterative linear finite element analysis and the non-linear algebraic analysis of the plate cross-section is used. Through independent analysis of the plate cross-section in different deformation states, a degraded plate stiffness can be obtained, which allows iterative update of displacements and rotations in the nodes of the FE model. Which in turn allows to update the deformation state and then correct translations and rotations in the nodes again. The results obtained from the full detailed 3D nonlinear FEM model and from the proposed approach are compared for different slab cross-sections. The obtained results from both models are consistent.

This study investigates the failure of a roof with steel truss construction of a factory building in Tekirdag in North-western part of Turkey. The failure occurred under hefty weather conditions including thunderbolt, lightning strikes, heavy rain and fierce winds. In order to interpret the reason for the failure, the effects of different combinations of factors on the design and dimensioning of the roof are checked. Therefore, finite element analysis is performed several times under different assumptions and considering different factors aiming to determine the dominant ones that are responsible for the failure using the commercial software Abaqus (Dassault Systèmes, Vélizy-Villacoublay, France). Each loading condition gives out a characteristic form of failure. The scenario with the most similar form of failure to the real collapse is considered as the most likely scenario of failure. Also, the factors included in this scenario are expected to be the responsible factors for the partial collapse of the steel truss structure.

To investigate the effect of Polyvinylpyrrolidone (PVP) addition and consequently porosity, two different sets of membranes are manufactured, since PVP is a widely used poring agent which has an impact on the mechanical properties of the membrane material. The first set (PAN 1) includes PAN and the necessary solvent while the second set (PAN 2) is made of PAN and PVP. These membranes are put through several characterisation processes including tensile testing. The obtained data are used to model the static behaviour of the membranes with different geometries, but similar loading and boundary conditions that represent their operating conditions. This modelling process is undertaken by using finite element method. The main idea is to investigate how geometry affects the load-carrying capacity of the membranes. Alongside membrane modelling, their materials are modelled with representative elements with hexagonal and rectangular pore arrays (RE) to understand the impact of porosity on the mechanical properties. Exploring the results, the best geometry is found as the elliptic membrane with the aspect ratio 4 and the better RE as the hexagonal array which can predict the elastic properties with an approximate error of 12%.

Aerated concrete, which is manufactured from binding material, sand, foaming agent and water, is currently being utilized in the construction industry because of its lightweight and durability. The binding material, cement, along with other materials used in the concrete produces huge carbon footprints during its fabrication. The utilization of natural aggregates name as coarse aggregates depletes the natural resources of the country. Therefore, huge amounts of agricultural wastes have led scholars to investigate the effectiveness of replacing conventional materials used in concrete with agricultural wastes. In the current study, rice husk ash (RHA) was used as supplementary cementing material, thereby reducing the amount of cement used in aerated concrete (AC) mixture will reduce carbon footprints. The experimental and numerical analysis were conducted to investigate structural behavior of reinforced RAC- B beams subjected to flexural load. Parametric study on structural performance of RAC- B beam under flexure were conducted using finite element analysis (FEA). From the experiment and FEA. Results from the parametric study showed that RAC-10%RHA-B with higher depth structurally performed better compared to RAC-B under flexure with greater load carrying capacity, lesser maximum deflection, and less cracks developing in the tension area.

Estimating the vehicle crashworthiness parameters experimentally is expensive and time consuming. For these reasons different modelling approaches are utilized to predict the vehicle behaviour and reduce the need for full-scale crash testing. The earlier numerical methods used for vehicle crashworthiness analysis were based on the use of lumped parameters models (LPM), a combination of masses and nonlinear springs interconnected in various configurations. Nowadays, the explicit nonlinear finite element analysis (FEA) is probably the most widely recognized modelling technique. Although informative, finite element models (FEM) of vehicle crash are expensive both in terms of man-hours put into assembling the model and related computational costs. A simpler analytical tool for early analysis of vehicle crashworthiness could greatly assist the modelling and save time. In this paper a simple piecewise LPM composed of a mass-spring-damper system, is used to estimate the vehicle crashworthiness parameters, focusing on the dynamic crush and the acceleration severity index (ASI). The model is first calibrated against a full-scale crash test and a FEM, post-processed with the LS-DYNA software, at an impact velocity of 56 km/h. The genetic algorithm is used to calibrate the model by estimating the piecewise lumped parameters (stiffness and damping of the front structure of the vehicle). After calibration, the LPM is applied to a range of velocities (40, 48, 64 and 72 km/h). The predictions for crashworthiness parameters from the LPM were compared with the predictions from the FEA and the results are much similar. It is shown that the LPM can assist in crash analysis, since LPM has some predictive capabilities and requires less computation time in comparison with the explicit nonlinear FEA.

Transient analysis of gas flow in pipeline was studied. Finite Element Method based on ANSYS FLOTRAN was used to account for changes in pressure, temperature and flow rate. Compressibility factor function of temperature and pressure was considered. For non-isothermal transient results, the pressure and the flowrate gave wave propagation as a result of slow transients created by demand condition at the outlet end of custody transfer. Results obtained were in agreement with the demand restrictions at the outlet end of custody transfer indicating that the predictions are accurate and reliable. The results demonstrated that the Finite Element Method gave accurate prediction of pressure, temperature and flowrate in transient gas studies. For steady state non-isothermal model, results showed that the magnitude of the average pressure drop was higher when pressure was predicted with constant compressibility factor, but the same average pressure drop was reduced when the pressure was predicted with variable compressibility factor, z. Since compressibility factor is a function of temperature and pressure, the above findings signifies that in the case when gas temperature does not stabilize, the prediction of pressure with isothermal model and constant compressibility factor will lead to significant errors.

Purpose: The aim of the present experimental 3D-finite element study was to evaluate the influence of an augmented sinus lift with an additional inserted bone graft block. The bone graft block stabilizes the implant in addition to conventional augmented bone. We placed the block in three different positions. The implants were loaded with axial force and forces secondary to laterotrusion and protrusion. Material and Methods: A simplified U-shaped 3D finite element model of the upper jaw and a more complex anatomical model of the left maxilla were created. The bone graft block was placed in three positions: in the lower third in contact with the sinus floor, the middle, and the upper third of the implant. Van Mises’ stress distribution was calculated and analyzed for the different models. We also compared the complex anatomical model with the simplified one. Results: The position of the bone graft block significantly influences the magnitude of stress distribution. A bone graft block positioned in the upper third or middle of the implant reduces the quantity of stress compared to the reference model without a bone graft block. The low bone graft block position is clearly associated with lower stress distribution in compact bone. We registered no significant differences in stress in compact bone with regard to laterotrusion or protrusion. Conclusions: Maximum values of von Mises stresses in compact bone can be reduced significantly by using a bone graft block. The reduction of stress is nearly the same for positions in the upper third and the middle of the implant. It is much more pronounced when the bone graft block is in the lower third of the implant near the sinus floor, which appeared to be the best position in the present study.

In this study, we design a new refractive index sensor based on a metal–insulator–metal waveguide coupled with a notched ring and vertical rectangular resonators. We use the finite element method to study the propagation characteristics of the sensor. According to the calculation results, the transmission spectrum exhibits a typical Fano resonance shape. The phenomenon of Fano resonance is caused by the coupling between the wide spectrum and narrow spectrum. In the design, the wide spectrum signal is generated by the vertical rectangular resonator, while the narrow spectrum signal is generated by the notched ring resonator. In addition, we varied the structural parameters of the resonators and filled the structure with media of different refractive indices to study the sensing properties. The maximum achieved sensitivity of the sensor reached 1071.4 nm/RIU. The results reveal potential applications of the coupled system in the field of sensors.

This paper presents an approach to the Euler-Bernoulli beam theory (EBBT) using the finite difference method (FDM). The EBBT covers the case of small deflections, and shear deformations are not considered. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected points of the structure. Generally, the finite difference approximations are derived based on fourth-order polynomial hypothesis (FOPH) and second-order polynomial hypothesis (SOPH) for the deflection curve; the FOPH is made for the fourth and third derivative of the deflection curve while the SOPH is made for its second and first derivative. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. In this paper, the FOPH was made for all of the derivatives of the deflection curve, and additional points were introduced at the beam’s ends and positions of discontinuity (concentrated loads or moments, supports, hinges, springs, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends and to satisfy the boundary and continuity conditions. Moreover, grid points with variable spacing were also considered, the grid being uniform within beam segments. First-order analysis, second-order analysis, and vibration analysis of structures were conducted with this model. Furthermore, tapered beams were analyzed (element stiffness matrix, second-order analysis). Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. The results obtained in this paper showed good agreement with those of other studies, and the accuracy was increased through a grid refinement. Especially in the first-order analysis of uniform beams, the results were exact for uniformly distributed and concentrated loads regardless of the grid. Further research will be needed to investigate polynomial refinements (higher-order polynomials such as fifth-order, sixth-order…) of the deflection curve; the polynomial refinements aimed to increase the accuracy, whereby non-centered finite difference approximations at beam’s ends and positions of discontinuity would be used.

Earthquakes present a significant risk to marine and underwater structures such as seawalls, piers, dolphins, breakwaters, buried pipelines, and sheet-piled constructions. These structures are particularly susceptible to the destabilizing effects of seismic events, which can result in severe damage. Furthermore, earthquakes can indirectly threaten these installations by destabilizing the underlying soil. Soil liquefaction, induced by seismic shaking, can undermine the stability and integrity of structures. This review discusses advanced simulation programs that have been systematically developed and validated using comparisons with commercial software and existing research . It extensively examines various factors that influence earthquake-induced dynamics in marine environments, including soil conditions, quality of structural components, tunnel length, mooring intervals, seismic wave propagation, and seaquake impacts . Over the past four decades, significant insights have been gained into the seismic design of marine structures, leading to the development of comprehensive guidelines. Recently, the simulation of marine infrastructure—such as breakwaters, platform jackets, and tunnels—has progressed considerably, driven by rapid economic growth and technological advancements in engineering. This review aims to evaluate and categorize recent studies on the numerical simulation of marine structures under seismic forces, with a focus on mitigating earthquake-related damage.

This study determines the equivalent stress intensity factor (SIF) model that best fits the experimental behavior of low-carbon steel under mixed mode (I and II). The study assessed Tanaka, Richard, and Pook's equivalent SIF models. The theoretical values used for comparison correspond to experimental results in a modified C(T) geometry by machining a hole ahead of the crack tip subjected to fatigue loads with a load ratio of R = 0.1. The comparison involved the SIF for six experimental points and the values computed through the numerical simulation. The Paris, Klesnil, and modified Forman-Newman crack growth models were used with each equivalent SIF to analyze the prediction in the estimated number of cycles. The Klesnil model showed the closest prediction since the error between the calculated and experimentally recorded number of cycles is the lowest. However, the material behavior reflects a reduced crack propagation rate attributed to plasticity in the crack tip. Results suggest that Asaro equivalent SIF conservatively estimates the element lifespan with increasing errors from 2.3% at the start of growth to 27 % at the end of the calculation.

Eddy current based approaches have been investigated for a wide range of inspection applications. Dodd-Deeds model and the truncated region eigenfunction expansion (TREE) method are widely applied in various occasions, mostly for the cases that the sample is relatively larger than the radius of the sensor coil. The TREE method converts the integral expressions to the summation of many terms in the truncated region. In a recent work, the impedance of the co-axial air-cored sensor due to a plate of finite radius was calculated by the modified Dodd-Deeds analytical approach proposed by authors. In this paper, combining the modified analytical solution and the TREE method, a new finite region eigenfunction expansion (FREE) method is proposed. This method involves modifying its initial summation point from the first zero of the Bessel function to a value related to the radius of the plate, therefore makes it suitable for plate with finite dimensions. Experiments and simulations have been carried out and compared for the verification of the proposed method. Further, the planar size measurements of the metallic circular plate can be achieved by utilising the measured peak frequency feature.

Functional graded materials (FGMs) have been successfully applied in improving the insulation performance of insulators, yet there is currently no relevant literature on their application in dry-type transformers. Therefore, it is meaningful to study on the feasibility of using FGMs to improve the insulation performance of dry-type transformers through research on electric field distribution. This paper sets FGMs in the dry-type transformer spacer only or both spacer and winding end to change the relative permittivity values in the corresponding areas. By comparing the electric field distribution characteristics on the inner and outer surfaces of the high-voltage coil and the lower surface of the spacer of the dry-type transformer under each control method with the case without setting FGMs, the study found that when FGMs were only set on the spacer, the control methods set in this paper could not improve the electric field distribution of the transformer. When FGMs are employed at the spacers and end windings, two modulation strategies can mitigate the electric field intensity on the outer surface of high-voltage windings. Specifically, a continuous increase in relative permittivity can reduce the electric field intensity by 34.02%, while an initial decrease followed by an increase can result in a reduction of 13.58%. This paper provides theoretical guidance for insulation design of transformers with high voltage levels of 66 kV and above, especially as transformers with high voltage levels face greater challenges in regulating surface electric fields.

This study presents the experimental and analytical evaluations on an innovative connection class under cyclic loading. The connection was an H- or I-shaped section as a short stub column (SSC) between the end-plate and the column. The stub section was bolted to the column flange on the one side and bolted to the beam along with the welded end-plate on the other side. Stiffeners were installed on the webs and flanges of the beam and stub to ensure the necessary stiffness of the connection. The variables included thickness of stub flange/web and height of stub web. The cyclic loading was applied on a sample of the proposed connection according to AISC ANSI/AISC 341-16 loading protocol. Experimental observations displayed that the proposed connection, one can either employ hot rolling sections or create the connection as a plate girder, in which case the flange thickness and stub web, as well as the distance between two flanges and their height, was obtained according to the demand and performance of the structure. The optimum cross-section of the stub was then obtained. Also, the analytical model was calibrated using the finite element method (FEM), and its performance was confirmed. The effective parameters of the proposed connection were determined by modeling in the FEM software, and the design method was presented.

According to the honeycomb structure, the preliminary design of bionic wheel hub for finite element analysis, in the working condition of a wheel by the force is through the spokes pointing to the center, the working condition of two wheel by the force is through the two spokes pointing to the center, and the ordinary wheel stress and displacement value results compared, observed the maximum stress value appears in the connection between the center of the wheel and the spokes, and predict the fatigue life of the preliminary design of bionic wheel. Using finite element analysis and response surface optimization analysis method for numerical simulation, through two optimization results show that, in the use of response surface analysis model optimization, the optimized bionic wheel maximum stress is 109.34MPa, compared with the ordinary wheel maximum stress 119.77MPa, the maximum stress reduced by 8.7%. The mass of the ordinary wheel is 34.02Kg, using the response surface optimization method, the optimized wheel mass is 29.89Kg, the mass reduced by 12.13%. The optimized wheel fatigue analysis, in the stress of 109.34MPa, load loading times at least 4.217e+005, in the position at the junction with the half shaft. In line with the national standard GB/T 5334-2021[16] for commercial vehicle wheel fatigue life requirements.

Background: Good reduction of the articular surface and restoration of the congruent joint space are crucial for promoting positive outcomes in ankle fracture treatment. A good understanding of the mechanism of ankle fractures is indispensable to achieve these goals. Unfortunately, until now, the mechanism has still not been fully elucidated, particularly regarding the posterior malleolar fractures (PMFs). This study aims to provide insight into the mechanism of PMFs, thereby grasping the key to proceeding with PMFs reduction and fixation. Methods: An FE model of a normal ankle in tiptoe posture with 9 collateral ligaments was constructed from CT images. On the model, lateral malleolar fractures were made below, through, and above the syndesmosis; the medial malleolus was broken at its base; and ligaments were suppressed to simulate ligament rupture. A force of 1000 Newtons was applied on the talus head along talus axis. The posterior malleolar fragment (PMF) displacement, whether the talus impacts the posterior malleolus, and the initial mechanism displacing the PMF, were recorded. Results: A total of 15 tests were conducted. The tests, 'all collateral ligaments rupture', 'LMFB', 'LMFC+AITFL rupture', 'ATFL+CFL+PTFL rupture', 'ATTL+PTTL+TCL rupture', 'MMF', and 'MMF+AITFL rupture', presented significant PMF displacement. In tests 'MMF' and 'ATTL+PTTL+TCL rupture', the talus hit the medial edge of the PMF; but in tests 'ATFL+CFL+PTFL rupture' and 'all collateral ligaments rupture', the talus impacted the center part of the PMF. Except tests 'all collateral ligaments rupture' and 'ATFL+CFL+PTFL rupture', where PMF displacement was initially caused by talus impact, the other tests demonstrated PMF displaced following PITFL drag. Conclusions: This research discovered that rotation is the initial and major mechanism causing PMF displacement. The subsequent impact of the talus on the posterior malleolus may be the reason for large/blocky PMF. The talus impacting the posterior malleolus is always combined with talus rotation. Without the medial ankle structure being broken, big/blocky PMF may not occur. These findings suggest that, in clinical practice, reconstructing both the medial and lateral ankle structures is critical for achieving ankle stability, and this procedure even takes precedence over posterior malleolar fracture fixation itself.

This study focuses on calibrating the linear Drucker-Prager yield criterion for polypropylene under multiaxial loading conditions, using the Arcana machine to achieve a wide range of triaxiality states. The Drucker-Prager criterion has been adapted for polymers, composites, and metals due to its ability to account for the influence of hydrostatic pressure on yield stress. This study employs the Arcana machine, which allows for testing flat plate samples under various angles, thereby simulating different stress states from pure shear to biaxial tension. Additionally, uniaxial compression tests were conducted to extend the range of triaxiality. The material used in this study is Sabic 83MF10, a polypropylene polymer. Samples were prepared by injection molding and cut to specific dimensions. Tests were performed using a universal testing machine with the Arcana handle, and the results were analyzed to determine the yield strength, pressure, and triaxiality for each sample. The results showed a linear relationship between von Mises stress and hydrostatic pressure, with a friction angle (β) of 20.65° and material cohesion (d) of 27.81. The numerical simulations in Abaqus confirmed the validity of the Drucker-Prager model, accurately reflecting the moment of yielding for the tested samples.

To enhance the uplift capacity and facilitate the installation of multi-row perfobond connectors at shallow burial depths, this study puts forward a novel notched T-perfobond connector. The design incorporates an integrated flange at the bottom of the connector and a notch at the edge of the hole. Through pull-out model tests on four notched T-perfobond connectors, this research investigates their failure mechanisms and pull-out capacities. Utilizing the explicit dynamics method in Abaqus, a finite element model of the pull-out resistance test for the notched T-perfobond connector is established and verified against experimental data. Furthermore, a detailed parametric analysis involving 54 models is conducted, examining crucial parameters such as rib dimensions, hole geometry, flange size, notch width, bar diameter, and material properties. Based on the combined experimental and numerical results, the paper assesses the suitability of current formulas for calculating the pull-out capacity of perfobond connectors and proposes a refined calculation method specifically for notched T-perfobond connectors. All the findings in this paper can serve as a reference in the design and construction of composite structures.

During the thinning process, in the semiconductor industry, residual stress from thin films can cause severe warpage in wafers, which can compromise their processability and packaging. This warpage can lead to bifurcation or buckling, making it important to understand this phenomenon for the manufacturing of large semiconductor wafers. We conducted a finite element analysis of the thinning process using ANSYS mechanical enterprise 2023/R2 to investigate the asymmetric curvatures resulting from bifurcation in a standard 300 mm Si (001) wafer. We induced an asymmetry in the system by applying a slight force on a point of the circumferential region perpendicularly to the substrate. We also simulated the compressive stress from the damaged layer formed during the thinning process and considered its influence on warpage and bifurcation in both linear and non-linear cases. We also included the influence of Earth's gravity induced deflection (GID) and investigated the influence of the compressive stress damaged layer on the grinding process of the whole wafer by adding a metal front layer to the ground back damaged wafer.

Traumatic Brain Injury represents a significant public health concern due to its role in traumatic death and disability, often caused by head impacts or rapid accelerations. Despite advancements in Finite Element Head Models, existing models primarily represent adult males, neglecting anatomical and physiological differences between genders and age groups. Additionally, the mechanistic roles of specific intracranial tissues, such as the dura mater, have conventionally been overlooked. To better understand the mechanical behavior of the brain and its injured regions in the context of TBI, including the dura mater, this study proposes the development of a Female Finite Element Head Model. This model, representative of a middle-aged female subject, was developed using medical image-derived geometry and finite element modeling techniques. The validation results demonstrate a similarity between the numerical displacement curves and those obtained experimentally. In the model with dura mater, the numerical results show analogous behavior to the experimental results, despite minimal variations in the amplitude of the curves. These advancements highlight the significance of including the dura mater in biomechanical brain modeling. Incorporating the dura mater in Finite Element Head Models produces more realistic results than models without it. This emphasizes the importance of enhancing biomechanical models to better represent anatomical complexity.

Magnetic levitation vibration isolation platform are essential for testing space telescopes. This paper presents the design of a Large Air-Gap Voice Coil Motor (LAG-VCM) specifically for use in such platforms. The electromagnetic design of the voice coil motor begins with the Equivalent Magnetic Circuit Method, followed by the analysis and optimization of the motor design using Finite Element Analysis (FEA). Subsequently, a prototype was constructed and tested. The test results confirm that the prototype meets the design specifications and requirements.

Abstract Background Context Pelvic ring fractures are becoming more common in the aging population and can prove to be fatal, having mortality rates between 10% and 16%. Stabilization of these fractures is challenging and often require immediate internal fixation. Therefore, it is necessary to have a biomechanical understanding of the different fixation techniques for pelvic ring fractures. Methods A previously validated three-dimensional finite element model of the lumbar spine, pelvis, and femur was used for this study. A unilateral pelvic ring fracture was simulated by resecting the left side of the sacrum and pelvis. Five different fixation techniques were used to stabilize the fracture. A compressive follower load and pure moment was applied to compare different biomechanical parameters including range of motion (contralateral sacroiliac joint, L1-S1 segment, L5-S1 segment), and stresses (L5-S1 nucleus stresses, instrument stresses) between different fixation techniques. Results Transiliac -trans sacral screw fixation at S1 and S2 showed the highest stabilization for horizontal and vertical displacement at the sacral fracture site and reduction of contralateral sacroiliac joint for bending and flexion range of motion by 165% and 121%, respectively. DTSF (Double transiliac rod and screw fixation) model showed highest stabilization in horizontal displacement at the pubic rami fracture site, while the L5_PF_W_CC (L5-Ilium posterior screw fixation with cross connectors) and L5_PF_WO_CC ( L5-Ilium posterior screw fixation without cross connectors) showed higher rod stresses, reduced L1-S1 (approximately 28%), and L5-S1 (approximately 90%) range of motion. Conclusions Longer sacral screw fixations were superior in stabilizing sacral and contralateral sacroiliac joint range of motion. Lumbopelvic fixations displayed a higher degree of stabilization in the horizontal displacement compared to vertical displacement of pubic rami fracture, while also indicating the highest rod stresses. When determining the surgical approach for pelvic ring fractures, patient-specific factors should be accounted for to weigh the advantages and disadvantages for each technique.

The fifth Industrial revolution (I5.0) prioritizes resilience and sustainability, integrating cognitive cyber-physical systems and advanced technologies to enhance machining processes. Numerous research studies have been conducted to optimize machining operations by identifying and reducing the sources of uncertainty and estimating the optimal cutting parameters. Virtual modeling and Tool Condition Monitoring (TCM) methodologies have been developed to assess the cutting states during machining processes. With a precise estimation of cutting states, the safety margin necessary to deal with the uncertainties can be reduced, resulting in improved process productivity. This paper reviews the recent advances in high-performance machining systems with a focus on the cyber-physical models developed for the cutting operation of difficult-to-cut materials using cemented carbide tools. An overview of the literature and background on the advances in the offline and online process optimization approaches have been presented. Process optimization objectives such as tool life utilization, dynamic stability, enhanced productivity, improved machined part quality, and reduced energy consumption and carbon emissions are independently investigated for the offline and online optimization methods. Addressing the critical objectives and constraints prevalent in industrial applications, the paper explores the challenges and opportunities inherent in developing a robust cyber-physical optimization system

Transcutaneous nerve stimulation has several neuromodulatory applications, including pain management, sensory restoration, and motor recovery. Transcutaneous stimulation delivered via electrode arrays offers ease of use and reconfigurability through non-invasive means. Also, employing electrode arrays allows regulation of the permeated electric fields for greater control over the stimulation performance. Field shaping involves the regulation of permeated electric fields and is a well-established practice for deep brain stimulation; however, its applicability to peripheral nerve stimulation still needs exploration. Adopting such techniques can help to overcome current limitations with transcutaneous stimulation, like poor targetability, early onset of fatigue, and perceived discomfort. This article provides a perspective on the applicability of field-shaping techniques in improving the overall stimulation outcome.

Microvascular self-healing composite materials have significant potential for application and their mechanical properties need in-depth investigation. In this paper, tensile and compressive properties of woven fabric carbon fiber reinforced polymers (CFRP) laminates containing three-dimensional microvascular channels were investigated experimentally. Several detailed finite element (FE) models were established to simulate the mechanical behavior of the laminate and the effectiveness of different models were examined. The damage propagation process of the microvascular laminates and the influence of microvascular parameters were studied by the validated models. The results show that microvascular channels arranged along the thickness direction (z-direction) of the laminates are critical locations under the loads. The channels have minimal effect on the stiffness of the laminates but causing a certain reduction in strength, which varies approximately linearly with the z-direction channel diameter within its common design range of 0.1~1mm. It is necessary to consider the resin-rich region formed around microvascular channels in the warp and weft fiber yarns of the woven fabric composite when establishing the FE model. And the layers in model should be assigned with equivalent unidirectional ply material in order to calculate the mechanical properties of laminates correctly.

This study investigates the elastic buckling behaviour of skew Functionally Graded Material (FGM) thin plates featuring a circular opening. FGM materials, known for their unique property gradients, have gained prominence in structural engineering due to their mechanical performance and durability. Including a circular opening introduces a critical geometric consideration, influencing the structural stability and load-carrying capacity of FGM plates. The study examines the effects of the skew angle, plate’s aspect ratio, opening position, and size on the critical buckling load, normalized buckling load, and various buckling failure modes through computer modelling and finite element analysis. The results offer valuable insights into the interplay between material heterogeneity, geometric configuration, and structural stability. For instance, the critical buckling load increases by 29%, 82%, and 194% with the increment of skew angle from 0° to 30°, 45°, and 60°, respectively. Moreover, as the opening shifts from the plate’s edge closer to the center, the critical buckling load decreases by 26%. This research contributes to the advancement of understanding FGM thin plates’ behaviour under skew loading conditions, with implications for the design and optimization of innovative structures. The findings presented provide a foundation for further exploration of advanced composite materials and their applications in structural engineering.

Digital twin, an innovative technology propelled by data and models, plays a seminal role in the digital transformation and intelligent upgrade of ships. This study introduces a digital twin methodology for real-time monitoring of ship structure deformation fields, based on finite discrete strain data, and a visualization tool framework is developed using virtual reality technology. First, the inverse Finite Element Method (iFEM) is employed to derive the deformation field of the ship structure in real time via sensor strain data. Secondly, the deformation field data obtained based on the iFEM algorithm is converted into general visualization data conducive to interpretation within VR applications. Lastly, a digital twin software tool is built to enable synchronous response and interaction between the virtual scene and the physical scene, directly superposing particular virtual objects (data acquired by sensors, CAD virtual models, and deformation field cloud images) onto the physical scene in real-time. The digital twin tool embodies a virtual reality visualization system framework integrating physical data measurement, reconstruction, analysis, expression, storage, rendering and interaction of deformation field data. Through practical application, the flexibility, effectiveness and compatibility of the developed prototype tool are verified. According to the results, the system can enhance the efficiency of scientific communication, model validation, and interdisciplinary sharing during the analysis and evaluation of mechanical properties of ship structures.

The paper presents the numerical and experimental results of research aimed at determining the influence of hardness in the range of 50-90 Shore A of layered tools made of elastomeric materials on the possibility of forming Inconel 625 nickel-based alloy sheets. A stamping die made of 90MnCrV8 steel (hardness 60HRC) was designed for forming embosses in drawpieces, ensuring various stress states on the cross-section of the formed element. The principle of operating the stamping die was based on the Guerin method. Finite element-based numerical modelling of the forming process for various configurations of polyurethane inserts was also carried out. The drawpieces obtained by sheet forming were subjected to geometry tests using optical 3D scanning. The results confirm that in the case of forming difficult-to-deform Inconel 625 Ni-based alloy sheets, the hardness of the polyurethane inserts significantly affects the geometric quality of the obtained drawpieces. The assumptions made in numerical simulations were verified in experimental studies. Based on the test results, it was concluded that the selection of polyurethane hardness should be determined by the shape of the formed element. Significant non-uniform sheet metal deformation was also found, which may pose a problem in the process of designing forming tools and technology of plastic forming of Inconel 625 Ni-based alloy sheets.

Modular high-rise buildings frequently adopt cast-in-place concrete to ensure building integrity. To reduce the material consumption and realize mold-free casting, a 30 mm thick light concrete mold shell wall panel ( MiC wall panel) is set as part of the shear wall. But the concrete casting process is prone to cracking and rising pattern. Its force performance and damage mechanism is not clear. To this end, this paper designs three groups of MiC wall panels with different material component types. The static load test is carried out in a graded partitioned loading mode. The effects of different material components on the mechanical properties of MiC wall panels are systematically investigated and parametrically analyzed by numerical simulation. The results show that the front cracks of the wall slab are mainly distributed under the truss tendons, and the back cracks are mainly distributed in the span position of the adjacent truss tendons. They are both along the span direction of the wall slab. Increasing the strength level of concrete and using steel fiber concrete can improve the mechanical properties of MiC wall panels to different degrees. Among them, increasing the concrete strength has a significant effect on improving the load-bearing capacity of wall panels, while incorporating steel fibers can significantly improve the deformation and crack resistance of wall panels. By building more finite element models, it was verified that the steel fiber admixture has limited effect on the mechanical properties of wall panels as the concrete grade increases. This study aims to summarize the types of wall panels with the best working properties. It provides a practical basis for optimizing the production and fabrication process of MiC wall panels.

Compared with ordinary asphalt pavement, the application of modified asphalt will improve the performance of asphalt pavement. However, using the same modification method for the whole road section can cause modification waste or result in some parts not meeting performance standards. This paper proposes a whole road section modification method. First of all, a finite element model was established using Abqus according to the actual road pavement material to analyze the mechanical response of each layer under different working conditions and establish the grading standard under comprehensive working conditions. Secondly, different modification methods were analyzed according to investigation and their modification effects were graded through clustering analysis method. Finally, different working conditions and different modification methods were reasonably matched and a section of highway was selected for verification. The results show that before modification, the calculated rutting is 6.23mm and after modification, it is reduced to 2.45mm which indicates that reasonable modification has a better effect and can achieve the purpose of rutting equalization control on the whole road section. It is of great significance to carry out reasonable modification on different conditions to control rutting equalization of asphalt pavements.

Nickel-based Superalloy Inconel 718 is widely used in the aerospace industry for its excellent high-temperature strength and thermal stability. However, milling Inconel 718 presents challenges due to significantly increased cutting force and vibration, which is a typical difficult-to-machine material. This paper focuses on the milling process of Inconel 718, establishing a milling force model to analyze the force trends under various processing parameters. Finite element analysis is employed to study the stress and temperature fields during milling. Dynamic equations for milling Inconel 718 are developed, and stability lobe diagrams are generated based on modal experiments. Milling experiments on Inconel 718 validate the milling force model and finite element analysis results. The fmincon optimization algorithm is utilized to identify the optimal machining parameters for Inconel 718. Through this research, valuable insights into enhancing the efficiency and quality of Inconel 718 machining are provided. This study contributes to a deeper understanding of Inconel 718 milling behavior, offering crucial guidance for more efficient machining processes.

The cutter structure and layout scheme of PDC (Polycrystalline Diamond Compact) bits are important factors in improving efficiency. To further improve the drilling efficiency of PDC bits, axe, triangular prism, and cylindrical PDC cutters were used as research objects. Based on the measured granite data, a finite element model of non-homogeneous granite was established and verified by uniaxial compression simulation. The rock-breaking process of different cutter combination schemes was compared using the finite element method, and the parameters in the best scheme were optimized using the Box-Behnken response surface method. The results show that the constructed model of non-homogeneous granite is consistent with the stress-strain relationship of real granite and is reliable. The axe PDC cutter is more aggressive than the other two cutters and is more suitable for the front row of the bit blade tooth arrangement, while the triangular prism cutter is the second most aggressive and is suitable for the rear row of tooth arrangement, and the best combination scheme is the front row of axe cutter and the rear row of triangular prism cutter arranged alternately with the axe cutter. The optimal transverse and longitudinal distances of the optimized triangular prism cutter from the front axe cutter are 10mm and 7mm, and the optimal transverse and longitudinal distances of the rear axe cutter from the front cutter are 10.06mm and 7mm. The drilling speed is more stable during drilling and the PDC bit with mixed tooth arrangement has 16.8% and 16.6% higher rate of penetration(ROP) compared with the bit with single axe cutter and triangular prism cutter, and the drilling speed is more stable during working, which can effectively improve the rock breaking efficiency of the PDC bit. The field application proves that the bit with mixed cutter arrangement is easier to break the complex formations, with more stable ROP and better efficiency. The study can provide theoretical support for the cutter layout of the PDC bit.

In the laboratory practice, several standards for testing the bending stiffness of corrugated board are used. There are often cases of tests where the results depend on the way the sample is placed on the supports. The problem arises when the board is five-ply (with two corrugated layers with different corrugation heights) or when the board has asymmetrically selected papers on the flat layers. In this article, we focus our attention on the problem related to boundary conditions, with particular attention to the local effects of the support of the sample. Since the cardboard layers, both flat and corrugated, have a small thickness, a slight deformation of the papers can always be observed at the point of contact between the sample and the support, which affects the readings of the measured stiffness. The paper presents theoretical and numerical analyzes showing how much the method of supporting the sample affects the measured bending stiffness of various samples. Numerical observations were also compared with the results of analyzes presented by other scientists as well as with experimental results.

Nickel superalloys occupy a reasonably broad niche in the industry. One of this group's most common and well-known alloys is Inconel 718. Parts from Inconel 718 are used in aerospace, energy, automotive, and some other vital industries, which creates a demand for research and improvement of the machining conditions of this alloy. This article is devoted to the systematization and generalization of the accumulated experience of machining of Inconel 718 at turning operation. The research methodology is based on the Pareto distribution. In the study, more than 50 articles devoted to turning Inconel 718 were analyzed. The tool materials, methods, and types of coating used in turning Inconel 718 were analyzed. Finite element method modeling of processing processes has been considered. The trends in the selection of cutting parameters and the geometry of the cutting tool and finite element method modeling of processing processes have been considered. Lacuna for promising future research was formulated.

The aim of this paper is to optimize thickness variation function of simply supported and cantilever beams, in terms of maximizing gaps between chosen neighboring frequencies and analyze obtained results. The optimization results are examined in terms of achieving the objective function (related to eigenproblem), but also in terms of their dynamics stiffness (forced vibrations excited by point harmonic load). In the optimization process, a Genetic Algorithm was used. Problems related to structural dynamics were solved by FEM implementation into the algorithm. Sample results were presented, and the developed algorithm was analyzed in terms of result convergence by examining several variable parameters. The authors demonstrated the validity of applying the described optimization tool to the presented problems. Conclusions were drawn regarding the correlation between stiffness and mass distribution in the optimized beams and the natural frequency modes in terms of which they were optimized.

The analysis concerned a thin-walled plate with a circular central hole made of CFRP (Carbon Fiber Reinforced Polymer) material. The test object was subjected to a tensile test to investigate the cracking mechanism of the composite structure. The tests were performed using two independent test methods experimentally and numerically. In the experiment, the force was recorded as a function of elongation of the tensile plate and the elongation of the plate was investigated using an Aramis non-contact optical measurement system. As part of the numerical analysis, an FEM model adequate to the experimental conditions was performed in the Abaqus system. The fracture process was modeled using the extended finite element method XFEM. This approach enabled the structure of the composite to be investigated over the full range of tensile loading. The behaviour of the circular hole plate before the damage symptoms was investigated, the damage initiation force was determined, then the initial cracking and delamination of the layers were analysed, as well as the crack propagation process leading to the complete failure of the composite structure (cracking of all laminate layers).

This study employs a high-fidelity numerical framework to determine the plastic material flow patterns and temperature distributions that lead to void formation during friction stir welding (FSW), and to relate the void morphologies to the underlying alloy material properties and process conditions. Three aluminum alloys, viz., 6061-T6, 7075-T6, and 5053-H18 were investigated under varying traverse speeds. The choice of aluminum alloys enables investigation of a wide range of thermal and mechanical properties. The numerical simulations were validated using experimental observations of void morphologies in these three alloys. Temperatures, plastic strain rates, and material flow patterns are considered. The key results from this study are: (1) The predicted stir zone and void morphology are in good agreement with the experimental observations, (2) The temperature and plastic strain-rate maps in the steady-state process conditions show a strong dependency on the alloy type and traverse speeds, (3) The material velocity contours provide a good insight into the material flow in the stir zone for the FSW process conditions that result in voids as well as those that do not result in voids. The numerical model and the ensuing parametric studies presented in this work provide a framework for understanding material flow under different process conditions in aluminum alloys, and potentially in other alloys. Furthermore, the utility of the numerical model for making quantitative predictions and investigating different process parameters to reduce void formation is demonstrated.

Microstructure evolution in metal additive manufacturing (AM) is a complex multi-physics and multi-scale problem. Understanding the impact of AM process conditions on the microstructure evolution and the resulting mechanical properties of the printed part is an active area of research. At the meltpool scale, the thermo-fluidic governing equations have been extensively modeled in the literature to understand the meltpool conditions and the thermal gradients in its vicinity. In many phenomena governed by partial differential equations, dimensional analysis and identification of important dimensionless numbers can provide significant insights into the process dynamics. In this context, a novel strategy using dimensional analysis and the method of linear least squares regression to numerically investigate the thermo-fluidic governing equations of the Laser Powder Bed Fusion AM process is presented in this work. First, the governing equations are solved using the Finite Element Method, and the model predictions are validated by comparing with experimentally estimated cooling rates, and with numerical results from the literature. Then, through dimensional analysis, an important dimensionless quantity - interpreted as a measure of heat ab- sorbed by the powdered material and the meltpool, is identified. This dimensionless measure of heat absorbed, along with classical dimensionless quantities such as Peclet, Marangoni, and Stefan numbers, is used to investigate advective transport in the meltpool for different alloys. Further, the framework is used to study the variations of thermal gradients and the solidification cooling rate. Important correlations linking meltpool morphology and microstructure evolution related variables with classical dimensionless numbers are the key contribution of this work.

The cuttlebone-like structure is a complex porous bionic structure with asymmetric sinusoidal S-shaped wall structure connecting laminar septa, and studies have shown that the cuttlebone-like structure has light weight, high strength and excellent energy absorption capability. It has become an important research to investigate the mechanical mechanism of cuttlebone-like structure and to design a bionic structure beyond nature structure based on it. In order to investigate the mechanical properties of the static compression and shearing processes of the cuttlebone-like structure, this paper establishes theoretical formulas and analyzes the influence of dimensionless parameters such as height-to-thickness ratio $\lambda$, period-to-thickness ratio $\xi$ and amplitude to period ratio $\eta$ on the compression and shearing processes with the help of ABAQUS finite element software.The parameter sensitivity analysis method was used to compare the important influence degree of each dimensionless parameter on the mechanical properties of the structure and to determine the relative optimal parameters of the structural mechanical properties under compression and shearing. Based on the relatively optimal parameter of the cuttlebone-like structure, the cross-sectional geometry of the structure is improved to build a new cuttlebone-like structure.The experimental results showed that the new cuttlebone-like structure with revised cross-section improved the mechanical properties of the cuttlebone-like structure.The study results are of theoretical guidance to enhance the mechanical properties and geometric design of cuttlebone-like structures.

Background: Cervical laminoplasty is a useful for treatment of cervical myelopathy. However, this procedure has limitations for kyphotic cervical alignments. We used the finite element (FE) analysis and investigated the biomechanical changes in an intact and laminoplasty models with lordosis, straight, and kyphosis cervical alignments. Methods: A three-dimensional FE model of the cervical spine (C2-C7) with ligaments was created from computed tomography. The model was modified with the following cobb angles and the C3-C6 laminoplasty was conducted; a) laminoplasty-lordotic model (LM-L; C2-C7 angle: -10°), b) laminoplasty-straight model (LM-S; C2-C7 angle: 0°), and c) laminoplasty-kyphotic model (LM-K; C2-C7 angle: 10°). A pure moment with a compressive follower load of 100N to represent the weight of the head/cranium and cervical muscle stabilization was applied to these models. The range of motion (ROM), annular stress, nucleus stress and facet forces were analyzed. Results: ROM of LM-K increased when compared to the other models except for flexion. The LM-K had the highest mobility with 49% increase in ROM observed under extension, compared to the intact model. In all motion except for flexion, LM-L models’ ROM decreased by more than 10%, and LM-S and LM-K models’ ROM increased by more than 10% at C2-C7 compared to the intact model. The annular stresses and nucleus stresses in LM-K were higher compared to the other models. The maximum increase in annular stresses was about 194% in LM-K compared to the intact model, observed at the C3-C4 segment. The facet contact forces were lowest in the LM-K, compared to the other models.Conclusions: Patients with a cervical kyphosis alignment are at a disadvantage of increased kyphosis compared to cases with lordosis or straight alignment and should be treated with caution.

Gol'denveizer's problem of a torus has been analyzed by Audoly and Pomeau (2002) and Sun (2021). However, all of the investigations of Gol'denveizer's problem of an elastic torus have been linear. In this paper, the finite element method is used to more accurately address this problem. Furthermore, Sun (2021) cannot be solved by nonlinear analysis. We research the nonlinear mechanical properties of Gol'denveizer's problem of circular and elliptic tori, and relevant nephograms are given. We study the buckling of Gol'denveizer's problem of an elastic torus, and propose failure patterns and force-displacement curves of tori in the nonlinear range. Investigations reveal that circular tori have more rich buckling phenomena as the parameter a increases. Gol'denveizer's problem of the buckling of an elliptic torus is analyzed, and we find a new buckling phenomenon called a "skirt." As a/b increases, the collapse load of an elliptic torus of the Gol'denveizer problem is enhanced gradually.

Linear elastic fracture modeling coupled with empirical material tension data result in good quantitative agreement with experimental measurements of fracture failure for both brittle and tough epoxy nanocomposites. The nanocomposites comprise diglycidyl ethers of bisphenol A cured with O,O’ bis (2-aminopropylpropylene glycol) (Jeffamine D230) and doped with rubber nanoparticles of varying concentrations. Toughness, critical load, and critical displacement in quasi-static single edge-notched three-point bending are predicted accurately using both surface-based cohesive zone (CZS) and extended finite element (XFEM) methods implemented in Abaqus software. Fracture initiation within a crack is taken at the yield stress from uniaxial tension data. Prediction of fracture processes using a generalized truncated linear traction-separation law was significantly improved by considering the case of a linear softening function. There are no adjustable parameters in the XFEM model. The CZS model requires only optimization of the element displacement at fracture parameter. Thus, these continuum methods describe these materials in mode I fracture with a minimum number of independent parameters.

Fatigue analysis plays a vital role in determining the structural integrity and life of a dental implant. With the use of such implants on the rise, there is a corresponding increase in the number of implant failures. As such, the aim of this research paper is to investigate the life of 3D-printed dental implants. The dental implants considered in this study were 3D printed according to the direct metal laser sintering (DMLS) method. Additionally, a finite element model was developed to study their performance, while fatigue life was predicted using Fe-Safe software®. The model was validated experimentally by performing fatigue tests. The life of the dental implants was analysed based on Normal strain and the Brown-Miller with Morrow mean correction factor algorithm. The model revealed that there was a strong correlation between the FEA and the experimental results. The clinical success of 3D-printed dental implant experimentally is 20.51 years and computationally under Normal strain is 19.89 years and Brown-Miller with Morrow mean correction factor is 26.82 years.

EIT-MESHER (https://github.com/EIT-team/Mesher) is C++ software, based on the CGAL library, which generates high quality Finite Element Model tetrahedral meshes from binary masks of 3D volume segmentations. Originally developed for biomedical applications in Electrical Impedance Tomography (EIT) to address the need for custom, non-linear refinement in certain areas (e.g. around electrodes), EIT-MESHER can also be used in other fields where custom FEM refinement is required, such as Diffuse Optical Tomography (DOT).

Lightning strike can cause a considerable damage in aircraft parts made from semiconducting materials such as Carbon Fiber Reinforced Plastics (CFRPs). Therefore, in recent years, the lightning strike phenomenon has attracted the interest of the academic community and the aircraft industry. Until now, the problem has been addressed mainly experimentally, while the reported numerical works are very limited. In the present work, a coupled electro-thermal FE model has been developed using the ANSYS commercial FE code to simulate the lightning strike damage in unidirectional CFRP laminates due to the Joule heat flux phenomenon. The model is based on the SOLID69 thermoelectric element and applies a non-linear, time-transient analysis. The main input to the model is the thermal-electrical properties of the composite material which vary with temperature. Using the model, a parametric study on the effect of mesh density and peak intensity on the thermal damage has been performed. Three electrical lightning strikes of low (10 kA), medium (30 kA) and high peak intensity (40 kA) have been applied according to the SAE ARP 5412 standard. The electro-thermal model has been validated against a numerical model from the literature. The numerical results reveal that the increase of peak intensity leads to the increase of the area and penetration depth of matrix thermal damage (pyrolysis) as well as to the increase of the area of fiber damage (deterioration and ablation). Through progressive damage modeling, the residual tensile strength of the CFRP plate after being subjected to lightning strike of different peak intensity has been predicted. Lightning strike initial damage has been simulated by translating the thermal field into degradation of elastic properties of the lamina. The results show an increase in the accumulated matrix damage and a decrease of tensile strength due to the initial lightning strike damage. For the maximum peak intensity of 40 kA, a decrease in tensile strength of 4.8% has been predicted

The purpose of this work is to optimize angular extrusion for sever plastic deformation of tubular specimens. A finite element model (FEM) was built in Deform-3D^® for three extrusion dies and analyzed. The dies were ECAP (135°), TCAP with external groove of 90° (denoted as TCAP-e) and TCAP with internal groove of 90° (denoted as TCAP-i). The analysis for process parameters (such as coefficient of friction-μ, die angles-ѱ and φ, temperature - T, and radius ratio-R) showed that TCAP-i was the optimal die in processing Al6063 tubes based on strain and load distribution of the model. The TCAP-i die model was further optimized for different parameters namely die angles, coefficient of friction, deformation ratio and temperature. The results showed that at constant process temperature of 25 °C, the optimal TCAP-i has the following parameters: φ₂=80⁰, Ѱ₁=30⁰, Ѱ₂=80⁰, Ѱ₃=30⁰, μ =0.4 and R =1.5.

This study aimed to establish and verify the validity of an acoustic simulation method during sustained phonation of the Japanese vowels /i/ and /u/. The study participants were six healthy adults. First, vocal tract models were constructed based on computed tomography (CT) data, such as the range from the frontal sinus to the glottis, during sustained phonation of /i/ and /u/. Next, cylindrical shapes virtually extended by 12 cm were added to the vocal tract models to imitate the trachea between the tracheal bifurcation and lower part of the glottis. The Kirchhoff–Helmholtz integral equation was formulated as the wave equation for sound propagation, and the boundary element method was used for discretization. As a result, the relative discrimination thresholds of the vowel formant frequencies for /i/ and /u/ against actual voice were 1.1%–10.2% and 0.4%–9.3% for the first formant and 3.9%–7.5% and 5.0%–12.5% for second formant, respectively. In the vocal tract model with nasal coupling, a pole–zero pair was observed at around 500 Hz, and for both /i/ and /u/, a pole–zero pair was observed at around 1000 Hz regardless of the presence or absence of nasal coupling. These findings demonstrated that /i/ and /u/ could be simulated with high validity in a vocal tract model constructed from CT data obtained during sustained phonation using the boundary element method.

The vibration response of soil is a key property in the field of agricultural soil tillage. Vibration components of tillage machinery are generally used to reduce tillage resistance and improve work efficiency, the pressure variation under low-frequency vibration will affect the fragmentation and dispersion of farmland soil. However, the gradient of pressure variation, frequency domain response, and effective transmission range are unclear. A new method based on DEM (Discrete Element Method) is presented to study the vibration response and pressure transmission under low-frequency vibration. Bench test results showed that peak pressure positively correlates with the vibration frequency and attenuates rapidly at the vibration distance of 100 to 250 mm. The results data were also selected to determine the simulation model parameters. Amplitude, vibration frequency, and soil depth were used as test factors in single-factor simulation tests, and their effects on the peak pressure, frequency domain response, and effective transmission distance were analyzed. The results showed a positive relationship between the peak pressure and the test factors. The peak pressure increases with a maximum gradient of 19.02 kPa/mm at a vibration distance of 50 mm. The amplitude, vibration frequency, and soil depth positively correlate with the dominant frequency amplitude. The main frequency is independent of amplitude and soil depth. At the vibration distance of 250 mm, the dominant frequency is approximately twice the vibration frequency at 7–11 Hz and approximately equal to the vibration frequency at 13–15 Hz. Multiple exponential functions were used to fit the peak pressure attenuation function, obtaining an effective transmission distance range of 347.15 to 550.37 mm for the 5kPa cut-off pressure. For a soil depth of 300 mm, the vertical shear wave diffusion angle is greater than the horizontal shear wave diffusion angle. The study clarifies the vibration response of soil under low-frequency vibration, which helps to design vibration type soil-engaging components of tillage machinery and match vibration parameters for energy saving and resistance reduction in soil tillage.

This paper reviews numerical methods used to simulate desiccation cracks in clayey soils. It examines five numerical approaches: Finite Element (FEM), Lattice Boltzmann (LBM), Discrete Element (DEM), Cellular Automaton (CAM), and Phase Field (PFM) Methods. The FEM is widely used to capture moisture diffusion, shrinkage, and cracking during drying. LBM is used to simulate fluid flow in clayey soils, while the DEM focuses on capturing the behavior of individual particles and their interactions. CAM simplifies crack evolution with computational efficiency, while PFM provides a continuous representation of crack formation and propagation. The author discusses the complexity of the problem, the continuum mechanics governing and constitutive equations that describe it, and the influence of various factors such as the multiphase nature of soils, heterogeneity, nonlinearities, coupling, scales of analysis, and computational aspects. The review highlights the characteristics, strengths, and limitations of each method. It emphasizes the importance of appropriate method selection for every problem depending on the aim of the analysis. The article concludes by reviewing the integration of multiple numerical methods to enhance the simulation of desiccation cracks in clayey soils.

Owing to the risk of atypical femoral fractures, total hip arthroplasty presents unique challenges for patients with ipsilateral osteoarthritis and localized periosteal thickening in the femoral subtrochanteric region. Stem length selection is critical for minimizing stress concentration in the thickened cortex to avoid such fractures. Herein, we report the case of a 78-year-old woman with ipsilateral hip osteoarthritis and localized subtrochanteric periosteal thickening. Preoperative planning included finite element analysis to assess stress distribution across various stem lengths. Simulation aimed to determine the optimal stem length to span the cortical thickening and reduce the risk of postoperative complications. The results of finite element analysis indicated that a stem length of >150 mm was required to effectively reduce the stress at the site of cortical thickening. A 175-mm stem was selected for total hip arthroplasty, which provided a favorable stress distribution and avoided the risk of stress concentration. Therefore, in cases of ipsilateral osteoarthritis with localized subtrochanteric periosteal thickening, finite element analysis can be valuable for preoperative planning to determine the optimal stem length, thereby reducing the risk of atypical femoral fractures. Further studies with multiple cases are recommended to validate these findings and improve surgical outcomes.

Corrugated boards are used for packaging because of their high strength-to-weight ratio, recyclability, and biodegradability. They have an orthotropic sandwich structure with unique characteristics for each direction owing to their flute shape. In this study, the bending behavior was qualitatively analyzed for various variables through an FE-based simulation, and the possibility of an alternative test method for the four-point bending test on a corrugated board was examined through a similarity analysis with the experimental results. The cross-machine direction (CD) bending stiffness through finite element analysis (FEA) simulation was closely related to the thickness of the corrugated board. In AB-flute-double-wall (AB/F-DW), the difference in CD bending stiffness based on the bending direction was approximately 17.2% in the FEA simulation and approximately 10.5% in the experiment. However, the differences in the bending behavior (bending force vs. deflection plot) and bending stiffness based on the phase shift between the two flutes constituting BB-flute-double-wall (BB/F-DW) were insignificant. Overall, the CD bending behavior of the corrugated board was simulated relatively well through FEA simulation. However, it was impossible to simulate the machine-direction (MD) bending behavior through FE-based simulation because of a decrease in convergence and a large error caused by the variability of the contact condition of the modeled test specimen with a non-uniform MD cross-section.

To improve the current density distribution and electromagnetic performance of air-core inductors, a structural optimization method combining BP neural network and genetic algorithm was proposed for the study of axial and radial spiral multi-winding inductors. The Monte Carlo method was used to extract the structural size samples of the inductors, and the training dataset was obtained through finite element calculation of electromagnetic fields. Based on BP neural networks, nonlinear mapping models between the inductance value, volumetric inductance density, current distribution non-uniformity coefficient, and inductor structural parameters were constructed. The sensitivity analysis of the inductor inductance value affected by structural parameters was conducted through the Sobol index calculation. Using the current distribution non-uniformity coefficient as the fitness function and the volumetric inductance density as the constraint condition, the genetic algorithm was applied to optimize the structural parameters of the inductor globally. The optimization results were verified through finite element comparison. The results show that, under the requirement of satisfying the volumetric inductance density, the current distribution non-uniformity coefficient of AHI type inductor was reduced by 4.57% compared with the best sample in the sampling, while that of RHI type inductor was reduced by 5.33%, demonstrating the practicality of the BP-GA joint algorithm in the structural optimization design of inductors.

Accurate cloud quantification is essential in climate change research. In this work, we construct an automated computer vision framework by synergistically incorporating deep neural networks and finite sector clustering to achieve robust whole sky image-based cloud classification, adaptive segmentation, and recognition under intricate illumination dynamics. A bespoke YOLOv8 architecture attains over 95% categorical precision across four archetypal cloud varieties curated from extensive annual observations (2020) at a Tibetan highland station. Tailor-made segmentation strategies adapted to distinct cloud configurations, allied with illumination-invariant image enhancement algorithms, effectively eliminate solar interference and substantially boost quantitative performance even in illumination-adverse analysis scenarios. Compared with the traditional threshold analysis method, the cloud quantification accuracy calculated within the framework of this paper is significantly improved. Collectively, the methodological innovations provide an advanced solution to markedly escalate cloud quantification precision levels imperative for climate change research, while offering a paradigm for cloud analytics transferable to various meteorological stations.

The paper delves into the critical issue of damage detection within the guy cable of a truss steel mast, a pivotal component in structural integrity. It introduces a model wherein damage manifests as a localized cross-section reduction in a single element of the cable. Employing Discrete Wavelet Transform (DWT), a pioneering methodology, the study scrutinizes the behavior of the affected elements through static structural analyses. Signal decomposition via the Mallat pyramid algorithm facilitates comprehensive examination. Static displacements at the connection point between the cable and the truss mast serve as the measured variables. Through systematic investigation, the paper evaluates the impact of damage size and location, external loading force, and cable tension force on the efficacy of the proposed approach. Utilizing data derived from Finite Element Method (FEM) computations for wavelet analysis the authors substantiate the findings with numerical examples, thus offering valuable insights into damage detection strategies for structural health monitoring and engineering applications.

This work aims at improving surface integrity and LCF life of machined Inconel 718. The correlation between the LCF life of Inconel 718 and various states of machined surface integrity is explored. In this paper, the surface integrity of Inconel 718 specimens is enhanced by low plasticity burnishing (LPB). The LCF life with different surface integrity specimens is predicted by micro-macro finite element method (FEM). Firstly, the LCF specimens with different surface integrity are machined by turning (turned specimen), polishing (matrix specimen) and LPB process (LPBed specimen). Secondly, the LCF experiment is carried out to reveal the effect of surface integrity on LCF life. Finally, the LCF micro-macro FEM model is proposed to predict the LCF of machined Inconel 718 specimens. The representative volume element (RVE) model is established based on the measured surface integrity and microstructure of Inconel 718 specimens. The effect of surface integrity on LCF life is transformed to equivalent load. The micro-macro FEM model combined with Tanaka-Mura dislocation crack initiation theory and extended finite element method (XFEM) is applied to predict the LCF life the machined specimens. The study results show that the LCF life of LPBed specimens can be improved by 90.5% and 36.1% compared with that of turned specimens and polished matrix specimens, respectively. The errors between FEM prediction results and experimental results are 13.1%, 9.2% and 12.2%, respectively. The proposed micro-macro FEM model could be utilized to predict the LCF life of Inconel 718 with different surface integrity, and further to apply the LCF life prediction in industry.

In order to further investigate the study of the effect of grouped stud effect on the force properties of stud connector, based on the premise that the correctness of the finite element simulation method in this paper, finite element modeling of grouped stud connectors was developed, the grouped stud effect and its sensitivity factors were analyzed, to validate the recommended formula for calculating the shear capacity of grouped stud connectors. Results show that the number of grouped stud rows and row spacing have a significant effect on the grouped stud effect, and the unevenness coefficient of grouped stud force Ⅰ is negatively correlated with the number of grouped stud rows as well as the grouped stud row spacing. For grouped stud connectors with commonly used concrete grades greater than C50 and height-to-diameter ratios greater than 4 in steel-concrete composite structural bridges, which is insensitive to changes in the concrete strength grades and the length of the studs. The direction of force transmission for the grouped stud changes with the change of loading angle, and the unevenness coefficient of force for the grouped stud will therefore be reduced. By comparing the results of the 62 existing groups of grouped stud connectors push-out tests, the mean value of the ratio of tested to calculated values is found to be 1.12, the variance is 0.023, the dispersion is small, it is shown that the recommended formula has a high degree of accuracy. The results of this paper can be used as a theoretical basis for the study of the shear stress performance of grouped stud connectors.

For enhancing the efficiency of cantilever casting the construction in the bridge, a novel rhombus traveling track of hanging basket was developed and optimized. The civil software titled Midas was applied to analyse the mechanical properties of the hanging basket with full load and control load. The strength and distortion of the walking mechanism were within the specified range when the maximum beam unit internal force of 149.69Mpa, which was less than the allowable stress of steel 175Mpa, and the anti-overturning safety factor of the hanging basket was 2.5 to meet the requirements. Through the comparative analysis of the key components and the finite element calculation, it was found that there was about 30% redundancy in the structure performance. Therefore, further optimization of each structure was carried out, and the front elevation of the main truss of 20° was achieved. It obtained the best performance in all aspects.

Nitrile butadiene rubber (NBR) and its various composites are extensively studied as matrix materials for oil-free lubricated mechanical friction surfaces. While numerous macroscopic friction experiments can compare the frictional performance of different NBR composites, the complex material preparation and low experimental efficiency pose challenges. Nanoscale molecular dynamics (MD) simulations can help develop materials with improved friction properties, but they cannot directly provide the friction coefficients of actual macroscopic friction surfaces. Therefore, it is important to integrate the advantages of macroscopic and nanoscale friction studies and perform a synergistic analysis to modify the friction properties of composite materials. In this study, a multi-scale approach is proposed to simulate the frictional characteristics of carbon nanotube (CNT)-reinforced NBR by combining the MD friction simulation， micromechanical bridging model, and macroscopic finite element analysis (FEA) methods. The results of the multi-scale friction simulation of copper (Cu)-CNT/NBR composites show that the addition of CNTs significantly improves the frictional properties of NBR, and the extent of the improvement is related to the mass fraction of CNTs. As the mass fraction of CNTs increases (0%, 1.25%, 2.5%, 5%), the coefficient of friction decreases from 0.50 to 0.38. The results are validated by ring-block friction experiments at the same scale as the simulations.

Functionally graded materials (FGM) have attracted considerable attention in the field of composite materials and rekindled interest in research on composite materials, due to their unique mechanical response achieved through material design and optimization. Compared to conventional composites, FGMs offer several advantages and exceptional properties, including improved deformation resistance, improved toughness, lightness properties, and excellent recoverability. This study focused on the production of functionally graded (FG) polymer materials by additive manufacturing (AM) method. FG structures were produced by the Fused Deposition Modeling (FDM) method using Acrylonitrile benzidine styrene (ABS) and Polylactic acid (PLA) materials and tensile tests were performed according to ASTM D638. The effects of different layer thicknesses, volume ratios, and total thicknesses on mechanical behavior were investigated. In addition, tensile tests were carried out by producing single-layer samples at certain volume ratios to create a numerical model with the Finite element method to verify the experimental data. As a result of this study, it is presented that the FG structure produced with FDM improves mechanical behavior.

The article presents a simplified method for determining the strength of a corrugated board packaging subjected to dynamic transport loads. The proposed algorithm consists of several calculation steps: (1) a static analysis of the compressive strength of the package, (2) an analysis of random vibrations in the frequency domain, used to determine the resonance frequencies, and (3) a dynamic analysis of the package loaded with computed resonant frequencies. For this purpose, numerical models of the static compression test of the packaging before and after the dynamic analysis of the package subjected to general transport loads were developed. In order to validate the model, the laboratory packaging compression tests were also performed for samples of boxes with 3-layer cardboard. Thanks to this, it was possible to verify the results of numerical simulations of compression tests for several box geometries. Which, in turn, allowed for the development of a method based on dynamic and post-dynamic (static) numerical analysis, allowing to determine with high accuracy the resistance of selected packaging to vibrations and dynamic loads. The results of the (validated experimentally) numerical analysis prove the usefulness of the simplified method presented here for precise estimation of the load capacity of various packages dynamically loaded during transport.

In this paper, we present the investigation results of radiation-induced effects in metal-oxide-semiconductor (MOS) doped with moderate amounts of Zinc Oxide (ZnO) as a potential candidate for space-borne application. The samples were fabricated via the sputtering method at a working pressure of 3mTorr and a deposition temperature of 300oC. The ZnO samples were exposed to 1.25-MeV gamma-ray utilizing Co60 source, and their electronic response was measured at ionizing doses ranging from 10 kGy to 300 kGy. A comparative work was performed through finite element method to simulate the electronic response of the PN junction diode due to ionizing radiation. The results indicate that the ideality factor of the MOS diode increases as the ionizing dose increases, rendering it unsuitable for use as a diode. The degradation of the electrical parameters was also simulated, showing the increase in hole concentration. These findings suggest that the ejection of electrons occurred, which agrees with the gamma radiation effects trend. Furthermore, as the intensity of radiation increases, the spatial charge that arises from the separation of hole-electron pairs results in a substantial reduction of the electric field in the central portion of the n-type region. These findings provide insights into the degradation of electrical parameters in MOS devices under gamma radiation and have implications for their use in space-borne applications.

To investigate the impact of a magnetic field on plaque development in a stenotic bifurcated artery, a finite element method is utilized. The blood flow is modelled as a stable, incompressible, Newtonian, biomagnetic, and laminar fluid. Furthermore, the arterial wall is assumed to be linear elastic. The Arbitrary Lagrangian Eulerian (ALE) method is employed to describe the hemodynamic flow in a bifurcated artery under the influence of an asymmetric magnetic field, taking into account two-way fluid-structure interaction coupling. A stable $P_2P_1$ finite element pair discretizes a nonlinear system of partial differential equations that requires a solution. The Newton-Raphson method is utilized to find a solution to the resulting nonlinear algebraic equation system. Numerical modelling is used to simulate the presence of magnetic fields, and the resulting displacement, velocity magnitude, pressure, and wall shear stresses are shown for a range of Reynolds numbers ($Re = 500$, $1000$, $1500$, and $2000$). The results of the numerical analysis demonstrate that the presence of a magnetic field has a significant effect not only on the magnitude of displacement but also on the velocity of the flow. The application of a magnetic field reduces flow separation, extends the recirculation area near the stenosis, and increases wall shear stress.

Alternating current field measurement (ACFM) testing is one of promising techniques in the field of non-destructive testing with advantages of the non-contact capability and the reduction of lift-off effects. In this paper, a novel crack detection approach is proposed to reduce the effect of the angled crack (cack orientation) by using rotated ACFM techniques. The sensor probe is composed of an excitation coil and two receiving coils. Two receiving coils are orthogonally placed in the centre of the excitation coil where the magnetic field is measured. It is found that the change of the x component and the peak value of the z component of the magnetic field when the sensor probe rotates around a crack follows a sine wave shape. A customised accelerated finite element method solver programmed in MATLAB is adopted to simulate the performance of the designed sensor probe which can significantly improve the computation efficiency due to the small crack perturbation. The experiments have also been carried out to validate the simulations. It is found that the ratio between the z and x components of the magnetic field remains stable under various rotation angles. It shows the potential to estimate the depth of the crack from the ratio detected by combining the magnetic fields from both receiving coils (i.e., the x and z components of the magnetic field) using the rotated ACFM technique.

Zero-rate output (ZRO) drift induces deteriorated micro-electromechanical system (MEMS) gy-roscope performances, severely limiting its practical applications. Hence, it is vital to explore an effective method toward ZRO drift reduction. In this work, we conduct an elaborate investigation on the impacts of the internal and packaging stresses on the ZRO drift at the thermal start-up stage, and propose a temperature-induced stress release method to reduce the duration and magnitude of ZRO drift. Self-developed high-Q dual mass tuning fork gyroscopes (TFGs) are adopted to study the correlations between temperature, frequency and ZRO drift. Furthermore, a rigorous finite element simulation model is built based on the actual device and packaging structure, revealing the temperature and stresses distribution inside TFGs. Meanwhile, the relationship between temperature and stresses are deeply explored. Moreover, we introduce a temperature-induced stress release process to generate thermal stresses and reduce the temperature-related device sensitivity. By this way, the ZRO drift duration is drastically reduced from ~2000 s to ~890 s, and the drift magnitude decreases from ~0.4 °/s to ~0.23 °/s. This stress release method achieves a small bias instability (BI) of 7.903 °/h and a low angle random walk (ARW) of 0.792 °/√h, and the long-term bias performance is significantly improved.

To understand the enhanced protection mechanism of CoCrNiAlY-YSZ-LaMgAl11O19 dou-ble-layer ceramic coating with aluminum plating, a finite element simulation method was used to simulate the distribution of thermal stress in the coating in all directions. The results show that in the air exposure of the un-aluminized coating, high temperature causes a large radial thermal stress on the surface of the LaMgAl11O19 (LMA) layer, and it increases with the increase in temperature, which is the main reason for the initiation of axial cracks. After arc aluminum plating, the aluminum plating layer effectively inhibited the volume shrinkage of the coating through good adhesion to the coating and internal diffusion, the thermal stress of the coating was considerably reduced, and the CoCrNiAlY-YSZ-LMA coating had an effective enhancement and protection effect; however, there was still a certain amount of shear thermal stress inside the LMA layer, the top of the crack, and the bottom of the crack. This thermal stress caused the initi-ation of radial microcracks in the LMA layer, which also becomes a risk point for the failure of the aluminum coating.

Variations in the implant thread shape and occlusal load behavior may result in significant changes in the biological and mechanical properties of dental implants and surrounding bone tissue. Most previous studies consider a single implant thread design, an isotropic bone structure, and a static occlusal load. However, the effects of different thread designs, bone material properties, and loading conditions are important concerns in clinical practice. Accordingly, the present study performs Finite Element Analysis (FEA) simulations to investigate the static, quasi-static and dynamic response of the implant and implanted bone material under various thread designs and occlusal loading directions (buccal-lingual, mesiodistal and apical). The simulations focus specifically on the von Mises stress, displacement, shear stress, compressive stress and tensile stress within the implant and the surrounding bone. The results show that the thread design and occlusal loading rate have a significant effect on the stress distribution and deformation of the implant and bone structure during clinical applications. Overall, the results provide a useful insight into the design of enhanced dental implants for an improved load transfer efficiency and success rate.

Corrosion is considered as one of the main factors in the structural performance deterioration of steel members. In this study, experimental and numerical methods were used to assess the reduction in compressive strength of short tubular steel columns with local corrosion damage. The corrosion damage was varied with different depths (0, 1.5, 2, 3, 4, 4.5, and 6 mm), height (0, 20, 40, 60, 80, 100, 120, 140, 160, and 180 mm), circumference (0, 90, 180, 270, and 360°), and location along the column. A parametric numerical study was performed to establish a correlation between the residual compressive strength and the severity of corrosion damage. The results showed that as the corrosion depth, height and circumference increased, the compressive strength decreased linearly. As for the corrosion height, the residual compressive strength became constant after decreasing linearly when the corrosion height was greater than the half-wavelength of buckling of the short columns. An equation is presented to evaluate the residual compressive strength of short columns with local corrosion wherein the volume of the corrosion damage was used as a reduction factor in calculating the compressive strength. The percentage error using the presented equation was found to be within 11.4%.

Mobile elevating work platforms (MEWPs) consist of a work platform, extending structure, and chassis, and are used to move persons to working positions. MEWPs are useful but are composed of pieces of equipment, and accidents do occur owing to equipment defects. Among these defects, accidents caused by the fracture of bolts fixed to the extension structure and swing system are increasing. This paper presents a failure analysis of the fixing bolts of MEWP. Standard procedure for failure analysis was employed in this investigation. Visual inspection, chemical analysis, tensile strength measurement, and finite element analysis (FEA) were used to analyze the failure of the fixing bolts. Using this failure analysis approach, we found the root cause of failure and proposed a means for solving this type of failure in the future. First, the chemical composition of the fixing bolt is obtained by a spectroscopy chemical analysis method, which determined that the chemical composition matched the required standard. The tensile test showed that the tensile and yield strengths were within the required capacity. The stress analysis was carried out at five different boom angles, and it was determined that the fixing bolt of MEWP can withstand the loads at all the boom angles. The outcomes of the fatigue analysis revealed that the fixing bolt fails before reaching the design requirements. The results of the fatigue analysis showed primarily that the failure of the fixing bolt was due to fatigue. A visual inspection of the fractured section of the fixing bolt also confirmed the fatigue failure. We propose a method to prevent failure of the fixing bolt of the MEWP from four different standpoints: the manufacturer, safety certification authority, safety inspection agency, and owner.

Globally 30% of landslides occur in the northeastern part of India [1]. One of the major earthquake events in Sikkim, India occurred on 18th September 2011 (Mw 6.9) led to over 300 landslides and 122 human deaths [2]. These landslides not only controlled by natural disasters but initiated due to human activities. The present study considered Lungchok landslide occurred in south district of Sikkim due to 2011 seismic event. The study focused on the failure mechanism of the landslide based on finite element analysis by adopting eight different cases. The deformation characteristic was investigated for dry and saturated slope conditions under static and dynamic behavior considering vehicle loads using GeoStudio software. The FEM analysis has been carried out using load deformation and linear elastic. The analysis shows that the failure of the slope was not sudden due to the 2011 earthquake event, but progressive failure was observed with time and construction activity. The paper demonstrates that, an increase in infrastructure development including construction by hill cutting increased the initiation of landslide with soil erosion. The cracks developed after 2011 earthquake event led to further deformations during future disasters required effective stabilization measures.

The finite volume method (FVM) is a kind of numerical analysis method that is simple and can take into account the topography as the finite element method (FEM). The adjoint state method (ASM) can easily and efficiently calculate the gradient vector of the objective function for the inversion. We developed a program to apply these methods to 3-D Magnetotelluric (MT) inversion and verified it using test models. In the forward modelling, it was confirmed that the program could perform analysis properly and efficiently considering topography and subsurface resistivity heterogeneity. Furthermore, we conducted the inversion using a synthetic open dataset and could reconstruct the true model of the dataset. The comparison with the previous work revealed that the computation time was comparable in spite of the fact that we calculated without MPI implementation, which the previous work used.

Our aim is to implement and test some less known numerical methods that might be used to solve an advection equation. The description is restricted to information necessary to implement and test the presented methods and it is given only for one-dimensional case.

Spacecraft recovery technology is crucial in the field of aerospace, in which the parachute plays a key role in slowing down the descent speed of the spacecraft and realizing a smooth landing. In order to construct a dynamically adjustable parachute deployment strategy, it is necessary to measure the parachute dynamic load accurately in real time. However, the existing sensor measurement scheme is difficult to meet the measurement requirements due to its large structure and complex wiring. In order to meet the current demand for real-time measurement of parachute cords dynamic load, a miniature measuring instrument is designed.According to the function and technical requirements of the miniature measuring instrument, the hardware modules of the acquisition system are selected and designed, and the integration debugging and performance optimization of the microcontroller module, A/D sampling module, signal acquisition circuit and power supply module are carried out. The software of the parachute cords tension acquisition system based on the miniature measuring instrument is developed. The elastomer is modeled by using SolidWorks and statically analyzed by using AnsysWorkbench finite element analysis software. Then the final structure of the elastomer and the pasting position of the strain gage are determined through the results analysis as well as experimental verification. The hardware module of the signal acquisition system for the miniature measuring instrument is then encapsulated. The force value of the miniature measuring instrument is calibrated and tested many times by using the microcomputer-controlled electronic universal testing machine. The experimental results show that the designed miniature measuring instrument has accurate data, strong stability and good real-time performance, which meets the demand for real-time accurate measurement of miniature measuring instrument, and can provide reliable data for parachute cords parameter validation and stepless unfolding design.

Cascading landslides represent a dynamic and hazardous geological phenomenon, char-15 acterized by a sequential chain of slope failures triggered by various factors such as heavy rainfall, 16 seismic activity, or anthropogenic activities. These events can amplify the damage caused by the 17 initial trigger and propagate instability along a slope, often resulting in significant environmental 18 and societal impacts. The Morino-Rendinara cascading landslide, situated in the Ernici mountains 19 along the border of Abruzzo and Lazio region (Italy), serves as a notable example of the complexities 20 and devastating consequences associated with such events. In March 2021, a substantial debris flow 21 event in Morino obstructed the Liri River, marking the latest step in a series of landslide events 22 characterized by a complexity far beyond initial expectations. This study employs a multidiscipli-23 nary approach, combining conventional techniques and advanced technologies, to understand the 24 complexities of the Morino-Rendinara landslide. Field activities, UAV image acquisition, SAR inter-25 ferometry based on SENTINEL-1 images and pixel offset analysis based on high-resolution Google 26 Earth images, offer insights into the geological and hydrogeological setting of the unstable slope, 27 landslide geometry, mechanism and kinematics. To supplement the analysis a specific FEM slope 28 stability analyses is used to reconstruct the deep geometry of the system emphasizing the modula-29 tion action of groundwater flow to the slope stability

This paper studies the flexural behavior of the bridge deck link slabs made with the Polyvinyl Alcohol-Strain Hardening Cementitious Composites (PVA-SHCC). Tensile and flexural properties of the PVA-SHCC with four volume fractions, i.e., 0%, 1%, 1.5%, and 2%, were evaluated first. Next, using the similarity theory, composite models with a geometric similarity ratio of 1: 5 were designed to represent the bridge deck with the link slabs. The models considered three materials for link slabs, including concrete, cement mortar and self-compacting PVA-SHCC; and two different curing ages at 7 and 56 days. Bending tests were performed to investigate the flexural behavior of the models. Based on the fractal theory, cracking characteristics of the models with different types of link slabs were compared, and the relationship between Df and the flexural behavior of the models were studied. Numerical models were built to correlate with the results from the bending tests. It is illustrated that the flexural behavior of the self-compacting PVA-SHCC link slab is better than that of concrete and cement mortar link slabs, where the crack initiation and propagation can be postponed.

Scaling up the production of power devices involves, in the semiconductor industry, the use of large semiconductor 8” and 12” wafers of Si or SiC. Their use implicates a successful an effective monitoring and control of wafer warpage, as the geometrical size of the substrate increases, as well as a command control on the degeneration of non-linear warpage into critical phenomena known as bifurcation or buckling.To mitigate the warpage, the use of taiko wafers, which consists of wafers having a thicker region on the rim region of the wafer, has been generally adopted in the production lines. In a previous work [1] we investigated the warpage of taiko wafers by introducing the concept of equivalent thickness, in the linear approximation, for the case of a front side metalized taiko wafer. However, as the thickness of the central semiconductor substrate region decreases the possibility that also the taiko wafer bifurcates emerges. For this reason, in this work we have investigated the occurrence of bifurcation or buckling in a taiko wafer. In particular, we implemented a finite element analysis method by using ANSYS Mechanical Enterprise 2023/R2 to investigate the phenomenon of bifurcation in a taiko wafer. The investigated system consisted of a 200 mm Si (001) front side metalized (FSM) taiko wafer, having a substrate thickness ranging from 100 µm to 400 µm, with a step height fixed to 450 µm. The FSM layer consisted of an Al layer of 4.5 µm thin. The system was subjected to a thermal load such that, getting cooled from a nominal temperature T to 25°C, the metal layer can develop a residual stress. To investigate the bifurcation, an asymmetry was induced in the system by applying a slight force on a point of the circumferential region perpendicularly to the substrate. The temperature investigated ranged from 30 °C to 170 °C. From this investigation, it results that also the bifurcation is mitigated by the taiko structure. Indeed, the bifurcation emerges at higher values of the stresses with respect to the case of flat wafers having the same size and thickness of the substrate. Moreover, the concept of equivalent thickness needs to be reformulated by taking into account of the bifurcation phenomenon. Indeed, from the warpage at the bifurcation point of the taiko wafer, we can determine an overall equivalent thickness of the taiko wafer.

This paper proposes an analytical model of the double-stator spoke-type permanent-magnet vernier machine (DSSTVM) using the subdomain model (SDM), which can be used to calculate the mag-netic field distribution and corresponding electromagnetic parameters of the DSSTVM. The whole field domain is divided into several subdomains according to the magnetic characteristics of each region, within which Laplace's and Poisson's equations are solved accordingly in terms of magnetic vector potential (MVP). Then, the magnetic flux density distribution, back electromotive force (EMF), and electromagnetic torque of the DSSTVM can be obtained correspondingly. Ultimately, finite element analysis (FEA) is adopted to validate the proposed analytical model's effectiveness for quickly predicting the no-load and on-load performances of the DSSTVM.

New or in-service truss bridges, with or without upper bracing systems, may display instability phenomena such as general lateral torsional buckling of the upper chord. The buckling of structural elements, particularly in the case of steel bridges, can be associated with the risk of collapse or temporary/permanent withdrawal from service. Buckling poses a significant danger as it often occurs at lower load values compared to those considered during the design phase. Additionally, this phenomenon can manifest suddenly, without prior warning, rendering intervention for its prevention impossible or futile. In contemporary times, most research and design calculation software offer the capability to establish preliminary values for buckling loads, even for highly intricate structures. This is typically achieved through linear eigenvalue buckling analyses, often followed by significantly more complex large displacement nonlinear analyses. However, interpreting the results for complex bridge structures can be challenging, and their accuracy is difficult to ascertain. Consequently, the objective of this paper is to introduce an original method for a more straightforward estimation of the buckling load of the upper chord in steel truss bridges. The method utilizes the theory of beams on an elastic foundation. The buckling load of the upper chord was determined using both the finite element method and the proposed methodology, yielding highly consistent results.

Thermoplastic composites are continuously replacing thermosetting composites in lightweight structures. However, the accomplished work on the fatigue behavior of thermoplastics is quite limited. In the present work, we propose a numerical modeling approach for simulating fatigue delamination growth and predicting the residual tensile strength of quasi-isotropic TC 1225 LM PAEK thermoplastic coupons. The approach is supported and validated by tension and fatigue (non-interrupted and interrupted) tests. Fatigue delamination growth is simulated using a mixed-mode fatigue crack growth model, which is based on the cohesive zone modeling method. All interfaces between the thermoplastic plies were modelled using cohesive elements. Quasi-static tension analyses on pristine and fatigued coupons were performed using a progressive damage model implementing a set of Hashin-type strain-based failure criteria and a damage mechanics-based material property degradation module. Using the fatigue model, delamination growth was predicted as a function of number of cycles at all interfaces. The results agree well with C-scan images taken on fatigued coupons during interruptions of fatigue tests. An unequal and unsymmetric delamination growth was predicted due to the quasi-isotropic lay-up. Moreover, the combined models capture the decrease in the residual tensile strength of the coupons. In the quasi-static tension analysis of the fatigued coupons, the driving failure mechanisms are the fast propagation of the pre-existing delamination and the severe matrix cracking.

Fabry-Pérto (FP) cavity is the essential component of ultra-stable laser (USL) for gravitational wave detection, which couples multi-physic (Optical/Thermal/Mechanics) fields and requires ultra-high precision. To satisfy the requirements of precise and efficient design, a multi-physic coupling optimization method for fixing cubic FP cavity based on data learning is proposed. A multi-physic model for the cubic FP cavity is established and the performance is obtained by finite element analysis. The key performance indices (V, wF, wF) and key design variables (d, l, F) are determined considering the features of the FP cavity. The 49 sets of data by orthogonal experiment are acquired for the establishment and comparison of different data learning models (NN, RSF, KRG). The result turns out that the neural network has the best performance. Based on NSGA-II, the Pareto optimal front is obtained and the optimal combination of design variables is finally determined as {5,32,250}. The performance after optimization has proven to be a great improvement, of which the displacement under the fixing force and vibration test are decreased by more than 60%. The optimization strategy can not only help the design of the FP cavity but also enlighten other optimization fields.

Graphyne is a material that has unique mechanical properties, but little is known about how these properties change when the material has holes. In this work, the effect of hole geometry, considering circular, triangle, and rhombus hole configurations, on the mechanical nonlinear response of γ-graphyne structures is studied. Graphyne, graphdiyne, graphyne-3, and graphyne-4 structures are under investigation. An efficient nonlinear finite element analysis (FEA) method is adequately implemented under large deformations for this purpose. The study varied the size and shape of the holes to understand how these changes affect the nanostructure's mechanical response. The results indicate that the hole geometry significantly impacts the mechanical nonlinear response of γ-graphyne structures. The holes' size and shape affect the structures' elastic behavior, deformation, and strength. The findings can be used to optimize the design of γ-graphyne structures for specific mechanical applications. The study highlights the importance of considering the hole geometries in the design and fabrication of these materials.

This work aims to describe the dielectric behaviours of four main tissues in the human forearm using mathematical modelling, including fat, muscle, blood, and bone. Determining the electrical response of multiple tissues can help develop physiological monitoring of an organ or a body’s section by filtering out the impedance contributions from the surrounding tissues to blood flow-induced impedance variations. Multi-frequency bioimpedance analysis (MF-BIA) was initially performed using the finite element method (FEM) with a 3D forearm model followed by a pilot study to characterise the response of actual forearm tissues from 1 kHz to 349 kHz. Both the simulation and experimental results were fitted to a single-dispersion Cole model (SDCM) and a multi-dispersion Cole model (MDCM) to determine the Cole parameters for each tissue. The correlation analysis and Bland-Altman plot indicated a good fit between raw and fitted impedance values using both SDCM and MDCM. Overall, MDCM exhibited better performance in fitting and estimation of the Cole parameters with correlation coefficient (R2) of 0.99 and 0.97, root mean squared error (RMSE) of 0.09 Ω and 0.14 Ω, and mean difference (mean±standard deviation) of 0.00±0.09 Ω and -0.03±0.14 Ω for the real part and imaginary part of impedance, respectively.

The ultrasonic vibration assisted pressing process can improve the fluidity, the uneven distribution of density and particle size of WC-Co powder. However, the microscopic mechanism of ultrasonic vibration on the powder remains unclear. In this paper, WC and Co particles are simulated by three-dimensional spherical models of two particle sizes with the aide of the secondary development of Python. At the same time, the forming process of powder at mesoscale is simulated in the virtue of the finite element analysis software ABAQUS. The influence of vibration amplitude on the fluidity, the filling density and the stress distribution of WC-Co powder is investigated when the ultrasonic vibration is applied to the conventional pressing process. From simulation results, the ultrasonic vibration amplitude has great influences on the density of the compact. With the increase of the ultrasonic amplitude, the compact density also increases gradually. At the initial stage of the compaction, the particles move violently under the action of ultrasonic vibration, the fluidity is obviously enhanced, the arch bridge effect between particles is destroyed, the fine particles quickly fill the pores between the large particles, which increase the density, and the residual stress in the billet decreases after the compaction. From the experimental results, the size distribution of the billet is more uniform, the elastic aftereffect is reduced, the dimensional instability is improved, and the density curves obtained by experiment and simulation are within a reasonable error range.

The inverse Finite Element Method (iFEM) is a model-based technique to compute the displacement (and then the strain) field of a structure from strain measurements and a geometrical discretization of the same. Different literature works exploit the error between the numerically reconstructed strains and the experimental measurements to perform damage identification in a Structural Health Monitoring framework. However, only damage detection and localization are performed, without attempting a proper damage size estimation. The latter could be based on machine learning techniques, however, an a priori definition of the damage conditions would be required. To overcome these limitations, the present work proposes a new approach in which the damage is systematically introduced in the iFEM model to minimize its discrepancy with respect to the physical structure. This is performed with a maximum likelihood estimation framework, where the most accurate damage scenario is selected among a series of different models. The proposed approach is experimentally verified on an aluminum plate subjected to fatigue crack propagation, which enables the creation of a Digital-Twin of the structure itself. The strain field fed to the iFEM routine is experimentally measured with an Optical Backscatter Reflectometry fiber and the methodology is validated with independent observations of lasers and the Digital Image Correlation.

The calculation of the collapse load of spherical domes is addressed using a semi-analytical approach under the hypotheses of small displacements and perfect plasticity. The procedure is based on the numerical approximation of the self-stress that represents the projection of the balance equilibrium null space on a finite dimensional manifold. The so obtained self-equilibrated stress span is superimposed to a finite element linear elastic solution to the prescribed loads yielding to the statically admissible set accordingly to Melan’s theorem. The compatibility of the stress with the constitutive law of the material has been enforced using linearized limit domain in terms of generalized stress, namely axial force and bending moment along the local spherical curvilinear coordinates. The procedure has been tested with reference to numerical and experimental data from the literature confirming the accuracy of the proposed method. The comparison with the literature confirms that the buckling load is much greater than the plastic collapse loads both calculated through the proposed procedure and reported in the quoted literature.

Cervical fusion has been a standard procedure for treating the abnormalities associated with the cervical spine. However, the reliability of anterior cervical discectomy and fusion (ACDF) has become arguable due to its adverse effects on the biomechanics of adjacent segments. One of the drawbacks associated with ACDF is adjacent segment degeneration (ASD) which has served as the base for the development of dynamic stabilization systems (DSS) and total disc replacement (TDR) devices for cervical spine. However, the hybrid surgical technique has also gained popularity recently but their effect on the biomechanics of cervical spine is not well researched. Thus, the objective of this FE study was to draw the comparison among single, bi-level and hybrid surgery with DCI implant with traditional fusion. Reduction in range of motion (ROM) for all the implanted models was observed for all the motions except extension, compared to intact model. The maximum increase in ROM of 42% was observed at C5-C6 level in Hybrid-DCI model. The maximum increase in adjacent segment’s ROM of 8.7% was observed in multilevel fusion model. The maximum von Mises stress in the implant was highest for the multilevel DCI model. Our study also showed that the shape of DCI implant permits flexion/extension relatively more compared to lateral bending and axial rotation.

A new two-dimensional nanomaterial – Titanium Carbide MXene (Ti₃C₂-MXene) – was reported in 2011. In this work, the microscale models of Ti₃C₂-MXene nanomaterial are considered with polymer matrix. The nanocomposites are modeled using nacre-mimetic brick-and-mortar assembly configurations due to enhanced mechanical properties and interlocking mechanism between the Ti₃C₂-MXene (brick) and polymer matrices (mortar). The polymer matrix material (Epoxy-resin) is modeled with elastic and viscoelastic behavior (Kelvin-Voigt Model). The Finite Element Method is used for numerical analysis of the microscale models with the multi-point constraint method to include Ti₃C₂-MXene fillers in the polymer matrix. Ti₃C₂-MXenes are considered as thick plate elements with transverse shear effects. The response of elastic and viscoelastic models of polymer matrix are studied. Finally, a tensile and compressive load is applied at the microscale and the effective load transfer due to nacre-mimetic configuration is discussed. This paper provides nacre-mimetic models to pre-design the nanocomposite for optimal performance with damage resistance and enhanced strength.

Background: Multifunctional nanoparticles are known for their wide range of biomedical applications. Controlling the magnetic properties of these nanoparticles is imperative for various applications, including therapeutic angiogenesis. The study was performed to evaluate the magnetic properties and their control mechanisms by the external magnetic field. Methods: A 100nm magnetic nanoparticle was placed in the magnetic field, and parametrically the magnet field strength and distance was evaluated. Various models of magnetic strength and disposition were evaluated. Magnetic flux density, force/weight, and magnetic gradient strength were the parameters evaluated in electromagnetic computational software. Results: The seven-coil method with three centrally placed coils as Halbach array, and each coil with a flux density of 7 Tesla, and with a coil dimension of 20cmx20cm (square model) of each coil showed a good magnetic strength and force/weight parameters in a distance of 15cm from the centrally placed coil. The particles were then evaluated for their motion characteristics in saline. It showed good displacement and acceleration properties. After that, the particles were theoretically assessed in a similar mathematical model after parametrically correcting the drag force. After application of high drag forces, the particles showed adequate motion characteristics. When the particle size was reduced further, the motion characteristics were preserved even with high drag forces. Conclusion: There is potential for a novel method of controlling multifunctional magnetic nanoparticles using high magnetic fields. Further studies are required to evaluate the motion characteristics of these particles in-vivo and invitro.

The formability and failure behavior of TRIP steel blanks were investigated through various stress states. The forming limit diagram (FLD) at fracture is constructed both experimentally and numerically. Numerical studies are performed to evaluate the applicability of different damage criteria in predicting the FLD as well as complex cross-die deep drawing process. The fracture surface and numerical results revealed that the material failed in a different mode for different strain path. Therefore, Tresca model which is based on shear stress accurately predicted the conditions where shear had the profound effect on the damage initiation, whereas Situ localized necking criterion was able to calculate the conditions which localization was dominant.

Long-term durability of the vacuum edge-seal plays a significant part in retrofitting triple vacuum glazing (TVG) to existing buildings in achieving towards zero-energy buildings (ZEB) target. Vacuum pressure decrement with respect to time between panes affect the thermal efficiency of TVG. This study reports a 3D finite element model, with validated mathematical methods and comparison, for the assessment of the influence of vacuum pressure diminution on the thermal transmittance (U value) of TVG. The centre-of-pane and total U values of TVG calculated to be 0.28 Wm⁻²K⁻¹ and 0.94 Wm⁻²K⁻¹ at the cavity vacuum pressure of 0.001 Pa. The results suggests that a rise in cavity pressure from 0.001 Pa to 100 kPa increases the centre-of-pane and total U values from 0.28 Wm⁻²K⁻¹ and 0.94 Wm⁻²K⁻¹ to 2.4 Wm⁻²K⁻¹ and 2.58 Wm⁻²K⁻¹, respectively. The temperature descent on the surfaces of TVG between hot and cold sides’ increases by decreasing the cavity vacuum pressure from 50 kPa to 0.001 Pa. To maintain the cavity vacuum pressure of 0.001 Pa for over 20 years of life span in the cavity of 10 mm wide edge sealed triple vacuum glazing, non-evaporable getters will maintain the cavity vacuum pressure that will enable the long-term durability to TVG.

Detection of early stage corrosion on slender steel members is crucial for preventing buckling failures of steel structures. An active photoacoustic fiber optic sensors (FOS) system has been reported for early stage steel corrosion detection of steel plates and rebars using surface ultrasonic waves. The objective of this paper is to investigate the surface corrosion/rust detection problem on steel rods using numerically simulated surface ultrasonic waves. The finite element method (FEM) is applied in simulating the propagation of ultrasonic waves on steel rod models. Transmission mode of damage detection is adopted, in which one source (transmitter) and one sensor (receiver) are considered. In this research, radial displacements at the receiver were simulated and analyzed by short-time Fourier transform (STFT) for detecting, locating, and quantifying a surface rust located between the transmitter and the receiver. From our time domain and frequency domain analyses, it is found that the presence, location, and dimensions (length, width, and depth) of surface rust can be estimated by ultrasonic waves propagating through the surface rust.

The Wind Energy Conversion System (WECS) plays an inevitable role across the world. In particular, the attention for Permanent Magnet Synchronous Generators (PMSGs) connected with wind farm is popular. This paper deals with the literature review that describes the recent advances, progresses and innovatory trends on PMSGs for WECS. Comparison between geared and direct-driven conversion systems and the classification of electrical machines used in WECS are discussed. A detailed analysis on the design aspects considering various topologies of PMSGs are encompassed in the literature. The PMSG design and optimization problems are solved by field computation techniques and optimized by using Soft Computing (SC) techniques .The three-dimensional, finite element software platform for the analysis and design of PMSGs is discussed. This paper also deals with the interdisciplinary modeling, analysis, and optimization of PMSG using Finite Element Analysis (FEM) and SC techniques. Finally, PMSGs are reviewed and compared for further exploration.

In the current study, the failure behavior of retrofitted steel structures was studied experimentally and theoretically with steel/CFRP double strap joints (DSJs) under quasi-static tensile loading. A series of DSJs with different bonding lengths are also considered and examined to experimentally assess the effective bond length. To predict the failure load values of the tested specimens, a new stress-based criterion, namely the point stress (PS) criterion is proposed. Although some theoretical predictive modelling for the strength between steel/CFRP joints under various loading conditions has been presented, in this work by using the new proposed approach, one can calculate rapidly and conveniently the failure loads of the steel/CFRP specimens. Furthermore, to assess the validity of the new proposed criterion, further experimental data on steel/CFRP DSJs available in the open literature are predicted using the PS criterion. Finally, it was found that a good agreement exists between the experimental results and the theoretical predictions based on the PS criterion.

An effective method to inspect inaccessible nuclear power facility by interface wave which propagate along the shrink fit boundary of reactor head is proposed in this study. Reactor head is relatively thick to inspect from the outside of reactor by conventional ultrasonic testing. The proposed interface wave can propagate a long distance from the fixed transducer position. The inside of nuclear reactor is limited to access due to the high radiation, so transducers are located at outside of nuclear facility and interface wave propagates into the nuclear reactor for defect detection. The numerical simulation and experiments were carried out to validate the effectiveness of the proposed interface wave inspection method. Various defect cases simulating field failures are also presented with satisfactory detectability by the proposed technique with the features for defect classification.

In this paper we introduce a computational model for the simulation of hypervelocity impact (HVI) phenomena which is based on the Discrete Element Method (DEM). Our paper constitutes the first application of DEM to the modeling and simulating of impact events for velocities beyond 5 kms⁻¹. We present here the results of a systematic numerical study on HVI of solids. For modeling the solids, we use discrete spherical particles that interact with each other via potentials. In our numerical investigations we are particularly interested in the dynamics of material fragmentation upon impact. We model a typical HVI experiment configuration where a sphere strikes a thin plate and investigate the properties of the resulting debris cloud. We provide a quantitative computational analysis of the resulting debris cloud caused by impact and a comprehensive parameter study by varying key parameters of our model. We compare our findings from the simulations with recent HVI experiments performed at our institute. Our findings are that the DEM method leads to very stable, energy–conserving simulations of HVI scenarios that map the experimental setup where a sphere strikes a thin plate at hypervelocity speed. Our chosen interaction model works particularly well in the velocity range where the local stresses caused by impact shock waves markedly exceed the ultimate material strength.

A computational methodology for the hydrodynamic analysis of horizontal axis marine current turbines is presented. The approach is based on a boundary integral equation method for inviscid flows originally developed for marine propellers and adapted here to describe the flow features that characterize hydrokinetic turbines. To this purpose, semi-analytical trailing wake and viscous-flow correction models are introduced. A validation study is performed by comparing hydrodynamic performance predictions with two experimental test cases and with results from other numerical models in the literature. The capability of the proposed methodology to correctly describe turbine thrust and power over a wide range of operating conditions is discussed. Viscosity effects associated to blade flow separation and stall are taken into account and predicted thrust and power are comparable with results of blade element methods that are largely used in the design of marine current turbines. The accuracy of numerical predictions tend to reduce in cases where turbine blades operate in off-design conditions.

Age-strengthened aluminum alloys as important structural lightweight materials have much lower fatigue properties than conventional steel. This study focuses on the effect of the precipitation-free zone (PFZ) on the mechanical properties of age-strengthened aluminum alloys. A polycrystalline model containing PFZ was created by secondary development of the traditional polycrystalline model and combined with the crystal plasticity finite element method to simulate the stress-strain response and local plastic deformation behaviors under different microstructures. The presence of PFZ significantly reduces the yield strength of the alloy and concentrates the stresses in the PFZ region, especially under cyclic loading, and the strain energy dissipation in the PFZ region is significantly higher than that in the intracrystalline region, and far exceeds that without the precipitation zone. region and far exceeds the model without PFZ, indicating that PFZ is the main region for fatigue crack initiation. In addition, the introduction of a rotation factor to simulate the inhomogeneous rotation within the grain reveals that the additional stress concentration in the PFZ introduced by the aluminum alloy forming process further increases the fatigue crack initiation driving force.

The automotive market is very competitive and presents a demand for the consistent improvement of the used technologies and reduction of the product cost and dimensions with each product model iteration. Hence, it is important to have access to well-defined reference designs of high quality products, in order to evaluate the new ideas and technologies. This paper provides the readers with a numerical model of such a high quality product, i.e., the IPMSM type traction motor from the 4th generation of the Toyota Prius hybrid transaxle. In the introductory section, the brief history of the Prius model is given together with the references to the reverse engineering reports. Afterwards, the machine dimensions, material properties and winding configuration are described. Then, the model is validated with the torque measurements at constant speed. In the next chapters the simulation results are presented and the numerical source data are supplied to the reader. Finally, the design philosophy of the Toyota drive is briefly discussed in comparison with the BMWi3 drive and the results are concluded in the last section.