Search

The potential of paper and paperboard as fiber-based materials capable of replacing conventional polymer-based materials has been widely investigated and evaluated. Due to paper's limited extensibility and inherent heterogeneity, local structural variations lead to unpredictable local mechanical behavior and instability during processing, such as mechanical forming. To gain a deeper understanding of the impact of mechanical behavior and heterogeneity on the paper forming process, the Finite Element Method (FEM) coupled with continuum modeling is being explored as a potential approach to enhance comprehension. To achieve this goal, utilizing experimentally derived material parameters alongside stochastic finite element methods allows for more precise modeling of material behavior, considering the local material properties. This work first introduces the approach of modeling heterogeneity or local material structure within continuum models, such as the Stochastic Finite Element Method (SFEM). A fundamental challenge lies in accurately measuring these local material properties. Experimental investigations are being conducted to numerically simulate mechanical behavior. An overview is provided of experimental methods for material characterization, as found in literature, with a specific focus on measuring local mechanical material structure. By doing so, it enables the characterization of the global material structure and mechanical behavior of paper and paperboard. [ABSTRACT FROM AUTHOR]

The present article provides a comprehensive analysis of the aging phenomenon of transformer insulation and the corresponding changes in insulation performance that occur over time. The present manuscript outlines diverse insulating materials that are employed in transformers, the fundamental aging mechanism of the insulation, the function of space charges in the aging of the insulation, and various simulation methodologies utilized for diagnosing the insulation status of the transformer. The primary focus of this paper is to examine the various numerical techniques employed in modeling and analyzing the aging characteristics of transformer insulation. The mechanisms underlying the accumulation of charge and distortion of electric fields in insulation are explained through the utilization of numerical modeling techniques. This paper focuses on understanding the mechanism of insulation aging, thereby facilitating the optimization and maintenance of insulation performance throughout its operational lifespan. [ABSTRACT FROM AUTHOR]

The Proton Exchange Membrane Fuel Cell (PEMFC) performances are strongly impacted by the compression of the Gas Diffusion Layer (GDL). Despite its fibrous microstructure, this material is usually considered as a continuous medium and characterized with uniform loading. However, the GDL is subjected to a heterogeneous compression onto rib/channel patterns in the fuel cell assembly. In the present study, a complex behavior of the GDL response is experimentally revealed when the material is loaded with a rib/channel pattern, compared to uniform compression. The tests are simulated by finite element modeling using a classical strain-dependent elastic law, using parameters fitted from uniform compression experiments. It is shown that the numerical results do not reproduce the effect of pattern observed experimentally. Hypotheses to interpret these results involve mechanisms at the fiber microscale including fiber fracture, cross-link breakage and fiber rearrangement, which are exacerbated by larger material deformation caused by the heterogeneous loading. • Unexpected impact of rib/channel patterns on GDL compression. • Experimental investigation of the effect of rib/channel width on GDL deformation. • Numerical study with classical model unable to reproduce experimental results. • Interpretations based on mechanisms at the fiber microscale. [ABSTRACT FROM AUTHOR]

Alumina-based multilayer laminates with layers of varied porosities and ceramic fiber reinforcements, which can be achieved by including paper-derived ceramic layers, are interesting candidates for chemically stable and thermal shock resistant materials. In order to investigate the complex sintering behavior of these laminates, the sintering of each alumina layer contained within a laminate, each one with its individual sintering anisotropy and sintering shrinkage, has to be considered separately. For this purpose, sintering curves are determined by dilatometry measurements, serving a consecutive examination by an analytical model. Modeling parameters specific to each layer, are applied in finite-element-type sintering simulations for symmetrical multilayer laminates. Differences in the sintering shrinkages may cause internal stresses at layer interfaces and, in turn, may influence the thermal and mechanical properties of multilayer laminates after sintering. Accordingly, the prediction of developing internal stresses within alumina-based multilayer laminates by the presented sintering simulations is a key aspect of this study. [Display omitted] • Individual laminates sintering dilatometry. • Anisotropic sintering assessment by analytic modeling. • Finite element simulation of multilayer laminates sintering. [ABSTRACT FROM AUTHOR]

Femoral deformities are highly prevalent in children with cerebral palsy (CP) and can have a severe impact on patients' gait abilities. While the mechanical stress regime within the distal femoral growth plate remains underexplored, understanding it is crucial given bone's adaptive response to mechanical stimuli. We quantified stresses at the distal femoral growth plate to deepen our understanding of the relationship between healthy and pathological gait patterns, internal loading, and femoral growth patterns. This study included three-dimensional motion capture data and magnetic resonance images of 13 typically developing children and twelve participants with cerebral palsy. Employing a multi-scale mechanobiological approach, integrating musculoskeletal simulations and subject-specific finite element analysis, we investigated the orientation of the distal femoral growth plate and the stresses within it. Limbs of participants with CP were grouped depending on their knee flexion kinematics during stance phase as this potentially changes the stresses induced by knee and patellofemoral joint contact forces. Despite similar growth plate orientation across groups, significant differences were observed in the shape and distribution of growth values. Higher growth rates were noted in the anterior compartment in CP limbs with high knee flexion while CP limbs with normal knee flexion showed high similarity to the group of healthy participants. Results indicate that the knee flexion angle during the stance phase is of high relevance for typical bone growth at the distal femur. The evaluated growth rates reveal plausible results, as long-term promoted growth in the anterior compartment leads to anterior bending of the femur which was confirmed for the group with high knee flexion through analyses of the femoral geometry. The framework for these multi-scale simulations has been made accessible on GitHub, empowering peers to conduct similar mechanobiological studies. Advancing our understanding of femoral bone development could ultimately support clinical decision-making. • Distal femoral growth plate stresses quantified with subject-specific simulations. • Large dataset of 50 femurs – 26 healthy, 24 femurs of participants with cerebral palsy. • Knee flexion angle is of high relevance for typical bone growth at the distal femur. • Workflow reveals plausible growth results supported by geometrical analysis. • Multi-scale workflow based on free software published on GitHub. [ABSTRACT FROM AUTHOR]

