Microstructural Topology Optimization of Constrained Layer Damping on Plates for Maximum Modal Loss Factor of Macrostructures

Fang, Zhanpeng; Yao, Lei; Tian, Shuxia; Hou, Junjian

doi:https://doi.org/10.1155/2020/8837610

Shock and Vibration

On this page

Abstract Introduction Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Nonlinear Vibration of Continuous Systems

View this Special Issue

Research Article | Open Access

Volume 2020 | Article ID 8837610 | https://doi.org/10.1155/2020/8837610

Microstructural Topology Optimization of Constrained Layer Damping on Plates for Maximum Modal Loss Factor of Macrostructures

Zhanpeng Fang,¹Lei Yao,¹Shuxia Tian,¹and Junjian Hou¹

Academic Editor: Konstantin V. Avramov

Received29 Apr 2020

Revised10 Sept 2020

Accepted24 Sept 2020

Published13 Oct 2020

Abstract

This paper presents microstructural topology optimization of viscoelastic materials for the plates with constrained layer damping (CLD) treatments. The design objective is to maximize modal loss factor of macrostructures, which is obtained by using the Modal Strain Energy (MSE) method. The microstructure of the viscoelastic damping layer is composed of 3D periodic unit cells. The effective elastic properties of the unit cell are obtained through the strain energy-based method. The density-based topology optimization is adopted to find optimal microstructures of viscoelastic materials. The design sensitivities of modal loss factor with respect to the design variables are analyzed and the design variables are updated by Method of Moving Asymptotes (MMA). Numerical examples are given to demonstrate the validity of the proposed optimization method. The effectiveness of the optimal design method is illustrated by comparing a solid and an optimized cellular viscoelastic material as applied to the plates with CLD treatments.

1. Introduction

Viscoelastic damping materials are desired in many engineering applications to reduce unwanted noise and vibration due to their favourable characteristics in dissipating dynamic energy. Damping treatment approaches including active, passive, and hybrid passive-active methods [1–3] have been developed to reduce noise and vibration and applied to many industrial fields. The constrained layer damping (CLD) treatments provide considerably more damping effect than unconstrained ones despite the addition of the additional weight of the constraining layer, and they are widely applied to attenuate the structural vibration. In order to design lightweight structures with high damping performance, the optimization for damping treatment has been extensively studied. Plunkett [4] proposed a method of cutting the constraining layer into appropriate lengths to increase the structure damping of the CLD structure. Ray and Baz [5] studied the optimization of energy dissipation of active constrained layer damping treatments of plates. Nakra [6] investigated how to design the layout of CLD treatments to improve the vibration performance. Zheng et al. [7] proposed the layout optimization of the CLD treatments to minimize vibration energy and sound radiation of cylindrical shells.

Topology optimization method is regarded as a powerful tool for seeking innovative structural design. Bendsøe and Kikuchi [8] firstly proposed the topology optimization method by employing homogenization theory. With the development of the topology optimization method, several well-known topology optimization methods include the solid isotropic material with penalization (SIMP) method [9], level set method (LSM) [10], the Evolutionary Structural Optimization (ESO) [11], and the bidirectional evolutionary structural optimization (BESO) [12]. A series of methods for structural topology optimization were proposed to address the layout problem of the damping layer in the past decades. Kim et al. [13] compared the modal loss factors of the structures with damping material treatment obtained by topology optimization to the ones obtained by other approaches and further demonstrated that topology optimization is the most effective way to find optimal damping layout in a shell structure. James and Waisman [14] proposed a time-dependent adjoint sensitivity analysis method to produce optimal designs of structures that exhibit viscoelastic creep deformation. Madeira et al. [15] optimized the distribution of CLD material on a laminated composite panel to minimize weight and maximize modal damping simultaneously.

The damping performance of CLD treatments mainly depends on the physical and geometric properties of the viscoelastic damping material. Therefore, there is a great desire to design the microstructures of viscoelastic materials with the demand property. Sigmund [16, 17] firstly proposed an inverse homogenization method to design materials with prescribed properties. Yi et al. [18] utilized the inverse homogenization method to optimize the microstructure of the damping material with the aim of improving the damping characteristics of viscoelastic composites. Huang et al. [19] utilized the BESO method to design the composite microstructure with the desirable viscoelastic properties. Liu et al. [20] proposed a topology optimization algorithm based on the BESO method to enhance the macroscopic modal damping and natural frequency of structures constructed by optimized viscoelastic materials.

