16+
Регистрация
РУС ENG
Расширенный поиск

Библиотека

(Всего документов: 1413)

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



РАО "ЕЭС России"



Открытое акционерное общество
"ВСЕРОССИЙСКИЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ
ГИДРОТЕХНИКИ им. Б.Е. ВЕДЕНЕЕВА"



ПОСОБИЕ ПО РАСЧЕТУ ДИНАМИЧЕСКОГО НАПРЯЖЕННОГО СОСТОЯНИЯ МЕТАЛЛИЧЕСКИХ ОБЛИЦОВОК КАМЕР РАБОЧИХ КОЛЕС ПОВОРОТНО-ЛОПАСТНЫХ ГИДРОТУРБИН С ОЦЕНКОЙ ИХ УСТАЛОСТНОЙ ПРОЧНОСТИ



П 77-2000
ВНИИГ

УДК 621.224

Срок введения в действие IV кв 2000 г

Внесены ОАО "ВНИИГ им. Б. Е. Веденеева"



Утверждены РАО "ЕЭС России" письмо № 02-1-03-4/632 от 03.07.98



Вводится впервые





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

Применение разработанных методик и программного обеспечения позволяет обосновать выбор плановых размеров и толщин облицовок камер рабочих колес (КРК) с помощью современных подходов к расчету нагрузок и динамического напряженного состояния облицовок КРК при обеспечении заданного уровня их надежности.

В приложениях к Пособию даны краткая инструкция для работы с пакетом программ по расчету гидродинамических нагрузок, динамическому напряженному состоянию облицовок КРК и прогнозу их усталостной прочности, а также пример расчета этих величин для конкретных ГЭС.

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

Пособие предназначено для работников проектных, конструкторских и научно-исследовательских организаций.

ПРЕДИСЛОВИЕ

"Пособие по расчету динамического напряженного состояния металлических облицовок камер рабочих колес поворотно-лопастных гидротурбин с оценкой их усталостной прочности" содержит материалы, обоснованные современным опытом строительства, проектирования и исследований вопросов, связанных с различными аспектами условий работы металлических облицовок (динамическое воздействие потока, адгезия к бетону и другие разрушающие факторы), и помогает в практике проектирования и реконструкции пропеллерных и поворотно-лопастных гидротурбин (ПЛ-гидротурбин) при оценке надежности облицовок камер рабочих колес (КРК). Пособие включает современное программное обеспечение для определения изменяющихся во времени гидродинамических нагрузок, динамического напряженного состояния и прогноза усталостной прочности облицовок КРК при различных схемах их работы.

Металлические облицовки относятся к объектам экономической категории, стоимость их составляет значительную часть в общей стоимости гидросооружений. Технические параметры таких объектов, в том числе надежность, подлежат, по возможности, наиболее точной и корректной оценке на уровне современного состояния науки в данной области с учетом практики строительства и эксплуатации КРК ПЛ-гидротурбин.

Настоящее Пособие позволяет обосновать выбор плановых размеров и толщин облицовок КРК с помощью современных подходов к расчету нагрузок и динамического напряженного состояния облицовок КРК при обеспечении заданного уровня их надежности.

Пособие разработано во ВНИИГ им Б.Е. Веденеева доктором физ.-мат. наук В.И.Климовичем и канд. техн. наук С.М. Левиной при участии ведущего инженера Н.Н. Короткова

1. ОБЩИЕ ПОЛОЖЕНИЯ

Анализ данных обследований разрушений облицовок КРК на действующих ГЭС позволяет выделить в качестве основных причин, приводящих к их непроектной работе, наличие пустот за КРК или отсутствие сцепления облицовок КРК со штрабным бетоном. Нарушение контакта с бетоном может происходить вследствие плохой проработки бетона, его усадки и выноса по фильтрационным ходам и т.п. Учет изменения проектных условий работы облицовок КРК должен быть произведен и в расчетах путем рассмотрения динамического напряженного состояния КРК с учетом эффекта гидроупругости и наличия подплитных областей. Кроме того, для расчета динамического напряженного состояния облицовок КРК необходимо также определить расчетным путем гидродинамические нагрузки, действующие на стенки КРК при различных режимах работы гидроагрегата, которые являются наиболее характерными для данной ГЭС. На основе указанных выше расчетов должен быть произведен прогноз усталостной прочности облицовки КРК и дана оценка его долговечности [1].

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

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

Полная нагрузка на облицовку КРК равна сумме нагрузок, действующих как со стороны проточной части, так и со стороны штрабного бетона или подплитной полости. При этом вследствие вращения РК гидродинамическая составляющая нагрузки со стороны проточной части переменна во времени.

Расчеты гидродинамической составляющей нагрузки на облицовку, местоположение и размеры которой заданы, производятся для заданных режимов работы гидроагрегата с помощью современных методов расчета гидродинамики проточной части [2,3]. Зазоры между камерой и лопастью РК полагаются малыми

При расчете гидродинамических нагрузок на облицовку со стороны камеры принимается гипотеза об идентичности спектра нагрузки на жесткую и упругую конструкции.

Режим работы гидроагрегата характеризуется расходом воды через гидротурбину, числом оборотов вращения РК, напором, высотой отсасывания, открытием направляющего аппарата (НА), углом установки лопастей РК.

Для гидродинамического расчета течений в проточной части гидротурбины требуется исходная информация (чертежи) о геометрии профиля НА, лопасти РК, нижнего и верхнего ободов проточной части.

Для расчета гидродинамической нагрузки на облицовку КРК необходимы чертежи КРК ПЛ-гидротурбин.

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

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

Облицовка рассчитывается как упругая пластина (схематизация облицовки КРК в виде плоских пластин для стандартных размеров элемента справедлива при диаметрах рабочего колеса больше 3 м), закрепленная на прямоугольном жестком контуре, граничащая с одной стороны с потоком в проточной части, а с другой - с подплитной полостью, наполненной водой.

При динамическом расчете учитывается гидроупругое взаимодействие пластины с окружающей ее водной средой.

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

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

При расчете вероятности отказа элемента стальной облицовки КРК принимается во внимание два критерия отказа:

превышение уровнем статических напряжений в элементе предела прочности Sст > RН (где Sст- расчетные напряжения от статических нагрузок; RН - предел прочности);

нарушение условия усталостной прочности с учетом коррозии Sуст > RN-1 где Sуст - расчетные усталостные напряжения с учетом коррозии; RN-1 - условный предел коррозионно-усталостной прочности.

Расчет по первому критерию проводится при воздействии на облицовку статической нагрузки, по второму критерию - на основе динамического прочностного расчета.

Динамический прочностной расчет состоит из двух этапов: расчет спектра гидродинамических нагрузок и спектра напряжений, соответствующих этим нагрузкам.

Расчет долговечности стальной облицовки для данной совокупности режимов предлагается проводить по двум теориям: детерминистической и вероятностной.

При работе агрегата на разных режимах для определения ресурса стальной облицовки КРК используется линейная модель суммирования повреждений.

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

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

Область применения

Настоящее Пособие распространяется на расчеты металлических облицовок (обечаек) камер рабочих колес поворотно-лопастных гидротурбин, проводимые как на стадиях проектирования и реконструкции КРК, так и на стадии ремонтно-восстановительных работ.

Пособие предназначено для работников проектных, конструкторских и научно-исследовательских организаций.

