風(fēng)電機(jī)組螺栓連接參數(shù)化有限元分析
——VBA與APDL聯(lián)合二次開發(fā)

張俊，沈雨虹，廖暉

（東方電氣風(fēng)電有限公司，四川 德陽，618000）

風(fēng)電機(jī)組螺栓連接參數(shù)化有限元分析
——VBA與APDL聯(lián)合二次開發(fā)

張俊，沈雨虹，廖暉

（東方電氣風(fēng)電有限公司，四川 德陽，618000）

在風(fēng)電機(jī)組的整體有限元分析中，重要且困難的部分是對(duì)主要部件之間的連接螺栓進(jìn)行極限強(qiáng)度和疲勞壽命有限元分析。文章介紹風(fēng)電機(jī)組螺栓連接參數(shù)化有限元分析流程，并以輪轂與葉片連接螺栓為例詳細(xì)介紹分析流程中的各個(gè)步驟，最后簡單介紹如何聯(lián)合VBA和APDL進(jìn)行分析流程中所需的二次開發(fā)。

風(fēng)電機(jī)組，螺栓連接，有限元分析，二次開發(fā)，參數(shù)化

0　引言

在風(fēng)力發(fā)電機(jī)組新機(jī)型的開發(fā)和認(rèn)證過程中，為保證機(jī)組的安全運(yùn)行和質(zhì)量，必須對(duì)主要部件（如輪轂、主軸、軸承座、機(jī)架等）及主要部件之間的連接螺栓進(jìn)行有限元分析，而且在風(fēng)電機(jī)組的整體有限元分析中，重要且困難的部分是對(duì)主要部件之間的連接螺栓進(jìn)行極限強(qiáng)度和疲勞壽命有限元分析。其重要性在于螺栓的失效將導(dǎo)致風(fēng)電機(jī)組的損壞，危及到人員安全；其困難在于螺栓的單元類型通常為梁單元 （因用實(shí)體單元模擬螺栓時(shí)，螺紋部位與連接件之間的接觸關(guān)系設(shè)為綁定，傳力面由螺栓桿部的等圓形截面突變?yōu)檫B接件的端部截面，會(huì)在關(guān)鍵螺紋部位產(chǎn)生不真實(shí)的應(yīng)力集中現(xiàn)象，而用梁單元模擬螺栓時(shí)，可將代表螺紋部位的各個(gè)梁單元的節(jié)點(diǎn)通過梁或桿單元連接至螺紋孔面上的各個(gè)節(jié)點(diǎn)，可較真實(shí)和精確地模擬螺栓連接中由多圈螺紋傳力的效果），且整個(gè)螺栓連接有限元模型中單元之間的接觸和連接關(guān)系很多，僅僅通過有限元分析軟件的界面操作無法或很難建立精確的螺栓連接有限元模型，然而通過命令流的形式能很好地解決該問題。

本文介紹了風(fēng)電機(jī)組螺栓連接參數(shù)化有限元分析流程 （包括極限強(qiáng)度和疲勞壽命），并以輪轂與葉片連接螺栓為例詳細(xì)介紹了分析流程中的各個(gè)步驟，最后結(jié)合VBA和APDL兩種語言的特點(diǎn)簡單介紹了如何聯(lián)合該兩種語言進(jìn)行分析流程中所需的二次開發(fā)。應(yīng)用本文介紹的方法，可簡化有限元分析的前后處理工作，并根據(jù)具體的機(jī)型和載荷高效合理地選取螺栓型號(hào)和數(shù)量，有利于機(jī)組優(yōu)化設(shè)計(jì)，具有很高的實(shí)用價(jià)值。

1　參數(shù)化有限元分析流程

圖1　螺栓連接參數(shù)化有限元分析流程圖