Topology optimization problems involve a large number of design variables. Hence, gradient-based methods are preferred over evolutionary- and population-based approaches. It has been shown that, based on current computing technology, the evolutionary methods are prohibitively time- and resource-consuming for most reasonably sized topology optimization problems [21]. The Method of Moving Asymptotes (MMA) [22] is a well-known gradient-based algorithm, and it is widely used in topology optimization of structures. Zheng et al. [23] used topology optimization with MMA approach to optimally design damping treatment by maximizing damping effect. Kang et al. [24] investigated the optimal distribution of damping material in vibrating structures which subject to harmonic excitations by using a topology optimization method. El-Sabbagh and Baz [25] conducted topology optimization for unconstrained damping treatment of plates. Takezawa et al. [26] proposed complex dynamic compliance as objective function for optimizing damping materials to reduce the resonance peak response. Fang et al. [27] proposed an efficient optimization procedure integrating the Pseudoexcitation Method (PEM) and the double complex modal superposition method to optimize the layout of the CLD structures subjected to stationary random excitation. Chen and Liu [28] investigated the optimal microstructural configuration of the viscoelastic material to improve the damping characteristics of the macrostructures. Yun and Youn [29] optimized the microstructures of viscoelastic materials to effectively attenuate transient responses of the viscoelastic damping structures. Zhao et al. [30] proposed an efficient concurrent topology optimization approach for minimizing the maximum dynamic response of two-scale hierarchical structures in the time domain.

Most of the present works concerning microstructural topology optimization of damping material has concentrated on the Sandwich beams. Microstructural topology optimization of damping material of a plate with CLD treatment is relatively limited. Therefore, the purpose of this paper is to optimize microstructure of viscoelastic material of the CLD plates with the aim of maximizing the modal loss factor of the macrostructure. The microstructure of the viscoelastic damping layer is represented by 3D periodical unit cell (PUC), and its effective shear modulus is obtained by using the strain energy-based method. The design sensitivities of modal loss factor with respect to the design variables are analyzed and the optimization problem is solved by the MMA approach.

2. Optimization Problem and Material Interpolation Scheme

A two-scale topology optimization problem of the CLD plate is studied in this paper. The macrostructure of the CLD plate is composed of the uniform cellular viscoelastic material. The microstructure of the viscoelastic damping layer is represented by 3D periodical unit cell. The optimization objective is to design the optimal layout of the viscoelastic material within the 3D periodical unit cell with the prescribed volume fraction so that the optimal CLD plate has a maximum modal loss factor for damping the vibration. Therefore, the optimization problem can be mathematically expressed as follows:where , are the design variables and is the total number of the elements within the 3D periodical unit cell. denotes that element is full of viscoelastic damping material and represents a void element. is the rth mode loss factor. is the volume of the ith element and is the prescribed volume fraction.

The viscoelastic damping material of an element in the 3D periodical unit cell can be treated to be isotropic and its physical properties are assumed to be a function of the elemental density, . To obtain relative clear 0-1 topological structure, the Solid Isotropic Material with Penalization (SIMP) model [31] with an exponential ‘power-law’ scheme is adopted to penalize the intermediate densities and the material interpolation scheme can be expressed as follows:where is the material loss factor and . and denote the density and the storage modulus of the viscoelastic damping material when . and are the penalty factors.

In the SIMP interpolation scheme, the reasonable penalty factor should be chosen to suppress the intermediate density elements in the final optimal topology. A small penalty factor will lead to a large number of intermediate density elements in the final optimal results as the punishment is insufficient in the process of solving. However, when the penalty parameter is too large, many high-relative-density elements (for example, greater than 0.5 and less than 1.0) may be irrational to delete. It will lead to checkerboard pattern numerical instability and even cause the incorrect optimization results [32]. The combination of and is always used to optimize modal loss factor and vibration response [26–28, 33, 34]. Therefore, the penalty factors and are 1 and 3 in this paper.

