Large Deformation Research Articles

Numerical Study on Large Deformation Characteristics of Tunnels Excavated in Strain-Softening Time-Dependent Rock Masses

Numerical Study on Large Deformation Characteristics of Tunnels Excavated in Strain-Softening Time-Dependent Rock Masses

Nonlinear analysis of three-dimensional hyperelastic problems using radial point interpolation method

Hyperelastic materials are primarily common in real life as well as in industry applications, and studying this kind of material is still an active research area. Naturally, the characteristic of hyperelastic material will be expressed when it undergoes large deformation, so the geometrical nonlinear effect should be considered. To analyze the behavior of hyperelastic material, the Neo-Hookean model is imposed in this study because of its simplicity. The model shows the nonlinear behavior when the deformation becomes large due to the nonlinear displacement-strain relation. This constitutive relation also gives a good correlation with experimental data. This study performs a meshless method, namely the radial point interpolation method (RPIM), to analyze the nonlinear behavior of Neo-Hookean hyperelastic material under a finite deformation state in three-dimensional space. The standard Newton-Raphson technique is applied to obtain the nonlinear solutions. Unlike mesh-based approaches, the meshless method shows its advantages in large deformation problems due to its mesh independence. In this paper, the numerical results of three-dimensional problems that undergo large deformation will be calculated and validated with solutions derived from previous studies. By investigating the obtained results, the superior ability of RPIM in hyperelastic problems can be proved.

Recyclable Chitosan-Modified Cellulose Fiber Porous Structure for Sensitive and Robust Moisture-Driven Actuators and Automatic Cooling Textiles.

Moisture-driven actuators featuring programmable stimuli-responsiveness and a rapid response have garnered substantial research attention. Cellulose-based actuators face challenges, including prolonged and unstable responsiveness, along with inadequate interfacial bonding. Herein, we developed a bilayer structured moisture actuator by integrating multiscale cellulose fibers with chitosan. The protonated chitosan forms strong electrostatic attractions with negatively charged cellulose nanofibrils (CNF), achieving a robust interfacial interaction. Leveraging the hierarchically porous structure and varying hygroscopicity of microfibrillated cellulose (MFC) and CNF, the film establishes an effective wettability gradient, enabling a stable and rapid moisture actuation performance. The bilayer film exhibits large deformation toward moisture with a bending angle of 60°, a short response time of 12 s, good stability over 50 wetting and drying cycles, and promising recyclability. Harnessing these advantageous properties, the bilayer film was demonstrated for its applications in automatic cooling textiles, contactless electrical switches, and artificial moisture-activated muscles, showing great potential for practical use.

Nonlinear seismic response analysis of long-span railway cable-stayed bridges crossing strike-slip faults

Active faults along railways in the mountainous regions of western China pose significant challenges to bridge safety. To ensure the safe operation of long-span railway bridges under complex geological conditions, this study investigates the synthesis of artificial ground motions for bridges crossing strike-slip faults and analyzes their nonlinear seismic response. First, we develop a theoretical method for simulating high- and low-frequency seismic motions using a finite fault and an equivalent velocity pulse model. Next, using a specific long-span railway cable-stayed bridge as a case study, we construct a nonlinear finite element model with OpenSees software. Finally, we assess the seismic response of key bridge components considering various crossing angles, seismic amplitudes, fault rupture directivity, and fling-step effects. The results show that the crossing angle significantly influences the seismic response, with longitudinal and transverse responses exhibiting opposite patterns. Additionally, the scaling factor of seismic motion significantly affects bridge response. For bridges crossing strike-slip faults, the longitudinal response exhibits a sudden increase in displacement due to instantaneous velocity pulses, while the transverse response shows notable residual displacement influenced by the fling-step effect. However, the critical section curvatures of bridge towers and piers remain within the elastic range across all crossing angles, indicating that controlling large displacement deformations is crucial for the seismic design of bridges crossing strike-slip faults.

Performance enhancement of magnetoactive elastomer soft robots through a coupled magnetoelastic topology optimization scheme

