# A Simple Estimation of Coupling Loss Factors for Two Flexible Subsystems Connected via Discrete Interfaces

Authors: Jun Zhang; Lin Ji; Zhenyu Huang; Pingping Zhang; Wei Wang
Journal: Shock and Vibration
 (2016)

---

## Abstract

A simple formula is proposed to estimate the Statistical Energy Analysis (SEA) coupling loss factors (CLFs) for two flexible subsystems connected via discrete interfaces. First, the dynamic interactions between two discretely connected subsystems are described as a set of intermodal coupling stiffness terms. It is then found that if both subsystems are of high modal density and meanwhile the interface points all act independently, the intermodal dynamic couplings become dominated by only those between different subsystem mode sets. If ensemble- and frequency-averaged, the intermodal coupling stiffness terms can simply reduce to a function of the characteristic dynamic properties of each subsystem and the subsystem mass, as well as the number of interface points. The results can thus be accommodated within the theoretical frame of conventional SEA theory to yield a simple CLF formula. Meanwhile, the approach allows the weak coupling region between the two SEA subsystems to be distinguished simply and explicitly. The consistency and difference of the present technique with and from the traditional wave-based SEA solutions are discussed. Finally, numerical examples are given to illustrate the good performance of the present technique.


---

## Body

## 1. Introduction

Statistical Energy Analysis (SEA) [1] is commonly used for higher frequency modelling of complex systems. In SEA, the whole structure is firstly partitioned into many individual SEA subsystems with joined interfaces, and then the power balance equation is set up for each subsystem by equating the power input into the subsystem to that of output, including both the power dissipated within the subsystem itself by the subsystem damping and that transmitted to other subsystems via the subsystem interfaces. A key assumption which underpins the set-up of these power balance equations is that the transmitted power between any two SEA subsystems is proportional to the energy difference between them, and the proportionality is determined by the so-called “coupling loss factors (CLFs)” between the two subsystems [1]. Although a rigorous proof is still under investigations [2, 3], it has been widely accepted that the coupling power proportional (CPP) assumption is generally valid if all the SEA subsystems are highly modal dense and weakly coupled [4–9]. It is seen therefore that a successful application of the SEA theory somehow largely relies on the following two aspects: (1) a clear distinguishing of the weak coupling region and (2) an accurate estimation of the coupling loss factors between coupled SEA subsystems. These actually have formed the two main issues of the SEA modelling techniques in the past decades and resulted in extensive research efforts from the SEA researchers, such as [10–16]. A good summary of the large quantity of the relevant publications can be found in [17].Regarding simple estimations of CLFs for discretely connected subsystems, the existing SEA literature shows that simple formulas of CLF estimations are only available for limited coupling cases so far, for example, plate to room/cavity, room to room/cavity, plate to plate with a line connection, and structure to structure with a single-point connection [1, 18]. For multipoint couplings between SEA subsystems, for example, the solar panels in a spacecraft [19], however, there is no such simple formula available. It is true that [1] did propose two simple formulas for CLF estimations between two subsystems connected via a single point by using wave approach (10.1.6) and then by a modal approach (10.1.18). But, as stated in [1] by Lyon and DeJong, the derivations were only rigorous for one-dimensional subsystems and incoherence waves. It thus becomes rather arguable whether the two formulas can be extended straightforward to two- and three-dimensional subsystems with multipoint couplings in that more implicit/explicit assumptions have to be involved.For multipoint couplings, on one hand, one can generally rely on the principle of power injection methods (PIMs) [20] to predict the SEA CLFs, either experimentally [21] or numerically [22, 23], despite such a fact that PIMs are generally expensive and also lack of the ability to measure the coupling strengths between SEA subsystems. On the other hand, many specialist methods are also available which were mainly proposed from finite elements using modal bases of uncoupled subsystems, for example [12–15]. These specialist approaches are found to be much more accurate and efficient than PIMs. However, these methods still need a certain level of extra knowledge of the modes of the uncoupled subsystems except for the gross modal properties. It is therefore very desirable if SEA CLFs can be determined only from the gross modal properties of SEA subsystems and the number of interface degrees of freedom, and, preferably, the predicting technique also allows the weak coupling region to be distinguished in a simple and explicit manner.With this in mind, the paper aims to propose a simple formula to estimate the SEA CLFs for multipoint coupling cases together with a simple measurement of the coupling strength level between SEA subsystems. The derivation procedure is closely in line with the theoretical frame of [1], but the intermodal stiffness coupling terms are employed in the present derivation instead of the wave transmission coefficient terms in [1]. By assuming each subsystem being of high modal density and the interface degrees of freedom act independently from each other, the statistical coupling relations between the two sets of subsystem modes can be put into a very simple form. It thus allows a simple CLF formula to be derived by accommodating with the results given in [1, 2]. In addition, the predicting procedure allows the weak coupling region between the two subsystems to be determined simply but explicitly. Compared with other forms of CLF solutions, the proposed technique is found to be extremely easy to apply.The paper is organized as follows. First, a subsystem modal analysis technique, as described in [24], is used to find out the theoretical relations (both analytical and statistical) of the intermodal coupling stiffness between two generic SEA subsystems with discrete interfaces. Then, in Section 3, the statistics of the intermodal coupling stiffness terms are accommodated within the theoretical frame of [1, 2] to yield two simple formulas: one for SEA CLF estimation and the other for coupling strength measurement. Finally, numerical examples of two plates connected via discrete points are provided in Section 4 to illustrate the good performance of the proposed technique.


## 2. Dynamic Stiffness Couplings between Two SEA Subsystems with Discrete Connections

In this section, a system model consisting of two generic SEA subsystems with discrete interfaces is set up, as shown in Figure1. The intermodal dynamic relations between the two sets of subsystem modes can be found by a subsystem modal analysis technique described in [24]. Similar approaches can also be found in [6, 25]. The main derivation procedure is briefly given as shown below.Figure 1
Two generic SEA subsystems connected via discrete couplings.


### 2.1. Intermodal Coupling Stiffness Terms for General Interfaces

