Наука к Образование
МГТУ им. Н.Э. Баумана
Сетевое научное издание
Наука и Образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2015. № 09. С. 336-351.
ISSN 1994-0408
Б01: 10.7463/0915.0802462
Представлена в редакцию: Исправлена:
© МГТУ им. Н.Э. Баумана
11.08.2015 09.09.2015
УДК 534.1
Уточнение конечно-элементной модели лопатки ГТД на основе результатов вибрационных испытаний с учётом разброса модальных параметров
Николаев С. М.1'*, Жулёв В. А.1, Киселёв И. А.1
1МГТУ им. Н.Э. Баумана, Москва, Россия
nikolaev-b gmail.com
В данной статье рассматривается вопрос уточнения конечно-элементной модели лопатки компрессора низкого давления, закреплённой в технологической оснастке с использованием стохастического подхода на основе результатов модальных испытаний. Цель уточнения модели - сведение к минимуму различия экспериментальных и расчётных собственных частот и форм колебаний. Предложенный алгоритм позволяет достичь указанной цели с помощью вариации приведённых модулей упругости элементов сборочной конструкции в области контакта с учётом разброса исследуемых частот. В работе выведены рекуррентные зависимости, лежащие в основе предлагаемого алгоритма. На примере уточнения модели лопатки компрессора газотурбинного двигателя показана эффективность предлагаемого метода для получения конечно-элементных моделей, адекватно описывающих динамику исследуемого изделия.
Ключевые слова: уточнение модели, экспериментальный модальный анализ, конечно-элементное моделирование, стохастический подход
Введение
Целью математического моделирования механических систем в динамике является исследование характеристик поведения объектов под действием различных нагрузок, переменных во времени. С помощью модели можно определить такие параметры системы как собственные частоты и формы колебаний, а также найти наибольшую амплитуду колебаний, распределение напряжений и деформаций при различных режимах эксплуатации конструкции. В настоящее время широкое распространение получил метод конечных элементов, который позволяет определить вышеупомянутые характеристики для изделий со сложной геометрией. Однако, для достоверного определения динамических параметров конечно-элементная модель должна с высокой точностью соответствовать объекту исследования. Выполнение такого требования в большинстве случаев затруднительно, так как, во многих случаях, объекты исследования является сборочными конструкциями, а точное моделирование мест соединений и стыков при
динамическом поведении конструкции не всегда возможно из-за большого количества неизвестных параметров.
Таким образом, возникает задача идентификации параметров конечно-элементной модели в местах соединений. Обычно, такими параметрами являются приведенные модули упругости элементов конструкции в местах соединений. В данном случае, приведенные модули упругости несут физический смысл параметров контактной жесткости, неизвестной заранее.
Как известно, собственные частоты и формы колебаний являются базовыми динамическими характеристиками любой механической системы. При этом методика испытаний по определению собственных частот колебаний являются достаточно хорошо изученной и простой для реализации. Методы определения собственных частот из экспериментальных данных подробно описаны в работах D. Ewins, P. Sas, R. Allemang, D. Brown и др. [1-3].
Разницу значений собственных частот, полученных экспериментально и численно, удобно использовать в качестве критерия точности расчётной модели. Наряду с собственными частотами, в качестве критерия соответствия модели и объекта, могут использоваться и собственные формы колебаний. Однако, получение экспериментальных собственных форм более трудоёмкий процесс и, в рамках данной статьи, на рассматривается.
Первые публикации в области уточнения и верификации модели стали появляться в конце 60-х - начале 70-х годов. Данные методы были улучшены и обобщены на случай, когда параметры механической системы известны с некоторой конечной точностью, то есть являются случайными величинами с известными статистическими характеристиками. Наиболее подробные математические выкладки проделаны в работе [4]. Несмотря на большое количество исследований, данная область продолжает развиваться в настоящее время [5-8]. Последние достижения в области уточнения расчетных моделей механических систем связаны с применением теорем теории вероятности и статистики, в частности с использованием теоремы Байеса, для определения функции плотности вероятности системы [9]. Стоит отметить, что вопрос автоматизированного уточнения конечно-элементных моделей сложных механических систем по результатам модальных испытаний является актуальной проблемой для внедрения автоматизации проектирования в различных отраслях высокотехнологичной промышлености.
Постановка задачи
В рамках работы решается задача уточнения конечно-элементной модели лопатки компрессора низкого давления, закрепленной в технологической оснастке. Динамические характеристики модели должны с заданной степенью точности соответствовать динамическим характеристикам партии изделий. Значения собственных частот колебаний партии лопаток являются случайными величинами, с известными из модальных испытаний статистическими характеристиками (математичекие ожидания и дисперсии). Такая постановка задачи требует применения стохастического подхода к уточнению конечно-элементной модели, который позволяет учесть случайную природу параметров модели и динамических характеристик объекта.
Получение рекуррентных соотношений для верификации конечно-элементной модели с учетом статистического разброса экспериментальных данных
Как было отмечено выше, в ходе проведения модальных испытаний собственные частоты и формы колебаний, а также коэффициенты демпфирования могут быть определены с некоторым статистическим разбросом, который необходимо учитывать при верификации модели. Ниже приведена последовательность получения рекуррентных соотношений для верификации и уточнения конечно-элементной модели по результатам модальных испытаний с учетом разброса параметров системы. При этом, выкладки приведены для общего случая, когда в качестве параметров модели могут быть выбраны любые величины.
Рассмотрим характеристическое уравнение механической системы п-го порядка без учета демпфирования
(К - р2М) q = 0, (1)
где К — матрица жёсткости системы, М — матрица масс системы, р — собственная частота колебаний системы, Ц — собственный вектор системы.
Можно ввести некоторую неизвестную функцию, которая описывает связь между собственными частотами системы и параметрами системы.
Данную функцию можно разложить в ряд Тейлора до первого члена [7], тогда соотношение между вектором собственных частот и параметрами системы может быть представлено в векторно-матричной форме:
Ш = Ш (хр ) + т-(х — ^ ). (2)
вектор идеальных параметров объекта, вектор параметров модели,
„ / 2 2 2 \ вектор, составленный из ур1,р2,...,рп),
матрица производных частот от параметров, то есть
где х х
Ш
Т
Т =
дщ дщ дщ
дх1 дх2 дхп
дщ дщ дсо2
дхх дх2 дхп
дщп дсоп дюп
дх1 дх2 дхп
Рассмотрим процедуру решения уравнения (2) с учетом нормального разброса экспериментальных данных. Перепишем (2) в виде
У = т • Г (3)
где
У = ю
- Ю ( х , )
Г = х - X.
(4)
(5)
Предполагается, что параметры модели имеют нормальное Гауссовское распределение. Тогда, математическое ожидание параметров объекта
я=ЕИ=х Р >
ковариационная матрица вектора х определяется как:
8хх = Е [( х-Ях )( х-Ях )Т ] = Е ( х - ^ )( х - ^ )Т
Математическое ожидание и ковариационная матрица вектора Г :
я = Е [г ] = Е [ х - х, ] = 0,
8„ = Е [ г • гТ ] = Е [(х - х ,)(х - х, )Т ] = 8хх, а математическое ожидание и ковариационная матрица вектора У :
Я = Е [У ] = 0, Е [у • уТ ] = Е [Т • г • гТ • ТТ ] = Т • Е [г • гТ ] • ТТ = Т • 8„ • Т
(6)
(7)
(8) (9)
(10) (11)
Рассмотрим систему (3) с учётом погрешности £ , которая возникает в результате измерений. Данная погрешность образуется как вследствие несовершенства аппаратуры и алгоритмов обработки экспериментальных данных, так и вследствие разброса параметров при переходе от одного изделия партии к другому.
Примем, что данная погрешность имеет нулевое математическое ожидание и ковариационную матрицу 8££ . Данная матрица получена из модальных испытаний и обычно является диагональной матрицей дисперсий.
у = Т • Г + £ . (12)
Для новой системы среднее значение у совпадает со значением у
я=Е [у ]=0
Матрица ковариации у принимает новый вид
8уу = Е [у • у] = Е [(Т • Г + £)(т • Г+£ )Т ]
Матрица ковариаций векторов у и Г определяется как
8уг = Е [у • гТ ] = Е [(Т • г + £ )• гТ ] = Е [Т • г • гТ + £ • гТ ] = Т • SГ
= Т • 8 • тт + 8„
(13)
(14)
(15)
Т
причем, вектора Г и 8 статистически независимы, следовательно,
£[е ■ гт ] = 0. (16)
Цель приведенных математических операций - получить наилучшее приближение параметров модели Г на основе измеренных параметров натурного объекта. Следовательно, необходимо решить обратную задачу:
Г * = о ■ у, (17)
где матрица О должна обеспечивать наименьшую разницу между действительным
вектором r и полученным вектором Г , то есть норма матрицы
\T'
Q = E
(Г* -r)(Г* -Г)
(18)
должна стремиться к минимальному значению.
Рассмотрим выражение (18) с учетом (17)
Q = E[(G• y - r)(G• y - r)T] = G• Syy • GT - G• Syr - S^ • GT + S„, (19)
и проварьируем (19) по G. Поскольку норма матрицы Q минимальна, то её вариация равна нулю
0 = • (Syy • GT - Syr ) + (G • Syy - Syr) • ^ . (20)
Отсюда следует, что
G = ^r • (21)
Г * = Sir •• y . (22)
Заметим, что матрица G является псевдообратной матрицей T, то есть
g = t+. (23)
Объединив уравнения (4), (5), (14) и (15) получаем
x* = x, + srr • tt (t • srr • tt + se )-1 • r, (24)
s*r = srr - srr • ttTt • srr • tt + se)-1 • t • srr. (25)
Таким образом, на основе работ [5-9] авторами статьи получены рекуррентные соотношения, которые позволяют верифицировать численную модель механической системы с учетом разброса экспериментальных данных. Если зависимость между параметрами модели и её частотами будет иметь линейный характер, то верификация будет осуществлена за одну итерацию. В действительности зависимость нелинейная и потребуется несколько итераций использования формул (24) и (25).
Получение зависимостей для определения коэффициентов
чувствительности
Следует отметить, что основная трудоёмкость алгоритма уточнения состоит в вычислении матрицы Т, состоящей из частных производных собственных частот по параметрам. В случае, когда в качестве параметров выступают модули упругости в конечных элементах, искомую производную оказывается возможно получить аналитически [10].
Для нахождения собственной частоты можно использовать соотношение Релея:
щ=< Кч,
чТ Мч,
(26)
J max
Покажем, как изменяется собственная частота щ при изменении механических
свойств материала в элементе п. Введем параметр х., который будет описывать
отклонение механических свойств (плотности и модуля упругости) от номинальных значений. Коэффициент х. может изменяться в некоторых пределах от х до х1
Р(^) = ^ 'Р,
Е (х ) = х • Е ,
где
р - номинальная плотность материала; Е - номинальный модуль упругости.
Определим производную собственной частоты по параметру х.. В соответствии с соотношением (26) получим:
(27)
йщ
1
йх 2щ чТ
2 ^ (К -щ М) ч, + чТ
йх
йк
йх
\ J
йМ
щ
йх
J У
(28)
Первое слагаемое может быть исключено в соответствии с (1)
йщ
1
йх 2щ Мд
йк
йх
\ J
щ
йм
йх
(29)
J У
В случае, если плотность и модуль упругости изменяются только в конечном элементе с номером е, производные от глобальной матрицы масс и жесткости могут быть выражены через изменения матрицы жесткости и масс конечного элемента с номером е :
-М =-Ме
й к = л. К е •
(30)
йх
йх
Используя формулы (29) и (30), получаем коэффициент чувствительности е-го элемента, соответствующий /-ой собственной частоте
dЩ! „ dщ 1 /
2 7 1 _ / лт« Л
_
е ул^е ^"""е <-у"л'е У
< = —_ = 2щ ■
dx„ dx„ т
1т
_=- (ч е)
ТЫI 4 7
2 дме ок щ2---
дх дх„
ч е, (31)
где Х^ - коэффициент чувствительности элемента е соответствующий собственной форме /,
т - модальная масса,
Ц е - вектор собственной формы /, соответствующий степеням свободы конечного элемента е .
Производные от глобальных матриц жесткости и масс в соотношении (31) будут равны соответствующим матрицам конечного элемента е
<=1 (д е )т (щ
тх 'У
хе —(яе) (щ-клIяе. (32)
В случае если изменяется только модуль упругости, либо только плотность материала конечного элемента, соотношение (32) принимает вид:
2
щ / \Т
е I _е »
Xе =Щ(це ) Меце, т р 4 '
1 / \т х=—(це) ке це.
т р
(33)
тЕ
Коэффициент чувствительности в виде (33) показывает, что уменьшение плотности (увеличение жесткости) в элементе с максимальным коэффициентом чувствительности приведет к наибольшему (по сравнению с изменением свойств в других КЭ) увеличению _-ой собственной частоты системы.
Следует обратить особое внимание на физический смысл коэффициента чувствительности. При варьировании модулей упругости конкретных конечных элементах, коэффициент чувствительности представляет собой масштабированную потенциальную энергию упругой деформации е-го конечного элемента на /-ой собственной частоте системы.
Таким образом, получены все необходимые соотношения для реализации алгоритма уточнения конечно-элементной модели на основе результатов испытаний с учетом разброса модальных параметров системы. Авторами статьи впервые предложен алгоритм, реализующий зависимости 23, 24, 32, 33 совместно. Использование зависимостей 32, 33 позволяет значительно сократить время работы алгоритма, заменяя численное вычисление матрицы чувствительности её вычислением по аналитическим зависимостям.
Алгоритм реализован в авторском программном обеспечении, позволяющем выполнять уточнение модели в автоматизированном режиме. Блок-схема данного алгоритма представлена на рис. 1.
Рисунок 1 - Алгоритм верификации модели с учётом разброса параметров
Этапы верификации модели
В качестве объекта верификации рассматривается лопатка компрессора низкого давления, закреплённая в технологической оснастке (рис. 2). Эксперимент проводился при помощи однокомпонентного датчика ускорения и виброударного молотка фирмы B&K, и аналого-цифрового преобразователя National Instruments. Полученные данные обрабатывались в программном комплексе LabView. Процесс идентификации
собственных частот и форм колебаний лопатки с помощью виброударных испытаний подробно описан в работах [1,2].
Рисунок 2 - Лопатка компрессора низкого давления, закреплённая в технологической оснастке
В результате эксперимента были получены первые три собственные частоты серии лопаток. Результаты испытаний представлены в таблице 1.
Таблица 1 - Результаты модальных испытаний лопатки ГТД
№ собственной формы Частота, Гц Дисперсия, Гц2
1 208.1 3.7
2 779.5 9.7
3 1104.9 19.3
При моделировании сборки контакт в области закрепления лопатки в оснастке и самой оснастки в станочных тисках не учитывается, т.е. конечно-элементная модель сборочной конструкции выполнена, как одна деталь, состоящая из разных материалов. Конечно-элементная модель представлена на рис. 3.
Рисунок 3 - Конечно-элементная модель лопатки
Первые три собственные формы колебаний изображены на рис. 4.
Рисунок 4 - Собственные формы колебаний лопатки
Значения соответствующих собственных частот представлены в таблице 2.
Таблица 2 - Результаты определения частот лопатки газо-турбинного двигателя методом конечных
элементов
№ собственной формы Частота, Гц
1 230.4
2 886.3
3 1124.7
Для автоматизированного уточнения модели была разработана программа, реализующая алгоритм, изображенный на рис. 1, с рекуррентными зависимостями (24, 25). Уточнение проводилось в области контакта лопатки и оснастки путем корректировки приведённых значений модуля упругости в каждом конечном элементе данной области. В качестве целевой функции, подлежащей минимизации была использована норма вектора ошибки расчетных собственных частот системы. Вектор ошибки рассчитывался следующим образом:
8 =
/ РЕЛ
Р
V
•100%
(34)
Результаты верификации представлены в таблице 3.
Таблица 3. Результаты уточнения лопатки в области закрепления
№ собственной формы Частота эксперимент альная, Гц Частота расчетная до уточнения, Гц Ошибка, % Частота расчетная после уточнения, Гц Ошибка, %
1 208.1 230.4 10.7 201.6 - 3.1
2 779.5 886.3 13.7 774.0 - 0.7
3 1104.9 1124.7 1.8 1104.9 0.0
На рис. 5 представлены графики сходимости собственных частот в процессе уточнения модели.
Как видно из приведённого графика, после уточнения модели погрешность по всем трём собственным частотам колебаний лопатки не превышает 4%. Таким образом, уточненную модель лопатки можно использовать для точного моделирования динамики изделия, в частотном диапазоне до 1150 Гц. На рис. 6 представлена уточненная модель лопатки. Цвет каждого конечного элемента модели соответствует значению приведенного модуля упругости в данном элементе после уточнения.
2 4 6 8 10 12 14 16 18 20 22 24
Номер итерации
Рисунок 5 - Графики сходимости собственных частот
Рисунок 6 - Изображение уточненной модели лопатки
Как видно из рис. 6 наибольшие изменения претерпели значения модулей упругости элементов, находящихся близко к контактной области лопатки и технологической оснастки. Это связано с тем, что в данной области, коэффициенты чувствительности собственных частот колебаний к изменению параметров упругости являются максимальными.
Стоит отметить, что для моделирования сложных переходных процессов, необходимо не только выполнять уточнение модели по собственным частотам колебаний, но и идентифицировать модель диссипативных сил, возникающих в системе.
Заключение
В рамках данной работы, по результатам модальных испытаний серии лопаток компрессора газотурбинного двигателя, было выполнено уточнение конечно-элементной модели данной лопатки с учетом статистического разброса параметров. Для уточнения модели был использован стохастический алгоритм, основанный на анализе чувствительности. Задача уточнения математической модели лопатки была сведена к минимизации различия между измеренными и расчетными значениями собственных частот колебаний с помощью варьирования приведенных модулей упругости в месте закрепления лопатки. До выполнения процедуры уточнения модели, собственные частоты лопатки, полученные в результате модальных испытаний отличались от результатов моделирования более чем на 13%. После процедуры уточнения модели максимальная погрешность по собственным частотам составила 3.79%. Таким образом, предложенный алгоритм может быть использован для уточнения моделей сложных сборочных конструкций, с учетом разброса параметров системы.
Уточненная конечно-элементная модель изделия может быть использована для прогнозирования динамического поведения под действием эксплуатационных нагрузок, точного назначения остаточного ресурса и оптимизации конструкции.
Список литературы
1. Ewins D.J. Modal Testing: Theory, Practice and Application. 2nd edition. Baldock: Research Studies Press LTD, 2000.
2. Heylen W., Lammens S., Sas P. Modal Analysis Theory and Testing. Leuven: KUL Press, 1997.
3. Brown D.L., Allemang R.J., Zimmerman R., Mergeay M. Parameter estimation techniques for modal analysis. SAE Technical Paper No. 790221. SAE, 1979. 19 p. DOI: 10.4271/790221
4. Collins J.D., Hart G.C., Hasselman T.K., Kennedy B. Statistical Identification of Structures // American Institute of Aeronautics and Astronautics Journal. 1974. Vol. 12, no.2. P. 185190. DOI: 10.2514/3.49190
5. Katafygiotis L.S., Papadimitriou C., Lam H.-F. A probabilistic approach to structural model updating // Soil Dynamics and Earthquake Engineering. 1998. Vol. 17, is. 7-8. P. 495-507. DOI: 10.1016/S0267-7261(98)00008-6
6. Mares C., Mottershead J.E., Friswell M.I. Stochastic model updating: Part 1 — theory and simulated example // Mechanical Systems and Signal Processing. 2006. Vol. 20, no. 7. P. 1674-1695. DOI: 10.1016/j.ymssp.2005.06.006
7. Mottershead J.E., Mares C., James S., Friswell M.I. Stochastic model updating: Part 2 - application to a set of physical structures // Mechanical Systems and Signal Processing. 2006. Vol. 20, no. 8. P. 2171-2185. DOI: 10.1016/j.ymssp.2005.06.007
8. Mottershead J.E., Link M., Friswell M.I. The sensitivity method in finite element model updating: A tutorial // Mechanical Systems and Signal Processing. 2011. Vol. 25, no. 7. P. 2275-2296. DOI: 10.1016/j.ymssp.2010.10.012
9. Behmanesh I., Moaveni B. Bayesian FE Model Updating in the Presence of Modeling Errors // Model Validation and Uncertainty Quantifications. Vol. 3. Springer International Publishing, 2014. P. 119-133. DOI: 10.1007/978-3-319-04552-8 12
10. Huang X., Zuo Z.H., Xie Y.M. Evolutionary topological optimization of vibrating continuum structures for natural frequencies // Computers & Structures. 2009. Vol. 88, no. 5. P. 357-364. DOI: 10.1016/j.compstruc.2009.11.011
Science and Education of the Bauman MSTU, 2015, no. 09, pp. 336-351.
DOI: 10.7463/0915.0802462
Received: Revised:
11.08.2015 09.09.2015
Science^Education
of the Bauman MSTU
ISS N 1994-0408 © Bauman Moscow State Technical Unversity
Updating Turbine Blade FE Model Considering the Modal Parameters Scatter
S.M. Nikolaev1'*, V.A. Zhulev1,1.A. Kiselev1
nikolaev.b oxig gmail.com :Bauman Moscow State Technical University, Moscow, Russia
Keywords: model updating, experimental modal analysis, finite-element modeling, stochastic
approach
The paper considers an issue of updating a turbine blade finite element model. The turbine blade is fixed in the special rig and the unknown contact conditions strongly influences on the blade dynamics. Hence, there is a difference between the real blade dynamic characteristics and the results of the finite element simulation. The task is to update the blade FE model on the basis of the modal test data to obtain fine dynamic model of the structure. The experimental modal parameters are known with some scatter, which is the result of measurement errors, and some differences from blade to blade within the kit of parts.
The proposed algorithm is capable to update the turbine blade finite element model by variation of elasticity modulus in contact areas taking into account the experimental modal parameters scatter. It is worth noting that elasticity modulus of finite elements near the contact area should be considered as unknown contact stiffness parameters.
The stochastic updating algorithm based on the sensitivity analysis is proposed. The paper presents equations to define sensitivity coefficients and provides the special recurrence relations, on which the proposed algorithm is based. The eigenvalue problem is resolved repeatedly with system parameters changing during the algorithm implementation. The system parameters (elasticity modulus in the contact areas) are recalculated according to the recurrence relations, mentioned above. To implement the proposed algorithm, a special software in C++ language has been developed.
The proposed algorithm has demonstrated its efficiency for updating the finite element model of the complex structure. Thus, the differences between the calculated eigenfrequencies and those of identified from modal testing have been reduced from 13% to 3%. The updated model of the turbine blade could be used for accurate simulation of the blade dynamic behavior under different operating conditions.
References
1. Ewins D.J. Modal Testing: Theory, Practice and Application. 2nd edition. Baldock, Research Studies Press LTD, 2000.
2. Heylen W., Lammens S., Sas P. Modal Analysis Theory and Testing. Leuven, KUL Press, 1997.
3. Brown D.L., Allemang R.J., Zimmerman R., Mergeay M. Parameter estimation techniques for modal analysis. SAE Technical Paper No. 790221. SAE, 1979. 19 p. DOI: 10.4271/790221
4. Collins J.D., Hart G.C., Hasselman T.K., Kennedy B. Statistical Identification of Structures. American Institute of Aeronautics and Astronautics Journal, 1974, vol. 12, no.2, pp. 185190. DOI: 10.2514/3.49190
5. Katafygiotis L.S., Papadimitriou C., Lam H.-F. A probabilistic approach to structural model updating. Soil Dynamics and Earthquake Engineering, 1998, vol. 17, is. 7-8, pp. 495-507. DOI: 10.1016/S0267-7261(98)00008-6
6. Mares C., Mottershead J.E., Friswell M.I. Stochastic model updating: Part 1 — theory and simulated example. Mechanical Systems and Signal Processing, 2006, vol. 20, no. 7, pp. 1674-1695. DOI: 10.1016/j.ymssp.2005.06.006
7. Mottershead J.E., Mares C., James S., Friswell M.I. Stochastic model updating: Part 2 - application to a set of physical structures. Mechanical Systems and Signal Processing, 2006, vol. 20, no. 8, pp. 2171-2185. DOI: 10.1016/j.ymssp.2005.06.007
8. Mottershead J.E., Link M., Friswell M.I. The sensitivity method in finite element model updating: A tutorial. Mechanical Systems and Signal Processing, 2011, vol. 25, no. 7, pp. 2275-2296. DOI: 10.1016/j.ymssp.2010.10.012
9. Behmanesh I., Moaveni B. Bayesian FE Model Updating in the Presence of Modeling Errors. Model Validation and Uncertainty Quantifications. Vol. 3. Springer International Publishing, 2014, pp. 119-133. DOI: 10.1007/978-3-319-04552-8 12
10. Huang X., Zuo Z.H., Xie Y.M. Evolutionary topological optimization of vibrating continuum structures for natural frequencies. Computers and Structures, 2009, vol. 88, no. 5, pp. 357-364. DOI: 10.1016/j.compstruc.2009.11.011