A Stochastic Dynamics Method for Time-Varying Damping Depending on Temperature/Frequency for Several Alloy Materials

In the field of aerospace and advanced equipment manufacturing, accurate response analysis has been paid more attention, requiring a more comprehensive study of the variation of mechanical parameters with the service environment. The damping variation characteristics of 304 aluminum alloy, Sa564 high-strength alloy, GW63K magnesium alloy, and Q235 steel were investigated in this paper, which plays a significant role in the dynamic responses of structures. Variable damping ratios were revealed by the damping tests based on a dynamic mechanical analysis (DMA). The numerical method of temperature/frequency-dependent damping parameters in stochastic dynamics was focused on. With a large variation in the damping ratio, a numerical constitutive relation for temperature-dependent damping was proposed, and an efficient stochastic dynamics method was derived to analyze the responses of structures based on the pseudo excitation method (PEM) and variable damping theory. The computational accuracy and validity of the proposed method are confirmed during the vibration tests and numerical analysis. Based on the comparison results of the two damping models and the experiments on GW63K alloy, we proved that the proposed method is more accurate to the real response of the actual engineering structure. The differences in dynamic responses between the constant damping and experiments are significant, and more attention should be paid to the numerical method of stochastic dynamic response of variable damping materials in the aviation and aerospace fields and high-temperature environments.


Introduction
Damping capacity plays a significant role in restraining vibration, absorbing sound, and reducing the noise of structures.A high-damping alloy material has increasing importance in fields with high requirements for vibration responses.Many studies [1,2] indicated that the damping characteristics of alloy materials are complicated and depend on temperature, frequency, strain amplitude etc.Several studies have investigated damping theories and test methods.Zhang [3] studied the damping of viscoelastic materials and indicated that their effective damping temperature range is usually extremely narrow and a comprehensive understanding of their temperature-dependent and frequency-dependent properties is lacking.The temperature damping capacities of Mg-3Al-1Zn-xSn alloys were investigated using a DMA under varied loading frequencies and Sn concentrations in [4].The addition of Sn resulted in a leftward shift of the P1 dislocation damping peak around 80 • C and a rightward shift of the P3 peak around 220 • C, which were attributed to the second phases in the vicinity of grain boundaries.Liu [5] prepared Fe-21Cr-4Al-based alloys by Sc addition, and their damping and mechanical properties were analyzed.The microstructure and phase composition were analyzed by XRD, SEM, and TEM.Ebrahimi [6] investigated the characterization of damping and its underlying mechanisms in a composite of CNTs/AZ91D, produced through cyclic extrusion.The study revealed an enhancement in high-temperature damping behavior compared to the original alloy.Ganguly [7] evaluated the damping behavior of SiCnp-reinforced AZ91 + 2.0Ca+ 0.3Sb (wt%) alloy in the temperature range of 25-350 • C and at frequencies of 1, 5, and 10 Hz; the research indicated that all nanocomposites exhibit higher storage modulus, loss modulus, and damping capacity than the alloy, and the damping capacity increased with an increase in temperature and decreased with frequency.Santosh's work [8] focused on exploring the damping properties of some copper-based ternary and quaternary SMAs using DMA by calculating internal friction.The results show that the peak value of damping depended on the temperature, and the addition of quaternary elements (Ni and Mn) decreases the transformation temperatures.Sakaguchi [9] revealed that a step-increase mode of damping capacity appears in the MnCuNiFe alloy when it is maintained at a temperature above 350 • C. The damping properties of Mg97Zn1Y2 alloy, Mg-10Gd-2Ye1Zn-0.5Zr-0.2Ndalloy, AZ91-Ca-Sb magnesium alloy, and magnesium alloy AZ61 were recently the focus of research, including the measurement and characterization methods of damping [7,[10][11][12][13][14].The damping characterization of high-strength high-manganese FeMn-based alloys at low strain amplitudes was studied.These FeMnCr-based composite alloys with ferritic FeCrMn layers showed higher damping capacity at low strain amplitudes [15].Fe-21Cr-4Al-based alloys and their damping and mechanical properties were analyzed.The microstructure and phase composition were analyzed by XRD, SEM, and TEM.The damping of alloys enables them to maintain high values of damping properties even over broad ranges of strain amplitude [5].
While many studies have concentrated on examining the impacts of alloy mechanical properties and damping capacity, the computation of random responses in alloy structures with genuine variable damping parameters has proven challenging due to a deficiency in both qualitative and quantitative analyses.Addressing dynamic analysis, especially for high-performance material structures in random conditions, the distinction between constant and variable damping computations could significantly influence structural assessment and optimization.Presently, there is a scarcity of studies addressing this complex issue.Hu [16] proposed a hybrid damping mechanical model to predict the hysteresis loops under different operating conditions.In this model, the total restoring torque of the hysteresis loop is divided into a high-order polynomial elastic torque, a hybrid damping torque, and an inertial torque.The superelastic response and damping capacity of quaternary Ni45•3Ti39.7Hf10Pd5 polycrystalline alloys were investigated in terms of temperature and loading frequency dependency [17].Fang [18] studied the vibration control of a rotating functionally gradient material beam under a thermal environment, and the difference was that the base beam is temperature-dependent.Jitender's study [19] focused on vibration control with variable damping and stiffness and presents a mathematical and Simulink model to analyze the performance of vibration absorbers, which was an effective model for variable damping systems.The study concluded that considering variable damping in practical applications is a very difficult task.A nonlinear geometric damping was studied in [20], in which a novel tuning methodology was proposed to analyze the shock loads and sustained vibrations induced by harmonic loads for arbitrary stiffness and damping.In the above studies, the characteristics of damping parameters of high-damping alloy materials changing with the service environment have been widely revealed, but the influence of such variable damping on the dynamic response of structure and equipment using high-damping alloy materials and the design method are less explored.
In this study, the characteristics of damping change were paid primary attention, and an efficient numerical method was derived to deal with the stochastic dynamics responses caused by these characteristics without exploring the microscopic mechanism that causes these changes.The temperature-dependent damping characteristics of several alloys (Sa564, GW63K, 304, and Q235) used in aerospace engineering were investigated using qualitative and quantitative analysis methods, and multiple evolution curves based on operating conditions were obtained to describe the damping change characteristics.A numerical constitutive relation was proposed to describe the temperature-dependent damping.Subsequently, a proposed quasi-non-stationary random vibration approach, relying on a classic pseudo excitation method (PEM) [21,22] and numerical constitutive relation, aims to address the dynamic responses of alloy structures by incorporating genuine variable material damping parameters.