First, describe each subsystem by the modes when the subsystem is uncoupled from the rest. For subsystema, for instance, let ϕja, ωja, and ηja represent the jth mass-normalized mode shape, the corresponding natural frequency, and modal loss factor, respectively. The modal receptance matrix of subsystem a, Ya, which is diagonal, can then be expressed as(1)Yjja=1ωja21+iηja-ω2.In (1) i2=-1. Clearly there exists Yjja=qja/fja, where qja andfja are the jth modal amplitude and modal force of the subsystem, respectively. For subsystem b, similar expressions can be established by replacing the superscript a with b and subscript j with k.Assume that the interface forces and displacements can be decomposed along the interface regionσI in terms of a set of interface basis functions ψσI, that is, (2)FIσI=∑pfI,pψpσI,(3)WIσI=∑pqI,pψpσI,where p is the order of the basis function, and (4)∫σIψpσIψp′σIdσI=δpp′where δ is Dirac delta function. Such an interface decomposition technique has been found to be particularly useful when dealing with continuous interfaces [25].If the modes in subsystemsa and b as well as the interface modes can be truncated into finite sets, by making use of the modal orthogonal properties, both the interface modal force, fI,p, and the interface modal amplitude, qI,p, can be expressed in terms of the sets of modal forces, fIa and fIb, as well as the sets of modal amplitudes qIa and qIb. By enforcing force equilibrium and displacement continuity boundary conditions along the interface between subsystems a and b, the modal amplitudes of each subsystem can finally be calculated as [24](5)qaqb=GaaGabGabTGbbfeafeb.In (5), the superscript T represents transpose and qa and fea are both column vectors whose jth elements qja and fe,ja represent, respectively, the jth modal amplitude and modal force of subsystem a corresponding to the external excitation Fea at σea, of which (6)fe,ja=∫σaϕjaσeaFeaσeadσea.Similar expressions can also be established for subsystemb by replacing the superscript a with b, and subscript j with k, as appropriate.G a a and Gab are, respectively, the Green function matrices which give the modal amplitudes of subsystem a upon per unit modal forces applied on the subsystem itself and the modal amplitudes of subsystem a upon per unit modal forces applied on subsystem b. These Green function matrices, after some algebraic, can be determined by(7)Gaa=YaI-αA+B-1αTYa,(8)Gbb=YbI-βA+B-1βTYb,(9)Gab=YaαA+B-1βTYb,where α and β are matrices whose j,pth and k,pth elements are given by(10)αjp=∫σIϕjaσIψpσIdσI,βkp=∫σIϕkbσIψpσIdσIwhile A and B are square matrices, given by(11)A=αTYaα,(12)B=βTYbβ.Physically,A (11) and B (12) represent, respectively, the receptance matrices of subsystem a and b at their interfaces. A+B-1 is therefore the so-called “interface dynamic stiffness matrix” of subsystems a and b, determined by the subsystem modal properties as well as their local mode shapes at the interfaces.From (5), a dynamic stiffness matrix can be obtained as (13)Daa-Dab-DabTDbb=GaaGabGabTGbb-1,where(14)Daa=Gaa-GabGbb-1GabT-1,(15)Dbb=Gbb-GabTGaa-1Gab-1,(16)Dab=DaaGabGbb-1.Physically, thej,kth entry in Dab (16) represents the coupling stiffness between the jth mode of subsystem a and the kth mode of subsystem b, while the j,j′th entry in Daa (14) represents that between the jth and j′th modes of subsystem a. For convenience, here, these intermodal couplings are termed, respectively, as “out-of-set” and “in-set” mode couplings [2]. More discussions on the distinctions of the out-of-set and in-set mode coupling terms can be found in Section 2.2.1.


### 2.2. Intermodal Coupling Stiffness Terms for Discrete Interfaces

Assume the discrete interfaces composed ofN coupling points. By the criterion of (4), the interface basis functions in (2) here can then be chosen as a unit diagonal matrix with a size of N. (Note that, in general, there should be many choices of the sets of basis functions provided they are complete and orthogonal and can represent the interface deformations between the two subsystems well [25]; usually the set of local modal shapes defined on the interface DOFs can be a better choice of the basis functions [24]). Consequently, (10) simply yields(17)αjp=ϕjaσI,p,p=1,…,N,βkp=ϕkbσI,p,p=1,…,N.Substituting (17) into (11) and (12), the p,p′th elements of matrices A and B can thus be obtained, respectively, as(18)App′=∑jϕjaσI,pϕjaσI,p′ωja21+iηja-ω2,Bpp′=∑kϕkbσI,pϕkbσI,p′ωkb21+iηkb-ω2.It is seen that App′ is exactly the receptance of subsystem a between its connecting points σI,p and σI,p′, that is, App′=Ypp′a. Similarly, Bpp′=Ypp′b. Therefore, A and B are actually the receptance matrices of subsystems a and b at their interfaces.If it is assumed that the interface degrees of freedom act independently so that the coherences among the interface points are ignorable, then by (18), A and B tend to be diagonally dominated so that (19)A+Bpp≈Yppa+Yppb.Note that (19) is usually true if both subsystems have high modal densities and the interface points are spaced at least one wavelength apart from each other [26]. As a result, Gaa and Gbb, by (7) and (8), become diagonally dominated as well with(20)Gjjaa≈Yjja1-∑p1Yppa+YppbϕjaσI,pϕjaσI,pωja21+iηja-ω2,Gkkbb≈Ykkb1-∑p1Yppa+YppbϕkbσI,pϕkbσI,pωkb21+iηkb-ω2.For SEA subsystems with high modal densities, the second terms on the brackets of the right sides in (20) tend to be much less than unity. It thus implies that the Green functions of each SEA subsystem can only be slightly affected by the other through the couplings. In which case, (9) leads to(21)Gjkab≈Yjja∑pϕjaσI,pϕkbσI,pA+BppYkkb.Substituting (20)-(21) into (14)–(16), the dynamic coupling stiffness between any two pair of subsystem modes can be simply determined as(22)Djkab≈∑pϕjaσI,pϕkbσI,pA+Bpp.Considering the diagonal properties of matrices A and B, (22) then becomes(23)Djkab≈∑pϕjaσI,pKppabϕkbσI,p,where,(24)Kppab=Yppa+Yppb-1≈Yin,∞a+Yin,∞b-1where Yin,∞a and Yin,∞b are the characteristic receptances of the two subsystems when the subsystems are extended to infinity. Equation (23) can then be further reduced to (25)Djkab≈Yin,∞a+Yin,∞b-1∑pϕjaσI,pϕkbσI,p.Equation (25) indicates that the intermodal coupling stiffness term Djkab is a frequency-dependent function, determined by the characteristic receptances of the two subsystems Yin,∞a and Yin,∞b as well as the mode shapes of the two subsystems ϕja and ϕkb at their interfaces σI. Moreover, as the mean of Djkab is mainly determined by Yin,∞a+Yin,∞b, the variance of Djkab is mainly determined by those of ϕja and ϕkb at the interfaces. In the contents below, (25) will be firstly used in a coupling case with single-point connection in order to demonstrate its application.