Нормативные ссылки

В настоящем Пособии приведены ссылки на следующие нормативные документы:

ГОСТ 19185-73 "Гидротехника. Основные понятия. Термины и определения";

ГОСТ 27.002-89 "Надежность в технике. Основные понятия и определения";

РТМ 108.243.101-83 "Турбины гидравлические поворотно-лопастные. Расчеты на усталостную прочность".

Терминология

Термины и определения, использованные в настоящем Пособии, приняты в трактовке ГОСТ 19185-73 , ГОСТ 27.002-89 и РТМ 108.243.101-83.

2. РАСЧЕТ ГИДРОДИНАМИЧЕСКОЙ НАГРУЗКИ НА ОБЛИЦОВКУ КРК СО СТОРОНЫ ПРОТОЧНОЙ ЧАСТИ ГИДРОТУРБИНЫ



2.1. Общие положения

Для расчета гидродинамической нагрузки на облицовку КРК со стороны проточной части гидротурбины необходимо для заданных режимов работы гидроагрегата рассчитать течение во всей проточной части, включая области НА, РК и свободных вихревых течений до НА, между НА и РК, за РК. Расчетно-теоретическое изучение течения жидкости затрудняется вследствие сложной геометрической формы проточной части, пространственного характера течения, необходимости учета влияния вязкости, наличие которой проявляется в образовании сильно завихренных областей пограничных слоев на стенках проточной части, а также вихревых следов за лопатками НА и за лопастями РК. В настоящее время существуют модели различной сложности для исследования потока в гидромашинах. При этом в основном используются двухмерные, квазитрехмерные и трехмерные модели. В настоящем Пособии используется квазитрехмерный подход к расчету параметров течения в проточной части гидротурбины. Целесообразность применения квазитрехмерной модели движения обуславливается возможностью учесть вихревой характер течения в проточных частях гидроагрегатов, а также взаимное влияние лопастных систем друг на друга. Практическая эффективность квазитрехмерных подходов объясняется удовлетворительным согласованием результатов квазитрехмерных расчетов с экспериментальными данными для проточных частей гидроагрегатов различного типа в широком диапазоне изменения их параметров.

2.2. Модель течения жидкости в проточной части гидротурбины. Расчет гидродинамической нагрузки на КРК

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

В областях лопастных систем уравнения квазитрехмерной модели течения имеют вид [3]:

Здесь H1, H2, R - коэффициенты Ламе, соответствующие криволинейным координатам q1, q2,φ (рис.2.1); V1, V2, Vu- компоненты абсолютной скорости , полученной путем осреднения скорости по окружной координате φ в пределах межлопастного канала; s, n - естественная система координат в меридиональной плоскости, координатными линиями которой являются линии тока осредненного течения; I - интеграл Бернулли в относительном движении;

φ - Ф(q1, q2) = 0 - уравнение средней поверхности лопасти.

Отметим, что в случае густых лопастных систем (ε = L/Т > 1) уравнения, описывающие квазитрехмерное течение, могут быть упрощены за счет отбрасывания слагаемых малого порядка А1, А2, А3 и учета связи осредненных параметров

где ω - угловая скорость вращения лопастной системы.



Рис. 2.1. Меридиональная проекция проточной части ПЛ-гидротурбины

В областях свободного вихревого течения уравнения квазитрехмерной модели течения имеют вид

где Е - интеграл Бернулли в абсолютном движении.

Отметим, что слагаемое А3 в формуле (2.6) в случае густых лопастных систем может быть отброшено, так как имеет малый порядок.

Численная схема решения уравнений (2.1) - (2.7) строится следующим образом. Уравнение (2.2) (см. также (2.6)), записанное для функции тока, аппроксимируется конечными разностями. Полученная нелинейная система разностных уравнений с учетом граничных условий решается методом установления. Решение задачи обтекания решеток профилей на осесимметричных поверхностях тока производится на основе программы, разработанной Б.С. Раухманом [4], в которую был добавлен расчет характеристик течения в канале между соседними профилями.

На основе полученного в результате расчета по уравнениям (2.1)-(2.7) поля скоростей из интеграла Бернулли определяется распределение давления в РК

где Евых = gPат, Рат - атмосферное давление (м вод.ст.); Z1 = Нs - Z, Нs - высота отсасывания; Z- координата, отсчитываемая по оси гидротурбины от оси поворота лопастей РК в направлении течения; RVu,вых - момент окружной скорости на выходе из РК; φ - окружная координата в относительном движении.

Из последнего соотношения в координатах относительного движения получаем следующее выражение для давления Р (м вод.ст.) в РК (P = p/(ρg))

В случае возникновения областей с давлением Р < Рп (что соответствует зоне кавитации, Рп - давление парообразования) в этих областях полагалось Р = Рп.

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

