На правах рукописи
Дубинский Сергей Иванович
Численное моделирование
ветровых воздействий
на высотные здания и комплексы
Специальность 05.13.18 –
математическое моделирование, численные методы и комплексы программ
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата технических наук
Москва – 2010
Работа выполнена в Государственном образовательном учреждении высшего профессионального образования Московский государственный строительный университет.
Научный руководитель:
|
доктор технических наук, профессор
Белостоцкий Александр Михайлович
|
Официальные оппоненты:
|
доктор технических наук, профессор
Травуш Владимир Ильич
|
кандидат технических наук,
старший научный сотрудник
Попов Николай Александрович
|
Ведущая организация:
|
НИИ Механики МГУ,
г. Москва
|
Защита диссертации состоится «___» _____________ 2010 г. в ______ часов на заседании диссертационного совета Д 212.138.12 при ГОУ ВПО Московский государственный строительный университет по адресу: 129337, г. Москва, Ярославское шоссе. д.26, ауд. ______.
С диссертацией можно ознакомиться в библиотеке ГОУ ВПО Московского государственного строительного университета.
Автореферат разослан «___»___________ 2010 г.
Ученый секретарь
диссертационного совета Анохин Н.Н.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы. Обоснование безопасности и комфортности проектируемых и уже эксплуатируемых высотных зданий современных архитектурных форм и оригинальных конструктивных решений и их комплексов в России (ММДЦ “Москва-Сити”, ЖК "Аквамарин" во Владивостоке и ряда других), для которых фактор ветрового воздействий является определяющим, сдерживается действующими нормативно-регламентированными методиками, которые не содержат рекомендаций по назначению аэродинамических коэффициентов для сложных по форме высотных зданий, не учитывают влияние интерференции, рельефа, изменения спектра набегающего потока в условиях плотной и изменяющейся застройки.
Нормами (СНиП 2.01.07-85* «Нагрузки и воздействия», МДС 20-1.2006 «Временные нормы по назначению нагрузок и воздействий ... на многофункциональные высотные здания и комплексы в Москве») предлагается использовать результаты испытаний крупномасштабных макетов в специализированных аэродинамических трубах, позволяющих воспроизвести атмосферный пограничный слой. Но такие испытания весьма трудоемки, причем для наиболее масштабных моделей их приходится выполнять за рубежом. Кроме того, методология экспериментального моделирования ветровых потоков и воздействий на высотные комплексы обладает собственными ограничениями и погрешностями.
В последние 10-15 лет бурно развивается вычислительная гидрогазодинамика (CFD), совершенствуются технологии расчетов ветровых воздействий на здания и сооружения при неуклонно возрастающей мощности компьютеров. Ведущие зарубежные исследовательские и проектные организации все чаще комбинируют испытания и "численные" эксперименты. В перспективе роль математического моделирования, как показал опыт в смежных отраслях (например, аэрокосмической) и задачах (строительная механика), будет только возрастать.
Разработка, программная реализация и применение верифицированной методики, основанной на численном решении трехмерных уравнений гидрогазодинамики, позволит отечественным специалистам определять ветровые воздействия на высотные комплексы с необходимыми достоверностью, полнотой и оперативностью.
Цели и задачи работы. Целью работы являлась разработка, программная реализация и верификация методики расчета ветровых воздействий на высотные здания и комплексы на основе численного решения задач гидрогазодинамики, позволяющей оперативное многовариантное решение большеразмерных практических задач строительного проектирования и мониторинга в форме вычислительного эксперимента.
Исходя из поставленной цели работы, решались следующие задачи:
Анализ результатов отечественных и зарубежных теоретических и экспериментальных исследований ветровых воздействий на высотные здания и их комплексы в условиях, прежде всего, их компактного расположения и интерференции, а также анализ применимости существующей нормативной базы;
Анализ и тестирование существующих методологий численного моделирования задач гидрогазодинамики и выбор программных средств с точки зрения их применимости к решению поставленных задач;
Разработка комплексной методики расчета ветровых воздействий на высотные здания-комплексы, основанной на численном моделировании уравнений гидрогазодинамики и обеспечивающей решение с единых позиций регламентируемых нормами задач: расчета средней и пульсационной составляющих ветровых нагрузок на несущие конструкции, определения пиковых давлений на ограждающие конструкции, оценки уровня комфортности в пешеходных зонах комплексов.
Bыбор наилучших (по соотношению “ресурсы-точность”) методологий построения расчетных сеток, моделей турбулентности, характеристик вычислительных алгоритмов применительно к данному классу задач и выбранному базовому программному комплексу (ПК);
Разработка "инженерной" методики оценки пульсационной составляющей давления и пиковых нагрузок на фасадные конструкции по результатам стационарных расчетов;
Разработка процедур передачи ветровых нагрузок в программы расчетов динамики и прочности конструкций, зданий и сооружений;
Реализация вышеупомянутых методик в форме единого программного модуля – надстройки над базовым ПК;
Верификация разрабатываемой методики и реализующего программного обеспечения на основе сравнения результатов расчетов с результатами испытаний в аэродинамических трубах и с данными натурных замеров на представительном наборе примеров;
Применение разработанной методики и программного обеспечения для определения ветровых воздействий на реальные проектируемые и строящиеся высотные комплексы;
Обучение студентов и аспирантов теоретическим основам и практическому применению разработанной методики.
Исследования проводились в рамках работ по Проекту 2.1.2/6414 аналитической ведомственной целевой программы "Развитие научного потенциала высшей школы (2009-2010 годы) ".
Научная новизна работы.
Разработана и верифицирована методика определения ветровых воздействий на высотные здания и их группы на основе решения уравнений Навье-Стокса в приближении Рейнольдса (модели турбулентности RANS, URANS и DES) с применением метода конечных объемов, позволяющая адекватно учитывать ключевые факторы – направления и профиль ветровых потоков, рельеф местности, окружающую застройку и поэтапность возведения комплексов – и обеспечивающая успешное практическое решение регламентируемых отечественными нормами задач (расчета средней и пульсационной составляющей нагрузок на несущие конструкции, определения пиковых давлений на ограждающие конструкции, оценки уровня пешеходной комфортности).
Разработана методика оценки пульсационной составляющей давления и пиковых нагрузок на фасадные конструкции по результатам стационарных расчетов осредненной энергии турбулентных пульсаций и с учетом осредненных коэффициентов обеспеченности и корреляции.
Разработана методология проблемно-ориентированного программного модуля для проведения вариантных вычислительных экспериментов, позволяющего в автоматизированном режиме формировать полный программный код для всех этапов расчета.
Практическая значимость работы.
Разработанные в диссертации методика математического моделирования и вычислительного эксперимента, а также реализующие ее программные средства обеспечивают приемлемую практическую точность, позволяют существенно уменьшить временные ресурсы на выполнение расчётов ветровых воздействий (средняя и пульсационная составляющие нагрузок на несущие конструкции, пиковые давления на фасадные конструкции, оценка уровня пешеходной комфортности) на высотные здания-комплексы;
Полученные с использованием разработанной методики результаты могут быть применены для верификации аналогичных технологий расчета и специализированных программных комплексов, а также для “настройки” вновь вводимых в эксплуатацию аэродинамических труб, для повышения эффективности, сокращения сроков испытаний и оптимального размещения датчиков;
Методика позволяет создание справочной базы данных по аэродинамическим параметрам для типовых форм высотных зданий, допускающей использование совместное с программными комплексами расчетов динамики и прочности конструкций;
Методика может быть использована в составе систем мониторинга перемещений и ускорений высотных зданий и комплексов для повышения эффективности их работы;
Методика может быть применена для определения ветровых воздействий на произвольные классы сооружений, включая транспортные.
Внедрение:
Методика и программное обеспечение применяются в Научно-образовательном центре компьютерного моделирования (НОЦ КМ) МГСУ и Научно-исследовательском центре “СтаДиО” для практических расчетов реальных высотных зданий и комплексов застроек;
Разработанная методика используется в практике обучения студентов кафедры “Информатика и прикладная математика” МГСУ по дисциплине “Вычислительная аэрогидромеханика”, а также аспирантов, прикрепленных к НОЦ КМ МГСУ для выполнения диссертационных исследований.
Личный вклад соискателя. Все исследования, изложенные в диссертационной работе, проведены лично соискателем в процессе научной деятельности. Из совместных публикаций в диссертацию включён лишь тот материал, который непосредственно принадлежит соискателю.
На защиту выносятся:
Методология численного моделирования ветровых воздействий на высотные здания и их группы, обеспечивающая успешное практическое решение регламентируемых отечественными нормами задач (расчета средней и пульсационной составляющей ветровых нагрузок, определения максимальных давлений на ограждающие конструкции, оценки уровня пешеходной комфортности);
Методика оценки пульсационной составляющей давлений и максимальных расчетных нагрузок на фасадные конструкции по результатам стационарных расчетов энергии турбулентных пульсаций и с учетом осредненных коэффициентов обеспеченности и корреляции;
Методология проблемно-ориентированного программного модуля для проведения вариантных вычислительных экспериментов, позволяющего в автоматизированном режиме формировать полный программный код для всех этапов расчета, для анализа и оформления результатов расчета, а также для передачи ветровых нагрузок в пакеты “прочностных” расчетов в качестве исходных данных;
Методология формирования расчетных аэродинамических моделей высотных сооружений и их комплексов в составе застройки.
Достоверность полученных результатов обеспечивается:
использованием апробированного математического аппарата (математические модели теории аэромеханики и механики деформируемого твёрдого тела) и численных методов решения;
применением апробированных в мировой практике технологий аэродинамических расчетов широкого круга задач машиностроения и строительства и верифицированного лицензионного программного комплекса ANSYS CFX, их реализующего;
успешным решением с использованием разработанной методики представительного набора верификационных тестовых задач;
согласованием получаемых результатов с данными расчетов авторитетных отечественных и зарубежных специалистов, использующих иные методики и программные средства;
согласованием с результатами экспериментальных исследований в аэродинамических трубах и данными натурных замеров для реальных зданий и застроек.
Апробация работы. Основные положения диссертации докладывались и обсуждались на симпозиумах Российской академии архитектуры и строительных наук (РААСН) «Актуальные проблемы компьютерного моделирования конструкций и сооружений» (Нижний Новгород, 2007 и Пермь, 2008), на конференциях «Теория и практика расчета зданий, сооружений и элементов конструкций. Аналитические и численные методы» (Москва, МГСУ, 2008, 2009 гг.), на научных семинарах НОЦ КМ МГСУ "Актуальные проблемы компьютерного моделирования зданий, сооружений и комплексов" (в декабре 2009 и феврале 2010), а также на конференциях и семинарах пользователей ПК ANSYS (2001-2002, 2004-2008 гг.) и ПК SCAD в г. Москва (2006-2009 гг.) и в г. Киев (2008 г.).
Публикации. По тематике диссертации опубликована 11 работ, в том числе 7 работ в изданиях, включенных ВАК в перечень рекомендуемых.
Структура и объем работы. Диссертационная работа состоит из введения, пяти глав (с выводами по каждой главе), заключения, списка литературы (194 наименования, в том числе – 110 на иностранных языках), 81 рисунка и 11 таблиц. Общий объём диссертации – 199 страниц.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении приводится обоснование актуальности работы, определены проблемы, цели и задачи исследований, перечислены основные научные и практические результаты, приведено краткое содержание диссертации по главам.
В первой главе дается аналитический обзор современного состояния вопросов определения ветровых воздействий (нормативного, натурно-метеорологического, экспериментального и численного) на высотные здания и их комплексы. Представлен анализ отечественной нормативной базы в ее развитии и опыта практических расчетов на ветровые нагрузки (Э.И. Реттер, А.Г. Соколов, Г.А. Савицкий, И.М. Беспрозванная, Ф.Л. Серебровский, А.И. Цейтлин, А.С. Бернштейн, М.Ф. Барштейн, Б.Г. Коренев, М.И. Казакевич, А.А. Петров, Н.А. Попов, В.И. Травуш, А.В. Перельмутер, Б.В. Остроумов и др.). Из зарубежных исследователей можно отметить вклад А.G. Davenport, E. Simiu, R. H. Scanlan, N.J. Cook, Y. Tamura, A. Kareem, J.D. Holmes, G. Solari и др.
Анализируются существующие отечественные и зарубежные экспериментальные оценки ветровых нагрузок на здания и комплексы (Р.А. Irwin, А.G. Davenport, Н.А. Попов, М.А. Березин, А.Б. Айрапетов, С.В. Гувернюк, В.Г.Гагарин, Б.В. Остроумов, Ю.В. Табунщиков, С.Г.Кузнецов и др.). Рассмотрены особенности аэродинамических испытаний и ветровой аэродинамики высотных зданий в условиях их компактного расположения и интерференции на примере сооружений ММДЦ “Москва-Сити”. Обсуждаются реальные ветровые режимы для гг. Москвы и Владивостока и нормативное районирование (Л.Е. Анапольская, М.М. Борисенко, Н. А. Дашко и др.).
Кратко описана история математического моделирования задач гидрогазодинамики. Рассмотрены как классические работы (И. Ньютон, Д. Бернулли, Ж.Л. Даламбер, О. Рейнольдс, Д. Стокс, Н.Е. Жуковский, Л. Прандтль, Т. фон Карман, А.Н. Колмогоров, Л.И. Седов, Л.Г.Лойцянский, Г. Шлихтинг и др.), так и современные отечественные (М.Х. Стрелец, Е.М. Смирнов, С.А. Исаев, А.Е. Усачов, И.К. Лифанов, О.М. Белоцерковский, С.В. Гувернюк, А.В. Атаманчук, В.Я. Шкадов, В.Н. Варапаев и др.).
Изучены существующие методические рекомендации по численному моделированию ветровой аэродинамики (J. Franke, С.P.W. Guerts, A.Mochida). Анализируются возможности существующих программных комплексов (ПК) и опыт расчетов ветровых воздействий (Stathopoulos, Blocken, Carmeliet, Bitsuamlak, Meroney, Franke, Kai Fan Liаw, F.Menter, R. Stangroom, С.В. Гувернюк, С.А. Исаев и др.). Отмечено, что применение наиболее мощных и признанных ПК (в частности, ANSYS CFX, ANSYS FLUENT, STAR-CD, PHOENICS) обеспечивает успешное решение практических задач.
Обосновывается необходимость разработки специализированного программного модуля, учитывающего особенности расчетов ветровых воздействий на высотные здания и комплексы.
Во второй главе описана разработка комплексной методики численного моделирования ветровых воздействий на высотные здания и их компактно расположенные группы.
Расчеты ветровых потоков и воздействий сводятся к численному решению трехмерных нестационарных нелинейных уравнений гидрогазодинамики в постановке Навье-Стокса:
 (1)
