Next Article in Journal
State-of-the-Art Review on the Seismic Performance Assessment of On-Ground Steel Cylindrical Tanks
Next Article in Special Issue
Efficient Modal Identification and Optimal Sensor Placement via Dynamic DIC Measurement and Feature-Based Data Compression
Previous Article in Journal
Ultrasonic Signal Time-Expansion Using DAC Frequency Modulation
Previous Article in Special Issue
Basic Study on Mechanical Vibration Suppression System Using 2-Degree-of-Freedom Vibration Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Approach to Optimize the Dynamic Behaviour of Structures Considering Structural Durability

Fraunhofer Institute for Structural Durability and System Reliability LBF, 64289 Darmstadt, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Vibration 2023, 6(3), 477-493; https://doi.org/10.3390/vibration6030030
Submission received: 31 May 2023 / Revised: 26 June 2023 / Accepted: 27 June 2023 / Published: 29 June 2023

Abstract

:
In the design of lightweight structures, both the dynamics and durability must be taken into account. In this paper, a methodology for the combined optimization of structural dynamics, lightweight design, and lifetime with discrete vibration engineering measures is developed and discussed using a demonstration structure. A two-sided welded bending beam is excited at the centre and optimal parameters for tuned mass dampers (TMD) are searched, satisfying the requirements for the dynamic behaviour, the overall mass, and the lifetime of the weldings. It is shown that the combination of a reduced order model with the implementation of the structural stress approach at critical welds enables an efficient evaluation of certain design concepts in the time domain. Using this approach, multi-criterial optimization methods are used to identify the best set of parameters of the TMD to reduce the structural vibrations and enhance the durability.

Graphical Abstract

1. Introduction

In various industrial sectors, such as mechanical engineering, cranes, intra-logistics, automotive, and agricultural engineering, the requirements on the vibration behaviour of components and structures are constantly increasing since it affects performance, comfort, and noise emission. At the same time, products must be optimized in terms of lightweight design to save resources during production and operation while they still have to satisfy various other requirements. For example, the overall stiffness or the maximum allowable deflection needs to be considered or sufficient, static and cyclic strength needs to be verified, and in some cases also the crash-worthiness needs to be addressed.
The optimization of such dynamically loaded structures is usually decoupled according to the state of the art. However, such a recursive approach may lead to many design cycles. All mechanical properties, e.g., structural dynamics and fatigue, are closely linked and interdependent. Those relations can be seen as a typical multidisciplinary optimization (MDO) problem. MDO is evolving with the growing complexity of modern engineering problems. It involves different engineering domains to evaluate interactions within a coupled optimization problem. Traditional design and domain optimization mainly focus on optimal properties on component level and engineering sub-domains. As complexity increases due to progress, even more interactions need to be considered. In addition, lightweight structures are not just a hot topic, they are state-of-the-art. However, with the goal of achieving the lightest possible structures, dynamic problems are steadily increasing. This can correlate with reduced durability. Thus, structural dynamics should be addressed together with fatigue investigations to enhance structural durability and system reliability. MDO is well established, theoretically covered, and, for different types of problems, investigated, e.g., [1,2,3]. Structural optimization with fatigue constraints can also be found, e.g., [4,5]. However, just a few contributions are known which address structural dynamic fatigue correlation in terms of vibration engineering for lightweight structures and none of them deals with parameter optimization of vibrational measures.
Zeiler demonstrated in [6] the correlation between structural dynamics and fatigue performance in an automotive design optimization task using fatigue analysis and sensitivity derivatives in the context of a MDO. Häusler and Albers provided a workflow in [7] for optimization of the shape of a hole in a dynamically loaded plate and included a fatigue analysis. The process involves a multi-body simulation whose results were applied to a finite element model as excitation in the time domain. The location of the inner nodes of the hole were the optimization design points and could be moved according to the stress distribution. A commercial software was then used for damage estimation. Miao et al. [8] presented a multidisciplinary fatigue optimization method and evaluated them on structural lightweight design of a railway carbody, including a fatigue assessment. The multibody structure is first investigated numerically, starting with a multibody simulation and afterwards evaluating the stress in a FE environment. They also compared numerical and experimental results and concluded that MDO in terms of structural dynamics and fatigue can improve lightweight design. Han et al. [9] propose a methodology for structural topology optimization with stress minimization considering non-linear continuum damage. Chen et al. [10] presented topology optimization methods with the penalization of elemental damage to achieve lightweight design by simultaneously restricting the maximal damage. Faes et al. [11] show a reliability-based design optimization, taking into account uncertainties. Martins et al. [12] provide a survey on different methods in multidisciplinary design optimization and classified them. They indicated that benchmarking of different approaches is difficult—which leads to the approach of evolving a process by starting the optimization procedure with a low order model. Additionally, they stated that large-scale problems in MDO have only slow convergence. Here, reduced order models can provide a benefit. For obtaining less computational effort in multidisciplinary design optimization problems, Meng et al. [13] used a surrogate model for uncertainty evaluations in MDO. Their investigations were made on a turbine blade considering aerodynamic performance and structural reliability. They concluded that uncertainty investigation of this coupled problem can be investigated more efficient using an MDO approach. An interesting multidisciplinary design optimization approach for a gear transmission is presented by Mahiddini et al. [14] with an analytical formulation based on product reliability and economic aspects like product costs and customer preference. However, aspects of dynamic behaviour and lightweight design are not addressed in that work.
Moreover, all these works do not focus on the durability of welds especially. Therefore, the aim of this work is to present a method for evaluating and optimizing both the structural dynamics and the durability of a lightweight structure with critical welds in the same computational environment via the use of vibration measures.

2. Materials and Methods

2.1. Modelling of Dynamic System Behaviour