風(fēng)電機(jī)組螺栓連接參數(shù)化有限元分析流程如圖1所示。第一步：在Microsoft Office Excel中輸入螺栓連接有限元模型所需的各種參數(shù)，包括螺栓型號(hào)和數(shù)量、連接件尺寸、各個(gè)極限工況下的載荷分量數(shù)據(jù) （用于生成與螺栓極限強(qiáng)度有限元模型對(duì)應(yīng)的宏文件 “Static.mac”）及按不同大小和方向施加的合成彎矩載荷 （用于生成與螺栓疲勞壽命有限元模型對(duì)應(yīng)的宏文件 “Fatigue.mac”）；第二步：將Excel中的部分參數(shù)，包括螺栓型號(hào)和數(shù)量及連接件尺寸傳入ANSYS Workbench的參數(shù)管理器 （WB Parameter Set）中，然后更新靜力結(jié)構(gòu)分析系統(tǒng) （WB Static Structural）的幾何 （Geometry）數(shù)據(jù)單元中的參數(shù)以更新幾何模型；第三步：將Excel中的其余參數(shù)傳入VBA工程，編寫VB過程生成與螺栓極限強(qiáng)度和疲勞壽命有限元模型分別對(duì)應(yīng)的宏文件"Static.mac"和"Fatigue.mac"，打開ANSYS Mechanical，將兩個(gè)文件的內(nèi)容分別復(fù)制到各自對(duì)應(yīng)的Command對(duì)象中；第四步：在ANSYS Mechanical中生成與螺栓極限強(qiáng)度和疲勞壽命有限元模型分別對(duì)應(yīng)的可用于求解的輸入文件"Static.inp"和"Fatigue.inp"，并將其依次導(dǎo)入ANSYS經(jīng)典環(huán)境 （ANSYS Classic）中進(jìn)行求解；第五步：在Excel的VBA工程中編寫VB過程生成與螺栓極限強(qiáng)度和疲勞壽命有限元模型分別對(duì)應(yīng)的用于提取螺栓應(yīng)力的宏文件"StressStatic.mac"和"StressFatigue.mac"，在ANSYS Classic中運(yùn)行這兩個(gè)宏，自動(dòng)提取螺栓的應(yīng)力數(shù)據(jù)并保存于txt文件中。本章后續(xù)內(nèi)容以輪轂和葉片連接螺栓為例詳細(xì)介紹了每步的具體操作過程。

1.1Excel中的參數(shù)輸入

圖2　有限元模型參數(shù) （螺栓型號(hào)和數(shù)量、連接件尺寸等）

輪轂與葉片之間的螺栓連接有限元模型所需的各種參數(shù)如圖2所示，參數(shù)定義表 （Parameter Definitions）包括參數(shù)名稱（Name）、參數(shù)值（Value）、單位 （Unit）及描述 （Description），參數(shù)值列中以紅色表示的數(shù)值為手動(dòng)輸入數(shù)據(jù)，以黑色表示的數(shù)值為自動(dòng)更新數(shù)據(jù) （輸入輪轂側(cè)和葉片側(cè)螺栓螺紋規(guī)格參數(shù)h_TS和b_TS后，單擊Update Data按鈕，即可完成數(shù)據(jù)的自動(dòng)更新）。圖2中未列出的參數(shù) （各個(gè)極限工況下的載荷分量數(shù)據(jù)和按不同大小和方向施加的合成彎矩載荷）在該Excel工作薄的其他工作表列出。

1.2利用WB Parameter Set傳遞和更新參數(shù)

ANSYS Workbench的參數(shù)管理器 （WB Parameter Set）用于管理各個(gè)分析系統(tǒng)（Analysis Systems）中參數(shù)的傳遞和更新，配合組件系統(tǒng)（Component Systems）中的 Microsoft Office Excel組件使用時(shí)，能更加方便、靈活地管理各個(gè)參數(shù)（如圖3所示）。將Excel中的部分參數(shù)，包括螺栓型號(hào)和數(shù)量及連接件尺寸傳入WB Parameter Set中，然后更新靜力結(jié)構(gòu)分析系統(tǒng)（WB Static Structural）的幾何 （Geometry）數(shù)據(jù)單元中的參數(shù)以更新幾何模型，如圖4所示。

圖3　ANSYS Workbench中參數(shù)的傳遞和更新

圖4　輪轂與葉片螺栓連接幾何模型

幾何模型中的軸對(duì)稱部分只建出1/N（N為一圈螺栓的個(gè)數(shù)），且盡可能地切割出規(guī)則的體，以便劃分出規(guī)則的網(wǎng)格，如圖5所示。若建出整個(gè)軸對(duì)稱部分并直接進(jìn)行自由網(wǎng)格劃分，則無法劃分出規(guī)則的網(wǎng)格，有可能導(dǎo)致有限元模型在求解時(shí)不收斂或求解的應(yīng)力結(jié)果不真實(shí)。因主要的分析對(duì)象為螺栓，故輪轂采用粗略的四面體單元?jiǎng)澐志W(wǎng)格，且在靠近螺栓的部分進(jìn)行了網(wǎng)格加密。

圖5　輪轂與葉片螺栓連接網(wǎng)格模型

1.3通過Excel VBA生成宏文件

