，

(1.大連理工大學 運載工程與力學學部 遼寧 大連 116024； 2.廣西大學 土木建筑工程學院 廣西 南寧 530004；3. 廣西大學 工程防災與結構安全教育部重點實驗室 廣西 南寧 530004)

0 引言

加肋板由于具有輕質高強的優(yōu)點而被廣泛應用在橋梁、海洋、船舶、高鐵和航空航天等領域中.在最初的研究中，主要利用梁格模型和正交各向異性板模型來模擬加肋板，而后發(fā)現(xiàn)這些模型僅適用于肋條布置較密的加肋板.之后的研究更傾向于將加肋板考慮為一個混合結構，肋條被看作是梁[1].在對加肋板的彎曲分析中，文獻[2]利用能量的標準對加肋板在受均布荷載作用下的穩(wěn)定性進行了分析；文獻[3]對有縱向肋和橫向肋的情況下加肋板受幾種不同的彎曲荷載進行了研究；文獻[4]利用一階剪切變形理論和移動最小二乘近似求解了矩形加肋板的線性彎曲問題.無網格法因利用一系列節(jié)點而非網格離散結構，相較有限元法更多了一些靈活性.在對加肋板進行計算時，肋條上的節(jié)點分布并不需要與板上的節(jié)點對齊，肋條位置變化也不需要對板節(jié)點重新分布.文獻[5]通過迭代優(yōu)化算法得到了實心變厚度薄板優(yōu)化問題的最優(yōu)解是一個具有無限細、無限密的加肋結構.本文基于無網格法，采用移動最小二乘近似和一階剪切變形理論得到加肋板的位移場函數(shù)，構造加肋板應變能，獲得計算板撓度的控制方程，再利用約束隨機方向法對加肋板的肋條位置進行優(yōu)化，使加肋板中心點的撓度最小.最后用算例驗證了本文方法的有效性，并對不同荷載作用下的加肋板肋條布置進行了優(yōu)化分析.

1 加肋板的無網格模型

加肋板的無網格模型如圖1所示.可以看出，x向肋條上的節(jié)點和y向肋條上的節(jié)點布置并沒有沿著板上節(jié)點的位置布置，而是分布在板上兩行或者兩列節(jié)點之間.根據移動最小二乘近似[6-7]，可得形函數(shù)的表達式為

圖1 加肋板的無網格模型Fig.1 The meshfree model of stiffened plate

(1)

式中：P(x)代表基函數(shù)矩陣；A(x)、B(x)分別是由權函數(shù)和基函數(shù)組成的矩陣.由移動最小二乘近似和一階剪切變形理論可以得到加肋板的位移場函數(shù)為

(2)

式中：u0pI、φpxI、v0pI、φpyI、wpI是板的節(jié)點參數(shù)；n代表板節(jié)點的個數(shù).

x向肋條的位移場函數(shù)可以表示為

(3)

式中：u0sI、φsxI、wsxI是x向肋條上的節(jié)點參數(shù)；N是肋條上的節(jié)點總數(shù)；φxI是肋條上節(jié)點由移動最小二乘近似計算出來的形函數(shù).

同樣地，y向肋條也有相似的表達式.

1.1 位移協(xié)調條件

圖2 加肋板的位移協(xié)調條件Fig.2 The displacement compatibility of stiffened plate

加肋板的位移協(xié)調條件如圖2所示.雖然加肋板在計算時將肋條考慮為梁，但是加肋板的肋條和板是結合在一起的一個整體，所以在肋條上的某一個節(jié)點必有與之相對應的板上的點，使得它們的位移有如下關系：

[wp]p=[wsx]s,

(4)

[φpx]p=[φsx]s,

(5)

[up]c=[usx]c,

(6)

式中：p點是與肋條上節(jié)點s相對應的板上的點；c點是肋條上節(jié)點s與板上的點p之間的連線與接觸面的交點.

將式(2)、 (3)帶入式(4)～(6)中，可得到用矩陣形式表示的位移協(xié)調關系為

Tpδpw=Tsxδsxw,

(7)

Tpδpφ=Tsxδsxφ,

(8)

Tpδpu+esTpδpxφ=Tsxδsu.

(9)

δsx=Tspxδp，

(10)

其中：δsx=[u0s1,wsx1,φsx1,…,u0sN,wsxN,φsxN]T，δp=[u0p1,v0p1,wp1,φpx1,φpy1,…,u0pn,v0pn,wpn,φpxn,φpyn]是一個3N×3n的矩陣，當加肋板上肋條位置改變時只需要重新計算T即可.

1.2 控制方程

根據一階剪切變形理論，加肋板的勢能可以表示為

.

(11)