Modelling and simulation of dynamic systems are fundamental steps in a product development process. In the mathematical description of these systems, balance equations for energies, masses, and momentum can be derived with the help of basic physical laws. This usually leads to ordinary differential equations (ODE). The ODEs are usually generated by analytical or numerical methods and can be solved by various algorithms. The general approach to the computational solution of higher order ODEs involves the transformation to first order ODE systems. This increases the dimension of the system matrices in proportion to the number of derivatives.
In structural mechanics, the representation form of a second-order system has proven to be useful, since the assignment of the system matrices to the physical properties of the structure is unambiguous. The second-order system of equations can be represented as a linear descriptor system or generalized control system (1).
M u ¨ ( t ) + B u ˙ ( t ) + C u ( t ) = B i n f ( t )
y ( t ) = B 1 o u t u ( t ) + B 2 o u t u ˙ ( t )
Here,  t R +  is the time variable and the vector functions  u : R R n f : R R k , and  y : R R l  describe the internal states, the inputs, and the outputs of the system. The system order n represents the dimension of the system matrices. The number of inputs or outputs is denoted by k and l. In structural systems, the real matrices  M B C R n × n  represent the mass, damping, and stiffness matrices, respectively. For systems with regular mass matrix  M , Equation (1) corresponds to a second order ODE system.
In practice, ODE systems with a huge amount of degrees of freedom (DOF) often arise in the modelling of dynamic systems. In particular, this applies to numerical discretization methods, e.g., Finite Element Method (FEM), and reduces the attractiveness for optimization tasks or time domain simulations. Models with approximately the same accuracy and considerably greater computational efficiency are often generated by model order reduction (MOR) methods. The reduced order models (ROM) are then available for efficient analysis or optimization tasks. For further details about numerical modelling, especially FEM and the implementation, the reader is referred to the literature (e.g., [15,16]).

2.1.1. Model Order Reduction Methods

A large part of the currently established MOR techniques are projective methods, where the output system of order n is projected into a subspace of order r. The prerequisites for this are that the dynamics between input and output of the dynamical system are preserved in a defined time or frequency range and that the methods are numerically stable and efficient. Basically, the methods are distinguished with respect to the system behaviour to be considered. The majority of the developed algorithms assume linear and time-invariant systems.
In structural dynamics, the modal methods for MOR are among the most frequently used approaches, since they enable a suitable approximation of the dynamic behaviour. The basic idea of all modal methods is the transformation of physical coordinates into the modal coordinate space. For this purpose, the eigenvectors of the FE model are calculated by solving the eigenvalue problem and are combined in the modal matrix. The transformation is then based on this same modal matrix and the calculated deformations of the overall structure arise from the superposition of eigenvectors. This procedure is called modal synthesis or modal superposition [17]. For systems with linear independent eigenvectors and nonzero eigenvalues, the modal stiffness and modal mass matrices are diagonal. Thus, the equations of the system in modal space are decoupled, i.e., the transfer function of the overall system corresponds to the superposition of the transfer functions of decoupled single-mass oscillators.
In the modal coordinate reduction the modal superposition is combined with modal truncation. With the FE model, a reduced amount of eigenvectors is calculated and the modal matrix is thus built from a reduced modal base. Since a transformation matrix with a limited number of dynamic eigenmodes is used for this MOR, the static behaviour of the ROM may exhibit an error [18], which corresponds to the static residuals of the modes not considered.

2.1.2. Modal Approach

In modal truncation, the ODE system (1) is transformed into generalized or modal coordinates. The transformation rule is
u ( t ) = Φ w ( t )
with the modal matrix  Φ  and the vector of modal degrees of freedom  w ( t ) . The time dependent representation is omitted in the following in favour of better readability.
By solving the generalized eigenvalue problem
C Φ i = λ i M Φ i
with the system matrices from Equation (1) and the eigenvalues  λ i , the modal matrix  Φ = [ Φ 1 , , Φ n ]  is calculated. Using this matrix as a projection matrix, a system of n decoupled differential equations is obtained:
Φ T M Φ w ¨ + Φ T B Φ w ˙ + Φ T C Φ w = Φ T B i n f
If the eigenvectors are mass-normalized and the damping matrix is approximated by a linear combination of the stiffness and mass matrices, the following equation
I w ¨ + diag ( 2 ϑ ω 0 ) w ˙ + diag ( ω 0 2 ) w = Φ T B i n f
with the global damping  ϑ  and the natural frequencies  ω 0  applies. Equations (5) and (6) describe the full model of dimension n with correspondingly many eigenvectors.
Since the system matrices are now diagonal and contain parameters to n decoupled equations, rows and columns of the modal transformed matrices associated to an eigenvalue can be truncated. The reduced modal basis  Φ r R n × r  is described only by the first r eigenmodes. Using the reduced modal basis as transformation matrix, the physical degrees of freedom
u = Φ r w ¯ + ϵ
are approximated with an approximation error  ϵ  that depends on the selected number of eigenmodes r. Based on Equation (1), the system description is derived with  r n  decoupled differential equations
M ¯ w ¯ ¨ + B ¯ w ¯ ˙ + C ¯ w ¯ = Φ r T B i n f
y = B ¯ 1 o u t w + B ¯ 2 o u t w ˙
with the reduced system matrices
M ¯ = Φ r T M Φ r ,
B ¯ = Φ r T B Φ r ,
C ¯ = Φ r T C Φ r R r × r
and
w ¯ R r .

2.1.3. System Formulation

Since each modal deformation corresponds to a modal stress state [19], the output matrix  B ¯ 1 o u t  can also be extended to stresses. The idea is based on the fundamental equation
j = 1 n σ j w j = σ
with modal stresses  σ j  in the particular mode j and the total stress  σ . The modal coordinates  w j  correspond to a weighting factor describing the influence of an individual modal stress to the total stress. Compared to the modal reduced model (Section 2.1.2) the nodal stresses are used instead of eigenvectors to build the output matrix. For the back transformation from modal to physical stress values, the following transformation matrix is used.
Φ σ = [ σ 1 , , σ r ]
This results in the final system representation according to Equations (8) and (9) with the extended output matrix
B ¯ 1 o u t = Φ r Φ σ
in which the stress states are now calculated in addition to the deformation states.

2.2. Fatigue Strength Assessment

2.2.1. Linear-Elastic Approaches

