1:1內(nèi)共振情況下輕質(zhì)材料層合板動力學(xué)的奇異性分析*

郭宇紅 張 偉 楊曉東

(北京工業(yè)大學(xué) 機電學(xué)院, 北京 100124)

引言

隨著高新技術(shù)的發(fā)展,人們已不再滿足于材料單純的輕質(zhì)化,而是尋找兼有輕質(zhì)化和其他某種或幾種優(yōu)良性能相結(jié)合的先進材料以適應(yīng)不同的需求,作為一種生物材料中普遍存在的典型結(jié)構(gòu),多孔介質(zhì)在自然界動植物體中發(fā)揮著不可替代的生理功能. 在這些構(gòu)型中,3D-Kagome結(jié)構(gòu)擁有更加稀少的桿件布局和更小的相對密度.3D-Kagome結(jié)構(gòu)是由兩個正四面體對頂連接而成的網(wǎng)架結(jié)構(gòu),將這種核心連接在面板上就形成了3D-Kagome點陣夾芯板[1-3]. 通過實驗和數(shù)值模擬得到結(jié)論,相對于其他構(gòu)型的點陣夾芯板,在同一相對密度情況下,3D-Kagome點陣夾芯板具有更高的強度和抗屈曲性能[4-6]. 此外,3D-Kagome點陣夾芯板在功能性方面的表現(xiàn)也比較突出,例如具有優(yōu)異的驅(qū)動和致動性能,能夠在受到很小的內(nèi)部抵抗力的情況下獲得很大范圍內(nèi)的整體變形. 點陣材料在航空航天等領(lǐng)域有著廣泛的應(yīng)用,可以減輕飛行器的重量,同時保證結(jié)構(gòu)的強度和剛度滿足要求.

自從上世紀(jì)70年代,Golubitsky、Schaeffer和Stewart等將奇異性理論和群論方法引入分叉問題的研究后,分叉理論得到了越來越多的關(guān)注,從而也推動了奇異性理論的發(fā)展.1981年,Golubitsky和Langford[7]研究了退化Hopf分叉問題的分類和開折.1982年,Golubitsky[8]和Martinet[9]等討論了在強等價下光滑映射芽的開折,給出了各種形式的通有開折定理.1985年,Golubitsky和Schaeffer[10]得到了余維數(shù)不大于3的分叉問題的分類,其中狀態(tài)變量只有一個且具有Z2對稱性. 1986年,Keyfitz[11]給出了余維數(shù)不大于7的分叉問題的分類,其中狀態(tài)變量只有一個但不具對稱性,給出了各種形式的通有開折定理,但是這些研究都是針對單狀態(tài)變量、單分叉參數(shù)的分叉系統(tǒng)而言. 隨后人們將奇異性理論推廣到了多狀態(tài)變量的分叉研究中,1986年,Golubitsky[12]研究了單參數(shù)兩狀態(tài)變量余維數(shù)不超過2的分叉問題的分類,其中狀態(tài)變量關(guān)于二面體群4D對稱,1988年,Melbourne[13]得到了單參數(shù)三個狀態(tài)變量余維數(shù)不超過1的分叉問題的分類.

在開始的研究中,研究工作沒有考慮分叉參數(shù)的對稱性. 1986年,Gaffney[14]將冪單代數(shù)群和冪零Lie代數(shù)應(yīng)用于多參數(shù)分叉問題中,給出了D(Γ)-等價. 1996年,Furter和Sitta[15]等考慮了分叉參數(shù)的對稱性,研究了余維數(shù)不大于1的分叉問題的分類.2003年,高守平和李養(yǎng)成等[16]討論了狀態(tài)變量和分叉參數(shù)具有不同對稱性的余維數(shù)不超過1的分叉問題的分類,其中狀態(tài)變量關(guān)于二面體群D4對稱,分叉參數(shù)關(guān)于S1對稱. 2006年,郭瑞芝[17]給出了狀態(tài)變量和分叉參數(shù)具不同對稱性的岔問題的分類及識別條件,其中狀態(tài)變量關(guān)于二面體群D3對稱,分叉參數(shù)關(guān)于O(2)對稱. 2007年,崔登蘭和李養(yǎng)成等[18]也研究了含有兩組狀態(tài)變量且參數(shù)具有對稱性的等變分叉問題.