Кроме того, должны удовлетворяться уравнения неразрывности (сохранения массы) и состояния:
 (2)
 (3)
Здесь u,v,w – искомые компоненты вектора скорости (по осям x,y,z), p –давления, t – время, –динамический коэффициент вязкости для воздуха, ρ – плотность, R – универсальная газовая постоянная, T –температура.
Для упрощения моделирования ветровые потоки предполагаются несжимаемыми и изотермическими, массовые силы не учитываются.
Прямое решение уравнений (1-3) с учетом вихрей всех масштабов (DNS, Direct Numerical Simulation) при современных возможностях ЭВМ практически реализуемо только для очень малых скоростей потока и чисто исследовательских задач. Поэтому, в современной расчетной практике доминирует полуэмпирический подход, основанный на разложении скорости на осредненную во времени и пульсационную составляющие  и переходе к решению т.н. "осредненных по Рейнольдсу уравнений Навье-Стокса" (Reynolds averaged Navier-Stoks Method, RANS):
 ,
  (4)
где  – средние давления, индексы i=1,2,3 и j=1,2,3 соответствуют координатам x,y,z. Сдвиговые (рейнольдсовы) напряжения  являются дополнительными шестью неизвестными к параметрам осредненного движения ( ,) и аппроксимируются, как правило, по гипотезе Буссинеска:  , (5)
где Т – дополнительная вязкость, вызванная пульсациями; k – осредненная энергия турбулентных пульсаций (TKE). Система является незамкнутой и требует дополнительных соглашений ("моделей турбулентности").
Важным аспектом решения уравнений Навье-Стокса является схема дискретизации. Наиболее эффективным в настоящее время признан "метод конечных объемов" (МКО), не требующий столь детального моделирования пограничного слоя, как метод конечных элементов (МКЭ), и более удобный при описании сложных расчетных областей реальных застроек, чем метод конечных разностей. Современные подходы (схемы адвекции второго порядка, применение пристеночных функций и увеличение числа точек интегрирования в ячейке) значительно снижают требования к расчетной сетке и ресурсам ЭВМ.
Решение уравнений (4) может проводиться как в стационарной, так и в нестационарной постановке. Стационарная задача сводится к итерационному решению редкозаполненной системы линейных алгебраических уравнений, с узловыми давлениями, компонентами скоростей и мерами локальной турбулентности (в частноcти, TKE) в качестве неизвестных. Определяются осредненные во времени величины, число которых для рассматриваемых объектов/задач может достигать сотен миллионов. Итерационный процесс решения системы проводится по связанной многосеточной схеме (coupled multigrid) и считается сошедшимся при малых относительных невязках решения системы. Для улучшения сходимости в качестве предобуславливателя используется неполное треугольное (ILU) разложение матрицы.
К сожалению, точность (и даже сходимость) стационарных расчетов RANS в ситуациях с развитым вихреобразованием не всегда удовлетворительна. Проведенное представительное тестирование на широком круге практических задач (включая сравнение с данными испытаний в аэродинамических трубах) позволило определить наиболее пригодную для данного класса задач модель турбулентности. Модель SST (Shear Stress Transport, перенос сдвиговых напряжений, Ф. Ментер,1993), комбинируя достоинства классических моделей “k-ε” и “k-ω”, обеспечивает достаточную точность результатов и эффективную сходимость итерационного процесса даже при относительно грубых сетках и умеренном разрешении пограничного слоя.
Выработаны рекомендуемые размеры расчетной области – минимум 10 высот "господствующего" задания во всех направлениях. Рассмотрены варианты построения сетки в расчетной области и типы ячеек-элементов.
В качестве характеристик набегающего потока используются профили скорости ветра и характеристик турбулентности (интенсивность и масштаб вихря), отвечающие заданным ветровым районам и типам местности по СНиП. Предлагается в качестве "внешних" граничных условий для расчетной области использовать "мягкие" условия по Нейману (равенство нулю производных).
Для учета шероховатости стенок рекомендуется вместо подбора пристеночных функций или коррекции их стандартных параметров применять по возможности "натуральное" моделирование (учет в модели рельефа местности, балконов и оконных проемов и т.п.).
Предложена и апробирована следующая схема расчетных исследований:
1) по результатам стационарных расчетов для всех направлений ветра определяются наиболее неблагоприятные направления по максимальным значениям средних нагрузок и максимальной энергии турбулентных пульсаций на поверхностях; для этих направлений выполняется уточненный расчет при сгущенной сетке и/или с применением схем автоматической адаптации;
2) для этих направлений выполняется нестационарный расчет; при существенном отличии средних значений может оказаться необходимым выполнить нестационарные расчеты для многих или "всех" направлений ветра. Нестационарные расчеты также следует выполнять для направлений, при которых стационарный расчет указывает на возникновение "стоячей волны".
Нестационарные расчеты RANS (Unsteady RANS, URANS) также не всегда позволяют правильно отследить срывные потоки. В этой связи наиболее активно в настоящее время развиваются "гибридные" подходы. При DES-подходе (моделировании “отсоединенных” вихрей, Detached Eddy Simulation, М.Х. Стрелец и P.R. Spalart,1997) комбинируются LES (моделирование крупных вихрей, Large Eddy Simulation) и URANS, что позволяет во много раз уменьшить трудоемкость расчетов по сравнению с “классическим” LES. Вихри в пограничном слое и вблизи него моделируются интегрально, отсоединенные вихри предлагается моделировать детально (аналогично LES), применяя достаточное разрежение сетки. При подходе SAS (Scalable Adaptive Simulation, Menter, Egorov, 2005) комбинируются уже DES и URANS, в процессе нестационарного расчета гибко изменяется предельный размер учитываемых вихрей.
Неявные схемы интегрирования позволяют применять временной шаг с числами Куранта Co>3 при условии отслеживания отклика системы. Для нестационарных расчетов необходима генерация нестационарного поля скоростей ветра со средними и спектральными характеристиками, отвечающими заданной высоте и ветровому району. Используется суммирование гармонических функций со случайными фазами согласно выражению:
 , (6)