In a fatigue assessment, linear elastically calculated stresses  σ S  in the structure are compared to endurable stresses  σ R  that depend on the material strength. For constant amplitude loading, the fatigue life N, representing the endurable number of cycles until a fracture is expected, can be calculated by the equation of the S–N curve
N = N R σ R σ S k
where k is the slope of the S–N curve and  N R  is the number of endurable cycles for the stress  S R .
In case of dynamically loaded structures, not a constant but a variable stress amplitude is expected. With an appropriate cycle-counting algorithm, e.g., such as described in ASTM E1049-85 [20], different stress cycles  n i  with individual amplitudes  σ a , i S  and mean values  σ m , i S  can be extracted. The damaging effect of these cycles be assessed by a linear damage accumulation
D = i D i = i n i / N i
in which the damaging effect of each counted cycle is evaluated against the endurable cycles  N i =  for a given stress amplitude and mean stress using Equation (17). The approach is visualized in Figure 1, where  Δ σ  is the stress used for the evaluation, e.g., noth stress  σ e , hot-spot stress  σ h s . According to IIW [21], the slope  k *  of the S–N curve in the very high cycle regime ( N > N k ) is decreased to:
k * = 2 k 1 .
In Figure 1 Δ σ  is the stress used for evaluation, e.g., hot-spot stress  Δ σ h s  or notch stress  Δ σ e . Further information on the assessment of variable amplitudes can be found in [22].
Furthermore, additional influencing factors, such as mean stresses, stress gradients, or roughness need to be taken into account [23]. These factors are used to decrease or increase the endurable fatigue strength by shifting the S–N curve up or down. The slope k of the S–N curve is typically set constant.

2.2.2. Application for Welded Joints

Various approaches are available for the fatigue assessment of structures with welded components [24]. The ones most often used are linear–elastic approaches, such as the nominal stress, the structural stress or the notch stress approach [21].
For a numerical assessment of a large and complex structure, the structural stress approach [25] is commonly used. A structural stress is used for the evaluation that is extracted by an extrapolation of the surface stress to the failure critical weld toe [21]. Stresses at two locations are used for the evaluation that have the distance of  0.4 · h  and  1.0 · h  from the weld toe, where h is the material thickness, as shown in Figure 2. An extrapolation to the weld toe
σ h s = 1.66 · σ ( x = 0.4 · h ) 0.67 · σ ( x = 1.0 · h )
leads to the so-called hot-spot stress  σ h s .
The derived acting structural stresses can directly be evaluated against the S–N curve.
For the structural stress approach an endurable stress range of  Δ σ h s R  = 100 MPa at  N R = 2 × 10 6  cycles can be used [21]. The slope for welded joints under stresses normal to the weld in as-welded condition depends on the sheet thickness. For thick ( h 7  mm) welded sheets a slope of  k = 3  and for thin sheets h > 7 mm a slope of  k = 5  is recommended [26]. The knee point of the S–N curve is typically expected to be located at  N k = 10 7  cycles.
Compared to other local approaches, such as the effective notch stress approach [21], the structural stress approach has the advantage that stresses at only two locations need to be evaluated. Of course, for real-world structures this evaluation has to be conducted for every relevant cross section. So, depending on the complexity, hundreds, or even thousands, of locations need to be considered.

2.3. Application Model

For the application of the proposed approach, a simple demonstrator is chosen. The geometry is inspired by a field sprayer which is typically used on tractors for, e.g., nutrient distribution, see Figure 3 and Figure 4. This low complexity model allows a perceivable interpretation of the results. For point masses were added to the structure in order to obtain low resonance frequencies while keeping the overall length at 1 m, considering a possible realization as lab demonstrator. The model is symmetric, while each side is distributed into three sections of different thicknesses. Construction steel material properties are considered ( E = 210   GPa  and  ν = 0.3 ). Additional specifications are noted in Table 1. The transition regions due to the thickness change represent the welds  w 1  and  w 2 , which are the focus of the following investigations. At this stage of process development, only two input locations for measures are considered per side. For the optimization of the dynamic behaviour of the structure, passive mass-stiffness-damping ( m T k t ϑ T ) absorber (tuned-mass-damper (TMD)) are included. Since the positions are fixed, the parameters of the TMD are of interest instead of the optimal positioning. The design space of each parameter is defined in Section 3.1.2.

2.3.1. Finite-Element Model

Figure 4 shows the FE model of the demonstrator. The model was built using the Software ANSYS [27]. The discretization of the model is based on a second-order approach with hexahedron elements. The specifications of the mesh at the welds are designed after recommendations in [21], as discussed in Section 2.2.2. The model consists of 4640 elements, 24,191 nodes, and 71,352 degrees of freedom. The outer and the inner sections are loaded with an additional mass which is equally distributed (weighted Multi-Point Constraint (MPC) interpolation (RBE3)) to all corresponding nodes of the area. The displacement of the middle section is constrained in the symmetry plane ( x = 0  and  y = 0 ), excluding the z-direction. No other displacement boundary conditions have been applied. As an excitation, a displacement base excitation  z i n ( t )  is chosen. It is considered that the fixing is achieved by a stiff spring support, thus a spring with a longitudinal stiffness of  k L = 1   ×   10 8   N / m  is implemented. The excitation is equal for all nodes of the lower surface of the middle area using a MPC resulting in one global excitation point, which is connected to the spring  k L .
At this point, no physical measurement data of a ground profile exists, so the input signal is synthetically constructed, as shown in Figure 5. The frequency content is dominated by low-frequency behaviour, representing, e.g., a farmland ground.

2.3.2. Verification

