Inelastic Model Research Articles

Multi-physical modeling and numerical simulation of the thermo-hygro-mechanical treatment of wood

The contribution at hand introduces computational modeling and realistic simulation concepts for a comprehensive description of the manufacturing and application of densified wood and wooden structures made from molded densified wood. Wood, as a natural material, is characterized by e.g. a very good mechanical load-bearing capacity related to its density. Nevertheless, the ratio between its mechanical properties and its density can be optimized by densification technology for an expanded use of wood in structural engineering. The wood densification process is not only a mechanical process with large and irreversible deformations, it is also denoted by temperature- and moisture-dependent treatments of the wooden specimens. Thus, the introduced approaches to predict the material and structural characteristics of compressed and molded wood consist of an inelastic and multi-physical constitutive modeling of wood at finite deformations as well as the computation of effective structural properties of wood after the thermo-hygro-mechanical densification process. A successful implementation of the modeling concepts into the finite element method (FEM) is presented, which is verified by numerical investigations. A validation of the numerical results is carried out by use of experimental data at beech wood (Fagus Sylvatica, L.), taken from literature.

Ab initio determination of interfacial thermal conductance for polymer-bonded explosive interfaces

Understanding the thermal transport in polymer-bonded explosives (PBXs) is critical for enhancing the safety and reliability during PBX design, especially in the absence of effective experimental measurements. In this work, we rigorously investigated the phonon properties of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) and polyvinylidene fluoride (PVDF) and calculated the interfacial thermal conductance using an ab initio approach. The diffuse mismatch model and anharmonic inelastic model were adopted to examine the interfacial thermal conductance as a function of temperature for the TATB–PVDF interface. Our calculation results indicate that low-frequency phonon modes and the two-phonon process play dominant roles in the thermal transport at interfaces. In contrast, high-order phonon processes involving three to eight phonons accounted for around 8% of the interfacial thermal conductance at the TATB–PVDF interface. Phonon properties, such as the velocity and degree of phonon density overlap, are discussed for the TATB–PVDF and 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX)–PVDF interfaces to estimate the interfacial thermal conductance of PBXs. This study provides a theoretical explanation for the establishment of a research method for PBX thermal transport.

One-dimensional compression of a saturated elastoviscoplastic medium

We explore a model for compaction of a two-phase medium in which the solid is allowed to deform elasto-viscoplastically. Previous models for compression of colloidal and cellulose fiber suspensions often prescribe the solid stress as a material function of the local solid concentration. Our current model distinguishes between purely elastic or plastic deformation and allows a solid viscous response. This is crucial for model comparison to experiment where a dramatic improvement over inelastic constitutive models is demonstrated.

SuSAv2 model for inelastic neutrino-nucleus scattering

