Научная статья на тему 'К вопросу динамического анализа статически неопределимых систем'

К вопросу динамического анализа статически неопределимых систем Текст научной статьи по специальности «Строительство и архитектура»

CC BY
61
20
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Кочегаров И.И., Жихарев К.В., Паршиков А.А., Таньков Г.В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «К вопросу динамического анализа статически неопределимых систем»

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

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

В соответствии с подходом Фойхта приведенный модуль упругость определяется выражением

п

еф = у • Е, (1)

у=1

где ^^ - относительное объемное содержание; Е - модуль упругости у-го компонента в структуре.

Применительно к каждому к -му слою гетерогенной структуры его модуль упругости в соответствии с (1) определяется выражением 7(к) _

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

• ^ • (3)

При использовании такого подхода в ранее предложенной модели при определении приведенного модуля и коэффициента ^ следует полагать:

E — E(3) . E — e(4) . G — G

Е3 — Епр ' e4 _ Епр ' G Gm

Епр ' Е4 Епр

Существенной особенностью при изменении рисунка

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

плате изменяются значения

Y

при

k — 2,3,4 ),

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

Следует заметить, что если защитный слой заполняет все «выемки» «А» (см. рис. 1), то это учитывается следующим образом: приведенные модули в слоях с к = 2 3 4 вычисляются с учетом

ЕПр ■ Ек ,

(2*

V

m ' k

Здесь V - объем k-

го слоя до

травления; V - объем к-го слоя после травления к

( V < V )• Следовательно Е(1) = Ел - для к к пР 1

основания

заполнения выемок «А»

Е(к) = Т • Ек + (1 )Е5 •

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

платы и Ejp — Е5 -

для защитного слоя.

ЛИТЕРАТУРА

1. Кочегаров И.И. Виды дефектов при распознавании повреждений печатных плат / И.И. Кочегаров, Е.А. Данилова // Современные наукоемкие технологии. 2014. № 5-1. С. 60-61.

2. Гришко А.К. Алгоритм верификации электромагнитной устойчивости радиоэлектронных плат / А.К. Гришко, И.И. Кочегаров, Е.С. Каракулов // Инновации на основе информационных и коммуникационных технологий. 2015. Т. 1. С. 301-304.

3. Кочегаров И.И. Межсистемное взаимодействие систем CAD и CAE при моделировании / И.И. Кочегаров, Н.К. Юрков, В.Б. Алмаметов // XXI век: итоги прошлого и проблемы настоящего плюс. 2015. № 4 (26). С. 161-166.

4. Кочегаров И.И. Выбор оптимального варианта построения электронных средств / И.И. Кочегаров, Н.В. Горячев, А.К. Гришко // Вестник Пензенского государственного университета. 2015. № 2 (10). С. 153-159.

5. Белов А.Г. Обеспечение влагозащитного покрытия печатных узлов датчика протечки / Белов А.Г., Баннов В.Я., Трусов В.А., Кочегаров И.И., Лысенко А.В., Горячев Н.В., Юрков Н.К. // Труды международного симпозиума Надежность и качество. 2014. Т. 2. С. 151-154.

6. Реута Н.С. Разновидности механических воздействий в радиоэлектронной аппаратуре / Реута Н.С., Горячев Н.В., Трусов В.А. // Молодой ученый. 2014. № 21 (80). С. 224-226.

7. Трусов В.А. Использование самоорганизующегося алгоритма для нелинейных технологических процессов / Трусов В.А. // Труды международного симпозиума Надежность и качество. 2012. Т. 2. С. 395.

УДК 004.94

Кочегаров И.И., Жихарев К. В., Паршиков А.А. Таньков Г.В.

К ВОПРОСУ ДИНАМИЧЕСКОГО АНАЛИЗА СТАТИЧЕСКИ НЕОПРЕДЕЛИМЫХ СИСТЕМ

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

В статье [1] на основе метода конечных разностей (МКР) разработан алгоритм расчета собственных частот и форм изгибных колебаний стержня. На основе этого алгоритма написана программа по расчету первой собственной формы и частоты изгибных колебаний стержня, где использован итерационный метод, который состоит в последовательном уточнении значений перемещений в узлах модели при резонансе и собственной частоты.

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

Проверка точности численного решения проводилась сравнением аналитического и численного решений.

Так для стержня сечением 3x6 мм, шарнирно закрепленного по концам (рис. 1), значение первой резонансной частоты определялось по формуле [2]:

где ш-1 = -

Д = <м1/2п,

Величины, входящие в формулу: Е =

I2

2,1 X 1011 Н/м2 - модуль упругости материала; / = 54 X 10-12м4 - момент инерции сечения относительно оси у; р = 7800 кг/м3 - плотность материала; 5 = 18 X 10-6 м2 - площадь стержня; I - длина стержня, м; 1 = 1 - волновое число.

Рисунок 1

длины I = 1 м.

Аналитическое решение дает для /1 = 14,11 Гц, для ¿ = 2 м,/1 = 3,53 Гц.

Численное решение дало для I = 1 м.Д = 14 Гц, для I = 2 м,/1 = 3,0 Гц.

Для этого же стержня, жестко закрепленного с одного края и шарнирно закрепленного на другом конце (рис. 2), использована аналитическая фор-

3,932 /Ёу I2 ^^

мула

Рисунок 2

Для 5,52 Гц.

для I = 2 м,/1 = для I = 1 м,/1 =

при длине пролета ¿>5 м. не превосходило

при длине пролета

(цитируется по

загруженных заданными пролетными нагрузками [3], то величину прогиба каждого пролета от действия статической нагрузки можно найти по формулам простых однопролетных статически определимых балок [3].

Выберем сначала профиль балки. Наибольший изгибающий момент здесь будет во втором пролете

(г = 10 м.).

ц12 о,бхю2

мЛ5

8

8

■ =7,5 тхм.

о)

Рг-Зг

А' Г^Н.

иТГПР I ПтТТП, ГГТГТп 1111/1 ТТТ777Г

.;. ю .1 * .|'Й1<ма1

длины I = 1 м,/1 = 22,1 Гц, а

Численное решение дает 22 Гц, а для ¿ = 2 м,/1 = 5,0 Гц.

Результаты расчетов показали достаточно высокую точность численного решения, поскольку расхождение численного решения с аналитическим составляет доли процента, поэтому предложенная модель и программа [1] могут быть использованы для анализа реальных конструкций.

В данной статье делается попытка применения указанной программы к расчету динамики статически неопределимых стержневых конструкций. Поскольку в конструировании РЭС подобных расчетов не проводилось, обратимся к опыту строительной механики и сопромата, где «статически неопределимая балка, имеющая более двух опор, называется неразрезной балкой. Такие балки находят широкое применение в различных конструкциях.» (цитируется по [3]). Но в балках большой длины, «подверженных переменным нагрузкам, упругие перемещения могут оказывать вредное влияние. Напряжения, возникающие при колебаниях как следствие больших перемещений, особенно в случае резонанса, могут возрасти до разрушающих значений. Поэтому балки должны удовлетворять условию прочности и условию жесткости. Например, нормы на проектирование стальных конструкций требуют, чтобы абсолютное значение наибольшего прогиба

а

зии

¿ = 7 м. составляло не более [4]).

