法蘭蓋螺栓襯墊

1、ANASYS有限元分析題目：螺栓預(yù)緊法蘭盤(pán)襯墊有限元分析姓名：老師：學(xué)號(hào)：學(xué)院：機(jī)電學(xué)院專(zhuān)業(yè)：機(jī)械工程螺栓預(yù)緊法蘭盤(pán)襯墊有限元分析西南石油大學(xué) ，成都 610500摘要：法蘭是使管子與管子或閥門(mén)相互連接的零件，連接于管端。法蘭上有孔眼，螺栓使兩法蘭緊連，法蘭間用襯墊密封。法蘭聯(lián)接的螺栓預(yù)緊力是保證法蘭聯(lián)接點(diǎn)不發(fā)生泄漏的重要影響因素。給定預(yù)緊力，通過(guò)anasys 分析得出螺栓最適用于法蘭螺栓連接系統(tǒng)的預(yù)緊力。關(guān)鍵詞：法蘭;螺栓預(yù)緊;襯墊;1. 引言1.1. 法蘭連接概念法連接就是把兩個(gè)管道、管件或器材，先各自固定在一個(gè)法蘭盤(pán)上，兩個(gè)法蘭盤(pán)之間，加上法蘭墊，用螺栓緊固在一起，完成了連接。有的管件和

2、器材已經(jīng)自帶法蘭盤(pán)，也是屬于法蘭連接【1】。法蘭蓋，就是蓋的邊緣有法蘭緣，可以用螺栓固定的法蘭連接是管道施工的重要連接方式。 1.2. 螺栓預(yù)緊法蘭蓋螺栓法蘭聯(lián)接是壓力容器，石油化工設(shè)備及管道應(yīng)用極為廣泛的一種可拆式靜密封聯(lián)接結(jié)構(gòu)。法蘭螺栓連接系統(tǒng)主要失效形式是泄露，而螺栓預(yù)緊和加裝襯墊是保證連接面不發(fā)生泄露的重要環(huán)節(jié)之一。在螺栓連接中，螺栓在安裝的時(shí)候必須擰緊，即在連接承受工作載荷之前，預(yù)受到力的作用，這個(gè)預(yù)加的力稱(chēng)為預(yù)緊力，預(yù)緊的目的在于增強(qiáng)連接的可靠性和緊密性，以防止受載后被連接件出現(xiàn)縫隙或發(fā)生相對(duì)滑移。螺栓預(yù)緊力不能過(guò)大或者過(guò)小。預(yù)緊力過(guò)小，達(dá)不到墊片壓緊并實(shí)現(xiàn)初始密封條件；預(yù)緊力過(guò)大

3、，使墊片壓壞或擠出。所以知道預(yù)緊力的大小值就變得尤其重要【2】。本文通過(guò)給定不同預(yù)緊力，通過(guò)anasys有限元分析，得出螺栓，襯墊，法蘭的應(yīng)力圖。2. 理論分析及其有限元分析2.1. 有限元分析的目的和意義法蘭為圓形盤(pán)狀，設(shè)計(jì)螺栓組布置時(shí)要滿(mǎn)足：在布置螺栓時(shí)要求螺栓的列應(yīng)有合適的間距，邊距；螺栓的布置應(yīng)使各螺栓的受力合理；分布在同一圓周上螺栓數(shù)目，應(yīng)取4,6,8等偶數(shù)；避免螺栓承受附加的彎曲載荷。螺栓工作時(shí)承受工作載荷和預(yù)緊力的聯(lián)合作用，這種螺栓連接在緊螺栓連接中比較常見(jiàn)，因此也是最重要的一種。襯墊是介于法蘭蓋與法蘭蓋之間的填封介質(zhì)。襯墊的材料以及預(yù)緊力不同，會(huì)對(duì)法蘭蓋的工作性能產(chǎn)生極大的影響