隨著研究的進一步深入,開始了多分叉參數(shù)的研究工作. 1993年,Lavassani等[19]通過奇異性理論研究了等變多參數(shù)分叉,給出了多參數(shù)分叉問題的有限確定定理和正規(guī)型,并討論了多參數(shù)分叉的穩(wěn)定性問題. 2000年,胡凡努和李養(yǎng)成[20]將狀態(tài)變量分為兩組,一組狀態(tài)變量可以獨立變化,而另一組狀態(tài)變量則依賴于前一組狀態(tài)變量,研究了該類分叉問題的通有開折.2003年,高守平和李養(yǎng)成[21]研究了多參數(shù)等變分歧問題及其開折,給出了通有的開折定理. 但是人們通常將分歧(分叉)問題中的狀態(tài)變量看作是“平等”的,并不加以區(qū)分. 2005年,郭瑞芝和李養(yǎng)成[22]研究了含有兩組狀態(tài)變量的多參數(shù)等變分歧問題在左右等價群下的開折,得出了通有開折的充要條件. 2010年,秦朝紅和陳予恕等[23-26]研究了含有兩個狀態(tài)變量和兩個分叉參數(shù)的分叉系統(tǒng)的奇異性理論,并給出了含有兩個分叉參數(shù)系統(tǒng)的轉(zhuǎn)遷集的計算方法.

本文以點陣夾芯板為切入點,計算1:1內(nèi)共振情況下點陣夾芯板的非線性動力學(xué)方程,推廣含有兩個狀態(tài)變量和三個及三個以上參數(shù)的一般非線性動力學(xué)分叉方程的奇異性理論,利用推廣的奇異性理論計算復(fù)合層合板的非線性動力學(xué)分叉方程的普適開折,同時計算普適開折的轉(zhuǎn)遷集,并進一步討論若干重要參數(shù)對穩(wěn)定性的影響.

1點陣夾芯板非線性動力學(xué)分叉方程的普適開折

1.1 點陣夾芯板的非線性動力學(xué)方程

考慮橫向激勵與面內(nèi)激勵聯(lián)合作用下四邊簡支點陣夾芯板,在板的中面上建立坐標(biāo)系oxy,板在x和y方向的長度分別為a和b,厚度為h,設(shè)夾芯板中面上任一點在x、y和z方向的位移分別為u、v和w,板承受沿z方向的橫向激勵f=F(x,y)cosΩ1t與沿y方向作用于x=0和x=a的面內(nèi)激勵p=p0+p1cosΩ2t聯(lián)合作用,這里,Ω1和Ω2分別為橫向激勵與面內(nèi)激勵的頻率. 點陣夾芯板由點陣夾芯層與覆蓋在上面的兩層蒙皮組成,蒙皮采用各向同性材料,其厚度為hf,芯層厚度為hc. 其動力學(xué)方程為:

(1a)

(1b)

其中,w1為第一階模態(tài)的振幅,w2為第二階模態(tài)的振幅,μ1和μ2表示對應(yīng)于兩階模態(tài)的阻尼,F1和F2表示對應(yīng)于兩階模態(tài)的橫向激勵幅值,βij(i=1,2;j=1…7)表示材料參數(shù). 方程(1)的詳細推導(dǎo)過程及其所有參數(shù)表達式參見文獻[27].

1.2 在1:1內(nèi)共振情況下點陣夾芯板動力學(xué)方程攝動分析

利用多尺度法進行研究,將方程(1)中的阻尼項、參數(shù)激勵項、熱激勵項和非線性項添加小擾動項ε,考慮點陣夾芯板的主參數(shù)共振-1:1內(nèi)共振的情況,共振關(guān)系如下:

2ω1=Ω1-εσ1,2ω2=Ω2-εσ2, Ω1=Ω2

(2)

式中ω1和ω2為相應(yīng)線性系統(tǒng)的第一階和第二階固有頻率,σ1和σ2為系統(tǒng)的調(diào)諧參數(shù),為了方便處理,令Ω1=1.

設(shè)方程(1)的一致漸近解為:

w(x,t,ε)=w0(x,T0,T1)+εw1(x,T0,T1)

(3)

其中T0=t,T1=εt.

則有微分算子:

(4a)

(4b)

將式(3)和(4)帶入方程(1)中,比較方程兩邊攝動參數(shù)ε同階次的系數(shù),得到如下方程:

ε0階:

(5a)

(5b)

ε1階:

(6a)

(6b)

方程(5a)和(5b)的解可以寫成如下復(fù)數(shù)形式:

(7)

將方程(7)代入方程(6a)和方程(6b)中得到:

(8a)

(8b)

其中cc和NST分別表示方程(8)右端函數(shù)的復(fù)數(shù)部分和長期項.

A1和A2可以表示為下列形式:

(9)

消除方程(8)中可以產(chǎn)生長期項的部分,將方程(9)代入方程(8), 將實部與虛部分離,得到極坐標(biāo)形式的四維平均方程為:

(10a)

(10b)

(10c)

(10d)

(11b)

展開方程(11),得到:

(12a)

(12b)

式中:

k19=-β13cos(φ1-φ2),k110=-β13β16cos(3φ1-φ2),

k111=-3β14β16cos(2φ1),k112=-3β14,

k116=2β16sin(2φ1),k117=2β16cos(2φ1),

k29=-β23cos(φ2-φ1),k210=-β23β26cos(3φ2-φ1),

k211=-3β24β26cos(2φ2),k212=-3β24,

k216=2β26sin(2φ2),k217=2β26cos(2φ2)

(12c)

1.3 在1:1內(nèi)共振情況下點陣夾芯板動力學(xué)方程分叉分析

令:

(13)

式中:

(14a)

(14b)

z=(a1,a2),λ=(σ1,σ2,p1)

(14c)

1.3.1限制切空間

定理1.1芽g(z,λ)的限制切空間RT(g,1)能夠表示為:

RT(g,1)=M2+M〈σ1,σ2,p1〉

(15)

證明:根據(jù)命題1.4(文獻[10]中,第二冊169頁),RT(g,1)的生成元有14個:

(g1,0),(g2,0),(0,g1),(0,g2),a1(g1,a1,g2,a1),

a2(g1,a1,g2,a1),σ1(g1,a1,g2,a1),σ2(g1,a1,g2,a1),

p1(g1,a1,g2,a1),a1(g1,a2,g2,a2),a2(g1,a2,g2,a2),

σ1(g1,a2,g2,a2),σ2(g1,a2,g2,a2),p1(g1,a2,g2,a2)

(16)

式中:

(17a)

(17b)

(17c)

(17d)

M2+M〈σ1,σ2,p1〉的生成元有18個:

(σ2a1,0), (p1a1,0), (σ1a2,0), (σ2a2,0),

(0,σ1a1), (0,σ2a1), (0,p1a1), (0,σ1a2),

(0,σ2a2), (0,p1a2)

(18)

下面證明,當(dāng)(a1,a2,σ1,σ2,p1)=(0,0,0,0,0)時,(16)和(18)之間存在一個可逆矩陣A,使得兩者可以互相表示:

(19)

式中:

(20a)