В строительной механике статически неопределимые балки после определения опорных моментов в промежуточных опорах считают с использованием уравнения трех моментов [3].

Для динамического анализа возьмем из [3] пример 9.4. (стр. 262), где для балки длиной 28 м. с тремя промежуточными опорами рассчитаны и построены эпюры перерезывающих сил и изгибающих моментов (от опорных моментов и моментов от нагрузки) - рис. 3а, б.

По эпюре моментов (рис. 3, б) создается представление о характере изогнутой оси балки (рис. 3, в) [3].

Так как рассматриваемая неразрезная балка эк-Бивалентна системе четырех однопролетных балок,

1' 1

Эпн>роМ] , I ! |

—I 1 I I

Рисунок 3

Необходимый момент сопротивления поперечного сечения [4]

ш = "ь™* = 750000 кг/см = 535,7см3 ^ 1400 кг/см2

где [о"в] - допускаемое напряжение (здесь взято 1400 кг/см2.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Выбираем двутавр №30а [3]. Его характеристики - ширина Ь=126 мм, высота 11=300 мм; площадь сечения S = 61,2 см2; момент инерции /2 = 8950см4; момент сопротивления поперечного стержня = 597см3.

Прогиб первого пролета от нагрузки Р1 = 3Г находим по формуле [4]

1 Аха2хЬ2 1х 3000 22 х 42

А1 = ^х"

х

■ г 6 мм,

3г 18795 х 103 3 х6

где Е = 2,1 х 1011 Н/м2 - модуль упругости.

Прогиб второго пролета от распределенной

нагрузки определим по формуле [4]

5г3 5 х 103

Д2 = ^ ™-= 0,6 х 10——^т———т г 41,4 мм.

2 384 х Я/ 384 х 18795 х 103

Прогиб четвертого пролета от двух сил Р2 и Р3

будет равен [4]

Дз =

¿7[1рй(3г2-4й2)]

1

[2 х 103

48 х 18795 х 103 х 2(3 х 82 - 4 х 22) + 2 х 103 х 2(3 х 82 - 4 х 22)] г 15 мм.

В масштабе прогибов изогнутая ось балки будет иметь вид (рис. 4):

Рисунок 4

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

Предлагаемая программа позволяет получить первую собственную форму такой неразрезной

балки. Введем в качестве исходных данных параметры двутавра №30а и длину балки ¿ = 28 м. Граничные условия: на левом конце жесткое защемление, на правом - жесткий шарнир. Три промежуточные опоры - жесткие шарниры. Ускорение внешнего воздействия - 1д. Результат расчета дан на рис. 5.

Рисунок 5

Рисунок 6

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

Но амплитуда прогибов уже при 1д внешнего воздействия во втором пролете равна 575,9 мм (в статике - 41 мм). Изгибающий момент в среднем

сечении второго пролета можно определить по разностной формуле

д2^ Шх + Ю- 2w(x) + = ^ = -^-).