4、。螺栓預(yù)緊法蘭蓋襯墊分析，選4個(gè)螺栓對(duì)稱(chēng)布置，1個(gè)襯墊，2兩法蘭蓋組成。螺栓和法蘭蓋由186非線(xiàn)性單元模擬，襯墊用194非線(xiàn)性單元模擬。約束邊界為約束法蘭內(nèi)圈的徑向和周向位移。與施加在4個(gè)螺栓上的力都為1000N。在anasys進(jìn)行分析室，在法蘭內(nèi)徑的徑向施加約束。2.2. 模型理論分析有限單元法的基本思想是將連續(xù)的求解區(qū)域離散為一組有限個(gè)、且按一定方式相互聯(lián)結(jié)在一起的單元的組合體。由于單元能按不同的聯(lián)結(jié)方式進(jìn)行組合，且單元本身又可以有不同形狀，因此可以對(duì)復(fù)雜的模型進(jìn)行求解。有限單元法作為數(shù)值分析方法的另一個(gè)重要特點(diǎn)是利用在每一個(gè)單元內(nèi)假設(shè)的近似函數(shù)來(lái)分片地表示全求解域上待求的未知場(chǎng)函數(shù)。單元

5、內(nèi)的近似函數(shù)通常由未知場(chǎng)函數(shù)或其導(dǎo)數(shù)在單元的各個(gè)節(jié)點(diǎn)的數(shù)值和其插值函數(shù)來(lái)表達(dá)。這樣一來(lái)一個(gè)問(wèn)題的有限元分析中，未知場(chǎng)函數(shù)或及其導(dǎo)數(shù)在各個(gè)節(jié)點(diǎn)上的數(shù)值就成為新的未知量，從而使一個(gè)連續(xù)的無(wú)限自由度間題變成離散的有限自由度問(wèn)題。一經(jīng)求解出這些未知量，就可以通過(guò)插值函數(shù)計(jì)算出各個(gè)單元內(nèi)場(chǎng)函數(shù)的近似值，從而得到整個(gè)求解域上的近似解。顯然隨著單元數(shù)目的增加，也即單元尺寸的縮小，或者隨著單元自由度的增加及插值函數(shù)精度的提高，解的近似程度將不斷改進(jìn)。如果單元是滿(mǎn)足收斂要求的，近似解最后將收斂于精確解。目前，有限單元法的應(yīng)用已由彈性力學(xué)平面間題擴(kuò)展到空間問(wèn)題、板殼問(wèn)題，由靜力平衡問(wèn)題擴(kuò)展到穩(wěn)定間題、動(dòng)力問(wèn)題和波

6、動(dòng)問(wèn)題。分析的對(duì)象從彈性材料擴(kuò)展到塑性、粘彈性、粘塑性和復(fù)合材料等，從固體力學(xué)擴(kuò)展到流體力學(xué)、傳熱學(xué)等連續(xù)介質(zhì)力學(xué)領(lǐng)域.在工程分析中的作用已從分析和校核擴(kuò)展到優(yōu)化設(shè)計(jì)并和計(jì)算機(jī)輔助設(shè)計(jì)技術(shù)相結(jié)合。2.3. Anasys分析步驟Anasys分析主要步驟：定義參數(shù)；選擇單元；定義材料屬性；定義材料壓縮曲線(xiàn)；定義襯墊的線(xiàn)性卸載曲線(xiàn)；創(chuàng)建幾何體；劃分網(wǎng)格；定義預(yù)拉伸單元；定義接觸單元節(jié)點(diǎn)組；定義目標(biāo)單元節(jié)點(diǎn)組；通過(guò)接觸向?qū)Ы⒔佑|對(duì)；約束法蘭盤(pán)內(nèi)徑的徑向位移；約束法蘭盤(pán)內(nèi)徑的周向位移；施加預(yù)拉伸力；求解；后期處理。3. ANSYS計(jì)算的具體流程3.1. 分析模型在開(kāi)始有限元計(jì)算前，需要對(duì)計(jì)算的工程問(wèn)題