Many plant materials in nature have the ability to change their shape to respond to external stimuli, such as humidity or moisture, to ensure their survival or safe seed release. A well-known example for this phenomenon is the pinecone, which is able to open its scales at low humidity due to the specific bilayer structures of the scale. Inspired by this, we developed a novel humidity-driven actuator based on paper. This was realized by the lamination of untreated paper made from eucalyptus fibers to a paper–carboxymethyl cellulose (CMC) composite. As observed, the hygroexpansion of the composite can be easily controlled by the amount of CMC in the impregnated paper sheet, which, thus, controls the morphologic deformation of the paper bilayer. For a more detailed understanding of these novel paper soft robots, we also studied the dynamic water vapor adsorption, polymer distribution and hygroexpansion of the paper–polymer composites. Finally, we applied a geometrically nonlinear finite element model to predict the bending behavior of paper bilayers and compared the results to experimental data. From this, we conclude that due to the complexity of structure of the paper composite, a universal prediction of the hygromorphic behavior is not a trivial matter. [ABSTRACT FROM AUTHOR]

The dynamic mechanical behavior and cushioning performance of honeycomb sandwich panels, which are extensively employed in product cushioning packaging due to their exceptional energy absorption capabilities, were examined using a combination of experimental and numerical methods. Several factors, such as maximum acceleration–static stress, cushioning coefficient–static stress, and other curves, were analyzed under various impact conditions. The simulated stress–strain, deformation modes, cushioning coefficients, and other parameters demonstrate consistency with the experimental results. The acceleration, maximum compression, and cushioning coefficient obtained from the experiment and simulation calculation were 30.68 g, 15.44 mm, and 2.65, and 31.96 g, 14.91 mm, and 2.79, respectively. The results indicate that all error values were less than 5%, confirming the precision and reliability of the model. Furthermore, the model was utilized to simulate and predict the cushioning performance of honeycomb sandwich panels with different cell structures and paper thicknesses. These results provide a solid basis for enhancing the design of subsequent honeycomb element structures. [ABSTRACT FROM AUTHOR]

The study aims to present the results of paper compression under an axial load. Different heights of samples subjected to compression were taken into account. The main goal of the analysis was to determine experimentally the maximum compression load. In addition, numerical models based on the finite element method (FEM) were validated to refer to empirical results. The performed numerical simulations were founded on Green–Lagrangian nonlinear equations for large displacements and strains. The progressive failure of the compressed orthotropic material after exceeding maximum stresses was based on Hill's anisotropy theory. Nonlinear calculations were conducted by using a typical Newton–Raphson algorithm for achieving a sequence convergence. The accuracy of the developed model was confirmed experimentally in compression tests. The technique of analysing the shape of the compressed paper sample on the basis of images recorded during the measurement was used. The obtained test results are directly applicable in practice, especially in the calculation of the mechanical properties of corrugated cardboard and in determining the load capacity of cardboard packaging. Knowing the maximum compressive stress that packaging paper can withstand allows packaging to be properly designed and its strength assessed in the context of the transport and storage of goods. [ABSTRACT FROM AUTHOR]

Ultrathin devices are rapidly developing for skin‐compatible medical applications and wearable electronics. Powering skin‐interfaced electronics requires thin and lightweight energy storage devices, where solution‐processing enables scalable fabrication. To attain such devices, a sequential deposition is employed to achieve all spray‐coated symmetric microsupercapacitors (μSCs) on ultrathin parylene C substrates, where both electrode and gel electrolyte are based on the cheap and abundant biopolymer, cellulose. The optimized spraying procedure allows an overall device thickness of ≈11 µm to be obtained with a 40% active material volume fraction and a resulting volumetric capacitance of 7 F cm−3. Long‐term operation capability (90% of capacitance retention after 104 cycles) and mechanical robustness are achieved (1000 cycles, capacitance retention of 98%) under extreme bending (rolling) conditions. Finite element analysis is utilized to simulate stresses and strains in real‐sized μSCs under different bending conditions. Moreover, an organic electrochromic display is printed and powered with two serially connected μ‐SCs as an example of a wearable, skin‐integrated, fully organic electronic application. [ABSTRACT FROM AUTHOR]

The aim of this discussion is to establish a connection between the model for sandwich beams presented in Schöftner (Acta Mech 234:2543–2560, 2023) and the so-called Krajcinovic model. Although the two models turn out to be quite similar, Krajcinovic model displays a richer kinematics and has been developed in the literature by following slightly different methodologies with respect to that adopted by Schöftner, the latter having the merit of providing equilibrium equations that allow the computation of the interfacial stresses. We also offer some remarks on the capabilities and accuracies of the models discussed. [ABSTRACT FROM AUTHOR]

This work deals with fracture mechanisms in quasi-brittle materials, focusing on the characterization of the Fracture Process Zone (FPZ) of specimens under tensile load. Particularly, paper was used as model material. Experiments were conducted on notched and unnotched specimens. Based on an image analysis of these observations, a stochastic finite element model was developed, using both a nonlocal stress-based approach and a discretized random field modeling of the Young's modulus. The proposed methodology allowed characterizing the damage zone and the size of the FPZ, analyzing the influence of the mesostructure, composed of flocs (fiber aggregates where the basis weight is larger than the average one) and antiflocs (complement of flocs). The area of the active FPZ and the normalized stress drop were linked using a surface energy dissipated in the active FPZ. The stress drop, until limiting value, increased with the width of the active FPZ. Finally, a relationship between the surface energy and the nonlocal internal length was established. [ABSTRACT FROM AUTHOR]

Today, the toilet paper market offers product types with varying number of plies, providing better mechanical strength and liquid absorption. Several tissue paper perforation systems exist, and the best commonly applied is a top-cutting mechanism that includes an oblique blade, a combined oblique blade, or a simple spiral blade. The perforation efficiency must be high to have an easy sheet separation from the roll of the toilet paper, which does not always occur. Hence, consumer satisfaction can depend on the perforation performance. To study this, a laboratory perforation system was used to perforate different commercial toilet papers (in brands and number of plies) and evaluate their perforation efficiency. A finite element method (FEM) was used to simulate the curve of the progression of perforation efficiency as a function of the cut distance. The main findings were a stabilization of the perforation efficiency from a cut distance of 6 mm and a 15% increase in the cut distance for the laboratory blade to match the industrial perforation efficiency. The FEM analysis confirmed the behavior of the evolution of perforation efficiency with the increase of the cut distance. [ABSTRACT FROM AUTHOR]

