基于PSO 算法的有限元仿真模型

張曉娜 常樂冉 吳煒 王一晴 山東科技大學(xué)

引言

標(biāo)槍運動是我國一貫較為重視的一項運動，標(biāo)槍競技的勝負(fù)裁決是根據(jù)投擲標(biāo)槍的距離所給出的。本文主要采用粒子群算法（PSO）進行標(biāo)槍投擲最佳距離的求解，通過迭代尋求全局最優(yōu)解，較為精確地獲得標(biāo)槍投擲的最佳距離以及相關(guān)角度、速度的關(guān)系。在了解標(biāo)槍幾何形狀和尺寸參數(shù)的基礎(chǔ)上，構(gòu)建標(biāo)槍的有限元模型，同時借助計算機仿真獲得最優(yōu)投擲距離和角度。

1 計算標(biāo)槍投擲距離

在無風(fēng)的前提下，排除標(biāo)槍飛行中所受空氣影響，從運動學(xué)角度對標(biāo)槍飛行軌跡進行分析。假設(shè)圖中投擲者出手高度為ho，初速度為vo，投槍角為透出器械，標(biāo)槍沿一條曲線運動，該曲線經(jīng)過與A 點同高的B 點并落于地面上的C 點。

圖1.標(biāo)槍投擲飛行軌跡

根據(jù)運動學(xué)可知，在不考慮空氣阻力的情況下，這條曲線近似為拋物線。根據(jù)拋物線軌跡來分析標(biāo)槍的運動，把運動軌跡分為AB 段和CD 段，即S1 和S2。標(biāo)槍由A 點投出，經(jīng)過最高點而落到與A 點同高的B 點。根據(jù)運動學(xué)的斜拋運動可知，水平投擲距離為：

標(biāo)槍由B 運動到D 的水平距離S2，可以近似的把曲線BC 看成是直線，其斜角近似等于出手角度，因此得：

所以，總的投擲距離為：

根據(jù)題目中給出的投擲出手速度為29.70m/s、出手角為36.6°、初始攻角為-0.9°，可以計算出投槍角a=35.7°，將以上數(shù)據(jù)帶入3-3式，即可得出標(biāo)槍投擲距離約為70.24m。

2 PSO 模型建立

通過粒子群算法，根據(jù)鳥類尋食的仿生學(xué)，結(jié)合個體的最優(yōu)解和群體的最優(yōu)解不斷地迭代來尋找區(qū)域的最優(yōu)解。位置限制[28,44]，取初始種群N 為50，迭代次數(shù)為100,在28~44 內(nèi)隨機生成一個50x1 的數(shù)據(jù)矩陣，每次搜尋都需要將當(dāng)前的適應(yīng)度和最優(yōu)解同歷史的記錄值進行對比，如果超過歷史最優(yōu)值，則更新個體和種群的歷史最優(yōu)位置和最優(yōu)解。

圖2.投擲距離變化圖

圖3.標(biāo)槍投擲最佳距離

3 有限元模型的建立

首先根據(jù)已知的標(biāo)槍尺寸確定標(biāo)槍的幾何形狀，標(biāo)槍的結(jié)構(gòu)分為槍頭、槍身及把手。根據(jù)其物理模型對每部分的尺寸進行設(shè)置，并確定標(biāo)槍質(zhì)量、形心位置等，利用有限元法對幾何模型進行網(wǎng)格式劃分，創(chuàng)建標(biāo)槍的有限元模型。

圖4.標(biāo)槍的有限元模型和空氣包模型

在建立了標(biāo)槍的有限元模型以后,在其周圍建立一個 2640 mm× 37mm×2mm 的長方體幾何模型；其次,對該空氣氣流場定義單元類型 (SolidElement)和材料屬性(Null)，假設(shè)標(biāo)槍為剛體，其材質(zhì)設(shè)定為合金材料；再次,對其進行有限元網(wǎng)格的劃分。同樣,所有的單元也為四邊形單元,共劃分了195360 個單元,模型如圖所示。

計算機仿真

將標(biāo)槍的有限元模型和空氣包模型建立完成以后,最后對這一系統(tǒng)模型進行第二次的參數(shù)設(shè)置。在參數(shù)設(shè)置的過程中,將對空氣采用Euler 單元，對標(biāo)槍采用Lagrange 單元,它們之間的耦合采用*CONSTRAINED_LAGRANGE_IN_SOLID 命令來實現(xiàn)。