The susperscaling model SuSAv2, already available for charged-current neutrino-nucleus cross sections in the quasielastic region, is extended to the full inelastic regime. In the model the resonance production and deep inelastic reactions are described through the extension to the neutrino sector of the SuSAv2 inelastic model developed for ($e$, ${e}^{\ensuremath{'}}$) reactions, which combines phenomenological structure functions with a nuclear scaling function. This work also compares two different descriptions of the $\mathrm{\ensuremath{\Delta}}$ resonance region, one based on a global scaling function for the full inelastic spectrum and the other on a semiphenomenological $\mathrm{\ensuremath{\Delta}}$ scaling function extracted from ($e$, ${e}^{\ensuremath{'}}$) data for this specific region and updated with respect to previous work. The results of the model are tested against ($e$, ${e}^{\ensuremath{'}}$) data on $^{12}\mathrm{C}$, $^{16}\mathrm{O}$, $^{40}\mathrm{Ca}$, and $^{40}\mathrm{Ar}$ and applied to the study of the charged current inclusive neutrino cross-section on $^{12}\mathrm{C}$ and $^{40}\mathrm{Ar}$ measured by the T2K, MicroBooNE, ArgoNEUT, and MINERvA experiments, thus covering several kinematical regions.

Effect of hull inelasticity on whipping responses by underwater explosions

ABSTRACT This paper aims to verify why the inelastic hull responses of buckling and plastic deformation should be considered in the underwater explosion (UNDEX) simulations. A battleship was discretized into the five longitudinal segments that were composed of Timoshenko beam elements. The number of elements was determined through convergence simulations of natural frequency. The inelastic nonlinear moment-curvature curves using Smith’s method were applied to the UNDEX models. The incompressible gas bubble model was chosen for the UNDEX simulations because the compressible one produced acceleration less conservatively than the incompressible one. The UNDEX simulations were performed for the far and less far stand-off distances. The bending moment in the inelastic model was less developed than those in the elastic model for two stand-off distances. More evident Inelastic effects were observed for the reduced stand-off distance. The inelastic effect should be considered in the UNDEX analyses for even far-field UNDEX cases.

Detection prospects for multi-GeV neutrinos from collisionally heated GRBs

Neutrinos with energies ranging from GeV to sub-TeV are expected to be produced in Gamma-Ray Bursts (GRBs) as a result of the dissipation of the jet kinetic energy through nuclear collisions occurring around or below the photosphere, where the jet is still optically thick to high-energy radiation. So far, the neutrino emission from the inelastic collisional model in GRBs has been poorly investigated from the experimental point of view. In the present work, we discuss prospects for identifying neutrinos produced in such collisionally heated GRBs with the large volume neutrino telescopes KM3NeT and IceCube, including their low-energy extensions, KM3NeT/ORCA and DeepCore, respectively. To this aim, we evaluate the detection sensitivity for neutrinos from both individual and stacked GRBs, exploring bulk Lorentz factor values ranging from 100 to 600. As a result of our analysis, individual searches appear feasible only for extreme sources, characterized by gamma-ray fluence values at the level of F$_{\gamma} \geq 10^{-2}$ erg cm$^{-2}$. In turn, it is possible to detect a significant flux of neutrinos from a stacking sample of ~ 900 long GRBs (that could be detected by current gamma-ray satellites in about five years) already with DeepCore and KM3NeT/ORCA. The detection sensitivity increases with the inclusion of data from the high-energy telescopes, IceCube and KM3NeT/ARCA, respectively.

Heavy long-lived coannihilation partner from inelastic Dark Matter model and its signatures at the LHC

The coannihilation mechanism is a well-motivated alternative to the simple thermal freeze-out mechanism, where the dark matter relic density can be obtained through the coannihilation with a partner particle of similar mass with dark matter. When the partner particle is neutral, the inelastic nature of dark matter can help it to escape the direct detection limits. In this work, we focus on the coannihilation scenario in which the annihilation cross section is dominated by the partner-partner pair annihilation. We pay special interest on the parameter space where the coannihilation partner is long-lived, which leads to displaced signatures at the collider. In such case, it opens the heavy mass parameter space for the coannihilation dark matter, comparing with those dominated by the partner-dark matter annihilation. Specifically, we study an inelastic scalar dark matter model with a specific parameter space, which realizes the domination of partner-partner pair annihilation. Then, we study two different realizations of the coannihilation partner decay and the existing constraints from the relic abundance, direct and indirect dark matter detection and the collider searches. We focus on the channel that the long-lived coannihilation partner decays to dark matter plus leptons. The high-luminosity LHC can reach good sensitivities for such heavy dark matter and coannihilation partner around 100–700 GeV.

Force of droplet impact on thin liquid films

Inertia-dominated droplet impact transfers momentum to a dry flat target within a short span of time t characterized by (droplet diameter D)/(impact speed U). We investigate experimentally how impact force dynamics change when a droplet hits a thin liquid film of thickness H, less than or approximately equal to the droplet diameter, atop the flat target. Impact force and morphology are recorded simultaneously by piezoelectric force sensing and high-speed video imaging. Compared with a dry surface, the force of droplet impact on a thin liquid film is found to follow the same initial tU/D scaling and reach a slightly higher peak value, but at a significantly later time. Modeling the impact process as a perfect inelastic collision between the droplet and a liquid column of height equal to the film thickness yields the proper timescale (H+D)/U to characterize temporal evolution of the impact force near the inertial peak and through its subsequent exponential decay. The impact crater penetration depth developing within the thin film over the same time span is also found to collapse to a self-similar form based on this characteristic timescale, which attests to the validity of the inelastic collision model in capturing the underlying impact flow physics.

Performance-based-plastic-design of damage-control steel MRFs equipped with self-centring energy dissipation bays

Driven by the demand of producing novel seismic resilient structures, this study focused on steel moment resisting frames equipped with superelastic shape memory alloy (SMA) connections in the self-centring energy dissipation bays, and the emphasis was given to a direct-iterative design method for performance-based-plastic-design of the system under earthquake motions. Utilising the equivalent single-degree-of-freedom (SDF) oscillator which may reasonably quantify the dynamic responses of a multi-storey structure, an inelastic structural seismic demand model motivated by nonlinear spectral analyses was developed, and the statistic features were accounted. As the theoretical basis of the design philosophy, the seismic energy equilibrium equation of the structure was correlated with structural arrangement and earthquake motion characteristics. Then, a stepwise plastic design methodology enabling practising engineers to search for a feasible design strategy was proposed. Three prototype structures were developed following the proposed method, and finite element models of the systems were developed. The adequacy of the modelling techniques was verified by the test database, and the sufficiency of the hysteretic model of the structure was also justified. The seismic responses of the prototype structures were estimated by nonlinear static and dynamic analyses. The analysis data pool confirmed that the prototype structures were able to achieve the damage-control objective by limiting the maximum interstorey drifts below the prescribed deformation threshold, and inelastic actions were concentrated in the self-centring energy dissipation bays.

Investigation of radiation damage using thermal spike model for SHI irradiation on Al2O3

The material exposed to energetic radiation may change its properties and performance. Ion beam irradiation is an established method to investigate radiation damage in materials. In the present study, SRIM code is employed as a tool to determine the electronic and nuclear energy losses and displacement per atom (DPA) associated with the damage induced in Al2O3 using 100 MeV Au ion. We have further studied the creation of 100 MeV Au ion induced molten zone along the ion track in Al2O3 in the framework of the thermal spike (TS) model. Using the in-elastic thermal spike model, elastic thermal spike model, and unified thermal spike model, the evolution of lattice temperature and latent track radius as a function of depth of Al2O3 for 100 MeV Au ion irradiation are investigated. The current results demonstrate that the latent track radius shows a U-shaped variation based on the combined effect of electronic stopping power and nuclear stopping power. In order to corroborate these results, an additional irradiation experiment was performed using swift heavy ion (SHI), 100 MeV Au ions with fluence 1×1014 ions/cm2 for Al2O3. To investigate the structural properties, surface morphology and to identify the defects in the samples, pristine and irradiated Al2O3 samples were characterized by Grazing Incident X-Ray Diffraction (GIXRD), Atomic Force Microscopy (AFM) and Photoluminescence Spectroscopy (PL). AFM study shows the increase in roughness after SHI irradiation. A decrease in intensity of PL spectra upon SHI irradiation indicates the formation of new defects and disordered structures in the material, which is quite in agreement with the theoretical results.

A new proposal for design of laterally unsupported monosymmetric steel I-section beams at elevated temperature

Fire causes one of the worst hazardous situations to structures. With the increase in temperature, there is a rapid decrease in steel material properties, which eventually leads to a reduction in the flexural capacity of the beam members. Lateral-torsional buckling (LTB) may predominate the design strength predictions for laterally unsupported beams. A rapid decrease in the capacity at high-temperature may result in premature failure of beams. Past studies have focused on predicting LTB design curves using a single equation, varying with reduction factors at elevated temperatures. These curves may result in a very conservative design that might change with the loading patterns and can sometimes predict unsafe design strengths. The present manuscript carries out a parametric study using a validated numerical model. The robust computational software ABAQUS is used for simulating a total of 2100 elastic and inelastic beam models. A new proposal is recommended for predicting the design capacity of monosymmetric I-section beams for the uniform moment, moment gradients, point load at mid-span, and uniformly distributed load (UDL) at elevated temperature. Further, a statistical study is conducted to ascertain the accuracy of the present proposal. The resulted data is also compared with the available literature. The results show that the current proposal provides economical, and safe predictions.

A spatial and temporal characterisation of single proton acoustic waves in proton beam cancer therapy.

An in vivo range verification technology for proton beam cancer therapy, preferably in real-time and with submillimeter resolution, is desired to reduce the present uncertainty in dose localization. Acoustical imaging technologies exploiting possible local interactions between protons and microbubbles or nanodroplets might be an interesting option. Unfortunately, a theoretical model capable of characterising the acoustical field generated by an individual proton on nanometer and micrometer scales is still missing. In this work, such a model is presented. The proton acoustic field is generated by the adiabatic expansion of a region that is locally heated by a passing proton. To model the proton heat deposition, secondary electron production due to protons has been quantified using a semi-empirical model based on Rutherford's scattering theory, which reproduces experimentally obtained electronic stopping power values for protons in water within 10% over the full energy range. The electrons transfer energy into heat via electron-phonon coupling to atoms along the proton track. The resulting temperature increase is calculated using an inelastic thermal spike model. Heat deposition can be regarded as instantaneous, thus, stress confinement is ensured and acoustical initial conditions are set. The resulting thermoacoustic field in the nanometer and micrometer range from the single proton track is computed by solving the thermoacoustic wave equation using k-space Green's functions, yielding the characteristic amplitudes and frequencies present in the acoustic signal generated by a single proton in an aqueous medium. Wavefield expansion and asymptotic approximations are used to extend the spatial and temporal ranges of the proton acoustic field.

Performance analysis of nonlinear vibration energy harvesting system with inelastic barrier under colored noise excitation

The vibro-impact structures are introduced into the piezoelectric vibration energy harvesting (VEH) systems, which can broaden the operating bandwidth of the piezoelectric VEH systems and thus improve the harvesting efficiency. The process of vibro-impact not only causes energy loss during the impact process but also causes the local deformation between two impact bodies. Therefore, this manuscript proposes a class of stochastic non-classical inelastic impact VEH models, whose non-classical impact processes are described by Hertz contact forces. Firstly, the relationship between the total energy of the unperturbed system and the potential energy of the system at the barrier position is used to determine whether the impact phenomena have occurred or not, then used to derive the period and frequency of the system. Secondly, using the quasi-conservative stochastic averaging method, the averaged drift term and the averaged diffusion term of the piecewise smooth averaged Ito^ equation are derived. The stationary probability density functions (PDFs), the mean square output (MSO) voltage and the root mean square (RMS) voltage of the system response can be obtained from solving the corresponding Fokker-Planck-Kolmogorov (FPK) equation for averaged Ito^ equation. Finally, an example is provided to verify the effectiveness of the proposed method, and the effects of the model parameters on the MSO voltage and power conversion efficiency (PCE) are analyzed. The relationship between the MSO voltage and the PCE is also analyzed. It is of great significance for the design and optimization of stochastic VEH systems.

Response spectrum of a reinforced concrete frame structure under various column removal scenarios

In a frame structure, sudden removal of a load-bearing member (e.g., a column) causes structural vibrations, often leading to a partial or total progressive collapse. Since the term ”sudden column removal” is imprecise, it is worth evaluating how the rate of column removal influences the maximum displacement and how this response refers to the one produced by the static column removal. Both linear and nonlinear, inelastic modeling of the reinforced concrete members are considered using the OpenSees finite element framework. The numerical simulations consist of a gradual increase in column removal time from near zero to asymptotically approaching the static response. Due to the asymmetry of the analyzed frame building, three column removal scenarios are investigated. It is observed that the impact of column removal rate depends on the ability of the remaining structure to withstand the lack of column and the natural period of the mode shape corresponding to the downward motion of the model without the column. Also, the column removal durations producing 95% of the maximum response are identified and their ratios to the respective natural periods are given. These investigations can be helpful when planning numerical simulations and physical experiments. Finally, a novel simplified single degree of freedom system with support removal is proposed and its response spectrum is developed. This simplified model helps to understand more complex multi-degree of freedom systems.

Comparison between recent implicit time integration methods with frequency dissipation for nonlinear structural applications

The present paper aims to test recent (Truly self-starting two sub-step method and three-parameter single-step implicit method) and classical (Generalized-α, HHT-α, and WBZ-α methods) time integration methods using the geometrically nonlinear Positional Finite Element Method (PFEM). The numerical formulation is based on the total Lagrangian approach and uses the Hessian matrix to obtain the response. The mixed hardening inelastic model applied to PFEM is also presented. Two examples validate the time integration algorithms and the inelastic model. In the first example, the mixed hardening inelastic model is compared with the the bilinear stress-strain model and the elastic-perfectly plastic hinge model, and aspects such as amplitude decay and period elongation are discussed. In the second example, the implemented algorithms are verified in a severe geometrically nonlinear example, considering the influence of numerical dissipation, time interval, and the number of elements in the response. Results show the relevance of numerical damping for numerical stabilization and the good performance of the Generalized-α algorithm.

An Inelastic Material Model for Lateral Torsional Buckling and Biaxial Bending of Steel W-Shapes

An Inelastic Material Model for Lateral Torsional Buckling and Biaxial Bending of Steel W-Shapes

Electron–phonon coupling factor and electron heat capacity of 6H-SiC

ABSTRACT For high-energy single ion events, the relationship between electronic energy loss mechanism and atomic processes is generally described by the inelastic thermal spike model. However, the parameters required for the model are not accurately known and are often estimated from a free electron gas model. To ensure this model is more reliable and predictive, a more accurate calculation for the parameters should be taken. In this paper, the temperature dependence of the electronic heat capacity and electron–phonon coupling factor (e-ph) has been calculated, using density functional theory (DFT). The calculation results demonstrate that the effect of electron thermal excitation on the thermodynamic parameters is sensitive to the bandgap, and these parameters vary immensely with the electron temperature. Besides, to verify the accuracy of obtained parameters, both the swift heavy ions (SHI) irradiation in crystalline 6H-SiC and amorphous 6H-SiC are simulated, using the inelastic thermal spike model. And the results are more consistent with the experiment conclusion, for the obtained parameters by DFT. Hence, the calculated temperature-dependent parameters are more suitable for the inelastic thermal spike model.

Profile extrusion die design: A comparative study between elastic and inelastic fluids

AbstractComputational modeling is widely used to support the design of profile extrusion dies. However, despite the viscoelastic nature of the polymer melts, the majority of the computational approaches resort to inelastic models. With the aim of assessing the accuracy of the approach usually employed on the modeling of profile extrusion dies, this work aims at comparing the behavior of profile extrusion dies when interrelated viscoelastic (elastic) and generalized Newtonian (inelastic) constitutive models are used. For this purpose, the polymer melt employed in the study was experimentally characterized, being the data collected used to fit a nonlinear multimode viscoelastic (elastic) Giesekus constitutive model. Subsequently, the fitted model was used to generate the material shear flow curve, which was then used to fit a Bird‐Carreau (inelastic) constitutive model. The numerical studies undertaken were focused on the extrusion die of two profiles: a simple one, with a rectangular cross‐section, and a complex swimming pool cover profile. The results obtained showed that in realistic case studies, the effect of elasticity might be relevant, even when modeling just the flow in the extrusion die flow channel and, thus, should be considered when designing profile extrusion dies.

Quantifying the non-Newtonian effects of pulsatile hemodynamics in tubes

• A new hemorheological constitutive TEVP model able to reproduce experimental data. • Pulsating hemodynamic conditions are studied. • At high pressure-gradient, RBC aggregation primarily affects hysteresis phenomena. • The streamwise normal stress is the dominant stress component. • Inelastic models underestimate the kinematic phase-lag compared to the TEVP model. We investigate in silico the pulsatile blood flow in straight, rigid tubes using an elasto-viscoplastic, constitutive model coupled with thixotropy (TEVP) (Varchanis et al., 2019; Giannokostas et al., 2020). In addition to blood viscoelasticity, our model accounts for RBC aggregation through a kinetic equation that describes the level of blood structure at any instance and the viscoplasticity at stasis conditions. We evaluate the model parameters using steady, simple shear, and transient multi-shear rheometric experiments for human physiological subjects with normal blood aggregation. Then, we accurately reproduce previous experimental results (Thurston, 1975; Bugliarello and Sevilla, 1970; Thurston, 1976; Womersley, 1955) and provide reliable predictions for the pressure, stress, and velocity fields. We also investigate blood flow under sinusoidal and experimentally determined waveforms of the pressure-gradient with different frequencies, amplitudes, and patterns, providing a thorough parametric study. We find that the streamwise normal stress is of considerable magnitude. This stress component stretches the RBCs and their aggregates in the flow direction. Commonly used inelastic hemorheological models (e.g., Casson and Newtonian) cannot predict this. Typical values of the time-averaged wall normal-stress are found in the range of 2 Pa − 30 Pa , which are an order of magnitude larger than the corresponding wall shear-stress for the same hemodynamic conditions. The inelastic models systematically underestimate the shear-stress and overestimate the mean velocity. Finally, the TEVP model accurately predicts the phase-lags between the pressure-gradient and the flow-rate, and between the pressure-gradient and the structure-parameter for the whole range of the frequencies expressed in terms of the Womersley number, indicating its physical completeness and robustness.

Approximation Solution for the Zener Impact Theory

Collisions can be classified as completely elastic or inelastic. Collision mechanics theory has gradually developed from elastic to inelastic collision theories. Based on the Hertz elastic collision contact theory and Zener inelastic collision theory model, we derive and explain the Hertz and Zener collision theory model equations in detail in this study and establish the Zener inelastic collision theory, which is a simple and fast calculation of the approximate solution to the nonlinear differential equations of motion. We propose an approximate formula to obtain the Zener nonlinear differential equation of motion in a simple manner. The approximate solution determines the relevant values of the collision force, material displacement, velocity, and contact time.