更新WB Static Structural的Geometry數(shù)據(jù)單元中的參數(shù)后，幾何和網(wǎng)格模型中的軸對(duì)稱部分只有1/N。接下來需要建立整個(gè)軸對(duì)稱部分的網(wǎng)格模型，若還是在Workbench的Mechanical模塊中進(jìn)行界面操作將無法完成后續(xù)工作，因整個(gè)軸對(duì)稱部分的網(wǎng)格模型需通過APDL中的EGEN命令對(duì)單元進(jìn)行旋轉(zhuǎn)陣列來生成，涉及到對(duì)單元和節(jié)點(diǎn)的操作，而Mechanical模塊對(duì)此方面的處理能力較弱。但幸運(yùn)的是Mechanical模塊中提供了使用APDL語言的接口，即可以在Mechanical模塊中插入Command對(duì)象 （見圖6），在該對(duì)象中可直接輸入APDL命令，以輔助完成復(fù)雜有限元模型的建立。

圖6　在Mechanical模塊中插入Command對(duì)象

在Mechanical模塊中插入Command對(duì)象很簡單，難點(diǎn)在于如何完成Command對(duì)象中的APDL代碼 （即"Static.mac"和"Fatigue.mac"兩個(gè)宏文件中的內(nèi)容）。既然所需的所有參數(shù)都列在Excel的工作薄中，便可通過Excel VBA生成宏文件。在Excel VBA中編寫與圖7中4個(gè)命令按鈕對(duì)應(yīng)的事件過程，生成"Static.mac"，"Fatigue.mac"，"Stre ssStatic.mac"，"StressFatigue.mac"4個(gè)宏文件，前兩個(gè)分別為建立完整的輪轂與葉片連接螺栓極限強(qiáng)度和疲勞壽命有限元模型所需的宏文件，后兩個(gè)為后處理時(shí)提取螺栓應(yīng)力所需的宏文件。

圖7　通過Excel VBA生成宏文件

1.4通過ANSYS Mechanical生成輸入文件

在Mechanical模塊中通過菜單項(xiàng)Tools＞W(wǎng)rite Input File...即可生成與極限強(qiáng)度和疲勞壽命有限元模型分別對(duì)應(yīng)的輸入文件"Static.inp"和"Fatigue. inp"。inp文件中包含了完整的有限元模型從前處理（網(wǎng)格劃分和約束、載荷施加）到求解 （求解設(shè)置和載荷步控制）所需的所有數(shù)據(jù)。

1.5提取螺栓應(yīng)力數(shù)據(jù)

將1.4節(jié)中提到的兩個(gè)輸入文件依次導(dǎo)入到ANSYS Classic中進(jìn)行求解，求解完成后通過1.3節(jié)中提到的宏文件"StressStatic.mac"和"StressFatigue.mac"來提取螺栓應(yīng)力，輪轂與葉片連接螺栓在某一極限工況下的應(yīng)力如圖8所示。對(duì)于極限強(qiáng)度有限元模型，提取的螺栓應(yīng)力可以直接用于螺栓極限強(qiáng)度的校核，而對(duì)于疲勞壽命有限元模型，提取的是合成彎矩的大小和方向?qū)β菟☉?yīng)力的二維影響矩陣 （見表1），不能直接用于螺栓疲勞壽命的計(jì)算。計(jì)算螺栓疲勞壽命的步驟如下：

（1）提取合成彎矩的大小和方向?qū)β菟☉?yīng)力的二維影響矩陣；

（2）將合成彎矩大小和方向的時(shí)間序列 （由Bladed后處理生成）在此二維影響矩陣的基礎(chǔ)上進(jìn)行二維插值，得出螺栓的應(yīng)力時(shí)間序列；

（3）結(jié)合SN曲線數(shù)據(jù)統(tǒng)計(jì)螺栓的年損傷值D，從而算出螺栓的壽命為1/D，單位為年。

要注意的是以上算出的壽命只針對(duì)單個(gè)螺栓，應(yīng)對(duì)所有的螺栓進(jìn)行以上3個(gè)步驟的疲勞壽命計(jì)算。

圖8　輪轂與葉片連接螺栓在極限工況下的應(yīng)力云圖

表1　合成彎矩的大小和方向?qū)蝹€(gè)螺栓應(yīng)力的二維影響矩陣

2　VBA與APDL聯(lián)合二次開發(fā)介紹

2.1VBA和APDL語言介紹

VBA（即Visual Basic for Application）是Visual Basic的一種宏語言，主要用來擴(kuò)展Windows的應(yīng)用程序功能，特別是Microsoft Office軟件。Excel5.0版本 （1994年發(fā)行）中就已具備VBA的宏功能。

