強聲法清除附著藤壺的力學仿真

易定和，歐陽清，周澤均

(海軍工程大學動力工程學院，湖北 武漢 430033)

0 引言

藤壺又稱為“馬牙”或“蚵沏仔”，屬節肢動物門甲殼綱圍胸目動物，是我國周邊海域一種主要的污損生物。藤壺大多生活在潮間帶，附著棲息于海水中固定或浮動的硬物上，例如船體、浮標、橋墩、碼頭、網箱及網具等［1］。

明確藤壺附著情況及機理，已找到一些清除藤壺的方法，包括:物理清污法、化學防污法和生物防污法。物理清污法主要指使用人工手段或機械設備對附著藤壺的船體表面進行擦除。化學防污法大多基于涂料本身形態或其釋放物質抑制藤壺生長。Schumacher等［2］發現增大涂層微結構的高寬比能顯著降低藤壺幼蟲和綠藻孢子的附著量。文獻 ［3］指出，對涂層引入親水的聚乙二醇 (PEGs)或者同時引入－CH2CH2O－片段與疏水片段能阻止污損生物附著。生物防污法是一種安全環保的防污方法，但當前直接提取海洋生物自身防污物質的效率低下以及制取仿生材料的昂貴成本都使該方法暫時無法大規模使用。

本文對強聲清除附著藤壺進行力學仿真研究，探討強聲清污的理論基礎，明確基于強聲作用使附著藤壺得到清除的2種方式。仿真結果表明，不論是基于強聲產生共振力還是產生“綜合力”都對附著藤壺有清除效果。

1 藤壺附著基本情況

藤壺粘接主要經歷介蟲形幼蟲暫時粘附、介蟲形幼蟲永久性粘接、幼體藤壺粘接和成體藤壺粘接4個過程。具體附著過程為:生活在海洋中的成體藤壺在浮游生物上分離出幼蟲形態的無節幼體，該幼體通過吸取浮游生物體養分最終蛻變成介蟲形幼蟲，幼體初期通過觸角附著接觸表面，該附著稱之為暫時粘接;隨后介蟲形幼蟲發現合適的表面，分泌聚合物并實現永久定居，此時暫時粘接轉變成永久粘接，而介蟲形幼蟲最終變為殼狀幼體，幼體藤壺繼續在基體表面分泌出幼體膠。幼體變為成體后，成體膠也分泌到基材上，使附著進一步牢固。

2 強聲三大基本效應

2.1 空化效應

空化是流體流動過程中局部壓力低于飽和蒸汽壓力時出現的空泡生成、長大、潰滅現象，其實質是在流體動力學及熱力學聯合作用下，液體介質局部液－氣相變［4］。強聲在液體內傳播，導致液體局部形成負壓區，產生空穴或氣泡并最終導致其潰滅?？张轁鐣a生微射流和沖擊波，強大的沖擊波對流體本身及內部物質產生強大的機械攪拌作用。

2.2 聲輻射力效應

在高強聲條件下，聲波的非線性效應顯著，會在聲壓中引起一個“直流”項，這一項時間的平均值具有固定的方向和大小，從而產生聲輻射壓力。聲波產生的聲輻射壓力是驚人的，它可以把常見固體懸浮于空中。長時間的輻射力對物體有“持續作用力”效果。

2.3 聲流效應

氣體或液體媒質中有強聲波傳播時，往往會引起一種非周期的運動，這種現象稱為聲流［5］。聲流可以產生巨大的剪切力，液體中的物體受聲流作用會發生運動或形變。聲流對破壞附面層，加速傳質傳熱，以及清除表面污垢、雜物都是非常有效的［6］。

3 強聲清污仿真分析

對藤壺的清除機理可以理解為對藤壺分泌膠體產生何種作用使其“失效”。當強聲作用在附著藤壺及其膠體表面時，分析強聲三大基本效應可知，其對藤壺分泌膠體可能產生共振力效果和綜合力效果。聲波本身是一種周期性機械波，其導致的空化、聲流效應都會有一定的周期 (頻率)，對附著藤壺施加一定頻率的力會使其產生振動等效果。此外，上述三大效應對流體中物體都會產生“力”效果，該效果不能用某單一變量衡量，暫稱其為“綜合力”效果。下面，對附著藤壺2種可能清除機理作相應研究。

