Introducing the Moving Least-Squares Aided Finite Element Method (MLS-FEM), we present a robust computational technique that effectively evaluates the development of mixing in polymeric mixing devices. Our approach uses the Moving Least-Squares (MLS) interpolant to enhance the classic finite element shape functions. As a result, it becomes possible to treat fixed or moving components, such as rotors of a mixer, as discontinuities within the flow domain. This allows analysis of the flow field with a background mesh independent of the rotors' shape, orientation, and position. We compared the MLS-FEM results with the classical boundary-fitted mesh FEM, which comprises the mesh regeneration steps by the movement of the rotors. We then showed that the new method provides accurate results for the flow field parameters (i.e., the velocity components and the pressure values) and the mixing evolution without requiring mesh regeneration and without compromising accuracy.

This paper introduces an innovative approach to numerically model Structure-Soil-Structure Interaction (SSSI) by integrating the Boundary Element Method (BEM) and the Finite Element Method (FEM) in a coupled manner. To assess the accuracy of the proposed method, a comparative study is undertaken, comparing its outcomes with those generated by the conventional FEM technique. Alongside accuracy, the computational efficiency aspect is crucial for the analysis of large-scale SSSI problems. Hence, the computational performance of the coupled BEM-FEM method undergoes a thorough examination and is compared with that of the standalone FEM method. The results from these comparisons illustrate the superior capabilities of the proposed method in comparison to the FEM method. The novel approach provides more reliable results compared to traditional FEM methods, serving as a valuable tool for engineers and researchers involved in structural analysis and design.

A complete model was developed to simulate the behavior of a circular clamped axisymmetric fluid-coupled Piezoelectric Micromachined Ultrasonic Transducer (PMUT). Combining Finite Difference and Boundary Element Matrix (FD-BEM), this model is based on the discretization of the partial differential equation used to translate the mechanical behavior of a PMUT. In the model, both the axial and the transverse displacements are preserved in the equation of motion and used to properly define the neutral line position. To introduce fluid coupling, a Green’s function dedicated to axisymmetric circular radiating sources is employed. The resolution of the behavioral equations is used to establish the equivalent electroacoustic circuit of a PMUT that preserves the average particular velocity, the mechanical power, and the acoustic power. Particular consideration is given to verifying the validity of certain assumptions that are usually made across various steps of previously reported analytical models. In this framework, the advantages of the membrane discretization performed in the FD-BEM model are highlighted through accurate simulations of the first vibration mode and especially the cutoff frequency that many other models do not predict. This high cutoff frequency corresponds to cases where the spatial average velocity of the plate is null and is of great importance for PMUT design because it defines the upper limit above which the device is considered to be mechanically blocked. These modeling results are compared with electrical and dynamic membrane displacement measurements of AlN-based PMUTs in air and fluid. This complete PMUT model using the FD-BEM approach is shown to be very efficient in terms of computation time and accuracy.

It is well known that the Continuous Galerkin Finite Element (CGFE) method is globally consistent with respect to the first law of thermodynamics. This means that, for any mesh, all obtained discrete solutions will conserve total energy. One might expect, that the method is, also, globally consistent with respect to the second law of thermodynamics. In this paper, we formally study if such conjecture is true. The heat conduction equation is used as the physical model for this analysis. In the present study it is proved that the conjecture is false: at least, for standard piecewise linear (1D and 2D) elements, the CGFE method is not always globally consistent with respect to the second law of thermodynamics. In other words, some obtained discrete solutions can violate the global postulate of the second law, which asserts that total entropy can never decrease.

The use of tuning forks to measure fluid density and viscosity is widely employed in fields such as food, medicine, textiles, automobiles, petrochemicals, and deep drilling. The explicit analytical model based on Euler-Bernoulli cantilever beam theory for the relationship between tuning fork resonance characteristics and the density and viscosity of fluid is only applicable to the situation where the fluid viscous effect is very small. In this paper, the finite element method is used to simulate the influence of large variations in fluid density and viscosity on the resonance characteristic parameters (resonant frequency and quality factor) of the tuning fork. The numerical simulation results are compared with the analytical analysis results and experimental measurement results. Then, the sensitivity of tuning fork resonance characteristic parameters to fluid density and viscosity is studied. The results show that compared with the analytical results, the numerical simulation results have a higher degree of agreement with the experimental measurement results. The relative difference in resonant frequency is less than 2%, and the relative difference in quality factor is less than 4%. This indicates that the finite element method includes the influence of fluid viscosity on tuning fork resonance parameters, which is more in line with the actual conditions than the analytical model. Simulating and analyzing the sensitivity of the tuning fork to fluid density and viscosity by the finite element method, it is possible to consider the situation where fluid density and viscosity vary over a large range. Compared with experimental measurements, this method has higher efficiency and can significantly save time and economic costs. This study can overcome the limitation of existing explicit analytical models, which are only applicable when the viscous effects of the fluid are very small. It enables a more accurate simulation of the coupling vibration between tuning forks and fluids, thereby providing theoretical references for further optimizing tuning fork structural parameters to enhance the accuracy of measuring fluid characteristic parameters.

Arthroplasty is commonly performed to treat advanced osteoarthritis or other degenerative joint conditions, but it can also be considered for young patients with severe joint damage that significantly limits their functionality and quality of life. Young patients are still at risk of aseptic mobilization and bone resorption due to the phenomenon of stress-shielding that causes an uneven distribution of tensions along the femoral contact surface prosthesis. This phenomenon can be limited by appropriately choosing the material of the prosthesis or by varying its stiffness, making sure that its mechanical behavior simulates that of the femur as much as possible. The aim of this study is to evaluate the mechanical behavior of a prosthesis optimized both in shape and material and in comparing the results with a standard titanium prosthesis. Methods: Through three-dimensional modeling and the use of software finite element method (FEM) such as ANSYS we simulate the mechanical behavior of traditional prosthesis and prosthesis optimized topologically respecting the ASTM F2996-13 standard Results: With topological optimization there is a stress reduction from 987 Mpa to 810 Mpa with a mass reduction of 30%. With the use of carbon fiber, it is possible to further reduce stress to 509 Mpa. Conclusions: The reduction of stress on the femoral stem allows an optimal distribution of the load on the cortical bone thus decreasing the problem of stress-shielding.

Understanding and measuring rock characteristics is crucial for scientific and practical applications. This study investigates stone quarries in Morelia, Michoacán, Mexico, evaluating ten quarries (both volcanic and crushed) using diverse analytical methods. Physical properties such as density were measured alongside mechanical tests like uniaxial compressive strength and ultrasonic pulse velocity. Results demonstrate that crushed materials exhibit superior physical and mechanical properties compared to volcanic stones. Advanced techniques including high-magnification scanning electron microscopy and X-ray fluorescence provided insights into surface morphology and elemental composition. Additionally, a novel finite element model was constructed, leveraging dynamic modulus of elasticity (Ed) derived from ultrasonic pulse velocity and maximum load data from uniaxial compressive strength tests. This model facilitated the classification of quarries based on their engineering characteristics. The findings underscore the model’s potential applicability to quarry assessments globally, offering valuable insights into optimizing stone resource utilization and enhancing structural integrity in construction projects.

This paper presents a steady-state and transient heat conduction analysis framework using the polygonal scaled boundary finite element method (PSBFEM) with polygon/quadtree meshes. The PSBFEM is implemented with commercial finite element code Abaqus by the User Element Subroutine (UEL) feature. The detailed implementation of the framework, defining the UEL element, and solving the stiffness/mass matrix by the eigenvalue decomposition are presented. Several benchmark problems from heat conduction are solved to validate the proposed implementation. Results show that the PSBFEM is reliable and accurate for solving heat conduction problems. Not only can the proposed implementation help engineering practitioners analyze the heat conduction problem using polygonal mesh in Abaqus, but it also provides a reference for developing the UEL to solve other problems using the scaled boundary finite element method.

The development of the finite element method (FEM) combined block polynomial interpolation with the concepts of finite difference formats and the variation principle. Because of this combination, FEM overcomes the shortcomings of traditional variation methods while maintaining the benefits of current variation methods and the flexibility of the finite difference method. As a result, FEM is an advancement above the traditional variation methods. The aim of the study is to experimentally determine the thermal behavior of two stomatognathic systems, one control and the other presenting orthodontic treatment by means of a fixed metallic orthodontic appliance, both being subjected to several thermal regimes. In order to carry out this experimental study, we examined the case of a female subject, who was diagnosed with Angle class I malocclusion. The patient underwent a bimaxillary CBCT investigation before initiating the orthodontic treatment. A three-dimensional model with fully closed surfaces was obtained by using the InVesalius and Geomagic programs. Like the tissues examined in the patient, bracket and tube-like components, as well as orthodontic wires, can be included to these models. Once it is finished and geometrically accurate, the model is exported to a FEM-using program, such as Ansys Workbench. The intention was to study the behavior of two stomatognathic systems (with and without a fixed metallic orthodontic appliance) subjected to very hot food (70 ⁰C) and very cold food (-18 ⁰C). From the analysis of the obtained data, it was concluded that, following the simulations carried out in the presence of the fixed metallic orthodontic appliance, significantly higher temperatures were generated in the dental pulp.

Permeability is a fundamental property of porous media, especially when studying resin transfer within a dense fibrous medium to manufacture composite materials. When the measurement of permeability tensor of composite reinforcements is difficult and time consuming, especially in thick, heterogeneous fibrous structures, a numerical approach based on the calculation of the tensor components on a 3D image of the material can be very advantageous. A new finite element based method proposed in this study allows solving very high-dimensional flow problems while limiting the biases associated with boundary conditions and the small size of the numerical samples addressed. The method is validated against academic test cases and against the results of a recent permeability benchmark exercise. The results emphasize that for heterogeneous and anisotropic microstructures, it is important to perform calculations on large samples using boundary conditions that do not suppress the transverse flows that occur when flow is forced out of the principal directions. Since these are not necessarily known in complex media, the permeability determination method must not introduce bias by generating non-physical flows.

The present work deals with the development of finite element methodology for obtaining the stress distributions in thick cylindrical HK40 stainless steel pipe that carry high temperature fluids. The material properties and loading are assumed to be random variables. Thermal stresses that are generated along radial, axial and tangential directions are computed generally using analytical expressions which are very complex. To circumvent such an issue, the probability theory and mathematical statistics have been applied to many engineering problems which allows to determine the safety both quantitatively and objectively based on the concepts of reliability. Monte Carlo simulation methodology is used to study the probabilistic characteristics of thermal stresses which is used for estimating the probabilistic distributions of stresses against the variations arising due to material properties and load. A 2-D Probabilistic finite element code is developed in MATLAB and the deterministic solution is compared with ABAQUS solutions. The values of stresses that are obtained from the variation of elastic modulus are found to be low as compared to the case where the load alone is varying. The probability of failure of the pipe structure is predicted against the variations in internal pressure and thermal gradient. These finite element framework developments are useful for the life estimation of piping structures in high temperature applications and subsequently quantifying the uncertainties in loading and material properties.

This study presents a comprehensive investigation of entropy generation during natural convection in a porous medium with Casson fluid. The governing equations, including the momentum, energy, and entropy balance equations, are solved numerically using finite element method. Through the analysis, various parameters affecting entropy generation are investigated, such as the Casson fluid parameter, Radiation, Prandtl number, and Rayleigh number. The results indicate that the Casson fluid parameter significantly influences the flow and heat transfer characteristics, while the Darcy number and Rayleigh number control the intensity of natural convection. Moreover, the Prandtl number determines the relative significance of heat transfer compared to viscous effects.

The inverse finite element method (iFEM) is a novel method for reconstructing the full-field displacement of structures by discrete measurement strain. In practical engineering applications, the accuracy of iFEM is reduced due to the positional offset of strain sensors during installation, and errors in structural installation. Therefore, coarse and fine two-stage calibration (CFTSC) method is proposed to enhance the accuracy of the reconstruction of structures. Firstly, the coarse calibration is based on a single-objective particle swarm optimization algorithm to optimize the displacement-strain transformation matrix related to the sensor position. Secondly, as selecting different training data can affect the training effect of self-constructed fuzzy networks, this paper proposes to screen the appropriate training data based on residual analysis. Finally, the experiments of the wing integrated antenna structure verify the efficient of the method on the reconstruction accuracy of the structural body displacement field.