7、，進(jìn)行認(rèn)真的分析。這其中包括：模型的簡(jiǎn)化，能否忽略幾何不規(guī)則性，能否把三維問(wèn)題簡(jiǎn)化為平面問(wèn)題；能否把一些載荷看作是集中載荷；能否把某些支撐看做固定的；判斷模型材料的應(yīng)力-應(yīng)變關(guān)系，即采用線(xiàn)性材料分析還是非線(xiàn)性材料分析等。3.2. 選擇單元根據(jù)第一步分析的結(jié)果選擇滿(mǎn)足條件的單元，ANSYS中有多種單元，點(diǎn)、面、實(shí)體等均有，而且每種單元又有較多類(lèi)型，因此在具體選擇時(shí)應(yīng)結(jié)合使用者對(duì)軟件的熟悉程度結(jié)合所分析研究模型進(jìn)行選擇。在單元選擇時(shí)，可結(jié)合具體選項(xiàng)適當(dāng)講解各種典型單元的用途和常用的一些單元【4】。螺栓法蘭蓋單元設(shè)置步驟:Preprocessor ElementTypeAdd/Edit/Delete

8、Addsolid20node 186OK如圖3-1。設(shè)置襯墊單元設(shè)置步驟：Preprocessor ElementTypeAdd/Edit/DeleteAddGasket 16node 194OKClose。兩種單元設(shè)置完以后如圖32所示。31設(shè)置186單元體32設(shè)置完兩種單元體3.3. 設(shè)置材料屬性材料屬性需要結(jié)合所分析研究對(duì)象的特性來(lái)設(shè)置，主要包括泊松比、楊氏模量，屬性材料包括屈服強(qiáng)度和切線(xiàn)模量（雙線(xiàn)型材料），疲勞分析則有應(yīng)力-工作次數(shù)關(guān)系曲線(xiàn)等，因此對(duì)構(gòu)件的材料掌握非常關(guān)鍵。本出采用非線(xiàn)性線(xiàn)型材料來(lái)研究材料的特性。1）螺栓和法蘭蓋材料設(shè)置步驟：PreprocessorMaterialMo

9、delsStructural ElasticIsotropic OK,然在在彈出的對(duì)話(huà)框里輸入彈性模量1.05E11，泊松比0.3。如圖3-3所示。2）襯墊材料的設(shè)置步驟：PreprocessorMaterialModelsStructuralSpecialized MaterialsGasketCmpresionAdd point,如圖3-4所示。把數(shù)值輸入表中得到襯墊材料店壓縮曲線(xiàn)如圖3-5，再點(diǎn)Linaer unloain得到卸載曲線(xiàn)，擊如圖3-6,7所示。33螺栓法蘭蓋材料設(shè)置34襯墊材料設(shè)置35襯墊壓縮曲線(xiàn) 36襯墊材料卸載數(shù)據(jù)點(diǎn) 37襯墊材料卸載特性曲線(xiàn)3.4. 建立模型ANSYS

10、中有四種建立模型方法，即實(shí)體建模，有限元建模，從其他CAD軟件中導(dǎo)入和參數(shù)化建模。前三種建模方法適合除優(yōu)化設(shè)計(jì)，可靠性分析的各種分析。建模過(guò)程：首先設(shè)置參數(shù)變量如圖3-8所示。H1=0.01 !法蘭盤(pán)1的厚度。H2=0.002 !襯墊厚度。H3=0.01 !法蘭盤(pán)2的厚度。H4=0.008 !螺帽的厚度。R1=0.05 !法蘭盤(pán)內(nèi)半徑。R2=0.08 !法蘭盤(pán)外半徑。R3=(R2-R1)/2+R1 !螺栓中心位置。R4=0.5*(R2-R1)/2!螺栓半徑。R5=0.7*(R2-R1)/2!螺帽半徑。38引入建模參數(shù)1）建立法蘭蓋和襯墊模型，通過(guò)圓柱命令建立不同高度的圓柱， 3-9所示圖，輸入

