CC BY
25
СТРОИТЕЛЬНЫЕ КОНСТРУКЦИИ, ЗДАНИЯ И СООРУЖЕНИЯ
УДК 515.2: 624.04
А. А. Абдюшев - старший научный сотрудник Тел.: (843) 272-42-45, e-mail: YZSRR@kfti.knc.ru
И.Х. Мифтахутдинов - кандидат технических наук, старший научный сотрудник Тел.: (843) 231-91-15, e-mail: mift1943@mail.ru
П.П. Осипов - доктор физико-математических наук, профессор, ведущий научный сотрудник Тел.: (843) 231-91-17, e-mail: ossipov@kfti.knc.ru
Институт механики и машиностроения Казанского научного центра Российской академии наук (ИММ КНЦРАН)
ПРОЕКТИРОВАНИЕ НЕПОЛОГИХ ОБОЛОЧЕК МИНИМАЛЬНОЙ ПОВЕРХНОСТИ
АННОТАЦИЯ
Разработан программный комплекс для построения минимальных поверхностей, широко применяемых в архитектуре. Для построения оболочек со сложной конфигурацией применена технология многоуровневых адаптивных сеток. Приведено сравнение результатов построения оболочек по двум моделям, описываемым уравнениями Лапласа и Эйлера-Лагранжа. Дан анализ сходимости итерационных алгоритмов вычисления и точности представления минимальных поверхностей по различным моделям.
КЛЮЧЕВЫЕ СЛОВА: непологие оболочки, минимальные поверхности, уравнения Эйлера-Лагранжа, Лапласа.
A.A. Abdiushev - senior staff scientist Tel.: (843) 272-42-45, e-mail: YZSRR@kfti.knc.ru
I.H. Miftakhutdinov - candidate of technical sciences, senior staff scientist Tel.: (843) 231-91-15, e-mail: mift.1943@ma.i1.ni
P.P. Osipov - doctor of physical-mathematical sciences, leading staff scientist Tel.: (843) 231-91-17, e-mail: ossipov@kfti.knc.ru
Institute Mechanical engineers and Machine building of Kazan Scientific Centre of Russia Academy of Sciences (IMM KSC RAS)
DESIGN OF STEEP SHELLS WITH MINIMAL SURFACE
ABSTRACT
The simulator generating widely used in architecture minimal surfaces is developed. The technology of multylevel adaptive grids for generation of shells having complex shape is applied. The comparison of numerical results for two models corresponding to Laplace and Euler-Lagrange equations is presented. The analysis of convergation of iteration algorithms and quality of minimal surface representation for two models is given.
KEYWORDS: steep shells, minimal surfaces, Euler-Lagrange, Laplace equations.
Введение
Формы минимальных поверхностей (МП) и уравнения, их описывающие, известны из работ Л. Эйлера (1755) и Ж.Л. Лагранжа (1760). В 1861 году бельгийский физик Ж. Плато показал всему миру эти поверхности в виде мыльной пленки, натянутой на проволочный каркас. В настоящее время широко известен ряд минимальных поверхностей [1, 2, 3, 4]. В экспериментальной работе [2] показано, что при описании МП можно приближенно пользоваться уравнением Лапласа. Условно назовем получаемые при этом поверхности поверхностями Лапласа (ПЛ). Широкому применению МП в архитектуре и
строительстве способствует их привлекательный архитектурно-выразительный внешний вид, поскольку они созданы самой природой. Их оптимальность заключается в минимальном расходе материалов, а также в уменьшении теплообмена с внешней средой. Начало применению МП в архитектуре и строительстве положил известный немецкий архитектор и инженер Ф. Отто. В лаборатории Штутгартского университета (ФРГ) под его руководством проведено большое количество экспериментов с мыльными пленками, тонкими резиновыми шнурами, тканевыми и металлическими сетями для нахождения МП. Эти
эксперименты нашли воплощение во многих тентовых и сетчатых (вантовых) конструкциях, запроектированных и построенных под руководством Ф.Отто. Среди них павильон ФРГ в Монреале (1967) из тентовой конструкции, светопрозрачное сетчатое покрытие трибун Олимпийского стадиона в Мюнхене (1972), прозрачная конструкция выставочного павильона в Мангейме (1978) и др.
Отечественными учеными также внесен определенный вклад в разработку теории МП и внедрению их в практику строительства: тентовое покрытие под куполом в молодежном лагере близ Казани (1967), тентовые и полимерные покрытия над рыночными павильонами Казани и Татарстана (Блинов Ю.И., Сладков В. А., Удлер ЕМ., Мифтахутдинов ИХ и др.).
В настоящее время появились дорогостоящие программные комплексы (МА8ТЯАМ, АМБУБ, Соєтоє’МзгкБ), позволяющие, в частности, рассчитывать геометрию оболочек, заданных на сложных в плане областях. В настоящей работе представлены результаты вычислений геометрии МП с использованием авторских программных разработок.
Математическая модель
Уравнение поверхности. МП характеризуются тем, что при заданных граничных условиях их средняя кривизна в любой точке равна нулю [2, 4]. Аппликаты этих поверхностей в декартовой системе координат удовлетворяют уравнению Эйлера-Лагранжа
¡1 +
дУ
+ ¡1 +
д2 z dz dz д2 z
„ —2— 2------------------+
дх дх ду дхду
дх
дУ2
или в дивергентной форме
д дх
( д z Ї д / д z"
m v дх 0 + дУ m V ду 0
= 0,
= 0,
2
m = ¡1+
¿z
дх
+
дУ
(1)
Для пологих оболочек подъемистости вдоль осей малы, т.е. дz / Эх << 1 и Эz / Эу << 1, и поэтому ц » 1. Для них хорошим приближением служит решение уравнения Лапласа, которое получается из (1) при
т = 1
д2;
- + -
д2;
дх2 ду2
0.
(2)
Для подъемистых оболочек величина ¡Л может значительно отличаться от единицы, и уравнение Лапласа неадекватно описывает МП.
Граничные условия. Основные граничные условия, применяемые в архитектуре и строительстве, сводятся
к условиям Дирихле. На границе G области задается проектная высота поверхности (аппликата)
z\G = w( x(t), y(t)), или ZG = w(t), (3)
где x(t) , y(t) - параметрические уравнения участков границ.
Алгоритм определения поверхности
Записав уравнение (1) в конечных разностях на девятиточечном шаблоне равномерной сетки (рис. 1)
и обозначив Tt j = z ( x t, y ■ ), получим:
Ti, j = g 1Ti+1, j + g 2 Ti, j+1 +
+ 73Ti-1, j + g 4 Ti, j-1 , (4)
где g k =----------—-----------, (k = 1,2,3,4),
m 1 + m 2 + m 3 + m 4
т ={1+ ((Т+1,,■ - Т,)/н )2 + ((т,+1 +
+ Т+1, ,■+1 - Т+1, ,■-1 - Т,, -1) /(4 н) )2 ¡-Л
т 2 = { + ((Т+1,, +1 + Тг+1,, - Тг-1,, + 1 -
- ,)/(4н))2 + ((Т,,.+1 - Ти,)/н)2
т з = {1+((Т„ - Т-1,,)/ н )2 + ((Т-1,,+1 + (5)
+ Ти+1 - Ти-1 - Т-и-1)/(4н))2 }-05,
т 4 = { + ((Т+1,, + Тг +1,,-1 - Тг-1,,1 -
- Т-1,,-1)/(4н))2 + ((Т,, - Ти,-1)/н)2 }-0 5.
Уравнение (4) решается с помощью метода последовательных приближений, т.е. итераций
(п = 0,1,2,...)
Тп+1 = ^ пТп + V пТп +
11,, П1, + 1, 1211,, + 1 "Г,
+ V зпТ-1,, + гт п,-1 (6)
0.5
^+1 гГ1
__________________________*¡+1^+1
1 ^2 7- ■ Мл у 1
^4
тил 71+1^-1
Рис. 1. Девятиточечный шаблон
где индексы /', ] пробегают по внутренним узлам сетки. Процесс устойчив, поскольку все ук > 0 и
У]” + У2 + у3 + У4 = 1. Устойчивость гарантирует
то, что итерация не расходится, но не гарантирует ее сходимости. Наши численные эксперименты неизменно показывали сходимость метода (6).
Для пологих оболочек все тк » 1 и уП » 1/4.
Заметим, что при численном решении уравнения Лапласа можно использовать итерацию (6), положив
все Ук = 1/4. В этом случае достаточно пятиточечного шаблона «крест».
Расчет МП на равномерных сетках
В качестве тестового примера была выбрана поверхность Шерка [1, 3]
2 = 1 1о§(соБ(ау)/ СОБ(ах)}. (7)
а
Был рассмотрен случай а = 1 для диапазона х, у е [- ат^(к), ат^(к)]; к > 0- Здесь к -
модуль градиента функции г(х, у) на границе
рассматриваемой области. Можно убедиться, что эта поверхность точно удовлетворяет уравнению (1). На рис. 2 показаны соответствующие участки поверхности
Шерка при разных значениях подъемистости к.
В таблице 1 приведены максимальные относительные погрешности расчета узловых аппликат поверхности Шерка (7) по уравнениям Лапласа (2) и Эйлера-Лагранжа (1). Погрешности отнесены к максимальному диапазону аппликат. Была использована равномерная квадратная сетка при трех вариантах числа разбиения N.
Расчет проводился до тех пор, пока модуль средней разности значений аппликат для двух последовательных
итераций не становился меньше 5 • 10 7 . Расчеты проводились с привлечением программных средств АРС ЭРА ПК2000 [5]. Анализ результатов показывает, что для построения МП применение уравнения Лапласа является приемлемым с практической точки
при подъемистости к £ 1. При к > 1 построение МП необходимо проводить с применением уравнения Эйлера-Лагранжа [4]. При этом следует отметить, что сходимость итерационного процесса к ПЛ, организованного по пятиточечному шаблону, происходила в среднем в два раза быстрее процесса итерирования к МП по девятиточечному шаблону.
Расчет МП на нерегулярных сетках
При построении сеток для произвольных областей использована технология многоуровневых сеток [6]. Такие сетки состоят из квадратных ячеек, получаемых многократной бисекцией крупной декартовой сетки.
Рис. 2. Участки поверхности Шерка. Слева: к = 0.5; справа: к = 10
Таблица 1
N Лаплас Эйлер-Лагранж
к = 0.2
10 0.000833 0
30 0.000904 0
50 0.000908 0
к = 0.5
10 0.004820 0.00001 1
30 0.005203 0.000001
50 0.005219 0.000000
к = 1.0
10 0.015577 0.000139
30 0.016607 0.000016
50 0.016621 0.000006
к = 5.0
10 0.083611 0.006885
30 0.084412 0.001285
50 0.084476 0.000509
к - = 10.0
10 0.119564 0.014674
30 0.120620 0.004707
50 0.120703 0.001925
Типичный пример ячейки и смежных ячеек изображен на рис. 3. Ранг ячейки характеризует число бисекций (размер ячейки). Чем выше ранг, тем больше ячейка. При этом правило рангов требует, чтобы ранг ячейки отличался не более чем на единицу от рангов смежных ячеек. В центре ячейки задана аппликата. Для получения девятиточечного шаблона определяются аппликаты на серединах сторон и в углах ячейки
(жирные точки на рис. 4 и 5). Для этого используется линейная интерполяция. В таблице 2 приведены
формулы вычисления аппликаты Л1 для середины первой стороны в зависимости от рангов и положения смежных ячеек. Для других сторон формулы
вычисления аппликат Л2, Л3, Л4 получаются аналогично.
^00 Х2
®01 а. Su
^00 ^10 X
S„‘ x4
Рис. 4. Варианты вычисления аппликаты в середине первой стороны
Таблица 2
1 -ячейка справа имеет тот же ранг 2 - две ячейки справа имеют меньший ранг 3 - ячейка справа имеет больший ранг, чем iSjj 4 - ячейка справа имеет больший ранг, чем Л’ю
, Aj + Я Лі= 2 , Я + Я,, + Я,, д _ а и_ 1 3 Ai = —Я + -Я ~ 6 3 6 Л _ 5 Î 1 1 1 1 -/V і — Æ H- Æа 6 3 6
Рис. 5. Варианты вычисления аппликаты для угла 12
Таблица 3
1 - угловая ячейка имеет тот же ранг 2 - угловая ячейка имеет меньший ранг 3 - угловая ячейка имеет больший ранг 4 - угловая ячейка имеет больший ранг и является одновременно верхней ячейкой
Л 1 Î 1 Î ■/Yi^ — Æ + Л*\ 2 2 1 3 Л„ =-Я + -Я, 3 3 2 1 к =-Я + -А, 12 3 3 1 Л .г = —Я + —Я, +-Я2 3 3 3
Рис. 6. Минимальная поверхность сложной конфигурации
Рис. 7. Многоуровневая сетка в двусвязной области на плоскости (х, у)
Рис. 8. Сравнение линий уровня а - МП, б - ПЛ
В таблице 3 приведены формулы вычисления аппликаты для угла 12 в зависимости от рангов и положения смежных ячеек (рис. 5). Для других углов формулы получаются аналогично.
Результаты расчета
На рис. 6 показана МП, заданная на многосвязной области, имеющей в плане сложную конфигурацию (рис. 7).
На рис. 8 приведены линии уровня минимальной поверхности и поверхности Лапласа при средней подъемистости поверхности порядка 10. Разница высот в характерных точках достигает 30 %.
Заключение
Разработанный авторами программный комплекс позволяет оперативно проводить многовариантные расчеты геометрии оболочек МП. На основе проведенных вычислений авторы пришли к выводу, что для оболочек МП с подъемистостью более чем 1 необходимо применять уравнение Эйлера-Лагранжа. Для построения пологих МП с подъемистостью меньшей, чем 1, допускается использование уравнения Лапласа. Технология многоуровневых сеток, предложенная в статье для построения сложных в плане МП, позволяет получать сетки для конечно-
элементного или конечно-разностного анализа напряженного состояния оболочек. Перспективными задачами для авторов является расширение разработанного программного комплекса для расчета напряженного состояния оболочек МП.
Литература
1. Бляшке В. Введение в дифференциальную геометрию. - М.: Изд-во технико-теоретической литературы, 1957. - 223 с.
2. Мифтахутдинов И.Х. Оболочки минимальной поверхности в природе и архитектуре. - Казань: ЗАО «Новое знание», 2007. - 144 с.
3. Мифтахутдинов. И.Х. Наглядная геометрия оболочек минимальной поверхности. - Казань: ЗАО «Новое знание», 2009. - 40 с.
4. Михайленко В.Е., Ковалев С.Н. Конструирование форм современных архитектурных сооружений. -Киев: Будивельник, 1978. - 112 с.
5. Лукашенко В.И., Абдюшев А.А., Доронин М.М., Нуриева Д.М., Сладков А.В. Экспертиза, расчет, анализ пространственных конструкций (АРС ЭРА-ПК2000). - Казань: КазГАСУ, 2006. - 320 с.
6. Simpson R.B. Automatic local refinement for irregular rectangular meshes. // Int. J. Num.Meth. Eng, 1979, vol.14. - P. 1665-1678.