There exists a great potential for application in soil mechanics with the use of mechanistic models such as the finite element model, and other non-local models, such as the peridynamic model, since the soil mass can be modelled as a group of particles that interact with each other. In this work, we determine the bearing capacity of deep foundations where stresses are transmitted to deep deposits, which generally present better characteristics in terms of compressibility and shear strength. One of the main elements are driven piles, which must be previously designed, made and tested before their final use. In this work, a comparison with different methods is presented. The finite element method and the peridynamic model are used. Because the effect of consolidation was not taken into account, in these examples, the assumption is made that the structures are cemented in inert soil. Likewise, it is assumed that the water table has no influence on the simulations. Failure envelopes were observed where bonds between particles present high gradients of deformation and fracture. Bearing capacities were estimated and compared with those obtained from the Terzaghi and Meyerhof methods.

Design and multi-objective optimization of auxetic sandwich panels (ASPs) are performed to enhance the blastworthiness of armored fighting vehicles (AFVs). Various metastructures in the form of four auxetic geometries are proposed as the sandwich core: re-entrant honeycomb (REH), double-arrow honeycomb (DAH), star honeycomb (SH), and tetra-chiral honeycomb (CH). This paper employs a combination of finite element (FE) and machine learning (ML) methodologies to evaluate blastworthiness performance. Optimization is carried out using the non-dominated sorting genetic algorithm II (NSGA-II) method. The optimization results show significant improvements in blastworthiness performance, with notable reductions in permanent displacement and enhancements in specific energy absorption (SEA). Global sensitivity analysis (GSA) using SHapley Additive exPlanations (SHAP) reveals that cell thickness is the most critical factor affecting blastworthiness performance, followed by the number of cells and corner angle or radius for CH. The application of optimized ASP on AFVs shows promising results, with no failure occurring in the occupant floor. Furthermore, AFVs equipped with the optimized ASP DAH significantly reduce maximum displacement and acceleration by 39.00% and 43.56%, respectively, and enhance SEA by 48.30% compared to optimized aluminum foam sandwich panels (AFSPs). This study concludes that ASPs have potential applications in broader engineering fields.

In order to examine the impact of geometric configuration of the human vocal tract on the redistribution of acoustic energy in the frequency spectrum, it is necessary to assemble computation models, which facilitate simple modification of the model’s geometric characteristics. This article describes how the human vocal tract model was constructed using so-called super-elements using the improved reduced system method. The mentioned procedure allows using such a model for demanding optimization calculations, for which the model is modified on the super-element level. Furthermore, the article presents an alternative approach, in which an isoparametric element of a higher degree polynomial was derived. Both approaches were compared and used to find the geometric configuration of the human vocal tract amplifying acoustic energy values on the frequency spectrum to around 3kHz.

We present a modified characteristic finite element method that exhibits second-order spatial accuracy for solving convection-reaction-diffusion equations on surfaces. The temporal direction adopted backward-Euler method, while the spatial direction employed surface finite element method. In contrast to regular domains, it is observed that the point in the characteristic direction traverses the surface only once within a brief time. Thus, the good approximation of the solution in the characteristic direction holds significant importance for the numerical scheme. In this regard, Taylor expansion is employed to reconstruct the solution beyond the surface in the characteristic direction. The stability of our scheme is then proved. A comparison is made with an existing characteristic finite element method based on face mesh. Numerical examples are provided to validate the effectiveness of our proposed method.

This study investigates the solution of a boundary value problem using two numerical methods: the Rayleigh-Ritz method and the finite element method (FEM). The aim is to compare the performance of these methods and assess the reliability of FEM as a generalization of the Rayleigh-Ritz approach for more complex problems. The Rayleigh-Ritz method and the linear element formulation of the finite element method were employed to solve the boundary value problem. A detailed comparison of the results obtained from both methods was performed. Graphical illustrations were used to present the solutions, and potential sources of error were analyzed, including element and domain approximation errors, round-off errors, and the impact of using linear rather than quadratic elements in FEM. The solutions generated by both methods were found to be in close agreement, demonstrating that FEM is a viable alternative to the Rayleigh-Ritz method for solving boundary value problems. The minor discrepancies observed can be attributed to approximation errors and the choice of linear elements in the finite element analysis. This work highlights the applicability and effectiveness of both the Rayleigh-Ritz method and FEM in solving boundary value problems. It underscores the finite element method’s flexibility, especially in handling more complex boundary conditions and geometries, and contributes to the understanding of the factors influencing the accuracy of numerical methods in structural analysis.

The variable-length arbitrary Lagrange-Euler (ALE)-ANCF finite element, which employ nonrational interpolating polynomials, cannot exactly describe the rational cubic Bezier curves such as conic and circular curves. The rational absolute nodal coordinate formulation (RANCF) finite element, whose reference length (undeformed length) is constant, can exactly represent the rational cubic Bezier curves. A new variable-length finite element called the ALE-RANCF finite element, which is capable of accurately describe the rational cubic Bezier curves, is proposed by combining the desirable features of the ALE-ANCF and RANCF finite element. In order to control the reference length of ALE-RANCF element within a suitable range, element segmentation and merging schemes are proposed. It is demonstrated that exact geometry and mechanic is maintained after the ALE-RANCF element is divided into two shorter ones, and compared with the ALE-ANCF elements, there are smaller deviations and oscillations after two ALE-RANCF elements are merged into a longer one. Numerical examples are presented and the feasibility and advantages of the ALE-RANCF finite element are demonstrated.

This paper investigates, experimentally and numerically, the buckling behavior of innovative sustainable Low-Cost Bamboo Composite (LCBC) structural columns under compressive loading. The LCBC columns are manufactured from bamboo culms in combination with bio-based resins to form composite structural columns. Different LCBC cross-sectional configurations are investigated in this study including the Russian doll (RD), Big Russian doll (BRD), Hawser (HAW), and Scrimber (SCR). Extra-large, large, medium, and small sizes of bamboo are employed to form the proposed configurations. Two bio-based resins including one bio-epoxy and one furan-based resin, in addition to a soft bio-based filler and a synthetic epoxy resin are applied. The bamboo species used as the cast-in-place giant bamboo for all configurations include Moso, Guadua, and Tali. In addition to experimental testing of the slender LCBC short columns, finite element analysis (FEA) software ABAQUS is applied to model the composite columns and predict their behavior. A buckling analysis using static Riks method is carried out to determine the LBCB’s response to axial compressive loading. A model is calibrated using test results with various imperfections, accounting for the variable behavior of each bamboo specimen, to accurately predict the behavior of LCBC columns with different resin contents. Slender LCBC columns showed maximum stress at buckling up to 60 MPa, highlighting the potential of bio-based resins for structural applications. The study found that the samples with bio-epoxy resin (BE2) exhibited enhanced material stiffness when compared to with synthetic epoxy (EPX) and furan-based resin (PF1), while PF1 specimens demonstrated increased ductility. The BE2 samples showed strain hardening in their stress-strain behavior, and the PF1 samples exhibited a long plateau phase before reaching the maximum load. Among the specimens with Moso bamboo and BE2 resin, those with SCR and HAW configurations achieved the highest compressive strengths.

Ming-style round-backed armchairs as a classic Chinese chair design are renowned for the blend of ergonomics, practicality, and elegance. Although the Chinese armchair furniture industry has made strides in sustainability, there is still a need for further improvements in material selection, manufacturing processes, and supply chain management to fully support the development of sustainable furniture. Laminated bamboo lumber as new sustainable materials are becoming a popular option for furniture designers in the era of "sustainable design". This research aims to investigate the feasibility and application of laminated bamboo lumber in the design of Chinese armchair furniture and proposes innovative ideas for optimizing the armchair's design. First, an empirical study was conducted to comprehensively understand the classification of Ming-style round-backed armchairs’ structure to develop a 3D model. Secondly, an experimental research method was then employed to explore the parameters of laminated bamboo lumber for future studies. Additionally, a validation experiment was conducted to compare real-world scenarios with simulations using finite element analysis in Ansys. Finally, finite element analysis was used to test the feasibility and application of laminated bamboo lumber in furniture design, focusing on the material’s mechanical properties. Results indicated that laminated bamboo lumber has excellent mechanical properties as furniture design materials, subsequently, using the optimized results, the use of 0.6319 kg of steel was successfully reduced by 23.68% to ensure the armchair’s stability. This provides a reference for the future compter-assist and standardized innovative design of Chinese armchair. By incorporating compter-assist design and lightweight optimization, it is possible to save materials and make more efficient use of resources, thereby contributing to sustainable development goals in the field of furniture design.

New families of direct serendipity and direct mixed finite elements on general planar, strictly convex polygons were recently defined by the authors. The finite elements of index r are H1 and H(div) conforming, respectively, and approximate optimally to order r+1 while using the minimal number of degrees of freedom. The shape function space consists of the full set of polynomials defined directly on the element and augmented with a space of supplemental functions. The supplemental functions were constructed as rational functions, which can be difficult to integrate accurately using numerical quadrature rules when the index is high. This can result in a loss of accuracy in certain cases. In this work, we propose alternative ways to construct the supplemental functions on the element as continuous piecewise polynomials. One approach results in supplemental functions that are in Hp for any p≥1. We prove the optimal approximation property for these new finite elements. We also perform numerical tests on them, comparing results for the original supplemental functions and the various alternatives. The new piecewise polynomial supplements can be integrated accurately, and therefore show better robustness with respect to the underlying meshes used.

To enable purposeful design and implementation of automated concrete technologies, precise assessment and prediction of the complex material flow at various stages of the process chain are necessary. This paper investigates the intermediate stage of the extrusion and deposition processes in extrusion-based 3D-concrete-printing, using a numerical model based on the Particle Finite Element Method (PFEM). In PFEM, due to the Lagrangian description of motion, remeshing algorithms and the alpha shape method are used to track the free surface during large deformation scenarios. The Bingham constitutive model was used for describing the rheological behaviour of fresh concrete. This model is validated by comparing the numerically predicted layer geometries with those obtained from laboratory 3D printing experiments. Extensive parametric studies were then conducted using the numerical simulation, delineating the influence of process and material parameters on the layer geometries, the dynamic surface forces generated under the extrusion nozzle and the inter-layer interactions.

There are different ways to analyse energy absorbance (EA) in the human auditory system. In previous research, we developed a complete finite element model (FEM) of the human auditory system. In this mentioned work, the external auditory canal (EAC), middle ear, and inner ear (spiral cochlea, vestibule, and semi-circular canals) were modelled based on human temporal bone histological sections. Multiple acoustic, structure and fluid-coupled analyses were conducted using the FEM to perform harmonic analyses in the 0.1–10 kHz range. Once the FEM had been validated with published experimental data, the FEM numerical results were used to calculate the EA or energy reflected (ER) by the tympanic membrane. This EA was also measured in clinical audiology tests which were used as a diagnostic parameter. A mathematical approach was developed to calculate the EA and ER, with numerical and experimental results showing adequate correlation up to 1 kHz. Another published FEM had adapted its boundary conditions to replicate experimental results. Here, we recalculated those numerical results by applying the natural boundary conditions of human hearing and found that the results almost totally agreed with our FEM. This boundary problem is frequent and problematic in experimental hearing test protocols: the more invasive they are, the more the results are affected. One of the main objectives of using FEMs is to explore how the experimental test conditions influence the results. Further work will still be required to uncover the relationship between the middle ear structure and EA to clarify how to best use FEMs. Moreover, the FEM boundary conditions must be more representative in future work to ensure their adequate interpretation.

This study investigates the magnetization of Cobalt (Co) nanoparticles with two different crystal structures: a centered cubic (fcc) and an hexagonal close packed (hcp), under the magnetic field of the Earth. For that reason, the hysteresis loops of the two material samples are exported and compared with the corresponding value by experimental procedure in the laboratory under the Earth’s magnetization. For both samples, the nanoparticle number variation and the nearest neighbor distance effect, in the Co hysteresis loop are investigated, structural information that is useful for the sample preparation. For the simulations the Finite Element Method (FEM) was used, while experiments are conducted using a measurement set-up based on fluxgate magnetometers. From the simulation, it is proven that, for both samples, the nanoparticle number effect does not change the hysteresis loop. In contrast, when the distance between the nanoparticles increases, the magnetization obtains smaller value. This study aims to build a novel, simple and economical measurement arrangement enabling in situ magnetic measurements for ferromagnetic nano-material characterization. As the need for continuous knowledge of the condition and the phase of ferromagnetic nanoparticles that is crucial in many applications, like catalytic processes, this work aims to introduce an alternative approach that will add another dimension to conventional chemical reaction engineering this of the magnetic state of the material resolved in time and space and to contribute in the field of the monitoring and the early fault diagnosis of nanoparticle materials.