To verify that the results of the ROM are valid, a verification the FE solution is computed using ANSYS. For this aim, a full transient simulation is conducted, which represents the FOM. While the FOM is solved using an implicit scheme, the ROM is solved with an explicit fourth-order Runge–Kutta algorithm [28] using MATLAB/SIMULINK software environment [29]. For comparison,  1   s  of the shown base excitation signal is taken into account. The FE FOM simulation, as well as the ROM system simulation are performed with a time step size of  t s = 1   ×   10 4   s , which ensures numerically stable solutions covering the resonance of the stiff spring (stability condition only for ROM). The ROM is established after the description in Section 2.1.2. The used eigenvectors are obtained from a modal analysis using ANSYS. The results, here the modal eigenvectors (displacements) and the modal stresses, are extracted and organized as a state-space representation, as described in Section 2.1.3. This first order system of linear equations can be evaluated in the time domain. This is equal to a modal superposition transient analysis composed of modal states. If the displacements  u  of the ROM are comparable to the displacements of the FOM  u ref , then the stress states are also valid  σ σ ref . This can be derived from the usual isotropic stress–strain relationship following Equation (21). The stress tensor  σ  depends on the strain tensor  ε , as well as the material tensor  C , and is calculated by the Frobenius inner product.
σ = C , ε F
The strains are the derivatives of the calculated displacements. The fourth order material tensor depends on the material description. In case of linear–elastic, isotropic behaviour, it is symmetric with 36 components while 12 entries are non-zero elements which depend on only two of three parameters ( E , G , ν ), see, e.g., [30,31].
As a criterion for the amount of modes to consider, the effective mass is investigated [32]. The applied displacement excitation is fully one-dimensional, thus only the fraction in the direction of the excitation is relevant. It has been detected that the first 16 (modal) eigenvectors are necessary to achieve a valid ROM stress result compared to the FOM solution. This amount of eigenmodes mobilize 98.8% of the total mass in z-direction (cumulative mass fraction in z-direction). Taking more eigenmodes into account does not lead to any further relevant change—they only enlarge the model size. This amount can be further reduced by neglecting all modes in this range without any contribution to the effective masses in z-direction. The total number of modes then is reduced to 10. Since a modal approach is used, the ROM consists of the modal eigenvectors. These are computed using the FOM, so there is no mismatch between the eigenvectors for the points of interest.
To ensure that the ROM is able to approximate FOM results a verification stage is necessary. The easiest and computationally most efficient way is to compare frequency domain data. Figure 6 shows this for a normal stress transfer function at one exemplary welding position up to  1000   H z  for magnitude and phase response. It can be seen that both solutions are very well comparable, only slight differences occur in the higher frequency domain. However, the following investigations are made within the time domain, thus a time domain verification is also completed. Figure 7 shows the results of the FOM and the ROM for the maximum principle stress at the weld, denoting the transition from the middle to the inner section, excited by the base excitation input signal shown in Figure 5. The shown  1 s  cutout verifies the ROM solution since the characteristics, as well as the amplitudes, are very well comparable. The ROM can therefore be used in the optimization process, presented in the following chapter. Again, slight differences are visible but comparing the simulation times after Table 2, the advantage is predominant; the FOM takes roughly half an hour (1907 s) to solve a single design point, optimization procedures are not expedient—only the ROM enables this kind of procedure. The overall process then includes parameter optimization using the ROM and concludes with a final verification solution using the discovered optimal parameters with the FOM.

3. Results and Discussion

3.1. Optimization

The overall goal of the optimization process is to find a set of parameters that minimizes both the dynamic amplitudes and the damage values of the critical welds with a minimum expense of additional mass. Since these demands are in conflict with each other and cannot be minimized at the same time, the optimization problem can be seen as a MDO task from a mathematical point of view, which can be denoted as a multi-criterial optimization since several criteria are considered. However, to gain a good understanding of the problem, the system is initially analysed in the three dimensions structural dynamics, durability and lightweight design independently.

3.1.1. Mono-Criterial Optimization

Input for the damage accumulation are the counted hot-spot stress ranges  Δ σ h s , i  using ASTM 1049-85 based on the simulated stress-time course and the design S–N curve with  σ R  = 100 MPa and the slopes  k = 3  and  k * = 5 . The damage D is evaluated with the linear damage accumulation.
Four values are defined to asses each parameter set:  A 1  is the maximum displacement range (peak to peak) at the outer tip within the given time frame,  A 2  is the mean displacement at the same point, M is the overall mass and D is the damage value of the most critical weld, that is the maximum of all damage values  d i . The equations for each value is given in (22) to (25) with u denoting the difference in displacement in z-direction (Equation (26)).
A 1 = max ( u ( t ) ) min ( u ( t ) )
A 2 = mean ( | u ( t ) mean ( u ( t ) ) | )
M = m 0 + 2 m T
D = max ( d i )
u ( t ) = z o u t ( t ) z i n ( t )
The calculation of  A 1  and  A 2  according to Equations (22) and (23) is visualized in Figure 8 for exemplary time data of u and a simulation time of 5 s.
Since each value should be as low as possible, various optimization algorithms can be used to find optimal design parameters in a given design space. However, for a simple structure as this a full factorial parameter variation can be performed with reasonable computational effort. This typically gives a good understanding of the problem and allows an illustrative 2D visualization in the case of only two design variables. Therefore, the simulation model was run 21 × 21 = 441 times, while the absorber mass  m T  was varied from 0 to 0.300 kg in steps of 0.015 kg and the absorber stiffness  k T  was varied from 100 N/m to 900 N/m in steps of 40 N/m. The damping ratio  ϑ T  was kept constant at 3%, which is a realistic average value for vibration absorbers with elastomer-based damping.
Figure 9 shows the relative values for  A 1 A 2 , M, and D for the full factorial parameter variation (relative to the corresponding values of the initial system, denoted as  A 1 , o A 2 , o M o , and  D o ). The optimal value is marked in each case and given in Table 3, together with the corresponding design parameters mass  m T , stiffness  k T , and the related undamped eigenfrequency of the absorber  f T .
As observed, the single parameter optimization leads to very different optimal design parameters. The value  A 1  can be reduced most effectively with an absorber mass of 0.255 kg tuned to 8.85 Hz, while the mean displacement  A 2  is lowest for an absorber mass of 0.135 kg tuned to 10.5 Hz. The best lightweight design has no absorber at all, and for the most durable design an absorber mass of 0.060 kg tuned to 10.98 Hz is needed. Interestingly, the design parameter set for the best dynamics differs from the one for the best durability, which underlines that evaluating the damage of the structure together with its dynamic behaviour is essential.
With the results shown in Figure 9 the engineer can choose a proper set of design parameters, taking into account the specific boundary conditions and requirements of the application, and see, at a glance, the implications of a certain design regarding dynamics, durability, and lightweight design.

3.1.2. Multi-Criterial Optimization