где t - время; sp = sp (q1, q2 - координата, отсчитываемая вдоль периферийного обода; φ0, q1, q2 - неподвижная система координат (φ0 = φ + ωt).

Вследствие осевой симметрии (между соседними лопастями РК распределение давления в канале одинаково) гидродинамическая нагрузка со стороны проточной части Р0 является периодической функцией по t с периодом, равным T = 2p/(N1ω), где N1 - число лопастей РК.

2.3. Расчет спектральных характеристик гидродинамической нагрузки на элемент облицовки КРК

После расчета параметров потока в проточной части ПЛ-гидротурбины и гидродинамической нагрузки на КРК определяется средняя по площади нагрузка Рs(t) на облицовку КРК площадью S

Спектральные характеристики гидродинамической нагрузки на облицовку Рs(t) находятся на основе разложения Рs(t) в ряд Фурье

При этом спектральная плотность gpp(Ω) динамической нагрузки Ps(t) имеет вид

Отметим, что спектральная плотность gpp(Ω) имеет вид (2.13) при ее вычислении с помощью соотношений Хинчина-Винера [5]

и учете разложения Ps(t) в виде (2.12).

Характеристики спектральной плотности гидродинамической нагрузки (коэффициенты c2k и Ωk), рассчитанные с помощью зависимостей (2.12), (2.13), автоматизированно передаются в блок расчета динамического напряженного состояния облицовки (см. раздел 3).

2.4 Определение режимных параметров ПЛ-гидротурбины, необходимых для расчета гидродинамической нагрузки

Определение режимных параметров для расчета гидродинамики проточной части и гидродинамической нагрузки на КРК должно осуществляться на основе эксплуатационной и универсальной характеристик ПЛ-гидротурбины. При этом, для заданного режима, т.е. при известных напоре Н, мощности N и оборотах РК п, необходимо установить открытие НА а0, угол установки лопастей РК φл, расход через турбину Q. Кроме того, для расчета распределения давления на КРК необходимо определить высоту отсасывания Hs.

Режимные параметры a0, φл, Q и Hs при заданных Н, N и n находятся следующим образом. По эксплуатационной характеристике гидроагрегата по заданным H и N вычисляются коэффициент полезного действия (кпд) η и высота отсасывания Нs. С помощью соотношения

N = ρgQHη / 1000 (кВт)     (2.15)

определяется расход жидкости Q (м3/с) через турбину.

Далее рассчитываются приведенный расход Q‘1 и приведенные обороты n‘1:

где D1 - диаметр РК.

С помощью универсальной характеристики по найденным значениям Q‘1 и n‘1 и с учетом поправки на масштаб модели находятся открытие НА модельной турбины a‘0 и угол установки лопастей РК φл. Окончательное открытие НА для условий натурной гидротурбины а0 пересчитывается следующим образом

a0 = a‘D1 / Dм,     (2.17)

где Dм - диаметр модели турбины.

Для выполнения гидродинамического расчета необходимо также задать энергию потока Е0 на входе в проточную часть и закрутку потока RVu,вх, формируемую спиральной камерой. Эти значения находятся из соотношений:

где b0 - высота НА; δсп - угол спирали.

Угол спирали δсп выбирается на основе данных [6, 7] в зависимости от угла охвата спиральной камеры, который известен для каждой гидротурбины.

2.5. Подготовка исходных данных для гидродинамического расчета проточной части ПЛ-гидротурбины

На основе полученных значений режимных параметров а0 и φл (см. п. 2.3) для натурных условий необходимо решить две вспомогательные геометрические задачи по определению требуемых для гидродинамического расчета геометрических характеристик лопастных систем: коэффициента стеснения τ(R,Z) и угловой координаты средней линии лопасти (лопатки) Ф(R,Z) (см. п. 2.2).

Для указанных геометрических задач требуется иметь теоретический чертеж лопасти РК и чертеж лопатки НА.

Для расчета функций Ф и τ в НА по известному открытию задаются:

число лопаток НА - N0,

координата R0N (м) цилиндрического радиуса оси поворота лопаток HA(R0N = D0 / 2);

высота НА - b0 (м);

координата Z(м) нижнего кольца НА (Z = 0 на оси поворота лопастей РК);

массив координат {xi, yi} (м) верхней и нижней стороны профиля (рис. 2.2);

координаты x0, y0 (м) оси поворота лопатки НА;

координаты xв, yв (м) центра вписанной окружности и ее радиус rв на входной кромке НА ;

координаты xE, yE (м) центра вписанной окружности и ее радиус rE (м) на выходной кромке НА (условно принимается, что радиус rE равен половине толщины выходной кромки НА);

открытие НА - а0 (мм).

Для расчета функций Ф и τ в РК по известному углу установки лопастей φл задаются:

число лопастей РК N1;

число цилиндрических сечений R = const, в которых задаются сечения лопасти РК;

массив цилиндрических радиусов R (м) этих сечений ;

расстояния до верхней b1 (м) и нижней баз b2 (м) от оси поворота лопасти в развертке цилиндрического сечения на плоскость (рис.2.3);

угол разворота лопастей φл (град);

массивы координат {xвi, уi} (м) точек верхней стороны профилей (от верхней базы) в сечениях R = const;

массивы координат {xHi, уi} (м) точек нижней стороны профилей (от нижней базы) в сечениях R = const;

массив координат хc, уc (м) центров вписанных окружностей на входных кромках профилей и их радиус rc (м) в сечениях R = const (хc отсчитывается от верхней базы);

массив координат хc1, уc1 (м) центров вписанных окружностей на выходных кромках профилей и их радиус rc1 (м) в сечениях R = const (хc1 отсчитывается от верхней базы, условно принимается, что rc1 равен половине толщины выходной кромки лопасти);

массив координата y0 (м) положения оси вращения лопасти в сечениях R = const.

Для расчета гидродинамики проточной части требуется также задание геометрии нижнего (CD) и верхнего (АВ) ободов проточной части (см. рис.2.1). При этом координаты R, Z линий АВ и CD должны выбираться в одинаковом числе точек и точки рекомендуется располагать приблизительно равномерно вдоль линий АВ и CD. Линии DA и ВС задаются соответственно как R = const и Z = const. Линия DA выбирается между статором и НА, линия ВС - за РК в конусе отсасывающей трубы.

Местоположение и размеры облицовки КРК (длина облицовки Ds по оси турбины, ширина облицовки DL в окружном направлении, координата Z положения нижнего края облицовки ZHO) находят с помощью чертежа КРК.

Помимо геометрических данных проточной части ПЛ-гидротурбины, в соответствии с п. 2.4, находят также обороты РК n (об/мин), расход Q (м3/сек), энергию потока Е0 и закрутку потока RVu,вх на входе в проточную часть, атмосферное давление Рат (м), высоту отсасывания Hs (м).

3. РАСЧЕТ ДИНАМИЧЕСКОГО НАПРЯЖЕННОГО СОСТОЯНИЯ СТАЛЬНОЙ ОБЛИЦОВКИ



3.1. Выбор расчетной схемы

В качестве расчетных схем для облицовок камер рабочих колес типичными, наиболее распространенными, принимаются две основные схемы (рис.3.1):

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

облицовка в пределах рассчитываемой панели, ограниченной ребрами жесткости, не омоноличена и имеет подплитную полость глубиной а (рис.3.1, б, в - а > 0; рис. 3.1, г - а = 0, что соответствует термину "отслоение"). При этом выделяются два частных случая: облицовка не имеет отверстий в панели, соединяющих воду в подплитной полости с транзитным потоком (рис.3.1, б); облицовка имеет отверстие в панели, соединяющее воду в подплитной полости с транзитным потоком (рис. 3.1, в, г).

3.2. Нагрузка на элемент облицовки

Полная нагрузка Рвн на облицовку КРК определяется соотношением

Pвн = Ps(t) + Pф,     (3.1)

где Ps(t) - давление в камере со стороны проточной части, состоящее из осредненной Р и пульсационной Рг составляющих; Рф - фильтрационное давление.

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



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

3.3. Основные предпосылки расчета на прочность элемента облицовки

Расчет на прочность элемента облицовки состоит из статического и динамического расчетов.

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

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

Pф = kpf* (H - Hs + Zp),     (3.2)

где Н - напор; Нs - высота отсасывания; Zp - координата центра облицовки, отсчитываемая от оси поворота лопасти РК; kpf - коэффициент, меняющийся от 0 до 1 (значение 1 соответствует наихудшему случаю, когда давление в подплитной полости равно напору).

Анализ влияния фильтрационного давления и остаточных температурных напряжений на упругую устойчивость и вынужденные колебания облицовок, проведенный на основе решения задачи гидроупругости в нелинейной постановке, показал, что достижение границ статической и динамической устойчивости для применяемого в КРК диапазона толщин облицовок (точнее диапазона их жесткости) возможно и при реальных скоростях транзитного потока в том случае, если уровень напряжений в срединной поверхности облицовки Nx, Ny ~ (1 - 3)107 Н/м2.

При проведении расчетов динамического напряженного состояния и усталостной прочности фрагментов облицовки КРК и сопрягающего пояса следует придерживаться геометрической схематизации: расчетным элементом облицовки является фрагмент, ограниченный вертикальными и горизонтальными ребрами жесткости КРК, размеры которого отсчитываются по образующей цилиндрической конструкции КРК (b) и по окружности сечения КРК перпендикулярно ее оси (L = 2 l). Фрагмент облицовки представляется как плоская, упругая панель постоянной толщины, радиусы кривизны R реальных фрагментов облицовки принимаются равными ∞.

Следует иметь в виду, что это предположение может вносить существенные погрешности в расчет спектра собственных частот панелей при расчете КРК гидротурбин с малыми (< 3 м) радиусами РК. В этом случае наличие цилиндричности панели приводит к уплотнению спектра собственных частот, подавлению одних и возбуждению других форм колебаний.

Учитывается наличие подплитной полости глубиной а и возможность хотя бы одного незаваренного инъекционного отверстия.

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

При динамическом расчете фрагмента облицовки необходим учет гидроупругих эффектов. С ростом скорости транзитного потока собственная частота панели облицовки падает. Значение скорости Ud, при которой она становится равной нулю, соответствует потере пластинкой статической устойчивости (началу дивергенции - выпучиванию). В первом приближении оценку минимального значения этой скорости можно проводить по графику (рис.3.2)



Рис. 3.2. Зависимость безразмерной скорости дивергенции Ud от безразмерного комплекса μ = ρ02l / ρmh (ρ0 - плотность воды; ρm - плотность металла; h - толщина; 2l - длина элемента облицовки).
___ - Ud/2lω01; --- - Uср/2lω01 = 30/(0,5ω01)

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

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

р = рн + рп.     (3.3)

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

Результаты экспериментальных исследований подтверждают априорные предположения о том, что податливость границы потока должна приводить к перестройке спектра пульсаций давления на ее поверхности за счет добавления индуцированных пульсаций рн. Это проявляется в появлении в спектре потока дополнительных, ярко выраженных пиков на частотах, близких к первым собственным частотам ω01 колебаний панелей в жидкости, причем стандарты пульсаций на этих частотах (для натурных облицовок) становятся сравнимы со стандартами на ведущих частотах основного потока. Кроме этого значительно увеличивается стандарт пульсаций давления в широкой полосе частот. Происходит существенное возрастание коэффициентов взаимной пространственной корреляции между пульсациями давления на всей упругой поверхности, что подтверждает справедливость используемого при динамическом расчете допущения возможности осреднения пульсационной нагрузки Рг(t) по площади фрагмента.

Элемент облицовки испытывает динамические нагрузки со стороны проточной части и пульсаций давления, проникающих в подплитную полость (при наличии непроваренного отверстия). Приближенно определить частоты последних можно, например, по формуле Росситера, которая предназначена для оценки резонансной частоты малой полости, соответствующей отверстию в элементе облицовки: , где m - номер моды колебаний пластинки; d0 - длина отверстия вдоль потока; k = C1/U≈0,5 - безразмерное значение скорости распространения возмущений; M0 - число Маха в заторможенном потоке.

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

Методика динамического расчета стальных облицовок КРК исходит из предположения отсутствия напряжений в срединной поверхности (удовлетворяется критерий жесткости - прогиб меньше 1/5 толщины) и что конструкция работает в области динамической устойчивости.

Согласно многочисленным расчетам максимальный прогиб элемента облицовки КРК колеблется в пределах 2 - 5 мм, поэтому при глубине подплитной полости > 1 см колебания конструкции не будут испытывать помех со стороны бетона, при меньших глубинах амплитуда колебаний будет ограничена со стороны бетона, но величина ее при вогнутости внутрь камеры не уменьшится. Таким образом, при глубине водной прослойки 1-2 мм стальная облицовка не работает совместно с бетоном. Глубина подплитной полости >1 см мало влияет на колебания и, соответственно, на напряжения, возникающие в теле элемента облицовки. Расчет долговечности элемента облицовки на большую глубину подплитной полости может иметь чисто проверочный характер.

При определении динамических напряжений в теле облицовки предполагается, что конструкция достаточно жесткая и колебания ее можно считать малыми, а саму систему линейной. Тогда напряжения σ будут выражаться через нагрузки в виде σ = L(Pвн), где L - линейный оператор, а спектр напряжений будет равен спектру внешней нагрузки Sрг, умноженному на квадрат передаточной функции системы

Передаточная функция Пσ(iω), представляющая собой ничто иное как реакцию упругой системы на единичную нагрузку вида Pвн = eiωt, находится теоретически путем решения задачи гидроупругости.

3.4. Постановка задачи гидроупругости

Задача решается при следующих предпосылках:

рассматривается потенциальный поток жидкости, сжимаемость которой определяется степенью равномерно распределенной аэрации потока;

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

краевые условия элемента облицовки принимаются либо в виде шарнирного опирания, либо в виде жесткого защемления;

рассматривается напорное течение жидкости;

подплитная полость глубиной а заполнена водой (возможна ее гидравлическая связь с транзитным потоком);

колебания элемента считаются малыми.

Математическая формулировка задачи:

Уравнение колебания пластины

где w - динамический прогиб пластины; Рвн и р определяются, соответственно, по формулам (3.1) и (3.3);

γ - коэффициент неупругого сопротивления; - цилиндрическая жесткость пластины; ρm - плотность материала пластины; Е - модуль упругости; ν - коэффициент Пуассона; h - толщина; l - половина длины пластины; t - время.

Граничные условия:

а) шарнирное опирание w = 0, , при х = -l и х = l;

