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

Полный текст

Введение. Решение задач оптимизации свойств упругих тел направлено на улучшение некоторых характеристик (максимизация первой резонансной частоты или критической силы) за счет выбора формы или переменности упругих модулей. Подробный исторический обзор приводится в монографии [1], где отмечается, что задачи оптимизации упругих конструкций берут свое начало с работы Г. Галилея по оптимизации профиля консольной балки, нагруженной силой на конце. Таким образом, первый тип задач оптимизации связан с нахождением законов изменения геометрических характеристик (в основном толщины для пластин и формы поперечного сечения для колонн и балок), обеспечивающих минимальный вес при заданной нагрузке, или, что эквивалентно, максимальную прочность при заданном весе. Второй тип задач связан с максимизацией критической силы в задачах устойчивости колонн [2, 3], где фундаментальные результаты принадлежат Ж. Лагранжу, Т. Клаузену и Е.Л. Николаи. Третий тип задач – динамические, как и задачи устойчивости, связаны со спектральными задачами [4]. Динамические задачи по оптимизации частот собственных колебаний струны и балки были впервые рассмотрены в работах М.Г. Крейна [5] и Ф. Ниордсона [6].

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

С ростом применения на практике армированных композитов возникали задачи по оптимизации внутренней структуры упругих тел, которые в большинстве случаев решались на основе моделей кусочно-однородных и анизотропных тел. Однако в работах [8, 9] решены задачи об оптимизации первой собственной частоты для продольно неоднородной балки. В обеих работах получены оптимальные законы изменения безразмерного модуля Юнга, причем в [8] учтены предварительные напряжения, а в [9] приведены результаты для всех основных типов граничных условий (консоль, шарнирное опирание, закрепленные торцы).

В настоящее время с развитием технологий производства материалов и конструкций из функционально-градиентных материалов (ФГМ) стала актуальной задача по отысканию непрерывных оптимальных распределений плотности, пористости, упругих модулей, позволяющих увеличить резонансную частоту. Так, например, в работе [10] при помощи использования генетического алгоритма решена задача об определении переменного модуля упругости для увеличения критической силы. В работе [11] рассмотрена задача увеличения критической силы для колонны, где в качестве управляющих функций фигурируют как модуль Юнга, так и форма поперечного сечения.

Большинство работ по максимизации первой резонансной частоты выполнено в рамках вариации геометрических характеристик. В работе [12] изучается задача о максимизации основной резонансной частоты башни за счет оптимизации формы ее поперечного сечения. Оптимальное решение для диаметра и толщины секций разыскивается в классе кусочно-линейных функций. Предложено использование МКЭ с кубической аппроксимацией прогиба и реализован алгоритм по отысканию оптимального решения. В работе [13] изучена задача об оптимизации первой собственной частоты для круглой анизотропной пластинки. Управляющей функцией является переменная толщина пластинки, а решение задачи строится при помощи метода разложения по малому параметру для определенного типа анизотропии. В работе [14] описан итерационный алгоритм решения оптимизационной задачи, основанный на методе возмущений. Решена задача об оптимизации частоты колебаний при помощи управления площадью поперечного сечения. В работе [15] в рамках градиентного подхода осуществляется максимизация первой резонансной частоты для трехмерного стержня с использованием МКЭ.

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

1. Общие условия оптимальности. Рассмотрим установившиеся колебания неоднородного тела, занимающего объем V и ограниченного гладкой поверхностью S. Уравнения колебаний имеют вид:

${\sigma }_{ij,j}^{}+{\mathrm{\rho \omega }}^2{u}_i^{}=\mathrm{0,}i=\mathrm{1,}\mathrm{2,}\mathrm{3,}$ (1.1)

определяющие соотношения запишем в виде закона Гука для анизотропной среды:

${\sigma }_{ij}=c_{ijkl}{u}_{k,l},$ (1.2)

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

Будем считать, что граница S состоит из двух частей S_u, $S_{\sigma }$, на которых заданы однородные граничные условия:

${\left.u_i^{}\right|}_{S_u}=\mathrm{0,}{\left.{\sigma }_{ij}^{}{n}_j\right|}_{S_{\sigma }}=0.$ (1.3)

Здесь n_j – компоненты вектора внешней нормали к поверхности S.

Умножим (1) на u_i и проинтегрируем по объему V. Получим:

$\underset{\mathrm{V}}{\int }{\sigma }_{ij,j}^{}{u}_i^{}dV+{\mathrm{\omega }}^2\underset{\mathrm{V}}{\int }\mathrm{\rho }{u}_i^{}{u}_i^{}dV=0.$

Используя теорему Гаусса–Остроградского и граничные условия (1.3), получим:

$J={\mathrm{\omega }}^2=\frac{{\displaystyle \underset{V}{\int }{c}_{ijkl}{u}_{i,j}{u}_{k,l}dV}}{{\displaystyle \underset{V}{\int }\mathrm{\rho }{u}_i^{}{u}_i^{}dV}}=\frac{\mathrm{П}}{K}.$ (1.4)

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

Будем считать, что среднее по объему значение модуля фиксировано:

$\frac{1}{V}\underset{V}{\int }E\left(x\right)dV=m_0,$ (1.5)

и введем функционал

$J_1=\frac{1}{V}\underset{V}{\int }E\left(x\right)dV-m_0.$ (1.6)

Затем сформируем функционал:

$J_{*}\left[E\right]=J+\mathrm{\Lambda }{J}_1,$ (1.7)

где $\mathrm{\Lambda }$ – множитель Лагранжа.

Проварьируем функционал (1.7) и найдем его первую вариацию:

$\begin{array}{l}\delta J_{*}={\left.\frac{\partial }{\partial t}{J}_{*}[E+t\delta E]\right|}_{t=0}=\\ =\frac{{\displaystyle \underset{V}{\int }\left(\delta E{R}_{ijkl}{u}_{i,j}{u}_{k,l}+2E{R}_{ijkl}\delta u_{i,j}{u}_{k,l}\right)dV}}{{\displaystyle \underset{V}{\int }\mathrm{\rho }{u}_i{u}_{i}dV}}-2{\mathrm{\omega }}^2\frac{{\displaystyle \underset{V}{\int }\left(\mathrm{\rho }\delta u_i{u}_i\right)dV}}{{\displaystyle \underset{V}{\int }\mathrm{\rho }{u}_i{u}_{i}dV}}+\frac{\mathrm{\Lambda }}{V}\underset{V}{\int }\delta EdV.\end{array}$

Условие стационарности функционала J_* дает

$\mathrm{\delta }{J}_{*}=\underset{V}{\int }[\mathrm{\Phi \delta }E+{\mathrm{\Phi }}_i\mathrm{\delta }{u}_i]dV=\mathrm{0,}$ (1.8)

где введены следующие обозначения:

$\begin{array}{l}\mathrm{\Phi }=R_{ijkl}{u}_{i,j}{u}_{k,l}+2\mathrm{\Lambda }{K}_0{V}^{-1},\\ {\mathrm{\Phi }}_i=2\left[{\left(E{R}_{ijkl}{u}_{k,l}\right)}_{,j}+{\mathrm{\rho \omega }}^2{u}_i^{}\right],\\ K_0=\frac{1}{2}\underset{\mathrm{V}}{\int }{\mathrm{\rho u}}_{\mathrm{i}}^{}{\mathrm{u}}_{\mathrm{i}}^{}\mathrm{dV}.\end{array}$  (1.9)

Учитывая независимость вариаций $\mathrm{\delta }E,\mathrm{\delta }{u}_i$, получим следующее условие оптимальности:

$2\mathrm{\Pi }=R_{ijkl}{u}_{i,j}{u}_{k,l}={\mathrm{\Lambda }}_{*},$ ${\mathrm{\Lambda }}_{*}=-2{\mathrm{\Lambda {\rm K}}}_0{V}^{-1}.$ (1.10)

Дополняя (1.10) уравнениями движения и условием нормировки (1.5), получаем полную систему:

 ${\left(E{R}_{ijkl}{u}_{k,l}\right)}_{,j}^{}+{\mathrm{\rho \omega }}^2{u}_i^{}=\mathrm{0,} i=\mathrm{1,}\mathrm{2,}\mathrm{3,}$  $\frac{1}{V}\underset{V}{\int }E\left(x\right)dV=m_0.$ 

Таким образом, задача об оптимальном законе изменения модуля E(x) сводится к нелинейной краевой задаче относительно E, u_i, которая может быть решена только численно, например итерационным способом, подобно изложенному в работе [6].