где: N – количество частот, для которых рассчитывается спектр S(f), Δf –заданное приращение (шаг) для частот; текущая частота fк =k⋅Δf, t – время, угол φ - случайная величина, равномерно распределенная в диапазоне [0;2π].
Ввиду крайней трудоемкости вариантных нестационарных расчетов для минимизации их объема предлагается практическая методика оценки и пиковых расчетных нагрузок на фасадные конструкции (P*max и P*min) по результатам стационарных расчетов энергии турбулентных пульсаций TKE и с учетом осредненных коэффициентов обеспеченности θmax и θmin. Пульсационная составляющая нагрузок на конструкцию может быть получена в результате интегрирования по поверхности пульсационной составляющей давлений Pпульс .
Приведем основные расчетные соотношения методики:
TKE =3/2 (I V)2 = 3 P/ ρ I2
I = (ρ TKE / abs (P)/3)1/2 , p = (I2+2 I) abs(P) (7)
P*max = P (1+ p θmax), P*min = P (1- p θmin)
Pпульс = ( P*max - P*min ) ν/2
Здесь: р – среднеквадратичное отклонение (стандарт) давлений P, I – интенсивность турбулентности (стандарт пульсаций скорости), V – локальная скорость ветра, ν – коэффициент корреляции давлений на поверхностях (по СНиП 2.01.07-85*, либо вычисленный из анализа нестационарных результатов).
Для предварительных расчетов можно принять θmin=6 и θmax=3 (по рекомендациям А.Davenport, 2005); уточнение значений выполняется по результатам аналогичных нестационарных расчетов по соотношениям:
θmax =(Pmax-P)/ р , θmin =(P - Pmin)/ р, (8)
где Pmax, Pmin и P – вычисленные значения минимумов, максимумов и среднего давлений в точке поверхности, р – стандарт давлений.
Высокие скорости ветра, образующиеся вокруг высотного здания, неблагоприятно воздействуют на человека, особенно при низких температурах, а в некоторых случаях являются опасными. В этой связи, нормативные документы МДС 20-1.2006 требуют, чтобы скорость ветра (в порывах на уровне 1.5 м) превышала Vcr1=6 м/с не чаще Kcr1 =1000 часов в год (1-й уровень комфортности); превышала Vcr2 =12 м/с не чаще Kcr2 =50 часов в год (2-й уровень комфортности); превышала Vcr3 =20 м/с не чаще Kcr3=5 часов в год (3-й уровень комфортности). После проведения численного моделирования для “всех” направлений ветра (как правило, j=1,…,24) "время дискомфорта l-го уровня" Kcrl (l =1,2,3) для представительного набора точек пешеходных зон определяется соотношениями:
Kcrl =Σ Sij Tij, Vij= Vi/V10 (Vj + θ . I) (9)
где: Vi=1,2,3,... Vimax – скорости в таблице метеоданных ("розе ветров"); Tij - продолжительность (по метеоданным, часов в год) ветрового воздействия направления j и средней скорости Vi; Vj – средняя скорость ветра в данной точке согласно расчету для направления j при скорости V10 на высоте 10 м; Vij – максимальная скорость в точке в порывах при скорости ветра Vi; θ – коэффициент обеспеченности по заданию на расчет (обычно в диапазоне от 1 до 3); Sij – признак (0 или 1) превышения локальной скоростью ветра в точке Vij критического значения Vcrl для данного уровня комфорта l.
В главе 3 описана структура и возможности специализированного программного модуля WINDLOAD/CFX, реализующего разработанную численную методику.
В качестве основного расчетного инструмента выбран лицензионный универсальный ПК ANSYS, установленный в МГСУ и допускающий также проведение (в развитие диссертационных исследований) прочностных и связанных аэроупругих расчетов.
Поскольку расчеты на ветровые воздействия требуют задания множества нестандартных (для универсальных расчетных программ) параметров расчета и опций, определено, что для повышения эффективности вариантных расчетов и минимизации ошибок ввода наиболее целесообразно программно формировать текстовый стартовый файл, содержащий все необходимые признаки и опции для расчетной области (домена) и для граничных условий на внешних границах расчетной области, а также параметры рекомендуемых моделей турбулентности, опции решения, точки мониторинга и записываемые результаты.
Разработанные расчетные процедуры были оформлены в специальный программный модуль WINDLOAD/CFX, включающий в настоящее время около 80 "макросов" (процедур) на языке программирования APDL.
В рамках ПК ANSYS создаются и файлы с сеткой (*.cdb) и программный код (*.ccl) для ПК ANSYS CFX. Макросы снабжены русифицированным интерфейсом и рассчитаны на специалистов в области прочности, не имеющих большого опыта численного моделирования задач газодинамики. Макросы составлены для различных типовых форм сооружений, для двумерной и трехмерной постановок, для различных вариантов учета симметрии, для различных подходов к формированию сетки. В рамках одного входного файла можно провести генерацию сетки в цикле (например, при вариациях геометрическиx параметров задачи и для вариантов застройки). Разработаны макросы для формирования стационарных и нестационарных профилей скорости ветра и параметров турбулентности, а также процедуры методик, описанных в главе 2. Для упрощения и повышения эффективности обработки и оформления результатов средствами программирования APDL создается и файл настроек постпроцессора (в формате *.cst). Формируются плоскости сечений, ракурсы рисунков, вспомогательные параметры и переменные. Определяются спектры воздействия (давлений и нагрузок на характерные зоны), а также числа Струхаля
St=Df/V, (10)
где D – характерный размер (как правило, поперек потока), V – средняя скорость потока, f – частота срыва вихрей.
Рассмотрены преимущества одновременного расчета модели застройки сразу при нескольких направлениях ветра для оптимизации использования многопроцессорной вычислительной техники и особенности применения кластеров к решению данного класса задач. Описана методика импорта геометрии зданий и передачи ветровых нагрузок (средней и пульсационной составляющей) из ANSYS CFX в программы расчетов НДС зданий (ANSYS, SCAD, Лира, MicroFE, СТАДИО).
|