11、數(shù)字，建立1襯墊，2個(gè)法蘭的模型。操作步驟：PreprocessorModelingCreateVolumesCylinderBy DimensionsOk 39法蘭盤(pán)與襯墊三維模型2）通過(guò)布爾減操作，建立法蘭盤(pán)和襯墊上的四個(gè)圓孔，如圖3-10所示。操作步驟：PreprocessorModelingOperateBooleansSubtract選擇建立的襯墊與法蘭蓋apply選擇4個(gè)小的圓柱OK.3）通過(guò)創(chuàng)建圓柱操作，對(duì)建立的3個(gè)圓柱體進(jìn)行布爾相加，得到螺栓，如圖3-11所示。操作步驟：PreprocessorModelingCreateVolumesCylinderBy Dimension

12、4）把的螺栓通過(guò)平移到法蘭蓋的孔處，如圖3-12。再把螺栓通過(guò)復(fù)制粘貼到其他三個(gè)法蘭孔的位置，如圖3-13。操作步驟：WorkplaneChange active CS toGloable Cylidric PreprocessorModelingCopyVolumes再在DY欄輸入粘貼角度90,180,270。 310創(chuàng)建孔 311創(chuàng)建螺栓 312螺栓偏移 313四個(gè)螺栓創(chuàng)建完成3.5. 網(wǎng)格劃分ANSYS中主要有兩種網(wǎng)格劃分方法，即自由網(wǎng)格和映射網(wǎng)格。自由網(wǎng)格劃分的成功率高，但在動(dòng)力學(xué)計(jì)算中精度稍差，用戶(hù)根據(jù)需要進(jìn)行選擇使用。1）劃分一個(gè)法蘭蓋網(wǎng)格。法蘭蓋網(wǎng)格采用映射網(wǎng)格。操作步驟:Pre

13、processorMeshingMeshtoolSet186單元Set網(wǎng)格尺寸0.005HexSweep選擇一個(gè)法蘭蓋，如圖3-14，3-15。 314法蘭蓋網(wǎng)格參數(shù)設(shè)置 2)劃分襯墊網(wǎng)格。襯墊單元在兩個(gè)體的中間，所以先劃分的一個(gè)法蘭，再襯墊，再劃分另外一個(gè)法蘭蓋。操作步驟: PreprocessorMeshingMeshtoolSet186單元 Meshing Mesh Interface3D Interface選擇41,42號(hào)面再在如圖力輸入值，完成網(wǎng)格劃分，如圖3-16,3-17所示。4）劃分另一個(gè)法蘭蓋 另一個(gè)法蘭蓋的劃分方法與第一法蘭蓋的劃分網(wǎng)格方法一樣，網(wǎng)格劃分完如圖3-18所示。

14、5)劃分螺栓網(wǎng)格。螺栓網(wǎng)格采用自由網(wǎng)格劃分。操作步驟:PreprocessorMeshingMeshtoolSet186單元TetFreeMesh選擇4個(gè)螺栓Ok，如圖3-19所示。 315 劃分法蘭蓋 316襯墊網(wǎng)格參數(shù)設(shè)置317 襯墊單元網(wǎng)格劃分完成法 318蘭螺栓劃分完成 319劃分螺栓3.6. 確定分析類(lèi)型這一步用戶(hù)可以選擇需要的分析類(lèi)型：靜力學(xué)分析，模態(tài)分析，諧響應(yīng)分析，瞬態(tài)動(dòng)力學(xué)分析，譜分析，特征值屈曲分析，子結(jié)構(gòu)分析等，3.7. 施加邊界條件這一步用戶(hù)根據(jù)模型的實(shí)際工況定義邊界條件，如對(duì)稱(chēng)邊界，完全約束等。1）定義預(yù)拉伸單元操作步驟：:PreprocessorMoleingEle

