Метод Ньютона-Канторовича для расчета параметров равновесия при растяжении с кручением специального образца

Размещено на http://www.allbest.ru/

Метод Ньютона-Канторовича для расчета параметров равновесия при растяжении с кручением специального образца

Привалова В.В.

Екатеринбург, Россия

Рассмотрим стержневую систему, осуществляющую растяжение с кручением специального образца (рис. 1). Стержни 1 и 2 выполнены из линейно упругого материала. Стержень 1 передаёт на образец 3 растягивающую нагрузку. Жёсткость стержня 1 при растяжении равна . Стержень 2 передаёт на образец закручивающий момент. Жёсткость стержня 2 при кручении ? . Полый образец 3 имеет такую геометрию, что удлинение по величине равно деформации растяжения , угол закручивания ? деформации сдвига , растягивающая сила ? растягивающему напряжению , а крутящий момент ? касательному напряжению . Система нагружается либо заданием перемещения u сечению AA стержня 1 и угла закручивания в сечении DD стержня 2 (жёсткое нагружение), либо заданием соответственно растягивающей силы и крутящего момента (мягкое нагружение). Возможны также комбинации кинематических и силовых нагружений (смешанное нагружение).

Размещено на http://www.allbest.ru/

Рис. 1. Механическая система

Свойства материала образца при растяжении с кручением заданы в общем случае невыпуклым потенциалом [2]:

где E=2*104 кГ/ммІ, G=7,7*10і кГ/ммІ, , a=100р/E. Этот потенциал описывает состояния материала образца на всех стадиях его деформирования (устойчивых и неустойчивых). Таким образом, учитывается как состояние упрочнения материала, так и состояние разупрочнения. Требуется найти предельные значения параметров нагружения и параметры равновесия данной механической системы.

При подборе определённых геометрических параметров системы [ссылка на предыдущую статью из Вестника, где всё было подробно описано] потенциальную функцию запишем в виде

Здесь первые два слагаемых ? потенциальные энергии упругих деформаций стержней 1 и 2, , ? параметры управления, , ? параметры состояния системы. Величины и фиксируем (=500 кГ/мм, =500 кГ*мм). Критические точки функции есть решения системы уравнений равновесия [3]

Уравнения равновесия являются неоднозначным отображением пространства параметров управления на пространство параметров состояния. Для определения параметров состояния по заданным параметрам управления выделим области однозначности следующим образом.

Из уравнений равновесия выразим параметры управления

(1)

Вырожденные критические точки, образующие в пространстве управлений сепаратрису, обращают в нуль детерминант матрицы Гессе функции (матрицы вторых производных). Приближённые значения координат вырожденных критических точек определяем следующим образом. Возьмём на множестве прямоугольник , и построим в нём сетку узлов с шагом 2*10-4 (по ) и 3,2*10-4 (по ). В каждом узле вычисляем значения детерминанта матрицы Гессе и выделяем те узлы, в которых детерминант близок к нулю с достаточной степенью точности (рис. 2). Координаты этих узлов подставляем в уравнения (1) и находим значения . В пространстве управления точки с координатами расположены на сепаратрисе функции (или в достаточной близости от неё).

На рис. 3 изображена полученная сепаратриса функции . Потеря устойчивости процесса деформирования происходит тогда, когда путь нагружения в пространстве управлений выходит из области, ограниченной кривыми сепаратрисы, т.е. пересекает кривую A'B' (рис. 3). Координаты точек этой кривой и определяют предельные значения параметров управления.

Размещено на http://www.allbest.ru/

Рис. 2. Нули детерминанта матрицы Гессе функции

Размещено на http://www.allbest.ru/

Рис. 3. Сепаратриса при жёстком нагружении

При отображении, описывающемся уравнениями равновесия пространства параметров управления на пространство параметров состояния получаем наложение значений, что даёт неоднозначность отображения. Таким образом получаем:

1) каждой паре значений параметров управления из области О'D'С' (т.е. области ограниченной координатными осями и кривой С'D') соответствует только одна пара значений параметров состояния из области OAB;

2) каждой паре значений параметров управления из области С'D'B'A' (т.е. области ограниченной координатными осями, кривой С'D' и кривой A'B') соответствует три пары значений параметров состояния: одно из области OAB, второе из области ABCD и третье из области, ограниченной координатными осями и лежащей за пределами кривой СD;

3) каждой паре значений параметров управления из области ограниченной координатными осями и лежащей за пределами кривой С'D' соответствует одна пары значений параметров состояния: одно из области, ограниченной координатными осями и лежащей за пределами кривой СD.

Таким образом получаем 3 области однозначности решений.

Определение параметров состояния по заданным параметрам управления проведём в следующих трёх этапах:

? Разбиваем сеткой плоскость параметров состояния разделяя точки сетки на 3 множества: точки, принадлежащие области OAB (первое множество), точки, принадлежащие области ABCD (второе множество) и точки, лежащие в первом квадранте плоскости за кривой CD (третье множество). Отображаем точки всех трёх множеств в пространство управлений и получаем, соответственно, 3 множества точек на плоскости . В результате отображения область OAB отобразилась на область O'A'B', а область ABCD - на область A'B'C'D'. Область, содержащая третье множество точек, лежит в первом квадранте плоскости за кривой C'D'. Получаем, что точкам-образам, лежащим в области A'B'C'D', отвечают 3 прообраза (из всех трёх множеств точек пространства ).