Determination of Temperature-Dependent Damping of Alloy Materials
A simplified model of damping is commonly used in the analysis of the dynamics of vibration systems [23].Additionally, the theories of viscous, coulomb, and hysteretic damping are commonly used in structural analysis.In recent studies, some advanced analytical models of hysteresis were proposed.A novel rate-independent hysteretic model was formulated in [24], which adopts closed-form expressions for evaluating the output variable with important benefits in terms of computational efficiency and implementation ease.A Vaiana-Rosati model for complex rate-independent mechanical hysteresis phenomena was devised in [25].The closed-form expressions provided by the analytical reformulation are expressed in rate form to foster its use, especially in nonlinear dynamics.
The magnitude of a parameter is indicated by the extent of the hysteresis curve formed through the stress-strain curve under cyclic loading.General expression forms include specific damping ψ, loss factor η, loss tangent tan ϕ, logarithmic decrement δ, reciprocal of quality factor Q −1 , and damping ζ.A dynamic mechanical analysis (DMA) can perform measurements over a wide range of temperature and frequency changes.Based on the sensitivity of damping to working conditions, the DMA can be used to measure and demarcate the damping ratio.Its principle can be described by the stress-strain relationship under cyclic loading: in which ε 0 represents the strain amplitude, σ 0 signifies the stress amplitude, ω denotes the angular frequency, and t corresponds to time; additionally, φ is the loss angle, indicating the phase lag of strain concerning the stress phase difference.The interplay between strain and stress produces a hysteresis loop, resulting in the dissipation of mechanical vibration energy.
The capacity to dissipate this vibration energy is recognized as the damping capacity [25].
It can be calculated by the following: where E ′′ is the loss modulus, E ′ is the storage modulus, φ is the loss angle that strain hysteresis stress, and tan φ is the loss tangent.In a DMA, η and tan φ are used to characterize the damping capacity.