(3.6)

б) жесткая заделка w = 0, , при х = -l и х = l.

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

Уравнение для потенциала скоростей, индуцированных колеблющимся элементом в транзитном потоке φн:

Давление, индуцированное колеблющимся элементом, определяется с помощью уравнения Лагранжа-Коши

Аналогичная задача ставится для определения реакции жидкости на колебания элемента облицовки со стороны подплитной полости, заполненной водой. Давление рn определяется через потенциал скоростей возмущений φn

который находится из уравнений гидродинамики

Граничные условия:

Решение поставленной задачи осуществляется следующим образом. Определяется реакция системы на единичную нагрузку вида Pвн = eiωt. Прогиб представляется в виде разложения в ряд по формам колебаний в воздухе

Если к уравнению (3.4) применить метод Галеркина получается система алгебраических уравнений для определения коэффициентов в разложении. Найдя коэффициенты Аk, определяется прогиб от единичной нагрузки (т.е. передаточная функция по прогибу),которая будет функцией от целого ряда параметров:

В выражение (3.14) входят безразмерные комплексы, составленные из ряда величин, которые при расчете считаются известными:

ρm - плотность; Е - модуль упругости; ν - коэффициент Пуассона материала облицовки (табличные данные);

h - толщина облицовки;

l - половина длины элемента облицовки;

а - глубина подплитной полости;

b - расстояние от первой по потоку кромки элемента облицовки до отверстия в нем;

d - диаметр отверстия в элементе облицовки, соединяющего транзитный поток и подплитную полость;

ω01 - первая собственная частота элемента облицовки в воздухе;

Нг - глубина потока (расстояние между стенкой камеры и втулкой рабочего колеса);

U0 - средняя скорость потока над рассматриваемым элементом облицовки;

С и С0n - скорость звука в аэрированном потоке над и под облицовкой;

ψн и ψп - коэффициент затухания, определяющий диссипацию упругих волн.

Значения С0 и ψ зависят от степени аэрации потока α [8] и определяются по таблице.

Зависимости величин скорости звука и коэффициента затухания упругих волн в воде от коэффициента аэрации

α 0 0,01 0,03 0,05 0,12
С0, м/сек 1200 440 160 155 111
ψ, сек/м 0,001 0,001 0,001 0,004 0,012

Соответственно передаточная функция по напряжению определяется

Следует отметить вполне удовлетворительное совпадение результатов расчета напряженного состояния панели облицовки по вышеизложенной методике с экспериментальными данными, полученными при моделировании облицовок КРК. Некоторые несоответствия в экспериментальных и расчетных значениях спектральных плотностей пульсаций напряжений в области ω>2ωвед и ω≈0,1ωведвед - ведущая частота процесса) следует отнести как за счет погрешностей экспериментального моделирования, так и за счет допущений, принимаемых в расчетных схемах.

