Учебное пособие по курсу лекций (1164067), страница 43

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

Текст из файла (страница 43)

(В нашемпримере с кривой, проходящей через две точки, это означает, что все онидолжны удовлетворять граничным условиям y (a ) = A, y (b) = B ).Удовлетворять условиям полноты – ряд, построенный на этих функцияхдолжен быть сходящимся. Обычно пользуются рядами, сходимостькоторых доказана – тригонометрическим и степенным.Тогда функционал J ( y ) принимает вид:⎤⎡∞J [ y ( x )] = J ⎢∑ aiϕ i ( x )⎥ .⎦⎣ i =1Поскольку координатные функции ϕ i (x) заранее заданы, тофункционал становится функцией коэффициентов ai 63:J = J (ai ) .Условие стационарности функционалаδJ = 0в этом случае преобразуется в систему бесконечного числа уравнений:∂J=0i = 1,2…, ∞∂aiНа практике поступают следующим образом. Выстраиваютпоследовательность решений, ограничиваясь сначала одним, затем двумя и62Искомыми функциями в задачах теории пластичности является функцииизменения скоростей материальных точек в объеме деформируемого тела.63Зависит от переменных, а не от функции.282т.д.

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

Поэтому эти функции получилиназвание подходящих функций.Последовательность решения задачи определения напряженнодеформированного состояния методом Ритца.Выбор подходящих функций и количества членов ряда.Определение компонентов поля скоростей в соответствии с условиемпостоянства объема.Определение компонентов тензора деформаций.Вычисление интегралов в формуле для функционала полной мощности.Точное интегрирование во многих случаях оказывается невозможным,тогда использую приближенные, а также численные методы.Составляют и решают системы уравнений условия стационарностифункционала ∂J ∂ai = 0 для определения неизвестных варьируемыхпараметров ai .По найденным значениям коэффициентов ai определяют полескоростей vi и, затем, полную мощность.

Приравняв полную мощностьмощности внешних сил на известных скоростях, определяют удельныесилы деформирования.По найденному полю скоростей при необходимости можно определитьконечное формоизменение (изменение формы деформируемого тела) иполе напряжений.4.7.5 Определение напряженного состояния по заданному полюскоростей.После определения поля скоростей по результатам решениявариационной задачи можно определить поле напряжений. Задача283определения поля напряжений по полю скоростей имеет общее значение,поскольку поле скоростей может быть определено и другими методами64.Пусть v( x, y, z ) - поле скоростей деформируемого тела.

Тогда полескоростей деформации можно определить по известным из теориидеформаций формулам:∂v ⎞∂v11 ⎛ ∂vε ij = (vi , j + v j ,i ), например ε xx = x , ε zy = 2γ zy = ⎜⎜ z + y ⎟⎟2 ⎝ ∂y∂x∂z ⎠2Теория пластического течения постулирует следующую взаимосвязьдевиаторов напряжений и скоростей пластической деформаций:Dε = dλ ⋅ Dσ ,где переменная характеристика среды dλ пропорциональна отношениюинтенсивности скоростей деформации ε i и интенсивности напряжений σ i :3 εidλ =2σiЗначение интенсивности скоростей деформации в произвольной точкеможет быть определено по компонентам тензора скоростей деформаций:22I 2 (Dε ) =eij eij =εi =33()2(ε xx − ε yy )2 + (ε yy − ε zz )2 + (ε zz − ε xx )2 + 6 ε xy2 + ε yz2 + ε zx23Согласно критерию пластичности Мизеса:σi =σsТаким образом, характеристика среды dλ определена в каждой точкесреды.

Если пренебречь упругой деформацией, то средняя скоростьдеформации ε cp = 1 3 ⋅ (ε xx + ε yy + ε zz ) = 0 и девиатор скоростей деформаций=равен тензору: eij = ε ij − δ ij ε cp = ε ij . Тогда компоненты девиатора напряжениймогут быть определены:2σi2σi2σis xx = σ xx − σ cp =ε xx , s yy = σ yy − σ cp =ε yy , s zz = σ zz − σ cp =ε zz ,3 εi3 εi3 εi2σi2σi2σis xy = σ xy = τ xy =ε xy , s yz = σ yz = τ yz =ε yz , s zx = σ zx = τ zx =ε zx .3 εi3 εi3 εiНедиагональные компоненты тензора и девиатора напряжений равны,следовательно, касательные напряжения могут быть определенынепосредственно.