Abstract We have found that considering local magnetic fields, large deformations, and magnetoelastic coupling simultaneously significantly affect the resulting shape in magnetoelastic topology optimization in a uniaxial actuator case. In contrast to the work presented here, other works incorporate magnetoelastic formulations that include simplifying assumptions on the local field, and subsequent effects on the magnetization response of the material, or the absence of large deformation mechanics, or both. These assumptions were shown to produce solutions that differ substantively from cases where local fields and large deformations are addressed concurrently. Magnetoelastic topology optimization schemes are needed to optimize magnetoactive elastomer (MAE) devices. MAE devices are magnetic particle-filled polymer matrices designed for specific actuations and controlled remotely by an external magnetic field. They garner considerable research interest as an emerging technology for actuators in soft robots or in applications requiring untethered actuation. The material properties of MAEs are dependent on the volume fraction of particles in the elastomer matrix, where a high-volume fraction increases relative permeability (for soft magnetic particles) but also increases elastic modulus. For optimal actuation, a tradeoff between low stiffness and high magnetic response must be made by adjusting volume fraction and controlling material placement. Using a topology optimization scheme that considers both the magnetic and mechanical properties of the material, the shape and material composition of the device can be tuned to best achieve the desired actuation displacement. In this work, a density-based magnetoelastic multimaterial topology optimization scheme for soft magnetic material is developed in COMSOL Multiphysics. The optimization scheme uses multiphysics coupling that considers local magnetic fields and large deformations at each iteration through a Maxwell stress tensor formulation. A simulated example is then considered to demonstrate the effectiveness and necessity of a coupled optimization. The effect of considering large deformations during optimization is also investigated. It was found that a coupled topology optimization scheme with large deformations produced shapes with modes of actuation not captured by schemes with simplifying assumptions, leading to better performance at lower material cost. Considering large deformations in the coupled scheme offered significantly better performance, with an increase of 81.3% in a side-by-side performance simulation when compared to uncoupled cases.

Performing Magnetic Boundary Modulation to Broaden the Operational Wind Speed Range of a Piezoelectric Cantilever-Type Wind Energy Harvester

Small piezoelectric wind-induced vibration energy harvesting systems have been widely studied to provide long-term sustainable green energy for a large number of wireless sensor network nodes. Piezoelectric materials are commonly utilized as transducers because of their ability to produce high output power density and their simple structure, but they are prone to material fracture under large deformation conditions. This paper proposes a magnetic boundary modulated stepped beam wind energy harvesting system. On the one hand, the design incorporates a composite stepped beam with both high- and low-stiffness components, allowing for efficient vibration and electrical energy output at low wind speeds. On the other hand, a magnetic boundary constraint mechanism is constructed to prevent the piezoelectric sheet from breaking due to excessive deformation. Experiments have confirmed that the effective operational wind speed range of the harvester with magnetic boundary constraints is doubled compared to that of the harvester without magnetic boundary constraints. Furthermore, by adjusting the magnetic pole spacing of the boundary, the harvesting system can generate sufficiently high output power under high-wind-speed conditions without damaging the piezoelectric sheet.

Cyclic warm rolling: A path to superior properties in Mo[sbnd]Cu composites

Molybdenum–copper (MoCu) composites have a low coefficient of thermal expansion, good electrical and thermal conductivity and mechanical properties, and are widely used in microelectronic packaging heat dissipation materials, aerospace and other fields. Due to the large difference in the properties of Mo and Cu, the deformation of MoCu composites is difficult. At present, there is a lack of research on the deformation process and property changes of MoCu composites with large deformation. In this study, 74 % MoCu30 composites with large deformation are prepared by cyclic warm rolling, and the deformed materials have excellent mechanical and physical properties, and the evolution of the microstructure of the composites during the deformation process is described.

Quantitative Macromolecular Modeling Assay of Biopolymer-Based Hydrogels

