Abstract

The aim of this paper is to develop a new frequency response function- (FRF-) based indirect inverse substructuring method without measuring system-level FRFs in the coupling DOFs for the analysis of the dynamic characteristics of a three-substructure coupled product transport system with rigid and flexible coupling. By enforcing the dynamic equilibrium conditions at the coupling coordinates and the displacement compatibility conditions, a closed-form analytical solution to inverse substructuring analysis of multisubstructure coupled product transport system is derived based on the relationship of easy-to-monitor component-level FRFs and the system-level FRFs at the coupling coordinates. The proposed method is validated by a lumped mass-spring-damper model, and the predicted coupling dynamic stiffness is compared with the direct computation, showing exact agreement. The method developed offers an approach to predict the unknown coupling dynamic stiffness from measured FRFs purely. The suggested method may help to obtain the main controlling factors and contributions from the various structure-borne paths for product transport system.

1. Introduction

Packaging dynamics for protection of the packaged product has attracted increasingly more attention during the last two decades [19]. In distribution process, the product, packaging, and vehicle constitute a complex built-up product transport system. The packaged product will be damaged when the dynamic response due to the environmental vibration and shock exceeds some limit value [10]. To protect the packaged product, one needs to isolate the packaged item with flexible cushions. Hence, it is of significant importance to identify the coupling dynamic stiffness between packaged product and vehicle system for cushioning packaging design [11, 12]. However, the interaction between the packaged product and vehicle is extremely complicated, including internal and external cushion, container, pallet, nail, and rope. It is difficult to predict the coupling dynamic stiffness accurately by reported identification methods [1214], as some of the frequency response functions (FRFs) required cannot be easily measured between the complex coupling interfaces for a product transport system. To address this issue, Zhen et al. [15] proposed an inverse formulation that the FRFs of individual components and dynamic characteristics of the coupling elements can be predicted directly from the system-level FRFs. Although this technique includes several mathematical operations that might be very sensitive to measurement errors and/or and inconsistencies, it has its unique advantages in various application fields since no component-level spectra response is needed. Then it was applied to study the dynamics of product transport system [16, 17], showing its great application prospect in industry. In spite of those promising successful applications, the two-component structure assumption of the inverse substructure theory proposed by Zhen restricts its application for a complex product transport system. The discretely connected system should be taken as a multisubstructure coupled system in order to perform the parametric analysis of each substructure and connector more efficiently. In our recent published work [18, 19], the inverse substructuring method of three-component coupled product transport system for both single coordinate coupling case and multicoordinate coupling case was proposed and applied to perform the parameter analysis of product transport system. However, in the previous work, the dynamic stiffness of the coupling interface was predicted purely from the system-level FRFs. Actually, the system-level FRFs from coupling degree of freedoms (DOFs) may not be measured accurately due to the difficulties of vibration excitation and measurement for the coupled interface between packaged product and vehicle within the limited space. In this paper, an indirect inverse substructuring method is proposed to identify the coupling dynamic stiffness. Instead of identification of the coupling dynamic stiffness only from system-level FRFs, the proposed method aims to obtain the packaging stiffness accurately from easy-to-monitor DOFs, which will avoid the difficulties of vibration excitation and measurement due to space limitation. The accuracy of the proposed theory is first validated by a lumped parameter model.

2. Indirect Inverse Substructuring Method for Multicomponent Coupled Product Transport System with Rigid and Flexible Coupling

A free substructure is shown in Figure 1, where , , and are the response, excitation, and coupling coordinates, respectively. At the coupling coordinate, the substructure is connected to the rest of the system. The governing equations of motion can be written as [16]where , , and are the mass, damping, and stiffness matrices, respectively, and the subscript refers to substructure . Assuming forced periodic excitation, the steady-state vibratory response in frequency domain can be expressed as where is the excitation frequency. This form is used to generate the necessary frequency response function (FRF) matrix of the substructure . Here, represents the frequency response function matrices of substructure (response at coordinate due to an excitation at coordinate); for example, for and , is the transfer response function between the excitation coordinate and response coordinate of substructure .

Consider a system with three substructures , , and coupled at interface by a rigid connector between and and at interface by a flexible connector defined as between and (see Figure 2).

By applying the substructuring method, the system-level FRFs can be predicted from the component-level FRFs. First, take the structure coupled by substructure and substructure as a system, and the FRFs of structure can be expressed by the FRFs of substructure and substructure as [18, 19]whereand .

Then, take the system as a coupled one by substructure and substructure (here, structure is treated as a component), and the FRFs of the system can be expressed by the FRFs of substructure and substructure aswhereand .

From (5) and (3), it can be seen that is expressed by and , and is expressed by and . Substituting (5) into (3) yields the FRFs of the system expressed by the FRFs of the three substructures, ,, and . Because the response at coordinate of substructure due to an excitation at coordinate of substructure is most concerned, we give here the explicit expression of the system-level FRF by the FRFs of the three substructures [18, 19]:

Similarly, one can obtain the system-level FRF for multicoordinate coupling case as (see Figure 3)

The above derivation develops a method for a three-substructure coupled system to obtain the system-level responses from the prior knowledge of FRFs of the three substructures and the dynamic stiffness of the coupling interfaces. However, for many complex structure systems, the substructure-level FRFs may not be easily obtained. Besides, the dynamic stiffness of the coupling interface may be unknown. If the physical system exists but is not conveniently separable into two or more substructures, it is desirable to express the problem in terms of measurable system-level FRFs.

