一種針對(duì)含非線性子結(jié)構(gòu)多點(diǎn)耦合運(yùn)輸包裝系統(tǒng)的解耦方法

王維凱，王軍，盧立新，潘嘹，侯雪

一種針對(duì)含非線性子結(jié)構(gòu)多點(diǎn)耦合運(yùn)輸包裝系統(tǒng)的解耦方法

王維凱1，王軍1，盧立新1，潘嘹1，侯雪2

（1.江南大學(xué)，無錫 214122；2.汕頭東風(fēng)印刷股份有限公司無錫分公司，無錫 214000）

考慮到運(yùn)輸包裝系統(tǒng)耦合形式復(fù)雜，包裝材料及包裝結(jié)構(gòu)具有非線性特性，不容易測量局部物理參數(shù)，需要對(duì)傳統(tǒng)逆向子結(jié)構(gòu)方法進(jìn)行優(yōu)化，使之能夠求解非線性多點(diǎn)耦合系統(tǒng)中子結(jié)構(gòu)的動(dòng)態(tài)響應(yīng)特性。使用描述函數(shù)法將非線性的運(yùn)輸包裝系統(tǒng)線性化，測量其在若干特定振動(dòng)幅值下的頻率響應(yīng)函數(shù)；之后，應(yīng)用逆向子結(jié)構(gòu)方法和參數(shù)識(shí)別方法，計(jì)算包裝件的模態(tài)參數(shù)；最后，擬合包裝件模態(tài)參數(shù)與振動(dòng)幅值之間的關(guān)系，構(gòu)建函數(shù)來描述包裝件的動(dòng)態(tài)響應(yīng)特性。在集總參數(shù)模型中，解耦預(yù)測值與實(shí)際值吻合；在有限元模型中，對(duì)響應(yīng)峰值的預(yù)測誤差小于5%，對(duì)響應(yīng)跳躍現(xiàn)象所在頻率的預(yù)測誤差小于3%。該研究將傳統(tǒng)逆向子結(jié)構(gòu)方法的應(yīng)用范圍拓展到了非線性多點(diǎn)耦合系統(tǒng)，對(duì)復(fù)雜運(yùn)輸包裝系統(tǒng)動(dòng)力學(xué)模型的構(gòu)建和防振包裝的設(shè)計(jì)具有指導(dǎo)意義。

逆向子結(jié)構(gòu)方法；運(yùn)輸包裝；振動(dòng)；非線性

在現(xiàn)代緩沖包裝設(shè)計(jì)流程中，設(shè)計(jì)師經(jīng)常借助動(dòng)力學(xué)計(jì)算和有限元仿真來確定緩沖材料和包裝結(jié)構(gòu)的參數(shù)。其中，包裝件動(dòng)態(tài)響應(yīng)特性的獲取是精準(zhǔn)有效建立其動(dòng)力學(xué)模型的關(guān)鍵[1]，然而，在實(shí)際運(yùn)輸過程中，產(chǎn)品、包裝與運(yùn)載車輛耦合聯(lián)結(jié)在一起，構(gòu)成復(fù)雜的運(yùn)輸包裝系統(tǒng)，致使包裝件的動(dòng)態(tài)響應(yīng)特性難以被直接測量[2]。此外，包裝件在耦合狀態(tài)與非耦合狀態(tài)下的響應(yīng)特性不同，對(duì)其單獨(dú)測量的結(jié)果不能代表其在耦合狀態(tài)下的特性。針對(duì)這一問題，王志偉、王軍等[3-6]將逆向子結(jié)構(gòu)方法應(yīng)用到了運(yùn)輸包裝領(lǐng)域，實(shí)現(xiàn)了對(duì)運(yùn)輸包裝系統(tǒng)中子結(jié)構(gòu)動(dòng)態(tài)響應(yīng)特性的間接預(yù)測。