u2=(σ2(g1,a1,g2,a1)p1(g1,a1,g2,a1)a1(g1,a2,g2,a2)

a2(g1,a2,g2,a2)σ1(g1,a2,g2,a2)σ2(g1,a2,g2,a2)

(20b)

(σ2a1,0) (p1a1,0) (σ1a2,0) (σ2a2,0)

(20c)

(0,σ2a1) (0,p1a1) (0,σ1a2) (0,σ2a2)

(20d)

(20e)

(20g)

(20h)

式中:

(21a)

(21b)

(21c)

(21d)

(21e)

k28μ2a1a2+k29σ2a1a2+k210p1a1a2+

(21f)

(21g)

2k213σ2a1a2+2k214p1a1a2

(21h)

2k113σ1a1a2+2k114p1a1a2

(21i)

(21j)

因為μ1>0,μ2>0,當(dāng)a1=a2=σ1=σ2=p1=0時,

(22)

這樣矩陣A等于:

(23)

式中:

(24a)

(24b)

(24c)

(24d)

在14×18的矩陣A中,容易發(fā)現(xiàn)第1行和第5行線性相關(guān),第4行和第11行線性相關(guān),第7、8、9、13、14、15列是空的,剔除第5、11行和第7、8、9、13、14、15列,剩下一個12×12的矩陣,很容易證明12×12的矩陣的行列式的值是非零的.

證明完成.

1.3.2 簡單識別

定理1.2令:

g(z,λ)=h(z,λ)+q(z,λ)

(25)

式中:

(26a)

(26b)

可以得到g和h是強等價的.

證明:根據(jù)定理1.1,M2+M〈σ1,σ2,p1〉的一個高階項是:

M3+M2〈σ1,σ2,p1〉

(27)

對于多項式g,我們觀察發(fā)現(xiàn):

(28a)

(28b)

根據(jù)定理4.1(文獻[10]中,第二冊185頁),可以得到g和h是強等價的.

證明完成.

在下面的討論中,g將被h取代.

引理1.1h的非退化條件是:

(29)

滿足非退化條件的h等價于:

(30)

證明:h是關(guān)于a1,a2的二次齊次多項式,根據(jù)方程(2.7) (文獻[10]中,第一冊402頁),

(31)

根據(jù)方程(2.8)(文獻[10]中,第一冊402頁),

(32)

因為μ1>0,μ2>0,

(33)

h的正規(guī)形式能夠表示為:

(34)

證明完成.

觀察表達式(34),h的普適開折需要表達式(34)補足h的線性項和常數(shù)項.

定理1.3令:

H(a1,a2,σ1,σ2,p1)=(h1(a1,a2,σ1,σ2,p1),h2(a1,a2,σ1,σ2,p1))是分叉問題h的一個4參數(shù)開折,如果滿足引理1.1,則H是h的普適開折,當(dāng)且僅當(dāng):

det(Q)=

(35)

式中γ是輔助參數(shù),這里:

(a1,a2,σ1,σ2,p1,γ)=(0,0,0,0,0,0)

證明: 矩陣Q能夠被表示為:

Q=(α1,α2,α3,α4,α5,α6)T

(36)

式中:

α1=(0,h1,a1a1,h1,a1a2,0,h2,a1a1,h2,a1a2)

(37a)

α2=(0,h1,a2a1,h1,a2a2,0,h2,a2a1,h2,a2a2)

(37b)

α3=(h1,σ1,h1,σ1a1,h1,σ1a2,h2,σ1,h2,σ1a1,h2,σ1a2)

(37c)

α4=(h1,σ2,h1,σ2a1,h1,σ2a2,h2,σ2,h2,σ2a1,h2,σ2a2)

(37d)

α5=(h1,p1,h1,p1a1,h1,p1a2,h2,p1,h2,p1a1,h2,p1a2)

(37e)

α6=(h1,γ,h1,γa1,h1,γa2,h2,γ,h2,γa1,h2,γa2)

(37f)

將h代入方程(37), 得到:

α1=(0,2,0,0,0,0),α2=(0,0,0,0,0,2),

α3=(0,0,0,0,0,0),α4=(0,0,0,0,0,0),

α5=(0,0,0,0,0,0),α6=(0,0,0,0,0,0)

(38)

在方程(38),僅有2個向量α1和α2是線性無關(guān)的,因此存在補足h的4個線性無關(guān)的向量:

(39)

式中λ,λ′和λ″的取值是σ1,σ2和p1三個分叉參數(shù)中的任意一個,且互相之間取值不同.

這樣,我們得到:

(40)

把方程(40)代入方程(35),得到:

det(Q)≠0.

(41)

反過來,在方程(38)中,僅有向量α1和α2是線性無關(guān),因為det(Q)≠0, 補足h的線性無關(guān)的向量需要4個:

(42)

這樣,得到:

證明完成.

1.3.3 多項式空間的維數(shù)

定理1.4多項式空間能夠簡化為:

R{(dh)z,λ(Y1),…,(dh)z,λ(Ym),hλ,λhλ,λ2hλ,…}

=R{(a1,0),(0,a2)}

(43)

證明: 推導(dǎo)方程(43),建立下面的矩陣關(guān)系:

(44)

式中:

(45a)

(45b)

(45c)

(45d)

(45e)

(45f)

(45g)

(45h)

顯然:

D1=2,D2=2

(46)

R{(dh)z,λ(Y1),…,(dh)z,λ(Ym),hλ,λhλ,λ2hλ,…}

=R{(a1,0),(0,a2)}

(47)

證明完成.

1.3.4 奇異性理論推廣

根據(jù)定理1.1,有下面的關(guān)系:

h∈M2?M2+M〈σ1,σ2,p1〉

(48)

得到:

(49)

因此,對奇異性理論作如下的推廣:

〈h〉⊕p(n)=T(h+1)

(50)

式中p(n)是一個有限維子空間.

同理,存在m個線性無關(guān)的向量補足T(h,1),使得

(51)

式中p′(m)是一個有限維子空間.

把方程(50)代入方程(51), 得到

(52)

證明完成.

推論1.1n=6-2-m.

證明:根據(jù)定理1.3,存在6個線性無關(guān)的向量,但是在方程(38)中,僅有兩個線性無關(guān)的向量α1和α2. 根據(jù)引理1.1, 可知2+m+n=6,即:n=6-2-m.

證明完成.

引理1.3讓輔助參數(shù)γ∈R作用在單項式λl1xl2上, 如果l1≠0,則|γ|能夠被剔除,sgn(γ)λl1xl2將被獲得.

證明:對于γλl1xl2,如果l1≠0,|γ|能夠嵌入在λl1中,得到sgn(γ)λl1,這是因為γ和λ都是參數(shù),有共同的輔助特性.

證明完成.

在后面的討論中,為了簡化符號,sgn(γ)被記作ε.

1.3.5 普適開折

ε11p1+γ)

(53a)

ε21p1+γ)

(53b)

ε31σ1+γ)

(53c)