APDL（即ANSYS參數(shù)化設(shè)計(jì)語言，ANSYS Parametric Design Language）是一種解釋性語言，可用來在ANSYS Classic（即ANSYS經(jīng)典環(huán)境）中完成一些通用性強(qiáng)的任務(wù)，也可用于根據(jù)參數(shù)來建立模型。APDL還包括其他許多特性，如重復(fù)執(zhí)行某條命令、宏、if-then-else分支、do循環(huán)、標(biāo)量、向量及矩陣操作等。

雖然只通過APDL即可完成參數(shù)化有限元建模，但ANSYS Workbench的組件系統(tǒng)中已集成了Microsoft Office Excel，在Excel中定義參數(shù)時(shí)可增加圖示，顯得更加直觀，而且可通過VBA宏實(shí)現(xiàn)根據(jù)某個(gè)輸入條件批量更新參數(shù)值的功能 （如圖2）。因此聯(lián)合VBA和APDL兩種語言進(jìn)行參數(shù)化有限元分析，可達(dá)到取長補(bǔ)短的效果。

2.2VBA與APDL聯(lián)合開發(fā)步驟

VBA與APDL聯(lián)合開發(fā)看似復(fù)雜，實(shí)則簡單，即通過參數(shù)化的Excel VBA宏生成參數(shù)化的APDL宏，參數(shù)化思想貫穿VBA與APDL聯(lián)合開發(fā)的始終。如圖9所示，第一步：通過VBA代碼創(chuàng)建文本文件"Static.mac"并返回用于該文件讀寫的文件流對(duì)象"apdl"；第二步：通過文件流對(duì)象"apdl"在文件"Static.mac"中寫入APDL命令，通過字符串連接運(yùn)算符"&"將VBA中的參數(shù)"Length"的值寫入APDL命令中；第三步：完成所有APDL命令的寫入后關(guān)閉文件流對(duì)象"apdl"。通過以上簡單的三步即可完成VBA與APDL的聯(lián)合，當(dāng)然，這只是萬里長征的第一步，剩下的路程是如何完成所有的APDL命令。

圖9　通過VBA宏生成APDL宏示例

3　總結(jié)

本文介紹了風(fēng)電機(jī)組螺栓連接參數(shù)化有限元分析流程 （包括極限強(qiáng)度和疲勞壽命），并以輪轂與葉片連接螺栓為例詳細(xì)介紹了分析流程中的各個(gè)步驟，最后結(jié)合VBA和APDL兩種語言的特點(diǎn)介紹了如何聯(lián)合兩種語言進(jìn)行分析流程中所需的二次開發(fā)。應(yīng)用本文介紹的方法，可簡化有限元分析的前后處理工作，并根據(jù)具體的機(jī)型和載荷高效合理地選取螺栓型號(hào)和數(shù)量，有利于機(jī)組優(yōu)化設(shè)計(jì)，具有很高的實(shí)用價(jià)值。

［1］Germanischer Lloyd.Guideline for the Certification of Wind Turbines［S］.2010

［2］VDI2230 Part 1，Systematic Calculation of High Duty Bolted Joints，Jointswith One Cylindrical Bolt［S］.

［3］ANSYS，ANSYS Workbench 14.0 Help Documentation［DB］.Workbench User'sGuide

［4］ANSYS，ANSYS Workbench 14.0 Help Documentation［DB］.Mechanical APDL ANSYS Parametric Design Language Guide

［5］Steve Saunders，Jeff Webb.Programming Excel with VBA and.NET［M］.O'Reilly，2006

Parametrized FEA of Bolted Connections of Wind Turbines—Second Development through Combining VBA and APDL

Zhang Jun，Shen Yuhong，Liao Hui

（Dongfang Electric Wind Power Co.，Ltd.，Deyang Sichuan，618000）

In finite element analysis of wind turbines，the important and difficult part is to do ultimate strength and fatigue life FEA of bolted connections between main components.In this article，the analysis flow of parametrized FEA for bolted connections of wind turbines is introduced，and the bolted connection of hub and blade is taken as an example to describe in detail various procedures in the analysis flow.At last，the way of combining VBA and APDL to do necessary second development in the analysis flow is described.

wind turbine，bolted connection，finite element analysis，second development，parameterization

TK83

B

1674-9987（2015）02-0030-05

10.13808/j.cnki.issn1674-9987.2015.02.006

張俊 （1983-），男，工學(xué)碩士，2007年3月畢業(yè)于華中科技大學(xué)機(jī)械工程學(xué)院機(jī)電系，現(xiàn)在東方電氣風(fēng)電有限公司從事結(jié)構(gòu)分析工作。