（福建農(nóng)林大學機電工程學院，福建 福州350002）

目前，農(nóng)用大棚在結(jié)構(gòu)上具有很強的區(qū)域性及相似性，在同一區(qū)域內(nèi)，大棚結(jié)構(gòu)系列化程度較高。國內(nèi)外許多學者[1-5]已經(jīng)開始對大棚結(jié)構(gòu)進行研究。對大棚的結(jié)構(gòu)穩(wěn)定性的研究有多種方法，如應(yīng)用材料力學原理、通過有限元算法或有限元軟件進行計算與分析等。然而，由于沒有同時將大棚參數(shù)化的有限元模型結(jié)構(gòu)穩(wěn)定性檢測過程與通用的編程語言結(jié)合起來，使得各類分析法并不適用于大棚產(chǎn)品質(zhì)量檢驗與鑒定或生產(chǎn)與安裝工作人員對大棚產(chǎn)品的質(zhì)量檢測或驗收。為此，本研究以整體的連棟型大棚參數(shù)化模型為基礎(chǔ)，運用有限元分析軟件ANSYS參數(shù)化設(shè)計語言與Visual Basic 6.0結(jié)合，開發(fā)了一種適用性較強的人機交互式大棚結(jié)構(gòu)穩(wěn)定性參數(shù)化仿真檢測系統(tǒng)。

1 大棚有限元模型參數(shù)化

參數(shù)化通常指的是參數(shù)化造型，是一種重要的幾何參數(shù)快速構(gòu)造和修改幾何模型的造型方法[6]。本研究以福建省農(nóng)業(yè)機械鑒定推廣總站推廣的SPG型和SPJ型大棚為研究對象，各系列大棚包括跨度8m和6m，又根據(jù)其使用的管材不同，分為Ⅰ、Ⅱ及Ⅲ級，是具有相同結(jié)構(gòu)，而材料參數(shù)不同的系列化大棚，適合運用參數(shù)化方法建模與分析。

考慮到大棚結(jié)構(gòu)主要以鋼管焊接結(jié)構(gòu)為主，系列化程度高，可對大棚中不影響穩(wěn)定性的附屬結(jié)構(gòu)進行簡化，將大棚的跨度、高度、管材材料截面、材料物理屬性和荷載等作為參數(shù)，并結(jié)合APDL命令流，以建立與調(diào)用ANSYS宏文件的方式，得到參數(shù)化的大棚有限元模型如圖1所示。

模擬大棚的管材采用Beam188單元，薄膜采用Shell63單元，各單元均采用自適應(yīng)方式進行網(wǎng)格劃分。荷載方面，根據(jù)標準[8-9]中所提供的數(shù)據(jù)對其進行標準值的計算，而后施加到大棚模型上。故荷載標準值各個分量，也可以實現(xiàn)參數(shù)化。以SPG-8Ⅲ型大棚模型受雪荷載情況為例，結(jié)構(gòu)參數(shù)取值如表1所示，該型大棚跨度8m，高度3.5m，彈性模量2×1011，泊松比0.3，薄膜厚度0.0002m。

圖1 參數(shù)化的大棚有限元模型

2 檢測系統(tǒng)的開發(fā)

以Visual Basic6.0軟件為開發(fā)平臺，開發(fā)基于Visual Basic與APDL命令結(jié)合的人機交互式農(nóng)用大棚模型結(jié)構(gòu)穩(wěn)定性檢測系統(tǒng)。該系統(tǒng)的技術(shù)路線如圖2所示。

表1 SPG-8Ⅲ級大棚材料參數(shù)

2.1 系統(tǒng)界面

系統(tǒng)主要包括用戶登錄界面、參數(shù)設(shè)置界面、載荷設(shè)置與求解界面以及結(jié)果顯示界面。用戶在登錄后進入?yún)?shù)設(shè)置界面，如圖3所示，可選擇要分析的大棚類型、跨度及型號，再對該種大棚結(jié)構(gòu)的管材截面尺寸進行設(shè)置。然后，進入荷載設(shè)置界面，如圖4所示，設(shè)置所要分析的受載情況，調(diào)用ANSYS軟件進行計算。最后，可進入結(jié)果顯示界面，查看計算分析后的結(jié)果。

圖2 檢測系統(tǒng)的技術(shù)路線

圖3 參數(shù)設(shè)置界面

圖4 荷載設(shè)置與求解界面