Необходимые для дальнейших расчетов математическое ожидание и стандарт напряжений в теле облицовки находятся по спектральной плотности напряжений по формулам:

4. ОЦЕНКА НАДЕЖНОСТИ СТАЛЬНЫХ ОБЛИЦОВОК КРК ПЛ-ГИДРОТУРБИН



4.1 Оценка вероятности безотказной работы элемента облицовки

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

Рассчитывается вероятность отказа (или вероятность безотказной работы) элемента стальной облицовки КРК. При этом принимаются во внимание два критерия отказа:

превышение предела прочности уровнем статических напряжений в элементе

σст > σвр;     (4.1)

нарушение условия усталостной прочности с учетом коррозии

σуст > σN-1,     (4.2)

где σN-1 - условный предел коррозионно-усталостной прочности; σуст - расчетные усталостные напряжения с учетом коррозии.

Оценка вероятности отказов ведется по методике, основанной на работах А.Р. Ржаницына [9].

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

S = R - σ > 0.     (4.3)

В расчетах на прочность: σ - наибольшее значение усилия или напряжения в элементе, выраженное через внешнюю нагрузку; R - несущая способность, выраженная в тех же единицах, что и σ, и отвечающая предельному состоянию конструкции по прочности; S - резерв прочности. В общем случае усилия и несущая способность являются случайными функциями времени, но в рассматриваемой постановке R и σ считаются случайными величинами с заданными законами распределения f(R) и f(σ).

Отказ элемента облицовки происходит во всех случаях, когда не выполняется условие (4.3). Под отказом в данном случае понимается достижение элементом предельного состояния по прочности.

При любых законах распределения R и σ

где S,R, σ - математические ожидания, а S‘, R‘, σ‘ - стандарты, соответственно, резерва прочности, несущей способности и напряжения в элементе.

Под характеристикой безопасности понимается число стандартов S, укладывающееся в интервале 0 ≤ S ≤ S

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

где k(R,σ) - корреляционный момент случайных величин R и σ.

В практических расчетах для характеристик прочности можно принимать нормальное распределение. Закон распределения напряжений в элементах облицовки КРК тоже считается нормальным.

В предположении, что R и σ подчинены нормальным законам распределения, проектная вероятность безотказной работы элемента конструкции выражается

Р = P (R > σ) = 0,5 + Ф(γ),     (4.7)

где - интеграл вероятности Гаусса.

Для оценки надежности элемента стальной облицовки необходимо знать математические ожидания и стандарты соответствующих характеристик прочностных свойств конструкции R и внешних воздействий σ(3.16).

Приводимые в ГОСТ нормативные значения предельных напряжений являются минимально допустимыми и их можно отнести к 98 - 99% вероятности неразрушения. Поэтому необходимо ввести коэффициент для перехода к математическому ожиданию предельных напряжений, т.е. к 50% вероятности неразрушения. Так как предел текучести σm имеет обеспеченность 98%, то предельное значение функции распределения равно Fпр = 1 - 0,98 = 0,02. Отсюда

где σm и σ‘m - математическое ожидание и стандарт предела текучести. Равенство выполнятся при

где - коэффициент вариации предела текучести.

В случае, если входящие в формулу (4.3) случайные величины R и σ являются функциями других случайных величин, их вероятностные характеристики можно определять следующим образом. Пусть Х = φ(Y1, ..., Yn) - известная дифференцируемая функция случайных аргументов Y1, ..., Yn. Из разложения функции φ в ряд Тэйлора вытекают зависимости для математического ожидания и стандарта величины X:

где - частная производная функции φ по фактору Yi, вычисленная при Yi = Y1,..., Yn = Yn. В практических расчетах некоторые из величин Yi можно считать детерминированными. В этом случае соответствующие им стандарты Y‘i = 0.

4.2 Оценка долговечности элемента облицовки

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

Вероятностная методика - методика, основанная на кумулятивной модели отказов и использующая линейную модель суммирования усталостных повреждений [10]. При этом характеристическое значение условной долговечности Tc, близкое к математическому ожиданию или наиболее вероятному значению условного ресурса (срока службы), определяется [10] по формуле

Плотность вероятности максимумов Р(σ) для стационарного гауссовского процесса σ(t) с математическим ожиданием σ выражается формулой Райса

Здесь β - коэффициент широкополостности процесса

За σmin принимается предел выносливости σ-1, а за σmax - предел текучести σm.

Коэффициент α согласно линейной теории равен 1.

Показатель долговечности Тс - случайная величина. Однако, чем больше число циклов до разрушения, тем меньше относительная изменчивость условной долговечности. Поэтому, в первом приближении, можно принять Т1 = Тс с вероятностью, равной 1.

Детерминистическая методика - широко используемая в настоящее время методика, где за долговечность элемента принимается время Т, за которое происходит снижение запаса прочности элемента nσ до его нормативного значения [n*σ] [11]:

где ;

Здесь σ-1N0 - предел выносливости при симметричном цикле нагружения при базовом числе циклов N0 = 107; (kσ)D - коэффициент, учитывающий масштабный фактор, концентрацию напряжений и влияние обработки; ψσ - коэффициент, учитывающий относительное уменьшение предельной амплитуды с ростом среднего напряжения цикла; р - степенной показатель (аналог коэффициента регрессии). По зависимости (4.14) находится число циклов колебаний N, за время которых запас прочности элемента снизится до нормативного.

Время работы (срок службы) определится как

T2 = N / Nгод,     (4.15)

где Nгод - заданное число колебаний нагрузки за год.

Величина n*σ может быть назначена, например, по рекомендациям НПО ЦКТИ как при расчете усталостной прочности лопастей ПЛ-гидротурбин, равной [n*σ] = 1,5.

При широкополосных спектрах нагрузок и напряжений для получения более надежных результатов предлагается также проводить оценку времени безотказной работы как величину, обратную сумме "элементарных мер повреждения" (Tj)-1

где для Тj с учетом зависимости (4.14) и моделирования спектра его разложением в ряд по гармоническим составляющим можно записать

В целом оценки по обеим теориям дают совпадение по граничным значениям параметров.

По изложенным выше методикам просчитывается условная долговечность для ряда возможных значений параметров (геометрических) элемента облицовки. За конечное значение долговечности принимается:

T = min (T1, T2).     (4.18)

Для аналитического описания кривой усталости рекомендуется пользоваться зависимостью вида

где базовое число циклов N0 и параметр кривой m берутся по справочникам; σ-1 - предел выносливости; N(σ) - число циклов до появления усталостного разрушения при действии переменных напряжений с амплитудой σ.

За время работы элемента, равное долговечности, вероятность его безотказной работы (надежность) близка к 1, но не равна ей, т.к. имеется разброс прочностных характеристик, воздействий и свойств элемента.

Оценка вероятности безотказной работы элемента облицовки при воздействии переменной нагрузки производится по зависимостям (4.5) или (4.6) и (4.7).

4.3. Оценка долговечности элемента облицовки с учетом реального графика работ

Турбина в процессе эксплуатации имеет разные режимы работы (пусть их будет n), каждому из которых соответствует свой нагрузочный режим на стальную облицовку. При оценке долговечности металлической облицовки возможны два подхода:

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