ε41σ1+γ)

(53d)

ε51σ2+γ)

(53e)

ε61σ2+γ)

(53f)

式中γ是輔助參數(shù),εij=+1,0,-1；i=1,…,6；j=1,2,3.

證明:根據(jù)方程(2.7) (文獻[10]中,第二冊211頁),可知:

R{(a1,0),(0,a2)}

(54)

ItrT(h,1)是包含在T(h,1)中的極大理想:

(55)

可知:

(σ2,0),(p1,0),(-1,0)

(0,a1),(0,a2),(0,σ1),

(0,σ2),(0,p1),(0,1)}

(56)

維數(shù)是12.

(a2,0), (0,a1), (σ1,0), (0,σ1), (σ2,0),

(0,σ2),(p1,0), (0,p1),(-1,0),(0,1)

(57)

根據(jù)定理1.3,簡化方程(57),

(a2,0), (0,a1), (σ1,σ1), (σ2,σ2),

(p1,p1), (-1,1)

(58)

(σ1a1,0), (σ2a1,0), (p1a1,0), (σ1a2,0),

(σ2a2,0), (p1a2,0), (0,σ1a1), (0,σ2a1),

(0,p1a1), (0,σ1a2), (0,σ2a2), (0,p1a2)

(59)

在方程(58)中,對于向量(a2,0)和(0,a1),存在兩個輔助參數(shù)β1和β2,使得(β1a2,0)和(0,β2a1)發(fā)生. 根據(jù)定理1. 3,(β1a2,0)和(σ1a2,0), (σ2a2,0), (p1a2,0)三個向量中的任意一個線性相關(guān)；(0,β2a1)和(0,σ1a1), (0,σ2a1), (0,p1a1)三個向量中的任意一個線性相關(guān). 因此,(a2,0)和(0,a1)將被剔除.

在h的普適開折中,存在4個輔助參數(shù)α1,α2,α3和γ,使得(α1p1,α1p1),(α2σ1,α2σ1), (α3σ2,α3σ2)和(-γ,γ)發(fā)生,根據(jù)引理1.3,(ε1p1,ε1p1),(ε2σ1,ε2σ1),(ε3σ2,ε3σ2)和(-γ,γ)發(fā)生. 因為p1,σ1,σ2和γ都是參數(shù),根據(jù)定理1.3,把四個向量(ε1p1,ε1p1),(ε2σ1,ε2σ1),(ε3σ2,ε3σ2)和(-γ,γ)代入方程(37)中,僅有兩個向量是線性無關(guān)的,根據(jù)定義1.1,m=2.

根據(jù)定理1.3, 在矩陣Q中,可知向量(-γ,γ)必然發(fā)生.

為了簡化符號,(ε2σ1,ε2σ1)記作(ε1σ1,ε1σ1),(ε3σ2,ε3σ2)記作(ε1σ2,ε1σ2).

根據(jù)推論1.1,n=2.

根據(jù)上面的分析,得到:

ε11p1+γ)

(60a)

ε21p1+γ)

(60b)

ε31σ1+γ)

(60c)

ε41σ1+γ)

(60d)

ε51σ2+γ)

(60e)

ε61σ2+γ)

(60f)

根據(jù)引理1.3,當(dāng)所有隱輔助參數(shù)和輔助參數(shù)γ等于零,即:εij=0 (i=1,…,6;j=1,2,3)和γ=0,我們可以得到:

Gi(a1,a2,p1,0,0,0,0)=h(a1,a2,p1)

(61)

根據(jù)定理1.3,當(dāng)εij≠0,對于方程(60)有:

(62a)

(62b)

(62c)

(62d)

(62e)

(62f)

證明完成.

1.3.6 轉(zhuǎn)遷集

在下面的分析中,我們討論方程(53)的轉(zhuǎn)遷集.

對于方程(53a),分叉集滿足的條件是:

(63a)

(63b)

(63c)

由方程(63c)的第三個等式,可得分叉集的表達式:

(64)

方程(63c)的第一和第二個等式,表明分叉發(fā)生時:

(a1,a2,σ1,σ2,p1,γ)=(0,0,0,0,0,0)

(65)

滯后集:

G=0, det(dG)=0和d2G(v,v)∈range(dG),

(66)

任意非零v∈ker(dG).

滯后集滿足的條件:

(67a)

(67b)

a1a2-ε12ε13σ1σ2=0

(67c)

(67d)

式中v=(v1,v2).

我們假設(shè)σ1≠0和σ2≠0,則:

v=(a1,ε12σ1)

(68)

(69)

計算方程(69),可得:

(70)

讓a1乘以方程(70)的兩端,然后代入方程(67c)中,得到:

(71)

由方程(67c)和方程(71),可得:

(72)

方程(67a)加方程(67b),可得:

(73)

把方程(71)和方程(72)代入方程(73)中,得到滯后集的表達式:

(74)

雙極限點集滿足的條件:

(75a)

(75b)

a1a2-ε12ε13σ1σ2=0

(75c)

(a11,a21)≠(a12,a22),當(dāng)p1=const

(75d)

根據(jù)方程(64)計算,可知:

σ1=0,p1≤0和σ2=0,p1≤0當(dāng)ε11=1,

(76a)

σ1=0,p1≥0和σ2=0,p1≥0當(dāng)ε11=-1

(76b)

我們定義ρ,θ和δ是fTσ2σ1-空間的圓柱坐標(biāo)的轉(zhuǎn)換,令:

σ1=ρ3cos3θ,σ2=ρ3sin3θ,p1=ρ6δ

(77)

把方程(77)代入方程(64)和方程(74)中,B,H和D有如下形式:

(78a)

(78b)

(78c)

同理,方程(53b),(53c), (53d),(53e)和(53f)的轉(zhuǎn)遷集∑也能夠被獲得.

2 數(shù)值模擬

本節(jié)利用Runge-Kutta法對方程(53a)進行數(shù)值計算,針對分叉表達式,滯后表達式和雙極限點集得到關(guān)于三個分叉參數(shù)σ1,σ2和p1的平面轉(zhuǎn)遷集和立體轉(zhuǎn)遷集,同時給出了相應(yīng)的力-幅圖.

根據(jù)以上分析結(jié)果,對方程(53a)進行數(shù)值計算,由方程(64)和方程(74),我們得到方程(53a)關(guān)于三個分叉參數(shù)σ1,σ2和p1的平面轉(zhuǎn)遷集和立體轉(zhuǎn)遷集,當(dāng)ε11=-1時,如圖1和圖2所示；當(dāng)ε11=1時,如圖5和圖6所示,B表示分叉集,H表示滯后集,D表示雙極限點集. 從方程(64)的表達式可以發(fā)現(xiàn),當(dāng)p1固定時,在圖2(a)中,分叉集是一個圓；當(dāng)σ1或σ2固定時,在圖2(a)中,分叉集是一個拋物線.

圖1和圖2將方程(53a)平衡點附近鄰域分為不同的區(qū)域,分別對應(yīng)點陣夾芯板不同的振動形式,在區(qū)域(1)中方程(53a)有一個零解；在區(qū)域(2)中方程(53a)有兩個零解；在區(qū)域(3)中方程(53a)沒有零解.

圖1 ε11=-1時方程(53a)的平面轉(zhuǎn)遷集Fig. 1 Plan transition variety for Equation (53a) when ε11=-1

圖2 ε11=-1時方程(53a)的三維轉(zhuǎn)遷集Fig.2 Three-dimensional sketch of Transitions variety for Equation (53a) when ε11=-1

圖3表示,當(dāng)ε11=-1,ε12=1和ε13=1時,方程(53a)的第一階模態(tài)的力-幅響應(yīng),其對應(yīng)參數(shù)值分別為:

(a)σ1=-3,σ2=-1,γ=0

(b)σ1=-1,σ2=-3,γ=0

(c)σ1=-3,σ2=-1,γ=1

(d)σ1=-1,σ2=-3,γ=1

圖4表示,當(dāng)ε11=-1,ε12=1和ε13=1時,方程(53a)的第二階模態(tài)的力-幅響應(yīng),其對應(yīng)參數(shù)值分別為:

(a)σ1=-1,σ2=-3,γ=0

(b)σ1=-3,σ2=-1,γ=0

(c)σ1=-1,σ2=-3,γ=1

(d)σ1=-3,σ2=-1,γ=1

圖3 ε11=-1時方程(53a)在不同參數(shù)下第一階模態(tài)的力-幅響應(yīng)曲線Fig.3 Force-amplitude response curves of the first-order mode with different parameters for Equation (53a) when ε11=-1

圖4 ε11=-1時方程(53a)在不同參數(shù)下第一階模態(tài)的力-幅響應(yīng)曲線Fig.4 Force-amplitude response curves of the second-order mode with different parameters for Equation (53a) when ε11=-1

同理,圖5和圖6將方程(53a)平衡點附近鄰域分為不同的區(qū)域,分別對應(yīng)點陣夾芯板不同的振動形式,在區(qū)域(1)中方程(53a)有一個零解；在區(qū)域(2)中方程(53a)有兩個零解；在區(qū)域(3)中方程(53a)沒有零解.