The rubber elasticity theory has been lengthily applied to several polymeric hydrogel substances and upgraded from idealistic models to consider imperfections in the polymer network. The theory relies solely on hyperelastic material models in order to provide a description of the elastic polymer network. While this is also applicable to polymer gels, such hydrogels are rather characterized by their water content and visco-elastic mechanical properties. In this work, we applied rubber elasticity constitutive models through hyperelastic parameter identification of hydrogels based on their stress–strain response to compression. We further performed swelling experiments and determined the intrinsic properties, i.e., density, of the specimens and their components. Additionally, we estimated their equilibrium swelling and employed it in the swelling-equilibrium theory in order to determine the polymer–solvent interaction parameter of each hydrogel with regard to cross-linking. Our results show that the average mesh size obtained from the rubber elasticity theory can be regarded as a concentration-dependent characteristic length of the hydrogel’s network and couples the non-linear elastic response to the specimens’ inherent visco-elasticity through hysteresis as a quantifier of energy dissipation under large deformation.

Computational design of 4D printed shape morphing lattices undergoing large deformation

Abstract In 4D Printing, active materials are embedded in structures such that the application of an external stimulus, usually coming from the environment, results in a structural response. To design structures that achieve a targeted shape change for a defined stimulus, also known as shape morphing, the material distribution and structure needs to be tuned. However, the computational design of such material distributions and structures is a challenging task and remains, despite recent advances, unable to fully leverage the entire design freedom offered by state-of-the-art 4D printing technology. Notable gaps concern the handling of large and complex deformations, the high computational cost, and the exploration of the design space by the generation of alternative solutions. In this article, a method is presented to fill this gap. First, an artificial neural net is trained that represents a deformation map that occurs during actuation. Then, a shape morphing truss is designed that achieves this deformation during actuation. The method is used to solve four shape morphing problems, where superior capabilities are demonstrated in terms of magnitude and complexity of deformations that can be handled, efficient generation of alternative solutions and versatility. Due to these capabilities, the method enables exploration of the full potential of 4D printing technology to create stimuli-responsive, multifunctional structures.

Comprehensive review of recent research trends in dielectric elastomer transducers: applications, materials, and fabrication

Abstract While rigid transducers have led to innovations in industrial automation and control systems,
their rigidity limits them to the controlled environment of the factory. Recently, soft transducers
have become an attractive area of research because their compliant nature presents potential to
be applied to soft-bodied robots operating in uncontrolled environments. Dielectric elastomer
transducers (DETs) demonstrate large electrically-driven deformations, high energy density, fast
response, long lifespans, and self-healing capabilities which makes them a favorable replacement
for rigid transducers. This review first introduces the working principles of DETs, modeling work,
and different transducer designs. Then the applications of DETs to robotics, generators, and
electronics are reviewed. Lastly, the challenges currently faced by DET technology are discussed
and potential approaches are explored.

Modelling Method and Meshing Characteristics of a Novel Curve-Surface Conjugate Internal Gear Drive

Space exploration has become a major focus in the field of technology, with gear transmissions in aerospace equipment playing a crucial role. In the extreme environment of space, gear transmissions face challenges like large temperature differentials, deformation and maintenance difficulties, which will severely impact transmission accuracy and service life. To meet the growing demands for high-performance gear transmissions with high transmission efficiency and error adaptability in the aerospace field, this paper proposes a novel curve-surface conjugate internal gear drive consisting of an involute internal gear and a curve-surface gear. The fundamental theory of curve-surface conjugation is introduced, and the construction method for curve-surface gear based on a selected contact path and meshing tube is presented. The analysis models including induced curvature, sliding ratio and tooth contact analysis with errors (ETCA) are simulated to evaluate the meshing characteristics. Additionally, prototypes are manufactured and experimental setups are established to validate the transmission performance. These results indicate that as the rotational speed increases, the transmission efficiency of the curve-surface conjugate internal gear drive improves, which is contrary to the trend observed in involute gear drives. And the transmission efficiency of the curve-surface conjugate internal gear drive surpasses that of the involute gear drive at higher rotational speeds. Moreover, this novel gear drive exhibits excellent error adaptability, maintaining intact contact paths and high transmission efficiency even in the presence of assembly errors. This study provides new ideas for the design and manufacture of high-performance gear transmissions from the perspective of spatial geometric elements.