為建立與實(shí)際情況更加接近的理論模型，逆向子結(jié)構(gòu)方法在近年來持續(xù)發(fā)展。主要趨勢有：為適應(yīng)集裝、堆碼等復(fù)雜耦合情況，從單點(diǎn)耦合系統(tǒng)向多點(diǎn)耦合系統(tǒng)發(fā)展[7-8]；為重點(diǎn)考慮產(chǎn)品中的易損部件，從二級(jí)耦合系統(tǒng)向多級(jí)耦合系統(tǒng)發(fā)展[9]；設(shè)法降低誤差，如利用奇異值分解方法抑制計(jì)算過程中矩陣求逆導(dǎo)致的誤差放大[10-11]。

在以往研究中，均假設(shè)運(yùn)輸包裝系統(tǒng)是線性的，而實(shí)際包裝材料與包裝結(jié)構(gòu)具有非線性特性[12-14]，在使用傳統(tǒng)逆向子結(jié)構(gòu)方法進(jìn)行分析時(shí)，將非線性結(jié)構(gòu)視作線性結(jié)構(gòu)會(huì)導(dǎo)致求解結(jié)果不準(zhǔn)確。對(duì)此，筆者在之前的研究中，結(jié)合逆向子結(jié)構(gòu)方法與參數(shù)識(shí)別技術(shù)，實(shí)現(xiàn)了對(duì)含非線性子結(jié)構(gòu)耦合系統(tǒng)的解耦，其中只考慮了單點(diǎn)耦合的情況[15]。文中結(jié)合運(yùn)輸包裝實(shí)際，針對(duì)含非線性子結(jié)構(gòu)的多點(diǎn)耦合系統(tǒng)做進(jìn)一步討論，拓展了逆向子結(jié)構(gòu)方法的應(yīng)用范圍，并分別在集總參數(shù)模型和有限元模型中進(jìn)行了驗(yàn)證。

1 理論方法

包裝件與運(yùn)載車輛之間的耦合關(guān)系見圖1a，該耦合系統(tǒng)可抽象化為圖1b所示的物理模型，文中以該模型為基礎(chǔ)進(jìn)行理論分析。

1.1 非線性系統(tǒng)線性化

逆向子結(jié)構(gòu)方法需要測量系統(tǒng)的頻率響應(yīng)函數(shù)來進(jìn)行解耦運(yùn)算，而非線性系統(tǒng)不具備傳統(tǒng)意義上的頻率響應(yīng)函數(shù)，因此需要將其線性化。

非線性結(jié)構(gòu)的動(dòng)力學(xué)方程可以寫成：

式中：、和為質(zhì)量矩陣、阻尼矩陣和剛度矩陣；和L分別為外部激勵(lì)力向量和內(nèi)部非線性力向量；為位移向量。

假設(shè)激勵(lì)是簡諧振動(dòng)，根據(jù)歐拉公式，在時(shí)域上能被寫成：

式中：為簡諧激勵(lì)的幅值；為激勵(lì)的角頻率；e為自然常數(shù)，i為虛數(shù)單位。

非線性結(jié)構(gòu)的穩(wěn)態(tài)響應(yīng)以級(jí)數(shù)形式表示為：

式中：為結(jié)構(gòu)的第階諧波響應(yīng)；為第階諧波響應(yīng)的幅值。通常，第1階諧波是響應(yīng)的主要成分，而其他階次諧波相對(duì)較小[16]。在忽略高階諧波的前提下，式（3）可寫成：

圖1 包裝件–運(yùn)載車輛耦合系統(tǒng)

式中：為第1階諧波響應(yīng)的幅值。與之對(duì)應(yīng)地，非線性內(nèi)力也只考慮第1階諧波，并設(shè)為第1階內(nèi)部非線性力的幅值，得到：

將式（2）、（4）和（5）代入到式（1）中，得到：

如果該結(jié)構(gòu)的非線性特征體現(xiàn)為：剛度隨響應(yīng)幅值的變化而變化。那么，根據(jù)描述函數(shù)法[16]，可以把非線性內(nèi)力的幅值寫成響應(yīng)幅值的一次函數(shù)，即：