15、mentspretensionPretensn MeshWith optionDivide at voluElements in volue。如圖3-20 320定義預(yù)拉伸單元2）定義接觸單元。螺栓底面與法蘭蓋表面由接觸。創(chuàng)法蘭表面的建節(jié)點(diǎn)組操作步驟：拾取法蘭表面節(jié)點(diǎn)SelectEntiliesAreasBy PickFrom AllSele AllApplyOk,如圖3-21所以。再SelectEntiliesNodesAttached toAreas AllFrom AllPlotOk,如圖3-22所示。創(chuàng)建法蘭表面節(jié)點(diǎn)組CN1，如圖322所示。用同樣創(chuàng)法蘭給另一面的節(jié)點(diǎn)組CN2。 321

16、創(chuàng)建節(jié)點(diǎn)組322創(chuàng)建節(jié)點(diǎn)組創(chuàng)建螺栓底面節(jié)點(diǎn)組操作步驟：拾取法蘭表面節(jié)點(diǎn)SelectEntiliesAreasBy PickFrom AllSele AllApply選擇四個(gè)螺栓底面Ok如圖3-23,再 Select Entilies Nodes Attached toAreas AllFrom AllPlotOk。 323拾取螺栓節(jié)點(diǎn)組 324螺栓和法蘭節(jié)點(diǎn)組 創(chuàng)建接觸步驟如下：Contact ManagerContace WizrdNodal ComponentTN1Nodal NexeCompentCN1Option SettingConcat OperitiesCreatFinsh.以同

17、樣的方法創(chuàng)建第二接觸。建立過(guò)程如圖3-25,3-26,3-27，3-28所示。325創(chuàng)建接觸326選擇接觸目標(biāo) 327設(shè)置不對(duì)稱(chēng) 328 接觸設(shè)置完成3.8. 求解用戶(hù)在完成以上操作后，就進(jìn)入求解階段，ANSYS有直接求解，載荷步求解和自適應(yīng)求解方法，根據(jù)需要適時(shí)選擇。通過(guò)約束法蘭內(nèi)徑的軸向和軸向位移。所以約束R1一周的節(jié)點(diǎn)。選取法蘭內(nèi)徑節(jié)點(diǎn)：WokplaneChange WorkplaneChange active CS toGloable Cylidric SelectEntiliesAreasBy locationR1OK.過(guò)程設(shè)置如圖。329選取法蘭內(nèi)環(huán)節(jié)點(diǎn)約束法蘭內(nèi)徑節(jié)點(diǎn)：Prepr

18、ocessorSolutionApplyStrucuralOn NodesPick all Ux,Uy。如圖3-30,3-31所示。 330設(shè)置步 331施加約束. 332約束后的顯示 定義預(yù)拉伸載荷SolutionDefine LoadapplyStructuralPretnsn secth. 333設(shè)置載荷 334載荷設(shè)置過(guò)程求解SolveCurrent Ls3.9. 后處理在該步用戶(hù)可以使用兩種后處理方法，觀(guān)察結(jié)果，校核計(jì)算結(jié)果是否符合工程結(jié)構(gòu)設(shè)計(jì)要求，如果滿(mǎn)足要求，則保存結(jié)果，如果不滿(mǎn)足要求，則需要修改模型重新計(jì)算直到符合要求為止。 觀(guān)察最大應(yīng)力點(diǎn)：General PostprocPo

19、lt ResulContour PoltElemnet Solution。如圖3-35所示。335整體變形圖觀(guān)察螺栓應(yīng)力步驟：SelectEntiliesVolumBy pick螺栓OkSelectEntiliesElementAttacth toVolumOk,如圖所示3-36所示。用同樣方法觀(guān)察法蘭，襯墊。如圖3-37,3-38所示。336螺栓應(yīng)力圖337法蘭應(yīng)力圖338襯墊應(yīng)力圖4. 結(jié)論通過(guò)對(duì)螺栓襯墊法蘭的的力學(xué)分學(xué)模型分析，再借助anasys的有限元計(jì)算，得出：襯墊在孔的位置處，應(yīng)力最大。螺栓，在有形狀急劇變的話(huà)的位置處，應(yīng)力集中，應(yīng)力最大。本模型還可以施加不同的預(yù)緊力，去觀(guān)察螺栓，