Evidence regarding the effect of the onlay preparation design for different CAD/CAM restorative materials considering the preservation of cusps is lacking. Molars were 3D modeled in four preparation designs for onlay restoration: traditional design with functional cusp coverage (TFC), non-retentive design with functional cusp coverage (NFC), traditional design with non-functional cusp coverage (TNFC), non-retentive design with non-functional cusp coverage (NNFC). Restorations were simulated with two CAD/CAM restorative materials: LD – lithium disilicate (IPS e.max CAD) and RC - resin composite (GrandioBloc). A 100 N axial load was applied to the occlusal surface simulating the centric contact point. Von Mises (VM) and maximum principal (Pmax) stresses were evaluated for restorations, cement layer and dental substrate. The non-retentive preparation design reduced the stress concentration in the tooth structure in comparison to the conventional retentive design. For LD onlays, the stress distribution on the restoration intaglio surface showed that the preparation design as well as the prepared cusp, influenced the stress magnitude. The non-retentive preparation design provided better load distribution in both restorative materials and more advantageous for tooth structure. The resin composite restoration on non-functional cusp is recommended when functional cusp is preserved, in order to associate conservative dentistry and low stress magnitude.

A reduced-dimension (RD) method based on the proper orthogonal decomposition (POD) technology and the linearized Crank-Nicolson mixed finite element (CNMFE) scheme for solving the 2D nonlinear Extended Fisher-Kolmogorov (EFK) equation is proposed. The method reduces CPU runtime and error accumulation by reducing the dimension of the unknown CNMFE solution coefficient vectors. For this purpose, The CNMFE scheme of the EFK equation is established and the uniqueness, stability and convergence of the CNMFE solutions are discussed. Subsequently, the matrix-based RDCNMFE scheme is derived by applying the POD method. Furthermore, the uniqueness, stability and error estimates of the linearized RDCNMFE solution are proved. Finally, numerical experiments are carried out to validate the theoretical findings. In addition, we contrast the RDCNMFE method with the CNMFE method, highlighting the advantages of the dimensionality reduction method.

Brushless Direct Current Motors are significantly improved for using many electrical applications. The design investigation is increased year by year and most of the motor design researchers find out the superior efficiency of the BLDC motor. The traditional motor consumed more electrical power than BLDC motor. The design of BLDC motors that are compact, small, easy to control, save energy and light weight. Due to this fact, the design of BLDC motor should be analyzed to improve the motor characteristics performance. The research focuses on the design of BLDC motor performance using Taguchi Method and Response Surface Methodology (RSM) for Finite Element Method (FEM). Firstly, the Taguchi Method is applied to identify the optimization point of design parameters and consider using the reference model of actual BLDC electric fan motor with a rated power of 26Watts design. Based on the Taguchi Method to create the RSM is used to two-step optimization design performance using in FEM analysis. The performances of motor design proposed approach to provide the effeteness for motor design investigator. In this investigation aims the better structural performances in output power, torque, efficiency and torque ripple. The final design of proposed model achieved an energy efficiency is 10% higher than the reference design of prototype modes, the output power result is 10W higher this reference design. We can reduce the loss of the proposed motor design. Furthermore, the maximum torque can also increase the 0.032Nm of that design. The design simulation results were produced and studied by the FEM analysis of JMAG software. From the evaluation results, the configurations of reference model design are compared with the proposed design by using the three methods of motor design optimization showed the best characteristics for performance of motor investigation design.

In recent years, there has been a growing interest in the preoperative modeling of fluid-structure interaction in the treatment of cerebral aneurysms. In this study, we investigate two cases involving laminar incompressible fluid flow interacting with hyperelastic materials (a vertical flap and aneurysm walls) under the effect of fluid flow. We present a finite element method (FEM) for these prototypical two-dimensional ($2D$) configurations, taking into account the complex flows and deformation of the variant's structure models. We use an Arbitrary Lagrangian-Eulerian (ALE) formulation in a continuum, fully monolithic coupled way, and discretize the fluid and solid domains using the quadratic LBB stable $P_2P_1$ finite element pair to approximate the displacement, velocity, and pressure spaces independently. The resulting discretized form of the nonlinear algebraic system is linearized using a variant of Newton's procedure, and the Jacobian matrices are approximated via the divided difference method. The resulting linear systems are solved using a direct solver MUMPS. We evaluate hydrodynamic forces such as drag and lift coefficients for nonlinear elastic material models, including Saint-Venant Kirchoff, Neo-Hookean, and Mooney-Rivlin (2 parameters) separately. To gain more physical insight into the problem, we verify the computed results by comparing the velocity, viscosity, and pressure fields. Our aim is to qualitatively analyze the changes in the mechanics of the wall behavior of the elasticity of the flap vs. the complexity of deformation, interface, details for prototypical flow situations, corresponding displacement, and deformation. We also provide a physical interpretation of blood flow magnitude and basic hemodynamic phenomena of wall shear stresses (WSS) in aneurysma and the range of deformation. Overall, this study presents a comprehensive approach to modeling fluid-structure interaction in cerebral aneurysms, which could help in the development of more effective treatment strategies.

To reduce weight and simplify maintenance, ship structures frequently include openings and cut-outs. While these features offer practical advantages, they can weaken the structural integrity of key components. This study explores the effects of these geometric discontinuities on the double-bottom floor plates of a Panamax-class container ship under axial and transverse loads. Through numerical simulations and experimental testing, we analyzed different cut-out configurations and stiffening strategies to assess their impact on stress distribution, plate thickness and fatigue behavior. The results reveal that side cut-outs significantly increase stress, particularly under transverse loads, while central openings have less impact. Additionally, increasing plate thickness consistently reduces stress levels across all models, improving structural durability. Fatigue analysis shows that certain stiffening configurations, particularly those with longitudinal stiffeners in the bottom plates, enhance fatigue life. These findings offer critical design insights for optimizing hull structures, balancing weight reduction with strength improvement. The study provides valuable recommendations for reducing stress concentrations and extending the fatigue life of ship components, contributing to more efficient and safer ship designs.

Objective: Despite long simulation times, recently developed finite element analysis (FEA) models of knee joints have demonstrated their suitability for predicting individual risk of onset and progression of knee osteoarthritis. Therefore, the objective of this study was to assess the feasibility of machine learning (ML) to replicate outcomes obtained from FEA when simulating mechanical responses and predicting cartilage degeneration within knee joint. Design: Two ML models based on the Gaussian Process Regression (GPR) algorithms were developed. The first model (GPR1) utilized age, weight, and anatomical joint dimensions as predictor variables to predict tissue mechanical responses and cartilage degeneration based on FEA data. The second model (GPR2) utilized age, weight, height, and gender to predict anatomical joint dimensions, which were then used as inputs in the GPR1 model. Finally, the GPR1 and combined GPR1+GPR2 models were used to investigate the importance of clinical imaging when making personalized predictions for knees from healthy subjects with no history of knee injuries. Results: In the GPR1 model, R2 of 0.9 was exceeded for most of the predicted mechanical parameters. The GPR2 model was able to predict knee shape with R2 of 0.67 - 0.9. Both GPR1 and combined GPR1+GPR2 models offered equally good performance (AUC = 0.73 - 0.74) at classifying patients at high risk for the onset and development of knee osteoarthritis. Conclusions: In the future, real-time and easy-to-use GPR models may provide a rapid technology to evaluate mechanical responses within knee for researchers or clinicians who have no former knowledge from FEA.

(1) Background: This study aimed to compare the effect of polyetheretherketone (PEEK) and titanium abutments on implant prosthetic systems and supporting bone using three-dimensional finite element analysis (FEA); (2) Methods: A three-dimensional finite element analysis was conducted using CATIA and Abaqus software to model mandibular first molar implant systems with titanium and PEEK abutments. Under external loading conditions, finite element analysis was conducted for the stresses in implant components and surrounding bones of each group; (3) Results: The implant fixture of the PEEK model exhibited the highest von Mises stress (VMS). The lowest VMS was observed in the abutment screw of the titanium model. Both implant systems demonstrated similar stress distribution and magnitudes in cortical and cancellous bones; (4) Conclusion: PEEK abutments show similar stress distribution in surrounding bone compared to titanium. However, PEEK absorbs the stresses within the implant system and exhibits the highest VMS values due to its low mechanical and physical properties. Therefore, PEEK abutments need improved mechanical properties for better clinical application.

This paper introduces a low-tech capacitive micromachined ultrasonic transducer (CMUT) designed with low environmental footprint materials. The fabrication process involves a copper plate as the fixed bottom electrode, a polymer-based adhesive as a dielectric material, and an aluminum foil as the top electrode.Finite element simulations include studies of displacement, mechanical stresses, and eigenfrequency. Experimental measurements validate the device's electromechanical behavior, showing an eigenfrequency of 88.6 kHz and a displacement of 22 pm. The low-tech CMUT demonstrates potential for applications such as ultrasonic actuation and energy harvesting, offering simplicity, biocompatibility, and low environmental impact. While not directly ready for applications, these transducers provide hands-on experience with technology similar to high-performance silicon-based implementations. These low-tech MUTs are perfect practical case studies for teaching purposes, combining simulation and experimental validation.

The ANSYS 2024R1 student edition was used to create a finite element model of a basketball rim and backboard. This finite element model included the use of steel for the rim and its mount, a tempered glass backboard, and an aluminum frame behind the backboard. After a mesh was created, fixed support boundary conditions were applied to the four corners of the aluminum frame. The subsequently calculated mode shapes and frequencies were compared to empirical modal analysis previously done at the United States Military Academy at West Point, New York. Four mode shapes and frequencies agreed rather well between the theoretical finite element analysis and previously published empirical modal analysis, specifically where the rim was vibrating in the vertical direction, which was the direction that the accelerometer was aligned for modal analysis. However, three mode shapes missed by the empirical modal analysis were found where the vibration of the rim was confined to the horizontal plane, which was orthogonal to the orientation of our accelerometer.

The design of water Vessel fender system can minimize impact energy, which could result in severe vessel or offshore structure damage during collision. In order to protect Offshore structures and vessels from damage after impact, a Fender design approach is proposed which uses Soft materials like rubber, sand, inflated rubber tube, etc. to reduce damage on impacted vessels as well as offshore structures like wind turbines, bridges, oil drilling platforms in water etc. The fender is designed in such a way that it can handle vessels made of different material, sizes, and weights.The proposed solution is a fender consisting of a number of layers of soft sandwich materials, each layer getting softer as it moves away from the structure i.e. Hard-Soft fender design which can absorb impact energy upon collision. To handle lighter vessels made of soft or brittle materials, the outer layer is made of inflatable or Very soft foam rubber; however, the inner layer covering the offshore structure is made of hard material to handle heavy vessels made of hard material such as steel. Furthermore, the use of sand as a natural fender along waterways such as narrow channels is also proposed to reduce the cost of fender installation and catastrophic vessel impact damage on rocky ground since it produces less deformation just like rubber fender as demonstrated in the simulation. Finite Element Analysis using ANSYS software, Explicit Dynamic Autodyn, was performed. The results showed a considerable decrease in the total deformation and internal energy of the vessel and fended structure when collision occurs. Furthermore, the result showed that Hard-Soft fender layer performs much better than Soft-Hard-Soft followed by Soft-Hard Fender layer design.

In this paper, we aim to address the challenge of airflow interference during fault detection in high-speed train bogies by introducing a flow field and investigating the characteristics of the sound field distribution of critical components under its influence. This approach overcomes the limitation observed in previous studies that ignored the impact of the flow field. Furthermore, we evaluate and compare various layouts for inter-track acoustic sensor arrays. The proposed sensing methods hold significant practical engineering value, and the collected data serves as a valuable reference for acquiring acoustic signals from key bogie components in future research.

An estimation of the electric sources in the heart were conducted using a novel method, based on Huygens’ Principle, aiming at a direct estimation of equivalent bioelectric sources over the heart surface in real time. The main scope of this work was to establish a new, fast approach for the solution of the inverse electrocardiography problem. The study was based on recorded electrocardiograms (ECGs). Based on the Huygens’ principle, measurements obtained from a patient’s thorax surface were interpolated over the surface of the employed volume conductor model and considered as a secondary Huygens’ sources. These sources, being non-zero only over the surface under study, were employed to determine the weighting factors of the eigen-functions expansion describing the generated voltage distribution over the whole volume conductor volume. With the availability of the potential distribution stemming from measurements, the electromagnetics reciprocity theorem is applied once again to yield the equivalent sources over the pericardium. The methodology is self validated, since the surface potentials calculated from these equivalent sources are in very good agreement with ECG measurements. The ultimate aim of this effort is to create a tool providing the equivalent epicardial voltage or current sources in real time, i.e. during the ECG measurements with multiple electrodes.