Шаг сетки в данном случае равен /гх = 280 мм. Значение амплитуд смещения в середине второго пролета равны: ^(х+я) = 575,9 мм, М(Х) = 570,6 мм, = 559,9 мм. Тогда МЬтах = 129,4 Г хм. При моменте

сопротивления = 597 см3, максимальное напряжение в сечении будет

12940000 кг/см

^Ьтох = -^-5- = 21675 кг/см2 » 1400 кг/см2

0т 597 см3

Момент сопротивления в этом случае будет равен Щ, = 8032 см3 » 597 см3.

Такая нагрузка для данной балки неприемлема по условию прочности.

Снизим амплитуду внешнего воздействия до 0,1д. Результат расчета - на рис. 6.

Здесь амплитуда формы приближается к статическим прогибам, но во втором пролете еще высока и для этой амплитуды момент сопротивления = 924 см3 > 597 см3, следовательно и изгибное напряжение будет больше допустимого.

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

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

Для анализа из [3] взят пример 9.6 (стр. 270), где для балки с тремя промежуточными опорами построена эпюра моментов и изогнутая ось балки от осадки опоры 2 (рис. 7)

Рисунок 7

Для примера возьмем длину балки 4/ = 8 м. и профиль сечения - двутавр №30а. Крепление балки в точке 1 - жесткое защемление, в точке 2, 3, 4, 5 - жесткий шарнир. Введем эти данные в программу и зададим внешнее возмущение амплитудой 1д с учетом подвижности опоры 2 в пролете 1-3 (рис. 7). Результат решения дан на рис. 8.

I

Рисунок 8

При амплитуде внешнего воздействия в 2д и 3д амплитуда прогиба пролета 1-3 в собственной форме будет соответственно 17,1 мм. и 26,6 мм., но момент сопротивления сечения под опорой 2 будет равен 154 см3 и 257 см3, что меньше ^ = 597 см3.

При внешнем воздействии 10д амплитуда прогиба пролета 1-3 возрастает до 85,4 мм. В этом случае изгибающий момент в сечении над опорой 2 будет равен МЬтах = 12,7 Г/м, а напряжение изгиба

1270000кгхсм 597 см3

= 2128 кг/см2,

что больше допустимого напряжения 1400 кг/см2 в 1,5 раза.

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

Статья подготовлена в рамках реализации проекта «Разработка методов и средств создания высоконадежных компонентов и систем бортовой радиоэлектронной аппаратуры ракетно-космической и транспортной техники нового поколения» (соглашение №15-19-10037 от 20 мая 2015 г.) при финансовой поддержке научного фонда.

ЛИТЕРАТУРА

1. Таньков Г.В. Разработка и исследование дискретных моделей несущих конструкций РЭС / Г.В. Таньков, В.А. Трусов, И.И. Кочегаров - Прикаспийский журнал: Управление и высокие технологии. -

№2, 2 016 г.

2. Бабаков И.М. Теория колебаний. - М.: Наука, 1965 - 560 с.

3. Рудицын М.Н. Справочное пособие по сопротивлению материалов / М.Н. Рудицын, П.Я. Артемов, М.И. Любошиц. - Изд-во «Высшая школа», Минск, 1970 - 628 с.

4. Дрейер Г. Учение о прочности. - М.: Машиностроение, 1964 - 416 с.

ь

i Надоели баннеры? Вы всегда можете отключить рекламу.