Materials and Methods
Dynamic mechanical analysis (DMA) uses viscosity damping theory, coulomb damping theory, and hysteresis damping theory to measure the material damping parameters over large variations in frequency and temperature [26].This study aims to characterize damping characterization for alloys resonance DMA methods, especially with the beam model.Here, a beam is excited and oscillates at its resonance.To overcome this drawback in this investigation, a forced frequency DMA setup was used.The measuring instrument is shown in Figure 1.
This paper uses typical alloy materials to test the damping properties, including 304 aluminum alloy, Sa564 high-strength alloy, GW63K magnesium alloy, and Q235 steel.The corresponding material parameters are shown in Table 1.The specimen is designed according to the specimen requirements of the DMA test, as shown in Figure 2. model.Here, a beam is excited and oscillates at its resonance.To overcome this drawb in this investigation, a forced frequency DMA setup was used.The measuring instrum is shown in Figure 1.This paper uses typical alloy materials to test the damping properties, including aluminum alloy, Sa564 high-strength alloy, GW63K magnesium alloy, and Q235 steel.T corresponding material parameters are shown in Table 1.The specimen is designed cording to the specimen requirements of the DMA test, as shown in Figure 2.  Here, the damping characteristics are focused on three aspects: the relationship tween damping and strain amplitude, the relationship between damping and tempe ture, and the relationship between damping and excitation frequency.A set of compari tests is first designed to guarantee the accuracy of the instrument's measurements.In t test, two calibrated devices were used to test GW63K magnesium alloy specimens, a the comparison results are shown in Figure 3.  model.Here, a beam is excited and oscillates at its resonance.To overcome this drawback in this investigation, a forced frequency DMA setup was used.The measuring instrument is shown in Figure 1.This paper uses typical alloy materials to test the damping properties, including 304 aluminum alloy, Sa564 high-strength alloy, GW63K magnesium alloy, and Q235 steel.The corresponding material parameters are shown in Table 1.The specimen is designed according to the specimen requirements of the DMA test, as shown in Figure 2.  Here, the damping characteristics are focused on three aspects: the relationship between damping and strain amplitude, the relationship between damping and temperature, and the relationship between damping and excitation frequency.A set of comparison tests is first designed to guarantee the accuracy of the instrument's measurements.In this test, two calibrated devices were used to test GW63K magnesium alloy specimens, and the comparison results are shown in Figure 3. Here, the damping characteristics are focused on three aspects: the relationship between damping and strain amplitude, the relationship between damping and temperature, and the relationship between damping and excitation frequency.A set of comparison tests is first designed to guarantee the accuracy of the instrument's measurements.In this test, two calibrated devices were used to test GW63K magnesium alloy specimens, and the comparison results are shown in Figure 3.
Figure 3 depicts the damping parameters of 304 alloy obtained through various test pieces at room temperature (25 • C) with a frequency of 1 Hz.The results exhibit consistent variation tendencies and convergence trends.Generally, the damping ratio of 304 alloys increases with the strain amplitude and stabilizes at tanφ = 0.0180, around 10 −1 % of the strain amplitude, often termed the 'critical unpinning strain amplitude point'.This specific point is commonly adopted as the constant damping ratio in engineering practices.To validate measurement accuracy, a cantilever beam model with length 4.8 cm, width 1.2 cm, and height 0.2 cm was devised for the vibrating reed method test, and the results were compared with a numerical simulation employing the same damping ratio.
The mechanism of the evolution of the damping characters of alloy materials with temperature and strain amplitude can be explained by optical microstructure, such as Figure 4.Under low strain amplitude at room temperature, the deformation of the material is twin-coordinated deformation of grain, with clear grain boundary and no dynamic recrystallization.With the increase in temperature, dynamic recrystallization occurs in the grain boundary and grain interior, which changes from basal slip to non-basal slip.The number of slip systems increases, resulting in dislocation movement, which improves the damping performance.At the same temperature, the increase in strain amplitude can increase the twinning tendency of the grain boundary, and the deformation mechanisms such as the cross-slip of screw dislocation and the migration of grain boundary cannot be carried out.The dislocation movement is hindered, and the damping increases greatly.Figure 3 depicts the damping parameters of 304 alloy obtained through various tes pieces at room temperature (25 °C) with a frequency of 1 Hz.The results exhibit consisten variation tendencies and convergence trends.Generally, the damping ratio of 304 alloy increases with the strain amplitude and stabilizes at tanφ = 0.0180, around 10 −1 % of th strain amplitude, often termed the 'critical unpinning strain amplitude point'.This spe cific point is commonly adopted as the constant damping ratio in engineering practices To validate measurement accuracy, a cantilever beam model with length 4.8 cm, width 1. cm, and height 0.2 cm was devised for the vibrating reed method test, and the results wer compared with a numerical simulation employing the same damping ratio.
The mechanism of the evolution of the damping characters of alloy materials wit temperature and strain amplitude can be explained by optical microstructure, such as Fig ure 4.Under low strain amplitude at room temperature, the deformation of the materia is twin-coordinated deformation of grain, with clear grain boundary and no dynamic re crystallization.With the increase in temperature, dynamic recrystallization occurs in th grain boundary and grain interior, which changes from basal slip to non-basal slip.Th number of slip systems increases, resulting in dislocation movement, which improves th damping performance.At the same temperature, the increase in strain amplitude can in crease the twinning tendency of the grain boundary, and the deformation mechanism such as the cross-slip of screw dislocation and the migration of grain boundary cannot b carried out.The dislocation movement is hindered, and the damping increases greatly.It has also been reported that the damping capacities of alloy are dependent on man working conditions, such as temperature, excitation frequency, physical dimensions, an so on [1,2].To reveal temperature-dependent characteristics, some qualitative and quan titative analysis experiments were designed.The results are shown in Figures 5-8.  Figure 3 depicts the damping parameters of 304 alloy obtained through various test pieces at room temperature (25 °C) with a frequency of 1 Hz.The results exhibit consistent variation tendencies and convergence trends.Generally, the damping ratio of 304 alloys increases with the strain amplitude and stabilizes at tanφ = 0.0180, around 10 −1 % of the strain amplitude, often termed the 'critical unpinning strain amplitude point'.This specific point is commonly adopted as the constant damping ratio in engineering practices.To validate measurement accuracy, a cantilever beam model with length 4.8 cm, width 1.2 cm, and height 0.2 cm was devised for the vibrating reed method test, and the results were compared with a numerical simulation employing the same damping ratio.
The mechanism of the evolution of the damping characters of alloy materials with temperature and strain amplitude can be explained by optical microstructure, such as Figure 4.Under low strain amplitude at room temperature, the deformation of the material is twin-coordinated deformation of grain, with clear grain boundary and no dynamic recrystallization.With the increase in temperature, dynamic recrystallization occurs in the grain boundary and grain interior, which changes from basal slip to non-basal slip.The number of slip systems increases, resulting in dislocation movement, which improves the damping performance.At the same temperature, the increase in strain amplitude can increase the twinning tendency of the grain boundary, and the deformation mechanisms such as the cross-slip of screw dislocation and the migration of grain boundary cannot be carried out.The dislocation movement is hindered, and the damping increases greatly.It has also been reported that the damping capacities of alloy are dependent on many working conditions, such as temperature, excitation frequency, physical dimensions, and so on [1,2].To reveal temperature-dependent characteristics, some qualitative and quantitative analysis experiments were designed.The results are shown in Figures 5-8.It has also been reported that the damping capacities of alloy are dependent on many working conditions, such as temperature, excitation frequency, physical dimensions, and so on [1,2].To reveal temperature-dependent characteristics, some qualitative and quantitative analysis experiments were designed.The results are shown in Figures 5-8.The variation characteristics of alloy temperature-dependent damping were investigated, as shown in Figures 5-8.The interaction of temperature with the damping properties of different alloys can be discussed.
The temperature has a great influence on the damping properties of the four alloy materials studied; the maximum change of damping performance is 218% and the minimum is 167% in the test range of 25 °C to 350 °C allowed by DMA equipment.
According to the test results, the damping performance of Sa564 decreases monotonically with the increase in temperature, and the damping performance of GW63K increases monotonically with the increase in temperature.The damping performance of Q235 decreases first and then increases, and the damping performance of 304 increases first and then decreases, both of which produce an inflexion point at around 200 °C.This indicates that the influence of temperature on the damping properties of different alloy materials is inconsistent and should be discussed separately based on the experiments.
In this paper, the phenomenon of damping as a function of temperature is focused on, and the research pays attention to the difficulties brought by the phenomenon of The variation characteristics of alloy temperature-dependent damping were investigated, as shown in Figures 5-8.The interaction of temperature with the damping properties of different alloys can be discussed.
The temperature has a great influence on the damping properties of the four alloy materials studied; the maximum change of damping performance is 218% and the minimum is 167% in the test range of 25 • C to 350 • C allowed by DMA equipment.
According to the test results, the damping performance of Sa564 decreases monotonically with the increase in temperature, and the damping performance of GW63K increases monotonically with the increase in temperature.The damping performance of Q235 decreases first and then increases, and the damping performance of 304 increases first and then decreases, both of which produce an inflexion point at around 200 • C.This indicates that the influence of temperature on the damping properties of different alloy materials is inconsistent and should be discussed separately based on the experiments.
In this paper, the phenomenon of damping as a function of temperature is focused on, and the research pays attention to the difficulties brought by the phenomenon of damping depending on temperature in the calculation and analysis of structural dynamic response and proposes computational mechanics methods.Therefore, the microscopic mechanism of this temperature-dependent damping change is not the focus of research.
In the test results, the experimental phenomena of some alloy materials have more significant divergence, contributed to by the blessing accuracy, the accuracy of the heating process, the dimensional accuracy of the specimen, and the environmental influence.It is still difficult to achieve complete consistency during experiments.
However, the damping parameters depended on temperature, and the variation characteristic can be obtained, which is of positive significance for subsequent research on computational mechanics methods.