#### 2.2.1. Application to a Single-Point Coupling Case

For a single-point coupling between subsystemsa and b, the interface function can be simply defined as ψpσI=1 with p=1. A and B (see (11) and (12)) are simply the input-receptances of the two subsystems at their connecting points σIa and σIb, as(26)A=YaσIa=∑jϕja2σIaYjja,B=YbσIb=∑kϕkb2σIbYkkb.It thus gives (by (24))(27)Kab=1YaσIa+YbσIb≈1Yin,∞a+Yin,∞b.Then, (7)–(9) become(28)Gaa=YaI-KabΦaΦaTYa,(29)Gbb=YbI-KabΦbΦbTYb,(30)Gab=KabYaΦaΦbTYb.Equations (28)–(30) can then be combined with (13) to determine both the in-set (Daa and Dbb) and out-of-set (Dab) intermodal coupling stiffness terms.In the sense of weak coupling, that is, the second terms in brackets of the right sides of (28) and (29) are both negligible [14], the coupling stiffness between the jth mode of subsystem a and the kth mode of subsystem b, by (16), thus gives(31)Djkab=ϕjaσIaKabϕkbσIb.In [2], the intermodal coupling stiffness, for two subsystems connected by a single spring K, is given by(32)Djkab≈ϕjaxIaKϕkbxIb.It is seen that (31) and (32) are in a very similar form, except that Kab is a frequency-dependent function in the former equation while K is a constant in the latter one. It thus allows the results in [2] for spring-coupled subsystems to be extended for the coupling cases with rigid-point connections in a way fairly straightforward. These are described in the section below.


## 2.1. Intermodal Coupling Stiffness Terms for General Interfaces

First, describe each subsystem by the modes when the subsystem is uncoupled from the rest. For subsystema, for instance, let ϕja, ωja, and ηja represent the jth mass-normalized mode shape, the corresponding natural frequency, and modal loss factor, respectively. The modal receptance matrix of subsystem a, Ya, which is diagonal, can then be expressed as(1)Yjja=1ωja21+iηja-ω2.In (1) i2=-1. Clearly there exists Yjja=qja/fja, where qja andfja are the jth modal amplitude and modal force of the subsystem, respectively. For subsystem b, similar expressions can be established by replacing the superscript a with b and subscript j with k.Assume that the interface forces and displacements can be decomposed along the interface regionσI in terms of a set of interface basis functions ψσI, that is, (2)FIσI=∑pfI,pψpσI,(3)WIσI=∑pqI,pψpσI,where p is the order of the basis function, and (4)∫σIψpσIψp′σIdσI=δpp′where δ is Dirac delta function. Such an interface decomposition technique has been found to be particularly useful when dealing with continuous interfaces [25].If the modes in subsystemsa and b as well as the interface modes can be truncated into finite sets, by making use of the modal orthogonal properties, both the interface modal force, fI,p, and the interface modal amplitude, qI,p, can be expressed in terms of the sets of modal forces, fIa and fIb, as well as the sets of modal amplitudes qIa and qIb. By enforcing force equilibrium and displacement continuity boundary conditions along the interface between subsystems a and b, the modal amplitudes of each subsystem can finally be calculated as [24](5)qaqb=GaaGabGabTGbbfeafeb.In (5), the superscript T represents transpose and qa and fea are both column vectors whose jth elements qja and fe,ja represent, respectively, the jth modal amplitude and modal force of subsystem a corresponding to the external excitation Fea at σea, of which (6)fe,ja=∫σaϕjaσeaFeaσeadσea.Similar expressions can also be established for subsystemb by replacing the superscript a with b, and subscript j with k, as appropriate.G a a and Gab are, respectively, the Green function matrices which give the modal amplitudes of subsystem a upon per unit modal forces applied on the subsystem itself and the modal amplitudes of subsystem a upon per unit modal forces applied on subsystem b. These Green function matrices, after some algebraic, can be determined by(7)Gaa=YaI-αA+B-1αTYa,(8)Gbb=YbI-βA+B-1βTYb,(9)Gab=YaαA+B-1βTYb,where α and β are matrices whose j,pth and k,pth elements are given by(10)αjp=∫σIϕjaσIψpσIdσI,βkp=∫σIϕkbσIψpσIdσIwhile A and B are square matrices, given by(11)A=αTYaα,(12)B=βTYbβ.Physically,A (11) and B (12) represent, respectively, the receptance matrices of subsystem a and b at their interfaces. A+B-1 is therefore the so-called “interface dynamic stiffness matrix” of subsystems a and b, determined by the subsystem modal properties as well as their local mode shapes at the interfaces.From (5), a dynamic stiffness matrix can be obtained as (13)Daa-Dab-DabTDbb=GaaGabGabTGbb-1,where(14)Daa=Gaa-GabGbb-1GabT-1,(15)Dbb=Gbb-GabTGaa-1Gab-1,(16)Dab=DaaGabGbb-1.Physically, thej,kth entry in Dab (16) represents the coupling stiffness between the jth mode of subsystem a and the kth mode of subsystem b, while the j,j′th entry in Daa (14) represents that between the jth and j′th modes of subsystem a. For convenience, here, these intermodal couplings are termed, respectively, as “out-of-set” and “in-set” mode couplings [2]. More discussions on the distinctions of the out-of-set and in-set mode coupling terms can be found in Section 2.2.1.


## 2.2. Intermodal Coupling Stiffness Terms for Discrete Interfaces