2.2 參數(shù)處理及調(diào)用

APDL命令中關(guān)于設(shè)置參數(shù)的代碼，通過Microsoft Access 7.0軟件建立ANSYS認可的APDL代碼參數(shù)表，再通過VB代碼對表中的數(shù)據(jù)進行實時的更新。即在參數(shù)設(shè)置界面下對各參數(shù)對應(yīng)的text控件的 “click”事件[10]利用以下語句實現(xiàn)參數(shù)表中的參數(shù)值實時更新。

SQL＝"UPDATe＃＃SET value＝＇" ＆load.Text ＆"＇WHERE prmt＝＇load.＇"

在參數(shù)設(shè)置完成后，將參數(shù)表以文本的形式與APDL命令代碼中的其余代碼重新編輯成一個新的APDL代碼文件，供VB程序通過代碼直接調(diào)用ANSYS程序進行后臺運行計算。調(diào)用ANSYS的關(guān)鍵代碼如下：

其中，“dy”只是一個用來存放shell函數(shù)返回值的變量。

2.3 后處理結(jié)果

ANSYS軟件可以提供多種的后處理方式，可利用APDL命令得到各種方式的結(jié)果。在VB調(diào)用ANSYS程序后臺運行完畢后，可以應(yīng)用VB進行后處理。如對APDL命令生成位移分布云圖，在VB程序中在新建窗口中以加載方式顯示圖片，以各大系列Ⅲ級型大棚結(jié)果為例，結(jié)果顯示界面如圖5、圖6所示。

圖5 Ⅲ級大棚受雪荷載結(jié)果顯示界面

圖6 Ⅲ級大棚受風荷載結(jié)果顯示界面

3 仿真結(jié)果分析

通過上述方法，可以得到大棚在受荷載后每個單元所發(fā)生的位移量，將得到的最大位移量值與允許位移量值比較，從而分析該大棚結(jié)構(gòu)是否失穩(wěn)，即不失穩(wěn)條件是：

式中，δmax為大棚結(jié)構(gòu)受荷載后得到的最大位移量值；[δ]為大棚結(jié)構(gòu)允許位移量值，為大棚跨度的1／200[11]，即：[δ]SPG8＝ [δ]SPJ8＝8／200＝0.04m；[δ]SPG6＝ [δ]SPJ6＝6／200＝0.03m。

也就是說，跨度為8m的大棚，其允許位移量值為0.04m，而跨度為6m的大棚，其允許位移量值為0.03m，每個大棚的位移量值若超過相應(yīng)的允許位移量值，則將發(fā)生失穩(wěn)現(xiàn)象，甚至倒塌、損壞。表2為各系列大棚在受載后的仿真檢測結(jié)果。

表2 仿真檢測結(jié)果 m

設(shè)最大位移量與允許位移量比值為η，即：

則大棚模型受荷載后情況如表3所示。

表3 各系列大棚受雪荷載η最大值 %

4 結(jié)論

（1）由各位移值分布云圖可知，所建立的整體連棟型大棚模型與單拱整體大棚模型或連棟型非整體大棚模型受荷載后變形的位移分布情況相仿，均在大棚的中部并逐漸向周圍及兩端遞減；受風荷載后的變形在迎風面的1／4拱與中間拱交接處以及背風面附近最大，大棚中間拱較小。而且，整體連棟型大棚模型受荷載后的最大位移量值大于其他類型的大棚模型。原因在于整體連棟型大棚同時考慮了薄膜受載變形、連棟型受荷載情況與骨架間等因素的相互作用而使得大棚骨架的位移增大。

（2）對于大棚模型受雪荷載情況，上述結(jié)果顯示，在對各系列大棚模型進行受雪荷載時的穩(wěn)定性仿真檢測分析時，各大棚均未失穩(wěn)，檢測到的結(jié)果是Ⅰ級的最為穩(wěn)定，Ⅱ級次之，Ⅲ級的穩(wěn)定性較差。原因在于SPG型大棚中，Ⅲ級的大棚骨架管選用的是鍍鋅板材管，Ⅱ級選用的則是熱鍍鋅管，而Ⅰ級選用的是熱鍍鋅管且管的壁厚有所增大；SPJ型大棚中，Ⅲ級的大棚骨架為鍍鋅板材管，Ⅱ級的為增加壁厚的鍍鋅板材管，而Ⅰ級選用與Ⅲ級同樣壁厚的熱鍍鋅管。但是，在造價方面也同樣的Ⅰ級最大，Ⅲ級的最經(jīng)濟。同時，SPJ型抗雪荷載能力明顯高于SPG型，原因在于SPJ型大棚在結(jié)構(gòu)上優(yōu)于SPG型大棚，在同樣的高度下，其拱形弧度較低，又有立柱等輔助支撐，但在造價方面，SPG型較為經(jīng)濟。