Numerical Stochastic Dynamics Method for Temperature-Dependent Damping
The PEM method is a mature, accurate, and efficient computational mechanics method for structural random vibration, which was proposed by Prof Jiahao Lin, and is widely used in aerospace, bridge and structural seismic engineering [21,22].In typical applications, the damping parameter is set to constant damping in stationary and non-stationary random vibrations.Based on the proposed PEM method, it is impossible to deal with variable damping during vibration.In this part, the work focuses on dealing with the dynamic responses of variable damping caused by the environmental influence of this alloy material based on PEM.This method applies to the above four alloy materials and has guiding significance for any variable damping that can establish a time function.
Generally, alloy structures are used as high-performance materials in complicated working conditions and are subject to random loads in a wide frequency range and harsh temperatures.Random dynamics analysis methods under a constant damping parameter have been developed, including the response spectrum method, the step-by-step integration method, and the PEM.However, if a variable damping parameter is considered instead of a constant damping parameter in a random dynamic system, the linear random vibration issue will evolve into a non-linear random dynamic issue, which is more difficult.The differences between linear and non-linear systems can be described by the following Equation (5) [21,22].
where M, K, and C comprise the non-time-varying matrix, P is the excitation amplitude, g(t) is the modulation function, and x(t) is the random excitation which varies over time.In general, the responses of Equation ( 5) can be obtained quickly and accurately based on the PEM.For the excitation of a stationary random process vector x(t), the pseudo excitation x = √ S xx e iωt can be constructed according to its self-spectral spectrum density S xx (ω); the structural response spectrum characteristics can be obtained by the following equations.
where H(ω) is the structural frequency response function matrix, * is a complex conjugate, S yy is the self-spectral density matrix of responses, and S xy is the cross-spectral density matrix of response.
The correlation matrix of y(t) can be expressed as follows: and Equation ( 7) can be calculated by the Wiener-Khintchine theorem.
when we combine Equations ( 7) and ( 8), we get the following: where For the time-varying auto-power spectrum matrix S yy (ω, t): According to the above derivation, if we construct the pseudo-excitation as follows S xx (ω)g(t)e iωt the response of time t is as follows: Therefore, the non-stationary random vibration can be solved by the pseudo-excitation method.The result is as shown in Equation (13).
For temperature/frequency-dependent damping model, the damping matrix C is variable parameter matrix.Unlike the constant damping Equation ( 5), the damping term C y, .y is not only related to the velocity response but also related to the displacement response (strain amplitude), excitation frequency, and working temperature.Thus, a linear random dynamic system is evolved into a complex non-linear stochastic dynamic system.
The problem of Equation ( 6) mainly involves a non-linear analysis.The non-linear parts of dynamics can be divided into two categories according to the damping change mechanism.One type is the damping variation caused by the strain amplitude associated with the structural response.This type needs to be solved based on the linear theory of the non-linear problem.The other can be attributed to time-dependent damping variation, such as temperature and excitation frequency.The second type of problem can transform the damping properties into a time-dependent process, which is the target of this paper for stochastic dynamics.Sidorov [27] proposed a non-local-in-time damping models, in which the material had separate parts of locality and nonlocality, and the damping properties were also divided into nonlocal and local parts.Within damping-with-memory models, the internal damping of a structure at the current time is assumed to be dependent not only on the instant strain-rate magnitude or displacement-velocity magnitude but also on the strain-rate or velocity values along the previous time history; the equation of motion could be solved using the modified Newmark method.The new model brings a new idea to the solution of variable damping dynamics.Unlike [27], temperature/frequency-dependent damping can be discretized in time steps and converted into a time-dependent function, since both the heating process and the excitation process are time-dependent.In this way, the equation of motion contains one type of damping, which is a function of time.According to a certain characteristic of loading, both can be converted into time-dependent parameters as follows: where T e , which is a set of vectors, represents the change characteristic of the structural working environment temperature with time.The recursive function can be performed according to the curve of the damping pair and the function relationship: In these equations, η T represents the evolution characteristic of material damping with temperature.In this way, the characteristic of variable damping coefficient was evolved into a curve that changes with time in which the method of solving such time-varying damping was given more attention.Although the change characteristics of damping evolve into a time-varying function, the orthogonal characteristics are still maintained.The Rayleigh damping theory can be used in a random dynamic system and can be expressed as follows: Here, α(t) and β(t) denote the Rayleigh damping ratio parameters associated with time-dependent damping, while ω i and ω j represent the modal circular frequencies for each order.For calculation convenience, the decisive first-and second-order modes are used to calculate.η t (t) signifies the component of the time-dependent material damping ratio vector, with the subscript 'T' indicating its variation with temperature.The original random dynamic equation corresponding to these parameters can be expressed as follows: The random dynamics problem should still be defined as a stationary random vibration problem.Discretization in the time domain is necessary to solve this dynamic problem because the damping matrix C of the system is time-dependent.Referring to the constant damping non-stationary random vibration solution method, Equation ( 13) is a quasi-nonstationary random dynamic system, which is solved based on the non-stationary method.The solution form of the pseudo excitation method is as follows.At a certain time t, the pseudo excitation is determined as follows: During modal reduced order to the system, taking the first q-order modal ϕ Its responses can be obtained through S yy (ω) = y * • y T .This deviates from a nonstationary system across the entire time range due to the time modulation of the pseudo excitation.Typically, a non-stationary random pseudo excitation is represented as follows: In the quasi-non-stationary process involving dynamic damping, the pseudo-excitation Equation ( 23) from the stationary process remains applicable.However, its time-varying features are manifested in the evolution of the damping matrix at the next time step t + ∆t: Dynamic analysis of this dynamic damping can be easily acquired through numerical methods within a linear framework.
A similar approach is proposed to analyze random vibration responses to address the features of frequency-dependent damping.Adhering to a specific frequency characteristic, damping variations can be translated into time-dependent parameters as follows: where F e , which is a set of vectors, represents the change characteristic of the structural working force frequency with time.The recursive function can be performed according to the curve of the damping pair and the function relationship:

Numerical Constitutive Relations for Temperature-Dependent Damping
Based on the evolution characteristic η T of material damping with temperature, the random vibration Equation ( 5) with constant damping can be drawn into the temperaturedependent damping stochastic dynamics Equation (13).A quasi-non-stationary method with temperature-dependent damping is proposed to calculate the random responses.However, damping constitutive relation is a key point, which determines the accuracy and computational efficiency of the proposed algorithm.The traditional physical constitutive relation is complicated and involves many factors, which makes it difficult to analyse the dynamic equation.Here, a numerical constitutive relation fitted according to the temperature curve is proposed to simulate the complex temperature and damping dependence and improve the efficiency of solving the dynamic equation.Take the temperature-dependent damping of GW63K alloy as an example, as shown in Figure 8.The polynomial interpolation method can be used to simulate the temperature-dependent damping curves.
In Figure 9, combining polynomials and least square method, the constitutive relation of temperature-dependent damping can be simulated as follows in Equation (27) or Equation ( 28): where a, b, c, d, f are fitting parameters.The order and function form of the numerical simulating model can be determined according to the requirements of accuracy and computational efficiency.In general, the time-dependent higher-order terms will cause higherorder nonlinearity.Based on the consideration of computational difficulty and accuracy, low-order nonlinear dynamic equations are adopted in this paper.

Dynamic Analysis and Experimental Validation of GW63K Alloy Considering Temperature-Dependent Damping
In the presented example, a numerical simulation and experimental verification model were implemented for the GW63K magnesium alloy in [28].The analysis involved a cantilever beam mode subjected to continuous sinusoidal excitation at the free end, with amplitude responses collected across the temperature range of 25-350 °C using Dynamic Mechanical Analysis (DMA).The amplitude responses were calculated using the constant damping ratio for a strain amplitude at the 10-1 level.Meanwhile, the proposed method employed the temperature-dependent damping ratio to compute amplitude responses.To mitigate high-order non-linearity in the dynamic system, linear interpolation was applied to linearly fit the variable damping.Figure 8 illustrates the numerical constitutive relation of simulation Curve I for temperature-dependent damping.The experimental setup is depicted in Figure 10, and detailed material and model parameters are provided in Table 2.The heating rate of the DMA is 0.194 °C/s.The unit of time was seconds, and the damping evolution can be expressed as a determined time function in conjunction with Figure 10.The response results obtained by the three methods are summarized in the following figures.
The examination of the test piece, through both experimental and numerical calculations, reveals a notable impact of the damping variable model on acceleration and dis-

Dynamic Analysis and Experimental Validation of GW63K Alloy Considering Temperature-Dependent Damping
In the presented example, a numerical simulation and experimental verification model were implemented for the GW63K magnesium alloy in [28].The analysis involved a cantilever beam mode subjected to continuous sinusoidal excitation at the free end, with amplitude responses collected across the temperature range of 25-350 • C using Dynamic Mechanical Analysis (DMA).The amplitude responses were calculated using the constant damping ratio for a strain amplitude at the 10-1 level.Meanwhile, the proposed method employed the temperature-dependent damping ratio to compute amplitude responses.To mitigate high-order non-linearity in the dynamic system, linear interpolation was applied to linearly fit the variable damping.Figure 8 illustrates the numerical constitutive relation of simulation Curve I for temperature-dependent damping.The experimental setup is depicted in Figure 10, and detailed material and model parameters are provided in Table 2.

Dynamic Analysis and Experimental Validation of GW63K Alloy Considering Temperat Dependent Damping
In the presented example, a numerical simulation and experimental verific model were implemented for the GW63K magnesium alloy in [28].The analysis invo a cantilever beam mode subjected to continuous sinusoidal excitation at the free end, amplitude responses collected across the temperature range of 25-350 °C using Dyn Mechanical Analysis (DMA).The amplitude responses were calculated using the con damping ratio for a strain amplitude at the 10-1 level.Meanwhile, the proposed me employed the temperature-dependent damping ratio to compute amplitude respo To mitigate high-order non-linearity in the dynamic system, linear interpolation wa plied to linearly fit the variable damping.Figure 8 illustrates the numerical constit relation of simulation Curve I for temperature-dependent damping.The experim setup is depicted in Figure 10, and detailed material and model parameters are prov in Table 2.The heating rate of the DMA is 0.194 °C/s.The unit of time was seconds, and damping evolution can be expressed as a determined time function in conjunction Figure 10.The response results obtained by the three methods are summarized in th lowing figures.
The examination of the test piece, through both experimental and numerical cal tions, reveals a notable impact of the damping variable model on acceleration and placement.The heating rate of the DMA is 0.194 • C/s.The unit of time was seconds, and the damping evolution can be expressed as a determined time function in conjunction with Figure 10.The response results obtained by the three methods are summarized in the following figures.
The examination of the test piece, through both experimental and numerical calculations, reveals a notable impact of the damping variable model on acceleration and displacement.Figures 11-13 present a comparison between the data obtained from the Dynamic Mechanical Analysis (DMA) experiment and the results derived from temperature-dependent damping and constant damping.Notably, constant damping maintains consistency in the energy region under continuous sinusoidal loading, ensuring a stable response region.The experimental findings align well with the outcomes of the temperature-dependent damping calculation; the maximum difference between the test results and the calculated results was only 4.32%, exhibiting a relatively consistent change and numerical correspondence.However, during the stable tightening period when the test fixture is excited, some divergence was observed.For displacement response, at a temperature of 350 • C, the variable damping model response is 0.00335 m, and the constant damping response is 0.00635 m; the discrepancy in response at elevated temperatures is substantial, reaching approximately 47.24%.Furthermore, the numerical results indicate a significant difference in displacement responses with increasing temperature.dependent damping and constant damping.Notably, constant damping maintains consistency in the energy region under continuous sinusoidal loading, ensuring a stable response region.The experimental findings align well with the outcomes of the temperature-dependent damping calculation; the maximum difference between the test results and the calculated results was only 4.32%, exhibiting a relatively consistent change and numerical correspondence.However, during the stable tightening period when the test fixture is excited, some divergence was observed.For displacement response, at a temperature of 350 °C, the variable damping model response is 0.00335 m, and the constant damping response is 0.00635 m; the discrepancy in response at elevated temperatures is substantial, reaching approximately 47.24%.Furthermore, the numerical results indicate a significant difference in displacement responses with increasing temperature.dependent damping and constant damping.Notably, constant damping maintains consistency in the energy region under continuous sinusoidal loading, ensuring a stable response region.The experimental findings align well with the outcomes of the temperature-dependent damping calculation; the maximum difference between the test results and the calculated results was only 4.32%, exhibiting a relatively consistent change and numerical correspondence.However, during the stable tightening period when the test fixture is excited, some divergence was observed.For displacement response, at a temperature of 350 °C, the variable damping model response is 0.00335 m, and the constant damping response is 0.00635 m; the discrepancy in response at elevated temperatures is substantial, reaching approximately 47.24%.Furthermore, the numerical results indicate a significant difference in displacement responses with increasing temperature.