Model-based approaches form the basis of the majority of damage detection and localization studies; however, built-in online detection technique can offer an alternative for future structural health monitoring technologies. Thus, in this paper, we present a dynamic localization algorithm for damage localization in functionally graded plates. The method depends on the creation of a grid matrix that contains the dynamic response of the structure over time. Subsequently, an optimization process is carried out via a linear equation representing the information contained within the grid to achieve accurate localization of damage. Due to the inherent sparsity of the localization nature, the FISTA algorithm has been employed as a problem solver. The approach is tested in the case of the functionally graded plate under clamped free boundary conditions. Several, damage scenarios were investigated for the case of damage signals on and off the grid. The results show that the proposed approach is capable of accurately predicting damage position and can be suitable for use in low size data systems.

A novel phase-field method based numerical approach to modeling the compaction process of sedi- ments is presented. Water-logged rock or soil sediments are deposited on the water basins over time and the increasing volume of the sediment compacts under its own weight and external pressure. Coupled evolution of mass conservation, Darcy flow, and the viscoelastic constitutive response, in conjunction with the evolution of porosity and permeability, make this problem highly non-linear and involve moving boundaries. We adopt a phase-field approach to represent the moving sediment interface, and the underlying multiphasic model permits for studying compaction in dynamically evolving sediments. This approach does not necessitate a change of domain sizes as is the case of existing traditional models of sedimentation but rather treats sedimentation growth as a problem of interface motion. We first model a classical compaction problem in 1D to compare with existing results in the literature, and then extend the framework to 2D to model the compaction process taking place under the influence of gravity. The model is then extensively applied to understand the effect of the sediment initial state and the sediment material properties on the compaction process and its spatiotemporal evolution. Lastly, a strong validation of our phase-field treatment results against predictions from traditional methods of solving open boundary sediment compaction problems, and a good agreement between the conventional understanding and the numerical predictions of porosity dependence on the fluid pore pressure for various cases of geological interest, are used to demonstrate the applicability and scalability of this novel numerical framework to model more advanced sediment compaction problems and geometries.

Polymer nanocomposites are composed of polymer materials reinforced with nano fillers. In the present study the effective thermal conductivity of the composites filled with nanofillers has been investigated using commercially available finite element software ABAQUS 6.11. The nanofillers used were alumina particles and multi-walled carbon nanotubes (MWNTs) and the matrix was considered to be made of epoxy. For the analysis 2D and 3D Representative Volume Elements (RVEs) were generated using Random Sequential Adsorption (RSA) algorithm using MATLAB and Python scripts. Thermal conductivity was found out for 2D and 3D RVEs for different area and weight fractions respectively. Two different shapes of the alumina nanoparticles were considered: spherical (circular) and ellipsoidal (elliptical) for 3D (2D) analysis. It was found that the thermal conductivity was increased with the addition of nanofillers. The increase in thermal conductivity was approximately same for both types of inclusions at corresponding area or weight fractions in 2D or 3D analysis. The results showed that addition of MWNTs to the composites lead to a significant increase in thermal conductivity than spherical or ellipsoidal inclusions.

This work presents a linear smoothing scheme over high-order triangular elements in the framework of a cell-based strain smoothed finite element method for two-dimensional nonlinear problems. The main idea behind the proposed linear smoothing scheme for strain-smoothed finite element method (S-FEM) is no subdivision of finite element cells to sub-cells while the classical S-FEM needs sub-cells. Since the linear smoothing function is employed, S-FEM is able to use quadratic triangular or quadrilateral elements. The modified smoothed matrix obtained node-wise is evaluated. In the same manner with the computation of the strain-displacement matrix, the smoothed stiffness matrix and deformation graident are obtained over smoothing domains. A series of benchmark tests are investigated to demonstrate validity and stability of the proposed scheme. The validity and accuracy are confirmed by comparing the obtained numerical results with the standard FEM using 2nd-order triangular element and the exact solutions.