3. Finite Element Analysis and Sensitivity Analysis

3.1. Two-Scale Structure Analysis

The finite element analysis should be conducted both for the CLD plate at the macroscale and the 3D periodical unit cell at the microscale. The CLD structure uses its damping mechanism based on vibratory energy dissipation through transverse shear strains induced in the viscoelastic layer, so viscoelastic layer carries only transverse shear and no normal stress when the finite element model of the CLD plate is established. In the macroscale, the viscoelastic layer is defined as an orthotropic material. Thus, the global stiffness matrix and mass matrix can be calculated by [23]

In these equations, , , , , , and are the eth element stiffness and mass matrix of the base plate, viscoelastic layer, and constrained layer, respectively. and are the strain matrix and the shape function matrix, respectively. , , , and are constitutive matrix and mass density of the base plate and constrained layer, respectively. is the average mass density of the viscoelastic layer. and are the effective shear modulus of the viscoelastic layer. and are the shape function matrix of the shear strains and in the viscoelastic layer, respectively.

The global stiffness matrix contains the real part and the imaginary part , which are defined as follows:where “Re” denotes “Real part of” and “Im” represents “Imaginary part of.”

The governing equation of the CLD structure for free vibration is written as follows:where is the displacement vector in the macrostructure.

According to the Modal Strain Energy Method, the modal loss factor at the rth mode can be found as follows:where is the eigenvector, which is obtained by the real modal analysis.

The effective elastic properties of the viscoelastic layer directly depend on its microstructure. In micromechanics, there are two basic theories: the homogenization theory [28–30] and the strain energy-based method [35–37]. The basic feature of the strain energy-based method is that both the unit cell of the periodic structure and the corresponding unit volume of the homogeneous solid undergo the same strain energy. Compared with the homogenization theory, the strain energy-based method can provide a relatively simple and efficient way in the estimation of the effective elastic properties and the sensitivity analysis [37].

As illustrated in Figure 1, the 3D unit cell is a cube, which does not change during the optimization process. The topology of viscoelastic damping material in a 3D unit cell is changed by varying the design variables . The 3D unit cell of the viscoelastic layer is defined as an orthotropic material. The generalized Hooke’s law can be written as follows:where , , and () are the effective stress tensor, strain tensor, and elastic tensor.

The effective shear modulus of the unit cell and can be determined by applying the following two unit strain fields:

The uniform strain boundary condition is replaced with equivalent Dirichlet boundary condition in numerical implementation and the corresponding boundary conditions are resumed in Tables 1 and 2 [37]. In the tables, , , and are the sizes of the 3D unit cell in the X, Y, and Z directions, respectively.

The strain energy of the orthotropic material under the unit strain can be calculated bywhere is the volume of the unit cell.

The total strain energy can also be obtained using the domain integration of the local strain energy density in the unit cell as follows:

Under the boundary condition described in Table 1, the finite element equation governing the unit cell can be expressed as follows:where and are the equivalent nodal displacement vector and force vector caused by the unit strain . is the real part of the global stiffness matrix of the unit cell, which is expressed as follows:where is the strain matrix of the ith element in the unit cell. is the real part of the stiffness matrix of ith element.

In finite element analysis, equation (12) can be rewritten as follows:where represents the equivalent nodal displacement vector of the ith element.

Keeping in mind that , the effective shear modulus of the unit cell can be expressed as follows:

Similarly, the effective shear modulus of the unit cell is expressed as follows:where is the equivalent nodal displacement vector of the ith element, which can be obtained through the statics analysis of the unit cell under the boundary condition described in Table 2.

The constitutive matrix of an orthotropic material has nine independent components, see equation (9). The procedure for the analysis of the effective shear modulus and is shown in the paper, and the procedure of the identification of the other seven components can be obtained from [37].

The effective density of the unit cell is evaluated through the following relationship:

3.2. Sensitivity Analysis

An important step in finding the optimum distribution of the unit cell is to carry out sensitivity analysis in which the effect of each design variable on the modal loss factor is estimated. The sensitivity of the modal loss factor with respect to the design variables can be formulated as follows:

Substituting equations (5) and (6) to equation (19), equation (20) can be obtained:

Considering the definition of matrix in equation (3), the sensitivities of and can be formulated as follows:

By analyzing the first partial derivative of equations (16) and (17) with respect to the design variables, the sensitivities of the effective shear modulus can be obtained:

3.3. Optimization Procedure

The entire topology optimization procedure for the design of the viscoelastic material microstructure is illustrated in Figure 2. Based on the initial microstructure, the effective shear modulus and density are computed by the homogenization process in equations (16), (17), and (18). The modal loss factor of the macrostructure and the corresponding sensitivities with respect to the design variables are calculated by a two-scale analysis in equations (8) and (20). Then, the MMA algorithm is utilized for the optimizer to update the new design variables. This process is repeated until a certain convergence criterion is satisfied. The convergence criterion is that the maximum difference of the design variables of two adjacent iteration steps is less than the given small value . The convergence parameter is .

The MMA algorithm was proposed to adjust the curvature of the convex linearization method. Giving the current design , the MMA approximation of the maximum modal loss factor problem in equation (1) yields to the following linear programming problem:

The lower and upper asymptotes and are iteratively updated to mitigate oscillation or improve convergence rate. The heuristic rule is as follows [22].

For and ,

For ,where

4. Numerical Examples

Two numerical examples are presented in this section to illustrate the microstructural design of the CLD plates and demonstrate the effectiveness of the proposed optimization approach. The first is a rectangular plate with two short edges clamped and the second is a cantilever rectangular plate. It is assumed that the material properties of both cases are the same, as shown in Table 3.

It is assumed that the 3D unit cell characteristic size is much smaller than the characteristic size of the macrostructure. Such hypotheses support the homogenization theory that is used here to compute equivalent mechanical properties of the periodic viscoelastic damping materials. Then, the homogenized material properties of viscoelastic damping materials are independent of the size of the 3D unit cell [20]. For simplicity, nondimensional sizes 1 × 1 × 1 for the 3D unit cell are assumed. The 3D unit cell of the microstructure is discretized with eight-node hexahedron elements. The initial designs are of a uniform distribution with the given volume fraction except the elements at the center. The density of the elements at the center is set to be 1.

4.1. The Plate Clamped with Two Short Edges

The plate clamped with two short edges is shown in Figure 3. The length and width of the cantilever rectangular plate are 0.4 m and 0.2 m, respectively. The thickness of the base plate, viscoelastic layer, and constrained layer are 0.002 m, 0.0005 m, and 0.0005 m, respectively.

Topology optimization of the microstructure is carried out first for four different meshes of the 3D unit cell: 15 × 15 × 15, 20 × 20 × 20, 25 × 25 × 25, and 30 × 30 × 30. The optimization objective is to maximize the modal loss factor of the first mode and the volume fraction is 0.3. All of the results are shown in Table 4. It can be clearly seen that the optimum objective values of the final topologies using different meshes are close to each other. With the increase of the number of elements in the 3D unit cell, the optimized configurations are more detailed, but the computation time is significantly increased. The number of elements in a 3D unit cell can be determined by comprehensively considering the detail of the optimized configurations and the computational cost. The 3D unit cell in the following examples is divided into 20 × 20 × 20. From Table 4, it can also be seen that the optimal configurations have very few intermediate density elements. It means that while the design variables can be fractions, most design variables are close to 0 or 1.

The design results of maximizing the modal loss factor of the second mode and the modal loss factor of the third mode are shown in Table 5. It is found that the optimal microstructure of maximizing the modal loss factor of the second mode is obviously different from those of maximizing the modal loss factor of the first mode and the third mode. That is because the first mode and the third mode are bending modes and the shear strains occur in the plane y-z, while the second mode is torsional mode and the shear strains occur in the planes y-z and x-z.

For each mode, there is a certain value of the shear modulus that maximizes the modal loss factor in the damping structure. The volume fraction varies with the optimization objective. Three different design objectives are considered. They are (1) maximizing the modal loss factor of the first mode, (2) maximizing the modal loss factor of the second mode, and (3) maximizing the modal loss factor of the third mode. The design results are shown in Table 6. It can be seen that the iterations are smooth. For the single mode design, the modal loss factors of optimal designs are increased by 48.28%, 14.72%, and 6.56% than the traditional design (the solid viscoelastic material) at modes 1, 2, and 3, respectively.