x向肋條的勢能可以表示為

.

(12)

則加肋板的總勢能可以表示為U=Up+Usx+Usy，根據虛功原理可得加肋板的控制方程為Kδp=f.

1.3 完全轉換法

由于形函數(shù)并不滿足N1J=δ1J的條件，所以并不能直接施加邊界條件,而是需要通過一些方法使得剛度方程中的位移表現(xiàn)出的是真實的位移，以便施加邊界條件.通過完全轉換法[8]之后的節(jié)點形函數(shù)滿足N1J=δ1J的條件，可以像有限元法那樣直接施加邊界條件.

2 約束隨機方向法在肋條布置優(yōu)化中的程序流程

約束隨機方向法是在滿足約束條件的可行范圍內隨機產生方向進行搜索的一種方法.這種方法主要依靠的是隨機生成的點與向量，并利用隨機生成的點與向量找到近似于目標函數(shù)的最速變化方向進行求導操作.因為肋條位置優(yōu)化的數(shù)學模型是選取剛度方程作為目標函數(shù)，而對剛度方程進行求導是非常耗費時間的一項操作，所以選取約束隨機方向法進行近似.

利用無網格法和約束隨機方向法對加肋板進行優(yōu)化計算,具體的程序流程如圖3所示.

圖3 基于無網格法和約束隨機方向法的程序流程Fig.3 Program flow based on meshfree method and constrained random direction method

3 算例分析

3.1 驗證分析

圖4 單根肋條加肋板受均布荷載作用(單位：m)Fig.4 Single rib stiffened plate under uniform loads (unit: m)

序號肋條位置y/m板中心點撓度/10-4m優(yōu)化結果的相對誤差/%10.512 9091.499 051.7420.489 8981.489 441.0930.498 8101.474 450.0740.442 5721.967 2233.5250.608 4482.906 1897.2560.493 4711.480 600.4970.498 0201.474 820.1080.446 2031.906 9729.4390.248 8178.741 17493.27100.495 2731.477 570.28

對約束隨機方向法進行了驗證分析.當加肋板受到均布荷載時，單根肋條在板中間時加肋板的撓度最小.取邊長L=1 m，板厚h=0.01 m的正方形加肋板,肋條高度hs=0.1 m，肋條寬度ws=0.01 m,彈性模量E=1.7×107Pa，泊松比為0.3，受z向均布荷載1 Pa的作用，單根肋條加肋板受均布荷載作用如圖4所示.

無網格模型采用9*9的均布節(jié)點離散，四邊固支邊界條件.采用二次基函數(shù)，權函數(shù)的支撐域大小取4倍的相鄰兩節(jié)點之間的距離.當兩次相鄰計算結果的差值小于0.01 m時，迭代終止.當肋條位于板中間時，加肋板中心點的撓度為1.473 38×10-4m.

為了驗證約束隨機方向法的尋優(yōu)效果，對加肋板受到均布荷載的情況進行了優(yōu)化分析.表1為單根肋條受到均布荷載時中心點撓度的優(yōu)化結果.

由表1可以看出，有些計算結果的相對誤差很大，這是因為約束隨機方向法的局部尋優(yōu)能力較強,隨機選取的初始點對于尋優(yōu)的整個過程以及尋優(yōu)的效率起到非常關鍵的作用.隨機初始點的值對于約束隨機方向法來說是很重要的，往往要進行很多次計算來確定一個適用的最優(yōu)設計方案.

選取結果較優(yōu)的一組數(shù)據和較差的一組數(shù)據對尋優(yōu)過程進行模擬，其中第7組和第9組的數(shù)據尋優(yōu)過程分別如圖5、圖6所示.圖中虛線為y=0.5 m的標準線，當初始點的取值很接近最優(yōu)點時，得出的優(yōu)化結果有極大概率為近似最優(yōu)解.因此，對于約束隨機方向法的優(yōu)化需要進行多次計算.

另外，通過圖形的中心點撓度曲線可以看出總趨勢是趨于下降的，這就證明了約束隨機方向法的迭代過程是具有優(yōu)化意義的，最終得到的點一定會比初始點的值更優(yōu)秀.但是，約束隨機方向法僅僅可以在局部進行尋優(yōu)，全局尋優(yōu)能力較弱，容易找到局部最優(yōu)點.

圖5 第7組數(shù)據尋優(yōu)過程Fig.5 The seventh set of data optimization process

圖6 第9組數(shù)據尋優(yōu)過程Fig.6 The ninth set of data optimization process

圖7 單根肋條加肋板受矩形局部荷載作用(單位：m)Fig.7 Single rib stiffened plate under rectangular partial loads (unit: m)