圖7表示,當(dāng)ε11=1,ε12=1和ε13=1方程(53a)的第一階模態(tài)的力-幅響應(yīng),其對應(yīng)參數(shù)值分別為:

(a)σ1=-3,σ2=-1,γ=0

(b)σ1=-1,σ2=-3,γ=0

(c)σ1=-3,σ2=-1,γ=1

(d)σ1=-1,σ2=-3,γ=1

圖5 ε11=1時方程(53a)的平面轉(zhuǎn)遷集Fig.5 Plan transition variety for Equation (53a) when ε11=1

圖6 ε11=1時方程(53a)的三維轉(zhuǎn)遷集Fig.6 Three-dimensional sketch of Transitions variety for Equation (53a) when ε11=1

圖7 ε11=1時方程(53a)在不同參數(shù)下第一階模態(tài)的力-幅響應(yīng)曲線Fig.7 Force-amplitude response curves of the first-order mode with different parameters for Equation (53a) when ε11=1

圖8表示,當(dāng)ε11=1,ε12=1和ε13=1時,方程(53a)的第二階模態(tài)的力-幅響應(yīng),其對應(yīng)參數(shù)值分別為:

(a)σ1=-1,σ2=-3,γ=0

(b)σ1=-3,σ2=-1,γ=0

(c)σ1=-1,σ2=-3,γ=1

(d)σ1=-3,σ2=-1,γ=1

圖8 ε11=1時方程(53a)在不同參數(shù)下第二階模態(tài)的力-幅響應(yīng)曲線Fig.8 Force-amplitude response curves of the second-order mode with different parameters for Equation (53a) when ε11=1

同樣的方法,也可以對方程(53b)、(53c)、 (53d)、 (53e)和(53f)進行數(shù)值計算,針對分叉表達式,滯后表達式和雙極限點集得到關(guān)于三個分叉參數(shù)σ1、σ2和p1的平面轉(zhuǎn)遷集和立體轉(zhuǎn)遷集.

3 結(jié)論

針對1:1內(nèi)共振情況下點陣夾芯板的非線性動力學(xué)分叉方程,推廣了對于含有兩個狀態(tài)變量和三個及三個以上分叉參數(shù)的一般非線性動力學(xué)方程的奇異性理論,得到點陣夾芯板的非線性力學(xué)分叉方程余維4的6個普適開折的表達式,從方程(64)的表達式可以發(fā)現(xiàn):當(dāng)p1固定時,分叉集是一個圓；當(dāng)σ1或σ2固定時,分叉集是一個拋物線,同時,在區(qū)域(1)中,方程(53)有一個零解；在區(qū)域(2)中方程(53)有兩個零解；在區(qū)域(3)中方程(53)沒有零解. 研究結(jié)果對理解點陣夾芯板結(jié)構(gòu)的穩(wěn)定性具有指導(dǎo)意義,并為該類構(gòu)件的工程應(yīng)用提供參考.

1Queheillalt D T, Murty Y,Wadley H N G. Mechanical properties of an extruded pyramidal lattice truss sandwich structure.ScriptaMaterialia, 2008,58(1):76～79

2Lim C H, Jeon I, Kang K J. A new type of sandwich panel with periodic cellular metal cores and its mechanical performances.MaterialsandDesign, 2009,30(8):3082～3039

3Park J S, Joo J H, Lee B C,et al. Mechanical behaviour of tube-woven Kagome truss cores under compression.InternationalJournalofMechanicalSciences, 2011,53(1):65～73

4Kim T,Hodson H P,Lu T J. Fluid-flow and end wall heat-transfer characteristics of an ultralight lattice-frame materia1.InternationalJournalofHeatandMassTransfer, 2004,47(6-7):1129～1140

5Kim T,Hodson H P,Lu T J. Contribution of vortex structures and flows eparation to local and overall pressure and heat transfer characteristics in an ultralight weight lattice material.InternationalJournalofHeatandMassTransfer, 2005,48(19-20):4243～4264

6Lu T J,Valdevit L,Evans A G. Active cooling by metallic sandwich structures with periodic cores.ProgressinMaterialsScience, 2005,50(7):789～815

7Golubitsky M, Langford W F. Classification and unfoldings of degenerate Hopf bifurcations.JournalofDifferentialEquations, 1981,41(3):375～415

8Golubitsky M, Guillemin V. Stable mapping and their singularities. New York:Graduate Texts in Mathematics, 1973

9Martinet J. Singularities of smooth functions and maps. Landon:Cambridge University Press, 1982

10 Golubistky M S, Schaeffer D G. Singularities and groups in bifurcation theory. New York :Spring-Verlag, 1988,220(4):1023～1074