Нормальные же напряжения могут быть определенытолько с точностью до среднего напряжения (гидростатического давления).Однакоопределенностькасательныхнапряженийпозволяетпроинтегрировать уравнения равновесия. Например:64Например, известен визиопластический метод экспериментальноаналитического определения напряженного состояния по результатамэкспериментального определения поля скоростей в деформируемом теле.284известноизвестно∂σ x ∂(τ yx ) + ∂ (τ zx ) = 0+∂x∂y∂zТогда:⎛ ∂τ∂τ ⎞σ x = − ∫ ⎜⎜ yx + zx ⎟⎟dx + ϕ1 ( y, z ) .∂z ⎠⎝ ∂yЗдесь ϕ1 ( y, z ) - произвольная функция, определяемая из граничныхусловий (обычно на свободном контуре). Аналогично могут быть получены идругие компоненты тензора напряжений.⎛ ∂τ∂τ ⎞σ x = − ∫ ⎜⎜ yx + zx ⎟⎟dx + ϕ1 ( y, z )∂z ⎠⎝ ∂y⎛ ∂τ xy ∂τ zy ⎞⎟⎟dy + ϕ 2 ( x, z )+∂∂xz⎝⎠∂τ ⎞⎛ ∂τσ z = − ∫ ⎜⎜ xz + yz ⎟⎟dz + ϕ 3 ( z, x )∂y ⎠⎝ ∂zВ.Л.Колмогоров предлагает решать вариационную задачу сиспользованием принципа виртуальных скоростей и напряжений.

Инымисловами и поле скоростей и поле напряжений считаются допустимыми иодновременно варьируются. В своих работах он предложил функционал,стационарность которого соответствует одновременно истинному полюскоростей и напряжений. Функционал имеет достаточно сложный вид ивыходит за рамки программы курса.σ y = − ∫ ⎜⎜4.7.6 Решениезадачиосадкицилиндрическойвариационным методом (по И.Я.Тарновскому).заготовки4.7.6.1 Постановка задачи и расчетная схема.Определить деформированное состояние и вычислить удельную силудеформирования при осадке цилиндрической заготовки. В решении учестьискажение формы поперечного сечения.Примем следующие допущения:материал – изотропное, идеальное жесткопластическое тело,контактное трение подчиняется закону Прандтля – Зибеля,искаженная форма поперечного сечения заготовки (форма бочки) близка кпараболической кривой,инструмент движется с постоянной скоростью, инерционнымипроцессами пренебрегаем.Ранее было показано, что при пренебрежении инерционнымипроцессами и одинаковом трении на обеих контактных поверхностях процессосадки одним деформирующим инструментом со скоростью v эквивалентеносадке двумя деформирующим инструментами, движущимися другнавстречу другу с одинаковой скоростью 0.5v .

Математически это285эквивалентно использованию подвижной системы координат, связанной ссечением, проходящим через центр высоты заготовки.z∆hh−∆ha0ϕ0Rρa1ϕ1a2ϕ2(z)a2ϕ2(ρ)При этих допущениях заготовка будет находиться в осесимметричномнапряженном состоянии.Решаем задачу методом Ритца.4.7.6.2 Выбор подходящих функцийЭкспериментальный анализ деформированного состояния заготовки,выполненный методом делительных сеток показывает, что горизонтальныесечения, плоские до деформации, после деформации получают прогиб,форма которого близка к параболической. Прогиб отсутствует в среднемсечении и в сечениях, непосредственно примыкающих к контактнойповерхности.Поскольку мы имеем экспериментальные данные о поле перемещений,а не о поле скоростей, то наиболее обоснованным подходом являетсяварьирование полем перемещений. В этом случае функционал, который надоминимизировать, чтобы получить истинное поле перемещений будет286соответствовать принципу минимума полной энергии, а не принципу полноймощности, как в случае виртуальных скоростей.Для осесимметричного напряженного состояния невозможнонезависимо задавать поля перемещений, поскольку они связаны друг сдругом условием неразрывности:∂u ρ u ρ ∂u zε ρ + εθ + ε z = 0 ,=0++ρ∂z∂ρПоэтому в качестве кинематически возможного независимо можновыбрать только одно поле перемещений, например u z , поле перемещений u ρдолжно быть определено из дифференциального уравнения неразрывности.Поле перемещений u z должно удовлетворять граничным условиям:u z z =0 = 0; u z z =h = − ∆hи воспроизводить искажение поперечных сечений.Для этого аппроксимирующая функция перемещений должнанелинейно зависеть как от координаты z , так и от координаты ρ .И.Я.Тарновскийпредложилаппроксимироватьполеперемещенийследующим степенным рядом, состоящим из трех членов:u z ( z, ρ ) = a0ϕ 0 ( z ) + a1ϕ1 ( z ) + a2ϕ 2 ( z, ρ )Первый член ряда подбирают таким образом, чтобы он отражалфункцию перемещений в условиях равномерной деформации, а два другихслагаемых как бы «уточняют» распределение перемещений.В условиях равномерной деформации функция ϕ 0 ( z ) должна бытьлинейной, а коэффициент a0 определится из граничных условий.