Assume the discrete interfaces composed ofN coupling points. By the criterion of (4), the interface basis functions in (2) here can then be chosen as a unit diagonal matrix with a size of N. (Note that, in general, there should be many choices of the sets of basis functions provided they are complete and orthogonal and can represent the interface deformations between the two subsystems well [25]; usually the set of local modal shapes defined on the interface DOFs can be a better choice of the basis functions [24]). Consequently, (10) simply yields(17)αjp=ϕjaσI,p,p=1,…,N,βkp=ϕkbσI,p,p=1,…,N.Substituting (17) into (11) and (12), the p,p′th elements of matrices A and B can thus be obtained, respectively, as(18)App′=∑jϕjaσI,pϕjaσI,p′ωja21+iηja-ω2,Bpp′=∑kϕkbσI,pϕkbσI,p′ωkb21+iηkb-ω2.It is seen that App′ is exactly the receptance of subsystem a between its connecting points σI,p and σI,p′, that is, App′=Ypp′a. Similarly, Bpp′=Ypp′b. Therefore, A and B are actually the receptance matrices of subsystems a and b at their interfaces.If it is assumed that the interface degrees of freedom act independently so that the coherences among the interface points are ignorable, then by (18), A and B tend to be diagonally dominated so that (19)A+Bpp≈Yppa+Yppb.Note that (19) is usually true if both subsystems have high modal densities and the interface points are spaced at least one wavelength apart from each other [26]. As a result, Gaa and Gbb, by (7) and (8), become diagonally dominated as well with(20)Gjjaa≈Yjja1-∑p1Yppa+YppbϕjaσI,pϕjaσI,pωja21+iηja-ω2,Gkkbb≈Ykkb1-∑p1Yppa+YppbϕkbσI,pϕkbσI,pωkb21+iηkb-ω2.For SEA subsystems with high modal densities, the second terms on the brackets of the right sides in (20) tend to be much less than unity. It thus implies that the Green functions of each SEA subsystem can only be slightly affected by the other through the couplings. In which case, (9) leads to(21)Gjkab≈Yjja∑pϕjaσI,pϕkbσI,pA+BppYkkb.Substituting (20)-(21) into (14)–(16), the dynamic coupling stiffness between any two pair of subsystem modes can be simply determined as(22)Djkab≈∑pϕjaσI,pϕkbσI,pA+Bpp.Considering the diagonal properties of matrices A and B, (22) then becomes(23)Djkab≈∑pϕjaσI,pKppabϕkbσI,p,where,(24)Kppab=Yppa+Yppb-1≈Yin,∞a+Yin,∞b-1where Yin,∞a and Yin,∞b are the characteristic receptances of the two subsystems when the subsystems are extended to infinity. Equation (23) can then be further reduced to (25)Djkab≈Yin,∞a+Yin,∞b-1∑pϕjaσI,pϕkbσI,p.Equation (25) indicates that the intermodal coupling stiffness term Djkab is a frequency-dependent function, determined by the characteristic receptances of the two subsystems Yin,∞a and Yin,∞b as well as the mode shapes of the two subsystems ϕja and ϕkb at their interfaces σI. Moreover, as the mean of Djkab is mainly determined by Yin,∞a+Yin,∞b, the variance of Djkab is mainly determined by those of ϕja and ϕkb at the interfaces. In the contents below, (25) will be firstly used in a coupling case with single-point connection in order to demonstrate its application.


### 2.2.1. Application to a Single-Point Coupling Case

For a single-point coupling between subsystemsa and b, the interface function can be simply defined as ψpσI=1 with p=1. A and B (see (11) and (12)) are simply the input-receptances of the two subsystems at their connecting points σIa and σIb, as(26)A=YaσIa=∑jϕja2σIaYjja,B=YbσIb=∑kϕkb2σIbYkkb.It thus gives (by (24))(27)Kab=1YaσIa+YbσIb≈1Yin,∞a+Yin,∞b.Then, (7)–(9) become(28)Gaa=YaI-KabΦaΦaTYa,(29)Gbb=YbI-KabΦbΦbTYb,(30)Gab=KabYaΦaΦbTYb.Equations (28)–(30) can then be combined with (13) to determine both the in-set (Daa and Dbb) and out-of-set (Dab) intermodal coupling stiffness terms.In the sense of weak coupling, that is, the second terms in brackets of the right sides of (28) and (29) are both negligible [14], the coupling stiffness between the jth mode of subsystem a and the kth mode of subsystem b, by (16), thus gives(31)Djkab=ϕjaσIaKabϕkbσIb.In [2], the intermodal coupling stiffness, for two subsystems connected by a single spring K, is given by(32)Djkab≈ϕjaxIaKϕkbxIb.It is seen that (31) and (32) are in a very similar form, except that Kab is a frequency-dependent function in the former equation while K is a constant in the latter one. It thus allows the results in [2] for spring-coupled subsystems to be extended for the coupling cases with rigid-point connections in a way fairly straightforward. These are described in the section below.


## 2.2.1. Application to a Single-Point Coupling Case

For a single-point coupling between subsystemsa and b, the interface function can be simply defined as ψpσI=1 with p=1. A and B (see (11) and (12)) are simply the input-receptances of the two subsystems at their connecting points σIa and σIb, as(26)A=YaσIa=∑jϕja2σIaYjja,B=YbσIb=∑kϕkb2σIbYkkb.It thus gives (by (24))(27)Kab=1YaσIa+YbσIb≈1Yin,∞a+Yin,∞b.Then, (7)–(9) become(28)Gaa=YaI-KabΦaΦaTYa,(29)Gbb=YbI-KabΦbΦbTYb,(30)Gab=KabYaΦaΦbTYb.Equations (28)–(30) can then be combined with (13) to determine both the in-set (Daa and Dbb) and out-of-set (Dab) intermodal coupling stiffness terms.In the sense of weak coupling, that is, the second terms in brackets of the right sides of (28) and (29) are both negligible [14], the coupling stiffness between the jth mode of subsystem a and the kth mode of subsystem b, by (16), thus gives(31)Djkab=ϕjaσIaKabϕkbσIb.In [2], the intermodal coupling stiffness, for two subsystems connected by a single spring K, is given by(32)Djkab≈ϕjaxIaKϕkbxIb.It is seen that (31) and (32) are in a very similar form, except that Kab is a frequency-dependent function in the former equation while K is a constant in the latter one. It thus allows the results in [2] for spring-coupled subsystems to be extended for the coupling cases with rigid-point connections in a way fairly straightforward. These are described in the section below.