3.1 基于共振力清污仿真分析

考慮藤壺膠體因共振而失效時，可以從膠體固有模態著手。固有模態是事物一種本真自然屬性，當外界作用力頻率與事物固有頻率接近或相同時，事物由于共振具有可能的最大失效破壞。一般說來，每種事物都有若干固有頻率，對應于每種模態，事物的失效形式是不盡相同的。分析藤壺分泌膠體各階模態振型，可以明確藤壺分泌膠體的失效機理。

根據藤壺分泌膠體基本屬性，建立膠體分析模型，設定膠體材料參數。仿真時假定膠體不存在氣孔等空隙，結構連接完好且結合面不存在缺陷。設置其材料模型為軟件內置的viscoelastic material model，查閱相關資料［7－8］，設定其主要參數為密度1190 kg/m3，楊氏模量3E9 Pa，泊松比0.1。然后對邊界條件進行相應約束。在本仿真實驗中，將膠體下邊界設定為固定約束以模擬膠體實際粘接狀態，利用三角形網格單元進行網格劃分，劃分后共得到402個單元。代入求解器進行求解運算，得到“模擬膠體”的前4階固有模態振型如圖1～圖4所示。

圖1 一階模態 (831 Hz)下膠體變形圖Fig.1 A modal of colloid deformation diagram

圖2 二階模態 (833 Hz)下膠體變形圖Fig.2 Two modals of colloid deformation diagram

圖3 三階模態 (899 Hz)下膠體變形圖Fig.3 Three modals of colloid deformation diagram

圖4 四階模態 (931 Hz)下膠體變形圖Fig.4 Four modals of colloid deformation diagram

從圖1～圖4得知，藤壺分泌膠體的前4階固有頻率接近850 Hz。4階模態振型是不盡相同的。在前2階模態中，膠體所受應力基本相當，其變形以偏移變形為主，區別是膠體偏移方向不相同，基于聲場作用使膠體在該頻率下振動，通過對膠體長時間施加“相同方向作用力”可使膠體的附著特性受到破壞，有利于附著藤壺的清除;在第三階模態，膠體的變形基本可以忽略，此時，共振的唯一效果是導致膠體內部不同部位出現不同程度的應力集中，對藤壺的清除效果有限;在第四階模態，該膠體上半部分各處受到近似對稱的應力集中并發生不一致的擠壓變形，該種變形使得藤壺與膠體的連接處出現破壞松動，該種頻率外力作用下可能出現的結果是藤壺被清除，其分泌膠體仍然附著在船體表面。

3.2 基于“綜合作用力”清污仿真分析

采用三維建模方法對單個附著藤壺進行清除力學仿真研究，建立模型如圖5所示。模型主要幾何參數為:下部矩形模擬船體鋼板部分厚度，幾何尺寸為20 mm×15 mm×3 mm;上部矩形模擬對象為靜態流體，幾何尺寸為15 mm×10 mm×7 mm;橢圓體表征單個成體藤壺，其三軸幾何尺寸分別為2 mm，2 mm，2.5 mm，藤壺膠體用三維液體橋［9］形式表征，其主要尺寸為底部半徑1.8 mm，高度1.8 mm，拉伸半角35°。

材料參數設定。對于下部矩形區域，模擬對象為船體，使用結構鋼材料，主要參數為密度7850 kg/m3，楊氏模量3E9 Pa，泊松比0.33;對于藤壺分泌膠體區域，如3.1設置其主要參數為密度1190 kg/m3，楊氏模量3E9 Pa，泊松比0.1;對于藤壺體區域，其外殼具有結構強度大，比較堅硬等特點，自定義材料參數，設定主要參數為密度5000 kg/m3，楊氏模量100E9 Pa，泊松比0.3。