20、襯墊，法蘭的變形與應(yīng)力狀況。因此可以對(duì)施加的預(yù)緊力進(jìn)行控制，得到理想的預(yù)緊力。參考文獻(xiàn)1 王平.襯墊的壓應(yīng)力分析.石油礦廠(chǎng)機(jī)械，20012年2月，20(1):5051 2 2張宏英，李東陽(yáng).螺栓預(yù)緊力分析， 機(jī)械裝備，2011，40(3)：2023.3 肖文生，ANASYS13.0基礎(chǔ)與實(shí)例教程.2011年11月.4 秦宇.ANSYS11.0基礎(chǔ)與實(shí)例教程，化學(xué)工業(yè)出版社，2009，21（2）：48.命令流通過(guò)對(duì)求解的模型導(dǎo)出命令流。通過(guò)修改命令流里面的參數(shù)，可以計(jì)算同模型下，不同問(wèn)題的求解。命令流如下：/BATCH ! /COM,ANSYS RELEASE 14.5 UP20120918 0

21、9:30:28 01/03/2014 /input,start145,ans,D:ANSYS Incv145ANSYSapdl !* /FILNAME,xieshuai2,0*SET,H1,0.01 ! *SET,H2,0.002 ! *SET,H3,0.01 ! *SET,H4,0.008 ! *SET,R1,0.05 ! *SET,R2,0.08 ! *SET,R3,(R2-R1)/2+R1 ! *SET,R4,0.5*(R2-R1)/2! *SET,R5,0.7*(R2-R1)/2! /PREP7 !* ET,1,SOLID186 !* ET,2,INTER194 !* !* MPTEM

22、P, MPTEMP,1,0 MPDATA,EX,1,1.05E11MPDATA,PRXY,1,0.3 TB,GASKET,2,13,COMP ! TBPT,0.508000E-04, 0.161226E+07TBPT,0.101600E-03, 0.520884E+07TBPT,0.152400E-03, 0.113134E+08TBPT,0.203200E-03, 0.200499E+08TBPT,0.254000E-03, 0.259960E+08TBPT,0.304800E-03, 0.290345E+08TBPT,0.355600E-03, 0.357453E+08TBPT,0.406

23、400E-03, 0.440064E+08TBPT,0.457200E-03, 0.563189E+08TBPT,0.508000E-03, 0.748254E+08TBPT,0.558800E-03, 0.967287E+08TBPT,0.609600E-03, 0.129001E+09TBPT,0.683260E-03, 0.157147E+09TB,GASKET,2,5,LUNL ! TBPT,0.152400E-03, 2.430000E+11TBPT,0.304800E-03, 3.565000E+11TBPT,0.406400E-03, 5.923000E+11TBPT,0.558

24、800E-03, 1.088000E+12TBPT,0.683260E-03, 1.490000E+12!* CYLIND,R1,R2,0,H1,0,360,CYLIND,R1,R2,H1,H1+H2,0,360,CYLIND,R1,R2,H1+H2,H1+H2+H3,0,360, NUMMRG,ALL CYLIND,R4,0,0,H1+H2+H3,0,360, ! VGEN,4,R3,1 ! CSYS,1 ! ! VGEN,2,4,90,0 VGEN,2,4,180,0 VGEN,2,4,270,0 ! FLST,2,3,6,ORDE,2 FITEM,2,1 FITEM,2,-3 VSBV,

25、P51X,4 FLST,2,3,6,ORDE,2 FITEM,2,8 FITEM,2,-10 VSBV,P51X,5 FLST,2,3,6,ORDE,2 FITEM,2,1 FITEM,2,-3 VSBV,P51X,6 FLST,2,3,6,ORDE,3 FITEM,2,4 FITEM,2,-5 FITEM,2,8 VSBV,P51X,7 CYLIND,R4*0.95,0,0,H1+H2+H3,0,360, CYLIND,R5,0,H1+H2+H3,H1+H2+H3+H4,0,360, CYLIND,R5,0,-H4,0,0,360,! /VIEW,1,1,1,1 ! /ANG,1 ! /RE

26、P,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /AN