## 3. Formulas of CLF Estimation and Coupling Strength Measurement

In [2], expressions of CLF estimations and coupling strength measurement were presented for two sets of oscillators coupled via a set of spring stiffness, as shown in Figure 2, based on the statistical power-energy relations between the two subsystems. These are first briefly reviewed in Section 3.1, and then extensions are made in Section 3.2 for the rigid-point coupling cases.Figure 2
Two sets of oscillators connected via a set of stiffness.


### 3.1. Formulas for Spring-Connected Subsystems

For simplicity, both subsystems are assumed to be uniform and homogenous so that both the modal masses (ma and mb) and modal loss factors (ηa and ηb) are constant.LetDjkab represent the coupling stiffness between the jth oscillator in set a and the kth oscillator in set b. The coupling strength between the two sets of spring-coupled oscillators can be measured by a so-called “coupling strength parameter” γ2 [1], which was given in [2] in the form of (33)γ2=EDjkab2mambηaηbω4,where ω is the centre frequency of the frequency band of interest, na and nb are the modal densities of subsystems a and b, and E· represents the expectation over the ensemble. The weak coupling region between the two subsystems corresponds to [1].(34)γ2≪1.In case of weak coupling (i.e. (34) is satisfied), the coupling loss factors ηab and ηba for the two subsystems a and b can then be estimated as [2](35)naηab=nbηba≈π2nanbmambω3EDjkab2.For a single-spring couplingK between a and b, for example, Djkab can be derived as (32). Substituting (32) into (35), ηab and ηba can then be calculated in a very simple manner [2]. In the subsection below, it will show how (33) and (35) can be extended readily to the coupling cases with multipoint interfaces.


### 3.2. Formulas for Rigid-Point Connections

For rigid-point couplings betweena and b, Djkab (given by (25)) becomes a frequency-dependent complex variable. Equations (33) and (35) thus need to be modified as (36)γ2=ED-jkab2mambηaηbω4,(37)naηab=nbηba≈π2nanbmambω3ED-jkab2,where ·- represents the expectation over the frequency band of interest.Considering the independence ofϕja, ϕkb, and Kppab, (23) leads to (38)ED-jkab2≈∑pEϕjaσpa2EK-ppab2Eϕkbσpb2.If ϕja and ϕkb are chosen to be mass-normalized, there exist(39)Eϕjaσpa2=1Ma,Eϕkbσpb2=1Mb,where Ma,b are the total masses of subsystems a and b. Substituting (39) into (38) and considering (24) yield(40)ED-jkab2≈NMaMb1Yin,∞a+Yin,∞b2.The above equation shows thatED-jkab2 is determined by the characteristic dynamic properties and the total mass of each subsystem as well as the number of coupling points between them regardless of the subsystem modal properties (modal mass, modal density, and modal damping). It thus indicates that the statistics of the dynamic coupling stiffness terms between two sets of subsystem modes are independent of the exact modal properties of each subsystem.Substituting (40) into (36) and considering ma=mb=1 (mass-normalized mode shapes), the coupling strength parameter can finally be expressed as(41)γ2≈NMaMbηaηbω41Yin,∞a+Yin,∞b2.The weak coupling region of two discretely connected subsystems can then be distinguished simply and explicitly byγ2≪1. As a result, the low bound of the frequency range of weak coupling region can be determined as (42)ωweak2≫1Yin,∞a+Yin,∞bNMaMbηaηb.For the frequency range where (42) is satisfied, the CLFs of two subsystems with discrete rigid-point couplings, by substituting (40) into (37), can be estimated as(43)naηab=nbηba≈Nπ2nanbMaMbω31Yin,∞a+Yin,∞b2.Equation (43) is expected to be valid if (42) (weak coupling) is satisfied and also the coherences among the coupling points are ignorable which commonly requires that coupling points are located at least one wavelength apart from each other [26]. A good supporting evidence is given here in the first instance to show the validity of (43).


#### 3.2.1. A Special Case of Single-Point Coupling

For a single-point connection between subsystemsa and b, (43) gives(44)ωnaηab≈π2nanbMaMbω21Yin,∞a+Yin,∞b2.Another well-known CLF expression is available in [1] (by substituting (10.1.3) into (10.1.6) in [1]), as(45)ωnaηab≈12π4ReZ∞aReZ∞bZ∞a+Z∞b2,where Z∞a represents the characteristic impedance of subsystem a, by [1](46)Z∞a=1iωYin,∞aReZ∞a=ωYin,∞a≈πna2Ma.Z∞b and ReZ∞b can be determined similarly by replacing the superscript a to b in (46). Substituting (46) into (45), it generates exactly the form of (44).Note that (45) was derived in [1] by using the wave approach and transmission coefficient; what (44) proposed here was derived by using the modal approach and intermodal couplings. The good consistency between (44) and (45) thus strongly supports the validity of the derivation procedure of (43).As stated previously in Section1 that, apart from (10.1.6), a CLF formula for the single-point connection proposed by a subsystem modal approach and the transmission coefficients is also available in [1] (i.e. (10.1.18), p. 188). However, because both (10.1.6) and (10.1.18) were developed by using one-dimensional subsystems and incoherent waves and then assumed (without being proven) to be applicable for two- and three-dimensional subsystems, one may become quite hesitated to extend directly (10.1.6) or (10.1.18) to coupling cases with many interface degrees of freedom in that more implicit and/or explicit assumptions have to be introduced and involved. Here the derivation of (43) has just proved that (10.1.6) or (10.1.18) can be extended to the multidegrees of freedom coupling cases but by using intermodal couplings rather than the wave transmission coefficients. Nevertheless, comparing (44) with (10.1.18) in [1], it can be seen clearly that the intermodal coupling-based solution has a much simpler form than the transmission-coefficient-based one, although the two expressions should be closely related to each other.


### 3.3. Result Discussion