（3）對于大棚模型受風荷載情況，上述結(jié)果也顯示，在對各系列大棚模型進行受風荷載時的穩(wěn)定性仿真檢測分析時，各大棚均未失穩(wěn)，SPG、SPJ型大棚各級結(jié)果與其受雪荷載時相類似，均為Ⅰ級最為穩(wěn)定，Ⅱ級次之，Ⅲ級較差。同時，SPJ型同樣在抗風荷載能力上高于SPG型。

（4）上述結(jié)果還說明，SPG型與SPJ型大棚在抗風荷載能力上要比抗雪荷載能力上弱些，在今后的生產(chǎn)與安全防范方面，應(yīng)采取適當?shù)拇胧嗉訌姼鞔笈锏目癸L能力。

（5）ANSYS參數(shù)化設(shè)計語言能夠?qū)崿F(xiàn)整體連棟型大棚建模的參數(shù)化驅(qū)動，可通過ANSYS經(jīng)典操作環(huán)境得到宏代碼，對利用Visual Basic 6.0對其宏代碼進行參數(shù)數(shù)據(jù)編輯、更新，從而得到參數(shù)化的結(jié)構(gòu)穩(wěn)定性檢測系統(tǒng)。

（6）采用Visual Basic與APDL命令結(jié)合，開發(fā)的人機互動式大棚穩(wěn)定性檢測系統(tǒng)，并采用參數(shù)化建立模型，可修改大棚的結(jié)構(gòu)及荷載參數(shù)，對不同結(jié)構(gòu)、受不同荷載的大棚穩(wěn)定性進行分析檢測。該系統(tǒng)操作人性化高、可操作性強，適用于專業(yè)與非專業(yè)大棚產(chǎn)品檢測與鑒定工作人員對大棚產(chǎn)品的結(jié)構(gòu)穩(wěn)定性進行檢驗。

[1]齊 飛，童根樹.連棟溫室鋼結(jié)構(gòu)框架穩(wěn)定設(shè)計方法 [J].農(nóng)業(yè)工程學報，2009，（9）：202-209.

[2]沈正炳，黃文彬.多聯(lián)棟溫室框架結(jié)構(gòu)的彈塑性計算 [J].農(nóng)業(yè)工程學報，2000，（3）：105-108.

[3]蘆永杰，田慧婧.日光溫室的荷載分析及結(jié)構(gòu)優(yōu)化計算 [J].農(nóng)業(yè)技術(shù)與裝備，2010，（4）：15-17.

[4]俞永華，盧如國.塑料大棚結(jié)構(gòu)承載特性的分析研究 [J].農(nóng)機化研究，2011，（1）：165-168.

[5]湯 慶，伍德林，朱世東，等.雙層拱架結(jié)構(gòu)塑料大棚透光率及穩(wěn)定性分析 [J].農(nóng)機化研究，2012，（9）：34-41.

[6]任會禮，李江波，高崇仁.基于ANSYS的塔式起重機臂架有限元參數(shù)化建模與分析 [J].起重運輸機械，2006，（9）：11-13.

[7]龔曙光，黃云清.有限元分析與ANSYS APDL編程及高級應(yīng)用 [M].北京：機械工業(yè)出版社，2009.

[8]中華人民共和國國家標準.建筑結(jié)構(gòu)荷載規(guī)范(06版）GB50009-2001[S].2006-07-25.

[9]中華人民共和國國家標準.溫室結(jié)構(gòu)設(shè)計荷載GB／T 18622-2002[S].2002-01-16.

[10]寧正元.Visual Basic程序設(shè)計教程 [M].北京：清華大學出版社，北京交通大學出版社，2004.

[11]程勤陽.溫室結(jié)構(gòu)設(shè)計的基本方法（一）——溫室結(jié)構(gòu)設(shè)計基本要求及構(gòu)件計算 [J].農(nóng)業(yè)工程技術(shù)-溫室園藝，2006，（9）：11-12.