For multi-criterial optimization problems in the engineering context many different techniques are known and described in the literature [33,34,35].
Here, the most straight forward approach with weighting functions for each target value and an arithmetic summation of all four components is chosen, which mathematically reduces the problem to the mono-criterial optimization of a cost function. The corresponding equations are given in Equations (27)–(30).
For the maximum displacement difference  A 1  a step function is chosen to guarantee that a certain threshold value  A 1 , m a x  is never exceeded, which may result from practical considerations in the application ((27), with H denoting the Heaviside step function). For the mean displacement  A 2 , a quadratic weighting function is chosen with the optimization parameter  ζ . This accounts for the assumption that the vibration of the system should be as low as possible at all times for a good performance (Equation (28)). For the overall mass M, a linear weighting function is chosen with the optimization parameter  ξ , assuming that any additional mass is seen as disadvantageous in the application (Equation (29)). For the damage D, again a step function is chosen with a threshold value  d m a x , taking into account that in many engineering problems a minimum lifetime must be guaranteed, but no further requirements regarding lifetime are made (Equation (29)). The formulation of these weighting functions and the specification of the parameters can be adjusted to each individual application scenario considering further requirements and boundary conditions of the system.
S A 1 = H ( A 1 A 1 , m a x )
S A 2 = ζ ( A 2 ) 2
S M = ξ ( M m 0 )
S D = H ( D d m a x )
The summation of the weighting functions results in the overall cost function S (31) that needs to be minimized.
S = S A 1 + S A 2 + S M + S D
The optimization parameters for this example are set to  A 1 , m a x = 30  mm,  ζ = 0.005  1/mm²,  ξ = 0.1  1/kg and  d m a x  = 0.001 as an initial guess. The determination of suitable values is not focused on here as the emphasis is on the methodology. Therefore, the values in this context were chosen exemplary. In practice, engineering expertise, experience, and economic considerations and design restrictions will be used in determining the values. The consideration of individual example cases or a pairwise comparison can also be helpful here. Since full factorial data are available, the goal function S can easily be visualized and the minimum value identified. The surface plot of S is shown in Figure 10. The white areas are invalid regions because the cost function is above 1. The parameters of the minimum are given in the first row of Table 4.
A standard optimization algorithm is used to find an optimum of S. In this case, the function fminsearch of Matlab is implemented, which is based on the simplex search method of Lagarias et al. [36]. Upper and lower limits are additionally implemented and the algorithm was confined to a maximum of 44 evaluations, that is 10% of the number of evaluations of the full factorial simulation. Using the central point of the design space as initial guess ( m T  = 0.150 kg,  k T  = 500 N/m), the algorithm yields a result that is worse compared to the minimum found by the full-factorial search. This is due to the fact that the function S is not very smooth and the algorithm is susceptible to become caught in local minima. Therefore, the minimum found by the full-factorial search was used as initial guess for another optimization run, yielding now an improved design point. The values of these two optimization runs are given in the second and third row of Table 4 and Table 5.
In the optimization runs, the damping ratio has been kept constant at 3%. Finally, an optimization run with the damping ratio as free design variable (within 0.5 and 6%, accounting for practical limitations) is performed, starting from the same design point as before. The results of this optimization process are given in the last row of Table 4 and Table 5.
Finally, all optimization results obtained in this study are visualized in Figure 11 in the three dimensional solution space regarding structural dynamics (in terms of mean displacement  A 2 ), lightweight design (in terms of overall mass M) and durability (in terms of damage D). On the left side, the conflict of goals between overall mass and displacement reduction can be seen. Up to an additional mass of 7% an increasing displacement reduction is possible with higher masses, clearly showing a typical Pareto frontier for these objectives with the optimum found by method I being part of the front line. The Pareto frontier is defined as the sum of all solutions for which any improvement in one objective can only take place if at least one other objective worsens [37]. The quick optimization run, starting from the centre of the design space (method II), does not reach the Pareto frontier. However, the successive optimization runs starting from the best full factorial value (method III and IV) exceed the Pareto frontier. On the right side, the relationship between displacement reduction and lifetime increase is comprehensible. By and large, a reduction in mean displacement corresponds with the reduction in the damage value; however, this is not true in general as, for some cases, a displacement reduction can also mean a reduction in durability. The dashed line indicates the threshold  d m a x .
The results of this simple example show the challenges of multi-criterial optimization. As it can be seen, the full factorial analysis gives a good understanding of the problem and is recommended if the required simulation capacity is available. A pure gradient-based optimization with an arbitrary starting point runs the risk of ending in an unfavourable local minimum due to the complexity of the cost function and should be used with care. A combination of a coarse full-factorial search with a limited number of design parameters and a subsequent gradient-based optimization starting from a previously found minimum can yield satisfactory results with justifiable numerical effort. However, more complex search algorithms that are designed to search for global minima like particle swarm, surrogate, or pattern search are also possible options [38]. As can be seen, numerical optimization processes need many evaluations and can rapidly become time-consuming for complicated structures in the time domain. In these cases, an evaluation of both the structural dynamics and the durability is recommended in the frequency domain, which is typically accompanied by a loss in accuracy in assessment [39], but which will lead to better optimization results within a limited computation time.

4. Conclusions

In this paper, a numerical approach for optimizing lightweight structures that are subjected to dynamic loads has been presented. This simulative approach takes into account both the structural dynamics and the lifetime assessment of the most critical welds. The focus has been on the method itself, which has been presented and discussed by means of a simple structure with parameterized absorbers. For this aim, we demonstrated how to establish a numerical valid reduced order model for stress evaluation in time from a full order model and what advantages can be achieved. Furthermore, a general optimization process was handled and specific results for mono-criteria and multi-criteria optimization approaches were presented.
From the investigations, the following main conclusions can be drawn:
  • Numerically efficient assessment of dynamically loaded structures with critical welds is possible using reduced order model methods and structural stress approach;
  • For a limited number of design variables, a full factorial simulation helps to understand the system and identify conflicting objectives between different criteria;
  • If a single cost function is defined combining the different objectives, a careful selection of the optimization algorithm and start values is recommended in order to avoid trusting optimization results arising from local minima;
  • Moreover, the robustness of the optimum should be analysed with respect to the optimization parameters. In case of doubt, a robust, less good design point is preferable to an optimal but very sensitive design point.
Current and future work is focused on applying this method to more practical structures and also to integrate techniques of parametric model order reduction in order to consider a variation of, e.g., sheet thicknesses. Moreover, more sophisticated multi-criterial optimization algorithms will be used which can also address and analyse the robustness of the identified optimal design points. Beyond this, a more holistic approach will be developed to simplify the initialization and optimization procedure for easy and safe industrial handling.

Author Contributions

Conceptualization, W.K.; durability methods, J.B., modelling, M.B.; methodology, C.T., M.B., J.B., W.K.; optimization W.K., writing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data will be available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations and symbols are used in this manuscript:
CADComputer-Aided Design
DOFDegree Of Freedom
FEFinite Element
FEMFinite Element Method
FOMFull Order Model
MDOMultidisciplinary Optimization
MORModel Order Reduction
MPCMulti-Point Constraint
ODEOrdinary Differential Equation
RBERigid Body Element
ROMReduced Order Model
TMDTuned Mass Damper
Nomenclature
  A 1 maximum displacement range
  A 2 mean displacement
Nnumber of cycles
Ddamage
  d m a x damage threshold