式中：為描述函數(shù)矩陣。矩陣內(nèi)部各位置的元素為：

式中：Δ為位于矩陣中第行列的元素；Δ為位于矩陣中第行列的元素；v為結(jié)構(gòu)中坐標(biāo)和地面之間的非線性力的描述函數(shù)；v為結(jié)構(gòu)中坐標(biāo)和坐標(biāo)之間非線性力的描述函數(shù)。v和v的值均與響應(yīng)幅值有關(guān)，因此，是隨變化的函數(shù)。在此基礎(chǔ)上，式（7）能被表示成式（9）的形式。

其中：

在文中被稱作準(zhǔn)線性頻率響應(yīng)函數(shù)矩陣，它的值與結(jié)構(gòu)的振動(dòng)幅值有關(guān)，只能表征結(jié)構(gòu)在特定振動(dòng)幅值下的響應(yīng)特性。

至此，實(shí)現(xiàn)了非線性結(jié)構(gòu)的線性化，得到了非線性結(jié)構(gòu)的準(zhǔn)線性頻率響應(yīng)函數(shù)。

1.2 逆向子結(jié)構(gòu)方法

將圖1所示的包裝件—運(yùn)載車輛耦合系統(tǒng)抽象化，考慮一個(gè)由非線性子結(jié)構(gòu)A和線性子結(jié)構(gòu)B組成的系統(tǒng)S，見圖2。2個(gè)子結(jié)構(gòu)之間通過對(duì)耦合點(diǎn)連接在一起。圖1中，a和b分別表示非線性子結(jié)構(gòu)A和線性子結(jié)構(gòu)B的內(nèi)部自由度，a和b分別表示它們的耦合自由度，C為耦合剛度矩陣。

圖2 包含非線性子結(jié)構(gòu)的多點(diǎn)耦合系統(tǒng)

其準(zhǔn)線性頻率響應(yīng)函數(shù)矩陣可表示為：

式中：S、S和S分別為系統(tǒng)S的質(zhì)量矩陣、阻尼矩陣和剛度矩陣；S為系統(tǒng)S中以描述函數(shù)法線性化的非線性內(nèi)力矩陣。S和S隨響應(yīng)幅值變化。

系統(tǒng)S的動(dòng)力學(xué)方程可寫作：

其中：

舉例說明式（13）中變量的物理含義：[S,iaia]表示耦合系統(tǒng)S中，自由度a的準(zhǔn)線性驅(qū)動(dòng)點(diǎn)頻率響應(yīng)函數(shù)；[S,iaca]表示耦合系統(tǒng)S中，從自由度a到自由度a的準(zhǔn)線性跨點(diǎn)頻率響應(yīng)函數(shù)。

另一方面，根據(jù)力的平衡條件，可得：

式中：[B,cbcb]為線性子結(jié)構(gòu)B中，耦合自由度b的驅(qū)動(dòng)點(diǎn)頻率響應(yīng)函數(shù)。

結(jié)合式（12）、（13）、（14）可得到二級(jí)耦合系統(tǒng)的解耦公式：

式（15）在本研究之前已有其他研究者得出，其含義是非線性子結(jié)構(gòu)A的準(zhǔn)線性頻率響應(yīng)函數(shù)可通過對(duì)耦合系統(tǒng)S進(jìn)行測量得出。

與單點(diǎn)耦合系統(tǒng)相比，多點(diǎn)耦合系統(tǒng)中的耦合坐標(biāo)有個(gè)，因此，頻率響應(yīng)函數(shù)矩陣中對(duì)應(yīng)位置的元素為維矩陣。以兩點(diǎn)耦合系統(tǒng)（耦合自由度為a1和a2）為例，其耦合自由度的準(zhǔn)線性頻率響應(yīng)函數(shù)矩陣表示為：

1.3 基于模態(tài)參數(shù)識(shí)別的響應(yīng)特性求解

