CC BY
56
Оптимизация толстостенной сферической оболочки на основе теории прочности Мора
Б.М. Языев, А.С. Чепурненко, А.В. Муханов
Толстостенные оболочки в настоящее время находят все более широкое применение в таких конструкциях, как радиационно-тепловые экраны ядерных реакторов, тепловые и биологические защиты и т.д. Напряжения в таких конструкциях распределяются неравномерно. Исчерпание несущей способности может происходить в небольшой области, т.е. материал конструкций часто используется нерационально.
Одним из способов оптимизации толстостенных оболочек является изменение модуля упругости материала. Идея данного метода принадлежит академику РААСН проф. В.И. Андрееву [1,2]. Сущность метода заключается в том, что если создать искусственную неоднородность, уменьшив модуль упругости там, где возникли наибольшие напряжения, можно добиться такого состояния, когда расчетные напряжения в каждой точке конструкции одинаковы, т.е. конструкция будет равнонапряженной. Если при изменении модуля упругости прочностные характеристики материала не меняются, то конструкция также будет равнопрочной.
Решения задач оптимизации толстостенных цилиндров и сфер, а также других конструкций рассматриваются во многих работах[1-10]. В основном в них авторы определяют эквивалентные напряжения по четырем классическим теориям прочности. В данной статье мы рассмотрим решение задачи оптимизации сферической оболочки при использовании теории прочности Мора.
Пусть имеется толстостенная сферическая оболочка с внутренним радиусом а и внешним радиусом Ь, нагруженная внешним давлением рЬ и внутренним давлениемра (рис.1).
Данная оболочка будет находиться в условиях центральносимметричной задачи теории упругости. Основное разрешающее уравнение относительно напряжений <5г при изменении модуля упругости по толщине имеет вид [2]:
Рис. 1. Расчетная схема оболочки.
// ,4 Е1. / 12(1 - 2v) E1 _ (1)
а" +(---------)а'г------Ь-----а г = 0. (1)
r Е r 1 -v Е
Штрих здесь и далее - дифференцирование по радиусу.
Эквивалентное напряжение по теории прочности Мора для центральносимметричной задачи принимает вид: аэкв = ae-kar, где k = [а ]/[ас ]-отношение допускаемых напряжений на растяжение и сжатие.
Найдем распределение модуля упругости, при котором аэкв = const (2). Условию (2) эквивалентно равенство:
аЭкв =ае - kai = 0. (3)
Напряжения ае можно выразить через аг из дифференциального
уравнения равновесия:
_/ . 2(аг-ае) , К-/
аг +---г----^ = 0 ^ ае = аг + - гаг. (4)
г 2
После подстановки (4) в (3) приходим к следующему
дифференциальному уравнению:
га// =(2^ - 3К. (5)
Общее решение уравнения (5) имеет вид:
г 2к-2
°г=С-21Т^+С2. (6)
Неизвестные константы С1 и С2 найдем из граничных условий аг (а) = -Ра, аг (Ь) = -рЬ. Ограничимся случаем, когда ра = р, рЬ = 0.
С = р(2к - 2) С = -рЬ2к-2
Ь2к-2 _ а2к-2 , 2 Ь2к-2 _ а2к-2. V '
(у Сф (у Сф
Подставим теперь выражение (6) в (1). В результате получим
дифференциальное уравнение с разделяющимися переменными:
йЕ С(2к + 1)г2 к-3 2(1 - 2у)
— =-------------------------, где т = —-------.
Е СГ к-2(1 + -^-) + тС2 1 (8)
1 v 2к - 2 2
Введем следующие обозначения:
А = С,(2к +1), В = С,(1 + -т-г), С = С2т.
2к - 2
Тогда искомая зависимость модуля упругости от радиуса будет иметь
вид:
Е = С0 (Вг2 к-2 + С)В (2к-2) (9)
С0 в формуле (9) - произвольная постоянная. Выражение (9)
справедливо при В ф 0 и к ф 1. Отдельно рассмотрим случай, когда В = 0:
п п т т V
С, (1 +-----------------) = 0 ^ к = 1-=-.
1 2к - 2 2 1 -V
Расчетное напряжение по II теории прочности:
= °е - Ч°г + %) = °е - Ч°г + °е) = К - 7-^—аг)(1 - V) (10)
1 -V
Как видно из выражения(10) расчетное напряжение по II теории прочности при В=0 отличается от расчетного напряжения по теории Мора только множителем (1 -V), то есть законы распределения Е(г) в обоих случаях должны быть одинаковыми. Таким образом, подставив в формулу (9)
V
к«------ , мы получим приближенное решение для II теории прочности.
1 -V
Первая и третья теория прочности являются частными случаями теории Мора соответственно при к = 0 и к = 1.
Проверить правильность решения мы можем, подставив полученную зависимость модуля упругости в исходное дифференциальное уравнение(1) и решив прямую задачу, например, методом конечных разностей.
На рис. 2 и 3 изображены кривые изменения модуля упругости и напряжений, построенные при V = 0,3; Ь / а = 1,5; р = 10МПа; к = 0.5.
Штриховой линией на рис. 3 показано решение для неоднородной сферы, а сплошной - для однородной.
20
15
10
-5
-10
-15
0эхе_
^экв
<7г
\<*г.
1.1
1.2
1.3
1.4
1.5
г/а
Рис.3. График изменения напряжений по толщине сферы.
Рис.2. График зависимости модуля упругости от радиуса (Е0=Е(а)) для неоднородной сферы.
Как видно из графиков, максимальные расчетные напряжения для неоднородной сферы при той же толщине уменьшились в 1.6 раза. График напряжений аэке при переменном модуле упругости представляет собой прямую линию, что свидетельствует о правильности решения обратной задачи.
Литература:
1. Андреев В.И. Потехин И. А. Оптимизация по прочности толстостенных оболочек: монография. М.:МГСУ, 2011. 86с.
2. Андреев В.И. Некоторые задачи и методы механики неоднородных тел: монография. - М.: Издательство АСВ, 2002. - 288 с.
3. Потехин И.А. Способ оптимизации конструкций на основе решения обратных задач теории упругости неоднородных тел. Дисс. канд. техн. наук. М., 2009. - 144 с.
4. Andreev V.I. Minaeva A.S. Creation on the basis of the first theory of strength model equal stressed cylinder exposed to power and temperature loads. International Journal for Computational Civil and Structural Engineering. Volume 7, Issue 1, 2011. p. 71-75
5. Андреев В.И., Булушев С.В. Оптимизация неоднородной толстостенной сферической оболочки, находящейся в температурном поле. Вестник МГСУ, 2012, №12, стр. 40-46.
6. Andreev V.I. Optimization of thick-walled shells based on solutions of inverse problems of the elastic theory for inhomogeneous bodies. Computer Aided Optimum Design in Engineering XII (OPTI XII). WIT Press. 2012, p.189-201
7. Чепурненко А.С., Андреев В.И., Языев Б.М. Построение модели
равнопрочной многопролетной балки [Электронный ресурс] // «Инженерный вестник Дона», 2013, №1. - Режим доступа:
http://ivdon.ru/magazine/archive/n1y2013/1571 (доступ свободный) - Загл. с экрана. - Яз. рус.
8. Чепурненко А.С., Языев Б.М. Оптимизация формы поперечного сечения сжатых стержней из условия устойчивости//Научное обозрение. 2012. № 6. — С. 202-204.
9. Чепурненко А.С., Андреев В.И., Языев Б.М. Построение модели равнонапряженного цилиндра на основе теории прочности Мора //Вестник МГСУ. №5 2013, с.56-61.
10. Козельская М.Ю., Чепурненко А. С., Литвинов С.В. Применение метода Галёркина при расчете на устойчивость сжатых стержней с учетом ползучести [Электронный ресурс] // «Инженерный вестник Дона», 2013, №2.
- Режим доступа: http://ivdon.ru/magazine/archive/n2y2013/1714 (доступ
свободный) - Загл. с экрана. - Яз. рус.
