Процедура вычисления

Описанная выше процедура позволяет получить ожидаемое значение глобального минимума функции качества, но еще не решает поставленной задачи — определения координат глобального экстремума, т. е. оптимальной комбинации параметров теплообменных аппаратов.

Описанная процедура позволяет установить требуемые уровни обслуживаемости, которые необходимо обеспечить, и связь их с вопросами надежности. Далее необходимо выбрать принцип обслуживания, удовлетворяющий предъявляемым требованиям к обслуживаемости.

Статистические решения для нескольких состояний. Выше были рассмотрены случаи, когда статистические решения принимались для различения двух состояний (дихотомия). Принципиально такая процедура позволяет провести разделение на п состояний, каждый раз объединяя результаты для состояния Dt и Dt. Здесь под Dt понимаются любые состояния, соответствующие условию «не ?>;». Однако в некоторых случаях представляет интерес рассмотреть вопрос и в прямой постановке — статистические решения для классификации п состояний. Рассмотрим, как и раньше, системы с одним диагностическим параметром х. Вначале остановимся на случае, когда проводится разделение на три состояния (рис. 9). Правило решения состоит в следующем:

Приведенные соображения подкрепляют ранее сделанный вывод о том, что двусторонняя процедура позволяет оценивать тепловое сопротивление (толщину) дефекта, но не его глубину залегания; это, в частности, делает невозможной двустороннюю тепловую томографию.

но-, так и двухстороннего ТК. Двухсторонняя процедура позволяет определить термическое сопротивление дефекта, а односторонняя - его глубину, причем в обоих случаях используют так называемый локально-одномерный подход, описанный в п. 3.1.2. Затем вводят концепцию пространственно-усредненного контраста Лапласа, который определяют по экспериментальным данным. Как показано в [28], значение этого контраста пропорционально Ъ х с. Данный алгоритм определения поперечных размеров дефектов является достаточно громоздким и представляет скорее теоретический интерес, однако и он может стать практическим по мере роста мощности компьютеров.

Такая процедура позволяет установить совокупность безразмерных чисел, характерных для изучаемого процесса. Эти числа подобия в общем случае являются мерой относительного влияния действующих сил и процессов переноса на течение жидкости и теплообмен.

Описанная выше процедура позволяет построить матрицы жесткости и этого элемента. Подобным же;образом определяют общую матрицу жесткости и находят напряжения. При аппроксимации поля перемещений полиномами второй степени (3.109) напряжения на границе элементов не будут иметь разрывов. Имеется возможность еще более уточнить решение, приняв на грани не одну, а две или более промежуточных точек. Матрицы жесткости элемента при этом становятся более громоздкими, но так как поле .перемещений и напряжений более точно, то можно значительно уменьшить число элементов во всей системе. Рассмотрим еще один двумерный конечный элемент применительно к задаче поперечного изгиба пластин (рис. 3.18). Деформация жест-

Рассмотренная итерационная процедура позволяет в нашем случае по заданным макродеформациям определить все неизвестные величины, относящиеся к микроуровню (структурные напряжения и деформации, значения функций микроповрежденности) и неизвестные

определение НДС упругих конструкций при неосесимметрич-ном нагружении и нагреве (при гармонической зависимости внешних нагрузок от времени процедура позволяет определить амплитудно-частотные характеристики конструкции);

определение НДС вязкоупругих конструкций при неосесим-метричном нагружении (при гармонической зависимости внешних нагрузок от времени процедура позволяет определить амплитудно-и фазочастотные характеристики конструкции);

определение НДС упругих конструкций при статическом на-гружении и нагреве (при гармонических внешних нагрузках процедура позволяет определять амплитудно-частотные характеристики конструкции);

F - fitfun(CC); /*fitfun - процедура вычисления значения целевой функции, реализует структурно-критериальную модель, СС - хромосома*/

Процедура вычисления вектора параметров состояния х -Do

Pi — процедура вычисления коэффициентов эквивалентного уравнения дискретной составляющей первого порядка с апериодическим переходным процессом; коэффициенты эквивалентного уравнения вычисляются по алгоритмам расчета динамических процессов в дискретных системах первого порядка;