因?yàn)榉蔷€性系統(tǒng)的準(zhǔn)線性頻率響應(yīng)函數(shù)是在唯一固定的響應(yīng)幅值下確定的，因此，由式（15）只能得到非線性子結(jié)構(gòu)在特定振動(dòng)幅值下的響應(yīng)特性，而無法預(yù)測其在任意大小激勵(lì)下的響應(yīng)特性。對(duì)此，提出一種基于模態(tài)參數(shù)識(shí)別的求解策略。

首先，使用激振器給系統(tǒng)施加種不同幅值的激勵(lì)，記錄非線性子結(jié)構(gòu)中某點(diǎn)的響應(yīng)幅值，測得系統(tǒng)在種不同響應(yīng)幅值下的準(zhǔn)線性頻率響應(yīng)函數(shù)S(1),S(2), …,S()。之后，應(yīng)用模態(tài)參數(shù)識(shí)別技術(shù)求解出對(duì)應(yīng)于每個(gè)響應(yīng)幅值的固有頻率和阻尼比。經(jīng)此過程，會(huì)得到與響應(yīng)幅值有關(guān)的一組固有頻率{n1,n2, ……,nm}和一組阻尼比{1,2, ……,}。最后，以作為自變量，分別擬合其與固有頻率n和阻尼比的關(guān)系，得到n()和()，兩者都是的函數(shù)。

對(duì)應(yīng)于黏性阻尼系統(tǒng)頻率響應(yīng)函數(shù)的模態(tài)表達(dá)式，將作為響應(yīng)幅值函數(shù)的模態(tài)參數(shù)代入到頻率響應(yīng)函數(shù)表達(dá)式中，可以重構(gòu)出非線性結(jié)構(gòu)關(guān)于響應(yīng)幅值的頻率響應(yīng)函數(shù)(,)，其表達(dá)式為：

式中：為所分析的模態(tài)數(shù)；nr為第階模態(tài)的固有頻率；ζ為第階模態(tài)的阻尼比；C()為第階模態(tài)頻率響應(yīng)函數(shù)表達(dá)式中的常數(shù)項(xiàng)。

根據(jù)頻率響應(yīng)函數(shù)的定義，有：

對(duì)式（18）進(jìn)行求解，即可求得測量點(diǎn)在任意大小激勵(lì)下的響應(yīng)。

2 模型驗(yàn)證

2.1 集總參數(shù)模型

該節(jié)中，參考一般運(yùn)載車輛的結(jié)構(gòu)，建立了一個(gè)含非線性部件的集總參數(shù)模型，見圖3。通過對(duì)這個(gè)模型的分析，來驗(yàn)證所提出解耦方法的可靠性。

該模型由非線性子結(jié)構(gòu)A和線性子結(jié)構(gòu)B兩部分組成。子結(jié)構(gòu)A代表具有非線性特性的包裝件，子結(jié)構(gòu)B代表運(yùn)載車輛。其中，為集中質(zhì)量，集中質(zhì)量之間存在線性剛度和線性阻尼，它們的下標(biāo)為不同位置參數(shù)的編號(hào)。在子結(jié)構(gòu)A中加入一個(gè)非線性彈簧NL，以體現(xiàn)包裝材料和包裝結(jié)構(gòu)的非線性，這使得子結(jié)構(gòu)A成為一個(gè)非線性子結(jié)構(gòu)，也使得整個(gè)系統(tǒng)成為一個(gè)非線性系統(tǒng)。2個(gè)子結(jié)構(gòu)之間通過彈簧4和5、阻尼4和5，形成兩點(diǎn)耦合。模型中參數(shù)的值見表1。

下面說明如何在不拆解結(jié)構(gòu)、只能在耦合狀態(tài)下激勵(lì)和測量的前提下，計(jì)算非線性子結(jié)構(gòu)A的振動(dòng)響應(yīng)特性。

假設(shè)非線性彈簧NL具有立方剛度的性質(zhì)，其非線性內(nèi)力L表示為

式中：kNL的值為1×107 N/m3。

表1 集總參數(shù)模型中參數(shù)的值