區域及邊界條件設置。考慮藤壺附著時受力情況進行如下屬性設置:設定船體處于靜止狀態，整體設置為固定約束;對于膠體與船體的接觸邊界面，為便于觀察仿真效果，設置其為固定約束。不考慮強聲波在介質中的傳播，設定聲源置于附著藤壺附近 (其幾何位置如圖5中黑點所示)。上部矩形邊界聲場條件設置需要保證邊界條件與實際情況的一致性即不限定聲波傳播范圍。因此，設定該矩形周邊聲場邊界條件為可允許聲波正常傳播。

劃分網格，進行計算求解。對該模型的網格劃分需要考慮2個重要因素。根據文獻 ［10］，在進行聲學方面的仿真實驗時，為了提高計算精度，需要保證每個波長內有6個網格單元的分析要求。此外，為準確描述膠體的物理性質，使用十分詳細的實體建模方法對膠體進行建模，在膠層自由邊緣必須至少劃分4個單元以保證能準確表現出膠層的變形效果，膠層表面內也應劃分盡可能多的網格單元，這樣劃分能正確描繪邊緣變形的衰減［11］。因此，對膠體等按文獻要求設定網格尺寸大小，用三角形單元進行網格劃分，得到如圖6所示的三維模型網格劃分示意圖，該網格模型共有21774個單元。最后代入求解器進行相關求解后得到應力應變示意圖。

圖5 附著藤壺三維模型Fig.5 Barnacle adhesion model

圖6 三維模型網格劃分示意圖Fig.6 3D mesh division

圖7 點聲源 (2.5 W)作用下藤壺及其膠體應力應變圖Fig.7 Under the action of point source(2.5 W)the barnacle stress strain diagram

圖8 點聲源 (5 W)作用下藤壺及其膠體應力應變圖Fig.8 Under the action of point source(5 W)the barnacle stress strain diagram

圖9 點聲源 (7.5 W)作用下藤壺及其膠體應力應變圖Fig.9 Under the action of point source(7.5 W)the barnacle stress strain diagram

圖10 點聲源 (10 W)作用下藤壺及其膠體應力應變圖Fig.10 Under the action of point source(10 W)the barnacle stress strain diagram

在不同功率聲源 (頻率為100 Hz)條件下，藤壺及其分泌膠體受到強聲“綜合力”產生應力應變，如圖7～圖10所示，4種情況的聲功率分別為2.5 W，5 W，7.5 W，10 W。由圖可知:點聲源產生強聲波使附著藤壺受到力作用，導致膠體出現形變。隨著聲源功率的改變，強聲波對附著藤壺的作用效果是有區別的，主要體現在2方面:一是膠體局部所受應力不同，隨聲源功率的增大，藤壺與分泌膠體的連接處出現應力增大、應力進一步集中等現象，當聲功率為10 W時，膠體最大應力達到10.1 MN/m2;聲功率為2.5 W時，膠體最大應力為5.1 MN/m2，二者的區別十分明顯;另一個區別在于膠體形變程度，隨著聲功率的增大，分泌膠體形變越來越顯著 (在仿真分析中不會出現膠體斷裂的現象)，反映藤壺被清除的可能性隨聲功率的增大而增加。

以膠體中心線處 (如圖5)所受應力為研究對象，觀察其應力變化趨勢，得到如圖11所示的數值變化圖。從圖中可以看到:膠體所受應力呈近似“V”型分布，在外力 (強聲波)作用方向受到大的拉應力，在擠壓方向受到大的擠壓應力，在中心區域所受應力相對小;該處膠體所受應力隨聲源功率的增加而增大，且變化趨勢基本保持一致。

圖11 膠體中心線處應力變化圖Fig.11 Colloidal center line of stress change chart

由于將膠體與船體接觸面設置為固定約束致使其在聲波作用時給出的仿真結果表現為所受應力為0和沒有應變，藤壺整體也沒出現清除脫離。但從膠體失效趨勢判斷，可以認為強聲“綜合力”作用在膠體表面使膠體失效部位及方式已有所體現。此外，根據有關海洋生物膠粘劑的綜述文獻 ［12］，科學家合成得到與藤壺分泌膠體性質相近的人工模擬膠，通過實驗得到這種人工模擬膠粘接Fe基質時的拉伸強度為1.52 MPa，剪切強度為1.92 MPa。上述強度與此時強聲作用下膠體 (即將脫離時)所受應力強度基本相當，這也反映強聲對附著藤壺有清除效果。