? Для заданных значений параметров управления и необходимо знать, какой области принадлежит точка с такими координатами в пространстве . Если точка принадлежит первой области, то ей будет соответствовать только один прообраз в пространстве состояний. Если точка принадлежит второй области, то ей будет соответствовать три прообраза в пространстве состояний (во всех трёх областях пространства состояний). Если же точка принадлежит третьей области, то ей будет соответствовать два прообраза в пространстве состояний (во второй и третьей областях). Итак, для заданных значений параметров и определяем, в какой из областей и сколько решений существует.

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

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

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

Опишем используемый модифицированный метод Ньютона-Канторовича для данной задачи более подробно. Модифицированный метод Ньютона-Канторовича в общем случае записывается в виде [4]

Введём обозначения:

где

Тогда модифицированный метод Ньютона-Канторовича записываем в матричном виде:

Пусть ? начальные приближения параметров состояния полученные из второго этапа решения задачи. Тогда метод можно записать в виде

Обозначим В результате получаем выражения для вычисления приближений

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

При мягком нагружении состояние механической системы описывает потенциальная функция

Здесь , ? параметры управления, , ? параметры состояния системы. Критические точки функции есть решения системы уравнений равновесия

Отсюда

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

На рис. 4 изображена полученная сепаратриса функции . Потеря устойчивости происходит после пересечения путём нагружения кривой A'B' (рис. 4). Координаты точек этой кривой и определяют предельные значения параметров управления.

Аналогично схеме проведённой для жёсткого нагружения, проводим определение параметров состояния по заданным параметрам управления проведём в трёх этапах:

? Разбиваем сеткой плоскость параметров состояния разделяя точки сетки на 3 множества. Отображаем точки всех трёх множеств в пространство управлений и получаем, соответственно, 3 множества точек на плоскости . В результате отображения область OAB отобразилась на область O'A'B', область ABCD - также на область O'A'B' со сгущением значений в точке O'. Область, содержащая третье множество точек, лежит в первом квадранте частично в области O'A'B', а также за её пределами. Таким образом, в случае мягкого нагружения, все 3 множества параметров состояния при отображении на пространство параметров управления будут пересекаться.

? Для заданных значений параметров управления и известно, что если точка принадлежит области O'A'B', то ей будет соответствовать три прообраза в пространстве состояний.

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

Размещено на http://www.allbest.ru/

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

Для случая мягкого нагружения функция определения последовательных приближений имеет вид

где

В остальном метод повторяет приведённый для случая жёсткого нагружения.

При смешанном нагружении состояние механической системы описывает потенциальная функция

Здесь , ? параметры управления, , ? параметры состояния системы. Критические точки функции есть решения системы уравнений равновесия

Отсюда

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

Размещено на http://www.allbest.ru/

На рис. 5 изображена полученная сепаратриса функции . Потеря устойчивости происходит после пересечения путём нагружения кривой A'B' (рис. 5). Координаты точек этой кривой и определяют предельные значения параметров управления.

Аналогично схеме проведённой для мягкого нагружения, проводим определение параметров состояния по заданным параметрам управления проведём в трёх этапах:

? Разбиваем сеткой плоскость параметров состояния разделяя точки сетки на 3 множества. Отображаем точки всех трёх множеств в пространство управлений и получаем, соответственно, 3 множества точек на плоскости . В результате отображения область OAB отобразилась на область O'A'B', область ABCD - также на область O'A'B' со сгущением значений на отрезке O'D'. Область, содержащая третье множество точек, лежит в первом квадранте частично в области O'A'B', а также за её пределами. Таким образом, в случае смешанного нагружения, все 3 множества параметров состояния при отображении на пространство параметров управления будут пересекаться.

? Для заданных значений параметров управления и известно, что если точка принадлежит области O'A'B', то ей будет соответствовать три прообраза в пространстве состояний.

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

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

Для случая мягкого нагружения функция определения последовательных приближений имеет вид

где

ньютон канторович нагружение растяжение

В остальном метод повторяет приведённый для случая жёсткого нагружения.

Работа выполнена при финансовой поддержке РФФИ (проект № 13-08-00186).

Литература

1. Седов Л. И. Механика сплошной среды. М. : Наука, 1970. Т. 1. 492 с.

2. Стружанов В. В., Просвиряков Е. Ю., Бурмашева Н. В. Об одном методе построения единого потенциала // Вычислительная механика сплошных сред. Т. 2. № 2. 2009. С. 96 - 107.

3. Гилмор Р. Прикладная теория катастроф: в 2-х книгах. М. : Мир, 1984. Кн. 1. 350 с.

4. Канторович Л.В., Акилов Г.П. Функциональный анализ. М. : Наука, 1977. 742 с.

Размещено на Allbest.ru