Fletcher-2-rus (1185919), страница 23
Текст из файла (страница 23)
(см. рис. 13.28), в результате чего находится набор значений Г, в зависимости от хл, где хл — локальная координата профиля, 0 < <хл <1. Полная длина поверхности между точками А и С определяется выражением ГАС, ааа«ГАВ + ХС вЂ” ХВ. (13.?2) Однако функция растяжения (!3.44) дает нормализованный па- РаМЕтР Га таКОй, Чта 0 ~ (Глс ~~ 1 при 0 ( $ ( 1. (13.73) 13.4.1. А(.ОЕМ: программа построения сеток около хороито обтекаемых тел Изображенная на рис. 13.25 поверхность АВ является аэродинамическим профилем семейства МАСА-00'1'.
Здесь ! — двухзначное число, представляющее толщину профиля в процентах. Так, профиль МАСА-00!2 является симметричным профилем, толщина которого составляет 12 о1о от единичной хорды. Профили семейства МАСА-00'У задаются уравнением у =1(а х Га+ а х + а ха+ а хз+ а х'), (13.70) 136 Гл. 13. Построение сеток Следовательно, физическая поверхностная координата глс, о вычисляется по формуле ГЛС, П ГЛСГЛС, Н 44 (13.74)4 где глс рассчитывается в подпрограмме 5ТГсЕСН (см. рис.!3.27), а глс,~ек определяется из (13.72). При глс и ( г„е физические координаты (х(1),у(1)) граничного сегмента АВ получаются в результате интерполяции гл и хл и использования (13.70) для расчета соответствующих значений у.
Это делается в подпрограмме Е01Е при !ЯТ = 1. При глс, о ) гли координаты поверхности получаются из линейной интерполяции хв ( х ( хс. Физические координаты поверхности АВС в программе АЕОЕМ (рис. 13.26) обозначены через Х5(1,7), т'Ь(1, У), поскольку они являются компонентами х и у поверхности Х, из п. 13.3.3. Аналогично в результате интерполяции нормализованной поверхностной координаты гго получаются значения поверхностной координаты гго, о и физические координаты Х5 (4, У) „ У5(4,7) граничной поверхности РЕР. Поверхностная координата гго, о отсчитывается вдоль дуги окружности РЕ и прямой линии ЕР.
Положение точек сетки на АР и СР (рис. 13.25) зависит от параметров растяжения злг и зоо. Однако физические значения координат точек сетки на АР и СР получаются явным образом в результате работы алгоритма многих поверхностей (13.76). В рассматриваемой программе в методе многих поверхностей (п. 13.3.3) используются четыре поверхности.
Граничные поверхности АВС и РЕР образуют поверхности Х, и Х4. Исходные промежуточные поверхности Хе и Хе получаются в результате линейной интерполяции между Х, и Х4 по формулам Хе = Х, + зе(Х4 — Х,), Х44 = Х, + з,(Х4 — Х,). (13.75) На данном этапе поверхности Хе и Хе правильно располагаются в физическом пространстве, однако корректное соответствие точек по координате г еще не установлено. Если алгоритм многих поверхностей (13.76) применить с промежуточными поверхностями Хе и Хз, определенными (13.75), то 41-линии полученной 4 сетки будут прямыми, соединяющими точки с одинаковыми значениями индекса 1 на АВС и РЕР (см. рис.
13.24 при У= 2). Следовательно, зависимость Хе и Хе от г должна быть изменена так, чтобы линии, проведенные через точки на Хе и соответствующие точкам (точки с тем же значением г) на Хь были перпендикулярны ХП линии, соединяющие Хе и Х,, должны быть перпендикулярны Х,. В этом случае алгоритм многих поверхно- 9 )Згл Численная реализация алгебраического отображения [37 ЛС ПС 3 С 4 С 5 С 6 7 '$ 9 С 10 11 12 13 14 15 16 17 18 19, 20 21 22 23 24 25 С 26 С 27 С 2$ 29 30 31 С 32 С 33 С 34 35 С Зб 37 С ЗЗ 39 С 40 41 С 42 34 45 46 47 48 49 50 С 51 С 52 С 53 54 С льбен АРРыез л но01$1$9 ноьтт[4)-ЗОРРЛСЕ тксйнломе ТО тие Окиеалт10И 03' А ОЛ10 АВОмт А ИАсл-00'т' леа0301ь Ат Екйо Хнс1юемсе.
тик МРРВВ НАЬР ОВ10 1$6ейейлтев ПЭИКИ$10И Х(53,5И,Т(51,51),ХВ(б),ТВ(6),1$(4,51),ТЗ(4,51) ВЭИКНЗХОИ йлс(51),ВРВ[51),ЗАР(51),ЗОВ(51),ЗН(4),ХА[51),йл(51) СОМНОИ ФА,ХЛ ОРЕМ(1,13ЬЕ АЬОЕН.ВЛТ ) ОРБИ(6,11ЬК 'ЛЬЗЕН.ООТ') йелн(1,1)днах,хнах,1РВ,1АРь,т,$2,$з,лм 1 Роанлт(415,бк10.3). ВЕЛО(3,2)Рлс,алс,РРВ,СРВ,РЛР,САР,Рсвдса 2 Роййлт(8$(о.з) МВ1ТЕ(6,3) ма1те(6,4)дйлх,ийх,1РВ.7,$2,$3,АМ Э РОЮВТ(' НМЬТ1$0ВРЛСЕ ОВ10 6ЕИЕВЛТЭОИ') 4 РОВИАТ(' дйлк ЕНЛХ ',213,' 1РВ ',12, 1' Т ',$10.3,5Х,' $2.$3"',ЭЕ30.3,' ЛМ ',РЗ.Э,IЛ . НА17$(6,5)РлсАлс,РРВ,ОРВ МВ1ТЕ В, 6) РАР. САР, РСВ, ЯСВ 5 Роанлт( Рлсй,кло,'3.
Олс ~',В10.3, ° РРВ ',Б10.3,' ВРВ '.Е10.3) б РОВНЛТ(' РЛР ',$10.3,' ПАР ',Б10.3,' РСВ~',ИЗО.Э,' ССПМ',Е10,3) ВБР1НВ соаикй $01нт$ ОР тне Вомнвлат ВАТА ХВ!2.00,3.00,5.00,5.00,2.25,0.00/ ВАТА ТВ/0.00,0.00,0.00,2.25,2.25,0.00/ Р1 Э.14159265 6КНШТЕ $ТЯЕТСВ1ИО РМНСТЭОИЗ СЛЬЬ Зтййсй(д)ВХ,РАС,САС,ВАС) сльь атакой(дилх,РРВ,ОРП,ВРВ» СЛЬЬ ЗТАБСИ(ХНАХ,РАР,СЛР,$АР) САЬЬ ЗТВБСН(ййХ,РСВ,ССВ,ЗСВ) МВ1ТЕ(6,7) ()ВС(д) .3~1,дйлй) МВЭТБ[6,8) (ЯРОВ),д Э~дйлк) МВ1ТЕ[6.9] (ЗАР(д) д 1.дМЛХ) МВ1ТЕ(6,10) (ЗСВЫ),д 3,0НАХ) 7 РОПВТ(( йЛС=',18$7.4) '8 РОВНЛТ[' ВРВ~',1817.4) 9 РОА3йТ(' ЗЛР ',1817,4) ЭО РОВВТ[' зов=','3817.4,3) овтлТИ $0атлск СООВВЭНАТБЗ оР вовт, лв Сльь Р01Ь[О,Т.ВАВ,ХВ,ТВ) Рис. [3.26.
Распечатка программы АсСЕМ (начало). ]38 Гл. )3. Построение сеток иснх - ив + хв<з> - хв<г> КГЕ 0.5*Г1*(УВ(5) — УВ(1]) КГОНХ КГЕ + ХВ(4) — ХВ(5) 1Г(198 .ЬО. о]сото 13 нкзте<6,11> <и(5],ь 1.5Л НК1ТЕ(6,12) (КЯ(1],Ь 1,51) Гокит(( хх ',18Г7.4] Гокнят< кя=,)ЗГ7.4] СЯ\ А Г01Ь(1,Т,КАСО,ХО,УО) СЕИЕКАТЕ 1ИТЕКХОН СЯ1О Рнс. (3.26 (продолжение). 55 56 57 58 59 60 61 62 63 С 64 С 65 С 66 67 68 69 70 71 72 С и 74 С 75 76 77 78 79 80 $1 82 83 Вб 85 86 Вт С 88 С 89 С ЗО С 91 92 С 93 94 95 96 97 98 99 1ОО 101 102 1ОЗ 104 С 105 С 106 С 11 12 13 14 15 16 17 1В 19 20 21 СЕИЕЯЯТХ ВООМО1МС ЗОКГАСЕЗ 1 АНО 4 ОО И д = 1.СНЯХ ИСВ = КЯС(д] НАСИХ 1Г(КАСО .Ет. КАВ)СОТО 14 Х5(1,Л ХВ(2) + (ИСО " ИВ)*(ХВ(3) - ХВ(2))/<ИСМХ-$АВ) УЗ(1,д] = О.
СЬТО и ХЗ(1,д) = ХВ*(ХВ(2]-ХВ(1)) + ХВ(1) Т5(1,Л = УО*(ХВ(2]-ХВ(1]) + УВ(1) КГОО АГО<д)*КГВНХ 1Г<ЯГОО .Ет. ЯГИСОТО 16 ХЗ(4,Л ХВ(5) + (НГОО-КГЕ]*(ХВ(4)-ХВ(5))/(КГОМХ-КГЕ] 75(4,д) = УВ(5) СОТО 17 КК УВ(5) — ТВ(1) тне = кГОО/ка ХЗ(4,Л ° ХВ(5] " КК*СОЗ(ТНЕ) УЗ(4,д) КК*$1И(ТНЕ] соитзное БОКСЕ СЕИЕКЯТЕЗ ЗОКГАСЕ5 2 Анв 3 $0 ТНАТ ТНЕ СК1О АВдАСЕИТ ТО 50КГАСЕ5 1 АИО 4 1$0ЯТНОСОИАЬ СААЬ ЗОКСВ<ОНАХ,$2,$3,ХЗ,УЗ) 00 21 5 1.4 НК1ТЕ(6,18)1 ГОКИАТ(1ОН ЗОЯГЯСЕ ,11) ОО 21 д = 1,дИАХ,18 дк д дВ дк+ 17 НК1ТЕ (б, 19] <ХЗ (Ь.дс),дс=дв,дв) нктте(6, го) <ув <Е.дс) .дс и,дв) ГОНИТ(' ХЗ ', 18Г7.4) ГОКНАТ(' И ',гвхт.б, и СОИТ1ИОЕ Гл.
(3. Построение сеток (40 ВУВВООТ1ИЕ ЭТВЕСВ(И,Р,Я,В) НТМЕН$10И $ (51) АИ И-1 ВЕТА 1./АИ ТЯ1 и:1./ТАИН(Я) 301 Ь 1,М .АП*~ Ь 1 ЕТА н АьаУЕТА ° УУН и Я" (1. - ЕТА) иук . - ВАМИ(ВОМ)*тя1 $(Ь) Р*ЕТА + (1.-Р)"РУМ 1 СОНТ1МУЕ ВЕТУВИ ° ЕИУ отей (и. 13.3.3) обеспечивает ортогональность сетки вблизи поверхностей ХФ и ХФ. Поскольку УФ(г) известно, достаточно провести через каждую точку сетки ХЗ(!,У), Ф'8(1,У) прямую линию, перпендикулярную ХФ, до пересечения с УВ.Точками пересечения будут ХБ(2,У), "Ф'8(2, У). Данное построение эквивалентно этапу предиктор, изображенному на рис. 13.17. Поскольку задание узла У определяет значение г на ХФ, точки ХЗ(2, У), т'3(2, У) имеют то же значение г и ортогональны к ХФ в точках ХЗ(1, У), Ф'8(1, У). Расчет ортогональных точек сетки на Ха осуществляется в подпрограмме 3()ВСН (рис.
13.29). Точки сетки ХЗ(З, У), (8(3, У) на поверхности Уа строятся так, что они оказываются ортогональными поверхности ХФ в точках ХЗ(4, У), т'8(4, У). Расчет этих точек также осуществляется в подпрограмме 8!)ЙСН. При заданных координатах на поверхностях от ХФ до 24 алгоритм многих поверхностей в случае четырех поверхностей осуществляется следующим образом !Е!зетап, 1979): 4 4 х(У, К)= ~., ЗН(У.)ХЗ(У,, У), У(У, К)= ~' ЗН(У.) Ф'8(У., У), (!3.76) где ЗН (1 ) = ( 1 — $)' ( 1 — а,з), ЗН (2) = $ ( 1 — $)' (а, + 2), ЗН(3)=$$(1 — $)$(аа+ 2), ЗН(4)=$$(1 — аа(1 — $)), (13.77) а, = 2/(За — 1), аа = 2/(2 — Зам).
(13.78) 1 Э С Э С 4 С $ С 6 С 7 В 9 10 11 С 11 13 16 16 16 17' 1$ 19 10 СОМРУТЕ$0ИЕ-У1НЕН$10ИАЬ ВТВЕТСН1ИО РУНСТ10И, $ Ратта + (1.-Р)*(1.-ТИФЕ(Я*(1."ЕТА))/ТМФЕ(Я))а РОВ $1теи соитаоь Ри(и(ете$5, Р Аиу я. Рис. (3,27, Распечатка программы 3Т((ЕСН. 1 2С эс 40 5С б С ТС ЯС 9 ЗО 11 12 зз 34 С 15 С 16 С 17 Ха 19 ЗО И 22 гэ 24 25 26' 27 ге 29 зо Эг 32 зз Эб 35 36 Зт зв ЗЗ С 4О С 41 С 41 С 43 44 4$ бб 47 46 69 $О $1 52 $3 54 5$ $6 $7 56 воаяоОтгие РОХьиит,т,ив,х,т) 50ВИСЕ РВОЕ1ЬЕ 13 ИАСЛ-00'Т' АхйОЕ01Ь 1Р 1ит ° О, %)меВ1сАььт 1итебВАтя то ОВта1М 5<Зисе СООВВ1МАТЕ 1Е 1МТ 1, 1ИТЕЯРОЬАТЕ 50ЯЕАСВ СООЯР1ИАТВ5 ТО ОВТАХИ СОАВЕЗРОИВ1НО (Х,Т) В1ИЯИ51ОИ А(5),ХЛ(ЗВ,ВА(51) СОННОМ ВА.ХА ВАТА Л<1.4779155,-0.624424,-1.727016,1.364037,-0.439769/ РХ 3.14159265 хе<Хит .ХО.
Иоото 2 ]пжея1сАььт 1мтеойатя то ОВта1и йх(ь) А5 А еоистХОМ Ое и(ь) ЯА(1) м О. ХА<1) О. ХА<2) = (А<1) l(1./Т " Ы2])1**2 ВЛ(2] О. $~РХ~ХА(2) ВОН 3.*А(4) + 4.*А(5)*ХЛ[21 ВШ< 2.*А(3) + ВОН*ХА(2) ВЯИ Т*(0.$<ЗОЯТ(ХА(21) + А(2] + ВШ<*И<2]) Еьт Явят<1. + ВОН*Вою ВО 1 Ь 2,50 ЛЬ Ь ЬР Ь+1 ха<ьи - 0.0253ь ВХ ХЫЬР) ХА(В ' Еь " ЕЬР, В<Ж ,З.*А(4] Ф 4.*Л($]*И(ЬР) РОМ 2.*А(З) + РОН'ХА<ЬИ ВОМ Т*(0<5ЩЯТ(И(ЬР)) + А(2) + РВИ*И(ЬР)) ЕЬР войт<1. + РОМ*ВОН) ВА(ЬР) ВА(Ь] + 0.5 ° (ЕЬ + ЕЬР)*РХ 1 СОИТПШЕ ВЛВ ВЛ(511 ВВТОВИ ХМТЕВРО<ЛТЕ ЯА(В ТО ОВТА1И Х СОЯВЕЗРОИВХИО ТО ЙАВ эцвэеоовитьт овтии т ейои лильттгс илса-оо т яюетье 2 ВО 3 Ь 2,$1 1ЕШАВ .ОТ. ВА(Ы)ООТО 3 ЬН Ь-1 Х И<ЬИ) + <И<Ш-И<ЛИ)) <ЯАВ-ВЛ(ЬИ)ИВА<В-и[ЛИ)) 1Е(Х .ЬТ. 1.0х-об)Х 1.0х-об В<Ж А(4) + Л<$)*х РШ< А(31 + ВШ<"Х В<)Н ° А(2) + РОМ*Х т та<Л<11*ВРВТ[11 + Ввиах) Ветойи 3 СОМТХИОН иазте[6,4)ив,И<1),йл<$31 4 ЕОВИАТ(' ВАВ ООТ51ВЕ ВАНОВ',51.' ВАР*',ХХО.З, 1' йа(11 '.Е10.3,' ВА(51) ° ',Е10.31 ВЕТ<ЖИ ЮВ Рас.
33.28. Распечатка программы РОП.. 142 Гл. 13. Построение сеток 1 2 С 3 С ЗС 5 С б 7 $ С 9 С 10 С 11 12 13 14 15 16 17 1$ 19 С 20 С 21 С 22 23 24 25 26 27 28 29 ЗО Э1 Э2 33 Э4 35 36 37 38 39 40 41 42 43 44 46 С 46 С 47 С аз 49 50 61 52 53 54 55 ЗОВЯООТ1НЕ ЗОЯСВ(диАХ,52,5Э,Х5,75) ОЕНЕЯАТ5$ ЗОЯРАСЕЗ 2 АИВ 3 ТО СВЕАТЕ ОЯТВ060НАЬ ВООНВАВУ 65105 01НЕИ510$ Х$2 (51],752 (51],ХБ3 (51),ТБ3 (51),ХБ (4,51),УЗ (4,51) ЗНАР 3НАХ - 1 Ряеьгнгияяу Окнкяят1ои от ЯОВРАскз 2 АНВ э ВО 1 д Э,днях ВХ5 ХБ(4,Л " ХЗП,Л ВТЗ 75[4,Л вЂ” УЗП,Л ХЗ(2,3) = ХЗ П,д) + $2"ОХЗ УЗ(2,д] = ТЗП,Л + 52*075 Х5(З,Л Х5 П,д) + $3*ВХ5 УЗ(З,д) " 75 П,д) + ЯЭ*ВУЯ 1 СОИТХИОК Ряодест оятноооияььу Раои БОВРАск 1 онто ЗОВРАск 2 ВО 9 д 2,3ИАР 1РНВЗ(ХЗП,д+1)-Х5П,У"1)1 .ОТ.
1.0Е"06)ООТО 2 ЕН1 1.0506'(УЗ (1,д+1) -ТЯП, д"1) ) сото 3 2 ЕН1 (УЗП,д+1]-ХЗП,3-1])([ХБП,З+1) "ХЗП,У-1)) 3 1Р(АВЗ(ХЗ[2,Л-Х5(2,У-1)) .ОТ. 1.0$-06)СОТО 4 ЕН2 1.05~06'[УЗ(2,Л-УЗ(2,3-1)] ООТО 5 4 кнг = Гуз(г,л-тз[г,д-и)((хзп,л-хзд„у-л> 5 Хг = (КНЭ*(тЗ(Э,Л-УЗП,Л+Ккг*ХЗП,Л)+ХЗП,Л)((1.+Вид Кнг> У2 УБ (2,3) + ЕН2" [Х2 - ХЗ (2, д) ] ЗТЛ[ 50ЯТ ((Х2-Х5 (2,д-1) ) **2 + (Уг-тз (2,д-1) ) **2) Зддн ЗОЯТ((ХБ(2,Л-ХБ(2,д-1))**2 + (75(2,3)-75[2,3-1))**2) 1Р(5ТЭН .ЬТ. Здди)ООТО 8 1Р(АВ5[Х5(2,д+1]"ХБ(2,3)) .ОТ.
1.0Е-ОЯ)60ТО б ЕН2 1.05+06*[У5[2,3+1>-75(2,Л ) 60ТО 7 б ЕН2 (75(2,д+1)-ТЗ[2,3!)((Х5(2,3+1)-ХЗ(2,3)) 7 Х2 (ЕН1*(75П,З)-У5(2,7)+ЕН2*ХЗ(2,д))+ХБП,Л)(П.+ЕК1*ЕН21 72 У5[2,3) + ЕИ2 (Х2-Х5(2,3)) 8 Х52(д! Х2 752(д) У2 9 СОИТ1ИОЕ РЯОдЕСТ ОЯТВОООНАЬЬУ РВОН 56$РАСЕ 4 ОНТО БПЯРАСВ 3 ВО 17 д 2.3НАР 1Р($85(Х5(4,д+1)-хз(4,3-1]) .ОТ.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