EYoung’s modulus
Gshear modulus
  ν Poisson’s ratio
kslope of S–N curve
hthickness of sheet metal
  σ maximum principal stress
Soverall cost function
  S A 1 structural dynamics cost function based on  A 1
  S A 2 structural dynamics cost function based on  A 2
  S M lightweight cost function
  S D durability cost function based on D
ttime
  t s time step
  T s i m simulation time
Mmass
Sdimensionless goal function
  m T absorber mass
  m 0 mass of the base structure
  k T absorber stiffness
  f T absorber eigenfrequency
  ϑ T damping ratio of absorber
  ζ optimization parameter in weighting function  S 2
  ξ optimization parameter in weighting function  S M
  z i n input displacement
  z o u t output displacement
urelative displacement
  w i welds
HHeaviside step function
  σ second order stress tensor
  ε second order strain tensor
  C fourth order material tensor
udisplacement
  u displacement vector
  w modal coordinate vector
  C stiffness matrix
  B damping matrix
  M mass matrix
  B o u t output matrix
  B i n input matrix
  Φ modal matrix
  Φ σ modal stress matrix

References

  1. Alexandrov, N.M.; Hussaini, M.Y. (Eds.) Multidisciplinary Design Optimization: State of the Art; SIAM Proceedings Series; SIAM: Philadelphia, PA, USA, 1997. [Google Scholar]
  2. Blachut, J. Emerging Methods for Multidisciplinary Optimization; CISM International Centre for Mechanical Sciences; Springer International Publishing: Cham, Switzerland, 2021; Volume 425. [Google Scholar]
  3. Salagame, R.R.; Ramu, P.; Narayanaswamy, I.; Saxena, D.K. Advances in Multidisciplinary Analysis and Optimization; Springer: Singapore, 2020. [Google Scholar] [CrossRef]
  4. El-Sayed, M.; Lund, E. Structural optimization with fatigue life constraints. Eng. Fract. Mech. 1990, 37, 1149–1156. [Google Scholar] [CrossRef]
  5. Zhao, L.; Xu, B.; Han, Y.; Xue, J.; Rong, J. Structural topological optimization with dynamic fatigue constraints subject to dynamic random loads. Eng. Struct. 2020, 205, 110089. [Google Scholar] [CrossRef]
  6. Zeiler, T. Use of structural dynamic and fatigue sensitivity derivatives in an automotive design optimization. Struct. Multidiscip. Optim. 2002, 23, 390–397. [Google Scholar] [CrossRef]
  7. Häußler, P.; Albers, A. Shape optimization of structural parts in dynamic mechanical systems based on fatigue calculations. Struct. Multidiscip. Optim. 2005, 29, 361–373. [Google Scholar] [CrossRef] [Green Version]
  8. Miao, B.; Luo, Y.; Peng, Q.; Qiu, Y.; Chen, H.; Yang, Z. Multidisciplinary design optimization of lightweight carbody for fatigue assessment. Mater. Des. 2020, 194, 108910. [Google Scholar] [CrossRef]
  9. Han, Y.; Xu, B.; Duan, Z.; Huang, X. Stress-based bi-directional evolutionary structural topology optimization considering nonlinear continuum damage. Comput. Methods Appl. Mech. Eng. 2022, 396, 115086. [Google Scholar] [CrossRef]
  10. Chen, Z.; Long, K.; Wen, P.; Nouman, S. Fatigue-resistance topology optimization of continuum structure by penalizing the cumulative fatigue damage. Adv. Eng. Softw. 2020, 150, 102924. [Google Scholar] [CrossRef]
  11. Faes, M.G.; Valdebenito, M.A. Fully decoupled reliability-based design optimization of structural systems subject to uncertain loads. Comput. Methods Appl. Mech. Eng. 2020, 371, 113313. [Google Scholar] [CrossRef]
  12. Martins, J.R.R.A.; Lambe, A.B. Multidisciplinary Design Optimization: A Survey of Architectures. AIAA J. 2013, 51, 2049–2075. [Google Scholar] [CrossRef] [Green Version]
  13. Meng, D.; Yang, S.; Zhang, Y.; Zhu, S.P. Structural reliability analysis and uncertainties-based collaborative design and optimization of turbine blades using surrogate model. Fatigue Fract. Eng. Mater. Struct. 2019, 42, 1219–1227. [Google Scholar] [CrossRef]
  14. Mahiddini, B.; Chettibi, T.; Benfriha, K.; Aoussat, A. Multidisciplinary design optimization of a gear train transmission. Concurr. Eng. 2019, 27, 268–281. [Google Scholar] [CrossRef]
  15. Zienkiewicz, O.C.; Zhu, J.; Taylor, R.L. The Finite Element Method: Its Basis and Fundamentals, 7th ed.; Elsevier: Amsterdam, The Netherlands; Boston, MA, USA; Heidelberg, Germany, 2013. [Google Scholar] [CrossRef]
  16. Pozrikidis, C. Introduction to Finite and Spectral Element Methods Using MATLAB®, 2nd ed.; CRC Press Taylor & Francis Group: Boca Raton, FL, USA, 2014. [Google Scholar]
  17. Breitbach, E. Modal Synthesis Modal Correction—Modal Coupling. In Identification of Vibrating Structures; Natke, H.G., Ed.; Springer: Vienna, Austria, 1982; pp. 321–348. [Google Scholar]
  18. Qu, Z.Q. Model Order Reduction Techniques: With Applications in Finite Element Analysis; Springer: London, UK, 2004. [Google Scholar]
  19. Pelayo, F.; Skafte, A.; Aenlle, M.L.; Brincker, R. Modal Analysis Based Stress Estimation for Structural Elements Subjected to Operational Dynamic Loadings. Exp. Mech. 2015, 55, 1791–1802. [Google Scholar] [CrossRef]
  20. ASTM E1049-85; Standard Practice for Cycle Counting in Fatigue Analysis. ASTM: West Conshohocken, PA, USA, 2005.
  21. Hobbacher, A.F. Recommendations for Fatigue Design of Welded Joints and Components; Springer International Publishing: Berlin/Heidelberg, Germany, 2016. [Google Scholar] [CrossRef]
  22. Baumgartner, J.; Waterkotte, R.; Hesseler, J. Fatigue assessment of a welded automotive differential under multiaxial and variable amplitude loading. Int. J. Fatigue 2021, 149, 106292. [Google Scholar] [CrossRef]
  23. Rennert, R.; Kullig, E.; Vormwald, M.; Esderts, A.; Siegele, D. Analytical strength assessment of components: FKM Guideline, 7th ed.; VDMA: Frankfurt, Germany, 2020. [Google Scholar]
  24. Radaj, D.; Sonsino, C.M.; Fricke, W. Fatigue Assessment of Welded Joints by Local Approaches; Woodhead Publishing Limited: Soston, UK, 2006. [Google Scholar] [CrossRef]
  25. Niemi, E.; Fricke, W.; Maddox, S.J. Structural Hot-Spot Stress Approach to Fatigue Analysis of Welded Components; Springer: Singapore, 2018. [Google Scholar] [CrossRef]
  26. Baumgartner, J.; Hobbacher, A.F.; Rennert, R. Fatigue assessment of welded thin sheets with the notch stress approach—Proposal for recommendations. Int. J. Fatigue 2020, 140, 105844. [Google Scholar] [CrossRef]
  27. ANSYS Inc. ANSYS Software Package Version 2021 R2; ANSYS Inc.: Canonsburg, PA, USA, 2021; Available online: https://www.ansys.com (accessed on 26 June 2023).
  28. Hoffman, J.; Frankel, S. Numerical Methods for Engineers and Scientists; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  29. The MathWorks Inc. Software, MATLAB Version: 9.9.0.1718557 (R2020b); The MathWorks Inc.: Natick, MA, USA, 2020; Available online: https://www.mathworks.com (accessed on 26 June 2023).
  30. Zienkiewicz, O.C. The Finite Element Method in Engineering Science, 2nd ed.; McGraw-Hill: London, UK, 1971. [Google Scholar]
  31. Haupt, P. Continuum Mechanics and Theory of Materials, 2nd ed.; Advanced Texts in Physics; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  32. Clough, R.; Penzien, J. Dynamics of Structures; International Student Edition; McGraw-Hill: New York, NY, USA, 1975. [Google Scholar]
  33. Gunantara, N. A review of multi-objective optimization: Methods and its applications. Cogent Eng. 2018, 5, 1502242. [Google Scholar] [CrossRef]
  34. Marler, R.T.; Arora, J.S. Survey of multi-objective optimization methods for engineering. Struct. Multidiscip. Optim. 2004, 26, 369–395. [Google Scholar] [CrossRef]
  35. Mirjalili, S.; Jangir, P.; Saremi, S. Multi-objective ant lion optimizer: A multi-objective optimization algorithm for solving engineering problems. Appl. Intell. 2017, 46, 79–95. [Google Scholar] [CrossRef]
  36. Lagarias, J.C.; Reeds, J.A.; Wright, M.H.; Wright, P.E. Convergence properties of the Nelder–Mead simplex method in low dimensions. SIAM J. Optim. 1998, 9, 112–147. [Google Scholar] [CrossRef] [Green Version]
  37. Sawaragi, Y.; Nakayama, H.; Tanino, T. Theory of Multiobjective Optimization; Elsevier: Amsterdam, The Netherlands, 1985. [Google Scholar]
  38. Koessler, E.; Almomani, A. Hybrid particle swarm optimization and pattern search algorithm. Optim. Eng. 2021, 22, 1539–1555. [Google Scholar] [CrossRef]
  39. Muñiz-Calvente, M.; Álvarez-Vázquez, A.; Pelayo, F.; Aenlle, M.; García-Fernández, N.; Lamela-Rey, M.J. A comparative review of time- and frequency-domain methods for fatigue damage assessment. Int. J. Fatigue 2022, 163, 107069. [Google Scholar] [CrossRef]