The influence of static versus dynamic pressure distribution strategies for modelling nonlinear waves generated by ships

A moving static pressure distribution is commonly used to simulate a travelling ship. However, the ship movement changes the fluid velocity around the hull, inducing pressures on the hull surface that are no longer equal to the static pressure. Therefore, we introduce a dynamic pressure correction strategy, which can accurately simulate the impact of the ship movement on the hull-surface pressure and preserve the desired hull shape under both stationary and transient conditions. The strategy is applied to a high-order spectral model and used to investigate ship-induced waves and wave resistance over a both flat and variable topography. We explore various parameters in our study, including the average water depth to ship draft ratio ( $h_0/d$ ), the channel width to ship width ratio ( $W/B$ ), the Froude number ( $Fr_0=U/\sqrt {gh_0}$ ) and variations in bathymetric slope. Compared with experiments on a flat bottom, the numerical results with dynamic correction show better accuracy in the simulation of ship-induced waves and wave resistance than those obtained using a static pressure distribution. The correlation coefficient for wake waves between the numerical and experimental results is improved by approximately 0.25 with the dynamic correction strategy. The amplitude and wavelength of ship-induced mini-tsunamis over a variable topography are found to be reduced when employing a dynamic correction compared with a static pressure distribution, and this effect becomes more pronounced with higher Froude number. The static pressure approach is shown to allow large deformations of the desired hull shape and changes in ship volume which are responsible for the different wave patterns from the two approaches.

Stacking fault energy during DRV-DRX competition in high-nitrogen austenitic stainless steel used in orthopedic implants

High-nitrogen (high-N) austenitic stainless steels (ASS) are used in orthopedic implants due to their mechanical properties, human biocompatibility, and affordable cost. However, the high content of (Ni–Nb–N)-rich precipitates can cause allergic reactions and significant challenges in the manufacturing of prostheses. This study investigates the competition between work hardening (WH), dynamic recovery (DRV), and dynamic recrystallization (DRX) during the physical simulation of an ASTM F-1586 steel by thermomechanical processing, using the Kocks-Mecking and Avrami constitutive models. The parameters were obtained through continuous isothermal torsion tests at the temperature range of 900–1200 °C, strain rates between 0.01 and 10 s−1 and a total deformation of 4.0. Determined by compositional-analytical methods, the stacking fault energy (γsfe) was correlated with the stress level (σi) according to the Uesugi and Dai models. Results indicated a high activation energy for hot deformation (Qdef = 587 kJ/mol), affecting the shape of the curves. The γsfe value varied along the curves, delaying the onset of the DRX (σc = 0.928σp). Scanning electron microscopy (SEM) and transmission electron microscopy (TEM) analyses revealed a direct competition between work-hardened and recrystallized grains in the early part of the curves, due to WH-DRV synergy. After the peak stress (σp), the progress of DRX was slow, with the Avrami exponent (n) between 1.1 and 1.9, completing only after large deformations (σss = 0.714σp). The moderate γsfe value (69 mJ/m2) and fine precipitates of the Z-phase (CrNbN) (<20 nm) influence the grain boundary mobility during the DRV-DRX competition, delineating the stress-strain curve shape.

Casein network formation at oil–water interfaces is reduced by [formula omitted]-casein and increased by Ca[formula omitted]