From the analysis given in the above subsection, it is seen that (1) subsystem damping loss factors have strong influence on the coupling strength (41) but have little effect on the CLFs (43) which is in good agreement with the existing SEA literature; (2) both coupling strength and coupling loss factor are proportional to N but inversely proportional to Ma, Mb, and Yin,∞a+Yin,∞b2; and (3) the coupling strength seems affected more by frequency than by the CLFs.If it is assumed thatMa≈Mb≈M, considering Yin,∞a≈πna/2Maω [1], (41) and (42) can be further simplified as (47)Δa+Δb≫2πN,where, Δa and Δb are, respectively, the modal overlap factors of subsystems a and b. Equation (47) thus explains well that weak couplings tend to occur if both subsystems have high modal overlaps. Although being derived under the condition of Ma≈Mb≈M, (47) can be reasonably taken as valid for more general coupling cases in that it makes a good sense of physical intuitive.It is worth noting here that (47) implies that weak coupling may not necessarily require each individual subsystem to be with a high modal overlap factor. Instead, it may only need the summation of the subsystem modal overlap factors to be big enough. It thus makes it possible to extend the conventional SEA modelling to the mid-frequency range where one or two subsystems might with a quite low modal overlap factor [24, 27]. Further relevant investigations are underway [28].


## 3.1. Formulas for Spring-Connected Subsystems

For simplicity, both subsystems are assumed to be uniform and homogenous so that both the modal masses (ma and mb) and modal loss factors (ηa and ηb) are constant.LetDjkab represent the coupling stiffness between the jth oscillator in set a and the kth oscillator in set b. The coupling strength between the two sets of spring-coupled oscillators can be measured by a so-called “coupling strength parameter” γ2 [1], which was given in [2] in the form of (33)γ2=EDjkab2mambηaηbω4,where ω is the centre frequency of the frequency band of interest, na and nb are the modal densities of subsystems a and b, and E· represents the expectation over the ensemble. The weak coupling region between the two subsystems corresponds to [1].(34)γ2≪1.In case of weak coupling (i.e. (34) is satisfied), the coupling loss factors ηab and ηba for the two subsystems a and b can then be estimated as [2](35)naηab=nbηba≈π2nanbmambω3EDjkab2.For a single-spring couplingK between a and b, for example, Djkab can be derived as (32). Substituting (32) into (35), ηab and ηba can then be calculated in a very simple manner [2]. In the subsection below, it will show how (33) and (35) can be extended readily to the coupling cases with multipoint interfaces.


## 3.2. Formulas for Rigid-Point Connections

For rigid-point couplings betweena and b, Djkab (given by (25)) becomes a frequency-dependent complex variable. Equations (33) and (35) thus need to be modified as (36)γ2=ED-jkab2mambηaηbω4,(37)naηab=nbηba≈π2nanbmambω3ED-jkab2,where ·- represents the expectation over the frequency band of interest.Considering the independence ofϕja, ϕkb, and Kppab, (23) leads to (38)ED-jkab2≈∑pEϕjaσpa2EK-ppab2Eϕkbσpb2.If ϕja and ϕkb are chosen to be mass-normalized, there exist(39)Eϕjaσpa2=1Ma,Eϕkbσpb2=1Mb,where Ma,b are the total masses of subsystems a and b. Substituting (39) into (38) and considering (24) yield(40)ED-jkab2≈NMaMb1Yin,∞a+Yin,∞b2.The above equation shows thatED-jkab2 is determined by the characteristic dynamic properties and the total mass of each subsystem as well as the number of coupling points between them regardless of the subsystem modal properties (modal mass, modal density, and modal damping). It thus indicates that the statistics of the dynamic coupling stiffness terms between two sets of subsystem modes are independent of the exact modal properties of each subsystem.Substituting (40) into (36) and considering ma=mb=1 (mass-normalized mode shapes), the coupling strength parameter can finally be expressed as(41)γ2≈NMaMbηaηbω41Yin,∞a+Yin,∞b2.The weak coupling region of two discretely connected subsystems can then be distinguished simply and explicitly byγ2≪1. As a result, the low bound of the frequency range of weak coupling region can be determined as (42)ωweak2≫1Yin,∞a+Yin,∞bNMaMbηaηb.For the frequency range where (42) is satisfied, the CLFs of two subsystems with discrete rigid-point couplings, by substituting (40) into (37), can be estimated as(43)naηab=nbηba≈Nπ2nanbMaMbω31Yin,∞a+Yin,∞b2.Equation (43) is expected to be valid if (42) (weak coupling) is satisfied and also the coherences among the coupling points are ignorable which commonly requires that coupling points are located at least one wavelength apart from each other [26]. A good supporting evidence is given here in the first instance to show the validity of (43).


### 3.2.1. A Special Case of Single-Point Coupling

For a single-point connection between subsystemsa and b, (43) gives(44)ωnaηab≈π2nanbMaMbω21Yin,∞a+Yin,∞b2.Another well-known CLF expression is available in [1] (by substituting (10.1.3) into (10.1.6) in [1]), as(45)ωnaηab≈12π4ReZ∞aReZ∞bZ∞a+Z∞b2,where Z∞a represents the characteristic impedance of subsystem a, by [1](46)Z∞a=1iωYin,∞aReZ∞a=ωYin,∞a≈πna2Ma.Z∞b and ReZ∞b can be determined similarly by replacing the superscript a to b in (46). Substituting (46) into (45), it generates exactly the form of (44).Note that (45) was derived in [1] by using the wave approach and transmission coefficient; what (44) proposed here was derived by using the modal approach and intermodal couplings. The good consistency between (44) and (45) thus strongly supports the validity of the derivation procedure of (43).As stated previously in Section1 that, apart from (10.1.6), a CLF formula for the single-point connection proposed by a subsystem modal approach and the transmission coefficients is also available in [1] (i.e. (10.1.18), p. 188). However, because both (10.1.6) and (10.1.18) were developed by using one-dimensional subsystems and incoherent waves and then assumed (without being proven) to be applicable for two- and three-dimensional subsystems, one may become quite hesitated to extend directly (10.1.6) or (10.1.18) to coupling cases with many interface degrees of freedom in that more implicit and/or explicit assumptions have to be introduced and involved. Here the derivation of (43) has just proved that (10.1.6) or (10.1.18) can be extended to the multidegrees of freedom coupling cases but by using intermodal couplings rather than the wave transmission coefficients. Nevertheless, comparing (44) with (10.1.18) in [1], it can be seen clearly that the intermodal coupling-based solution has a much simpler form than the transmission-coefficient-based one, although the two expressions should be closely related to each other.