Figure 1. Design S–N curve, example stress spectrum and damage accumulation procedure graphically exemplified for the load level  i = 3 .
Figure 1. Design S–N curve, example stress spectrum and damage accumulation procedure graphically exemplified for the load level  i = 3 .
Vibration 06 00030 g001
Figure 2. Determination of structural stresses at a weld.
Figure 2. Determination of structural stresses at a weld.
Vibration 06 00030 g002
Figure 3. Sketch of the demonstration model with boundary conditions and defined designations. The length of the model is 1 m.
Figure 3. Sketch of the demonstration model with boundary conditions and defined designations. The length of the model is 1 m.
Vibration 06 00030 g003
Figure 4. Demonstrator model: CAD representation and FE model with zoom to welding  w 1 . Symmetrical welding joints  w 1 , 2  and connection points for tuned mass dampers I, II.
Figure 4. Demonstrator model: CAD representation and FE model with zoom to welding  w 1 . Symmetrical welding joints  w 1 , 2  and connection points for tuned mass dampers I, II.
Vibration 06 00030 g004
Figure 5. Synthetic base excitation signal  z i n ( t ) .
Figure 5. Synthetic base excitation signal  z i n ( t ) .
Vibration 06 00030 g005
Figure 6. Normal stress  σ x x  transfer function comparison for magnitude and phase response of the FOM and the ROM using unit excitation for verification.
Figure 6. Normal stress  σ x x  transfer function comparison for magnitude and phase response of the FOM and the ROM using unit excitation for verification.
Vibration 06 00030 g006
Figure 7. (Left) Comparison of the maximum principal stress solution of the inner node for the calculation at welding  w 1 . (Right) Detail view of left image. Excitation with signal shown in Figure 5.
Figure 7. (Left) Comparison of the maximum principal stress solution of the inner node for the calculation at welding  w 1 . (Right) Detail view of left image. Excitation with signal shown in Figure 5.
Vibration 06 00030 g007
Figure 8. Visualization of the computation of  A 1  and  A 2  for an exemplary time data of  u ( t ) .
Figure 8. Visualization of the computation of  A 1  and  A 2  for an exemplary time data of  u ( t ) .
Vibration 06 00030 g008
Figure 9. Full factorial parameter variation for mono-criterial optimization with two variable parameters and 21 × 21 evaluations. The white circles indicate the respective minima.
Figure 9. Full factorial parameter variation for mono-criterial optimization with two variable parameters and 21 × 21 evaluations. The white circles indicate the respective minima.
Vibration 06 00030 g009
Figure 10. Surface plot of objective function S based on full factorial parameter variation with two variable parameters and 21 × 21 evaluations. The white circles indicate the minima found with method I and III.
Figure 10. Surface plot of objective function S based on full factorial parameter variation with two variable parameters and 21 × 21 evaluations. The white circles indicate the minima found with method I and III.
Vibration 06 00030 g010
Figure 11. Visualization of all optimization results in the three dimensional solution space. The legend applies to both graphs. The arrows indicate the optimization runs for method II, III, and IV.
Figure 11. Visualization of all optimization results in the three dimensional solution space. The legend applies to both graphs. The arrows indicate the optimization runs for method II, III, and IV.
Vibration 06 00030 g011
Table 1. Additional model specifications.
Table 1. Additional model specifications.
SectionLengthWidthThicknessAdd. MassBase Stiff.
mmmkgN/m
Outer0.200.072 × 10−30.5-
Inner0.250.074 × 10−30.5-
Mid0.100.0710 × 10−3-1 × 10−8
Table 2. Comparison of computing times for simulation of  T s i m = 1   s  with  t s = 1 × 10 4   s . FOM solved with ANSYS. ROM solved with MATLAB/SIMULINK. ANSYS used 16 cores in distributed memory parallel mode. MATLAB used four cores in shared memory parallel mode. Before the ROM solution, a modal analysis must be performed once with ANSYS, which takes  9   s . The solvable system formulation must also be established once and takes  23   s . The required time per iteration is the last denoted time specification ( + 2   s ). In addition, the computing time for  T s i m = 30   s  using the ROM (=ROM  30   s ) is noted in the last column to underline the ROM advantage.
Table 2. Comparison of computing times for simulation of  T s i m = 1   s  with  t s = 1 × 10 4   s . FOM solved with ANSYS. ROM solved with MATLAB/SIMULINK. ANSYS used 16 cores in distributed memory parallel mode. MATLAB used four cores in shared memory parallel mode. Before the ROM solution, a modal analysis must be performed once with ANSYS, which takes  9   s . The solvable system formulation must also be established once and takes  23   s . The required time per iteration is the last denoted time specification ( + 2   s ). In addition, the computing time for  T s i m = 30   s  using the ROM (=ROM  30   s ) is noted in the last column to underline the ROM advantage.
TypeFOM 1 sROM ANSYS 1 sROM 1 sROM 30 s
Time   1907 s   812 s   9 s + 23 s + 2 s   9 s + 23 s + 19 s
Table 3. Results of the mono-criterial optimization gained by full factorial calculation (21 × 21 evaluations).
Table 3. Results of the mono-criterial optimization gained by full factorial calculation (21 × 21 evaluations).
Goal   m T   k T   f T   A 1 / A 1 , o   A 2 / A 2 , o   M / M o   D / D o
kgN/mHz
  A 1 0.2557808.850.42970.41761.12750.3688
  A 2 0.13558010.500.52240.39591.06750.2930