27、G,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,

28、1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST FLST,2,3,6,ORDE,2 FITEM,2,4 FITEM,2,-6 VADD,P51X CSYS,0 ! VGEN,7,R3,1 CSYS,1 VGEN,2,7,90,0 VGEN,2,7,180,0 VGEN,2,7,270,0 ! SAVE, xieshuai2,db, ! SAVE, xieshuai2,db, ! SAVE, xieshuai2,db, ! SAVE, 36,db,C:UsersAdministrator ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAS

29、T ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST TYPE, 1 MAT, 1REAL, ESYS, 0 SECNUM, !* ESIZE,0.005,0, /UI,MESH,OFFESIZE,0.005 TYPE,1 MAT,1 VSWEEP,2TYPE,2 MAT,2 IMESH,AREA,19,25, ,0,0,H2,0.0005TYPE,1 MAT,1 VSWEEP,3MSHKEY,0MSHAPE,1,3D VMESH,4 VMESH,5 VMESH,6 VMESH,7 PSMESH,

30、 ,P1, ,VOLU,4, 0,Z,(H1+H2+H3)/2, PSMESH, ,P2, ,VOLU,5, 0,Z,(H1+H2+H3)/2, , , , , PSMESH, ,P3, ,VOLU,6, 0,Z,(H1+H2+H3)/2, , , , , PSMESH, ,P4, ,VOLU,7, 0,Z,(H1+H2+H3)/2, , , , , ! /DIST,1,0.924021086472,1! /REP,FAST ASEL,ALLASEL,S, , , 48 NSEL,ALLNSLA,S,1! NPLOT NSLA,S,1CM,CN1,NODE ! /DIST,1,0.924021

31、086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08

32、222638492,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! vlist, all ! elist,all,0,0 ! /DIST,1,0.92402108

33、6472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.0822

34、2638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.

35、924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ASEL,S, , ,28

36、 NSLA,S,1CM,CN2,NODE ALLSEL,ALL ! /MREP,EPLOT ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ASEL,ALLFLST,5,4,5,ORDE,4 FITEM,5,50 FITEM,5,54 FITEM,5,64 FITEM,5,74 ASEL,S, , ,P51X NSEL,ALLNSLA,S,1! NPLOT ! NPLOT NSLA,S,1! /MREP,EPLOT

37、CM,TN1,NODE FLST,5,4,5,ORDE,4 FITEM,5,70 FITEM,5,80 FITEM,5,53 FITEM,5,60 ASEL,S, , ,P51X NSEL,ALLNSLA,S,1CM,TN2,NODE ALLSEL,ALL ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.0822263

38、8492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /MREP,EPLOT ! /COM, CONTACT PAIR CREATION - START C

39、M,_NODECM,NODE CM,_ELEMCM,ELEM CM,_KPCM,KP CM,_LINECM,LINE CM,_AREACM,AREA CM,_VOLUCM,VOLU ! /GSAV,cwz,gsav,tempMP,MU,1,0.25MAT,1 MP,EMIS,1,7.88860905221e-031R,3 REAL,3 ET,4,170ET,5,174R,3,1.0,0.1,0,RMORE,1.0E20,0.0,1.0, RMORE,0.0,0,1.0,1.0,0.5RMORE,0,1.0,1.0,0.0,1.0KEYOPT,5,4,0KEYOPT,5,5,0NROPT,UNS

40、YM KEYOPT,5,7,0KEYOPT,5,8,0KEYOPT,5,9,0KEYOPT,5,10,2 KEYOPT,5,11,0 KEYOPT,5,12,0 KEYOPT,5,2,0KEYOPT,4,5,0! Generate the target surface NSEL,S,TN1CM,_TARGET,NODE TYPE,4 ESLN,S,0ESURF CMSEL,S,_ELEMCM ! Generate the contact surface NSEL,S,CN1CM,_CONTACT,NODETYPE,5 ESLN,S,0ESURF ALLSEL ESEL,ALLESEL,S,TY