In order to further verify the effectiveness of the optimal designs, the CLD plate is excited by a central time-harmonic force with the amplitude F = 1N. The location of the excitation force is shown in Figure 3 and the response is evaluated at the same point. The amplitude of frequency response function of the traditional design and the optimal designs are compared in Figure 4. It can be seen that when maximizing the first modal loss factor (Case 1), the amplitude of frequency response function at the first mode is minimum. Case 2 and Case 3 show the same tendency with Case 1.

4.2. The Cantilever Rectangular Plate

In order to further verify the effectiveness of the proposed optimization algorithm, another structure with different boundary conditions and size is discussed. A cantilever rectangular plate is shown in Figure 5. The length and width of the plate are 0.2 m and 0.1 m, respectively. Other physical and geometrical parameters are the same as the first example. The volume fraction is 0.3. Three different design objectives are considered.

The design results are shown in Table 7. Similar to the first example, the optimal microstructure and the effective shear modulus of the first mode are similar to those of the third mode. The optimal microstructure and the effective shear modulus of the second mode are obviously different from those of the first mode. That is also because the first mode and the third mode are bending modes and the shear strains occur in the plane y-z, while the second mode is the torsional mode and the shear strains occur in the planes y-z and x-z.

The volume fraction varies with the optimization objective. The design results are shown in Table 8. For the single mode design, the modal loss factors of optimal designs are increased by 160.32%, 1.55%, and 4.67% than the traditional design at modes 1, 2, and 3, respectively.

The loading condition of frequency response analysis is similar with the first example and the location of the frequency response point is shown in Figure 5. The amplitude of frequency response function of the traditional design and the optimal designs are compared in Figure 6. For mode 1, the resonance peak of the optimal design at the first mode is obviously lower than the traditional design. For mode 2 and mode 3, the resonance peaks of the optimal design for the corresponding mode are slightly higher than the traditional design. This is because the stiffness of the optimal design is naturally lower than the traditional design [29].

5. Conclusions

A microstructural topology optimization of viscoelastic materials of CLD plates is presented for maximizing the modal loss factor of the macrostructure. The strain energy-based method is employed to calculate effective properties of the viscoelastic materials composed of 3D periodic unit cells. Based on the homogenized properties, the modal loss factor of the macrostructure is performed by the Modal Strain Energy Method. The design sensitivities of the modal loss factor with respect to the design variables are analyzed. The density method is adopted to find the optimum layout on microscale of the viscoelastic damping layer. The numerical examples demonstrated the effectiveness of the proposed method. The lightweight damped plates with high damping capabilities can be designed through the proposed method.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This project was supported by the Natural Science Foundation of China (Grant no. 51805491), Plan of Key Research Projects of Higher Education of Henan Province (Grant no. 18A460035), and Doctoral Research Fund Project of Zhengzhou University of Light Industry (Grant no. 2015BSJJ031).

References