Mixtures of bovine caseins can serve as a benchmark for understanding the functionality of microbial-based recombinant caseins at oil–water interfaces. In this work we show that, in the presence of Ca2+, the individual casein fractions form viscoelastic networks at the oil–water interface with comparable stiffness. In the absence of Ca2+, αs2- and β-casein interfacial network formation was strongly inhibited over the full deformation regime. For αs1-casein, the network stiffness was increased in the absence of Ca2+ at small deformations (<15%), but at large deformations (>50%) it was completely disrupted, to a similar stiffness as αs2- and β-casein. The interfacial structure formed by κ-casein was largely unaffected by Ca2+ due to limited phosphorylation. We hypothesize that the differences between calcium-sensitive caseins lie in the conformation they assume at the interface. Both αs2- and β-casein adsorb in a train-tail conformation with a tail extending into the aqueous bulk phase, whereas αs1-casein adsorbs in a loop-train conformation, with a loop that extends less into the bulk phase. The tail-train configuration is hypothesized to increase the inter-molecular Ca2+ bridging thereby increasing the interfacial stiffness of αs2- and β-casein.Blending the casein fractions revealed a strong negative effect of β-casein on the interfacial modulus, which was more pronounced at a higher concentration. The presence of Ca2+ remained important for interfacial network formation of a casein blend. Without Ca2+, the interfacial network was less stiff, more viscous, and behaved like a 2d polymer solution.With this work we showed that casein interfacial network formation at oil–water interfaces is mediated by Ca2+ bridging. Blending the different casein fractions decreased the interfacial viscoelastic properties through the presence of β-casein. These results indicate that future work on recombinant caseins should focus on single genetic variants, since a blend of variants will likely decrease interfacial functionality.

High-Conductivity and Ultrastretchable Self-Healing Hydrogels for Flexible Zinc-Ion Batteries.

Aqueous zinc-ion batteries are promising candidates for flexible energy storage devices due to their safety, economic efficiency, and environmental friendliness. However, the uncontrollable dendrite growth and side reactions at the zinc anode hinder their commercial application. Herein, we designed and synthesized a dual network self-healing hydrogel electrolyte with zwitterionic groups (PAM-PAAS-QCS), which can be used for the large deformations of flexible devices due to its excellent stretchability (ε = 5100%). The incorporation of zwitterionic groups into the PAM-PAAS-QCS hydrogel electrolyte endows it with high ionic conductivity (33.61 mS/cm), a wide electrochemical stability window, and the ability to suppress zinc dendrite formation and side reactions. Besides, the Zn//Zn symmetric cell with PAM-PAAS-QCS can stably plate and strip zinc for 1500 h at 0.5 mA/cm2, and the Zn//Polyaniline full cell retains 82.4% of its capacity after 1500 cycles at 1 A/g. Additionally, flexible batteries based on both the original and self-healed PAM-PAAS-QCS hydrogel electrolytes demonstrate good cycling stability and stable charge-discharge performance under various bending conditions. This self-healing hydrogel electrolyte with excellent stretchability and high ionic conductivity is expected to pave the way for the development of high-performance flexible energy storage and wearable devices.

Simulation analysis of different types of balloon dilatation catheters for the treatment of intracranial arterial stenosis.

The application of balloon dilation catheters in the management of intracranial arterial stenosis has been gradually increasing. However, studies on the feasibility and effectiveness of different types of balloons remain relatively scarce. In this study, catheter models of three different materials were designed to simulate balloon crimping,splitting, and dilatation processes. A compliant balloon produces large deformations with poor dilatation and a stress concentration phenomenon. During dilatation, the shear stress generated in the intima and lesion area by the semi-compliant balloon was smaller than that generated by the non-compliant balloon. These results demonstrate the feasibility of using semi-compatible balloons.

Finite element analysis of large deformation for progressive failure of soil based on Cosserat continuum theory

The large deformation accompanied by strain softening and strain localization of the material is usually encountered in geotechnical engineering. Integrating the RITSS method and the Cosserat continuum theory, a large deformation Cosserat-RITSS finite element method is developed to tackle this issue. Based on the second development interfaces of UEL and SDVINI offered by the finite element software ABAQUS, the Cosserat-RITSS method adopting the MUEM interpolation and mapping method is realized automatically and used to simulate the problem of progressive failure and large deformation characterized by strain localization. The accuracy and effectiveness of the proposed method are verified by the comparison of the small deformation and the large deformation analysis of the penetration process of a rigid strip foundation. Further, the Cosserat-RITSS method is used to analyze the large deformation of a retaining wall resisting on the soil with non-associated plastic flow and strain-softening properties. Compared with the classical RITSS method, it is demonstrated that the proposed method has better performancy of simulating large deformation and overcoming the pathological mesh-dependent problem usually suffered in the classical FEA of large deformation accompanied by strain localization.