This paper presents an approach to the Timoshenko beam theory (TBT) using the finite difference method (FDM). The Timoshenko beam theory covers cases associated with small deflections based on shear deformation and rotary inertia considerations. The FDM is an approximate method for solving problems described with differential equations. It does not involve solving differential equations; equations are formulated with values at selected points of the structure. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. The model developed in this paper consisted of formulating differential equations with finite differences and introducing additional points at the beam’s ends and at positions of discontinuity (concentrated loads or moments, supports, hinges, springs, brutal change of stiffness, spring-mass system, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends. Moreover, grid points with variable spacing were considered, the grid being uniform within beam segments. First-order, second-order, and vibration analyses of structures were conducted with this model. Furthermore, tapered beams were analyzed (element stiffness matrix, second-order analysis, vibration 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.

The article presents the results of experimental studies on wound models made in ballistic gel casts and dressed with a vacuum dressing. Compared to previously published articles, the thesis presents quantitative parameters describing the behavior of wounds under the influence of vacuum. A significant reduction in wound volume was noted. This effect is independent of the size of the wound. The experimental studies results were supported by numerical calculations using the finite element method.

In this paper, we integrate the Finite Element Method (FEM) with orthogonal polynomial approximation in high-dimensional spaces to innovatively model the Moon’s surface gravity anomaly. Our aim is to approximate solutions to Laplace’s classical differential equations of gravity, employing classical Chebyshev polynomials as basis functions. Using classical Chebyshev polynomials as basis functions, the least-squares approximation was used to approximate discrete samples of the approximation function. These test functions provide an understanding of errors in approximation and corresponding errors due to differentiation and integration. These test functions provide an understanding of errors in approximation and corresponding errors due to differentiation and integration. The first application of this project is to substitute the globally-valid classical spherical harmonic series of approximations with locally-valid series of orthogonal polynomial approximations (i.e. using the FEM approach). With an error tolerance set at 10−9ms−2, we use this method to adapt the gravity model radially upwards from the lunar surface. The results showcase a need for a higher degree of approximation on and near the lunar surface, with the necessity decreasing as the radius increases. Notably, this method achieves a computational speedup of five orders of magnitude when applying the method to radial adaptation. More intrinsically, our second application involves using the methodology as an effective tool in solving boundary value problems. Specifically, we implement this approach to solve classical differential equations involved with high precision, long-term orbit propagation. This application provides a four-order of magnitude speedup in computational time while maintaining an error within the 10−10ms−2 error range for various orbit propagation tests. Alongside the advancements in orthogonal approximation theory, the FEM enables revolutionary speedups in orbit propagation without compromising accuracy.

Boeing and Airbus developed a special testing procedure to investigate the compressive response of laminates that have been impacted (following standards ASTM D 7137 and DIN 65561). This study focuses on both experimental and numerical analysis of Kevlar plates subjected to compression after impact. To ensure high quality and appropriate mechanical properties, the composite plates were manufactured using autoclaving. The DIN 65561 protocol was followed for all three test systems. Initially, ultrasonic C-scanning was performed on all plates before testing to confirm they were free of any significant defects arising from the manufacturing process. Subsequently, low-energy impact testing was conducted at levels ranging from 0 to 8 Joules. Three specimens were tested at each energy level. After the impact, all specimens underwent ultrasonic C-scanning again to assess the internal delamination damage caused by the impactor. Finally, both pristine and impacted specimens were subjected to compressive testing using the special jig specified in DIN 65561. The compressive impact strength results obtained from these tests were plotted against the delamination area measured by C-scanning. This data was then compared to results obtained from specimens with artificial damage. Semi-empirical equations were used to fit both sets of curves. The same procedure (impact testing, C-scanning, and data analysis) was repeated to investigate the relationship between impact energy and total delamination area. Lastly, finite element modeling was employed to predict the buckling stresses that develop under compression in the impacted systems studied. These modelling approaches have demonstrated good accuracy in reproducing experimental results for CAI tests.

Graphene foam prepared by chemical vapor deposition method is a promising thermal interfacial material. However, the thermal properties of graphene foam highly depend on the experimental fabrication conditions during the chemical vapor deposition process. Aiming to reveal how to prepare the appropriate graphene foam for the various thermal management scenarios, the influence of experimental conditions on thermal properties of graphene foam is investigated. Furthermore, the contribution of thermal conductivity and thermal radiation to the effective thermal coefficient of graphene foam is carried out for comparison. Research results show that, the porosity and the cross-section shape of struts of growth template are two critical factors to affect the thermal transport of graphene foam, especially with the increase of temperature. In addition, the deposition time of graphene determines the wall thickness and affects the thermal conductivity directly. The thermal radiation contributes more than thermal conductivity when the temperature climbs continuously. Comparatively, the effective thermal coefficient of graphene foam composite with high porosity and circular shape struts is much superior to others at high temperature. The research findings can provide important guidance for the graphene foam fabrication and its applications in the field of thermal management.

This study investigates the impact of specimen size on the mechanical properties of wood, specifically focusing on compression strength, elastic modulus, and Poisson's ratio. Compression tests were conducted using three different specimen sizes (20mm×20mm×30mm, 40mm×40mm×60mm, 60mm×90mm×90mm) in the longitudinal, radial, and tangential directions. Mechanical parameters, failure mechanisms, load-displacement curves, and stress-strain relationships were systematically analyzed for each size. The study also evaluated the influence of specimen size on the accuracy of finite element numerical simulations by utilizing the obtained mechanical parameters. The results reveal a significant correlation between compressive strength and specimen size, indicating a decrease in compressive strength with increasing specimen size. Conversely, elastic modulus and Poisson's ratio exhibit less sensitivity to variations in specimen size. Notably, parameters derived from small-sized specimens (20mm×20mm×30mm) exhibited substantial errors, while those obtained from medium-sized (40mm×40mm×60mm) and large-sized specimens (60mm×90mm×90mm) demonstrated greater reliability, providing precise results in finite element numerical simulations.

Lightweight component design is effectively achievable through the sandwich structures, many past research studies in the aerospace and racing sectors (since 1920s) has proven it. To extend their application into automotive and other transport industries, the manufacturing cycle times must be reduced. This can be achieved by thermoplastic sandwich materials of continuous fibre-reinforced thermoplastic (CFRTP) cover layers and thermoplastic honeycomb cores. To widespread the application of flat thermoplastic-based sandwich panels into complex parts, a novel forming technology was developed by the Fraunhofer Institute of Microstructure of Materials and Systems (IMWS). For an increase in reproducibility of the sandwich component production and a decrease in manufacturing defects, numerical modelling methods and experimental validation are needed. The finite element (FE) modelling of this thermoforming process is a thermo-mechanical assignment with all three: material, geometry, and boundary non-linearities. In the scope of this work, the material behavior of different sandwich components (cover layers, honeycomb core, interface) under thermo-mechanical loads is characterized in detail, and a complex numerical model is developed to describe the thermoforming process of thermoplastic-based sandwich laminate. The critical forming condition like wrinkling that occurs due to friction in a sandwich or poor intra-ply shear in cross-ply studied, validated, and optimized. A thermoplastic sandwich made of polypropylene (PP) honeycomb core and polypropylene glass fibre (PP/GF) cross-ply as cover layers was used for this study, and numerical modelling was executed in LS-DYNA software.

The insufficient toughness of concrete can lead to the occurrence of crack expansion, ultimately resulting in a decline in the physical and mechanical properties such as bearing capacity and deformation characteristics of structural components. In severe cases, this may lead to brittle failure, where the entire structure and system fail suddenly without significant deformation after being stressed. Additionally, during the process of shear and bending failure of reinforced concrete beams, the cracks generated on the side of the beam are mostly I-II composite cracks. Considering the material nonlinearity and geometric irregularity of reinforced concrete, it is of great significance to study the dynamic evolution law of crack expansion in reinforced concrete members. Based on experimental data, this paper aims to utilize the nonlinear finite element simulation software Abaqus to conduct finite element simulations of crack expansion in simple supported reinforced concrete beams. By accurately controlling the constitutive curves of concrete and reinforcement materials, the design of seismic performance of structural members can be enhanced, thereby making essential contributions to the preliminary analysis of actual crack control engineering. The nonlinear finite element results indicate that the joint cracks of the frame with the floor slab primarily manifest as flexural cracks, with the concrete damage degree of the floor slab being higher than that of the beam. This finding underscores the importance of considering the behavior of concrete and reinforcement materials in designing for seismic performance and crack control in actual engineering applications.

In this study, the autofrettage method, which is used to increase the durability of cylindrical parts, was applied to a barrel made of steel material according to the AISI 4340 standard. Assuming that the barrel works continuously under the pressure of 400 MPa, optimum autofrettage was performed. The stresses occurring on the barrel before and after autofretage were calculated. Additionally, analytical equivalent stress values (Von-Mises) formed on the part's wall were compared with numerical analysis results obtained with a commercial finite element program ANSYS®, and their compatibility was examined. As a result of the calculations and analysis performed, it was observed that autofretage has positive effect on increasing the pressure carrying capacity of the barrel. After autofrettage the cylindrical part, it was determined that the Von-Mises stress values formed in the part wall changed and these stresses took the lowest value at a point (elastic-plastic interface). It has been evaluated that this point gives the optimal autofretage pressure value for a certain operating pressure. After the autofrettage application, it was determined that the Von-Mises Stress value in the barrel wall decreased by 25.27 % and the pressure carrying capacity of the barrel increased without any design or dimensional change. It has been seen that the analytical and numerical analysis results are consistent with each other, and it is expected that this study will guide and contribute to the studies to be conducted on increasing the strength or lightening the weight of the parts.

In recent years, with the increase in maritime trade, the necessity of increasing the capacities of the ports has emerged. However, while it is planned to increase the capacities of the ports, it is important that the port continues to operate at the same time. In this respect, the old port structures should not be damaged during the capacity increase. In this study, the strengthening of a port in Guinea is discussed as a case study. In the study, the existing quay wall was evaluated, and geotechnical and structural alternatives of the new structure to be built for capacity increase were evaluated. A combined system was designed as a pile foundation and a reinforced foundation with plastic piles so as not to damage the existing quay wall. The pile capacities obtained as a result of the analyses were verified by loading tests.

This paper presents a finite element analysis of the bi-directional orthogonal model, which incorporates individual crack strain separation and tracking. The objective of the research is to expand the current shear friction model to manage bi-directional cracking at any angle, allowing for a more universal model that can be applied to intricate structures and non-proportional loading cases. The proposed model was initially developed as a total strain-based model, with the assumption that crack strains are equivalent to total strains, but was subsequently recalculated to improve accuracy by separating crack strains from total strains. Furthermore, a separate crack strain formulation was created to account for strains in the concrete's uncracked portions and locked-in crack strains. The article then discusses the testing of various convergence methods and loading programs to achieve high convergence. Comparative analyses of the generalized shear friction model with other models for crack orientation and loading cases similar to those of a reinforced concrete membrane are also presented. The MATLAB program successfully applied the bi-directional cracking model for one finite element under a uniform cyclical strain state, using a secant stiffness formulation.

Crystallization often occurs in the laser welding of amorphous alloys, reducing the amorphous alloys’ properties. So, it is necessary to simulate the temperature field curve of the welding area by computer and then determine the power and laser moving speed of laser welding. In this paper, the temperature field curve of Zr41.2Ti13.8Cu12.5Ni10Be22.5(Vit1) amorphous alloy in laser welding is obtained by finite element analysis. The continuous heating curve (CHT) of Vit1 is fitted by the Vogel-Fulcher-Tammann (VFT) equation and the Kissinger equation. If the temperature field curve intersects with the CHT curve, crystallization occurs. The experiment results show that the VFT equation can be used to predict the crystallization of Vit1 better in laser welding. The tem-perature and welding time are increased by using low welding speed. Therefore, the temperature of the weld zone cannot fall in time, resulting intersection of the temperature field curve and the CHT curve. So, crystallization can be avoided if the welding speed is controlled within a rea-sonable range, and the highest temperature is kept under the CHT curve.

Wearable thermoelectric generators (w-TEGs) convert thermal energy into electrical energy to realize self-powering of intelligent electronic devices, thus reducing the burden of battery replacement and charging, and improving the usage time and efficiency of electronic devices. Through finite element simulation, this study successfully designed high-performance thermoelectric generator and made it into wearable thermoelectric module by adopting “rigid device - flexible connection” method. It was found that higher convective heat transfer coefficient on cold-end leads to larger effective temperature difference and better power generation performance of device in typical wearable scenario. Meanwhile, at same convective heat transfer coefficient on the cold-end, longer TE leg length leads to larger temperature difference established at both ends of device, larger device output power and open-circuit voltage. However, when the convective heat transfer coefficient increases to a certain level, optimization effect of increasing TE leg length on device power generation performance will gradually diminish. For devices with fixed temperature difference between two ends, longer TE leg length leads to higher resistance of TEG, resulting in lower device output power but slight increase in open-circuit voltage. Finally, sixteen 16×4×2 mm2 TEGs (L=1.38 mm, W=0.6 mm) and two modules were fabricated and tested. At hot end temperature Th=33 ℃ and cold end temperature Tc=30 ℃, the actual maximum output power Pout of TEG is about 0.2 mW, and the actual maximum output power Pout of TEG module is about 1.602 mW, which is highly consistent with the simulated value. This work brings great convenience to research and development of wearable thermoelectric modules and provides new, environmentally friendly and efficient power solution for wearable devices.

The link of fracture analysis of large ships is an essential part of the structural safety of ships. However, the finite element analysis of large structure according to the traditional method not only has high requirements on the computer, but also the influence of the mesh encryption at special locations makes further requirements on the accuracy of the ship. In this paper, based on the extended finite element method (XFEM) and the multi-point constraint method (MPC), the dynamic crack expansion process of the rectangular plate containing cracks under the multi-scale model is comparatively studied by using ABAQUS simulation software Finally, the applicability of the method is verified by applying the multi-scale model of ship nodes. The results show that the computational effort of the multi-scale model is significantly reduced with the same accuracy and allowable error, which proves the accuracy and efficiency of the method.

This paper presents the experimental and numerical studies of the vertical stress distribution beneath the centerline of circular footings resting on the sand. Vertical stress values due to circular load in medium-dense sand were measured using soil stress transducers. Reduced-scale laboratory model tests were carried out in a square-shaped test box and vertical stresses were measured at the centerline of the circular footing at different depths of the sand. In order to determine the effect of footing width on vertical stress distribution, laboratory experiments were performed on circular footings of different diameters. In addition, a series of finite element analyses was carried out on the test model to validate the findings of the experimental study. In the numerical analysis, the sand was modeled as a non-linear elastoplastic material. The analyses were performed using the finite element method with two-dimensional axisymmetric and three-dimensional conditions. The results obtained from laboratory model tests and numerical analysis were compared with the vertical stress values obtained from the theoretical solution in the literature. Vertical stress distribution in horizontal planes at any depth from the footing base shows the same characteristic form in all methods. However, some differences were observed between experimental, numerical, and theoretical results and discussed.

This article presents a case study on the structural assessment of a reinforced concrete (RC) foundation exposed to low temperatures. The foundation supports a 19,500 m³-capacity tank with low-temperature (-89°C) ethane. Icing and bubbling were observed on the tank’s surface soon after it started operations. Condensation was also observed at the bottom of the 0.8 m-depth RC slab, which raised concerns about the structural condition of the concrete. This study provides details of the field and analytical investigations conducted to assess the structural condition of the foundation. Heat transfer finite element (FE) analyses were performed to examine the concrete sections subjected to low temperatures. It was found that the ethane leakage produced a low temperature on the top side of the concrete foundation of +9.7°C. Overall, the temperatures calculated by the FE analyses were in good agreement with actual field measurements, within a ±5% accuracy. The simplified heat transfer equation for porous media used in this study was sufficiently accurate to model the effects of the ethane leakage in the concrete foundation, provided the ambient temperature at the site is taken into account in the analysis. The results also confirm that reinforcing bars can be neglected in the thermal analysis of massive concrete slabs. The results from the field measurements and FE analyses confirmed that the structural integrity of the RC foundation was never compromised. The approaches, methods and techniques discussed in this article are deemed suitable to solve the practical and scientific challenges involved in the structural assessment and repairs of large special structures. Accordingly, they can serve as useful reference and guidance for engineers and practitioners working in the field of forensic engineering.

The Gol'denveizer problem of a torus was studied analytically by Audoly and Pomeau (2002), and the accuracy of the Audoly and Pomeau linear law was confirmed numerically by Sun (2021). However, the law does not include the major radius R of the torus. To find the influence of the major radius, we used finite element numerical simulation to simulate different cases, and we propose a modified Audoly and Pomeau linear law for vertical deformation, which includes R. A linear law of horizontal deformation is presented as well. Our studies show that the Audoly and Pomeau linear law has high accuracy. With modified vertical and horizontal deformation, a displacement-compatible relation between them is formulated.

Orofacial injuries are common occurrences during contact sports activities; however, there is an absence of data regarding the performance of hybrid occlusal splint mouthguards, especially during compressive loading. To evaluate the total deformation and stress concentration, a skull model was selected and duplicated to receive two different designs of mouthguard devices: one model received a conventional custom-made mouthguard (MG) with 4-mm thickness and the other received a novel hybrid occlusal splint-mouthguard (HMG) with the same thickness. Both models were subdivided into finite elements. The frictionless contacts were used, and a nonlinear analysis was performed simulating the compressive loading in occlusion. The results were presented in von-Mises stress maps (MPa) and Total Deformation (mm). A higher stress concentration in teeth was observed for the model with the conventional MG, while the HMG design displayed a promising mechanical response with lower stress magnitude. The HMG de-sign displayed a higher magnitude of stress on its occlusal portion than the MG design. The hybrid mouthguard (HMG) reduced (1) jaw displacement during chewing and (2) the generated stresses in maxil-lary and mandibular teeth.

Energy efficiency and leakage magnetic field (LMF) are two important issues in electric vehicle inductive chargers. In this work, the maximum achievable coil efficiency and the corresponding LMF strength are formulated as functions of hardware parameters, and figure of merits (FOM) are proposed for assessing the efficiency and LMF performance of the coil assembly pair. The impacts of the coil assemblies’ geometric parameters on both FOMs are examined with the aid of finite element analysis (FEA), and measures to improve the FOMs are extracted from FEA results. A coil assembly pair is manually optimized within given dimensional limits. Compared with the initial design, the optimized one achieves higher efficiency and lower LMF strength while consuming less copper. The performance improvement is verified by FEA results and experimental data measured on an 85 kHz electric vehicle inductive charger prototype. The key measures for coil assembly optimization are summarized.

The effect of the different dental implants positioning region on the stress performance of the implant-supported prosthesis is not yet clear. This study evaluated the dental treatment with six dental implants in three different models and three different occlusal loading conditions, in terms of the biomechanical response of implants, prosthetic screw and maxilla, using three-dimensional finite element analysis. The finite element models were modelled containing external hexagon implants, as well as a Cobalt-Chromium superstructure. Three types of loads were applied: in the area of ​​the central incisors, first premolar and in the second molars. For the finite element simulations, the von-Mises stress peaks in the implant and in the surrounding cortical bone were analyzed. All recorded results reported higher values ​​for the implant-supported prosthesis in group C compared to the groups A and B. The highest stress values, ​​regardless the evaluated model, was in the prosthesis in group C and in screws, the smallest were in group A.

In the process of metal cutting, the anti-wear performance of the tool determines the life of the tool and affects the surface quality of the workpiece. The finite element simulation method can directly show the tool wear state and morphology, but due to the limitations of the simulation time and complex boundary conditions, it has not been commonly used in tool life prediction. Based on this, a tool wear model was established on the platform of a finite element simulation software for the cutting process of titanium alloy TC4 by end milling. The key technique is to embed different types of tool wear models into the finite element model in combination with the consequent development technology. The effectiveness of the tool wear model is validated by comparing the experimental results with the simulation results. At the same time, in order to quickly predict the tool life, an empirical prediction formula of tool wear was established, and the change course of tool wear under time change was obtained.

Since large power transformers are custom-made, and their design process is a labor-intensive task, their design process is split into different parts. In tendering, the price calculation is based on the preliminary design of the transformer. Due to the complexity of this task, it belongs to the most general branch of discrete, non-linear mathematical optimization problems. Most of the published algorithms are using a copper filling factor based winding model to calculate the main dimensions of the transformer during this first, preliminary design step. Therefore, these cost optimization methods are not considering the detailed winding layout and the conductor dimensions. However, the knowledge of the exact conductor dimensions is essential to calculate the thermal behaviour of the windings and make a more accurate stray loss calculation. The paper presents a novel, evolutionary algorithm-based transformer optimization method which can determine the optimal conductor shape for the windings during this examined preliminary design stage. The accuracy of the presented FEM method was tested on an existing transformer design. Then the results of the proposed optimization method have been compared with a validated transformer design optimization algorithm.

This study presents the development of a number of finite element (FE) models of the pelvis using different continuum and structural modelling approaches. Four FE models were developed using different modelling approaches: continuum isotropic, continuum orthotropic, hybrid isotropic and hybrid orthotropic. The models were subjected to an iterative adaptation process based on the Mechanostat principle. Each model was adapted to a number of common daily living activities (walking, stair ascent, stair descent, sit-to-stand and stand-to-sit) by applying onto it joint and muscle loads derived using a musculoskeletal modelling framework. The resulting models, along with a structural model previously developed by the authors, were compared visually in terms of bone architecture, and their response to a single load case was compared to a continuum FE model derived from CT imaging data. The main findings of this study were that the continuum orthotropic model was the closest to the CT derived model in terms of load response albeit having less total bone volume, suggesting that the role of material directionality in influencing the maximum orthotropic Young's modulus should be included in continuum bone adaptation models. In addition, the hybrid models, where trabecular and cortical bone were distinguished, had similar outcomes, suggesting that the approach to modelling trabecular bone is less influential when the cortex is modelled separately.

Concrete-filled steel columns have been extensively used in the world due to having all suitable characteristics of concrete and steel, more ductility, increasing concrete confinement using steel wall, large energy-absorption capacity and appropriate fire behavior. In present paper, concrete-filled steel square columns have been simulated under the influence of blast load using ABAQUS software. These responses will be compared for scaled distances based on the distance to source and weight of explosive material. As result, it can be seen that although concrete deformation has been restricted using steel tube, but inner layer of concrete has been seriously damaged and column displacement will be decreased by increasing scaled distance. We also concluded that concrete-filled steel columns have high ductility and blast resistance.

Rapid post-earthquake damage diagnosis of bridges can guide decision-making for emergency response management and recovery. This can be facilitated using digital technologies to remove the barriers of manual post-event inspections. Prior mechanics-based Finite Element (FE) models can be used for post-event response simulation using the measured ground motions at nearby stations; however, the damage assessment outcomes would suffer from uncertainties in structural and soil material properties, input excitations, etc. For instrumented bridges, these uncertainties can be reduced by integrating sensory data with prior models through a model updating approach. This study presents a sequential Bayesian model updating technique, through which a linear/nonlinear FE model, including soil-structure interaction effects, and the foundation input motions are jointly identified from measured acceleration responses. The efficacy of the presented model updating technique is first examined through a numerical verification study. Then, seismic data recorded from the San Rogue Canyon Bridge in California are used for a real-world case study. Comparison between the free-field and the foundation input motions reveals valuable information regarding the soil-structure interaction effects at the bridge site. Moreover, the reasonable agreement between the recorded and estimated bridge responses shows the potentials of the presented model updating technique for real-world applications. The updated FE model is considered as the digital twin of the bridge and can be used to analyze the bridge and monitor the structural response at element, section, and fiber levels to diagnose the location and severity of any potential damage mechanism.

H-adaptivity is an effective tool to introduce local mesh refinement in FEM-based numerical simulation of crack propagation. The implementation of h-adaptivity could benefit the numerical simulation of fatigue or accidental load scenarios involving large structures such as ship hulls. In engineering applications, the element deletion method is frequently used to represent cracks. However, the element deletion method has some drawbacks such as strong mesh dependency and loss of mass or energy. In order to mitigate this problem, the element splitting method could be applied. In this study, a numerical method called ‘h-adaptive element splitting’ (h-AES) is introduced. The h-AES method is applied in FEM programs by combining h-adaptivity with the element splitting method. Two examples using the h-AES method to simulate cracks in large structures under linear-elastic fracture mechanics scenario are presented. The numerical results are verified against analytical solutions. Based on the examples, the h-AES method is proven to be able to introduce mesh refinement in large-scale numerical models that consist of structured coarse meshes. By employing the mesh refinement introduced in this paper, very small cracks are well represented in large structures.

The problems of modeling the human bladder and its stress-strain state under external static influence were considered. A method for identification of the anisotropic biomechanical characteristics of bladder tissue was proposed. A FEM model was created, which accepts into regard that the bladder is surrounded by fiber and affected by surrounding organs, and partially protected by pelvic bones. The model considered the presence of constant hydrostatic pressure on the walls of the bladder when it is full. It has been shown that isotropic mechanical characteristics of biological tissue can be used for studying the deformed state of a filled bladder if the filled bladder 300 ml is considered as initial non-deformed stage. This was shown by modeling and verification the effect of the external static force on the bladder. Numerical experiments were conducted based on the constructed model. To validate the results obtained, a series of natural experiments on the effect of external pressure on the bladder under ultrasound control was conducted. In the future, it is planned to use the constructed model to study rupture deformations of the bladder under the influence of static and dynamic loads.

The influence of applied loads on free vibrations analysis of plates and shells panels made of functionally graded materials is analyzed in the present work. Formulations for the static analysis considering geometrically nonlinear behavior, as well as linear buckling and free vibrations analyses are considered. The calculation of the through-thickness stress distribution is also performed. A finite element model based on a higher order shear deformation theory and using a non-conforming triangular flat plate/shell element with 3 nodes, and 8 degrees of freedom per node is used in the numerical implementation. The results obtained with this numerical model are presented, discussed, and compared with alternative solutions published by other authors in some benchmark applications.

Cycloid drives are generally used in precision machinery requiring high reduction ratios, such as robot joint (RV) reducers, in which the contact load and stress property of the internal meshing pairs greatly affect the dynamic transmission performance of the reducers. Therefore, in this study, Isogeometric Analysis (IGA) was employed to explore the multi-tooth contact problem of the cycloid pinwheel drive. Based on the NURBS curved surface generation method, the NURBS tooth profile of the cycloid gear was reconstructed. In addition, the NURBS surface of the cycloid gear - pin tooth - output pin was generated via the element splicing method. A geometrical analysis model of cycloid pinwheel drive was established to solve the contact force of the meshing pair under different input angles, and compared with the finite element method in terms of convergence, resultant accuracy, and solving timeliness. The results show that the Isogeometric analysis has higher accuracy and efficiency than the finite element method in calculating the contact stress and contact force. This method can improve the contact simulation accuracy of the cycloid drive and provides a reference for the design evaluation of RV reducer.

As a new type of gear, bevoloid gears have the advantages of compensating axial error, smooth transmission, and eliminating turning error, and they are widely used in applications requiring high transmission accuracy and stability. However, research on calculating the time-varying mesh stiffness of bevoloid gears is still limited, and there is an urgent need to propose a method that can calculate the stiffness of bevoloid gears quickly and accurately. This paper first establishes the bevoloid gear tooth profile expressions, then assumes a pair of bevoloid gears meshing with the same rack and derives the contact line equations of parallel axis bevoloid gear pairs, and analyze the contact process of bevoloid gears. We propose an analytical algorithm that uses the slicing method to calculate the stiffness of helical gears, straight bevoloid gears, and helical bevoloid gears, change the parameters of helical bevoloid gears respectively, and analyze the influence of different parameters on stiffness. Finally, the finite element method is used to verify the analytical method, and the correctness of the analytical calculation results is verified, and the errors are analyzed.

One of the most important problems of geotechnical engineering is the design of pile foundations. Piles are generally designed for vertical and horizontal static loads. However, it is important to design piles against dynamic loads in regions with high seismicity. This study verified a centrifuge experiment obtained from the literature with a finite element program. With the kinematic interaction analyses made for this purpose, the effect of parameters such as layered soil, groundwater, pile stiffness, and earthquake acceleration were investigated numerically. According to the results obtained, it was understood that the earthquake magnitude is the most important parameter in kinematic interaction analyses. In addition, in layered soil conditions, the fact that 75% of the pile length is in the weak soil creates the most unfavorable situation. It has been observed that pile stiffness and groundwater level also have certain effects on the kinematic interaction.

The flavor, aroma and color of coffee may change due to mechanical damage, reducing its qual-ity. To measure the mechanical response of materials, compression tests can be performed, de-termining the effect of force application on the mechanical behavior of the fruit at different stages of ripeness. In this context, this study aimed to analyze the deformation, strain energy and equivalent von Mises stress of coffee fruits at mature, semi-mature and immature maturation stages and compressed by collapse forces. Compression in three directions (x, y and z) was sim-ulated on coffee fruit models using the finite element method, with a compression support in the direction opposite to the force application axis. The numerical simulation of the compression process allowed to verify that the more mature the fruit, the greater the associated mean defor-mation (2.20 mm mm-1, 0.78 mm mm-1 and 0.88 mm mm-1), the lower the mean strain energy (0.07 mJ, 0.21 mJ and 0. 34 mJ) and the lower the mean equivalent von Mises stress (0.25 MPa, 1.03 MPa and 1.25 MPa), corresponding to ripe, semi-ripe and immature fruits; and these proper-ties may be used to analyze how processes in the production chain affect the amount of mechan-ical damage.

Emerging soft actuators with various soft materials and traveling-wave motion have been frequently discussed. Various configurations and the resulting performances have been proposed, and it is very challenging for soft actuators in realizing large strokes. This study presents an experimentally validated nonlinear finite element model to predict the deformation produced by a spiral pneumatic rubber actuator to generate traveling-wave motion. The actuator consists of a membrane mounted on a rubber substrate with three air chambers in a spiral configuration, simplifying the actuator's structure. The sequential deformations of the successive chambers interact with each other and produce radial traveling-waves on the membrane surface, driving the objects placed on the actuator. Finite element analysis using ANSYS computer software analyzed the elastic movement by considering the influence of different initial structural types. The structure with three different length ratios was introduced using the four parameters: membrane thickness, chamber width, chamber depth, and wall thickness. The simulation results indicated an optimal structure with specific ratios. A reasonable correlation was obtained during experimental validation; the predicted displacement values were approximately 17% smaller than the experimental values. Finally, the transportation performance of the prototype was tested and a velocity of 2.28 mm/s in the desired direction was achieved. We expect that our demonstration will expand the range of applications of the spiral pneumatic rubber actuator to include conveying or worm-like robots.

Advanced manufacturing techniques, aimed at implants with high dependability, flexibility, and low manufacturing costs, are crucial in meeting the growing demand for high-quality products like biomedical implants. Incremental sheet forming is a promising flexible manufacturing approach for rapidly prototyping sheet metal components using low-cost tools. Titanium and its alloys are used to shape most biomedical implants because of its superior mechanical qualities, biocompatibility, low weight, and great structural strength. The poor formability of titanium sheets at room temperature, however, limits their widespread use. The goal of this research is to show that gradual sheet formation of a titanium biomedical implant is possible. The possibility of creative and cost-effective concepts for the manufacturing of such complicated shapes with significant wall angles is explored in this study. A numerical simulation based on finite element modeling as well as a design process tailored to metal forming is used to complete the development. The mean of uniaxial tensile tests with a constant strain rate was used to study the flow behavior of the studied material. To forecast the crack, the obtained flow behavior was modeled using the behavior model and failure model.

A connection between the phenomenon of wafer bifurcation and that of spontaneous symmetry breaking (SSB) has been established. Indeed, by developing an analytical approximation of the elastic energy of a wafer coated with a thin layer (e.g. a metal), it is shown as the elastic potential energy, interpreted as a quantity contributing to the thermodynamic free energy, can be investigated within the framework of the Landau theory of the second order phase transitions. The elastic energy of a bifurcated wafer is a complex function of the stress and the curvatures in the two perpendicular directions. In this work, it is shown as a translation of the coordinate system of the curvatures allows to gain a potential which has a “Mexican hat” shape. This is a distinctive trait of the phenomenon of spontaneous symmetry breaking (SSB). Moreover, it is shown as the values of the coordinates at the minimum of the SSB potential agrees with those provided by the theory. Bifurcation is hence a phenomenon that can also be interpreted as a spontaneous symmetry breaking where the rotational symmetry of a disc shaped wafer is broken. It occurs because of a lowering of the wafer energy in the SSB broken symmetry configuration.

With the increasingly stringent airworthiness standards, the noise generated during the rotorcraft flight is gradually attracting people’s attention. It widely operated helicopters at low altitudes because of their maneuverability. The way to reduce the noise caused by the complex airflow of the helicopter rotor system has progressively become a hot topic for researchers. Using a hybrid acoustic analysis method, this paper investigates the improvement of the noise and thrust of the helicopter’s tail rotor through the tail rotor structural parameters. For the basic model, the turbulence simulation is performed using an incompressible detached eddy simulation (DES) method, and the Lighthill acoustic analog equation is calculated using the finite element method (FEM). We verified the accuracy of the method through wind tunnel tests. We chose a series of structural parameters for sound simulation and fluid simulation calculations. The results indicate that the modified tail rotor noise reduced by 16.5 dBA and the total thrust increased by 19.9% from the prototype model. This work can enhance the duct tail rotor design to improve aerodynamic and aeroacoustic performance.

A micromechanical model is implemented to indicate the progressive fatigue problem of a laminated composite with a central circular hole under fatigue loading based on a finite element model. The mechanical properties of the composite lamina are represented based on the characteristics of the fiber and the matrix through a micromechanics model. An appropriate algorithm is then adopted to simulate fatigue damage development in the composite lamina. According to this algorithm, the stress field of the composite subjected to fatigue load is initially obtained using the finite element method. Finally, the predicted results of the stresses in the constituents i.e. fiber and matrix are determined according to the micromechanical bridging model. Finally applying proper damage driving relations leads to damage degree in each element. The proposed model is proven to be successful in the observation of the fatigue behavior with stiffness degradation in each element of the composite in each cycle. Results are reported and validated using those micromechanical models and experimental data available in the literature.

In this paper, the results of finite element analyses of a single-layer cylindrical latticed shell under severe earthquake is presented. A 3D Finite Element model using fiber beam elements were used to investigate the collapse mechanism of this type of shell. The failure criteria of structural members are simulated based on the theory of damage accumulation. Severe earthquake of peak ground acceleration (PGA) of 500 gal was applied to the shell. The stress and defeomation of the shell were studied in detail. A three-stages collapse mechanism “double-diagonal -members-failure-belt” of this type of structure was discovered. Based on the analysis results, the measures to mitigate collapse of this type of structure is recomanded.

The paper aims in assessing risks of mandible fractures consequent to impacts or sport accidents. The role of the structural stiffness of mandible, related to disocclusion state, is evaluated through numerical simulations using the finite element method (FEM). It has been assumed that the quasi-static stress field, due to distributed forces developed during accidents, could explain the common types of mandibular fractures. Geometric model of adolescent mandible was built, upon the basis of medical imaging, in CAD software with distinction between cortical layer and inner spongy bone. The finite element model of disoccluded mandible was next created. Mandibular condyles were supposed jammed in the maxillary fossae. The total force of 700 N, simulating an impact on mandible, has been sequentially applied in three distinct areas: centrally, at canine zone and at the mandibular angle. Clinically most frequent fractures of mandible were recognized through the analysis of maximal principal stress and maximal principal strain fields. Mandibular fracture during accidents can be analyzed at satisfactory level using linear quasi-static FE models for designing protections in sport and transport. The proposed approach can be improved by introducing more realistic interactions between condylar processes and fossae.

Due to the demands of personal travels and entertainments, light airplanes and small business aircrafts are developing rapidly. Light airplane structure is simple; however, it lacks crashworthiness design, especially the considerations on the impact of energy absorption. Therefore, in an event of accident, significant damage to passengers will be usually incurred. Airplanes made of composite materials structurally have high specific strength and good aerodynamic configuration. These materials have become the primary choice for new airplane development. This study mainly explores the topology optimization analysis of the light aircraft’s cockpit made of carbon fiber reinforced composites. This paper compares the compression amounts in the original models of composite material and aluminum alloy fuselages with the models after optimization during the crash-landing, in order to investigate the safety of fuselages made of different materials after structural optimization under the dynamic crashing. This study found that the energy absorbed by the aluminum alloy fuselage during crash-landing is still higher than that by the carbon fiber reinforced composites fuselage. On the other hand, the aluminum alloy fuselage after topology optimization could have an energy absorption capability enhanced by 40%, as compared to the that of the original model of aluminum alloy fuselage.

This paper investigates the classical but yet intriguing question of modeling orthotropic RC floor systems endowed with lightening elements, by means of a uniform orthotropic slab, in order to describe accurately the building response under seismic loads. The simplified method, commonly adopted by engineers, based on the equivalence between the transverse stiffness of the RC elements of the real floor and those of the orthotropic slab is presented. A case study in which this simplified method is used is also provided. Then, an advanced finite elements (FE) based method to determine the elastic properties of the equivalent homogenized orthotropic slab is proposed. Based on the results obtained from the application of this method to a case study, a discussion on the adequacy of the simplified method is also provided.

This study examines stress distributions in adhesive joints under various loading and temperature conditions. Finite element analysis (FEA) was employed to compute the peel and shear stresses at the adhesive interface and bondline mid-section. Dependency analysis shows that mid-section peel stress significantly impacts the experimental shear strength of SLJs more than shear stress. This insight highlights the need to carefully analyze peel stress and bending moment factors. The analytical solutions proposed by Goland and Reissner were analyzed with modifications by Hart-Smith and Zhao. Hart-Smith’s approach performed more effectively, especially when the adhesive layer thickness (t_a) is 0.5 mm and the overlap length to thickness ratio (c/t_a) is ≥20. FEA revealed stress distributions at the adhesive/adherend interface and bondline mid-section. DP490 adhesive joints exhibited lower stresses than EA9696. Temperature variations significantly affected joint behavior, particularly above the adhesive’s glass transition temperature (T_g). Both EA9696 and DP490 adhesive joints displayed distinct responses to stress and temperature changes. The parabolic and biquadratic solutions for functionally graded adhesive (FGA) joints were compared. The biquadratic solution consistently yielded higher shear and peel stress values, with an increase ranging from 15% to 71% compared to the parabolic solution at various temperatures because of larger gradient of the Young's modulus distribution near the overlap ends. The ratio of peak peel stress to peak shear stress suggests selecting an adhesive with a superior peel strength or primarily reducing the peak peel stress by functionally grading is advisable, particularly if the adhesive is brittle. Comparison of stress distributions emphasize the importance of selecting adhesives based on stress type, temperature, and solution methods in optimizing adhesive bonding applications. These findings provide valuable insights for thermomechanical applications where thermal stimuli may be used for controlled debonding.

Dental implantology has evolved greatly as a result of the application of biomechanical engineering concepts, notably the utilization of additive manufacturing technologies that recreate the intricate architecture of real bone. This technique intends to improve functional compatibility with bone tissue while also addressing continuing problems of stress distribution inside dental implants. The research examines two types of dental implants: fully gyroid latticed and hybrid gyroid latticed with a solid neck, comprising three different cell sizes—FI-111, FI-222, FI-333, and their hybrid equivalents HI-111, HI-222, and HI-333. Traditional solid implants often cause stress shielding, which reduces long-term osseointegration and implant durability. Developed using nTopology for lattice configuration and meshing, this design undergoes a comprehensive stress analysis using finite element analysis (FEA). The study seeks to determine the most efficient implant design by maximizing stress distribution and mechanical stiffness. Analyses suggest that hybrid latticed implants, notably the HI-222 type, possess superior mechanical characteristics. These implants successfully balance stiffness and flexibility, reducing stress concentrations and improving stress distribution throughout the implant structure. This balance is critical for decreasing micromotions at the bone-implant interface, boosting osseointegration, and, ultimately, increasing implant longevity and success rates in clinical settings. The paper describes a substantial redesign of dental implants by combining biomechanical engineering and modern additive manufacturing. The novel proposed hybrid structure not only addresses typical challenges such as stress shielding, but it also establishes new standards in implantology, potentially leading to better clinical results and patient quality of life. This complete approach emphasizes the ability of personalized implant designs to fulfill individual clinical demands, indicating a significant improvement in dental treatment technology.

Wave propagation in buildings analyzed from its response to earthquakes can be applied for tracking the dynamic characteristics changes based on impulse response functions. The velocity of traveling shear-waves and intrinsic attenuation in buildings can be retrieved for buildings. It is an approach of importance for system identification. The Factor Building at the University of California, Los Angeles campus (henceforth referred as the UCLA Factor Building), an instrumented 15-story steel-moment frame structure, are selected for dynamic response characteristics identification. They are inferred from wave propagation using seismic interferometry for recorded motions. The deconvolved waves are used to compute shear wave travel time and wave attenuation. The natural logarithm for the envelopes of the waveforms deconvolved with the signal at the basement constrains the attenuation. The waveforms from deconvolution with the motion at the basement suggest the fundamental mode of the building. The frequency and decay with time constrain the shear velocity and attenuation as well. The velocity are measured using arrival times picked from these deconvolved waves and the average velocity is 147.1 m/s. The quality factor is equal to 10.8 with the damping ratio at 5%. The shear wave velocity and damping ratio estimated from the deconvolved waves have a good match with these inferred from the waveforms deconvolved with the basement recording. This consistency suggests that wave deconvolution can be applied to extract the building response from the excitation and ground coupling. Based on the resonant frequencies and damping ratio, the finite element model of the building is updated. This study implies that wave deconvolution is powerful to extract the building dynamic characteristics. It helps to better the understanding of the dynamic response of buildings to earthquakes.

The paper presents the problem of damage detection in steel girders. Static displacements at the selected point of the structure play the role of measured variables. Structural response signal decomposition is performed according to the Mallat pyramid algorithm, which is used to perform the discrete wavelet transform (DWT). This procedure allows us to quite well determine the location of structural damage. The geometry and the placement of any defective part of the structure may have a random character. It can be assumed that the random processes occurring in the broadly understood structure mechanics are Gaussian in nature. The first four probabilistic moments are estimated using three approaches independent: semi-analytical (SAM), perturbative (SPT), and Monte-Carlo simulations (MCS). The semi-analytical random approach seems to be the most optimal due to the necessary computation time. The incorporation of the mathematical stochastic apparatus into the classical (deterministic) analysis of the statics of the structure makes it possible to estimate the reliability measures of the analyzed girder.

Leaf springs are suspension components constituted of several structural and assembly components. Once the leaves are assembled, a stress state due to the assembly process is developed. Thus, the quantification of the assembly stress level should be known in order to be accounted for in the structural design. The aim of this article is precisely to quantify the assembly stress level in parabolic leaf springs for railway freight wagons. An experimental method with the aid of extensometry is considered to evaluate elastic assembly stresses. Next, a numerical model of the leaf spring assembly process with the required boundary conditions is developed to support the experimental results and allow obtaining a distribution of surface stresses along the master leaf length. The development of this work resulted in a surface stress profile along the length of the master leaf generated by the assembly process, where the deviations of the numerical outcomes regarding the analysed spots with strain gauges are mostly of little relevance.

Modal testing is a common step in aerostructure design, serving to validate the predicted natural frequencies and mode shapes obtained through computational methods. The strategic placement of sensors during testing is crucial to accurately measuring the intended natural frequencies. However, conventional methodologies for sensor placement are often time-consuming and involve iterative processes. This study explores the potential of machine learning techniques to enhance sensor selection methodologies. Three machine learning-based approaches are introduced and assessed, comparing their efficiency with established techniques. The evaluation of these methodologies is conducted using a numerical model of a beam to simulate real-world scenarios. The results offer insights into the efficacy of machine learning in optimizing sensor placement, presenting an innovative perspective on enhancing the efficiency and precision of modal testing procedures in aerostructure design.

Currently, in some domestic and foreign mines, the backfill drilling pipeline experience a rupture phenomenon even when the wear degree is low. This results in a delay in production due to the filling becoming 'sick'. This paper presents, for the first time, the damage mechanism from a mechanical perspective and re-derives the anti-extrusion strength model of the backfill drilling pipeline. We investigate the influence of law on the anti-extrusion strength of pipelines from the perspective of stratum and cement ring. We then verify the theoretical and simulation results through engineering examples. The results demonstrate that the Mises stress criterion is a suitable modification principle for the anti-extrusion strength model of the backfill drilling pipeline. The anti-extrusion strength of the pipeline is related to the elastic modulus and Poisson's ratio of the stratum, and the thickness of the cement ring. It is negatively affected by the depth of the stratum. For hard strata, a cement ring with a smaller elastic modulus is suitable, while for soft stratum, a cement ring with a larger elastic modulus is recommended. When the missing angle of the cement ring is less than 60°, the stress concentration factor increases up to 2.2. The stress unloading capacity of the cement ring ranges from 32.7% to 37.8%, and optimal performance of the cement ring is achieved when it has high strength and low rigidity. The backfill filling pipeline of a copper mine abroad was destroyed due to external extrusion force exceeding its anti-extrusion strength value. The modified pipeline anti-extrusion strength model is 18.2% higher than the casing API strength value. This finding can inform the design of the backfill filling pipeline for China's kilometer-deep wells in the future.

It is commonly recognized that a ship’s ability to maintain stability is dependent upon the way its loads are distributed. This means making sure the center of gravity is sufficiently low and evenly distributed. Yet when a ship encounters turbulent water with extremely large wave amplitudes, load distribution is a critical factor in determining how stressed the ship becomes. In particular, this article examines several methods for increasing the ship's longitudinal strength to prevent ship structural damage during hogging and sagging. This research presents a number of proposed techniques to enhance the longitudinal strength of cargo ships in order to guarantee minimal stress and deformation during sagging and hogging. The first method is the stern-to-bow use of curved support structures to increase the longitudinal strength of a cargo ship. And another way to stiffen the lateral plate is the proposed use of different designs of the side plate, such as wavy curves, stars, etc. To further boost the strength, the lateral plate upper part is suggested to curve upward. A numerical simulation using ANSYS of the stiffened plate when subjected to uniform load distribution and fixed support on either end shows that there is a decrease in deformation of the plate when we use wavy curved plates compared with those of the other plate stiffening methods used, and furthermore, curving the plate on top increases its strength during sagging as well as hogging. Finally, the use of a curved support structure further enhanced the overall longitudinal strength of the vessel.

In this study, a new structure is proposed based on the body-centered cubic (BCC) lattice structure: adding a cubic truss in the center of the BCC structure, which is denoted by TLC (Truss-lattice-cube). The different dimensions of the central cube can affect the mechanical properties of the lattice structure notably. With a fixed length (15mm) of a unit cell, the optimal size for the central cube is determined to be 5 mm. Quasi-static compressive tests are performed on specimens made of Polylactic Acid (PLA) using additive manufacturing technology. The deformation characteristics of the new structure is analyzed in detail by experiments and numerical simulations. Compared to the BCC structure, the mechanical properties of the TLC structure are significantly improved. The initial flow stress of the TLC has been increased by 122% at a strain of 0.1; the specific strength has been enhanced by 293% at a strain of 0.5; and the specific energy absorption has been improved by 312% at the strain of 0.6. Printing defects in the lattice structure could deteriorate the mechanical property remarkably. In this work, incorporation of microcracks into the finite element model makes the simulation capture the influence of printing defects and improve the predictive accuracy of simulation significantly.

This paper presents a numerical study on the laser shock wave propagation in a 3D woven CFRP material by means of detailed and homogenized FE models. The aim of the study is to numerically characterize the shock wave response of the 3D woven CFRP in terms of back face velocity profiles and the induced damage and to investigate whether the detailed FE models could be replaced by homogenized FE models. The 3D woven geometry was designed using the TexGen software, while the numerical analyses were executed using the LS-Dyna explicit FE software. A high-strain rate behavior was considered for the matrix. The fiber bundles in the detailed models were modeled as a high fiber content unidirectional composite laminate, with its mechanical properties calculated by micromechanical equations. Both for the fiber bundles of the detailed model and the homogenized models a progressive damage material model was applied. The results of the detailed model reveal a considerable effect of material’s architecture on the shock wave propagation and the sensitivity of the back face velocity profile to the spot location. Consequently, the homogenized model is not capable of accurately simulating the shock wave response of the 3D woven composite. Moreover, the detailed predicts with high accuracy matrix cracking in the resin-rich areas and in the bundles as well as fiber failure. On the contrary, the homogenized model predicts matrix cracking at the same areas and no fiber failure.

We investigated the biomechanical behaviour of different fixations of the tibial posterior malleolus (TPM), simulating distinct situations of involvement of the tibiotalar articular surface (TTAS) through a finite element model (FEM). A 3D computer-aided design model of the left ankle was obtained. The materials used were divided according to their characteristics into ductile and non-ductile, and all materials were assumed to be linear elastic, isotropic, and homogenous. Three different fracture lines of the TPM were defined, with sagittal angles of 10º, 25º, and 45º. For biomechanical comparison, different constructions using a trans-syndesmotic screw (TSS) only (Group T), a one-third tubular plate only with (Group PT) and without (Group PS) a TSS, and a locked compression plate with (Group LCPT) and without (Group LCPS) a TSS were tested. FEM was used to simulate the boundary conditions of vertical loading. Load application regions were selected in the direction of the 700N Z axis, 90% on the tibia and 10% on the fibula. Data on the displacement and stress in the FEM were collected, including the total principal maximum (MaxT) and total principal minimum (MinT) for non-ductile materials, total displacement (desT), localized displacement at the fragment (desL), localized displacement at syndesmosis (desS), and Von Mises equivalent stress for ductile materials. The data were analyzed using ANOVA and multiple comparison LSD tests were used. For TPM fractures with a sagittal angles 10º and 25º, desL in the PT and LCP groups was significantly lower, as well as Von Mises stress in Group LCPT in 10º, and PT and LCPT groups in 25º. TPM fractures with a sagittal angle of 45º, desL in the LCP group and Von Mises stress in Group LCPS and LCPT was significantly lower. We found that any TPM fracture may indicate instability of the distal tibiofibular syndesmosis, even when the fragment is small. Our study showed that in fragments involving 10% of the TTAS, the use of a TSS is sufficient, but when the involvement is greater than 25% of the TTAS, either non-locked or locked plate must be used to buttress the TPM. In posterior fragments affecting 45% or more of the TTAS, the use of a locking plate is recommended.

Precipitates are the primary source of strength for the Al-Mg-Si alloy. Aluminum alloy in the peak-aged state mainly contain β'' and β' precipitates. In the past, most of the literature only considered the strengthening effect of β''. Here we develop a single crystal intensity model including both precipitate enhancement effects for the first time. This model was subsequently implemented into a crystal plastic finite element method to model the uniaxial tensile process of a polycrystalline aggregate model of Al-Mg-Si alloy. The simulation results for uniaxial stretching are in good agreement with the experimental results, confirming that the constitutive parameters used for the single crystal strength model with two precipitates are based on realistic physical implications. Furthermore, by comparing the uniaxial tensile simulation results of a peak-aged alloy considering the actual precipitated phase composition of the alloy with those assuming that the precipitated phase is only the β'' phase, the predicted tensile strength of the former is around 5.65% lower than that of the latter, suggesting that the two kinds of precipitation should be considered separately when simulating the mechanical response of Al-Mg-Si alloy. It is highly expected that the present simulation strategy is not limited to the Al-Mg-Si alloys, and it can be equally applied to the other age-enhanced alloys.

The purpose of this study was to investigate the functional load capacity of the periodontal ligament (PDL) in a full arch maxilla and mandible model using numerical simulation. The goal of was to determine the functional load pattern in multi- and single rooted teeth with full and reduced periodontal support. CBCT data were used to create 3D-models of a maxilla and mandible. The DICOM dataset was used to create a CAD-model. For a precise description of the surfaces of each structure (enamel, dentin, cementum, pulp, PDL, gingiva, bone), each tooth was segmented separately, and the biomechanical characteristics were considered. A Finite-Element-Analysis (FEA) software computed the biomechanical behavior of stepwise increased force of 700N in cranial and 350N in ventral direction of the muscle approach of the Masseter muscle. The periodontal attachment (cementum-PDL-bone contact) was subsequently reduced in 1mm increments and the simulation repeated. Quantitative (pressure, tension, and deformation) and qualitative (color-coded images) data were recorded and descriptively analyzed. The teeth with the highest load capacities were the upper and lower molars (0.4-0.6MPa), followed by the premolars (0.4-0.5MPa) and canines (0.3-0.4MPa) when vertically loaded. Qualitative data showed that the area with the highest stress in the PDL were for single rooted teeth the cervical and apical area and for molars additionally the furcation roof. In both, single- and multi rooted teeth the gradual reduction of the bone levels caused an increase of the load on the remaining PDL. Cervical and apical areas as well as the furcation roof are the zones with the highest functional stress. The more bone loss, the higher the mechanical load on the residual periodontal supporting structures.

According to the development needs of marine seismic station networking, it is necessary to develop a marine seismometer with a networking function to realize quasi-real-time and large-scale long-term observation of undersea seismic signals. This can lay a solid foundation for marine seismic research and research on earth tectonics, activities, and processes. Based on the advantages and disadvantages of the existing observation equipment, this study proposed to develop a new type of the Autonomous Glide Marine Seismometer(AGMS) that can be applied to the global sea area, focusing on the design and study of the main structure of the carrier of the AGMS. The ultimate strength of the ellipsoidal pressure-resistant spherical shell is investigated by using the finite element analysis method. A reasonable finite element analysis model is proposed through comparing the finite element parameter results with the theoretical formulas and relevant specifications. Finally, the design process and strength checking of the pressure-resistant shell of the AGMS are established. These relevant results can also provide a reference for the design of other deep-sea pressure-resistant chambers.

The article reports on a patient who required a cranial protection system.Using additive manufacturing techniques and surgical planning with the help of bio-models, a patient-specific bone implant solution was proposed that allows aesthetic restoration of the affected area and provides an adequate level of protection. In addition, through a comparative analysis with finite elements, the mechanical response to external actions of the medical device, printed with two materials: polymethylmethacrylate (PMMA) and polyether-ether-ketone (PEEK), is simulated. The tested materials have recognized biocompatibility properties, but their costs on the market differ significantly. The results obtained demonstrate the similarities in the responses of both materials. It offers the possibility that low-income people can access these devices, guaranteeing adequate biomechanical safety, considering that PMMA is a much cheaper material than PEEK.

We have much to learn from historic buildings. The fact that they are still standing is obvious, but unfortunately often not appreciated. In this paper the results obtained from a series of parametric analyses, where the influence that geometric and mechanical parameters have in the structural response of existing cob walls within an Irish context, are presented. A design of experiments using central composite designs was implemented along with analysis of variance following two computational approaches, namely, finite element method and kinematic limit analysis. As results, a series of response surfaces and parametric equations with which it is possible to compute safety factors and collapse multipliers (within the range of values studied) are provided. Based on the results obtained, it could be concluded that traditional cob walls in Ireland are very robust. Relatively high acceleration values, unlikely to happen in a low seismic hazard region such as Ireland, would be needed to start the collapse mechanisms studied or cause yielding in typical vernacular cob walls. Furthermore, the equations generated with the refined regression models, can be used by practitioners as a first approach to estimate the safety levels of existing cob buildings with similar characteristics.

Multi-chiplet technique is expected to be a promising solution to achieve high-density system integration with low power consumption and high usage ratio. This technique can be integrated with glass interposer to accomplish a competitive low fabrication cost compared with the silicon-based interposer architecture. In this study, process-oriented stress simulation is performed by the element activation and de activation technique in finite element analysis architecture. Submodeling technique is also utilized to mainly conquer the scale mismatch and difficulty in mesh gridding design. It is used as well to analyze the thermomechanical responses of glass interposer with chiplet arrangement and capped epoxy molding compound (EMC) during curing. A three-factor, three-level full factorial design is applied using the analysis of variance method to explore the significance of various structural design parameters on stress generation. Analytic results reveal that the maximum first principal stresses of 130.75 and 17.18 MPa are introduced on the sidewall of Cu-filled via and the bottom of glass interposer, respectively. Moreover, the EMC thickness and through glass via pitch are the dominant factors in the adopted vehicle. They significantly influence the stress magnitude during heating and cooling.

Air conditioners are indispensable necessities in daily life and industrial production. However, the heat generation of their internal core power components can limit their performance release. In this work, we explore the application of thermoelectric coolers (TECs) in the field of power device heat dissipation through finite element simulation. First, we geometrically modeled the structure and typical operating conditions of core power devices in air conditioners. We compared the temperature fields in air cooling and TEC active cooling modes for high power consumption power devices in a 319 K operating environment. The simulation results show that in the single air cooling mode, the maximum temperature of the 173.8 W power device reaches 394.4 K and the average temperature reaches 373.9 K, which exceeds its rated operating temperature of 368.1 K. However, after adding TEC, the maximum temperature of the power device drops to 331.8 K at an operating current of 7.5 A and the average temperature of the device is 326.5 K. It indicates that TEC active cooling has a significant effect on the temperature control of the power device. At the same time, we also studied the effect of TEC working current on the temperature control effect of power devices. For TEC, there is a minimum working current of 5 A; when it is less than 5 A, TEC has no cooling effect; the cooling effect of TEC increases with an increase in working current. When the TEC working current is 10 A, the average temperature of the power device can be reduced to 292.9 K. This study has made a meaningful exploration of the application of TEC in chip temperature control and heat dissipation, providing a new solution for chip thermal management and accurate temperature control.