Paper has always been of great importance as a printing and packaging material. Due to its potential for recyclability and reusability, which result in a low environmental impact, it is expected that this importance will soar over the coming years. Nevertheless, even though paper has been in use for decades, many research questions are still open, especially in the field of mechanics. One of these questions is the appropriate size of a representative volume element (RVE) of the microstructure of paper. Therefore, a finite element network model of interconnected fibers was developed which represents the microstructure of paper. In order to create the network model realistically, geometry and material quantities were taken from experimental investigations. Additionally, single fibers and fiber bonds were investigated beforehand in order to incorporate their proper behavior into the network. Based on this, the networks were generated and subsequently loaded by applying mixed displacement and traction boundary conditions. The resulting effective stresses were compared to determine the necessary RVE size. Two approaches were considered for modeling the RVE: (i) one single "large" network, and (ii) a certain number of smaller networks whose average yielded a statistically representative effective response. As a result, the second approach turned out to be less computationally intensive due to smaller models and the possibility to run calculations simultaneously. Concluding, it was possible to determine the size of the RVE which will be used to investigate various individual microstructural effects on the macrostructure in the future. [ABSTRACT FROM AUTHOR]

To explore the application of cold-formed thin-walled steel-paper straw board(CTSPSB) composite wall in practical engineering and further meet people's living requirements, it was proposed to open holes in the composite wall to simulate the doors and windows in practical applications. Two composite wall specimens were tested to study the shear performance of the CTSPSB composite wall. Through the analysis of specimens' damage forms and experimental data, the characteristic values of bearing capacity and lateral stiffness were obtained. And then, the model of the composite wall was built by ANSYS, and finite element analysis (FEA) results were consistent with the experimental results, which could verify the feasibility of the finite element model. Moreover, the model needed to open holes and extensive parameter analysis was carried out. The FEA results indicate the most reasonable distance between screws around the opening is 150 mm; the most suitable spacing between the small studs is 400 mm; the position of the opening has the least influence on the shear performance, and the difference between the results of the five groups of models is within 5%; while the width of the opening has the greatest impact on the shear performance. Compared with the wall without opening, the bearing capacity of the wall with an opening width of 600 mm, 1200 mm and 1800 mm decreases by 38%, 46% and 52% respectively. Besides, the calculation formula of shear capacity of CTSPSB composite wall with openings was improved, which could be used as experience for practical engineering. [ABSTRACT FROM AUTHOR]

In this contribution, the causes leading to the undesired unevenness on the exterior surface of a paper‐based sandwich structure during the manufacturing process are investigated. This provides insight into the manufacturing parameters that need to be adjusted to reduce this waviness. A finite element model of the pressing process is developed in the commercial software ABAQUS and then used to analyze the formation of the surface unevenness. The sandwich structure adopted in the model consists of a decorative layer, an intermediate layer, a corrugated core and a bottom layer. In particular, the core has a sinusoidal structure with various geometry parameters, which are adjusted in the model to investigate the effect on the waviness of the top surface. The results of the quasi‐static process simulation correlate well with the experimental results. A subsequent parameter study based on the design of experiments (DoE) shows the relationship between the material parameters of the intermediate/bottom layer and the resulting waviness. [ABSTRACT FROM AUTHOR]

A square spiral incision was introduced into a 15 × 15 cm regular quadrilateral conductive film to study the modulation of electrothermal distribution by kirigami structure. Based on the simulation of finite element method, the transition from thermal focusing to thermal defocusing could be achieved by increasing the length of connection between adjacent modules. The effective electrothermal area can be increased by nearly two times through increasing the length increment of the notch length between adjacent modules. The maximum temperature was decreased from 254 °C to 44 °C, and the temperature difference was adjusted from 76 °C to 4 °C by increasing the number of kirigami module. Through the contact temperature measurement experiment on graphite paper kirigami structure, the trends could be verified. By creating an equivalent circuit model and considering the influence of heat conduction, the peak temperature at the center of the structure and the decrease of the overall temperature could be explained. By deriving the relationship between the temperature of conductive film of kirigami structure with time and position, the electrothermal characteristics of previously reported graphene kirigami structure could be predicted. [ABSTRACT FROM AUTHOR]

A square spiral incision was introduced into a 15 × 15 cm regular quadrilateral conductive film to study the modulation of electrothermal distribution by kirigami structure. Based on the simulation of finite element method, the transition from thermal focusing to thermal defocusing could be achieved by increasing the length of connection between adjacent modules. The effective electrothermal area can be increased by nearly two times through increasing the length increment of the notch length between adjacent modules. The maximum temperature was decreased from 254 °C to 44 °C, and the temperature difference was adjusted from 76 °C to 4 °C by increasing the number of kirigami module. Through the contact temperature measurement experiment on graphite paper kirigami structure, the trends could be verified. By creating an equivalent circuit model and considering the influence of heat conduction, the peak temperature at the center of the structure and the decrease of the overall temperature could be explained. By deriving the relationship between the temperature of conductive film of kirigami structure with time and position, the electrothermal characteristics of previously reported graphene kirigami structure could be predicted. [ABSTRACT FROM AUTHOR]

Paper, a porous-fibrous network, is made up of hydrophilic fibres, which are notably susceptible to deformations due to variations in moisture content. The response of paper sheets to full or partial wetting, based on the water-induced fibre swelling, has mostly been predicted in the past using highly idealised two-dimensional network geometries. The purpose of the present paper is to establish how precisely the two-stage process of determining effective hygro-expansivity under uniform wetting by a two-dimensional model followed by bending analysis using a continuum model predicts the curl of a fully three-dimensional network. In the present multi-scale modelling work, the moisture-induced homogenised sheet scale curl, derived from hygro-elastic deformations using non-uniform moisture profiles, is analysed based on an idealised non-woven three-dimensional fibrous network finite element model. First, the effective paper sheet swelling properties, derived computationally for the uniform wetting condition, are compared with analytically homogenised hygro-elastic properties. Whereas the analytically predicted elastic modulus has an inaccuracy of only ∼ 5%, the hygro-expansivity coefficient is over-predicted by roughly a third for a network with straight fibres. The difference is smaller if waviness of the fibres and wrap-around in the fibre bonds are taken into account. Fibre alignment turns out to have little influence. Subsequently, paper sheet curl predictions made by the analytical model are compared with that computed from the three-dimensional network model. The inaccuracy of the fully analytical model, based on the two-dimensional network description, relative to the fully detailed three-dimensional computational model with straight fibres is on the order of 40%–50%. The fidelity of the simple model is primarily constrained by the inaccuracy of the effective hygro-expansion coefficient. If the expansivity extracted from the three-dimensional simulations under uniform wetting is inserted, the inaccuracy drops to less than 10%. The continuum bending description hence performs adequately, provided it is based on a reliable estimation of the uniform expansion behaviour of the material. • Hygroscopic curl of paper sheets may be predicted using homogenised models. • We assess the accuracy of such predictions. • As a reference simulations of an idealised three-dimensional fibre network are used. • The predictions are accurate provided an accurate hygro-expansivity is used. • Two-dimensional network models may not be sufficient for this. [ABSTRACT FROM AUTHOR]