在對標(biāo)槍一空氣包建立有限元模型以及對其賦予材料、單元屬性、邊界條件、初始條件以及風(fēng)對標(biāo)槍的作用力以后,使用ANSYS/LS-DYNA 求解器對其進行求解。

在本文的研究中,對標(biāo)槍飛行的計算機模擬均是在考慮空氣阻力的情況下進行的。通過給定空氣的密度、粘性系數(shù)和狀態(tài)方程,利用流一固耦合的方法來模擬標(biāo)槍在空氣內(nèi)的運動情況。

由以下模擬結(jié)果可以看出,風(fēng)對標(biāo)槍的飛行遠(yuǎn)度的影響比較顯著。在其它初始條件`定的情況下,隨著逆風(fēng)風(fēng)速的增大,標(biāo)槍的飛行遠(yuǎn)度呈現(xiàn)出減小的趨勢。但是,在同樣的逆風(fēng)風(fēng)速下,當(dāng)初始條件中的出手角度不斷增大,標(biāo)槍的飛行遠(yuǎn)度又不斷增加,當(dāng)出手角度增大到42 度時,它的飛行遠(yuǎn)度又有所減小。由此也可說明,在空氣阻力的作用下,標(biāo)槍的最佳出手角度并不是45 度。

標(biāo)槍投擲最優(yōu)解

表1

4 標(biāo)槍投擲規(guī)律的總結(jié)

圖5. 標(biāo)槍投擲影響因素

圖6.標(biāo)槍初速度對投擲距離的影響

分析影響標(biāo)槍投擲距離的相關(guān)因素，主要從運動員、標(biāo)槍參數(shù)、比賽環(huán)境進行考慮，基于標(biāo)槍投擲距離模型的構(gòu)建，得出影響標(biāo)槍飛行距離的因素具體有以下幾點：

(1）在其他初始條件不變的情況下，標(biāo)槍飛行距離隨初速度vo增加而增加。圖中標(biāo)槍的出手角度為44°，初始俯仰角速度，初始攻角均為零，可見標(biāo)槍的初速度增加1m/s,標(biāo)槍的飛行距離增加4-5m，當(dāng)初速度繼續(xù)增加時，飛行距離有繼續(xù)增大的趨勢。

(2）在其他初始條件不變的情況下，標(biāo)槍在出手角度為0°~44°的范圍內(nèi)，隨著出手角度的增加，標(biāo)槍的飛行距離基本上是不斷增加的。但是，從模擬結(jié)果上來看，當(dāng)出手角度大于44°時，在空氣阻力的影響下，標(biāo)槍的飛行距離會由最大開始下降，所以標(biāo)槍的最佳出手角度大約在44°左右。

表2

表2.標(biāo)槍的初始攻角與初始俯仰角組合對標(biāo)槍飛行距離的影響

(3）在其他初始條件不變的情況下，且出手角度一定時，初始攻角和初始俯仰角速度的組合會影響標(biāo)槍飛行距離，其中最佳組合為“，當(dāng)初始條件偏離這一組合式，標(biāo)槍飛行距離有所減小。

(4）在其他初始條件一定的情況下，初始攻角的變化對其飛行距離的影響并不顯著，但是，在問題三種所討論的初始攻角與出手角度的不同組合中，其中最佳組合為“”時，標(biāo)槍的飛行距離最遠(yuǎn)，達到91.72m，當(dāng)初始攻角和出手角度的組合偏離這一組合時，飛行距離有所減小。

(5）風(fēng)對標(biāo)槍的飛行距離的影響較為顯著。根據(jù)第四問的計算結(jié)果進行分析，對于逆風(fēng)而言，在其他初始條件一定的情況下，隨著逆風(fēng)風(fēng)速的增大，標(biāo)槍的飛行距離呈現(xiàn)出減小的趨勢。但是，在同樣的逆風(fēng)風(fēng)速下，當(dāng)初始條件中的出手角度不斷增大，標(biāo)槍的飛行距離有不斷增大，當(dāng)出手角度AO增大到44°時，它的飛行遠(yuǎn)度有所減小。

5 結(jié)論