Tab.1 Parameters in the lumped parameter model

該系統(tǒng)的準(zhǔn)線性頻率響應(yīng)函數(shù)可由式（11）給出，其中：

控制1處的振幅從2 mm增加到20 mm，計(jì)算耦合自由度ca在不同振幅下的準(zhǔn)線性頻率響應(yīng)函數(shù)。以S,ca1ca1和S,ca1ca2為例，其部分響應(yīng)特性見圖4，可見該結(jié)構(gòu)的響應(yīng)特性會(huì)隨響應(yīng)幅值而改變，具有明顯的非線性特征。

類似地，分別計(jì)算出S,cacb、S,cbcb和S,cbca，將它們代入到式（14）中，可解耦得到非線性子結(jié)構(gòu)A的準(zhǔn)線性頻率響應(yīng)函數(shù)A,caca，見圖5?？梢钥闯觯咏Y(jié)構(gòu)A的一階固有頻率會(huì)隨響應(yīng)幅值的增加而遞增，表現(xiàn)出漸硬的特性。

為檢驗(yàn)該解耦方法的準(zhǔn)確性，以非耦合狀態(tài)下的子結(jié)構(gòu)A為對(duì)象，應(yīng)用諧波平衡法直接計(jì)算1處驅(qū)動(dòng)點(diǎn)頻率響應(yīng)函數(shù)在若干振幅下的實(shí)際值，將其與預(yù)測值進(jìn)行比較，結(jié)果見圖6，兩者保持高度一致。

圖4 子結(jié)構(gòu)A耦合自由度頻率響應(yīng)函數(shù)在不同響應(yīng)幅值下的測量結(jié)果

圖5 子結(jié)構(gòu)A耦合自由度頻率響應(yīng)函數(shù)在不同響應(yīng)幅值下的解耦結(jié)果

圖6 m1處驅(qū)動(dòng)點(diǎn)頻率響應(yīng)函數(shù)在不同響應(yīng)幅值下的預(yù)測值與實(shí)際值對(duì)比

基于非線性子結(jié)構(gòu)A的1處的驅(qū)動(dòng)點(diǎn)頻率響應(yīng)函數(shù)，應(yīng)用多項(xiàng)式擬合法識(shí)別其模態(tài)參數(shù)，并擬合模態(tài)參數(shù)與非線性彈簧響應(yīng)幅值之間的函數(shù)關(guān)系。將各模態(tài)參數(shù)作為振幅的函數(shù)代入到式（17）中，可以重構(gòu)出式（21）函數(shù)來描述的1處的振動(dòng)響應(yīng)特性：

假設(shè)外部激勵(lì)F為20 N，而激勵(lì)的位移幅值不再固定，作用于m1處對(duì)子結(jié)構(gòu)A進(jìn)行正向掃頻，根據(jù)式（18）所給出的關(guān)系，即可對(duì)位移響應(yīng)X進(jìn)行求解，這里使用的求解方法為牛頓迭代法。解得m1處的驅(qū)動(dòng)點(diǎn)頻率響應(yīng)函數(shù)，并將其與理論值進(jìn)行比較，結(jié)果見圖7?？梢妰烧吒叨纫恢?，且精準(zhǔn)預(yù)測了對(duì)非線性結(jié)構(gòu)進(jìn)行掃頻時(shí)出現(xiàn)的響應(yīng)跳躍現(xiàn)象。

2.2 有限元模型

在此節(jié)中，建立了一個(gè)包含非線性子結(jié)構(gòu)的多點(diǎn)耦合系統(tǒng)有限元模型，并通過對(duì)其進(jìn)行頻率響應(yīng)分析來檢驗(yàn)文中所提出理論在連續(xù)體模型中的應(yīng)用效果。該模型由一個(gè)具有線性剛度特性的懸臂梁和一個(gè)具有非線性剛度特性的T型梁通過2根彈簧連接而成，模型裝配效果及幾何參數(shù)見圖8。

圖8 包含非線性部件的多點(diǎn)耦合系統(tǒng)有限元模型及其幾何參數(shù)