Embossing is a functional and strategic process for creating high-quality multi-sensory tissue-paper products. Embossing modifies the sheet surface by generating hill and/or valley designs, changing the third-dimension z with a compressive die. This research work specifically concerns the impact study of the engraving finishing geometry on the final properties of tissue paper. This work led us to conclude that, even though the sheets individually present a higher hand-feel (HF) value for the straight finishing geometry, the highest softness was obtained in the two-ply prototype for the round finishing geometry. Moreover, this study confirmed that the HF value reduces with the increase of the bulk, being more accentuated for the micropattern. Relevant differences could not be seen in the spreading kinetics of the liquid droplets over time. Thus, the finishing geometry of the 3D plates did not impact the absorption kinetics on these samples. The finite element model allows us to understand the effect of the plate pattern and its finishing geometry on the paper, and the simulation results were in accordance with the experimental results, showing the same trend where patterns with a round finishing geometry marked the tissue-paper sheet more than patterns with a straight finishing did. [ABSTRACT FROM AUTHOR]

Embossing is a converting process in which the surface of a tissue paper sheet is changed under high pressure, allowing different functions. In this work, the authors intend to study how the embossing pressure affects the main properties of tissue paper, using a laboratory embossing system. An optimum pressure was achieved at 2.8 bar to this embossing laboratory set-up. The effect of pressure when densifying the paper sheet gives it a gain in mechanical strength but no differences in terms of liquid absorbency. The two embossing patterns present different behaviors but both evidence losses in mechanical and softness properties. On the other hand, the finite element method (FEM) does not show clear evidence of how the pressure affects the paper strength. For the deco die, it is possible to observe that the amount of yielding is slightly higher for lower pressure (2.4 bar), but this plasticity state parameter is very similar for 2.8 bar and 3.2 bar. For the micro die, FEM simulations of the manufacturing pressure do not show a considerable impact on the amount of plasticity state of the material; only for 3.2 bar, it shows a change in the pattern of the plasticity state of the paper during the embossing processes. In the end, to achieve a final product with excellent quality, it is important to make a compromise between the various properties. [ABSTRACT FROM AUTHOR]

The embossing operation is one of the processes of tissue paper converting. The embossing parameters influence the final properties of tissue products, such as mechanical, softness, and bulk. In this study, the influence of the rubber hardness used against the embossing steel rolls with a pattern created by intaglio engraving was studied. Three different configurations of rubber plates stacking, each plate with different hardness, were studied. After embossing, mechanical properties, softness, and bulk were evaluated to analyze the effect of rubbers hardness on these properties. Furthermore, a Finite Element Model of the embossing operation was used that considered the same rubber plates stacking configurations used in experiments, and it was able to replicate the experimental results. This work led us to conclude that the configuration where two rubber plates with different hardness, where the rubber plate with higher hardness is in contact with the tissue paper sheet, has shown to be the best solution to obtain higher softness. These findings support the use of embossing operations rubber rolls with a low hardness internal layer and a high hardness external layer in industry. Thus, finite element models were also shown to be reliable tools to virtually test other configurations, such as, for example, three or more rubber plates with different hardness. Since embossing is one of the tissue paper transformation operations with the greatest impact on the key properties of the final product, this study allows the producer to optimize them by varying the hardness of the rubber roll, as well as its configuration. [ABSTRACT FROM AUTHOR]

Alternating current (ac) losses in conductor-on-rounded-core (CORC) cables of YBCO high-temperature superconducting (HTS) tapes are a significant challenge in HTS power applications. This study employs two finite element analysis (FEA) models to investigate the contributions from different ac loss components and provide approaches for reducing ac losses in cables. An FEA model based on the T-A formulation treats the cross section of thin superconducting layers as 1D lines and, therefore, only can predict the ac loss generated by the perpendicular magnetic field. In contrast, the model based on H-formulation can be performed on the actual 2D rectangular cross section HTS tapes to provide the total ac losses generated by magnetic fluxes penetrating from both the edges and surfaces of HTS tapes, although this model requires more computing time and memory. The 1D and 2D simulation models were validated by cross comparing the results from both models and by comparing sub-section and full cross section models. Subsequently, two models relate cable design and operational parameters to the surface and edge losses of a two-layer CORC cable by considering the (1) relative contributions of edge and surface losses to the overall ac losses; (2) effect of the current distribution between inner and outer HTS layers on ac losses; (3) impact of the tape alignment on ac losses in each HTS layer; (4) influence of the thickness of HTS layers on ac losses; (5) effect of size and number of inter-tape gaps on ac losses; and (6) contribution frequency on the ac losses. The research results given in this paper are therefore not only valuable to suggest strategies for reducing ac loss in multi-layer cables but also for developing more accurate and effective methods to calculate ac loss in CORC HTS cables. [ABSTRACT FROM AUTHOR]

The influence of nonuniform aging has not been considered in the prediction of the frequency domain spectroscopy (FDS) of the transformer oil-paper insulation (OPI) when using the finite element method (FEM). Given this issue, the modified simulation model was constructed based on the XY model and FEM. In the current work, the dielectric response property of the pressboards with different aging conditions in the modified model was characterized by its dielectric response parameters measured in the lab. Then the modified simulation model constructed in a 3-D system would be applied to study the effect of both nonuniform aging and the axial geometry on the FDS of OPI. Finally, the feasibility of the proposed models was verified by comparing the measured FDS and the predicted FDS. This article is expected to provide a model for predicting the FDS of the transformer OPI including the effect of nonuniform aging. [ABSTRACT FROM AUTHOR]