если можно оценить рабочие режимы турбины в течение года, то рекомендуется проводить расчет долговечности с учетом продолжительности работы турбины на каждом режиме.

Пусть в течение года соотношение между режимами определяется коэффициентами ki. Время работы на первом режиме примем за единицу, следовательно k1 = 1, тогда годовая наработка будет

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

Процесс для каждого из режимов считается стационарным.

Пусть j - номер блока, причем j = 0 соответствует начальному нагрузочному режиму, начиная с t = 0. Значение меры повреждений непосредственно после окончания j-го блока обозначим ψj. Тогда имеем [10] уравнение

ψj - ψj-1 = w(ψj-1, qj),.....(j = 1,2...N),     (4.21)

где N - номер блока, соответствующего исчерпанию ресурса.

При этом в начале первого блока мера повреждений равна нулю, а после работы N блока нагружений равна единице, т.е. ψ0 = 0; ψN = 1.

Применяется модель линейного суммирования повреждений, тогда, если N - разрушающий блок, то

Физическое время непрерывно и уравнение для меры повреждения для каждого блока

В дальнейшем принимается кумулятивная модель: считается, что повреждения, вносимые каждым циклом или блоком, не зависят от состояния конструкции (ее повреждения) на предыдущем цикле. Тогда в формулах (4.21) и (4.23) правые части зависят только от вектора нагрузки.

Решение уравнения (4.23) для j-го блока нагружения



где - ресурс при стационарном режиме нагружения с заданным значением вектора qi.

Суммарная мера повреждений за N блоков равна 1, тогда с учетом того, что ψ0 = 0, имеем

При этом Δtj = kjΔt1 и тогда

где N1 + N2 +...... + Nn = N - сумма числа блоков нагружения qj.(j = 1, 2....n),

Суммарная долговечность равна

Если известен график годовой нагрузки, то N находится сравнительно легко из зависимости (4.26).

Для нормального стационарного процесса

4.4. Оценка надежности облицовки как системы

Конструктивно облицовка представляет собой стальные листы с ребрами жесткости, к которым приварены анкеры, т.е. ее можно рассматривать как некоторую конструкцию, состоящую из большого числа однотипных элементов. Отказ одного из них ведет за собой резкое ухудшение работы и отказ рядом стоящего и т. д. Для такой системы вероятность безотказной работы определяется вероятностью безотказной работы наиболее слабого звена, т.е. если облицовка КРК состоит из m элементов и Рkj - вероятность безотказной работы j-го элемента j = 1,..., m, в k-ом режиме, а (Pkj)min - минимальная из них, то вероятность безотказной работы всей облицовки равна (Pkj)min. Долговечность всей облицовки также определяется долговечностью наиболее слабого элемента.

Приложение 1

ИНСТРУКЦИЯ ПО РАБОТЕ С ПАКЕТОМ ПРОГРАММ ПО РАСЧЕТУ ГИДРОДИНАМИЧЕСКИХ НАГРУЗОК, ДИНАМИЧЕСКОМУ НАПРЯЖЕННОМУ СОСТОЯНИЮ ОБЛИЦОВОК КРК И ПРОГНОЗУ ИХ УСТАЛОСТНОЙ ПРОЧНОСТИ



П1.1. Запуск пакета

Для запуска пакета необходимо войти в директорию, где находится пакет, и набрать в командной строке:

a.bat или можно набрать а <имя>, где имя - имя проекта. В качестве имени проекта можно взять любую комбинацию из трех символов (можно использовать латинские буквы, цифры и знаки).

Если Вы не задали имя проекта, задача предложит его задать:

Введите название ПРОЕКТА    n O

После задания имени проекта (или при изменении имени в командной строке), например "XXX", появляется меню:

Гидродинамика проточной части Долговечность облицовки Выход

(П1.1)

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

Исходные данные Расчет Проект Результаты Помощь Удалить Выход

(П1.2)

Прежде чем начать работать с этим меню необходимо внимательно прочесть описание в пункте "Помощь". Пункт "Помощь" есть почти во всех меню и подменю пакета. Подробно работа с меню (П1.2) описана в п. П1.3.

Для расчета динамического напряженного состояния облицовки и оценки ее долговечности выбирается второй пункт меню (П1.1).

После выбора расчета долговечности облицовки пакет предложит выбрать вид облицовки:

Металлическая облицовка Композитная облицовка

После выбора вида облицовки требуется уточнить вид граничных условий

Шарнирное крепление или Двухслойная конструкция

Жесткая заделка Трехслойная конструкция


После этого появляется меню:

Исходные данные Расчет Проект Результаты Помощь Удалить Выход

(П1.3)

Подробно работа с этим меню описана в п. П1.2.

П1.2. Расчет надежности и долговечности облицовки КРК

Основное меню блока расчета надежности и долговечности облицовки имеет вид (П1.3). С помощью пункта "Исходные данные" этого меню вводятся и корректируются исходные данные. Подробно работа этого пункта меню описана в п. П1.2.1. Для проведения расчета долговечности облицовки выбирается пункт "Расчет". Пункт "Проект" позволяет выбрать новый проект. Пункт "Результаты" позволяет просматривать и печатать результаты расчета. Подробно работа этого пункта меню описана в п. П1.2.2. Пункт "Удалить" позволяет удалить с диска файлы данных этого блока.

П1.2.1. Подготовка исходных данных

При выборе пункта "Исходные данные" меню (П1.3) появляется следующее подменю:

Ввод
Расчет*
Вставить
Просмотр
Перекинуть
Помощь


Прежде, чем выбрать пункт "Ввод", необходимо ознакомиться с инструкцией по вводу исходных данных, приведенной в пункте "Помощь" данного подменю. Если для данного проекта отсутствуют данные гидродинамических характеристик или их надо заново пересчитать, то необходимо выбрать пункт "Расчет". Тогда получим меню (П1.2), описанное в п. П1.3. Если данные о гидродинамических нагрузках на облицовку КРК рассчитаны в соответствующем блоке данного пакета, то с помощью пункта "Вставить" эти данные будут вставлены в файл исходных данных для расчета надежности облицовки. Пункт "Просмотр" позволяет просмотреть файл исходных данных для данного блока. При необходимости организации файла исходных данных из существующего файла другого проекта необходимо выбрать пункт "Перекинуть", тогда появляется таблица существующих проектов. Перемещением курсора надо выбрать нужный проект и нажать "Enter". В результате исходные данные выбранного проекта перебросятся в текущий проект. Если файл исходных данных проекта формируется впервые, то все исходные параметры расчета задаются из стандартного файла для тестового примера.

П1.2.2. Просмотр результатов расчета

При выборе пункта "Результаты" меню (П1.3) появится следующее подменю:

(П1.4)

Просмотр
Печать
Графики
Сводка рез.


Пункт подменю "Просмотр" позволяет просмотреть результаты расчета.

Пункт подменю "Печать" вызывает следующее подменю:

Принтер
Широкий
Узкий


После выбора типа принтера печатается выходной файл (для узкого принтера печать плотная - 17 символов на дюйм).

Пункт "Графики" подменю (П1.4) вызывает следующее подменю:

Распр.спектра давл.над облиц.
Завис.перед.функц. от частоты
Распр. спектра напр. в облиц.
Помощь
Выход


Работа с графиками описана в пункте "Помощь".

