Диссертация (1025976), страница 16
Текст из файла (страница 16)
A thesis submitted for the degreeof Doctor of Philosophy. Adelaide. 2012. 226 p.112. Hempsell M., Bond A., Bond R. Progress on the SKYLON and SABREDevelopment Programme, 62nd International Astronautical Congress, CapeTown, IAC-11.D 2.4.2, 2011. 7 p.113. Hempsell M. Progress on the SKYLON and SABRE. 64th InternationalAstronautical Congress, Beijing, China, IAC-13,D2.4,6x19609, 2014. 14 p.114. Holland W.
SSME/SIDE Loads Analysis for Flight Configuration. NASA TMX-64841 Revision A, Sep. 1974. 90 p.115. Huzel D.K., Huang D.H. Design of Liquid Propellant Rocket Engines,Washington, D.C., NASA, 1967. 460 p.116. Jankovsky R.S., Arya V.K., Kazaroff J.A. Structurally Compliant RocketEngine Combustion Chamber—Experimental and Analytical Validation. NASATP-3431, March 1994. 12 p.117. Kasper H.J.
Thrust Chamber Life Prediction. NASA CP-2372, 1984, pp. 36-43.118. Korotkaya O.V. Substructure Method for Thermal-Stress Analysis of LiquidPropellant Rocket Engine Combustion Chamber. World Academy of Science,Engineering and Technology, International Science Index 88, InternationalJournal of Mechanical, Aerospace, Industrial and Mechatronics Engineering,2014, no.
8(4), pp. 734 - 737.163119. Kuhl D., Riccius J., Haidn O.J. Thermomechanical Analysis and Optimizationof Cryogenic Liquid Rocket Engines// Journal of propulsion and power, Vol. 18,No. 4, July–August 2002. pp. 835-846.120. Lal S., Asraff A.K., Thomas S.E. Stress strain analysis of rocket engine thrustchamber. Indian journal of applied research, Vol. 5, Issue 10, October 2015,pp. 225-228.121. Linne D.L.
A Rocket Engine for Mars Sample Return Using In Situ Propellants.35th Aerospace Sciences Meeting and Exhibit, Reno, Nevada, January 6-10,1997. 12 p.122. Luscher W.P., Mellisch J.A. Advanced High Pressure Engine Study for MixedMode Vehicle Applications. NASA CR-135141, Jan. 1977. 304 p.123. Meyer M.L. Design Issues for Lunar In Situ Aluminum/Oxygen PropellantRocket Engines. 1992 Aerospace Design Conference, Irvine, California,February 3-6, 1992.
21 p.124. Meyer M.L., Arrington L.A., Kleinhenz J.E. Testing of a Liquid Oxygen/LiquidMethane Reaction Control Thruster in a New Altitude Rocket Engine TestFacility. Space Propulsion 2012, Bordeaux, France, May 7–10, 2012. 14 p.125. Naraghi M.H.N., Armstrong E.S. Three Dimensional Thermal Analysis ofRocket Thrust Chambers. Thermophysics, Plasmadynamics and LasersConference, San Antonio, Texas, June 27-29, 1988. 32 p.126.
Naraghi M.H.N. A Computer Code for Three-Dimensional Rocket ThermalEvaluation. New York, Tara Technologies, 2002. 106 p.127. Naraghi M.H., Dunn S, Coats D. A Model for Design and Analysis ofRegeneratively Cooled Rocket Engines. AIAA-2004-3852, Joint PropulsionConference, Fort Lauderdale, 2004. 13 p.128. Oberth H. Die Rakete zu den Planetenräumen. München und Berlin, R.Oldenbourg, 1923. 92 p.129. Page B.R. Rocket Experiments of Robert H. Goddard, 1911-1930. The PhysicsTeacher, Nov.
1991, pp. 490-496164130. Pelaccio D., Jacobs M., Scheil C. Engine System Assessment Study UsingMartian Propellants. NASA CR-189188, June 1992. 339 p.131. Price H.G. Cooling of High Pressure Rocket Thrust Chambers with LiquidOxygen. 16th Joint Propulsion Conference, Hartford, Connecticut, June 30-July2 1980. 15 p.132. Price H.G., Masters P.A. Liquid Oxygen Cooling of High PressureLOX/Hydrocarbon Rocket Thrust Chambers.
NASA Technical Memorandum88805. 1986. 26 p.133. pSeven software. URL. https://www.datadvance.net/ (дата обращения:14.06.2015).134. Riccius J.R., Zametaev E.B. Stationary and Dynamic Thermal Analyses ofCryogenicLiquidRocketAIAA/ASME/SAE/ASEEJointCombustionChamberPropulsionConferenceWalls.and38thExhibit,Indianapolis, Indiana, July 7–10, 2002. 12 p.135.
Roncace E.A. In-Situ Propellant Rocket Engines for Mars Mission AscentVehicle. 27th Joint Propulsion Conference, Sacramento, California, AIAA Paper91-2445, 1991. 8 p.136. Rousar D.C., Miller F. Cooling with Supercritical Oxygen. AIAA paper 751248, Sept. 1975. 6 p.137. Ruf J.H., McDaniels D.M., Brown A.M. Nozzle Side Load Testing and Analysisat Marshall Space Flight Center. 45th AIAA/ASME/SAE/ASEE Joint PropulsionConference and Exhibit, Denver, Colorado, Aug.
2-5, 2009. 14 p.138. Schmucker R.H. Flow Processes in Overexpanded Chemical Rocket NozzlesPart 3: Methods for the Aimed Flow Separation and Side Load Reduction.NASA TM-77048, March 1983. 23 p.139. Schmucker R.H. Flow Processes in Overexpanded Chemical Rocket NozzlesPart 2: Side Loads due to Asymmetric Separation. NASA TM-77395, Feb. 1984.57 p.165140. Shi J.J. Rocket engine nozzle side load transient analysis methodology - apractical approach. 46th AIAA/ASME/ASCE/AHS/ASC Structures, StructuralDynamics and Materials Conference, Austin, Texas, 18-21 April 2005.
11 p.141. Smalley K.B., Brown A., Ruf J. Flow Separation Side Loads Excitation ofRocket Nozzle FEM. 48th AIAA/ASME/ASCE/AHS/ASC Structures, StructuralDynamics and Materials Conference, Honolulu, Hawaii, Apr. 23-26, 2007. 32 p.142. Spencer R.G., Rousar D.C. Supercritical Oxygen Heat Transfer. NASA CR135339, Nov. 1977. 101 p.143. Stubbs V.R. A graphite-lined regeneratively cooled thrust chamber. NASA CR120853, Jan. 1972. 260 p.144. Sutton G.P., Biblarz O. Rocket Propulsion Elements. New York, John Wileyand Sons, 2001. 751 p.145.
Von Braun W. History of Rocketry and Space Travel. New York, Thomas Y.Crowell Co., 1966. 244 p.146. Wall N. Characterisation of multiple concentric vortices in hybrid rocketcombustion chambers. A thesis submitted for the degree of Doctor ofPhilosophy. Sheffield. 2013. 269 p.147. Wang Ten-See. Transient Two-Dimensional Analysis of Side Load in LiquidRocket Engine Nozzles. 40th AIAA/ASME/SAE/ASEE Joint PropulsionConference and Exhibit, Fort Lauderdale, Florida, 11-14 July 2004. 11 p.148. Wang Ten-See. Transient Three-Dimensional Startup Side Load Analysis of aRegeneratively Cooled Nozzle. 38th Fluid Dynamics Conference and Exhibit,Seattle, Washington, June 23-26, 2008. 22 p.149.
Wang Ten-See, Lin J., Ruf J. Transient Three-Dimensional Side Load Analysisof Out-of-Round Film Cooled Nozzles. 46th AIAA/ASME/SAE/ASEE JointPropulsion Conference and Exhibit, Nashville, Tennessee, July 25-28, 2010.57 p.150. Wang Ten-See, Lin J., Ruf J. Transient Side Load Analysis of Out-of-RoundFilm-Cooled Nozzle Extensions. 48th AIAA/ASME/SAE/ASEE Joint PropulsionConference and Exhibit, Atlanta, Georgia, Jul. 30-Aug. 1, 2012.
29 p.166151. Zienkiewicz O. C., Taylor R. L. The finite element method: Solid mechanics.Oxford: Butterworth-heinemann, 2000. 459 p.167ПриложениеП.1. Тексты авторских программВ данном разделе приложения представлены тексты авторских программ,описанных в Главе 3.П. 1.1. ОптимизацияТекст программы, решающей последовательно стационарную задачутеплопроводности и прочностную задачу с заданными граничными условиямидля расчета оптимальных геометрических параметров тракта охлажденияпредставлен ниже./BATCH !Запуск ANSYS в BATCH-режимеK,8,r2,0,0!Геометрические параметрыK,9,r2,beta,0n=300K,10,r1+t,beta,0r1=352.2e-3K,11,r1,alpha,0t=1.2e-3K,12,r2,alpha,0h=1.3e-3h1=1.2e-3L,1,2alpha=0.882L,1,11beta=360/nL,11,3r2=r1+t+h+h1L,4,10/PREP7L,3,10!Построение геометрической моделиL,5,10CSYS,1L,5,6K,100,0,0,0L,6,7K,1,r1,0,0L,7,8K,2,r1+t,0,0L,8,12K,3,r1,beta,0L,5,9K,4,r1+t,alpha,0L,4,10K,5,r1+t+h,beta,0L,4,2K,6,r1+t+h,alpha,0CSYS,0K,7,r1+t+h,0,0L,4,6168CSYS,1CM,_Y,AREAL,9,12ASEL, , , ,P51XCSYS,0CM,_Y1,AREAL,4,11CMSEL,S,_YL,6,12CMSEL,S,_Y1AL,1,2,15,12AATT,AL,3,5,15,4CMSEL,S,_YAL,4,6,7,13CMDELE,_YAL,7,11,14,16CMDELE,_Y1AL,8,9,10,16!КЭ разбиение3, , ,0,FLST,5,10,4,ORDE,7!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!FITEM,5,1*Блок описания моделей всехFITEM,5,-2используемых материалов с учетомFITEM,5,5пластического поведения*FITEM,5,8!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!FITEM,5,-12!Присвоение материалов площадямFITEM,5,15!Силовая оболочкаFITEM,5,-16FLST,5,2,5,ORDE,2CM,_Y,LINEFITEM,5,4LSEL, , , ,P51XFITEM,5,-5CM,_Y1,LINECM,_Y,AREACMSEL,,_YASEL, , , ,P51X!*CM,_Y1,AREALESIZE,_Y1,0.000125, , , , , , ,1CMSEL,S,_YFLST,5,6,4,ORDE,6!*FITEM,5,3CMSEL,S,_Y1FITEM,5,-4AATT,5, , ,0,FITEM,5,6CMSEL,S,_YFITEM,5,-7CMDELE,_YFITEM,5,13CMDELE,_Y1FITEM,5,-14!Огневая стенкаCM,_Y,LINEFLST,5,3,5,ORDE,2LSEL, , , ,P51XFITEM,5,1CM,_Y1,LINEFITEM,5,-3CMSEL,,_Y169!*FLST,2,1,4,ORDE,1LESIZE,_Y1,0.00009, , , , , , ,1FITEM,2,12!*DL,P51X, ,TEMP,120ET,1,PLANE77,0,0,1 !Тип КЭMSHAPE,0,2DFLST,2,1,4,ORDE,1MSHKEY,1FITEM,2,8!*DL,P51X, ,TEMP,30FLST,5,5,5,ORDE,2FITEM,5,1!Линейное распределение температурыFITEM,5,-5*DIM,temp1,TABLE,2,1,1,X,,CM,_Y,AREA*SET,TEMP1(1,0,1) , r1+tASEL, , , ,P51X*SET,TEMP1(1,1,1) , 120CM,_Y1,AREA*SET,TEMP1(2,0,1) , r1+t+hCHKMSH,'AREA'*SET,TEMP1(2,1,1) , 30CMSEL,S,_Y!*FLST,2,1,4,ORDE,1AMESH,_Y1FITEM,2,13!*!*CMDELE,_YDL,P51X, ,TEMP, %TEMP1%CMDELE,_Y1CMDELE,_Y2FLST,2,1,4,ORDE,1FITEM,2,6TOFFST,273!*DL,P51X, ,TEMP, %TEMP1%!Плоская деформацияKEYOPT,1,1,0!Решение температурной задачиKEYOPT,1,3,0FINISH/SOL!Температура/STATUS,SOLU!Равномерное распределение температурыSOLVEFLST,2,2,4,ORDE,2FITEM,2,2/PREP7FITEM,2,-3!Изменение типа КЭ на structuralDL,P51X, ,TEMP,180ETCHG,TTS170!Плоская деформацияFITEM,2,13KEYOPT,1,3,2FITEM,2,6KEYOPT,1,5,0/GOKEYOPT,1,6,0!*!*SFL,P51X,PRES,19.5e6,!Закрепления!Температурная нагрузкаFLST,2,2,4,ORDE,2LDREAD,TEMP,,, , ,'sech2pldef','rth',' 'FITEM,2,1FITEM,2,9!Решение!*FINISH/GO/SOLDL,P51X, ,UY,/STATUS,SOLUSOLVEFLST,2,2,4,ORDE,2FINISHFITEM,2,2FITEM,2,3/POST1!*/GOSET,LASTDL,P51X, ,UX,/CYCEXPAND, ,ON!Циклическая симметрия!Получение максимального эквивалентногоCYCLIC, , , ,'CYCLIC'напряжения и максимальной пластическойдеформации!ДавлениеETABLE,seqv,S,EQVFLST,2,2,4,ORDE,2ESORT,ETAB,seqv,0,0FITEM,2,2*GET,max_seqv,SORT,0,MAXFITEM,2,3max_seqv=max_seqv/1000000/GO!*ETABLE,defpl,EPPL,EQVSFL,P51X,PRES,0.05e6,ESORT,ETAB,defpl,0,0*GET,max_defpl,SORT,0,MAXFLST,2,4,4,ORDE,4max_defpl=max_defpl*100FITEM,2,8FITEM,2,12/EXIT, MODEL171П.1.2.
Коэффициенты анизотропииПараметрическое задание геометрии на языке APDL:n=90!число каналов охлажденияt=1.3e-3!толщина огневой стенкиh=1.8e-3!толщина ребер охлажденияh1=2.8e-3 !толщина силовой оболочкиalpha=2.09 !угол, характеризующий ширину канала охлаждения в окружном направленииr1=31.6e-3 !внутренний радиусgamma=75 !угол закручивания каналов охлажденияbeta=360/n !угол сектораB=(r1+t+h)*(beta*3.14159/180)*tan(gamma*3.14159/180) !высота представительного элементаc=tan(gamma*3.14159/180)*((beta-alpha)/2)*3.14159/180*(r1+t+h) !вспомогательная величина/PREP7K,29,r1+t+h+h1,0,0K,22,r1+t+h,0,B-c!Построение ключевых точекK,32,r1+t+h+h1,beta,0K,24,r1+t+h,beta,B-cCSYS,1!цилиндрическаяK,23,r1+t,beta,B-cсистема координатK,11,r1,0,BK,25,r1+t,0,cK,100,0,0,0K,12,r1,beta,BK,26,r1+t+h,0,cK,13,r1+t,0,BK,28,r1+t+h,beta,cK,1,r1,0,0K,14,r1+t+h,0,BK,27,r1+t,beta,cK,2,r1,beta,0K,15,r1+t+h,(beta-alpha)/2,B!Построение линийK,3,r1+t,0,0K,16,r1+t,(beta-alpha)/2,BL,1,3K,4,r1+t+h,0,0K,17,r1+t,(beta+alpha)/2,BL,3,4K,5,r1+t+h,(beta-alpha)/2,0K,18,r1+t+h,(beta+alpha)/2,BL,4,5K,6,r1+t,(beta-alpha)/2,0K,19,r1+t+h,beta,BL,5,6K,7,r1+t,(beta+alpha)/2,0K,20,r1+t,beta,BL,3,6K,8,r1+t+h,(beta+alpha)/2,0K,30,r1+t+h+h1,0,BL,6,7K,9,r1+t+h,beta,0K,31,r1+t+h+h1,beta,BL,7,8K,10,r1+t,beta,0K,21,r1+t,0,B-cL,8,9172L,9,10L,13,21AL,45,30,47,12L,7,10L,21,22AL,1,45,19,39,46,33L,10,2L,22,14AL,11,47,29,44,48,38L,1,2L,19,24AL,39,23,53L,4,29L,24,23AL,46,53,24,55L,29,32L,23,20AL,33,55,28,44,57,5L,9,32L,1,11AL,6,57,48,59L,14,30L,25,21AL,59,38,10L,30,31L,2,12AL,39,20,41,40L,31,19L,27,23AL,33,34,35,2L,11,13L,26,22AL,44,27,42,43L,13,14L,28,24AL,38,37,36,9L,14,15L,29,30AL,53,22,54,40L,15,16L,32,31AL,7,59,37,60L,13,16L,21,16AL,34,55,25,56L,16,17L,22,15AL,4,57,43,58L,17,18L,25,17AL,35,49,41,16,51,13L,18,19L,26,18AL,14,52,17,51L,19,20L,6,23AL,15,36,50,42,18,52L,17,20L,5,24AL,8,60,36L,20,12L,7,27AL,41,21,54L,11,12L,8,28AL,49,54,32,56L,5,8!Построение площадейAL,3,35,56,26,42,58L,15,18AL,12,11,10,6,5,1AL,31,58,50,60L,3,25AL,5,2,3,4!Построение объемовL,25,26AL,10,9,8,7VA,9,5,10,1,11,16,15,14,13,12L,4,26AL,3,31,8,15,14,13VA,25,8,26,27,4,29,30,31,32,28L,10,9AL,19,30,29,28,24,23VA,17,6,29,21,12L,9,28AL,23,22,21,20VA,28,20,16,3,22L,28,27AL,28,27,26,25VA,2,7,14,18,19,23,24,31L,27,10AL,21,32,26,18,17,16173Текставторскойпрограммыдлявычислениякоэффициентованизотропии представлен ниже:!ТемператураTEMP=20!Параметры геометрииn=90 !число каналов охлажденияt=1.3e-3 !толщина огневой стенкиh=1.8e-3 !толщина ребер охлажденияh1=2.8e-3 !толщина силовой оболочкиalpha=2.09 !угол, характеризующий ширину канала охлаждения в окружном направленииr1=31.6e-3 !внутренний радиусgamma=75 !угол закручивания каналов охлажденияbeta=360/n !угол сектораB=(r1+t+h)*(beta*3.14159/180)*tan(gamma*3.14159/180)!высотапредставительногоэлемента!Модуль Юнга приведенного материала в меридиональном направлении, далее онизменяется в циклеEpriv=5e+10flag=0 !Если 0 - продолжать нагружать, если 1 - не нагружать, сразу следующий шаг цикла/PREP7K,12,r1,beta,B!Построение ключевых точекK,13,r1+t,0,BCSYS,1!цилиндрическаясистемаK,14,r1+t+h,0,BкоординатK,19,r1+t+h,beta,BK,100,0,0,0K,20,r1+t,beta,BK,1,r1,0,0K,29,r1+t+h+h1,0,0K,2,r1,beta,0K,30,r1+t+h+h1,0,BK,3,r1+t,0,0K,31,r1+t+h+h1,beta,BK,4,r1+t+h,0,0K,32,r1+t+h+h1,beta,0K,9,r1+t+h,beta,0K,102,r1,beta/2,0K,10,r1+t,beta,0K,103,r1+t+h+h1,beta/2,0K,11,r1,0,BK,104,r1,beta/2,B174K,105,r1+t+h+h1,beta/2,BL,4,9K,106,r1,0,B/2L,3,13K,107,r1,beta,B/2L,10,20K,108,r1+t+h+h1,0,B/2L,4,14K,109,r1+t+h+h1,beta,B/2L,9,19K,110,r1+(t+h+h1)/2,0,0L,102,103K,111,r1+(t+h+h1)/2,0,BL,103,105K,112,r1+(t+h+h1)/2,beta,BL,104,102K,113,r1+(t+h+h1)/2,beta,0L,105,104!Построение линийL,106,107L,1,3L,107,109L,3,4L,109,108L,9,10L,108,106L,10,2L,113,112L,1,2L,112,111L,4,29L,111,110L,29,32L,110,113L,9,32!Построение площадейL,14,30AL,1,5,4,23L,30,31AL,2,23,3,24L,31,19AL,6,24,8,7L,11,13AL,12,16,15,21L,13,14AL,13,21,14,22L,19,20AL,9,22,11,10L,20,12AL,5,17,16,18L,11,12AL,19,10,20,7L,10,9AL,1,17,12,25L,1,11AL,18,15,26,4L,2,12AL,2,25,13,27L,29,30AL,3,26,14,28L,32,31AL,6,27,9,19L,13,20AL,8,28,11,20L,14,19AL,23,25,21,26L,3,10AL,24,27,22,28175AL,31,32,30,29MPTEMP,5,300AL,33,34,35,36MPTEMP,6,400AL,37,38,39,40MPTEMP,7,500!Построение объемовMPTEMP,8,600VA,1,4,7,9,10,15MPDATA,KXX,3,,2.2VA,2,5,11,12,15,16MPDATA,KXX,3,,8.7VA,3,6,8,13,14,16MPDATA,KXX,3,,16.3!Дополнительные построенияMPDATA,KXX,3,,18.6FLST,2,3,6,ORDE,2MPDATA,KXX,3,,22FITEM,2,1MPDATA,KXX,3,,22.1FITEM,2,-3MPDATA,KXX,3,,24.9VSBA,P51X,17MPDATA,KXX,3,,26.6FLST,2,6,6,ORDE,2!Модель материала 12Х18Н10Т (модельFITEM,2,4№3)FITEM,2,-9MPTEMP,,,,,,,,VSBA,P51X,18MPTEMP,1,-253FLST,2,12,6,ORDE,4MPTEMP,2,-196FITEM,2,1MPTEMP,3,-70FITEM,2,-3MPTEMP,4,20FITEM,2,10MPTEMP,5,300FITEM,2,-18MPTEMP,6,500VSBA,P51X,19MPTEMP,7,700!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!MPDATA,EX,3,,2.2301E+011!Модели используемых в конструкцииMPDATA,EX,3,,2.11E+011материаловMPDATA,EX,3,,2.085E+011MPDATA,EX,3,,2.02E+011!Таблица коэффициентовMPDATA,EX,3,,1.82e+11теплопроводности 12Х18Н10Т (модельMPDATA,EX,3,,1.67e+11№3) коллекторMPDATA,EX,3,,1.53e+11MPTEMP,,,,,,,,MPDATA,PRXY,3,,0.3MPTEMP,1,-250MPDATA,PRXY,3,,0.3MPTEMP,2,-190MPDATA,PRXY,3,,0.3MPTEMP,3,100MPDATA,PRXY,3,,0.3MPTEMP,4,200MPDATA,PRXY,3,,0.3176MPDATA,PRXY,3,,0.3MPTEMP,2,200MPDATA,PRXY,3,,0.3MPTEMP,3,300MPTEMP,,,,,,,,MPTEMP,4,400MPTEMP,1,-253MPTEMP,5,500MPTEMP,2,-196MPTEMP,6,600MPTEMP,3,-100MPTEMP,7,700MPTEMP,4,-50MPTEMP,8,800MPTEMP,5,100MPDATA,KXX,5,,17.598MPTEMP,6,200MPDATA,KXX,5,,18.017MPTEMP,7,300MPDATA,KXX,5,,19.274MPTEMP,8,400MPDATA,KXX,5,,20.112MPTEMP,9,500MPDATA,KXX,5,,23.045MPTEMP,10,600MPDATA,KXX,5,,27.235MPTEMP,11,700MPDATA,KXX,5,,30.168MPTEMP,12,800MPDATA,KXX,5,,35.196UIMP,1,REFT,,,20!Модель материала (модель №5)MPDATA,ALPX,3,,11.2e-612Х21Н5Т-ВДMPDATA,ALPX,3,,13.0e-6MPTEMP,,,,,,,,MPDATA,ALPX,3,,15.0e-6MPTEMP,1,20MPDATA,ALPX,3,,15.2e-6MPTEMP,2,300MPDATA,ALPX,3,,16.6e-6MPTEMP,3,400MPDATA,ALPX,3,,17.0e-6MPDATA,EX,5,,1.76E+011MPDATA,ALPX,3,,17.4e-6MPDATA,EX,5,,1.53E+011MPDATA,ALPX,3,,17.8e-6MPDATA,EX,5,,1.45E+011MPDATA,ALPX,3,,18.5e-6MPDATA,PRXY,5,,0.3MPDATA,ALPX,3,,18.7e-6MPDATA,PRXY,5,,0.3MPDATA,ALPX,3,,18.7e-6MPDATA,PRXY,5,,0.3MPDATA,ALPX,3,,18.8e-6MPTEMP,,,,,,,,!Таблица коэффициентовMPTEMP,1,100теплопроводности (модель №5)MPTEMP,2,20012Х21Н5Т-ВД форсуночная головкаMPTEMP,3,300MPTEMP,,,,,,,,MPTEMP,4,400MPTEMP,1,100MPTEMP,5,500177UIMP,1,REFT,,,20MPTEMP,3,200MPDATA,ALPX,5,,11.1e-6MPTEMP,4,300MPDATA,ALPX,5,,14.0e-6MPTEMP,5,400MPDATA,ALPX,5,,16.8e-6MPTEMP,6,500MPDATA,ALPX,5,,16.8e-6MPTEMP,7,550MPDATA,ALPX,5,,17.4e-6MPTEMP,8,600!Таблица коэффициентовMPTEMP,9,650теплопроводности для ЭП750-ШMPTEMP,10,700(модель №1) силовая оболочка камерыMPDATA,KXX,2,,170MPTEMP,,,,,,,,MPDATA,KXX,2,,195MPTEMP,1,-253MPDATA,KXX,2,,220MPTEMP,2,-193MPDATA,KXX,2,,235MPTEMP,3,20MPDATA,KXX,2,,245MPTEMP,4,100MPDATA,KXX,2,,285MPTEMP,5,200MPDATA,KXX,2,,300MPTEMP,6,300MPDATA,KXX,2,,300MPTEMP,7,400MPDATA,KXX,2,,295MPTEMP,8,500MPDATA,KXX,2,,295MPTEMP,9,600!Таблица коэффициентовMPDATA,KXX,1,,2.9теплопроводности (модель №7)MPDATA,KXX,1,,7.5ХН55МБЮMPDATA,KXX,1,,18.0!внутренняя проставка соплаMPDATA,KXX,1,,19.5MPTEMP,,,,,,,,MPDATA,KXX,1,,21.5MPTEMP,1,-253MPDATA,KXX,1,,23.5MPTEMP,2,-223MPDATA,KXX,1,,25.5MPTEMP,3,-173MPDATA,KXX,1,,27.5MPTEMP,4,-113MPDATA,KXX,1,,29.5MPTEMP,5,-73!Таблица коэффициентовMPTEMP,6,-13темплопроводности БрХЦрТ (модель №2)MPTEMP,7,20! внутренняя оболочка камерыMPTEMP,8,400MPTEMP,,,,,,,,MPTEMP,9,500MPTEMP,1,20MPTEMP,10,600MPTEMP,2,100MPTEMP,11,700178MPTEMP,12,800MPDATA,EX,1,,1.71E+011MPTEMP,13,900MPDATA,EX,1,,1.605E+011MPDATA,KXX,7,,1.9MPDATA,EX,1,,1.512E+011MPDATA,KXX,7,,5.6MPDATA,EX,1,,1.45E+011MPDATA,KXX,7,,9.8MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,12.0MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,13.0MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,14.0MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,14.5MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,14.5MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,17.4MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,18.9MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,20.0MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,22.1MPDATA,PRXY,1,,0.3MPDATA,KXX,7,,23.5MPTEMP,,,,,,,,!Модель материала ЭП750-Ш (модельMPTEMP,1,-253№1)MPTEMP,2,-193MPTEMP,,,,,,,,MPTEMP,3,0MPTEMP,1,-253MPTEMP,4,200MPTEMP,2,-196MPTEMP,5,300MPTEMP,3,20MPTEMP,6,400MPTEMP,4,100MPTEMP,7,500MPTEMP,5,200MPTEMP,8,600MPTEMP,6,300MPTEMP,9,700MPTEMP,7,400MPTEMP,10,800MPTEMP,8,500MPTEMP,11,900MPTEMP,9,600UIMP,1,REFT,,,20MPTEMP,10,700MPDATA,ALPX,1,,10.6e-6MPDATA,EX,1,,2.09E+011MPDATA,ALPX,1,,12.9e-6MPDATA,EX,1,,2.07E+011MPDATA,ALPX,1,,16.4e-6MPDATA,EX,1,,1.94E+011MPDATA,ALPX,1,,16.5e-6MPDATA,EX,1,,1.89E+011MPDATA,ALPX,1,,16.8e-6MPDATA,EX,1,,1.855E+011MPDATA,ALPX,1,,17.5e-6MPDATA,EX,1,,1.8E+011MPDATA,ALPX,1,,18.1e-6179MPDATA,ALPX,1,,18.4e-6MPTEMP,10,500MPDATA,ALPX,1,,18.5e-6MPTEMP,11,600MPDATA,ALPX,1,,18.7e-6MPTEMP,12,700MPDATA,ALPX,1,,18.9e-6MPTEMP,13,800!Модель материала БрХЦрТ (модель №2)UIMP,1,REFT,,,20MPTEMP,,,,,,,,MPDATA,ALPX,2,,0.5e-6MPTEMP,1,-70MPDATA,ALPX,2,,8e-6MPTEMP,2,20MPDATA,ALPX,2,,12.8e-6MPTEMP,3,400MPDATA,ALPX,2,,15.4e-6MPTEMP,4,500MPDATA,ALPX,2,,16.4e-6MPTEMP,5,600MPDATA,ALPX,2,,17.4e-6MPTEMP,6,700MPDATA,ALPX,2,,18.0e-6MPDATA,EX,2,,1.4489E+011MPDATA,ALPX,2,,18.4e-6MPDATA,EX,2,,1.4E+011MPDATA,ALPX,2,,18.5e-6MPDATA,EX,2,,1.246E+011MPDATA,ALPX,2,,18.8e-6MPDATA,EX,2,,1.17E+011MPDATA,ALPX,2,,19.0e-6MPDATA,EX,2,,1.18E+011MPDATA,ALPX,2,,19.2e-6MPDATA,EX,2,,1.06E+011MPDATA,ALPX,2,,19.6e-6MPDATA,PRXY,2,,0.3!Модель материала (модель №7)MPDATA,PRXY,2,,0.3ХН55МБЮMPDATA,PRXY,2,,0.3MPTEMP,,,,,,,,MPDATA,PRXY,2,,0.3MPTEMP,1,-253MPDATA,PRXY,2,,0.3MPTEMP,2,-196MPDATA,PRXY,2,,0.3MPTEMP,3,20MPTEMP,,,,,,,,MPTEMP,4,500MPTEMP,1,-253MPTEMP,5,600MPTEMP,2,-196MPTEMP,6,700MPTEMP,3,-153MPDATA,EX,7,,2.185E+011MPTEMP,4,-73MPDATA,EX,7,,2.13E+011MPTEMP,5,20MPDATA,EX,7,,2.01E+011MPTEMP,6,100MPDATA,EX,7,,1.708E+011MPTEMP,7,200MPDATA,EX,7,,1.644E+011MPTEMP,8,300MPDATA,EX,7,,1.56E+011MPTEMP,9,400MPDATA,PRXY,7,,0.3180MPDATA,PRXY,7,,0.3FLST,5,4,6,ORDE,3MPDATA,PRXY,7,,0.3FITEM,5,1MPDATA,PRXY,7,,0.3FITEM,5,-3MPDATA,PRXY,7,,0.3FITEM,5,10MPDATA,PRXY,7,,0.3CM,_Y,VOLUMPTEMP,,,,,,,,VSEL, , , ,P51XMPTEMP,1,-253CM,_Y1,VOLUMPTEMP,2,-203CMSEL,S,_YMPTEMP,3,-193!*MPTEMP,4,100CMSEL,S,_Y1MPTEMP,5,200VATT,MPTEMP,6,300CMSEL,S,_YMPTEMP,7,400CMDELE,_YMPTEMP,8,500CMDELE,_Y1MPTEMP,9,600!*MPTEMP,10,700!Силовая оболочкаMPTEMP,11,800FLST,5,4,6,ORDE,2MPTEMP,12,900FITEM,5,15UIMP,1,REFT,,,20FITEM,5,-18MPDATA,ALPX,7,,10.0e-6CM,_Y,VOLUMPDATA,ALPX,7,,11.7e-6VSEL, , , ,P51XMPDATA,ALPX,7,,12.0e-6CM,_Y1,VOLUMPDATA,ALPX,7,,13.9e-6CMSEL,S,_YMPDATA,ALPX,7,,14.1e-6!*MPDATA,ALPX,7,,14.1e-6CMSEL,S,_Y1MPDATA,ALPX,7,,14.3e-6VATT,MPDATA,ALPX,7,,14.5e-6CMSEL,S,_YMPDATA,ALPX,7,,15.1e-6CMDELE,_YMPDATA,ALPX,7,,15.5e-6CMDELE,_Y1MPDATA,ALPX,7,,16.3e-6!*MPDATA,ALPX,7,,16.3e-6!Модель материала (модель №702) -102, , 1,101, , 1,00!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!приведенный материал № 2 БрХЦрТ!Присвоение материалов объемамMPTEMP,,,,,,,,!Огневая стенкаMPTEMP,1,0181MPDATA,EX,702,,EprivMPTEMP,9,400MPDATA,PRXY,702,,0.2MPTEMP,10,500MPTEMP,11,600MPTEMP,,,,,,,,MPTEMP,12,700MPTEMP,1,20MPTEMP,13,800MPTEMP,2,100UIMP,1,REFT,,,20MPTEMP,3,200MPDATA,ALPX,702,,0.5e-6MPTEMP,4,300MPDATA,ALPX,702,,8e-6MPTEMP,5,400MPDATA,ALPX,702,,12.8e-6MPTEMP,6,500MPDATA,ALPX,702,,15.4e-6MPTEMP,7,550MPDATA,ALPX,702,,16.4e-6MPTEMP,8,600MPDATA,ALPX,702,,17.4e-6MPTEMP,9,650MPDATA,ALPX,702,,18.0e-6MPTEMP,10,700MPDATA,ALPX,702,,18.4e-6MPDATA,KXX,702,,170MPDATA,ALPX,702,,18.5e-6MPDATA,KXX,702,,195MPDATA,ALPX,702,,18.8e-6MPDATA,KXX,702,,220MPDATA,ALPX,702,,19.0e-6MPDATA,KXX,702,,235MPDATA,ALPX,702,,19.2e-6MPDATA,KXX,702,,245MPDATA,ALPX,702,,19.6e-6MPDATA,KXX,702,,285MPDATA,KXX,702,,300!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!MPDATA,KXX,702,,300!НАЧАЛО ЦИКЛАMPDATA,KXX,702,,295*DO,Epriv, 5e+9, 5e+10, 2.5e+9MPDATA,KXX,702,,295!Замена модуля Юнга в моделиMPTEMP,,,,,,,,приведенного материала на новый,MPTEMP,1,-253соответствующий шагу циклаMPTEMP,2,-196MPDATA,EX,702,,EprivMPTEMP,3,-153MPTEMP,4,-73!Присвоение объему "размазанного" слояMPTEMP,5,20приведенного материалаMPTEMP,6,100FLST,5,8,6,ORDE,4MPTEMP,7,200FITEM,5,4MPTEMP,8,300FITEM,5,-9182FITEM,5,19AESIZE,P51X,raz2,FITEM,5,-20CM,_Y,VOLU!РазбиениеVSEL, , , ,P51XMSHAPE,1,3DCM,_Y1,VOLUMSHKEY,0CMSEL,S,_Y!*!*FLST,5,16,6,ORDE,4CMSEL,S,_Y1FITEM,5,1VATT,702, , 1,0FITEM,5,-10CMSEL,S,_YFITEM,5,15CMDELE,_YFITEM,5,-20CMDELE,_Y1CM,_Y,VOLU!*VSEL, , , ,P51X!Тип конечного элементаCM,_Y1,VOLUET,1,SOLID187CHKMSH,'VOLU'!Настройка КЭ сеткиCMSEL,S,_Yraz1=0.00026!*raz2=0.00036VMESH,_Y1!*!Размеры КЭCMDELE,_YFLST,2,68,5,ORDE,13CMDELE,_Y1FITEM,2,1CMDELE,_Y2FITEM,2,-18FITEM,2,20!Если flag=0 - продолжать, flag=1 -FITEM,2,-28решение, очистка сетки и след.