A. Baz, “Dynamic boundary control of beams using active constrained layer damping,” Mechanical Systems and Signal Processing, vol. 11, no. 6, pp. 811–825, 1997.
View at: Publisher Site | Google Scholar
A. L. Araújo, C. M. Mota Soares, and C. A. M. Soares, “A viscoelastic sandwich finite element model for the analysis of passive, active and hybrid structures,” Applied Composite Materials, vol. 17, no. 5, pp. 529–542, 2010.
View at: Publisher Site | Google Scholar
A. L. Araújo, P. Martins, C. M. Mota Soares, C. A. Mota Soares, and J. Herskovits, “Damping optimisation of hybrid active-passive sandwich composite structures,” Advances in Engineering Software, vol. 46, no. 1, pp. 69–74, 2012.
View at: Publisher Site | Google Scholar
R. Plunkett and C. T. Lee, “Length optimization for constrained viscoelastic layer damping,” The Journal of the Acoustical Society of America, vol. 48, no. 1B, pp. 150–161, 1970.
View at: Publisher Site | Google Scholar
M. C. Ray and A. Baz, “Optimization of energy dissipation of active constrained layer damping treatments of plates,” Journal of Sound and Vibration, vol. 208, no. 3, pp. 391–406, 1997.
View at: Publisher Site | Google Scholar
B. C. Nakra, “Vibration control in machines and structures using viscoelastic damping,” Journal of Sound and Vibration, vol. 211, no. 3, pp. 449–466, 1998.
View at: Publisher Site | Google Scholar
H. Zheng, C. Cai, G. S. H. Pau et al., “Minimizing vibration response of cylindrical shells through layout optimization of passive constrained layer damping treatments,” Journal of Sound and Vibration, vol. 279, no. 3–5, pp. 739–756, 2005.
View at: Publisher Site | Google Scholar
M. P. Bendsøe and N. Kikuchi, “Generating optimal topologies in structural design using a homogenization method,” Computer Methods in Applied Mechanics and Engineering, vol. 71, no. 2, pp. 197–224, 1988.
View at: Google Scholar
M. P. Bendsoe, “Optimal shape design as a material distribution problem,” Structural Optimization, vol. 1, pp. 193–202, 1989.
View at: Google Scholar
M. Y. Wang, X. Wang, and D. Guo, “A level set method for structural topology optimization,” Computer Methods in Applied Mechanics and Engineering, vol. 192, no. 1-2, pp. 227–246, 2003.
View at: Publisher Site | Google Scholar
Y. M. Xie and G. P. Steven, “A simple evolutionary procedure for structural optimization,” Computers & Structures, vol. 49, no. 5, pp. 885–896, 1993.
View at: Publisher Site | Google Scholar
X. Y. Yang, Y. M. Xie, G. P. Steven, and O. M. Querin, “Bidirectional evolutionary method for stiffness optimization,” AIAA Journal, vol. 37, no. 11, pp. 1483–1488, 1999.
View at: Publisher Site | Google Scholar
S. Y. Kim, C. K. Mechefske, and I. Y. Kim, “Optimal damping layout in a shell structure using topology optimization,” Journal of Sound and Vibration, vol. 332, no. 12, pp. 2873–2883, 2013.
View at: Publisher Site | Google Scholar
K. A. James and H. Waisman, “Topology optimization of viscoelastic structures using a time-dependent adjoint method,” Computer Methods in Applied Mechanics and Engineering, vol. 285, pp. 166–187, 2015.
View at: Publisher Site | Google Scholar
J. F. A. Madeira, A. L. Araújo, and C. M. Mota Soares, “Multiobjective optimization of constrained layer damping treatments in composite plate structures,” Mechanics of Advanced Materials and Structures, vol. 24, no. 5, pp. 427–436, 2017.
View at: Publisher Site | Google Scholar
O. Sigmund, “Materials with prescribed constitutive parameters: an inverse homogenization problem,” International Journal of Solids and Structures, vol. 31, no. 17, pp. 2313–2329, 1994.
View at: Publisher Site | Google Scholar
O. Sigmund, “Tailoring materials with prescribed elastic properties,” Mechanics of Materials, vol. 20, no. 4, pp. 351–368, 1995.
View at: Publisher Site | Google Scholar
Y.-M. Yi, S.-H. Park, and S.-K. Youn, “Design of microstructures of viscoelastic composites for optimal damping characteristics,” International Journal of Solids and Structures, vol. 37, no. 35, pp. 4791–4810, 2000.
View at: Publisher Site | Google Scholar
X. Huang, S. Zhou, G. Sun, G. Li, and Y. M. Xie, “Topology optimization for microstructures of viscoelastic composite materials,” Computer Methods in Applied Mechanics and Engineering, vol. 283, no. 253, pp. 503–516, 2015.
View at: Publisher Site | Google Scholar
Q. Liu, D. Ruan, and X. Huang, “Topology optimization of viscoelastic materials on damping and frequency of macrostructures,” Computer Methods in Applied Mechanics and Engineering, vol. 337, pp. 305–323, 2018.
View at: Publisher Site | Google Scholar
O. Sigmund, “On the usefulness of non-gradient approaches in topology optimization,” Structural and Multidisciplinary Optimization, vol. 43, no. 5, pp. 589–596, 2011.
View at: Publisher Site | Google Scholar
K. Svanberg, “The method of moving asymptotes-a new method for structural optimization,” International Journal for Numerical Methods in Engineering, vol. 24, no. 2, pp. 359–373, 1987.
View at: Publisher Site | Google Scholar
L. Zheng, R. Xie, Y. Wang, and E.-S. Adel, “Topology optimization of constrained layer damping on plates using method of moving asymptote (MMA) approach,” Shock and Vibration, vol. 18, no. 1-2, pp. 221–244, 2011.
View at: Google Scholar
Z. Kang, X. Zhang, S. Jiang, and G. Cheng, “On topology optimization of damping layer in shell structures under harmonic excitations,” Structural and Multidisciplinary Optimization, vol. 46, no. 1, pp. 51–67, 2012.
View at: Publisher Site | Google Scholar
A. El-Sabbagh and A. Baz, “Topology optimization of unconstrained damping treatments for plates,” Engineering Optimization, vol. 46, no. 9, pp. 1153–1168, 2014.
View at: Publisher Site | Google Scholar
A. Takezawa, M. Daifuku, Y. Nakano, K. Nakagawa, T. Yamamoto, and M. Kitamura, “Topology optimization of damping material for reducing resonance response based on complex dynamic compliance,” Journal of Sound and Vibration, vol. 365, pp. 230–243, 2016.
View at: Publisher Site | Google Scholar
Z. Fang, J. Hou, and H. Zhai, “Topology optimization of constrained layer damping structures subjected to stationary random excitation,” Shock and Vibration, vol. 2018, Article ID 7849153, 9 pages, 2018.
View at: Publisher Site | Google Scholar
W. Chen and S. Liu, “Microstructural topology optimization of viscoelastic materials for maximum modal loss factor of macrostructures,” Structural and Multidisciplinary Optimization, vol. 53, no. 1, pp. 1–14, 2016.
View at: Publisher Site | Google Scholar
K.-S. Yun and S.-K. Youn, “Microstructural topology optimization of viscoelastic materials of damped structures subjected to dynamic loads,” International Journal of Solids and Structures, vol. 147, pp. 67–79, 2018.
View at: Publisher Site | Google Scholar
J. Zhao, H. Yoon, and B. D. Youn, “Concurrent topology optimization with uniform microstructure for minimizing dynamic response in the time domain,” Computers & Structures, vol. 222, pp. 98–117, 2019.
View at: Publisher Site | Google Scholar
M. P. Bendsøe and O. Sigmund, “Material interpolation schemes in topology optimization,” Archive of Applied Mechanics, vol. 69, pp. 635–654, 1999.
View at: Google Scholar
Y. Du, S. Yan, Y. Zhang, H. Xie, and Q. Tian, “A modified interpolation approach for topology optimization,” Acta Mechanica Solida Sinica, vol. 28, no. 4, pp. 420–430, 2015.
View at: Publisher Site | Google Scholar
Y. Kun, C. G. Dong, and W. B. Ping, “Topology optimization of damping layers in shell structures subject to impact loads for minimum residual vibration,” Journal of Sound and Vibration, vol. 431, pp. 226–247, 2018.
View at: Google Scholar
H. Zhang, X. Ding, H. Li, and M. Xiong, “Multi-scale structural topology optimization of free-layer damping structures with damping composite materials,” Composite Structures, vol. 212, pp. 609–624, 2019.
View at: Publisher Site | Google Scholar
E. J. Barbero, Finite Element Analysis of Composite Material, Taylor and Francis Group, Milton Park, UK, 2008.
W. Zhang, G. Dai, F. Wang, S. Sun, and H. Bassir, “Using strain energy-based prediction of effective elastic properties in topology optimization of material microstructures,” Acta Mechanica Sinica, vol. 23, no. 1, pp. 77–89, 2007.
View at: Publisher Site | Google Scholar
A. Catapano and M. Montemurro, “A multi-scale approach for the optimum design of sandwich plates with honeycomb core. Part I: homogenisation of core properties,” Composite Structures, vol. 118, pp. 664–676, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Zhanpeng Fang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

435

Downloads

786

Citations