Пункт "Сводка рез." подменю (П1.4) позволяет просмотреть сводные результаты расчета.

П1.3. Расчет гидродинамической нагрузки на облицовку КРК

Основное меню блока расчета гидродинамической нагрузки имеет вид (П1.2). С помощью пункта "Исходные данные" этого меню вводятся и корректируются исходные данные. Подробно работа этого пункта меню описана в п. П1.3.1. Для проведения расчета гидродинамической нагрузки выбирается пункт "Расчет". Пункт "Проект" позволяет выбрать новый проект. Пункт "Результаты" позволяет просматривать и печатать результаты расчета. Подробно работа этого пункта меню описана в п. П1.3.2. Пункт "Удалить" позволяет удалить с диска файлы данных этого блока.

П1.3.1. Подготовка исходных данных

При выборе пункта "Исходные данные" меню (П1.2) появляется следующее подменю:

(П1.5)

Ввод
Расчет*
Просмотр
Графики
Перекинуть
Помощь


Прежде, чем выбрать пункт "Ввод", необходимо ознакомиться с инструкцией по вводу исходных данных, приведенной в пункте "Помощь" данного подменю. Если для данного проекта отсутствуют данные по направляющему аппарату или по рабочему колесу или их нужно заново пересчитать, то необходимо выбрать пункт "Расчет". Тогда получим меню, описанное в пп. П1.4 и П1.5. Пункт "Просмотр" позволяет просмотреть файл исходных данных для данного блока. Пункт "Графики" вызывает следующее подменю:

Геометрия прот. части
Энергия на вх. в пр. часть
Закрутка на вх. в пр. часть
Закрутка на вх. в р. колесо
Закрутка на вых. из р.колеса
Помощь
Выход


Работа с графиками описана в пункте "Помощь".

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

П1.3.2. Просмотр результатов расчета

При выборе пункта "Результаты" меню (П1.2) появится следующее подменю:

(П1.6)

Просмотр
Печать
Графики 1
Графики 2


Пункт подменю "Просмотр" позволяет просмотреть результаты расчета. Пункт подменю "Печать" позволяет распечатать результаты расчета. Пункт "Графики 1" подменю (П1.6) вызывает следующее подменю:

Осесимметричная задача

Мерид.линии тока в пр.части
Распр.мерид.скор.вдоль лин.тока
Распределение RCu в Н.А.
Распределение RCu в Р.К.
Распределение RCu на вх. в Н.А.
Распределение RCu на вх. в Р.К.
Распределение RCu на вых. из Р.К.
Помощь
Выход


Работа с графиками описана в пункте "Помощь".

Пункт подменю (П 1.6) "Графики 2" вызывает следующее подменю:

"Решеточная" задача

Pacnp.RCu на вх. в Р.К.
Pacnp.RCu на вых. из Р.К.
Распр.проф.потерь в Н.А.
Распр.проф.потерь в Р.К.
Распр.уд.потерь в Н.А.
Распр.уд.потерь в P.К.
Распределение RCu в Н.А.
Распределение RCu в Р.К.
Pacпp.сред.давл.над обл.
Распр.сп.ср.давл.иад обл.
Помощь
Выход


Работа с графиками описана в пункте "Помощь".

П1.4. Расчет геометрических характеристик направляющего аппарата

Основное меню блока расчета геометрических характеристик направляющего аппарата имеет вид:

Исходные данные Расчет Результаты Помощь Удалить Выход

(П1.7)

С помощью пункта "Исходные данные" этого меню вводятся и корректируются исходные данные. Подробно работа этого пункта меню описана в п. П1.4.1. Для проведения расчета геометрических характеристик направляющего аппарата выбирается пункт "Расчет". Пункт "Результаты" позволяет просматривать и печатать результаты расчета. Подробно работа этого пункта меню описана в п. П1.4.2. Пункт "Удалить" позволяет удалить с диска файлы данных этого блока.

П1.4.1. Подготовка исходных данных

При выборе пункта "Исходные данные" меню (П1.7) появляется следующее подменю:

(П1.8)

Ввод
Просмотр
Графики
Перекинуть
Помощь


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

Профиль Н. А.
Помощь
Выход


Работа с графиками описана в пункте "Помощь". При необходимости организации файла исходных данных из существующего файла другого проекта необходимо выбрать пункт "Перекинуть" меню (П1.8), тогда появляется таблица существующих проектов. Перемещением курсора надо выбрать нужный проект и нажать "Enter". В результате исходные данные выбранного проекта перебросятся в текущий проект. Если файл исходных данных проекта формируется впервые, то все исходные параметры расчета задаются из стандартного файла для тестового примера.

П1.4.2. Просмотр результатов расчета

При выборе пункта "Результаты" меню (П1.7) появится следующее подменю:

Просмотр
Печать
Графики


Пункт подменю "Просмотр" позволяет просмотреть результаты расчета. Пункт подменю "Печать" позволяет распечатать результаты расчета. Пункт подменю "Графики" вызывает следующее подменю:

Положение лопаток Н.А.
Помощь
Выход


Работа с графиками описана в пункте "Помощь".

П1.5. Расчет геометрических характеристик рабочего колеса

Основное меню блока расчета геометрических характеристик рабочего колеса имеет вид:

Исходные данные Расчет Результаты Помощь Удалить Выход

(П1.9)

С помощью пункта "Исходные данные" этого меню вводятся и корректируются исходные данные. Подробно работа этого пункта меню описана в п. П1.5.1. Для проведения расчета геометрических характеристик рабочего колеса выбирается пункт "Расчет". Пункт "Результаты" позволяет просматривать и печатать результаты расчета. Подробно работа этого пункта меню описана в п. П1.5.2. Пункт "Удалить" позволяет удалить с диска файлы данных этого блока.

П1.5.1. Подготовка исходных данных

При выборе пункта "Исходные данные" меню (П1.9) появляется следующее подменю:

(П1.10)

Ввод
Просмотр
Графики
Перекинуть
Помощь


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

Профиль Р.К.
Помощь
Выход


Работа с графиками описана в пункте "Помощь".

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

П1.5.2. Просмотр результатов расчета

При выборе пункта "Результаты" меню (П1.9) появится следующее подменю:

(П1.8)

Просмотр
Печать
Графики


Пункт подменю "Просмотр" позволяет просмотреть результаты расчета. Пункт подменю "Печать" позволяет распечатать результаты расчета. Пункт подменю "Графики" вызывает следующее подменю:

Линии Fi=Const
Помощь
Выход


Работа с графиками описана в пункте "Помощь".

П1.6. Завершение работы пакета

Для завершения работы в любом меню и перехода в меню более высокого уровня достаточно нажать клавишу Esc или выбрать пункт меню "Выход". Выход из самого верхнего меню приводит к завершению работы с пакетом.

Приложение 2

ПРИМЕР РАСЧЕТА ГИДРОДИНАМИЧЕСКИХ НАГРУЗОК, НАПРЯЖЕННОГО СОСТОЯНИЯ И ДОЛГОВЕЧНОСТИ ОБЛИЦОВОК КРК