Р\ — процедура вычисления коэффициентов эквивалентного уравнения дискретной составляющей первого порядка с колебательным переходным процессом; коэффициенты эквивалентного уравнения вычисляются по алгоритмам расчета динамических процессов в дискретных системах первого порядка;

Р3 — процедура вычисления коэффициентов эквивалентного уравнения дискретной составляющей второго порядка; коэффициенты эквивалентного уравнения вычисляются по алгоритмам расчета динамических процессов в дискретных системах второго порядка; в случае, когда в дискретной составляющей второго порядка, сформированной из пары дискретных звеньев, процесс оказывается апериодическим, в процедуре предусмотрена разбивка этой пары, при этом первое звено является составляющей первого порядка и рассматривается в процедуре Ръ а второе звено рассматривается совместно со следующим звеном;

Р4 — процедура вычисления коэффициентов эквивалентного уравнения дискретной составляющей первого порядка с однопи-ковым переходным процессом; коэффициенты эквивалентного уравнения вычисляются по алгоритмам расчета динамических процессов в дискретных системах первого порядка при ненулевых начальных условиях;

1) на экране дисплея САП УП выводит: шифр детали, формулу базы, номер базы, код поискового предписания (например, формула базы имеет вид Б1 = Л, Hi, II2, где Л — код вида базы; III — процедура вычисления диаметра базы; 112 — процедура вычисления расстояния от базы до правого торца детали;

В [2, 5, 10] содержатся приблизительные оценки средней, вызванной болезнью нетрудоспособности и среднего сокращения продолжительности жизни от наиболее обычных генетически детерминированных болезней в развитой стране. Для категорий болезней, индуцируемых излучением, оценки НКДАР ООН [2] приведены в табл. 2. Оценки генетического радиационного риска Международной комиссии по радиологической защите (МК.РЗ) (Публикация 45 [7] ) основаны на значениях, приведенных в табл. 2. Следует, однако, отметить, что процедура вычисления утраченных лет жизни как разности между средней продолжительностью жизни (в данном случае 70 лет) и возрастом в момент смерти является неверной, так как ожидаемая предстоящая жизнь человека в возрасте t больше разности между средней

К третьей группе относится критерий ортогональности плана, который требует такого выбора точек факторного пространства для измерений отклика, чтобы процедура вычисления оценок параметров модели данного вида была самой простой. В таком случае значения каждого параметра модели вычисляются независимо от значений других параметров.

Программа состоит из девяти процедур. 1. Процедура интегрирования — Procedure INT (х, у). Здесь х — подынтегральное выражение и у — значение интеграла. 2. Процедура решения интегрального уравнения методом последовательных приближений — procedure FUNCTION (х, у). Здесь обозначено х — функция FI; у — решение Ф(-. 3. Процедура итерполяции свойств материала по фактическим температурам диска — procedure INTERPOLATION. 4. Процедура нахождения нового секущего модуля и коэффициента Пуассона при упругопла-стическом решении — procedure PLASTICITY. 5. Процедура нахождения граничных условий иа и Nra по известным силовым, кинематическим или смешанным граничным условиям — procedure граничные условия. Процедура по известным граничным условиям вычисляет неизвестные граничные условия Nra или иа по формулам" § 4 гл. 1 по шести вариантам. В исходной информации задается условное число, которое определяет вид граничных условий. 6. Процедура вычисления запасов прочности и веса — procedure ЗАПАС. Процедура по формулам главы 4 вычисляет запасы местной статической прочности &Мг, &М6' запа по разрушающим оборотам k/ц, и запас по цилиндрическому разрушению kbz- В той же процедуре выбираются минимальные запасы и соответ-

интегральная функция, т — значение интеграла. Интегрирование ведется по формуле Симпсона по семи расчетным сечениям. 3. Процедура определения фундаментальных функций уравнения растяжения — procedure FUNCT1 (х, у); здесь х — функция F1{, у — решение Ф1(-. 4. Процедура определения фундаментальных функций уравнения изгиба — procedure FUNCT2 (х, у); здесь х — функция F2i у — решение Фае-. 5. Процедура интерполяции свойств материала по фактическим температурам диска — procedure INTERPOLATION. 6. Процедура нахождения нового секущего модуля при упругопластическом решении — procedure PLASTICITY. 7. Процедура печати исходной информации procedure INFORMATION. Процедура печатает исходную информацию о диске (геометрические параметры, свойства материала, нагрузки) в виде таблицы. 8. Процедура нахождения граничных условий уравнения растяжения иа и Nra по известным силовым, кинематическим или смешанным граничным условиям — procedure STRETCHING. 9. Процедура нахождения граничных условий уравнения изгиба da, Msa по известным силовым, кинематическим или смешанным граничным условиям -— procedure BEND. 10. Процедура печати результатов расчета — procedure RESULTS. Процедура печатает результаты упругого и упругопластического расчета в виде таблицы (см. текст программы). 11. Процедура вычисления запасов прочности — procedure ЗАПАС. Процедура вычисляет запасы местной статической прочности k^r &мв> запас по разрушающим оборотам kb и запас по цилиндрическому разрушению kb2- Далее выбираются минимальные запасы и соответствующий им радиус. 12. Процедура вычисления массы диска — REAL PROCEDURE F2 (х), 13. Процедура вычисления Л/0 и MQ — procedure DISTRICT. 14. Процедура вычисления угла подъема •& (г) в результате деформации—procedure TETA. 15. Процедура расчета диска на растяжение и изгиб с учетом восстанавливающего эффекта в упругопластической области — procedure DISK- Эта управляющая процедура с использованием всех описанных выше процедур производит расчет диска и печатает результаты расчета в упругой и упругопластической области. Составление исходной информации. Для расчета необходимо составить исходную информацию в следующем порядке: 1) Л/1 — число разбиений диска по радиусу: 2) N2 — число разбиений по кривой деформирования (включая ноль). 3) Л/3 — число табличных температур, при которых заданы характеристики материала (а, сгдл, Е), кривые деформирования. 4) Л/4 — число радиусов, по которым приложены сосредоточенное силы, включая радиус опоры диска. 5) ГАММА (г/см3) — плотность материала (у); 6) N — частота вращения, об/мин («). 7) DELTA — задаваемая точность расчета. 8) L — условное число, определяющее вид граничных условий (L = 1-нЮ): L = 1 —Nsa, Nsb> Msa, Msb; L = 2 — ua, uh> Msa, MSb\ L =3 — ua, ub, -flvj. $b\ L = 4 — Nsa> A/s&> #«, •&&; L = 5 — Nsa> ub, Msa, ®ь', L, = 6 — ua, NSb, fta, М&ь; L — 1 — — Nst>, Ms!>; L = 8 — ub, Msb; L = 9 — ub, -&b; L = 10 — Nsb, •&*,; L = l-f-6— диск с отверстием, L= 7-=-10—диск без отверстия. 9) иа, см. 10) иь, см. 11) ®а, рад. 12) $й, рад. 13) NSa> кгс/см. 14) Wsb> кгс/см. 15) MSa кгс. 16) М$ь, кгс. Неизвестные граничные условия задаются нулями. 17) /?1, см — радиус опоры диска. 18) R, см — радиусы расчетных сечений в количестве ЛП. 19) Н, см — толщины диска по соответствующим радиусам. 20) LF, см — расстояния от произвольно выбранной плоскости до правой стороны диска, 21) О [1 : : N1, 1 : 4] —двухмерный массив нагрузок на диск; О [1 : ./VI, 1] — температура диска в расчетных радиусах на левой стороне диска; G [1 : А/1,2] — то же самое на правой стороне диска; О [1 : N1,3] — qzr — распределенная осевая нагрузка (кгс/см2) по расчетным радиусам; G [1 : ЛПД] — ms — моменты, распределенные по радиусу. 22) Е — массив деформаций по кривой деформирования, включая ноль в количестве, равном Л/2. 23) El [1 : Л/3, 1:4] — массив табличных свойств материала, где ?1 [1 : Л/3,1 ] —табличные температуры, при которых известны свойства материала, причем табличные температуры должны охватывать диапазон температур, до которых нагрет диск; ?1 [1 : : Л/3,2 ] — с* СУград) — коэффициент линейного расширения при заданных табличных температурах; ?1 [1 : Л/3,3] —0дл (кгс/сма) при тех же температурах; ?1 [1 : Л/3,4]—? (кгс/см2)—модуль Юнга при тех же температурах. 24) ?2 [3 : Л/2, 1 : Л/3] — массив напряжений по кривым деформирования для