It is pointed out that the Levinson-Reddy beam theory should not be used as a higher-order model for the core, although it considers cross-sectional warping: this beam model is applicable for zero shear stress conditions at the core-skin-interfaces only, which is usually not the case for sandwich structures. The resulting sandwich beam model in [[2]] is a special solution from the Krajcinovic-Bardella derivation, which allows for richer kinematics because it considers zig-zag kinematics. Author's response to "Discussion on the paper by Schoeftner, J., "A verified analytical sandwich beam model for soft and hard cores: comparison to existing analytical models and finite element calculations", Acta Mech, 234, 2543-2560 (2023)" by Lorenzo Bardella. [Extracted from the article]

Thin-walled tubular deployable composite boom (DCB) has gained significant attention due to its lightweight, simple structure and high packaging efficiency. The experimental and numerical simulation methods were employed in this paper to investigate the folding behavior of the DCB. Using T700/epoxy unidirectional reinforced prepreg as raw material, DCB specimens were prepared by the vacuum bag method. To complete the folding experiments of DCB specimens, a folding mechanism was designed and manufactured. Folding experiments of DCB specimens were conducted and folding moment versus rotational displacement curves were measured. In addition, a Finite Element Model (FEM) was established to predict the folding behavior of the DCB. Different failure criteria were considered in the numerical analysis, including the Tsai‐Hill criterion and maximum stress criterion. Prediction results using the FEM were compared with experimental results, and both sets of results were in good agreement. It is shown that the DCB can achieve the folding function without failure. The research results in this paper are of great significance to the practical engineering application of the DCB. [Display omitted] • Six tubular deployable composite boom specimens and a precise folding mechanism were designed and prepared. • Folding experiments were conducted to measure the folding moment and rotational displacement curves. • A finite element model for predicting the folding behavior of the tubular deployable composite boom was established. [ABSTRACT FROM AUTHOR]

The thermal response in paper has been studied by thermography. It was observed that an inhomogeneous deformation pattern arose in the paper samples during tensile testing. In the plastic regime a pattern of warmer streaks could be observed in the samples. On the same samples digital image correlation (DIC) was used to study local strain fields. It was concluded that the heat patterns observed by thermography coincided with the deformation patterns observed by DIC. Because of its fibrous network structure, paper has an inhomogeneous micro-structure, which is called formation. It could be shown that the formation was the cause of the inhomogeneous deformations in paper. Finite element simulations was used to show how papers with different degrees of heterogeneity would deform. Creped papers, where the strain at break has been increased, were analysed. For these paper it was seen that an overlaid compaction of the paper was created during the creping process. During tensile testing this was recovered as the paper network structure was strained. [ABSTRACT FROM AUTHOR]

Optical eigenstates with a high quality (Q) factor provide substantial advantages for a broad spectrum of optical devices, particularly those demanding strong light–matter interactions. Recently, it has been demonstrated that ultrahigh- Q resonances can be realized in planar photonic structures by merging multiple bound states in the continuum (BICs) in the momentum space. Photonic lattices with thin-film geometry are known to support abundant TE-polarized and TM-polarized BICs. While prior research has explored the merging of TE-polarized BICs, this paper presents analytical and numerical results concerning the merging of TM-polarized BICs in laterally periodic one-dimensional photonic lattices. As the thickness of photonic lattices increases, TM-polarized accidental BICs descend along the dispersion curves and eventually merge at the upper edge of the second stop band. Employing coupled-mode analysis, we calculate the analytical merging thickness at which multiple TM-polarized BICs come together at the second-order Γ point. We confirm the merging of TM-polarized BICs through finite-element method simulations. Our results can be beneficial for achieving ultrahigh- Q resonances through the merging of BICs. [ABSTRACT FROM AUTHOR]

Wrinkling, creasing, and folding are frequent phenomena encountered in biological and man-made bilayers made by thin films bonded to thicker and softer substrates often containing fibers. Paradigmatic examples of the latter are the skin, the brain, and arterial walls, for which wiggly cross sections are detected. Although experimental investigations on corrugation of these and analog bilayers would greatly benefit from scaling laws for prompt comparison of the wrinkling features, neither are they available nor have systematic approaches yielding to such laws ever been provided before. This gap is filled in this paper, where a uniaxially compressed bilayer formed by a thin elastic film bonded on a hyperelastic fiber-reinforced substrate is considered. The force balance at the film–substrate interface is here analytically and numerically investigated for highly mismatched film–substrates. The onset of wrinkling is then characterized in terms of both the critical strain and its corresponding wavenumber. Inspired by the asymptotic laws available for neo-Hookean bilayers, the paper then provides a systematic way to achieve novel scaling laws for the wrinkling features for fiber-reinforced highly mismatched hyperelastic bilayers. Such novel scaling laws shed light on the key contributions defining the response of the bilayer, as it is characterized by a fiber-induced complex anisotropy. Results are compared with finite element analyses and also with outcomes of both existing linear models and available ad hoc scalings. Furthermore, the amplitude, the global maximum and minimum of ruga occurring under increasing compression spanning the wrinkling, period doubling, and folding regimes are also obtained. [ABSTRACT FROM AUTHOR]

The ability to determine the condition of paper (e.g. its brittleness) and its permanence is a need in libraries and archives. This study investigates various bend tests and applies these tests to aged paper with the goal of finding easy, non-destructive tests to determine the mechanical properties of paper. The cantilever test was previously shown to accurately assess mechanical properties of paper-based materials, such as elastic and bending moduli as well as the bending length. This work investigates the hanging pear loop and clamped fold tests and compares the results with those found with the cantilever test. The results show that the strain and curvature induced by the two tests are much larger than those experienced in the cantilever test. This large strain and curvature induce plastic behaviour and make the hanging pear loop and clamped fold tests inappropriate for use on paper-based materials. [ABSTRACT FROM AUTHOR]

A computationally efficient method to study the in-plane and out-of-plane dimensional instability of thin paper sheets under the influence of moisture changes is presented. The method explicitly resolves the bonded and the free segments of fibers in the sheet, capturing the effect of anisotropic hygroexpansion at the fiber level. The method is verified against a volumetric model. The importance of longitudinal fiber hygroexpansion is demonstrated in spite of the absolute value of longitudinal hygroexpansion being an order of magnitude lower than the transverse hygroexpansion component. Finally, the method is used to demonstrate the formation of macroscopic sheet curl due to a moisture gradient in structurally uniform sheets in the absence of viscoelastic or plastic constitutive behavior and through-thickness residual stress profiles. [ABSTRACT FROM AUTHOR]