## 3.2.1. A Special Case of Single-Point Coupling

For a single-point connection between subsystemsa and b, (43) gives(44)ωnaηab≈π2nanbMaMbω21Yin,∞a+Yin,∞b2.Another well-known CLF expression is available in [1] (by substituting (10.1.3) into (10.1.6) in [1]), as(45)ωnaηab≈12π4ReZ∞aReZ∞bZ∞a+Z∞b2,where Z∞a represents the characteristic impedance of subsystem a, by [1](46)Z∞a=1iωYin,∞aReZ∞a=ωYin,∞a≈πna2Ma.Z∞b and ReZ∞b can be determined similarly by replacing the superscript a to b in (46). Substituting (46) into (45), it generates exactly the form of (44).Note that (45) was derived in [1] by using the wave approach and transmission coefficient; what (44) proposed here was derived by using the modal approach and intermodal couplings. The good consistency between (44) and (45) thus strongly supports the validity of the derivation procedure of (43).As stated previously in Section1 that, apart from (10.1.6), a CLF formula for the single-point connection proposed by a subsystem modal approach and the transmission coefficients is also available in [1] (i.e. (10.1.18), p. 188). However, because both (10.1.6) and (10.1.18) were developed by using one-dimensional subsystems and incoherent waves and then assumed (without being proven) to be applicable for two- and three-dimensional subsystems, one may become quite hesitated to extend directly (10.1.6) or (10.1.18) to coupling cases with many interface degrees of freedom in that more implicit and/or explicit assumptions have to be introduced and involved. Here the derivation of (43) has just proved that (10.1.6) or (10.1.18) can be extended to the multidegrees of freedom coupling cases but by using intermodal couplings rather than the wave transmission coefficients. Nevertheless, comparing (44) with (10.1.18) in [1], it can be seen clearly that the intermodal coupling-based solution has a much simpler form than the transmission-coefficient-based one, although the two expressions should be closely related to each other.


## 3.3. Result Discussion

From the analysis given in the above subsection, it is seen that (1) subsystem damping loss factors have strong influence on the coupling strength (41) but have little effect on the CLFs (43) which is in good agreement with the existing SEA literature; (2) both coupling strength and coupling loss factor are proportional to N but inversely proportional to Ma, Mb, and Yin,∞a+Yin,∞b2; and (3) the coupling strength seems affected more by frequency than by the CLFs.If it is assumed thatMa≈Mb≈M, considering Yin,∞a≈πna/2Maω [1], (41) and (42) can be further simplified as (47)Δa+Δb≫2πN,where, Δa and Δb are, respectively, the modal overlap factors of subsystems a and b. Equation (47) thus explains well that weak couplings tend to occur if both subsystems have high modal overlaps. Although being derived under the condition of Ma≈Mb≈M, (47) can be reasonably taken as valid for more general coupling cases in that it makes a good sense of physical intuitive.It is worth noting here that (47) implies that weak coupling may not necessarily require each individual subsystem to be with a high modal overlap factor. Instead, it may only need the summation of the subsystem modal overlap factors to be big enough. It thus makes it possible to extend the conventional SEA modelling to the mid-frequency range where one or two subsystems might with a quite low modal overlap factor [24, 27]. Further relevant investigations are underway [28].


## 4. Numerical Examples

In this section, numerical examples are given to illustrate the performance of the two new formulas (see (41) and (43)) for discretely connected subsystems.The numerical model employed is shown in Figure3, which is comprised of two plates connected by rigid points. The number of coupling points considered is varied between 1 and 4, and both “large”- and “small”-spacing cases are involved. (Large spacing means the points are spaced at least one wavelength apart for the frequency range of interest, while small spacing is less than one wavelength apart). For each coupling case, coupling loss factors between the two plates are first calculated by (43) and then compared with the results generated from the PIM method. Here, the PIM-based results are obtained from the ensemble average of the CLFs of 30 samples of the plate-plate models. For each sample, both plates have fixed modal densities and total masses, while their natural frequencies and mode shapes may differ from one to another. This is achieved by varying the ratio of the length and width of each plate model while keeping its area unchanged (±5% of the nominal value). Frequency-averaged energy response and power transmission for each sample of the ensemble under randomly located force excitations are first calculated by exact modal approach, and then ensemble averages are taken over the whole 30 samples. Coupling loss factors can finally be estimated, by its definition given in [1], from the frequency- and ensemble-averaged power-energy relations between the two plates. The performance of the two new formulas (see (41) and (43)) can then be evaluated by comparing with the PIM-based results. Only the bending motions of the numerical model are considered in the following calculations.Figure 3
Two plates connected via discrete points.


### 4.1. System Model Description

In Figure3, both plates 1 and 2 are rectangular with simply supported edges. Plate 1 has an area of 0.63 m2 while plate 2 has an area of 0.8 m2, and both are with the same thickness of 2 mm. The material of the plates is chosen to be steel with Young’s modulus of 2.1×1011 N/m2, density of 7.85×103 kg/m3, Poisson ratio of 0.3, and damping loss factor of 0.01. The corresponding modal densities of plates 1 and 2 are 0.016 modes/rad/s and 0.02 modes/rad/s, respectively, and the characteristic impedances Z∞1=Z∞2=393.2 N/m/s.The frequency range considered is up to 2000 Hz which comprises approximately 200 resonant modes of plate 1 and 250 resonant modes of plate 2. The plate wavelength at 500 Hz is about 0.2 m while the wavelength at 2000 Hz is about 0.1 m. The coupling points on plates 1 and 2 are all randomly chosen over surface of each plate but deliberately kept to be spaced at least 0.2 m apart from all the plate edges in order to minimize the boundary effects for the frequencies above 500 Hz.


## 4.1. System Model Description