文中所使用的有限元建模與分析軟件為ABAQUS。在軟件環(huán)境中完成建模后，設(shè)置如下材料參數(shù)與約束條件：構(gòu)成懸臂梁和T型梁的材料均為鋼，其彈性模量為200 GPa，泊松比為0.3，密度為7 800 kg/m3。上方懸臂梁為線性子結(jié)構(gòu)，在其固定端設(shè)置完全固定約束；下方的T形梁為非線性子結(jié)構(gòu)，由一根懸臂梁和一根兩端固定薄梁剛性連接而成，在連接位置設(shè)置綁定約束，在其固定端設(shè)置完全固定約束。點(diǎn)與點(diǎn)、點(diǎn)與點(diǎn)之間各有一根剛度均為5 000 N/m的彈簧連接，通過設(shè)置彈簧連接器實(shí)現(xiàn)。

應(yīng)用文中所提出的解耦策略，可以在不對(duì)模型進(jìn)行拆解的前提下，通過測量該系統(tǒng)在4個(gè)耦合點(diǎn)的頻率響應(yīng)函數(shù)，來預(yù)測非線性子結(jié)構(gòu)（即T型梁）的頻率響應(yīng)函數(shù)。利用有限元軟件對(duì)該模型進(jìn)行模態(tài)分析，得到其一階模態(tài)的有效質(zhì)量占結(jié)構(gòu)實(shí)際質(zhì)量的83.2%，因此可以判斷該模型的響應(yīng)主要受一階模態(tài)影響，文中考慮其一階模態(tài)。

對(duì)系統(tǒng)進(jìn)行激振，控制點(diǎn)的響應(yīng)幅值分別為0.5、1、1.5、2、2.5、3 mm，測量每個(gè)響應(yīng)幅值下各耦合點(diǎn)的準(zhǔn)線性頻率響應(yīng)函數(shù)，它們組成了系統(tǒng)耦合自由度的準(zhǔn)線性頻率響應(yīng)函數(shù)矩陣。將它們代入到式（15）中，即可求得非線性子結(jié)構(gòu)在不同響應(yīng)幅值下的準(zhǔn)線性頻率響應(yīng)函數(shù)。以點(diǎn)為例，其驅(qū)動(dòng)點(diǎn)頻率響應(yīng)函數(shù)的求解結(jié)果見圖9。對(duì)該頻率響應(yīng)函數(shù)進(jìn)行參數(shù)識(shí)別，并擬合模態(tài)參數(shù)與響應(yīng)幅值的關(guān)系。

圖9 a點(diǎn)的驅(qū)動(dòng)點(diǎn)頻率響應(yīng)函數(shù)在不同響應(yīng)幅值下的解耦結(jié)果

只要給定激勵(lì)的幅值，即可代入式（18）進(jìn)行求解。為驗(yàn)證求解結(jié)果，使用ABAQUS顯式動(dòng)力學(xué)模塊對(duì)激勵(lì)為1 N時(shí)若干頻率點(diǎn)的響應(yīng)進(jìn)行計(jì)算，并將計(jì)算值與預(yù)測值進(jìn)行對(duì)比，結(jié)果見圖10，兩者保持高度一致。

圖10 在a點(diǎn)進(jìn)行掃頻激勵(lì)的頻率響應(yīng)函數(shù)預(yù)測值實(shí)際值對(duì)比

3 結(jié)語

考慮到包裝結(jié)構(gòu)和包裝材料的非線性與運(yùn)輸包裝系統(tǒng)耦合形式的復(fù)雜性，文中提出了一種針對(duì)非線性多點(diǎn)耦合系統(tǒng)的解耦方法，能夠在不拆解系統(tǒng)的前提下，通過測量運(yùn)輸包裝系統(tǒng)的頻率響應(yīng)函數(shù)來求解包裝件的動(dòng)態(tài)響應(yīng)特性，具體步驟如下。