Cantilevers are one of the most utilized mechanical elements for acoustic sensing. In comparison to the edge clamped diaphragms of different shapes, a single edge clamped cantilever makes an acoustic sensor mechanically sensitive for detection of lower pressure. The aspect ratio of cantilevers is one of the most important parameters which affect sensitivity. Herein, we present a mathematical, finite element method and experimental analysis to determine the effect of the aspect ratio on the resonant frequency, response time, mechanical sensitivity, and capacitive sensitivity of a cantilever-based acoustic pressure sensor. Three cantilevers of different aspect ratios (0.67, 1, and 1.5) have been chosen for sound pressure application to detect capacitance change. The cantilever with the smallest aspect ratio (0.67) has the highest response time (206 ms), mechanical sensitivity, and capacitive sensitivity (22 fF), which reduce after increasing the aspect ratio. The resonant frequency of the cantilever was also analyzed by applying sweep in sound frequency. It was found to be minimum for the cantilever with the smallest aspect ratio (510 Hz) and increases with an increase in the aspect ratio. We have applied the garage fabrication process using low cost, recyclable, and easily available materials such as metal coated polymer sheets, mounting tapes and glass slides as alternative materials for expensive materials. [ABSTRACT FROM AUTHOR]

Oil-immersed paper is a weak link in internal insulation of oil immersed power transformer, whose aging state distributes inhomogeneously. However, most of the current assessment methods fail to consider the spatial distribution of aging areas inside transformer. Aiming at this problem, this paper proposes the inversion detection method to obtain the resistivity of oil-immersed paper in various regions nondestructively, and uses the resistivity of oil-immersed paper to assess directly the insulation state of oil immersed power transformers. Finite element method is applied to establish the mapping relationship between the resistivity of oil-immersed paper in different regions, and the dielectric loss factor obtained from the port of transformer. Back propagation neural network is used to learn inversely this mapping relationship so as to calculate the resistivity. Applying the proposed method on a 10-kV transformer, the calculated results are close to the measured values. The study shows that the inversion method can effectively calculate the oil-immersed resistivity in various regions, which represent the aging state in different regions of the transformer. [ABSTRACT FROM AUTHOR]

This article describes a new method for evaluation of radial crush strength of paper cores. Based on the theory of multilayered composites, a mathematical model has been developed to calculate the value of destructive force during radial crushing of spirally shaped paper cores whose layers are made of different materials. To start calculations, it is necessary to have geometry parameters of the cores, Young's modulus and SCT of an individual component layers. [ABSTRACT FROM AUTHOR]

This paper proposes a layer-by-layer acoustic travel time approximation method based on ray theory for total focusing method (TFM) imaging in carbon fiber reinforced polymer (CFRP) laminates. The method considers the anisotropy in every monolayer and heterogeneity of CFRP, which approximates the path of propagation as straight in the whole material. The application of this method to TFM imaging is called TravelTimeAppro-TFM. In comparison to isotropic-TFM and Dijkstra-TFM, the experimental results indicated that TravelTimeAppro-TFM outperforms isotropic-TFM in terms of imaging amplitude gain with a maximum gain of 4.67 dB. On the other hand, this approach reduces the computational work compared to Dijkstra-TFM. The proposed method demonstrates significant improvements in both focusing performance and the speed of calculation. This paper also investigates the effective angular range of the layer-by-layer acoustic travel time approximation method through experimental and finite element simulation analysis. [ABSTRACT FROM AUTHOR]

An epiretinal membrane (ERM) is a fibrocellular proliferation on the inner surface of the retina causing blurred and distorted central vision. Surgery is the only effective method for ERM removal. This paper investigated the mechanical properties of ERM peeling using the finite element (FE) method. A FE model of ERM formation on the retina surface was constructed. The failure criterion was applied to the attachment pegs to represent the adhesive force between the ERM and retina. The simulation results were consistent with the experimental data in published research. The maximum peeling force was 4.1 mN at a peeling velocity of 2 mm/s and an angle of 30°. The peeling force was minimum at the peeling angle of 45° and increased with the increase in peeling velocity and Young's modulus of the membrane. The outcome of this paper can improve the safety and efficiency of ERM removal. [ABSTRACT FROM AUTHOR]

Electroporation is a complex, iterative, and nonlinear phenomenon often studied through numerical simulations. In recent years, simulations of tissue electroporation have been conducted with static models. However, the results of a static model simulation are restricted to a fixed protocol signature of the pulsed electric field. In this paper, we describe a novel dynamic model of tissue electroporation that also accounts for tissue dispersion and temperature to allow time-domain simulations. We have implemented the biological dispersion of potato tubers and thermal analysis in a commercial finite-element method software. A cell electroporation model was adapted to account for the increase in tissue conductivity. The model yielded 12 parameters divided into three dynamic states of electroporation. The thermal analysis describes the dependence of tissue conductivity on temperature. The model parameters were evaluated using experiments with vegetal tissue (Solanum tuberosum) under electrochemotherapy protocols. The proposed model can accurately predict the conductivity of tissue under electroporation from 100 to 1000 V/cm. A negligible thermal effect was observed at 1000 V/cm, with a temperature increase of 0.89 ° C. We believe that the proposed model is suitable to describe the electroporation at the tissue level and provides a hint of the effects on the cell membrane. [ABSTRACT FROM AUTHOR]

• Numerical model for solution of heat diffusion equation. • Thermographic investigation in optically semi-transparent materials. • Non-destructive detection of graphic features buried in cultural heritage artefacts. • Selective configurations for optical and thermal properties investigation in paper. In this paper, infrared thermography is applied to the quantitative characterization of optically semi-transparent media by using a newly introduced mathematical model. The method is specifically aimed to the detection and evaluations in graphic features buried beneath the surface of some particular kind of cultural heritage artefacts, such as books and documents. A numerical approach based on the finite element method, is proposed for solving the model equations. The computed thermographic signal depends on different inspected material physical properties, which have been selectively evaluated by comparing the theoretical results with the experimental ones in specifically designed test samples. Thereafter, the model is applied to the analysis of the results obtained in the case of graphic elements buried within samples constituted by the same previously characterized materials. [ABSTRACT FROM AUTHOR]