In Figure3, both plates 1 and 2 are rectangular with simply supported edges. Plate 1 has an area of 0.63 m2 while plate 2 has an area of 0.8 m2, and both are with the same thickness of 2 mm. The material of the plates is chosen to be steel with Young’s modulus of 2.1×1011 N/m2, density of 7.85×103 kg/m3, Poisson ratio of 0.3, and damping loss factor of 0.01. The corresponding modal densities of plates 1 and 2 are 0.016 modes/rad/s and 0.02 modes/rad/s, respectively, and the characteristic impedances Z∞1=Z∞2=393.2 N/m/s.The frequency range considered is up to 2000 Hz which comprises approximately 200 resonant modes of plate 1 and 250 resonant modes of plate 2. The plate wavelength at 500 Hz is about 0.2 m while the wavelength at 2000 Hz is about 0.1 m. The coupling points on plates 1 and 2 are all randomly chosen over surface of each plate but deliberately kept to be spaced at least 0.2 m apart from all the plate edges in order to minimize the boundary effects for the frequencies above 500 Hz.


## 5. Results and Discussions

Figure4 shows the coupling strength parameter γ2 calculated by (41) when the plates are connected by 1 point, 2 points, and 4 points, respectively. It is seen that the frequencies corresponding to γ2=1 for the three coupling cases are about 280 Hz, 400 Hz, and 560 Hz, respectively. The weak coupling regions can then be distinguished explicitly by (34), or equivalently, by (42).Figure 4
Coupling strength parameterγ2 between two coupled plates: N=1 (dotted line); N=2 (real line); N=4 (dash-dotted line).Figures5–7 compare the coupling loss factor results calculated by (43) and by the PIM-based simulations when N=1, 2, and 4, respectively. Here it should be mentioned that, in Figures 6 and 7, the PIM-based CLF solutions correspond to the “widely-spaced” coupling case; that is, the coupling points are located more than 0.2 m apart from each other so that the coherences between them can be neglected for the frequencies above 500 Hz.Figure 5
Comparison ofη12 when N=1: PIM simulations (real thin lines); ensemble average of the PIM solutions (thick real line); (43) (dotted line).Figure 6
Comparison ofη12 when N=2 (widely-spaced): PIM simulations (real thin lines); ensemble average of the PIM solutions (thick real line); (43) (dotted line).Figure 7
Comparison ofη12 when N=4 (widely-spaced): PIM simulations (real thin lines); ensemble average of the PIM solutions (thick real line); (43) (dotted line).It is seen clearly from Figures5–7 that the results of (43) tend to agree well with the PIM solutions from the frequency ranges of ω≥550 Hz for N=1 (Figure 5), ω≥650 Hz for N=2 (Figure 6), and ω≥900 Hz for N=4 (Figure 7). It thus indicates that the weak coupling regions start from about 550 Hz for the case of N=1, about 650 Hz for the case of N=2, and about 900 Hz for the case of N=4. It is very interesting to see in Figure 4 that all these weak coupling frequencies correspond to a coupling strength parameter of γ2≈0.3. It thus implies that (43) is valid for the frequencies where γ2≤0.3 is satisfied, at least for the coupling cases considered in the present investigation.The apparent fluctuations of the PIM solutions occurring at lower frequencies (below 500 Hz) can be explained as the result of insufficient statistical overlap of the ensemble examples [27]. They gradually improve as frequency increases. The poor agreements at lower frequencies are due to the strong coupling effects between the two plates.Figures8 and 9 show the comparisons of η12 calculated by (43) and that by the PIM simulations for the coupling cases of N=2 and N=4, respectively, of which the coupling points are spaced deliberately at less than 0.1 m apart from each other, that is, the “narrowly-spaced” couplings. It is seen that the agreements between the two sets of results in Figures 8 and 9 are relatively poor compared to those observed in Figures 6 and 7 for the “widely-spaced” couplings. Clearly this is because (43) was derived by ignoring the coherences among different coupling points, which does not meet when the coupling points are “narrowly-spaced” in Figures 8 and 9. However, it is also seen that (43) works well for certain frequency regions in Figures 8 and 9, such as for the frequency ranges of 600–1000 Hz and 1200–1700 Hz in Figure 8 and 900–1700 Hz in Figure 9, even when the coupling points are closely located with each other. It thus suggests that the coherences between the coupling points may be not so important for certain frequency ranges, at least in the sense of ensemble average.Figure 8
Comparison ofη12 when N=2 (narrowly-spaced): PIM simulations (real thin lines); ensemble average of the PIM solutions (thick real line); (43) (dotted line).Figure 9
Comparison ofη12 when N=4 (narrowly-spaced): PIM simulations (real thin lines); ensemble average of the PIM solutions (thick real line); (43) (dotted line).


## 6. Conclusions

This paper mainly concerns the estimation of coupling loss factors for SEA subsystems connected via discrete interfaces. First of all, the dynamic interactions between two subsystems connected via discrete interfaces are expressed explicitly as an intermodal coupling stiffness matrix by a subsystem modal approach given in [24]. If each subsystem is highly modal dense and each interface degree of freedom is assumed to act independently, the dynamic interactions are found to be strongly dominated by the modal couplings between the two sets of subsystem modes while the couplings within the same set of modes can be neglected. In which case, the modal stiffness coupling terms simply reduce to a function of the characteristic receptances of the two subsystems as well as their mode shapes at the interfaces. If further ensemble- and frequency-averaged, these modal coupling terms can be reduced to a function which depends only on the characteristic dynamic properties and the total mass of each subsystem as well as the number of interface degrees of freedom. The results can then be accommodated within the theoretical lines of [1, 2] to yield two simple formulas for both coupling strength measurement and coupling loss factor estimation, both depending only on the gross modal properties and the number of coupling points but regardless of any other detailed subsystem modal properties. The consistency and difference of the present CLF formula with and from the traditional CLF solutions in [1] are discussed. Meanwhile, the present theory shows clearly that (1) both coupling strength and CLFs are proportional to the coupling points N and inversely proportional to the total subsystem masses (Ma and Mb) and Yin,∞a+Yin,∞b2; (2) the coupling strength is strongly determined by the subsystem damping loss factors while the CLFs are independent of subsystem damping, which are in good agreement with the existing SEA literature; and (3) the coupling strength tends to be more sensitive to the frequencies than the coupling loss factors.To illustrate the validity and accuracy of the present technique, a few of numerical examples are presented for the plate-plate coupling models with discrete point connections. A very good performance of the present CLF formula is observed when compared with the commonly used power injection method (PIM).



---

*Source: 3636401-2016-11-22.xml*