MSIreg: an R package for unsupervised coregistration of mass spectrometry and H&E images.

Joint analysis of mass spectrometry images (MS images) and microscopy images of hematoxylin and eosin (H&E) stained tissues assists pathologists in characterizing the morphological structure of the tissues, and in performing diagnosis. Unfortunately, the analysis is undermined by substantial differences between these modalities in terms of aspect ratios, spatial resolution, number of channels in each image, as well as by large global or small local elastic spatial deformations of one image with respect to the other. Therefore, accurate coregistration of the images is a critical pre-requisite for their joint interpretation. We introduce MSIreg, an open-source R package for coregistration of MSI and H&E images. MSIreg is designed for high-dimensional MSI experiments where each spatial location is represented by thousands of mass features. Unlike most existing coregistration methods, MSIreg implements a landmark free workflow, and quantitative metrics for performance evaluation. We evaluate the performance of MSIreg on six case studies, including coregistration of contiguous tissues with large deformations, as well as simultaneous coregistration of 29 tissue microarray cores. The R package, installation instructions, and fully reproducible vignettes describing methods and Case Studies are are available open-source under the GPL-3.0 license at https://github.com/sslakkimsetty/msireg/.

Galerkin’s Method and Multivariate Newton’s Method for the Nonlinear Deformation of the Magneto-Electro-Elastic Bi-Layered Laminates

In this paper, mathematical modelling for the large deformation of a magneto-electro-elastic rectangular bi-layered laminate with general boundary conditions is presented. Constitutive equations involving the magneto-electro-elastic (MEE) material properties are introduced, Maxwell equations accounts for the electric and magnetic effects are also utilized. First-order shear deformation theory (FSDT) considering the von Karman nonlinear strain is adopted, and the plain strain/stress assumption applicable for thin plate analysis is used. A rather compact set of governing equations related to kinematical variables, electric/magnetic potentials and the Airy stress function is obtained as a consequence of the preliminary condensation for the electro-magnetic state to the plate kinematics. Semi-analytic solution for a bi-layered BaTiO3-CoFe2O4 laminate with specified boundary conditions subjected to various external applied loads is performed. By employing the Bubnov-Galerkin method, the set of nonlinear partial differential equations is transformed to a set of third-order nonlinear algebraic equations for the static deformation due to applied load. Numerical results are carried out by using the multivariate Newton's method with respect to various volume fractions indicating the volume ratio between piezoelectric BaTiO3 layer and piezomagnetic CoFe2O4. From the result, the nonlinearity of the von Karman strain appears to enhance system rigidity as smaller deformations will be detected when external load is applied. Also, some other interesting results are obtained which could be useful to future analysis and design of multiphase composite plates.

Failure Mechanism and Control Technology of the Primary Support for Asymmetric Deformation Disaster of Deep Tunnel

Asymmetric large deformation disasters inevitably occur when highway tunnels cross active fault zones. Relying on the Nanlangshan No.1 Tunnel Project, the mechanism of asymmetric deformation is analyzed based on the in-situ engineering geology and tunnel monitoring data, and numerical calculations are used for verification. Based on the results of asymmetric large deformation disaster analysis, the control scheme of double-layer primary support is proposed, and the design is optimized with numerical calculation before engineering practice. The conclusions are as follows: 1. The original primary support program cannot withstand the asymmetric surrounding rock pressure, and the slip of the weak interlayer in the engineering geology is the main reason for the asymmetric deformation disaster in the tunnel. 2. According to the results of the numerical calculations, the tilted weak interlayer greatly reduces the stability of the surrounding rock, which leads to the local bending damage in the left shoulder of the tunnel, and the local shear deformation damage in the right arch girdle. 3. When using the double-layer primary support program, the appropriate delay of support timing for the second layer of primary support can effectively decrease the tunnel asymmetric deformation. However, if the second layer of primary support lag distance is too large, it will also lead to tunnel deformation convergence larger, or support failure and other issues. So, it is recommended that the second layer of primary support lags behind the first layer of primary support by 2 m to be applied.