4 結語

借鑒強聲在除塵、去垢等方面的應用，將強聲引入船舶清污領域并進行仿真研究。強聲具有的三大基本效應是其能夠清污的理論基礎，通過仿真分析可以看到基于強聲產生的共振力或“綜合力”都對附著藤壺有清除效果。下一步準備設計簡易實驗裝置，通過真實實驗驗證強聲的清污效果。

［1］趙梅英，陳立僑，禹娜.藤壺－高度適應附著的海洋污損生物［J］.生物學教學，2004，29(5):56－57.ZHAO Mei-ying，CHEN Li-qiao，YU Na.Barnacle-marine fouling organism of highly adapted to attachment［J］.Biology Teaching，2004，29(5):56 －57.

［2］SCHUMACHER J F，ALDRED N，CALLOW M E，et al.Species-specific engineered antifouling topographies:correlations between the settlement of algal zoospores and barnacle cyprids［J］.Biofouling，2007，23:1.

［3］GUDIPATI C S，FINLAY J A，et al.The antifouling and fouling-release performance of hyperbranched fluoropolymer(HBFP)－poly(ethylene glycol)(PEG)composite coatings evaluated by adsorption of biomacromolecules and the green fouling alga ulva［J］.Langmuir，2005，21(7):3044.

［4］王福霞，趙麗付，顯威，等.空化初生研究綜述［J］.中國石油大學勝利學院，2009，20(3):3 －6.WANG Fu-xia，ZHAO Li-fu，XIAN Wei，et al.Study on cavitation inception［J］.Journal of Shengli College China University of Petroleum，2009，20(3):3 －6.

［5］錢祖文.非線性聲學［M］.北京:科學出版社，1992.QIAN Zu-wen.Nonlinear acoustics［M］.Beijing:Science Press，1992.

［6］馬大猷.現代聲學理論基礎［M］.北京:科學出版社，2004.345 －346.MA Da-you.Modern acoustics theory［M］.Beijing:Science Press，2004.345 －346.

［7］TIAN Y，PESIKA N，ZENG H，et al.Adhesion and friction in gecko toe attachment and detachment.Proceedings of the National Academy of Sciences，2006，103(51):19320 －19325.

［8］GAO H J，WANG X，YAO H M，et al.Mechanics of hierarchical adhesion structures of geckos［J］.Mechanics of Materials，2005，37(2 －3)275 －285.

［9］周群，于鵬，何斌，等.螞蟻附著力的測試及ANSYS分析［J］.同濟大學學報(自然科學版)，2008，36(5):670 －673.ZHOU Qun，YU Peng，HE Bin，et al.Testing of wet adhesive forces of ants and ANSYS analysis［J］.Journal of Tongji University(Natural science)，2008，36(5):670 －673.

［10］馬飛天，林逸彭，彥紅，等.轎車車內低頻噪聲的仿真計算及實驗研究［J］.中國機械工程，2005，16(16):1489 －1492.MA Fei-tian，LIN Yi-peng，YAN hong，et al.Simulation and experimental study of car interior low frequency noise［J］.China Mechanical Engineering，2005，16(16):1489 －1492.

［11］關英俊.膠粘結固緊調焦鏡有限元分析［J］.光電工程，2010，37(10):1 －5.GUAN Ying-jun.Finite element analysis of the focusing mirror fastened by adhesion［J］.Opto-Electronic Engineering，2010，37(10):1 －5.

［12］賀紅彬，王曉光，宋陽，等.天然海洋生物膠粘劑的研究進展［J］.中國膠粘劑，2006，15(5):42－45.HE Hong-bin，WANG Xiao-guang，SONG Yang，et al.Research progress of marine bioadhesives［J］.China Adhesives，2006，15(5):42 －45.