∂ u 1ε ij= u i j+ u j i+ u k iu, , , k , j; x i∈ V ;∂t2( 0)σ = D ε − c ) ; σ n = p + Φ ( t); x i∈ S2; (7)2iεij= eij+ c ;ij ij,j+ fi= ρ2ijσ ; ( )ijkl(kl klijui= u i; x i∈ S1; u ( x ,0) = c ( x ,0) = 0 ,i iij iгде n – единичная нормаль к границе тела, j = 1,2,3; {D} – тензор упругихсвойств материала; u – не изменяющиеся во времени известные значенияimaxперемещений точек поверхности S 1 ; Φ () = p A sin( Ω t + )ISSN 2078-9130. Вісник НТУ «ХПІ». 2012. № 55 (961) 27jiiii∑ ∞k=1t β , р i max – амплитудысоответствующих компонент поверхностных нагружений;k2k2kk/A = a + b ; Ω = 2πkT . Рассматриваем прямоугольные циклы изменениятемпературы от T max до T min и обратно. Считаем, что все точки тела в каждыймомент нагреты до одной и той же температуры. В работе [11] показано,что влиянием переходных изменений температурного поля на осредненнуюскорость ползучести при смене уровня нагрева во многих случаях можнопренебречь.Асимптотические разложения (4), (5) и усреднения на периоде позволяютразделить систему (7) на две: первая получается из соотношений, остающихсяпосле усреднения, и описывает процессы, происходящие в медленномосновном движении, а вторая – в результате вычитания из (7) первой системы.Таким образом, для определения напряженно-деформированного состояниятел при периодическом нагружении и изменении температуры решаетсяследующая система уравнений, в которой все неизвестные изменяются вмасштабе медленного времени:σ ij , j+ f i= 0; x i∈ V ; σ ijn j = pi; x i∈ S2;1ε ij = C ijklσkl+ cij; ε ij = ( u i,j + u j,i + uk,iuk,j ); x i∈ V;(8)2ui= u i; x i∈ S 1; ui( xi,0) = cij( xi,0) = 0с уравнениями состояния в виде (6).Для конкретизации предложенных в работе уравнений состояния (6)предварительно определяются компоненты тензоров амплитудных значенийнапряжений по решению краевой задачи в масштабе быстрого времени(k = 1,2,…):akak( ) 2ak maxσ = −ρΩ u ; x V ;ij, jk ii∈ σ ij n j = piAk; x i∈ S2; (9)ak akak( u + u ) C σak 1akε = = ; x V ;iji,j j,i ijmn mn i∈ ui= 0 , x i∈ S1.2Таким образом, решением системы (9) для каждой гармоники k находятсякомпоненты тензоров амплитудных значений напряжений, по которымkkk
определяются коэффициенты g σ i( ), g σ e( )n M kr M k, конкретизирующие уравнениясостояния (6).В работе использована конечноэлементная формулировка задач ползучеститонких неосесимметрично нагруженных оболочек вращения. Оболочкавращения с произвольной образующей при неосесимметричном напряженномсостоянии разбивается на конечные элементы, которые образованыдвумя поперечными сечениями, проходящими через ось вращения коническойоболочки. Таким образом, в результате использования метода конечныхэлементов решается система алгебраических уравнений[ K ]{ Δ } = { P } { } { } { }v+ Pc+ Pp+ Pn, (10)где [K] – глобальная матрица жесткости; { Δ } – вектор узловых перемещений;{P ν vT} – обобщенный вектор внешних узловых сил, { P } [ ] { p}dS ; {P с }∑ ∫ Φ28 ISSN 2078-9130. Вісник НТУ «ХПІ». 2012. № 55 (961)=edS , {P p } – обобщенный вектор узловых усилий– обобщенный вектор узловых сил, обусловленный необратимыми деформа-cTциями, { P } = ∑ ∫ [ B] [ R]{ c}e eSpот проекции обобщенных сил на нормаль, { } ∑ ∫ [ ΦT p] { p }=eeSeSP dS ; {P n } –обобщенный вектор узловых усилий, обусловленный нелинейной состав-nTляющей упругих деформаций, { } [ B][ R]{ }P ε dS .= ∑ ∫eСистема уравнений (10) дополняется уравнениями состояния (6).Для определения амплитудных значений компонент напряженнодеформированногосостояния для каждой гармоники предварительно необходиморешить систему уравнений (9), конечноэлементная формулировкакоторой принимает следующий вид:ak akK − Ω2 M Δ = P ;([ ] [ ]){ } { }kTгде [M] – матрица масс системы, [ ] = [ Φ] [ Φ]eS∑∫eeSnM ρ dS ; Ω k – частота внешнегонагружения; Δ ak – вектор амплитудных значений при полигармоническомнагружении в пределах цикла.Пример расчета. Рассмотрим задачу о ползучести цилиндрической оболочкис шарнирно опертыми краями, нагруженной внутренним периодическиизменяющимся по прямоугольному циклу давлением.Оболочка изготовлена из сплава S-321. Физико-механические характеристикисплава определены по диаграммам деформирования, кривым ползучестии справочным данным: E = 1,67×10 5 МПа; ρ = 7,8×10 3 кг/м 3 ; ν = 0,3;b = 2,07×10 6 МПа -n /час; n = 4,054; d = 0,219 МПа -r /час; r = 7,62; m = l = 2,07;Q = Q = 4,09×10 4 К. В кинетическом уравнении для параметра повреждаемо-
- Page 1 and 2: ISSN 2078-9130ВІСНИКНАЦІ
- Page 3: Вісник Національно
- Page 6 and 7: принят кандидатом
- Page 9 and 10: А. М. Журавлевой и О.
- Page 11 and 12: ции (1976 г.), орденом
- Page 13 and 14: ук.- Х.: 1955. - 12 с. 4. Бо
- Page 15 and 16: следующие формулы:
- Page 18 and 19: m+ 1 tε ω+( ) ( ) .1 0,5∫ gmt
- Page 20 and 21: ты количества движ
- Page 22 and 23: туды импульсов U с ,
- Page 24 and 25: напряжения на конд
- Page 26 and 27: напряжений построе
- Page 30 and 31: сти (6) выбираем зна
- Page 32 and 33: Введение. Одним из
- Page 34 and 35: ∂ψ ∂ψx y ∂ψ ∂ψψ xy = +
- Page 36 and 37: На рис. 2, а показано
- Page 38 and 39: c = 10 5 Н/м (кривая 5) е
- Page 40 and 41: которой величина н
- Page 42 and 43: Рисунок 11 - Одиннад
- Page 44 and 45: симметричной конст
- Page 46 and 47: Рисунок 24 - Третий в
- Page 48 and 49: Рисунок 36 - Шестой в
- Page 50 and 51: Рисунок 48 - Девятый
- Page 52 and 53: Рисунок 60 - Двенадц
- Page 54 and 55: Рисунок 72 - Пятнадц
- Page 56 and 57: variables functions and ANSYS. By t
- Page 58 and 59: числения были повт
- Page 60 and 61: Ψ - угол между осью x
- Page 62 and 63: Вычисление микро н
- Page 64 and 65: Табл. 4 и 5 показываю
- Page 66 and 67: вен разности двух д
- Page 68 and 69: безопасную работу
- Page 70 and 71: Из предыдущего опы
- Page 72 and 73: Однако, изменения т
- Page 74 and 75: УДК 539.1С. Н. ИСАКОВ,
- Page 76 and 77: абРисунок 1 - График
- Page 78 and 79:
k1 = -0,1; k2 = 8,0167; k3 = -13,75
- Page 80 and 81:
Рисунок 5 - Распреде
- Page 82 and 83:
УДК 519:539:534С.В. КРАСН
- Page 84 and 85:
жидкости. Схема при
- Page 86 and 87:
Максимальные велич
- Page 88 and 89:
2( x,t)∂ uu( x,t)= 0, = 0 при
- Page 90 and 91:
Тогда для изображе
- Page 92 and 93:
УДК 539.3О. О. ЛАРІН, к
- Page 94 and 95:
Рисунок 1 - Схема ко
- Page 96 and 97:
Під час досліджень
- Page 98 and 99:
Окрім СКЗ вібропри
- Page 100 and 101:
Список литературы:
- Page 102 and 103:
стях деформации (пр
- Page 104 and 105:
гдеcrεij - тензор ско
- Page 106 and 107:
абРисунок 5 - Перера
- Page 108 and 109:
Выводы. Для изучени
- Page 110 and 111:
нутого образца. В с
- Page 112 and 113:
pгде εi- интенсивнос
- Page 114 and 115:
= H11εr + H12ε θ ;θ H ε 21 r+
- Page 116 and 117:
Кинетический закон
- Page 118 and 119:
личение давления а
- Page 120 and 121:
Особенности модели
- Page 122 and 123:
расчетных точек в э
- Page 124 and 125:
сеточной дискретиз
- Page 126 and 127:
Рисунок 4 - Поперечн
- Page 128 and 129:
абРисунок 8 - Амплит
- Page 130 and 131:
ции такого вида пол
- Page 132 and 133:
моделирующих навес
- Page 134 and 135:
J Ψ + c ϕazy∗2⎡ 1my + byy+ c
- Page 136 and 137:
Структурная схема
- Page 138 and 139:
ключить влияние ве
- Page 140 and 141:
Общие соотношения
- Page 142 and 143:
Колебания прямоуго
- Page 144 and 145:
s K - положительные к
- Page 146 and 147:
2 ⎧ 4G( 0,0, − sK) = f2( δ )
- Page 148 and 149:
ленном крае пласти
- Page 150 and 151:
чия жесткого тела н
- Page 152 and 153:
Чтобы определить н
- Page 154 and 155:
( n+1)1εK= arctg, (21)( n)χ ⎡*
- Page 156 and 157:
Рисунок 1 - Графики
- Page 158 and 159:
С целью дальнейшей
- Page 160 and 161:
упругой анизотропи
- Page 162 and 163:
тигло 56 %, а σ Т упал
- Page 164 and 165:
УДК 539.3С. Ю. СОТРИХИ
- Page 166 and 167:
Рисунок 2 - Блок-схе
- Page 168 and 169:
УДК 539.1А. В. СТЕПУК,
- Page 170 and 171:
сти деформирования
- Page 172 and 173:
УДК 534.1:539.3А. Н. ШУПИ
- Page 174 and 175:
Рассмотрим цилиндр
- Page 176 and 177:
Рисунок 2 - Кристалл
- Page 178 and 179:
стержни отжигались
- Page 180 and 181:
с. 12. Шупиков А. Н. Не
- Page 182 and 183:
- неравномерность п
- Page 184 and 185:
Меридиональное сеч
- Page 186 and 187:
Рисунок 4 - Расчетны
- Page 188 and 189:
ЗМІСТКедровская О.
- Page 190:
НАУКОВЕ ВИДАННЯВІС