1）使用激振器對(duì)系統(tǒng)施加簡諧激勵(lì)，測量系統(tǒng)耦合位置在不同響應(yīng)幅值下的準(zhǔn)線性頻率響應(yīng)函數(shù)。

2）使用逆向子結(jié)構(gòu)方法求解非線性子結(jié)構(gòu)在不同響應(yīng)幅值下的準(zhǔn)線性頻率響應(yīng)函數(shù)。

3）識(shí)別非線性子結(jié)構(gòu)的模態(tài)參數(shù)并擬合其模態(tài)參數(shù)與響應(yīng)幅值之間的關(guān)系。

4）構(gòu)建激勵(lì)幅值與響應(yīng)幅值之間的關(guān)系方程并進(jìn)行求解。

在集總參數(shù)模型和有限元連續(xù)模型中，文中方法對(duì)非線性子結(jié)構(gòu)振動(dòng)響應(yīng)特性的預(yù)測結(jié)果與實(shí)際值基本吻合，證明了其有效性，在運(yùn)輸包裝系統(tǒng)參數(shù)獲取與模型構(gòu)建、運(yùn)輸過程產(chǎn)品狀態(tài)監(jiān)測等方面具備應(yīng)用前景。

[1] WANG zhi‐wei, ZHONG lin‐lin. Finite Element Analysis and Experimental Investigation of Beer Bottle‐Turnover Boxes Transport Unit under Random Vibration Excitation[J]. Packaging Technology and Science, 2020, 33(6): 197-214.

[2] 高德, 吳朝武, 陸俊杰. 運(yùn)輸包裝系統(tǒng)振動(dòng)行為研究與發(fā)展趨勢分析[J]. 包裝工程, 2020, 41(15): 51-58.

GAO De, WU Chao-wu, LU Jun-jie. Dynamic Behavior Research and Dynamic Development of Transportation Packaging System[J]. Packaging Engineering, 2020, 41(15): 51-58.

[3] WANG Zhi-wei, ZHANG Yuan-biao. Dynamic Characteristic Analysis of Refrigerator-Truck Transport System by Using Inverse Substructure Method[J]. Packaging Technology and Science, 2014, 27(11): 883-900.

[4] ZHANG Yuan-biao, WANG Zhi-wei. Investigation of Frequency Response Function of Product-Transport System Based on Multi-Coordinate Coupled Inverse Substructure Method[J]. Packaging Technology and Science, 2014, 27(5): 364-375.

[5] WANG Jun, SUN Guo-hua, LU Li-xin, et al. Indirect Inverse Substructuring Theory for Coupling Dynamic Stiffness Identification of Complex Interface between Packaged Product and Vehicle Transport System[J]. Packaging Technology and Science, 2015, 28(2): 141-155.

[6] WANG Jun, WANG Qi-li, LU Li-xin, et al. Inverse Sub-Structuring Theory for Coupled Product Transport System Based on the Dummy Masses Method[J]. Packaging Technology and Science, 2016, 29(3): 189-200.

[7] WANG Jun, HONG Xiang, QIAN Yi, et al. Inverse Sub-Structuring Method for Multi-Coordinate Coupled Product Transport System[J]. Packaging Technology and Science, 2014, 27(5): 385-408.

[8] MENG Tian-ya, WANG Jun, PU Guang-yi, et al. Inverse Sub-Structuring Method for Multi-Coordinate Rigidly Coupled Product Transport System Based on a Novel Shearing Probe Technique[J]. Packaging Technology and Science, 2017, 30(9): 601-618.

[9] WANG Jun, WANG Zhi-wei, LU Li-xin. Step-by-Step Decoupling Method for Inverse Substructuring Analysis of a Three-Component Coupled Packaging System[J]. Journal of Vibration and Control, 2015, 21(4): 676-683.

[10] 傅苗苗, 王軍, 盧立新, 等. 多重門限奇異值分解方法的實(shí)例驗(yàn)證[J]. 包裝工程, 2021, 42(21): 176-180.