2. Условие оптимальности в изотропном случае для параметров Ламе. В рамках такого же подхода можно получить условие оптимальности в изотропном случае для двух функций – параметров Ламе $\lambda \left(x\right)$ и $\mu \left(x\right)$, характеризующих упругие свойства неоднородного материала. Для этого можно использовать полученный выше функционал (1.4), причем для изотропного случая:

$\text{Π}=\frac{1}{2}\underset{V}{\int }\left(\mathrm{\lambda }{\left(u_{k,k}\right)}^2+\mathrm{\mu }(u_{i,j}+u_{j,i})(u_{i,j}+u_{j,i})\right)dV,$

а условие оптимальности будем отыскивать при дополнительных условиях:

$J_1=\frac{1}{V}\underset{V}{\int }\mathrm{\lambda }dV-{\mathrm{\lambda }}_0=\mathrm{0,}$ $J_2=\frac{1}{V}\underset{V}{\int }\mathrm{\mu }dV-{\mathrm{\mu }}_0=0.$ 

В этом случае сформируем функционал:

$J_{*}=J+{\mathrm{\Lambda }}_1{J}_1+{\mathrm{\Lambda }}_2{J}_2,$

где ${\mathrm{\Lambda }}_{1,}{\text{ Λ}}_2$ – множители Лагранжа, и из условия $\mathrm{\delta }{J}_{*}=0$, приравнивая нулю коэффициенты при независимых вариациях, находим следующие условия оптимальности:

$\begin{array}{l}{\left(u_{k,k}\right)}^2+{\mathrm{\Lambda }}_1{\left(2K_{0}V\right)}^{-1}=\mathrm{0,}\\ \left(u_{i,j}+u_{j,i}\right)\left(u_{i,j}+u_{j,i}\right)+{\mathrm{\Lambda }}_2{\left(2K_{0}V\right)}^{-1}=\mathrm{0,}\end{array}$  (2.1)

${\left(\mathrm{\lambda }{u}_{k,k}\right)}_{,j}+{\left(\mathrm{\mu }\left(u_{i,j}+u_{j,i}\right)\right)}_{,j}+{\mathrm{\rho \omega }}^2{u}_i^{}=\mathrm{0,}i=\mathrm{1,}\mathrm{2,}3.$ (2.2)

Таким образом, имеем:

 $\begin{array}{l}{\left(u_{k,k}\right)}^2={\mathrm{\lambda }}_1^2=\text{const},\\ (u_{i,j}+u_{j,i})(u_{i,j}+u_{j,i})={\mathrm{\lambda }}_2^2=\text{const}.\end{array}$ (2.3)

Эти условия характеризуют постоянство некоторых инвариантов тензора деформаций.

Из них можно определить u_i, j (один из простейших вариантов u_i = a_ijx_j + b_i, a_ij, b_i – постоянные), а затем из соотношений (2.2) получить систему дифференциальных уравнений в частных производных первого порядка для нахождения функций $\mathrm{\lambda },\mathrm{\mu }$:

${\left(\mathrm{\lambda }{a}_{kk}\right)}_{,i}+{\left(\mathrm{\mu }\left(a_{ik}+a_{ki}\right)\right)}_{,k}+{\mathrm{\rho \omega }}^2\left(a_{ik}{x}_k+b_i\right)=\mathrm{0,}i=\mathrm{1,}\mathrm{2,}3.$ (2.4)

Возвращаясь к условию (1.10), отметим, что оно содержит выражение для условной потенциальной энергии деформаций. Для одномерных задач (изгиб, растяжение, кручение стержня) это выражение может быть упрощено, откуда сразу получается условие на соответствующую деформационную характеристику, позволяющее достаточно просто получить аналитическое решение задачи. Так, например, для стержня V = S × [0, l] (S – поперечное сечение) в случае продольных колебаний (заложены гипотезы u₁ = u₂ = 0, u₃ = u(x), x = x₃) имеем $2\mathrm{\Pi }=F{\left(u\text{'}\right)}^2$, где u – характеризует продольное смещение, F – площадь поперечного сечения. Далее, находя u из условия ${\left(u\text{'}\right)}^2={\lambda }_{1*}^2$, получаем уравнение первого порядка относительно E(x). Наряду с гладким решением u(x) = x существует негладкое решение с изломом $u\left(x\right)=\left\{x,x<x_0,2x_0-x,x\ge x_0\right\}.$ Как показали расчеты для случая жесткой заделки на левом конце первое (гладкое) решение отвечает первой собственной частоте, а второе (негладкое) решение отвечает второй собственной частоте.