序號肋條位置y/m板中心點撓度/10-5m10.425 1391.509 9320.421 5411.524 5530.379 7652.451 3340.380 8942.408 5850.530 4805.815 5660.443 6591.609 5570.432 9791.515 9480.426 8171.506 8290.422 7311.518 52100.358 7373.417 39

3.2 加肋板肋條位置優(yōu)化設計

3.2.1局部荷載下肋條的最佳位置 利用約束隨機方向法對單根肋條加肋板受到局部荷載的狀態(tài)進行優(yōu)化設計.加肋板受到的矩形均布荷載仍為1 Pa，方向為z向，分布在左下角的正方形區(qū)域中，其他條件不變.單根肋條加肋板受矩形局部荷載作用如圖7所示.

同樣對受到局部荷載的加肋板進行了優(yōu)化分析，分別取其中10次計算數(shù)據進行分析對比，結果如表2所示.

由表2可以看出，當肋條位置在0.43 m附近時，板中心點的撓度是最小的，最小值為1.506 82×10-5m，是由第8組數(shù)據得出的.第1組、第2組、第7組和第9組數(shù)據都很接近第8組的計算結果，也從側面印證了第8組計算結果的準確性.由此可以得出，當單根肋條加肋板受到左下角局部荷載作用時，肋條的最佳擺放位置為0.43 m，此時加肋板中心點的撓度值為1.506 82×10-5m.

3.2.2雙平行肋條加肋板位置優(yōu)化 當平板中有兩根肋條時，約束隨機方向法的設計變量就變成了2個，且兩根肋條相互平行，間距為0.3 m，受到正方形局部荷載的作用.雙平行肋條加肋板受局部荷載作用如圖8所示,得到的10組結果見表3.

由表3可以看出，有兩根平行肋條的加肋板當固定肋條之間的距離為0.3 m，優(yōu)化得到的結果是兩根肋條的位置分別位于0.4 m和0.7 m時，加肋板中心點的撓度最小.

3.2.3雙垂直肋條加肋板位置優(yōu)化 當有兩根相互垂直的肋條共同作用在板上時，同樣只是對固支邊界條件的加肋板進行了分析，雙垂直肋條加肋板受局部荷載作用如圖9所示.

圖8 雙平行肋條加肋板受局部荷載作用(單位：m)Fig.8 Double parallel ribs stiffened plate under partial loads (unit: m)

圖9 雙垂直肋條加肋板受局部荷載作用(單位：m)Fig.9 Double vertical ribs stiffened plate under partial loads (unit: m)

加肋板為1 m×1 m×0.01 m的正方形薄板，兩根肋條的尺寸為1 m×0.01 m×0.1 m.受力占整個板面的四分之一，位置在加肋板的左下角.計算得到的10組結果如表4所示.

由表4可以看出，當x向肋條位于0.78 m左右、y向肋條位于0.43 m左右時，得到撓度最小的第5組數(shù)據，第9組數(shù)據也證明了在0.7 m和0.4 m附近擁有近似精確解的可能性大.所以，當雙垂直肋條加肋板受到矩形局部荷載作用時，兩根肋條的擺放位置在0.43 m和0.78 m處為佳.

表3 雙平行肋條加肋板優(yōu)化結果

表4 雙垂直肋條加肋板優(yōu)化結果

4 結論

用約束隨機方向法研究了加肋板肋條位置優(yōu)化的問題.對3個算例進行計算，得出了肋條擺放位置的近似最優(yōu)解.利用無網格法進行優(yōu)化設計，由于約束隨機方向法的局部性需要進行多次計算，但是增加的計算時間多是優(yōu)化方法的計算時間；而利用有限元法進行優(yōu)化計算，在重復計算的過程中增加的不僅僅是優(yōu)化的時間，還有因網格重構而增加的有限元計算時間.因此，利用無網格法進行計算可以提高計算效率，效率提高程度依賴于初始點的選取.約束隨機方向法計算方法簡單，不需要太大的計算量，是進行小規(guī)模優(yōu)化的有效計算方法.

算例中加肋板被當作板與梁的結合體，板與梁均為考慮剪切變形影響的Midlin板模型與Timoshenko梁模型.通過算例計算可以得到，在以固支邊界為前提的條件下，雙平行肋條加肋板在固定肋間距離為0.3 m時，兩根肋條的擺放位置分別在0.4 m和0.7 m處；雙垂直肋條中x向肋條的擺放位置在0.78 m、y向肋條的擺放位置在0.43 m時，加肋板中心點的撓度值最小.此方法計算簡單且容易操作，但是比較容易陷入局部死點，需要通過大量的重復計算來確定優(yōu)化的最終方案.