The voltage on the valve side windings of the converter transformer is non-sinusoidal, which contains many high-order harmonic components. To calculate the nonlinear electric field distribution under the non-sinusoidal voltage, the harmonic-balanced decomposition finite element method (HBDFEM) is proposed in this paper. By using the fixed-point method, the finite element stiffness matrix of each iteration is decomposed according to the harmonic order, by which the memory consumption can be remarkably reduced. To verify the effectiveness of the proposed algorithm, a typical model of the oil-paper insulation structure is taken for instance, and the results obtained by the HBDFEM is compared with two algorithms, the traditional HBFEM and the time step method. The calculation results show that the HBDFEM can effectively reduce the computation cost for the nonlinear electric field problems with the non-sinusoidal periodic voltage. And it can be applied to analyze the compound electric field problems of the actual converter transformer. The electric field distribution of a ±500 kV converter transformer is calculated and analyzed considering the nonlinear relationship between the electric field and the medium conductivity. [ABSTRACT FROM AUTHOR]

In this paper, the multi-objective optimization problem (MOOP) of selecting suitable materials from four materials for phononic crystals layouts is investigated based on an optimization method combined nondominated sorting-based genetic algorithm II (NSGA-II) and finite element method (FEM). Driven by different optimization objectives, the changes in the optimized structure compared to the seed structure are the change of the bandgap mechanism from a local resonance mechanism to a Bragg scattering mechanism, the change of the included materials from four to three, and the significant change in the location and extent of the pores and other materials, respectively. The obtained nondominated Pareto solution of MOOP can balance the bandgap and the structural mass, which provides the decision-maker trade-off to select the appropriate optimized solution. Compared with the single-objective optimization problem (SOOP), the MOOP not only obtains a nondominated solution close to the result of SOOP, but also obtains other nondominated solutions of MOOP, which proves the effectiveness of the optimization algorithm in this paper. The design method in this paper can be easily extended to select suitable materials from a wider variety of materials for bandgap optimization design of PnCs, which has very promising applications. [ABSTRACT FROM AUTHOR]

This paper presents a numerical comparison of finite-element methods resulting in local mass conservation at the element level for Poisson's problem, namely the primal hybrid and mixed methods. These formulations result in an indefinite system. Alternative formulations yielding a positive-definite system are obtained after hybridizing each method. The choice of approximation spaces yields methods with enhanced accuracy for the pressure variable, and results in systems with identical size and structure after static condensation. A regular pressure precision mixed formulation is also considered based on the classical RT k space. The simulations are accelerated using Open multi-processing (OMP) and Thread-Building Blocks (TBB) multithreading paradigms alongside either a coloring strategy or atomic operations ensuring a thread-safe execution. An additional parallel strategy is developed using C++ threads, which is based on the producer-consumer paradigm, and uses locks and semaphores as synchronization primitives. Numerical tests show the optimal parallel strategy for these finite-element formulations, and the computational performance of the methods are compared in terms of simulation time and approximation errors. Additional results are developed during the process. Numerical solvers often fail to find an accurate solution to the highly indefinite systems arising from finite-element formulations, and this paper documents a matrix ordering strategy to stabilize the resolution. A procedure to enable static condensation based on the introduction of piecewise constant functions that also fulfills Neumann's compatibility condition, and yet computes an average pressure per element is presented. • Numerical comparison of the mixed, primal hybrid, and other locally conservatives finite element methods. • Introduction of a twice hybridized primal hybrid method, yielding a positive definite matrix after static condensation. • Comparison of multi-threaded strategies and thread-safety mechanisms. • A numerical procedure to stabilize the linear system resolution for indefinite matrices is described. • A general procedure that enables static condensation and fulfills Neumann's compatibility condition. [ABSTRACT FROM AUTHOR]

Research laboratories worldwide use the spark plasma sintering (SPS) technique for sintering metallic and ceramic matrix composites. Typically, the SPS sintered parts are low circular cylinders (discs). The upscaling and industrialization of SPS technology face two challenges: the limitation of shaping non-cylindrical parts and the inhomogeneous sintering of large-sized products. This paper addresses these challenges in sintering axisymmetric shells within an SPS setup with improved thermal insulation. The paper considers the sintering of a 75 mm shell as a case study. The shell material is a zirconia-carbon nanotube composite. Finite element analysis revealed that homogeneous sintering of such a shell requires careful thermal insulation of the entire setup. Investigation of phases, microstructure, density, microhardness, and fracture toughness in the sintered shell showed a homogeneous structure and uniform properties. [ABSTRACT FROM AUTHOR]

ABSTRACT In digital healthcare engineering (DHE) for aging monopile‐type offshore wind turbines, predictive health analysis is essential for robust future maintenance planning. This paper presents the development of a DHE module for the predictive health analysis of corroded monopile‐type offshore wind turbines, focusing on wind and rotor blade rotation effects. An empirical formula for predicting time‐variant corrosion wastage is derived from a statistical analysis of a decade‐long corrosion wastage database and applied to predict the corrosion depth of a 5 MW monopile‐type offshore wind turbine, serving as an illustrative example. Nonlinear finite element analyses using LS‐DYNA are performed on the corroded turbine tower under combined wind‐induced loads and rotor‐induced thrust forces. Two types of corrosion wastage, pitting, and uniform (general) corrosion, are considered. The health condition of the corroded tower is evaluated based on serviceability limit state and ultimate limit state criteria. The methodology developed in this paper will be integrated into the DHE system currently being developed by the authors for aging monopile‐type offshore wind turbines. [ABSTRACT FROM AUTHOR]