Responses Analysis of Engineering Structures with Temperature-Dependent Damping by Abaqus Software
In aeronautical engineering, the GW63K magnesium alloy engineering structure is widely used for weight reduction and vibration reduction.In this example, the responses of the magnesium alloy support used in aircraft were calculated based on both constant (the first development) and variable (the second development) damping.Take the same damping in 25 °C-350 °C of GW63K alloy as an example.As shown in Figure 9, the polynomial interpolation method can be used to simulate the temperature-dependent damping curves.
Based on ABAQUS 6.10.1 software, for this bearing structure, the meshing method of tetrahedral element and mixed hexahedral element is selected according to its geometric characteristics and the characteristics of meshing.The finite element mesh in the process analysis can be set as the first-order element to improve the efficiency of finite element simulation and optimization.In the final result verification, it can be set as the secondorder element to improve the analysis and verification accuracy.In this model, 78,952 hexahedral elements and 12,874 tetrahedral elements were adopted.
Figure 14 shows the three-dimensional model of the magnesium alloy support.Table 1 shows the material parameters.Table 3 shows the natural frequencies of the magnesium alloy support.The duration of the sinusoidal excitation is 10 s.The fixed constraint is applied at the four corner screw holes.The sinusoidal excitation is applied at the two corner screw holes in the mid-span position.The dynamic response of a point at the mid-span position is

Responses Analysis of Engineering Structures with Temperature-Dependent Damping by Abaqus Software
In aeronautical engineering, the GW63K magnesium alloy engineering structure is widely used for weight reduction and vibration reduction.In this example, the responses of the magnesium alloy support used in aircraft were calculated based on both constant (the first development) and variable (the second development) damping.Take the same damping in 25 • C-350 • C of GW63K alloy as an example.As shown in Figure 9, the polynomial interpolation method can be used to simulate the temperature-dependent damping curves.
Based on ABAQUS 6.10.1 software, for this bearing structure, the meshing method of tetrahedral element and mixed hexahedral element is selected according to its geometric characteristics and the characteristics of meshing.The finite element mesh in the process analysis can be set as the first-order element to improve the efficiency of finite element simulation and optimization.In the final result verification, it can be set as the second-order element to improve the analysis and verification accuracy.In this model, 78,952 hexahedral elements and 12,874 tetrahedral elements were adopted.
Figure 14 shows the three-dimensional model of the magnesium alloy support.Table 1 shows the material parameters.Table 3 shows the natural frequencies of the magnesium alloy support.

