Оптимизация несущих конструкций блоков РЭС Текст научной статьи по специальности «Физика»

Селиванов В. Ф. Таньков Г. В.

ОПТИМИЗАЦИЯ НЕСУЩИХ КОНСТРУКЦИЙ БЛОКОВ РЭС

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

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

Физическая модель представляется в виде дискретной модели - сетки, расчетная схема является матема-

тическим описанием процессов в модели - сетке.

Модель - сетка сложной конструкции представляет собой сочетание дискретных элементов прямоугольной формы с размерами hx, hy, ^ в направлении трех координатных осей. От размеров элементов зависит точность расчета, поэтому относительные размеры элементов должны быть достаточно малыми. С другой стороны завышение числа элементов в модели сверх требуемого из условия точности приводит к большим затратам времени на вычисление и снижает надежность расчета.

Каждый элемент в модели - сетке находится под действием сил, таковыми являются силы упругости,

действующие на стороны соседних элементов, силы инерции и силы тяжести.

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

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

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

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

Оборудования и устройства учитываются в модели в виде сосредоточенных масс, добавленных к массе

соответствующих узлов модели сетки.

Рис. 1.

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

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

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

у р, = т

у р т э2

(1.1)

где Е^- силы упругости примыкающих стержней; 1, j - номер данного и соседнего узла соответственно; Ш1 - масса узла (элемента), т. е. масса входящих в элемент половинок стержней с возможными дополнениями в виде других навесных масс; «1 - перемещение узла по оси z.

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

Еч=оч^-01), (1.2)

где с^ - жесткость стержня связи между узлами 1, j.

При перемещении узлов модели по оси z вертикальные стержни будут растягиваться, горизонтальные изгибаться, а наклонные и растягиваться и изгибаться. Чтобы вывести общее выражение упругих сил для

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

В і - ом узле стержня введем локальную систему координат, направив ось х по оси стержня, а ось г' - перпендикулярно к ней.

Рис. 2.

Воздействуя на стержень в і - том узле растягивающей силой Ех'х' и изгибающей силой Г2'х , выразив их согласно (1.2) через перемещения узлов в локальной системе:

Ех'х'—Ср(иі — иі ) (1.3)

Р^'х'_Си("^ — )

где Ср и Си - жесткость стержня на растяжение и изгиб соответственно. Величина каждой жесткости определяется по формулам:

С' = Т ’

У

12 К7„

cu ' т3

1 У

Поскольку расчет осуществляется в системе координат xyz , спроектируем действующие силы из локальной системы в принятую систему координат (рис. 2)

Fzz=Fx'x' sin а (1.4)

Fzx=Fz'x' cos а

Также выразим перемещения локальной системы через перемещения узлов модели в системе xyz:

Ui'=Wi sin а; (1.9)

Uj'=Wj sin а ;

Wi'= Wi cos а;

Wj'= Wj cos а.

Подставив эти перемещения в (1.3), из (1.4) получим выражение сил, действующих в общем случае на любой узел модели:

Fzz= Cp sin2 а (Wj - Wi); (1.5)

Fzx= Cu cos2 а (Wj - Wi)

С учетом (1.5) уравнение динамического равновесия узла запишется как

(Q

'LicPiisin2 a(mj -Ti) + cUij cos2 a(Tj - T)] = mi

v ^ N -T

или LC0(j)(Tj -T) =m^~^2~ ' (1.6)

где c0(j) = cpiJ sin2 a + c\(1 - sin2 a) (1.7)

Учитывая, что при резонансных колебаниях конструкции на первой частоте все точки ее колеблются синхронно и синфазно, по гармоническому закону T = WsinT/ , (1.8)

где W- амплитуда колебаний узла; о - собственная частота ВЧj)(Wj - W)] = -miT2Wi (1.10)

j

Это и будет исходная система уравнений для расчета низшей резонансной частоты конструкции. Амплитуды Wi (собственная форма колебаний) определяется с точностью до постоянного коэффициента.

Для учета вибраций и ударов необходимо учесть диссипативные силы, и тогда система (1.6) приводится к виду:

с) с2 со

L[c0(j)(Tj-T)]+%%щ[c0(j)(T -Tj)] =m-Tt, (1.11)

где в левой части обозначены суммы сил, действующих на массу узла (Щ - коэффициент вязкости).

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

Тогда (1.6) можно записать

где Е - модуль упругости материала стержня в Н/м ; Б - площадь поперечного сечения стержня в м ;1

длина стержня в м ; ^ - момент инерции сечения стержня относительно оси, перпендикулярной стерж

4

ню, в м

т,(1 + г) - 2т,(I) +1о,( - г)

а = — ----------'-г2--- ---- , (1.12)

г г2

где о - перемещение узла в три момента времени.

Отсюда получим уравнение для расчета перемещений в узлах модели: ^+т); ^ ^-т)

т$ + г) = г2а2 + 2т,.(0 -т-г) . (1.13)

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

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

Щ = АЖС ,

где А - масштабный множитель.

Опуская промежуточные выкладки пишем формулу для определения А |8|

тЩ 2

' " (1.14)

п 2

где Ио - амплитуда внешнего гармонического воздействия.

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

Ж(о2

Ч,=-!—, (1-15)

§

где д- ускорение свободного падения.

Отметим, в заключении, особенности расчета статического напряжения, в этом случае система уравнений имеет вид:

2)(Щ -Щ) =Р,, (1.16)

/

где 1=1,2,...,Ы; Ы-число узлов модели; - сосредоточенная сила в 1-том узле.

Система (1.16) решается итерационным методом, в результате чего по описанным выше формулам определяются силы и напряжения.

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

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

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

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

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

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

[си+Дй-Мй = С

|м(0) = м|-0-|;м( 0) = и1-0-1

С, Д, М - соответственно матрицы: жесткости, демпфирования, инерции.

В общем виде задача оптимального проектирования формируется следующим образом: штФ(^,Ъ,ы)

Ь <е В <^Кп

\си+ Дй-Ми=С

|м(0) = м|-°'|;м( 0) = и1-0-1

где Ф(^, Ъ,ы) - целевая функция (масса), В - множество допустимых значений вектора Ь, Яп - п - мер-

ное действительное пространство, точкой которого является вектор Ь.

\£(* )|= К(-^доп;

И* )| = к*-°доп;

|в(0|=кв-адоп,

где £, о, 0 - соответственно деформация, напряжение, ускорение; К^, Ко, Ко - соответствующие допустимые коэффициенты запаса.

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

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

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

В зависимости от расчетного значения коэффициента запаса прочности

г/* &

Кпр =----- , (^2°)

&доп

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

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

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