Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета, 2007, №27
Покупка
Основная коллекция
Издательство:
Кубанский государственный аграрный университет
Наименование: Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета
Год издания: 2007
Кол-во страниц: 441
Дополнительно
Вид издания:
Журнал
Артикул: 640611.0001.99
ББК:
УДК:
ГРНТИ:
Скопировать запись
Фрагмент текстового слоя документа размещен для индексирующих роботов.
Для полноценной работы с документом, пожалуйста, перейдите в
ридер.
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 1 УДК 539.3 UDC 539.3 РАСЧЕТ ТЕПЛООБРАЗОВАНИЯ В ПРЕДВАРИТЕЛЬНО ДЕФОРМИРОВАННОМ РЕЗИНОВОМ ЦИЛИНДРЕ ПРИ ГАРМОНИЧЕСКОМ НАГРУЖЕНИИ CALCULATION OF WARMING IN PRELIMINARY DEFORMED RUBBER CYLINDER AT HARMONIOUS LOADING Старостенко Игорь Николаевич старший преподаватель Starostenko Igor Nikolaevich Senior lecturer Краснодарский университет МВД России, Краснодар, Россия Krasnodar university of the Ministry of Internal Affairs of Russia ,Krasnodar, Russia В настоящей статье на основании определяющих соотношений, полученных в работе, решена задача о гармоническом нагружении предварительно сжатого резинового цилиндра In the present article the problem about harmonious loading of preliminary compressed rubber cylinder is solved on the basis of the defining parities received during work Ключевые слова: МЕХАНИКА ДЕФОРМИРУЕМОГО ТВЕРДОГО ТЕЛА, ГАРМОНИЧЕСКОЕ НАГРУЖЕНИЕ, СТАТИЧЕСКАЯ ДЕФОРМАЦИЯ, МЕТОД КОНЕЧНЫХ ЭЛЕМЕНТОВ, ФУНКЦИЯ ПРИРАЩЕНИЯ ТЕМПЕРАТУРЫ Keywords: MECHANICS OF THE DEFORMABLE FIRM BODY, HARMONIOUS LOADING, STATIC DEFORMATION, THE METHOD OF FINAL ELEMENTS, FUNCTION OF THE INCREMENT OF TEMPERATURE 1. Постановка задачи. Постановка задачи включает: - вариационные уравнения движения 0 ) ( ] ) ( [ 0 0 0 0 0 , ∫ ∫ Σ = Σ − − + σ σ δ δ ρ δ d u P u k dV u F u u t j j j j j V i j ij r && ; (1) - связь между компонентами ковариантного метрического тензора и градиентами деформации 0 , , , V X x x g m j m i ij ∈ = r ; (2) - «закон сжимаемости» − − − = − 2 0 1 0 2 3 ) ( 2 1 1 χ β χ σ β F I g , 3 2 0 / χ χ β − = ; (3) - определяющие выражения для напряжений [ ] − − + + = − ij g ij ij g ij I g T a T f T f c kT G I S δ δ 1 0 2 0 5 ,0 3 3 1 ) ( 2 ) ( + − −∫ t g ij d I g t R 0 1 )) ( 3 1 ) ( )( , ( τ τ τ τ
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 2 [ ] [ ] + − − + + ij g ij g g I T I T a T f T f c 3 1 1 1 2 ) ( 3 2 3 ) ( 3 ) ( σ δ γ , (4) энтропии [ ]× + + − + = T a T f T f c I k G T T C S S g э 0 2 1 0 0 0 ) ( 4 ) ( 5,0 ln − − − × ∫ ∫ t ij ij t ij ij d g t R g d g t R g 0 0 ) ( ) , ( ) ( ) , ( η η η τ τ τ + − − − ∫ ∫ t g g t g g d I t R I d I t R I 0 1 1 0 1 1 ) ( ) , ( ) ( ) , ( 3 1 η η η τ τ τ [ ] [ ] − − + Θ − − + + ) ( 3 2 3 ) ( 3 ) ( 2 3 2 1 1 1 2 F I I T a T f T f c g g σ β − − − − − − − − − 3 2 0 1 3 0 3 2 2 0 1 2 0 2 ) ( 2 1 1 3 1 ) ( 2 1 1 5,0 χ β χ σ β χ χ β χ σ β χ [ ]× − + + − ∫ t ij ij g d g t R g T a T f T f c I T kT G 0 0 2 1 0 ) ( ) , ( ) ( 4 ) ( 5,0 τ τ τ ∂ ∂ × − − − × ∫ ∫ t g g t ij ij d I t R I d g t R g 0 1 1 0 ) ( ) , ( 3 1 ) ( ) , ( τ τ τ η η η [ ] [ ] + Θ − − + + − × ∫ 2 1 1 1 2 0 1 1 3 2 3 ) ( 3 ) ( ) ( ) , ( β η η η g t g g I T a T f T f c d I t R I − − − − − − + 2 2 0 1 2 0 2 2 3 ) ( 2 1 1 5,0 ) ( χ β χ σ β χ σ F I g 3 2 0 1 3 0 3 ) ( 2 1 1 3 1 − − − − χ β χ σ β χ , (5) диссипативной функции [ ]+ ∂ ∂ + = ∫ t ij ij d g t R t g T a T f T f c kT G w 0 0 2 0 * ) ( ) , ( ) ( 8 ) ( τ τ τ
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 3 ∫∫ ∫ − ∂ ∂ − + t t ij ij ij t d d g g t R t R t d g t R 0 0 0 ) ( ) ( ) , ( ) , ( ) ( ) , ( η τ η τ η τ η η η + + ∂ ∂ − ∫ ∫ t g t g g d I t R d I t R t I 0 1 0 1 1 ) ( ) , ( ) ( ) , ( 3 1 η η η τ τ τ ∂ ∂ + ∫∫ t t g g d d I I t R t R t 0 0 1 1 ) ( ) ( ) , ( ) , ( 3 1 η τ η τ η τ ; (6) - вариационное уравнение притока тепла [ ] ∫ + Θ Λ + Θ − 0 0 , , 3 *) ( V m n mn g dV T g I w S T δ δ & 0 ) ( ) ( 0 0 ∫ Σ = Σ Θ − + ф Ф ф T d T T k u k δ r ; (7) - условия закрепления u X X u X u 0 0 ), ( ) ( Σ ∈ = r r r r r ; (8) - начальные условия по перемещениям ) ( ) 0, ( 1 X v X u r r r r = ; (9) скоростям ) ( ) 0, ( 2 X v X u r r r &r = ; (10) температуре ) ( ) 0, ( 0 X T X T r r = ; (11) - условие теплообмена на внешней поверхности ф 0 Σ ) ( ф T T T k q − = ⋅νr r . (12) В формулах (1)-(12): 0 V - область, занимаемая телом в недеформированном состоянии; u 0 0 0 Σ ∪ Σ = Σ σ - ее поверхность, состоящая из двух кусков, на одном из которых заданы распределенные нагрузки Pr , на другом - перемещения 0 ur ; Xr - координаты Лагранжа; u X x r r r + = -
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 4 координаты Эйлера; us - вектор перемещения; 0 ρ - плотность материала в исходном состоянии; Fr - вектор объемной силы; ) (u k r - кратность изменения элементарной площадки; δ - символ вариации функции; ijt - компоненты тензора условных напряжений; ij g - компоненты метрического тензора с инвариантами gi I ; ij δ - символ Кронекера; 0 G - общее число активных молекулярных цепей в единице объема; k - постоянная Больцмана; 2 c , 0 a , 1a , 0 β , 2 χ - структурно-механические параметры резины; ] 1 [β , э С - теплофизические параметры резины; ) (T F , ) (T f , ) ( ] 1 [ T γ - известные функции температуры; σ - функция типа гидростатического давления; ) , ( τ t R - функция релаксации; 0 T T − = Θ - приращение температуры; Λ - коэффициент теплопроводности; Tk - коэффициент теплоотдачи; ) ( 1 X v r r , ) ( 2 X v r r , ) ( 0 X T r - известные функции координат; ф T - температура внешнего тела. Система (1)-(12) замкнута [2]. Задача заключается в нахождении вектор-функции ) , ( t X u r r и скалярных функций ) , ( t Xr σ , ) , ( t X T r , удовлетворяющих системе (1)-(12). 2. Уточнение внешних воздействий. Граничные условия имеют вид ≥ − + + ≤ = ,0 , sin ) , ( cos ) , ( ) , ( ) , ( ; ), , ( ) , ( 1 1 1 0 1 1 0 1 0 * ) 0 ( ) 0 ( 0 t t t X u t t X u t X u t X u t t t X u t X u s c N ω ω r r r r r r r r r r r r (13) ≥ − + + ≤ = ,0 , sin ) , ( cos ) , ( ) , ( ) , ( ; ), , ( ) , ( 1 1 1 0 1 1 0 1 0 * ) 0 ( ) 0 ( t t t X P t t X P t X P t X P t t t X P t X P s c N ω ω r r r r r r r r r r r r (14) где ) 0 ( ) 0 ( , P u r r - медленные функции времени, задающие процесс предварительного нагружения; 0 *ur , 0 cur , 0 sur , 0 *Pr , 0 cPr , 0 sPr - амплитуды
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 5 догружающих воздействий, Nt t t − = 1 , Nt - длительность предварительного нагружения. Условия (13), (14) позволяют линеаризовать задачу в окрестности промежуточного состояния и рассматривать малые моногармонические колебания предварительно деформированного вязкоупругого тела. 3. Построение приближенного решения. Приближенное решение задачи (1)-(12) ищем в виде ≥ − + + ≤ = ,0 , sin ) , ( cos ) , ( ) , ( ) , ( ; ), , ( ) , ( 1 1 1 1 1 1 * ) 0 ( ) 0 ( t t t X u t t X u t X u t X u t t t X u t X u s c N ω ω r r r r r r r r r r r r (15) ≥ − + + ≤ = ,0 , sin ) , ( cos ) , ( ) , ( ) , ( ; ), , ( ) , ( 1 1 1 1 1 1 * ) 0 ( ) 0 ( t t t X t t X t X t X t t t X t X s c N ω σ ω σ σ σ σ σ r r r r r r (16) ) , ( t Xr Θ = Θ . (17) где ) , ( ), , ( ) 0 ( ) 0 ( t X t X u r r r σ - части решения, отвечающие предварительному нагружению; ) , ( 1 * t X u r r , ) , ( 1t X uc r r , ) , ( 1t X us r r , ) , ( t Xr Θ - медленные функции времени. Для исключения из системы «быстрых» функций используем метод усреднения [3]. Линеаризованную задачу разделим на две подзадачи, первую из которых задает система квазилинейных дифференциальных уравнений движения относительно амплитудных переменных ), , ( t X uc r r ), , ( t X us r r ), , ( t X c r σ ) , ( t X s r σ с соответствующими начальными условиями, а также краевую задача статики относительно функций ) , ( * t X u r r , ) , ( * t Xr σ , вторую - определяет уравнение энергии (7). Особенностью изучаемой задачи в целом является присутствие в ней термомеханического
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 6 сопряжения, а также зависимости от параметров предварительного нагружения. В конечноэлементной формулировке задача имеет вид = − = + ,0 , * ) 0 ( * * * * σ α σ N i i N M j j M N i ij MN u D P D u A (18) = − = + + + ,0 , ) 0 ( ) ( ) ( ) ( c N ci i N M cj N si c ij MN N si c ij MN c j M N ci c ij MN u D P u C u B D u A σ α σ & (19) = − = + + + ,0 , ) 0 ( ) ( ) ( ) ( s N si i N M sj N ci s ij MN N ci s ij MN s j M N si s ij MN u D P u C u B D u A σ α σ & (20) ) ( ) ( ) ( ) ( ) ( e e e e e A B E D dt d C r r r r r − = Θ + Θ + Θ , (21) где C , D , E - соответственно матрицы теплоемкости, теплопроводности и теплоотдачи, ) (e Br - вектор обобщенной диссипации, ) (e Ar - вектор, обусловленные различием начальных температур тела и окружающей среды. Коэффициенты матричных уравнений (18)-(21) зависят от уровня предварительной деформации и текущего значения температуры. Дискретизация по времени выполним в виде простейшей шаговой процедуры. Применительно к шагу с номером )1 ( + r алгоритм имеет вид: 1) По известному распределению температуры ) , ( rt X T r в момент времени rt t = решаем задачу механики, определяем поля амплитуд перемещений ) , ( * rt X u r r , ) , ( r c t X u r r , ) , ( r s t X u r r и амплитудные значения функции типа гидростатического давления ) , ( * rt Xr σ , ) , ( r c t Xr σ , ) , ( r s t Xr σ ;
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 7 2) Формируем и решаем систему уравнений притока тепла, в результате находим распределение температуры ) , ( 1 + rt X T r на )1 ( + r - ом временном слое; 3) С учетом модифицированного температурного поля вновь переходим к пункту 1) алгоритма и процесс замыкается. 5. Саморазогрев предварительно сжатого резинового цилиндра. Решим задачу для кубообразного цилиндра = = h d 0,06 м с неоднородными кинематическими граничными условиями. На торцах, контактирующих с металлической арматурой, заданы встречные перемещения с одинаковыми амплитудами, а свободная от нагрузок боковая поверхность цилиндра находится в условиях конвективного теплообмена с воздухом. Систему координат связываем с главными центральными осями цилиндра: для определенности координатную плоскость 2 1Ox x считаем горизонтальной, ось 3x направляем вертикально вверх и совмещаем с осью цилиндра. Амплитуды догружающих вертикальных перемещений принимаем равными 0,003м (10% от начальной полувысоты). Интегрирование задачи выполняем при однородных механических начальных условиях, кроме того, полагаем, что в начальный момент времени температура во всех точках тела одинакова и равна температуре окружающей среды. Использованы призматические конечные элементы с билинейной аппроксимацией перемещений, температуры и постоянной функцией типа гидростатического давления. Расчеты выполнены при следующих исходных данных: Т T X T э ф 295 ) ( 0 = = = r ; Гц 40 = ω ; 1 4 10 5,1 − − ⋅ = К э α ; [ ] 1 4 1 10 0,1 − − ⋅ = К β ; ) /( 22 ,0 К м Вт ⋅ ⋅ = Λ ; ) /( 1153 2 ) ( К м Вт k м T = ;
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 8 ) /( 8,8 2 ) ( К м Вт k в T = ; /( 2000 3К м КДж Cэ = ; 3 0 / 1200 м кг = ρ , где ) ( м Tk - коэффициент теплоотдачи на границе контакта резины и металла, ) (в Tk - соответственно резины и воздуха [4,5]. Фактически решены следующие три задачи: 1) Предварительной деформации нет. Контакт между резиной и металлом идеальный (без трения); 2) Предварительной деформации нет. На границе контакта резина сцеплена с металлом; 3) Цилиндр поджимаем на 20% от первоначальной высоты, контакт идеальный. В практических расчетах было реализовано 24 внешних цикла по механической задаче, внутри каждого из которых за 10 шагов осуществлялось интегрирование уравнения энергии. Решение построено на временном интервале ] 12000 ,0 [ 1 ∈ t с. На рис. 1-4 даны функции избыточной температуры по двум горизонтальным сечениям цилиндра плоскостями 0 3 = X , h X 4375 ,0 3 = в момент времени 12000 1 = t с. Видно, что во всех случаях максимум температуры реализуется в центре области. Изображенные фигуры являются поверхностями вращения. На рис. 5 показаны графики роста температуры в центре цилиндра для каждой задачи. Графики носят «мягкий» характер, по истечению определенного времени температура стабилизируется. Влияние предварительного поджатия неоднозначно. С одной стороны поджатие способствует росту некомпенсированной теплоты, с другой - приводит к количественной перестройке матриц
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 9 теплопроводности и теплоотдачи и, в конечном итоге, к снижению избыточной температуры в рассматриваемой точке. В итоге температура в предварительно деформированном цилиндре может оказаться несколько ниже, чем в изначально недеформированном теле. Наибольшая стартовая скорость разогрева зарегистрирована в задаче 2 и обусловлена наличием неоднородного напряженно-деформированного состояния. По мере приближения к торцам цилиндра, где происходит интенсивный теплообмен с металлической арматурой, температура падает. Для исследования влияния числа шагов по параметру изменения внешних воздействий на интервале [ ] с t 500 ,0 1 ∈ первая задача была решена при t t ∆ = ∆ 1,0 1 , t t ∆ = ∆ 01 ,0 1 . Десятикратное увеличение числа шагов вносит поправку ∼2%. В критериях разрушения температура играет определяющую роль, это доказывает актуальность полученных решений. Рис. 1. Функция приращения температуры по сечению 0 3 = X (задача 1)
Научный журнал КубГАУ, №27(3), март 2007 года http://ej.kubagro.ru/2007/03/pdf/23.pdf 10 Рис. 2. Функция приращения температуры по сечению H X 4375 ,0 3 = (задача 1) Рис. 3. Функция приращения температуры по сечению 0 3 = X (задача 2)