Responses Analysis of Engineering Structures with Temperature-Dependent Damping Abaqus Software
In aeronautical engineering, the GW63K magnesium alloy engineering struct widely used for weight reduction and vibration reduction.In this example, the resp of the magnesium alloy support used in aircraft were calculated based on both con (the first development) and variable (the second development) damping.Take the damping in 25 °C-350 °C of GW63K alloy as an example.As shown in Figure 9, the nomial interpolation method can be used to simulate the temperature-dependent d ing curves.
Based on ABAQUS 6.10.1 software, for this bearing structure, the meshing meth tetrahedral element and mixed hexahedral element is selected according to its geom characteristics and the characteristics of meshing.The finite element mesh in the pr analysis can be set as the first-order element to improve the efficiency of finite ele simulation and optimization.In the final result verification, it can be set as the se order element to improve the analysis and verification accuracy.In this model, 78,95 ahedral elements and 12,874 tetrahedral elements were adopted.
Figure 14 shows the three-dimensional model of the magnesium alloy support. 1 shows the material parameters.Table 3 shows the natural frequencies of the magne alloy support.The duration of the sinusoidal excitation is 10 s.The fixed constraint is applied four corner screw holes.The sinusoidal excitation is applied at the two corner screw in the mid-span position.The dynamic response of a point at the mid-span posit  The duration of the sinusoidal excitation is 10 s.The fixed constraint is applied at the four corner screw holes.The sinusoidal excitation is applied at the two corner screw holes in the mid-span position.The dynamic response of a point at the mid-span position is extracted and compared with the corresponding dynamic response obtained with constant damping.The envelope curves of the dynamic responses are shown in Figures 15-20     It can be seen in the actual engineering structure that the dynamic response amplitude of the support gradually decreases with an increase in the damping parameter.For the constant damping system, the energy input and dissipation are in an equilibrium state.It can be seen in the actual engineering structure that the dynamic response amplitude of the support gradually decreases with an increase in the damping parameter.For the constant damping system, the energy input and dissipation are in an equilibrium state.It can be seen in the actual engineering structure that the dynamic response amplitude of the support gradually decreases with an increase in the damping parameter.For the constant damping system, the energy input and dissipation are in an equilibrium state.It can be seen in the actual engineering structure that the dynamic response amplitude of the support gradually decreases with an increase in the damping parameter.For the constant damping system, the energy input and dissipation are in an equilibrium state.In a variable damping system, the energy input is less than the dissipation.Therefore, finally, the amplitude of the displacement is reduced to 47.72% of the maximum, the amplitude of the speed is reduced to 47.73% of the maximum, and the amplitude of the acceleration is reduced to 47.74% of the maximum.This demonstrates that with an increase in the material damping performance of the structure, the energy dissipation of the structure increases, the kinetic energy decreases, the damping force increases, and the vibration reduction effects are strengthened.
finally, the amplitude of the displacement is reduced to 47.72% of the maximum, the amplitude of the speed is reduced to 47.73% of the maximum, and the amplitude of the acceleration is reduced to 47.74% of the maximum.This demonstrates that with an increase in the material damping performance of the structure, the energy dissipation of the structure increases, the kinetic energy decreases, the damping force increases, and the vibration reduction effects are strengthened.The dynamic responses of variable damping and constant damping in the middle of the support have obvious differences.Furthermore, at the two ends of the support, the dynamic responses were more similar because of the fixed constraint and because they were far from the sinusoidal excitation.For example, Figures 18 and 19 show the displacement response of the support based on the second development of ABAQUS software.The displacement responses of the variable damping in the middle section were concentrated in the range 7.798 × 10 −2 mm-8.789× 10 −2 mm.The displacement responses of the constant damping in the middle section were concentrated in the range 1.077 × 10 −1 mm-1.176× 10 −2 mm.At two ends, displacement responses of the constant damping and the variable damping are both concentrated in the dark blue area.These results show that there is a significant difference between the constant damping and variable damping, and this difference should not be ignored.The method presented in this paper can be used to calculate the random problem with variable damping, but its analytical capability needs Figure 17 shows the displacement response curves of the structure with constant damping.It can be seen that the dynamic response amplitudes exhibited during constant damping are stabilized after transient vibration.This characteristic is different from that of variable damping, wherein the amplitudes of the dynamic response vary with time.
The dynamic responses of variable damping and constant damping in the middle of the support have obvious differences.Furthermore, at the two ends of the support, the dynamic responses were more similar because of the fixed constraint and because they were far from the sinusoidal excitation.For example, Figures 18 and 19 show the displacement response of the support based on the second development of ABAQUS software.The displacement responses of the variable damping in the middle section were concentrated in the range 7.798 × 10 −2 mm-8.789× 10 −2 mm.The displacement responses of the constant damping in the middle section were concentrated in the range 1.077 × 10 −1 mm-1.176× 10 −2 mm.At two ends, displacement responses of the constant damping and the variable damping are both concentrated in the dark blue area.These results show that there is a significant difference between the constant damping and variable damping, and this difference should not be ignored.The method presented in this paper can be used to calculate the random problem with variable damping, but its analytical capability needs to be further verified by experiments.

The Dynamic Analysis of GW63K Alloy with Frequency-Dependent Damping
In this section, the focus was on investigating a magnesium alloy bearing within an engineering structure, specifically examining the response characteristics under wide-band random load excitation for both frequency-dependent damping and constant damping in [28].The bearing model is depicted in Figures 21 and 22. Utilizing two load input spectra detailed in Table 4, dynamic experiments were conducted to obtain time-domain and frequency-domain responses under random excitations.The numerical simulation involved two types of damping parameters, constant and frequency-dependent, established through a linear relationship of frequency-dependent damping, illustrated in Figure 23, specifically tailored for the GW63K alloy.The simulation computations were performed using the Hyperworks 14.0 software through secondary development.The dynamic responses of variable damping and constant damping in the middle of the support have obvious differences.Furthermore, at the two ends of the support, the dynamic responses were more similar because of the fixed constraint and because they were far from the sinusoidal excitation.For example, Figures 18 and 19 show the displacement response of the support based on the second development of ABAQUS software.The displacement responses of the variable damping in the middle section were concentrated in the range 7.798 × 10 −2 mm-8.789× 10 −2 mm.The displacement responses of the constant damping in the middle section were concentrated in the range 1.077 × 10 −1 mm-1.176× 10 −2 mm.At two ends, displacement responses of the constant damping and the variable damping are both concentrated in the dark blue area.These results show that there is a significant difference between the constant damping and variable damping, and this difference should not be ignored.The method presented in this paper can be used to calculate the random problem with variable damping, but its analytical capability needs to be further verified by experiments.