Тогда:∆ha0ϕ 0 ( z ) = −z = −εzhДве других функции должны быть подобраны таким образом, чтобыотражать нелинейную зависимость перемещения u z от z , ρ и не нарушатьграничных условий. И.Я. Тарновский предложил следующие подходящиефункции:⎛ρ 2 ⎛ z3 ⎞z3 ⎞⎟⎜ϕ1 ( z ) = − z⎜1 − 3 ⎟;ϕ 2 ( z, ρ ) = − z 2 ⎜⎜1 − 3 ⎟⎟R ⎝ h ⎠⎝ h ⎠Окончательно:⎡⎛ρ 2 ⎛ z 3 ⎞⎤z3 ⎞⎟⎜u z = − ⎢εz + a1 z ⎜1 − 3 ⎟ + a2 z 2 ⎜⎜1 − 3 ⎟⎟⎥R ⎝ h ⎠⎦⎝ h ⎠⎣65Тогда65d ⎡ ⎛z 3 ⎞⎤ ⎛z3 ⎞ ⎛z2 ⎞ ⎛z3 ⎞⎢ z ⎜1 − ⎟⎥ = ⎜1 − ⎟ + z ⎜ − 2 3 ⎟⎟ = ⎜⎜1 − 4 3 ⎟⎟dz ⎣ ⎜⎝ h 3 ⎟⎠⎦ ⎜⎝ h 3 ⎟⎠ ⎜⎝h ⎠ ⎝h ⎠287⎡ ⎛⎡⎛ρ2 ⎛ρ 2 ⎞⎛z3 ⎞z 3 ⎞⎤z 3 ⎞⎤ε z = − ⎢ε + a1 ⎜⎜1 − 4 3 ⎟⎟ + a2 2 ⎜⎜1 − 4 3 ⎟⎟⎥ = − ⎢ε + ⎜⎜ a1 + a2 2 ⎟⎟⎜⎜1 − 4 3 ⎟⎟⎥h ⎠R ⎝h ⎠⎦R ⎠⎝h ⎠⎦⎝⎣ ⎝⎣4.7.6.3 Определение компонент поля перемещений и тензора деформацийперемещение u ρ найдем из условия непрерывности, которое можнопривести к виду:1 ∂(u ρ ) = ε zρ ∂ρ ρ⎡ ⎛z 3 ⎞⎤ρ 2 ⎞⎛u ρ ρ = ∫ ⎢ε + ⎜⎜ a1 + a2 2 ⎟⎟⎜⎜1 − 4 3 ⎟⎟⎥ρdρ + C =R ⎠⎝h ⎠⎦⎣ ⎝⎛⎛z 3 ⎞⎤ ρ 4z3 ⎞ρ2 ⎡a2 ⎜1 − 4 3 ⎟⎟ + C=⎢ε + a1 ⎜⎜1 − 4 3 ⎟⎟⎥ +2h4R 2 ⎜h⎠⎝⎠⎦⎝⎣Постоянную C найдем из граничных условиях u ρρ =0= 0 .

Характеристики

Список файлов книги