11 Keyfitz B L. Classification of one state variable bifurcation problem up to codimension seven.DynamicsandStabilityofSystems, 1986,1(1):1～41

12 Golubitsky M,Stewart I. Hopf bifurcation with dihedral group symmetry: coupled nonlinear oscillators.MultiparameterBifurcationTheory, 1986:131～173

13 Melbourne I. The classification up to low codimension of bifurcation problems with octahedral symmetry[Ph.D Thesis]. Coventr:University of Warwick, 1988

14 Gaffney T. New methods in the classification theory of bifurcation problem.MultiparameterBifurcationTheory, 1986,56:97～116

15 Furter J E,Sitta A M,Stewart I. Singularity theory and equivariant bifurcation problems with parameter symmetry.MathematicalProcEedingsoftheCambridgePhiloSophicalSocietyR.M.Loynes, 1996,120(3):547～578

16 Gao S P, Li Y C. Classification of (D4,S1)-equivariant bifurcation problems up to topological condition 2.ScienceinChina.Ser, 2003,46(6):862～871

17 郭瑞芝. 等變分歧問題研究[博士學(xué)位論文]. 長沙:中南大學(xué),2006 (Guo R Z. Study on equivariant bifurcation problems[Ph.D Thesis]. Changsha: Central South University, 2006 (in Chinese))

18 崔登蘭,李養(yǎng)成. 含兩組狀態(tài)變量且參數(shù)具有對稱性的等變分歧問題及其開折的穩(wěn)定性. 應(yīng)用數(shù)學(xué)和力學(xué), 2007,28(2):209～215 (Cui D L, Li Y C. Equivariant bifurcation problems and the stability of open fold under contains two sets of state variables and parameters which have symmetry.AppliedMathematicsandMechanics, 2007,28(2):209～215 (in Chinese))

19 Lari-Lavassani A H, Lu Y C. Equivariant multiparameter bifurcation via singularity theory.JournalofDynamicsandDifferentialEquations, 1993,5(2):189～218

20 胡凡努,李養(yǎng)成. 關(guān)于兩狀態(tài)變量組的等變分歧問題的通用開折. 數(shù)學(xué)理論與應(yīng)用, 2000,20(3):50～57 (Hu F N, Li Y C. Versal unfolding of equivariant bifurcation problems about two sets of state variables.MathematicalTheoryandApplications, 2000,20(3):50～57 (in Chinese))

21 高守平,李養(yǎng)成. 多參數(shù)等變分歧問題關(guān)于左右等價的開折. 數(shù)學(xué)年刊A輯, 2003,24(3):341～348 (Gao S P, Li Y C. Open fold of equivariant bifurcation problems with multiparameter under the left and right equivalent group.ChineseAnnalsofMathematics,SeriesA, 2003,24(3):341～348 (in Chinese))

22 郭瑞芝,李養(yǎng)成. 含兩組狀態(tài)變量的等變分歧問題在左右等價群下的開折. 應(yīng)用數(shù)學(xué)和力學(xué), 2005,26(4):489～496 (Guo R Z, Li Y C. Open fold of equivariant bifurcation problems with two sets of state variables under the left and right equivalent group.AppliedMathematicsandMechanics, 2005,26(4):489～496 (in Chinese))

23 Qin Z H, Chen Y S. Singular analysis of bifurcation systems with two parameters.ActaMechanicaSinica, 2010,26(3):501～507

24 Qin Z H, Chen Y S. Singular analysis of a two-dimensional bifurcation system.ScienceChinaTechnologicalSciences, 2010,53(3):608～611

25 Qin Z H,Chen Y S, Li J. Singular analysis of two-dimensional elastic cable with 1:1 internal resonance.AppliedMathematicsandMechanics, 2010,31(2):143～150

26 秦朝紅. 兩狀態(tài)變量、兩分叉參數(shù)系統(tǒng)的分叉分析及其工程應(yīng)用[博士學(xué)位論文]. 哈爾濱:哈爾濱工業(yè)大學(xué), 2010 (Qin Z H. Singularity method for nonlinear dynamical analysis of systems with two parameters and its application in engineering[Ph.D Thesis]. Harbin: Harbin Institute of Technology, 2010(in Chinese))

27 陳建恩. 輕質(zhì)材料層合板的非線性動力學(xué)理論分析與實驗研究[博士學(xué)位論文]. 北京:北京工業(yè)大學(xué),2013 (Chen J E. Theoretical and experimental investigations on nonlinear dynamics of light-weight sandwich plate[Ph.D Thesis]. Beijing: Beijing University of Technology, 2013(in Chinese))