M0.000--1.00001.00001.00001.0000
D0.06030010.980.62420.43401.03000.2843
Table 4. Results of the multi-criterial optimization based on the cost function S with different optimization methods: I: full factorial search (2 variable parameters), II: fminsearch (2 variable parameters, start values  m T  = 0.150 kg,  k T  = 500 N/m), III: fminsearch (2 variable parameters, start values  m T  = 0.06 kg,  k T  = 340 N/m), IV: fminsearch (3 variable parameters, start values  m T  = 0.06 kg,  k T  = 340 N/m,  ϑ T  = 3%).
Table 4. Results of the multi-criterial optimization based on the cost function S with different optimization methods: I: full factorial search (2 variable parameters), II: fminsearch (2 variable parameters, start values  m T  = 0.150 kg,  k T  = 500 N/m), III: fminsearch (2 variable parameters, start values  m T  = 0.06 kg,  k T  = 340 N/m), IV: fminsearch (3 variable parameters, start values  m T  = 0.06 kg,  k T  = 340 N/m,  ϑ T  = 3%).
Method   m T   k T   ϑ T   f T SNumb.
kgN/m%Hz-Eval.
I0.0600340.03.011.980.03443441
II0.0976412.03.010.340.0426044
III0.0448239.23.011.630.03096485
IV0.0581300.45.911.440.02582485
Table 5. Mono-criterial values of the results of the multi-criterial optimization based on the cost function S with different optimization methods: I—full factorial search (2 variable parameters); II—fminsearch (2 variable parameters, start values  m T  = 0.150 kg,  k T  = 500 N/m); III—fminsearch (2 variable parameters, start values  m T  = 0.06 kg,  k T  = 340 N/m); and IV—fminsearch (3 variable parameters, start values  m T  = 0.06 kg,  k T  = 340 N/m,  ϑ T  = 3%).
Table 5. Mono-criterial values of the results of the multi-criterial optimization based on the cost function S with different optimization methods: I—full factorial search (2 variable parameters); II—fminsearch (2 variable parameters, start values  m T  = 0.150 kg,  k T  = 500 N/m); III—fminsearch (2 variable parameters, start values  m T  = 0.06 kg,  k T  = 340 N/m); and IV—fminsearch (3 variable parameters, start values  m T  = 0.06 kg,  k T  = 340 N/m,  ϑ T  = 3%).
Method   A 1 / A 1 , o   A 2 / A 2 , o   M / M o   D / D o
I0.63090.42671.03000.2925
II0.54240.43291.04880.3009
III0.55680.42271.02240.2666
IV0.46700.33981.02910.2090
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kaal, W.; Baumgartner, J.; Budnik, M.; Tamm, C. Numerical Approach to Optimize the Dynamic Behaviour of Structures Considering Structural Durability. Vibration 2023, 6, 477-493. https://doi.org/10.3390/vibration6030030

AMA Style

Kaal W, Baumgartner J, Budnik M, Tamm C. Numerical Approach to Optimize the Dynamic Behaviour of Structures Considering Structural Durability. Vibration. 2023; 6(3):477-493. https://doi.org/10.3390/vibration6030030

Chicago/Turabian Style

Kaal, William, Jörg Baumgartner, Maximilian Budnik, and Christoph Tamm. 2023. "Numerical Approach to Optimize the Dynamic Behaviour of Structures Considering Structural Durability" Vibration 6, no. 3: 477-493. https://doi.org/10.3390/vibration6030030

APA Style

Kaal, W., Baumgartner, J., Budnik, M., & Tamm, C. (2023). Numerical Approach to Optimize the Dynamic Behaviour of Structures Considering Structural Durability. Vibration, 6(3), 477-493. https://doi.org/10.3390/vibration6030030

Article Metrics

Back to TopTop