2.1. Closed-Form Solution for Single Coupling

The coupled system is first decoupled as two substructures, and , as illustrated in Figure 2. The force compatibility of the coupling coordinates between substructure and substructure leads to the following equation:The response vector of the coupling coordinates can then be solved as where

The coefficient matrix in (10) is the transfer matrix between the coupling coordinates on substructure and substructure , and it can be written in terms of system-level FRFs as

Comparing (10) and (12) and recalling (5) reveal a series of equations relating the component-level FRFs to system-level FRFs, from which we can obtain the expression of the coupling dynamic stiffness as depicted in

Equation (13) provides an inverse method to predict the coupling dynamic stiffness for product transport system. But the problem remaining is that some system-level FRFs cannot be easily measured from the coupling DOFs such as . To overcome this shortcoming, we recall back (7), make an inverse formulation, and then get

From (14), the coupling dynamic stiffness of the interface between packaged product and vehicle can be predicted from both system-level FRFs and component-level FRFs; all of the FRFs from DOFs of coupling c2 are measured in component level, which avoids the difficulties of vibration excitation and pickup. As the coupling dynamic stiffness is predicted from both system-level FRFs and component-level FRFs, we call the method proposed indirect inverse substructuring method.

2.2. Closed-Form Solution for Multicoordinate Coupling

As depicted in Figure 3, the components are connected by a set of [] of dimension . In many cases, [] is nearly diagonal especially in product transport system. Similar to the single coupling problem, the dynamic equilibrium conditions are applied to the coupling coordinates on both sides of the connecting springs to obtainwhich is inverted to express the response vector of the coupling coordinates aswhere .

The coefficient matrix in (16) is the transfer matrix [] between the coupling coordinates on substructure and substructure , and it can be written in terms of system-level FRFs as Comparing (16) and (17) and recalling (8) reveal a series of equations relating the component-level FRFs to system-level FRFs, from which we can obtain the expression of the coupling dynamic stiffness between and as in (18) after a lengthy derivation:

Equation (18) provides an inverse method to predict the coupling dynamic stiffness for product transport system. But the remaining problem is that some system-level FRFs cannot be easily measured from the coupling DOFs such as . To overcome this shortcoming, we recall back (8) and make an inverse formulation:where .

Equation (19) can be used to predict the coupling dynamic stiffness from both system-level FRFs and component-level FRFs, avoiding the difficulties of vibration excitation and pickup.

3. Numerical Validation Using a Lumped Parameter Model

To verify the indirect inverse substructure proposed above, a lumped mass-spring-damper model with three substructures, shown in Figure 4, is taken as an example. The specific parameters of the system are listed in Table 1. The necessary FRFs of the free substructures and the coupled system are generated from (2), and the computed system response functions as well as the necessary component response functions are used to predict the coupling dynamic stiffness applying (19) for packaging interface. Then, the results are compared to direct calculations, as shown in Figure 5. Specifically, for , , , and , , and , , , , and , and , , and , , and . All necessary FRFs are calculated from (2). The matrices for components , , , and are expressed explicitly as follows:

component :

component :

component :

system :

From Figure 5, we can see that the predicted dynamic stiffness is in exact agreement with the direct computation. Hence, the proposed indirect inverse substructuring method demonstrates its validity for identifying the dynamic stiffness at coupling interfaces. It should be noted that the formulation may be sensitive to input random errors. As illustrated in Figure 6, the predicted system-level frequency response function is affected remarkably by the input random errors of measured system-level frequency response functions, implying that the formulation for prediction of coupling stiffness is extremely sensitive to the experiment errors. A further study on the error propagation as well as the elimination approach will be conducted in the future.

4. Conclusions

Three sets of equations for multicomponent coupled product transport system are obtained, including the dynamic equilibrium conditions, the displacement compatibility conditions in addition to the relationship of the system-level FRFs, and the selected easy-to-monitor component-level FRFs. Then, the closed-form analytical solution to inverse substructuring analysis of multisubstructure coupled product transport system with rigid and flexible coupling is derived. The method developed offers an approach to predict the unknown coupling dynamic stiffness from system-level FRFs and easy-to-monitor component-level FRFs, avoiding the difficulties of vibration excitation and/or measurement in the coupling interface. Then, the proposed method is validated by a lumped mass-spring-damper model, and the coupling dynamic stiffness is compared with the direct computation, showing exact agreement. The method developed offers an approach to predict the unknown coupling dynamic stiffness from system-level FRFs and component-level FRFs which can be easily measured. The suggested method may help to obtain the main controlling factors and contribution from the various structure-borne paths for product transport system, which may certainly facilitate the cushioning packaging design.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to appreciate the financial support of National Natural Science Foundation of China (Grant no. 51205167), Research Fund for the Doctoral Program of Higher Education of China (Grant no. 20120093120014), and Fundamental Research Funds for the Central Universities (Grant no. JUSRP51403A), the project was supported by the Foundation (Grant no. 2013010) of Tianjin Key Laboratory of Pulp & Paper (Tianjin University of Science & Technology) and by the Foundation of Jiangsu Key Laboratory of Advanced Food Manufacturing Equipment and Technology (Grant no. FM-201403), and the research was also sponsored by Qing Lan Project of Jiangsu Province (2014).