Purpose: The paper's goal is to examine and illustrate the useful uses of submodeling in finite element modeling for topology optimization and stress analysis. The goal of the study is to demonstrate how submodeling – more especially, a 1D approach – can reliably and effectively produce ideal solutions for challenging structural issues. The paper aims to demonstrate the usefulness of submodeling in obtaining converged solutions for stress analysis and optimized geometry for improved fatigue life by studying a cantilever beam case and using beam formulations. In order to guarantee the precision and dependability of the optimization process, the developed approach will also be validated through experimental testing, such as 3-point bending tests and 3D printing. Using 3D finite element models, the 1D submodeling approach is further validated in the final step, showing a strong correlation with experimental data for deflection calculations. Design/methodology/approach: The authors conducted a literature review to understand the existing research on submodeling and its practical applications in finite element modeling. They selected a cantilever beam case as a test subject to demonstrate stress analysis and topology optimization through submodeling. They developed a 1D submodeling approach to streamline the optimization process and ensure result validity. The authors utilized beam formulations to optimize and validate the outcomes of the submodeling approach. They 3D-printed the optimized models and subjected them to a 3-point bending test to confirm the accuracy of the developed approach. They employed 3D finite element models for submodeling to validate the 1D approach, focusing on specific finite elements for deflection calculations and analyzed the results to demonstrate a strong correlation between the theoretical models and experimental data, showcasing the effectiveness of the submodeling methodology in achieving optimal solutions efficiently and accurately. Findings: The findings of the paper are as follows: 1. The use of submodeling, specifically a 1D submodeling approach, proved to be effective in achieving optimal solutions more efficiently and accurately in finite element modeling. 2. The study conducted on a cantilever beam case demonstrated successful stress analysis and topology optimization through submodeling, resulting in optimized geometry for enhanced fatigue life. 3. Beam formulations were utilized to optimize and validate the outcomes of the submodeling approach, leading to the successful 3D printing and testing of the optimized models through a 3-point bending test. 4. Experimental results confirmed the accuracy and validity of the developed submodeling approach in streamlining the optimization process. 5. The use of 3D finite element models for submodeling further validated the 1D approach, with specific finite elements showing a strong correlation with experimental data in deflection calculations. Overall, the findings highlight the effectiveness of submodeling techniques in achieving optimal solutions and validating results in finite element modeling, stress analysis and optimization processes. Originality/value: The originality and value of the paper lie in its innovative approach to utilizing submodeling techniques in finite element modeling for structural analysis and optimization. By focusing on the reduction of finite element models and the creation of smaller, more manageable models through submodeling, the paper offers designers a more efficient and accurate way to achieve optimal solutions for complex problems. The study's use of a cantilever beam case to demonstrate stress analysis and topology optimization showcases the practical applications of submodeling in real-world scenarios. The development of a 1D submodeling approach, along with the utilization of beam formulations and 3D printing for experimental validation, adds a novel dimension to the research. Furthermore, the paper's integration of 1D and 3D submodeling techniques for deflection calculations and validation highlights the thoroughness and rigor of the study. The strong correlation between the finite element models and experimental data underscores the reliability and accuracy of the developed approach. Overall, the originality and value of this paper lie in its comprehensive exploration of submodeling techniques, its practical applications in structural analysis and optimization and its successful validation through experimental testing. [ABSTRACT FROM AUTHOR]

The influence evaluation of the additive fibre on concrete structures is a concerned problem especially in strengthening and repairing engineering. Current theoretical studies on the damage mechanisms of tunnel lining are generally based on the tensile and compressive characteristics of the ordinary concrete. Little attention is paid to the toughness and damage development of the fibre reinforced concrete, as well as the influences of high-speed train dynamic loading. In this paper, both tensile and compressive experiments of fibre reinforced concrete with different lengths and volume contents are carried out considering three types of the glass, PVA (polyvinyl alcohol) and hybrid fibres. Afterwards, the damage evolution of fibre reinforced concrete is explored by application of the digital image correlation (DIC) technique. The hypothesis of strain equivalence and damage theory of Najar are employed to determine the damage parameters of the fibre reinforced concrete so as to establish three-dimensional finite element model of surrounding rock-tunnel. Finally, the parameters that affect damage development of concrete tunnel lining incorporated fibre under train loading, including the vehicle speed, fibre type, fibre length, and volume content, are performed. The results show that adding fibre into concrete can effectively enhance the damage resistance capacity of concrete and significantly improve the developing trend of multistep damage of ordinary concrete under train excitation loading. Moreover, fibre reinforced concrete can reduce the final damage values of tunnel lining subjected train loading. The paper contributes to the understanding of damage behaviours of fibre reinforced concrete lining for high-speed trains loading as a potential risk in the course of the operation. [ABSTRACT FROM AUTHOR]

The deterministic methods are efficient for solving the neutron transport equation (NTE), but suffer from discretization errors. The increasingly complex geometric models make spatial discretization errors the primary source of discretization errors. Considering that spatial discretization errors come from inaccurate geometric shape descriptions and low-accuracy numerical schemes, this paper develops a Discontinuous Galerkin Finite Element Method for the NTE on unstructured polygonal meshes to reduce spatial discretization errors. In this method, the physical modal basis is adopted to handle the polygonal mesh and to achieve high-order accuracy in a uniform and efficient way. The numerical results of various fixed-source and k-eigenvalue benchmarks demonstrate that the method developed in this paper can give accurate solutions on polygonal meshes with high convergence rates. [ABSTRACT FROM AUTHOR]

• Nonlinear model for a bended piezoelectric semiconductor beam is established. • Differential quadrature method based iteration procedure is developed. • Nonlinearity of constitutive relation and non-uniformity of structure are focused. • Potential barriers and wells are realized through designing structures. In this paper, the effects of nonlinear constitutive relation on the electro-elastic coupling behaviors in non-uniform piezoelectric semiconductor beams are investigated. On the basis of three-dimensional theory for piezoelectric semiconductor and double power series, one-dimensional bending model is developed. In order to deal with the nonlinear problem, a differential quadrature method based iteration procedure is proposed. Comparing with the results from finite element method, the proposed method has good convergence and high accuracy. In numerical analysis, the effects of nonlinearity and non-uniformity in piezoelectric semiconductor are discussed thoroughly. It is found the introduction of nonlinearity induces the loss of symmetric for all electrical field quantities but little changes the mechanical quantities. While the reduction of width along the length direction produces larger mechanical and electrical responses. Importantly, the effects of nonlinearity on the electrical field quantities are more evident in a piezoelectric semiconductor beam with rapidly varied cross section. Besides, we designed the piezoelectric semiconductor beams subject to double shear forces and a beam with locally varied cross section. In the designed structures, the potential barriers and wells are realized, which provide new approaches adjusting the electro-elastic behaviors in piezoelectric semiconductors. The study in this paper could be the guidance designing high performance electronic devices. [ABSTRACT FROM AUTHOR]