FU Miao-miao, WANG Jun, LU Li-xin, et al. Example Verification of Multi-Threshold Singular Value Decomposition Method[J]. Packaging Engineering, 2021, 42(21): 176-180.

[11] 傅苗苗, 王軍, 盧立新, 等. 基于頻率響應(yīng)函數(shù)的逆子結(jié)構(gòu)方法誤差分析[J]. 包裝工程, 2021, 42(23): 141-145

FU Miao-miao, WANG Jun, LU Li-xin, et al. Error Analysis of the Inverse-Substructuring Method Based on Frequency Response Function[J]. Packaging Engineering, 2021, 42(23): 141-145

[12] JOODAKY A, BATT G S, GIBERT J M. Prediction of Cushion Curves of Polymer Foams Using a Nonlinear Distributed Parameter Model[J]. Packaging Technology and Science, 2020, 33(1): 3-14.

[13] YANG Song-ping, WANG Zhi-wei. Acceleration Spectrum Analysis of Hyperbolic Tangent Package under Random Excitation[J]. Packaging Technology and Science, 2021, 34(9): 579-587.

[14] 應(yīng)玉萍, 王花蘭. 基于幾何非線性的果蔬運(yùn)輸隔振系統(tǒng)設(shè)計(jì)[J]. 包裝工程, 2020, 41(19): 153-158.

YING Yu-ping, WANG Hua-lan. Design of Vibration Isolation System for Fruit and Vegetable Transportation Based on Geometric Nonlinearity[J]. Packaging Engineering, 2020, 41(19): 153-158.

[15] WANG Peng-tao, WANG Jun, LIM T C, et al. A Strategy for Decoupling of Nonlinear Systems Using the Inverse Sub-Structuring Method and the Parametric Modal Identification Technique[J]. Mechanical Systems and Signal Processing, 2020, 140: 106695.

[16] KARAA?A?L? T, OZGUVEN H N. Experimental Modal Analysis of Nonlinear Systems by Using Response-Controlled Stepped-Sine Testing[J]. Mechanical Systems and Signal Processing, 2021, 146: 107023.

A Decoupling Method for Multi-coordinate Coupled Transport Packaging System Containing Nonlinear Sub-structure

WANG Wei-kai1,WANG Jun1,LU Li-xin1,PAN Liao1,HOU Xue2

(1. Jiangnan University, Wuxi 214122, China; 2. Wuxi Branch of Shantou Dongfeng Printing Co., Ltd., Wuxi 214000, China)

The work aims to optimize the traditional inverse sub-structuring method, to get the dynamic response characteristic of sub-structure in nonlinear multi-coordinate coupled system, which is helpful to resolve the issue that complex coupling forms and nonlinear characteristic of material make it difficult to measure the physical parameters in transport packaging system. The description function method is used to linearize the nonlinear system and the frequency response functions (FRFs) at several response amplitudes need to be measured. Then, the modal parameters of the sub-structure can be identified by inverse sub-structuring method. Lastly, the relationship between the modal parameters and vibration amplitudes is fitted and a function is constructed to describe the dynamic response characteristic of the sub-structure. In the lumped parameter model, the predicted response was consistent with the actual. In the finite element model, the prediction error of the response peak was less than 5%, and the prediction error of the jumping frequency was less than 3%. The application of traditional inverse sub-structuring method was extended to nonlinear multi-coordinate coupled system, which had guiding significance for the construction of dynamic model of complex transport packaging system and the design of anti-vibration packaging.

inverse sub-structuring method; transport packaging; vibration; nonlinear

TB485.3

A

1001-3563(2022)23-0252-07

10.19554/j.cnki.1001-3563.2022.23.030

2022?02?28

國家一流學(xué)科建設(shè)輕工技術(shù)與工程（LITE 2018–29）

王維凱（1997—），男，碩士生，主攻運(yùn)輸包裝。

王軍（1982—），男，博士，教授，主要研究方向?yàn)檫\(yùn)輸包裝、包裝新材料。

責(zé)任編輯：曾鈺嬋