41、PE,4 ESEL,A,TYPE,5 ESEL,R,REAL,3 ! /PSYMB,ESYS,1 ! /PNUM,TYPE,1! /NUM,1 ! EPLOT ESEL,ALLESEL,S,TYPE,4 ESEL,A,TYPE,5 ESEL,R,REAL,3 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACM CMSEL,S,_VOLUCM CMDEL,_V

42、OLUCM ! /GRES,cwz,gsav CMDEL,_TARGET CMDEL,_CONTACT ! /COM, CONTACT PAIR CREATION - END ! /COM, CONTACT PAIR CREATION - START CM,_NODECM,NODE CM,_ELEMCM,ELEM CM,_KPCM,KP CM,_LINECM,LINE CM,_AREACM,AREA CM,_VOLUCM,VOLU ! /GSAV,cwz,gsav,tempMP,MU,1,0.25MAT,1 MP,EMIS,1,7.88860905221e-031R,4 REAL,4 ET,6

43、,170ET,7,174R,4,1.0,0.1,0,RMORE,1.0E20,0.0,1.0, RMORE,0.0,0,1.0,1.0,0.5RMORE,0,1.0,1.0,0.0,1.0KEYOPT,7,4,0KEYOPT,7,5,0NROPT,UNSYM KEYOPT,7,7,0KEYOPT,7,8,0KEYOPT,7,9,0KEYOPT,7,10,2 KEYOPT,7,11,0 KEYOPT,7,12,0 KEYOPT,7,2,0KEYOPT,6,5,0! Generate the target surface NSEL,S,TN2CM,_TARGET,NODE TYPE,6 ESLN,

44、S,0ESURF CMSEL,S,_ELEMCM ! Generate the contact surface NSEL,S,CN2CM,_CONTACT,NODETYPE,7 ESLN,S,0ESURF ALLSEL ESEL,ALLESEL,S,TYPE,6 ESEL,A,TYPE,7 ESEL,R,REAL,4 ! /PSYMB,ESYS,1 ! /PNUM,TYPE,1! /NUM,1 ! EPLOT ESEL,ALLESEL,S,TYPE,6 ESEL,A,TYPE,7 ESEL,R,REAL,4 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEM

45、CM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACM CMSEL,S,_VOLUCM CMDEL,_VOLUCM ! /GRES,cwz,gsav CMDEL,_TARGET CMDEL,_CONTACT ! /COM, CONTACT PAIR CREATION - END ! /MREP,EPLOT ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST CSYS,1 NSEL,ALLNSEL,S,LOC,X,R1 ! N

46、PLOT NSEL,S,LOC,X,R1 ! /REPLOT,RESIZE FLST,2,1024,1,ORDE,12 FITEM,2,2473FITEM,2,-2536 FITEM,2,2809FITEM,2,-3124 FITEM,2,3753FITEM,2,-3884 FITEM,2,6297FITEM,2,-6424 FITEM,2,9017FITEM,2,-9332 FITEM,2,9913FITEM,2,-9980 NROTAT,P51X ! /REPLOT,RESIZE FLST,2,1024,1,ORDE,12 FITEM,2,2473FITEM,2,-2536 FITEM,2

47、,2809FITEM,2,-3124 FITEM,2,3753FITEM,2,-3884 FITEM,2,6297FITEM,2,-6424 FITEM,2,9017FITEM,2,-9332 FITEM,2,9913FITEM,2,-9980 NROTAT,P51X ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST FINISH /SOLFLST,2,1024,1,ORDE,12 FITEM,2,2473FITEM,2,-2536 FITEM,2,2809FITEM,2,-3124 FITEM,2,3753FITEM,2,-3884 FITEM,2,6297FITEM

48、,2,-6424 FITEM,2,9017FITEM,2,-9332 FITEM,2,9913FITEM,2,-9980 !* /GO D,P51X, , , , , ,UX,UY, , , , ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ALLSEL,ALL CSYS,0 ! vlist, all ! elist,all,0,0 ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /D

49、IST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST

50、! vlist, all ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,