The Dynamic Analysis of GW63K Alloy with Frequency-Dependent Damping
In this section, the focus was on investigating a magnesium alloy bearing within an engineering structure, specifically examining the response characteristics under wideband random load excitation for both frequency-dependent damping and constant damping in [28].The bearing model is depicted in Figures 21 and 22. Utilizing two load input spectra detailed in Table 4, dynamic experiments were conducted to obtain time-domain and frequency-domain responses under random excitations.The numerical simulation involved two types of damping parameters, constant and frequency-dependent, established through a linear relationship of frequency-dependent damping, illustrated in Figure 23, specifically tailored for the GW63K alloy.The simulation computations were performed using the Hyperworks 14.0 software through secondary development.Following the table's loading guidelines, it was configured to execute the 5-580 Hz cycle within a 30 s timeframe.In Figure 23, combining the linear equation method, the constitutive relation of frequency-dependent damping can be simulated as Black Curve based on 3 Specimen DMA test curve.The linear equation is as follows (29): η F = 5.9 −5 T( f ) + 0.02808 (29)    Following the table's loading guidelines, it was configured to execute the 5-580 loading cycle within a 30 s timeframe.In Figure 23  The acceleration time domain and frequency domain experimental curves of the shaking table are illustrated in Figure 24.Comparison between the results obtained from the frequency-dependent damping analysis method and the constant damping analysis method was conducted against the experimental data.The comparative outcomes are illustrated in Figure 25, and a detailed summary is provided in Table 5.Comparison between the results obtained from the frequency-dependent damping analysis method and the constant damping analysis method was conducted against the experimental data.The comparative outcomes are illustrated in Figure 25, and a detailed summary is provided in Table 5.
Observing the results, it becomes evident that the computed responses of constant damping and variable damping are relatively synchronized in the lower excitation frequency range, demonstrating a high level of consistency with the experimental results.However, as the excitation frequency increases, the impact of dynamic damping starts to diminish, leading to delayed attenuation.The disparity between the dynamic damping response and the constant damping calculation gradually amplifies, accompanied by a change in phase.Comparison between the results obtained from the frequency-dependent damping analysis method and the constant damping analysis method was conducted against the experimental data.The comparative outcomes are illustrated in Figure 25, and a detailed summary is provided in Table 5.  Observing the results, it becomes evident that the computed responses of constant damping and variable damping are relatively synchronized in the lower excitation frequency range, demonstrating a high level of consistency with the experimental results.However, as the excitation frequency increases, the impact of dynamic damping starts to diminish, leading to delayed attenuation.The disparity between the dynamic damping response and the constant damping calculation gradually amplifies, accompanied by a change in phase.
It is crucial to highlight that the experimental results and the calculated results of variable damping exhibit a greater consistency in both amplitude and phase compared to the results obtained with constant damping.The divergence between the outcomes of dynamic damping and constant damping becomes more pronounced in the higher frequency range.The numerical analysis reveals similarities in phase.The difference of RMS  It is crucial to highlight that the experimental results and the calculated results of variable damping exhibit a greater consistency in both amplitude and phase compared to the results obtained with constant damping.The divergence between the outcomes of dynamic damping and constant damping becomes more pronounced in the higher frequency range.The numerical analysis reveals similarities in phase.The difference of RMS between the test responses and numerical calculated results based on frequencydependent damping was only 2.7%, and the difference between the test responses and numerical calculated results based on constant damping reached 18.5%.The third example investigates the random dynamic responses based on three methods.It can be seen that the proposed variable damping method is closer to the experimental results in the response calculation at each order frequency response.The error between the constant damping and the experimental results is obvious.For test data, the first three orders are reliable to the experimental test results, and the high-order response is greatly affected by the excitation environment, which is also the reason for the large difference in the high-order responses.Such discrepancies demand significant consideration, especially in application fields with stringent precision requirements for the dynamic response of load-bearing structures.This underscores the importance of recognizing and addressing the impact of frequency sensitivity in alloy material damping on structural dynamics responses in high-performance damping materials and precision application domains.

Conclusions
In this study, the DMA test was used to measure the temperature-dependent damping characteristics of typical alloys.Furthermore, the variable damping dynamic response of the GW63K magnesium alloy support used in spacecraft was analyzed by ABAQUS.The conclusions drawn from the results of these tests are as follows.
(1) The DMA experiments revealed that the alloy damping parameters have high sensitivity to environment temperature.This variable damping characteristic is particularly important in fields with complex service environments and that require a high degree of precision from dynamic responses, such as the aeronautical industry.More test results also show that the damping of the four alloy materials is also affected by the excitation frequency of the environment, which is a frequency-dependent damping characteristic.

Materials 2024 ,
17, x FOR PEER REVIEW 4 o

Figure 3 .
Figure 3. Calibration results for alloy damping test.

Figure 3 .
Figure 3. Calibration results for alloy damping test.

Figure 9 .
Figure 9. Numerical constitutive relation of temperature-dependent damping of GW63K alloy
present a comparison between the data obtained from the namic Mechanical Analysis (DMA) experiment and the results derived from tempera

Figure 14 .
Figure 14.Three-dimensional model of magnesium alloy support.

Figure 14 .
Figure 14.Three-dimensional model of magnesium alloy support.
extracted and compared with the corresponding dynamic response obtained with constant damping.The envelope curves of the dynamic responses are shown in Figures15- 20.
extracted and compared with the corresponding dynamic response obtained with constant damping.The envelope curves of the dynamic responses are shown in Figures15- 20.

Figure 18 .
Figure 18.Displacement responses for the magnesium alloy support (constant damping).

Figure 17
Figure17shows the displacement response curves of the structure with constant damping.It can be seen that the dynamic response amplitudes exhibited during constant damping are stabilized after transient vibration.This characteristic is different from that of variable damping, wherein the amplitudes of the dynamic response vary with time.

Figure 18 .
Figure 18.Displacement responses for the magnesium alloy support (constant damping).

Figure 18 .
Figure 18.Displacement responses for the magnesium alloy support (constant damping).

Figure 17
Figure17shows the displacement response curves of the structure with constant damping.It can be seen that the dynamic response amplitudes exhibited during constant damping are stabilized after transient vibration.This characteristic is different from that of variable damping, wherein the amplitudes of the dynamic response vary with time.
, combining the linear equation met the constitutive relation of frequency-dependent damping can be simulated as B Curve based on 3 Specimen DMA test curve.The linear equation is as follows (29):The acceleration time domain and frequency domain experimental curves of shaking table are illustrated in Figure24.

Figure 24 .
Figure 24.The excitation signal in time domain of vibration test bench.

Figure 24 .
Figure 24.The excitation signal in time domain of vibration test bench.

Figure 24 .
Figure 24.The excitation signal in time domain of vibration test bench.

Figure 25 .
Figure 25.Acceleration response comparison of magnesium alloy bearing.

Figure 25 .
Figure 25.Acceleration response comparison of magnesium alloy bearing.

Table 1 .
The alloys material parameters.

Table 1 .
The alloys material parameters.

Table 1 .
The alloys material parameters.

Table 3 .
The natural frequencies of magnesium alloy support.

Table 3 .
The natural frequencies of magnesium alloy support.

Table 3 .
The natural frequencies of magnesium alloy support.

Table 4 .
Input working condition of random excitation.

Table 4 .
Input working condition of random excitation.

Table 4 .
Input working condition of random excitation.

Table 5 .
Random condition acceleration root-mean-square of bearing.

Table 5 .
Random condition acceleration root-mean-square of bearing.