В случае изгибных колебаний (заложены гипотезы модели Эйлера–Бернулли) $u_1=-x_{2}w\text{'}\left(x\right),u_2=\mathrm{0,}{u}_3=w\left(x\right),x=x_3$ имеем $2\mathrm{\Pi }=J{\left(w\text{'}\text{'}\right)}^2$, где J – момент инерции, w – прогиб нейтральной оси стержня, и условие оптимальности ${\left(w\text{'}\text{'}\right)}^2={\mathrm{\lambda }}_{2*}^2$ позволяет сразу определить прогиб и получить дифференциальное уравнение второго порядка для E(x).

В качестве частных случаев рассмотрим задачи о продольных и изгибных колебаниях стержня, с различными типами граничных условий. Будем считать для простоты, что плотность постоянна $\mathrm{\rho }={\mathrm{\rho }}_0$.

3. Продольные колебания. Пример 1. Обезразмерив задачу о продольных колебаниях неоднородного стержня, получим:

${\left(g{u}^{\text{'}}\right)}^{\text{'}}+\mathrm{\kappa }u=\mathrm{0,}$ $u\left(0\right)=\mathrm{0,}$ $g\left(1\right)u\text{'}\left(1\right)=\mathrm{0,}$ $x\in \left[\mathrm{0,1}\right],$ (3.1)

где g – безразмерный модуль Юнга, k – безразмерный спектральный параметр, пропорциональный квадрату частоты колебаний, u = u(x) — в данном случае безразмерная функция, описывающая продольные смещения. Будем искать такую функцию g = g(x), которая удовлетворяет условию

$\underset{0}{\overset{1}{\int }}g\left(x\right)dx=1$ (3.2)

и для которой первое собственное значение k примет максимальное значение. Используя условие оптимальности, которое имеет вид ${\left(u\text{'}\right)}^2=\mathrm{const}$, с точностью до амплитудного множителя получаем u(x) = x. Решая дифференциальное уравнение (3.1), получаем:

$g=\frac{3}{2}(1-x^2).$

На основе соотношения Релея или на основе численных расчетов получаем k = 3. Заметим, что первое собственное значение для однородной балки $\mathrm{k}=0.25{\mathrm{\pi }}^2\approx 2.467$ и выигрыш составляет порядка 20%.

Найдено также негладкое решение задачи:

$u\left(x\right)=\left\{\begin{array}{l}x,x<1/\mathrm{3,}\\ 2/3-x,x\ge 1/\mathrm{3,}\end{array}\right.$

которому отвечает функция

$g=27\left\{\begin{array}{l}\left(1-9x^2\right)/\mathrm{18,}x<1/\mathrm{3,}\\ \left(4x-3x^2-1\right)/\mathrm{6,}x\ge 1/3\end{array}\right.$

и собственное значение k = 27, которое является вторым.

4. Изгибные колебания. Пример 2. Рассмотрим обезразмеренное уравнение изгибных колебаний неоднородной балки:

$\left(g\left(x\right)w\text{'}\text{'}\left(x,\kappa \right)\right)\text{'}\text{'}-\kappa w\left(x,\kappa \right)=\mathrm{0,}$ $x\in \left[\mathrm{0,}1\right],$  (4.1)

где $g\left(x\right)=E\left(xl\right)E_0^{-1},$ $\mathrm{k}={\mathrm{\rho }}_{0}F{l}^4{\mathrm{\omega }}^2{E}_0^{-1}{J}^{-1}$, E₀ – характерное значение модуля упругости (например, среднее значение). Будем искать такую функцию g = g(x), которая удовлетворяет условию (3.2) и для которой первое собственное значение k принимает максимальное значение. При этом можно рассмотреть следующие типы граничных условий:

1. шарнирное опирание
   ${\left.w\right|}_{x=0}={\left.w\right|}_{x=1}={\left.\left(g{w}^{\text{'}\text{'}}\right)\right|}_{x=0}={\left.\left(g{w}^{\text{'}\text{'}}\right)\right|}_{x=1}=0;$ (4.2)
2. консоль
   ${\left.w\right|}_{x=0}={\left.w^{\text{'}}\right|}_{x=0}={\left.\left(g{w}^{\text{'}\text{'}}\right)\right|}_{x=1}={\left.{\left(g{w}^{\text{'}\text{'}}\right)}^{\text{'}}\right|}_{x=1}=0;$ (4.3)
3. заделка обоих концов
   ${\left.w\right|}_{x=0}={\left.w^{\text{'}}\right|}_{x=0}={\left.w\right|}_{x=1}={\left.w^{\text{'}}\right|}_{x=1}=0;$ (4.4)
4. консоль при наличии упругой заделки на конце x = 0
   

   ${\left.\left({\left(g{w}^{\text{'}\text{'}}\right)}^{\text{'}}+c_{1}w\right)\right|}_{x=0}={\left.\left(g{w}^{\text{'}\text{'}}-c_2{w}^{\text{'}}\right)\right|}_{x=0}={\left.\left(g{w}^{\text{'}\text{'}}\right)\right|}_{x=1}={\left.{\left(g{w}^{\text{'}\text{'}}\right)}^{\text{'}}\right|}_{x=1}=0.$(4.5)

Здесь c₁, c₂ – положительные параметры, характеризующие характеристики упругости в заделке [16].

Были построены оптимальные решения для граничных условий (4.2)–(4.4), которые совпали с рассмотренными ранее в работе [9]. Заметим, что при условиях (4.2) и (4.3) существуют, соответственно, гладкие решения в виде полиномов 4-го порядка:

$w\left(x\right)=x\left(1-x\right),g=5\left(x^4-2x^3+x\right),$

$w\left(x\right)=x^2,g=20\mathrm{\phi }\left(x\right),\mathrm{\phi }\left(x\right)=\left(\frac{x^4}{24}-\frac{x}{6}+\frac{1}{8}\right),$  (4.6)

а при условиях (4.4) функция w(x) имеет два разрыва второй производной, а g(x) – разрывы первой производной. Эти решения можно представить в виде:

$w\left(x\right)=\left\{x^2\text{},x<\mathrm{0.25,}0.125-{\left(0.5-x\right)}^2\text{},0.25\le x\le \mathrm{0.75,}{\left(1-x\right)}^2\text{},x>0.75\right.;$

можно записать с учетом симметрии относительно середины отрезка:

$g\left(x\right)=\left\{\begin{array}{l}{g}_1\left(x\right),x<\mathrm{0.5,}\\ g_1\left(1-x\right),x\ge \mathrm{0.5,}\end{array}\right.$

$g_1\left(x\right)=\frac{15360}{23}\left\{\begin{array}{l}\frac{23}{6144}+\frac{x^4}{24}-\frac{x}{64},x<\frac{1}{4},\\ \frac{x^4}{24}-\frac{x^3}{12}+\frac{x^2}{32}+\frac{x}{96}-\frac{7}{2048},\frac{1}{4}\le x\le \frac{1}{2}.\end{array}\right.$(4.7)

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

Таблица 1. Первое собственное значение для трех вариантов граничных условий

Оказалось, что по сравнению с однородным случаем первое собственное значение k удается увеличить соответственно на 23, 61 и 33%.

Обратимся теперь к задаче с граничными условиями (4.5). Получим соотношение Релея. Для этого умножим уравнение (4.1) на w и проинтегрируем по отрезку [0, 1]:

$\underset{0}{\overset{1}{\int }}\left({\left(g{w}^{\text{'}\text{'}}\right)}^{\text{'}\text{}\text{'}}w-\kappa w^2\right)dx=0.$

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

$\begin{array}{c}\kappa =\frac{{\displaystyle \underset{0}{\overset{1}{\int }}{\left(g{w}^{\text{'}\text{'}}\right)}^{\text{'}\text{}\text{'}}wdx}}{{\displaystyle \underset{0}{\overset{1}{\int }}{w}^{2}dx}}=\frac{{\left.{\left(g{w}^{\text{'}\text{'}}\right)}^{\text{'}}w\right|}_0^1-{\left.\left(g{w}^{\text{'}\text{'}}\right)w^{\text{'}}\right|}_0^1+{\displaystyle \underset{0}{\overset{1}{\int }}g{\left(w^{\text{'}\text{'}}\right)}^{2}dx}}{{\displaystyle \underset{0}{\overset{1}{\int }}{w}^{2}dx}}=\\ =\frac{c_1{w}^2\left(0\right)+c_2{\left(w^{\text{'}}\right)}^2\left(0\right)+{\displaystyle \underset{0}{\overset{1}{\int }}g{\left(w^{\text{'}\text{'}}\right)}^{2}dx}}{{\displaystyle \underset{0}{\overset{1}{\int }}{w}^{2}dx}}.\end{array}$

Условие оптимальности принимает тот же вид ${\left(w\text{'}\text{'}\right)}^2=$const. Решение с точностью до амплитудного множителя будем искать в виде w(x) = x² + kx + b. Заметим, что в зависимости от граничных условий решение задачи может иметь или не иметь разрывов второй производной, что демонстрируют случаи (4.6) и (4.7). Следовательно, не для всех значений параметров c₁, c₂ гладкое решение вида w(x) = x² + kx + b существует.

Решая (4.1) относительно g(x), находим:

$g=\frac{\mathrm{\kappa }}{2}\left(\frac{x^4}{12}+k\frac{x^3}{6}+b\frac{x^2}{2}\right)+mx+p,$ (4.8)

где $k=\mathrm{\kappa }\left(3c_1-\mathrm{\kappa }\right)/\mathrm{\Delta },$ $b=\mathrm{\kappa }\left(\mathrm{\kappa }+24c_2\right)/\left(6\mathrm{\Delta }\right),$ $m=-c_{1}b/2,p=c_{2}k/2,$ $\mathrm{\Delta }={\mathrm{\kappa }}^2-4\mathrm{\kappa }{c}_1-12\mathrm{\kappa }{c}_2+12c_2{c}_1.$

В результате решение запишется в виде:

$\begin{array}{l}w\left(x\right)=x^2+\mathrm{\kappa }x\frac{\left(3{\mathrm{c}}_1-\mathrm{\kappa }\right)}{{\mathrm{\Delta }}_1{\mathrm{\Delta }}_2-3{\mathrm{\kappa }}^2}+\mathrm{\kappa }\frac{1}{6}\frac{\left(24c_2+\mathrm{\kappa }\right)}{{\mathrm{\Delta }}_1{\mathrm{\Delta }}_2-3{\mathrm{\kappa }}^2},\\ g=\mathrm{\kappa }\left(\frac{x^4}{24}+\frac{\left(3c_1-\mathrm{\kappa }\right)}{2\mathrm{\Delta }}\left(\mathrm{\kappa }\frac{x^3}{6}+c_2\right)+\frac{\left(\mathrm{\kappa }+24c_2\right)}{12\mathrm{\Delta }}\left(\mathrm{\kappa }\frac{x^2}{2}-c_{1}x\right)\right).\end{array}$

Параметр k определяется из условия (3.2). Возникающее кубическое уравнение

$\frac{\kappa }{720}\frac{\left({\kappa }^2-9\kappa c_1-192\kappa c_2+432c_1{c}_2\right)}{\left({\kappa }^2-4\kappa c_1-12\kappa c_2+12c_1{c}_2\right)}=1$

решается численно. Для анализа полученного результата рассмотрим 2 частных случая:

Случай 1 $\left(c_2\to \infty \right)$: ${\left.({\left(g{w}^{\text{'}\text{'}}\right)}^{\text{'}}+c_{1}w)\right|}_{x=0}=\mathrm{0,}{w}^{\text{'}}\left(0\right)=0.$

Параметры в (4.8) принимают вид k = 0, b = k/3${∆}_1$, m = -c₁k/6${∆}_1$, p = k(3c₁ – k)/24${∆}_1$, _{${∆}_1$}= c₁ – k, что позволяет записать решение оптимизационной задачи в виде:

$\begin{array}{l}w\left(x\right)=x^2+\kappa /\left(3c_1-3\kappa \right),\\ g=\kappa \left(\frac{x^4}{24}+\frac{\kappa x^2}{12{\Delta }_1}-\frac{c_{1}x}{6{\Delta }_1}+\frac{\left(3c_1-\kappa \right)}{24{\Delta }_1}\right).\end{array}$ (4.9)

Параметр k определим из условия (3.2), которое сводится теперь к квадратному уравнению, имеющему решение следующего вида:

 $\kappa =\frac{9}{8}{c}_1+\frac{45}{2}-\frac{3}{8}\sqrt{9c_1^2+40c_1+3600}.$  (4.10)

Случай 2 $\left(c_2\to \infty \right)$: w(0) = 0, ${\left.\left(g{w}^{\text{'}\text{'}}-c_2{w}^{\text{'}}\right)\right|}_{x=0}=0.$

Аналогично имеем:

$b=\mathrm{0,}k=3\kappa /\left(4{\Delta }_2\right),m=-\kappa \left(24c_2+\kappa \right)/\left(48{\Delta }_2\right),p=3\kappa c_2/\left(8{\Delta }_2\right),{\Delta }_2=3c_2-\kappa ,$

$w\left(x\right)=x^2+3x\kappa /\left(12c_2-4\kappa \right),g=\kappa \left(\frac{x^4}{24}+\frac{\kappa x^3}{16{\Delta }_2}-\frac{\left(24c_2+\kappa \right)x}{48{\Delta }_2}+\frac{3c_2}{8{\Delta }_2}\right),$ (4.11)

$\kappa =160+24c_2-8\sqrt{9c_2^2+105c_2+400}.$ (4.12)

Асимптотические представления решений (4.9), (4.11) при больших c₁, c₂ имеют вид:

$g=\kappa \phi \left(x\right)+{\kappa }^2\left(\frac{x^2}{12}-\frac{x}{6}+\frac{1}{12}\right)\frac{1}{c_1}+o\left(\frac{1}{c_1}\right)$ (Случай 1), (4.13)

$g=\kappa \phi \left(x\right)+{\kappa }^2\left(\frac{x^3}{48}-\frac{x}{16}+\frac{1}{24}\right)\frac{1}{c_2}+o\left(\frac{1}{c_2}\right)$ (Случай 2). (4.14)

На основе приведенных асимптотических представлений параметр k также может быть определен как один из корней квадратного уравнения, возникающего при удовлетворении условию (3.2):

$\kappa =-\frac{9}{10}{c}_1+\frac{3}{10}\sqrt{9c_1^2+400c_1},$ (4.15)

$\kappa =-\frac{8}{5}{c}_2+\frac{8}{5}\sqrt{c_2^2+25c_2}.$ (4.16)

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

$\kappa =\frac{40}{1+\sqrt{1+{\left(\frac{15}{3}\right)}^2\frac{1}{c_2}+{\left(\frac{20}{3}\right)}^2\frac{1}{c_1}}}.$ (4.17)

5. Результаты численных расчетов. Выполним расчеты первого собственного значения для однородной и оптимальной балок. В случае конечных значений c₁, c₂ будем находить k численно и при помощи формулы (4.17). В случае, когда $c_2\to \infty$, будем сравнивать точное k, найденное из (4.10), и асимптотическое из (4.15). Аналогично сравним (4.12) и (4.16) для случая $c_2\to \infty$.

Замечание 2. У возникающих квадратного и кубического уравнений относительно k соответственно 2 и 3 корня. В (4.10), (4.12), (4.17) и в табл. 2 указан только один из них, поскольку остальные приводят к знакопеременной функции g(x); соответственно, в (4.15), (4.16) указаны положительные корни, а отрицательные опущены.

Таблица 2. Первое собственное значение для различных параметров закрепления

Рис. 1. Собственное значение k: по (4.10) – сплошная линия, по (4.15) – пунктир

Рис. 2. Собственное значение k: по (4.12) – сплошная линия, по (4.16) – пунктир

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

Работа выполнена за счет гранта РНФ № 22-11-00265, https://rscf.ru/project/22-11-00265/ в Южном федеральном университете

Об авторах

А. О. Ватульян

Автор, ответственный за переписку.
Email: aovatulyan@sfedu.ru
Россия, Ростов-на-Дону

В. О. Юров

Email: vitja.jurov@yandex.ru
Россия, Владикавказ

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

Дополнительные файлы

Доп. файлы

Действие

1. JATS XML

Скачать

2. Рис. 1. Собственное значение k: по (4.10) – сплошная линия, по (4.15) – пунктир

Скачать (45KB)

Метаданные

3. Рис. 2. Собственное значение k: по (4.12) – сплошная линия, по (4.16) – пунктир

Скачать (44KB)

Метаданные