В качестве примера рассмотрим расчет гидродинамических нагрузок, динамического напряженного состояния и долговечности облицовок третьего пояса КРК ПЛ-гидротурбины Миатлинской ГЭС (ПЛ-60-В600/1075) для режима максимального напора (Н = 59 м) и номинальной мощности (N = 113 Мвт). На рис. П2.1 - П2.6 приведены исходные данные о геометрии РК, НА, обводов, чертеж КРК, эксплуатационная и универсальная характеристики гидротурбины. Угол спирали δсп был равен 32°, n =125 об/мин., D1 = 6 м, D0 = 7,2 м, N0 =24 м, N1 = 7, b0 = 2,1 м.

С помощью эксплуатационной характеристики (рис.П2.5) определяем η = 93,5% и Нs = -1,6 м. С помощью формулы (2.11) определяем Q = 208,8 м3/с. Рассчитываем приведенный расход и обороты (см.(2.12)): Q1 = 0,765 м3/c, n1 = 97,64 об/мин. По универсальной характеристике (рис.П2.6) определяем φл = 1,6° и а0 = 26 мм. Открытие НА для натурной турбины определяется по (2.13) и равно а0 = 339,2 мм. Используя зависимость (2.14), определяем Е0 = 578,8 м22, RVи.вых = 23,3 м2/c. По чертежу КРК (рис. П2.4) определяем Ds = 0,464 м, DL = 0,7645 м, Zu0=1,028 м. Атмосферное давление Рат = 10 м.

Все необходимые геометрические размеры РК, НА и обводов приведены на рис.П2.1-П2.3.

Результаты расчетов гидродинамической нагрузки на облицовку 3-го пояса КРК и ее спектральные характеристики показаны на рис. П2.7 (F = Ω/2π).

Далее проводится расчет напряженного состояния облицовки КРК. Рассматривается цельнометаллическая конструкция, изготовленная из биметаллического листа (ст. 3 и 08Х13). Соответствующие прочностные характеристики приведены в табл. П2.1.

Таблица П.2.1

Марка стали Предел прочности σвр, МПА Предел текучести σт, МПА Предел выносливости σ-1, МПА
Ст.3 350 215 122,5
08Х13 600 420 210

В расчете принято, что был нарушен контакт металла с бетоном и между ними имелась щель 2 мм. Фильтрационное давление принималось равным гидростатическому. Толщина облицовки варьировалась и, соответственно, менялись прочностные характеристики.

Результаты расчетов в виде максимальных средних напряжений (σ) и стандартов напряжений (σ‘) приведены в сводной табл. П2.2. Используя их, можно оценить надежность и долговечность облицовки (см. табл. П.2.2 и рис.П.2.8).

Таблица П.2.2

Расчеты проводятся до тех пор, пока долговечность и вероятность безотказной работы не достигнут, соответственно, значений 100 лет и 1. Некоторый запас по долговечности взят из-за того, что по применяемой методике вычисляются средние значения долговечности, тогда как из-за разброса прочностных характеристик металлов и случайности других, входящих в расчет факторов, полученные в результате расчета долговечности имеют разброс и являются случайными величинами.



Рис. П.2.1. Очертания проточной части гидротурбины ПЛ60-В-600/1075
(D1 = 6 м, D0 = 7,2 м, b0 = 2,1 м);
кромки НА и РК показаны для режима Н = 46 м, N =113 МВт





Рис. П.2.2. Профиль лопатки НА (N0 = 24) для ПЛ60-В-600/1075



а)



б)



в)



г)



д)



е)



Рис. П.2.3. Чертеж лопасти РК ПЛ60/1075 (для модели D1 = 0,35 м, N1 = 7):
а - цилиндрическое сечение R = 0,0753 м; б - R = 0,0951 м; в - R = 0,1141 м;
г - R = 0,1331 м; д - R = 0,1522 м; е - R = 0,1712 м





Рис. П.2.4. Геометрия облицовок КРК ПЛ60-В-600/1075:
а - поперечный разрез камеры рабочего колеса; б - развертка 1/4 части камеры РК









Рис. П.2.6. Универсальная характеристика гидротурбины ПЛ 60/1075-В-46

а)



б)



Рис. П.2.7. Результаты расчета гидродинамической нагрузки на облицовку 3-го пояса КРК ПЛ60-В-600/1075 для режима H = 59 м, N = 113 МВт:
а - распределение давления под облицовкой; б - коэффициенты сk для спектральной плотности gpp нагрузки на облицовку КРК





СПИСОК ЛИТЕРАТУРЫ

1. Дзюбанов Е.М., Дмитриев Н.Ю., Климович В.И., Левина С.М., Штильман В.Б. О прогнозировании надежности камер рабочих колес диагональных и поворотно-лопастных гидротурбин // Гидротехническое строительство, 1996. № 2. С. 8-15.

2. Климович В.И. Расчет течений жидкости в проточной части гидромашин на основе последовательного решения осесимметричной задачи и задачи обтекания решеток профилей // Известия ВНИИГ им. Б.Б. Веденеева. 1990. Т. 220. С. 88-93.

3. Климович В.И. Квазитрехмерный расчет течений жидкости в проточных частях гидромашин//Изв. АН СССР. МЖГ. 1991. № 2. С. 69-76.

4. Этинберг И.Э., Раухман Б.С. Гидродинамика гидравлических турбин. Л: Машиностроение, 1978.

5 Корн Г., Корн Т. Справочник по математике. М.: Наука, 1968.

6 Справочник по гидротурбинам. Л.: Машиностроение, 1984.

7 Гидроэнергетическое и вспомогательное оборудование гидроэлектростанций. Справочное пособие. М.: Энергоатомиздат, 1988 (т. 1). 1990 (т. 2).

8. Мошков Л.В. Собственные частоты колебаний воды в трубе постоянного сечения при наличии аэрации //Известия ВНИИГ им. Б.Е. Веденеева. 1969. Т. 88. С. 48-55

9. Ржаницин А.Р. Теория расчета строительных конструкций на надежность. М.: Стройиздат, 1978.

10. Болотин В.В. Прогнозирование ресурса машин и конструкций. М.: Машиностроение, 1984

11 Серенсен С.В.., Когаев B.П., Шнейдерович Р.М. Несущая способность и расчеты деталей машин на прочность/Справочное пособие. М: Машиностроение, 1975.

СОДЕРЖАНИЕ

Предисловие

1. Общие положения

2. Расчет гидродинамической нагрузки на облицовку КРК со стороны проточной части гидротурбины


2.1. Общие положения

2.2. Модель течения жидкости в проточной части гидротурбины. Расчет гидродинамической нагрузки на КРК

2.3. Расчет спектральных характеристик гидродинамической нагрузки на элемент облицовки КРК

2.4. Определение режимных параметров ПЛ-гидротурбины, необходимых для расчета гидродинамической нагрузки

2.5. Подготовка исходных данных для гидродинамического расчета проточной части ПЛ-гидротурбины

3. Расчет динамического напряженного состояния стальной облицовки

3.1. Выбор расчетной схемы

3.2. Нагрузка на элемент облицовки

3.3. Основные предпосылки для прочностного расчета элемента облицовки

3.4. Постановка задачи гидроупругости

4. Оценка надежности стальных облицовок КРК ПЛ-гидротурбин

4.1. Оценка вероятности безотказной работы элемента облицовки

4.2. Оценка долговечности элемента облицовки

4.3. Оценка долговечности элемента облицовки с учетом реального графика работ

4.4. Оценка надежности облицовки КРК как системы

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

Приложение 2. Пример расчета гидродинамических нагрузок, напряженного состояния и долговечности облицовок КРК

Список литературы