40573

ВАРИАНТ ОПРЕДЕЛЕНИЯ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ УПРУГОГО ТЕЛА КОНЕЧНЫХ РАЗМЕРОВ С ТРЕЩИНОЙ

Курсовая

Производство и промышленные технологии

Рассмотрена задача о нахождении напряженно-деформированного состояния (НДС) в поврежденном трещиной теле конечных размеров. Трещина моделируется физическим разрезом с характерной толщиной и материальным слоем на его продолжении. Напряженное состояние слоя описывается средними по толщине и граничными напряжениями, связанными условиями равновесия

Русский

2014-12-28

1.75 MB

43 чел.

ВАРИАНТ ОПРЕДЕЛЕНИЯ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ УПРУГОГО ТЕЛА КОНЕЧНЫХ РАЗМЕРОВ С ТРЕЩИНОЙ1

Аннотация

Рассмотрена задача о нахождении напряженно-деформированного состояния (НДС) в поврежденном трещиной теле конечных размеров. Трещина моделируется физическим разрезом с характерной толщиной  и материальным слоем на его продолжении. Напряженное состояние слоя описывается средними по толщине и граничными напряжениями, связанными условиями равновесия. Для смежных со слоем тел получены вариационные уравнения равновесия, решение которых определяет НДС в рассматриваемом слое. В случае линейно упругого поведения материала рассмотрены случаи деформирования тела нормального отрыва и поперечного сдвига. Из конечно-элементного решения задачи выявлены тенденции к началу пластического деформирования.

Ключевые слова: трещина, характерный размер, нормальный отрыв, чистый сдвиг, пластическое состояние.

Введение.

Исследование проблем прочности и разрушения твердых тел в настоящее время является актуальной задачей, как в теоретическом, так и в прикладном плане. Под разрушением понимается макроскопическое нарушение сплошности тела в результате воздействия на него внешнего воздействия. В настоящее время разрушение рассматривают на разных масштабных уровнях [18,22,29,31,34-36,47,51]. И в этом случае модель трещины определяет соответствующий математический аппарат для её исследования.

Существуют два основных подхода  для описания модели трещины. Первый,  заключается в  представлении трещины в виде математического разреза. Недостатком данного подхода является сингулярность поля напряжения в вершине трещины.  Введение в модель сил сцепления позволяет решить эту проблему, однако эти внешние нагрузки устанавливаются a priori, без решения соответствующей граничной задачи. Здесь стоит отметить работы Баренблатта [3,50], Гольдштейна [14-20,52], Лавита [26,27], который обобщил соответствующий подход на случай упругопластического деформирования.

Основы описания трещины как математического разреза были заложены английским ученым А. А. Гриффитсом.[54-55]  Он предложил энергетический подход для описания разрушения. Суть данного подхода состоит в том, что для  роста трещины необходимо, чтобы количество высвобождающейся потенциальной энергии должно превышать поверхностную энергию необходимую для преодоления сил взаимодействия слоев атомов. Этот подход называют энергетическим критерием разрушения.

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

Орован обобщил подход Гриффитса для материалов, при разрушении которых в кончике трещины развиваются необратимые пластические деформации. Было установлено, что пластические деформации сосредоточены в малой зоне вблизи кончика трещины. Из этого было сделано предположение, что затраты энергии в процессе создания новых поверхностей при развитии трещины главным образом связаны с работой пластической деформации объемов материала, расположенных перед фронтом трещины. Ирвин [31,56-58] установил, что процесс разрушения материала при распространении трещины обуславливается напряженно-деформированным состоянием в окрестности вершины трещины, которое в свою очередь, в линейно упругом теле определяется коэффициентом интенсивности напряжений. Поэтому предполагается, что трещина распространяется при достижении коэффициентом интенсивности напряжений некоторого критического значения. Критические значения коэффициентов интенсивности [38-44] напряжений являются постоянными материала, характеризующими его трещиностойкость при заданной температуре, внешней среде и т.п. Этот критерий разрушения получил название силового критерия разрушения. Ирвином была показана эквивалентность силового критерия разрушения энергетическому критерию Гриффитса.

В настоящее время интенсивно разрабатываются так называемые нелокальные критерии прочности [23-25], в частности интегральный критерий, или критерий средних напряжений. Этот критерий обычно связывают с именами Г. Нейбера [41] и В.В. Новожилова [42-43].

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

В настоящее время активно развивается моментная теория упругости и ее применение в теории трещин [1,2,4,38]. Однако, отсутствие оценок введенного линейного параметра для  конкретных материалов не дает право применению соответствующей теории в практике расчетов поврежденных конструкционных материалов.

Второй подход представляет трещину как физический разрез с характерной толщиной. Здесь отметим работы Прандтля, Гудьера, Канинена, где на продолжении физического разреза вводится слой со связями типа пружин, работающих на растяжение и сжатие. В работе Прандтля [60] рассматриваются два упругих тела (балки), скрепленные по всей длине за исключением конечного отрезка (трещины) поперечными упругими связями, хрупко разрушающимися по достижении некоторого удлинения, и принимаются условия предельного равновесия такой системы под действием равномерно распределенных отрывающих усилий. В модели отсутствует сингулярность, но ее применимость ограничена нормальным отрывом.

Отметим, что для трещин в виде математического и физического разрезов может быть применен  критерий J-интеграла, предлагаемый Черепановым [48] и Райсом [61-63] для хрупких материалов и обобщенный в работах [5,9,10] в случае упругопластического деформирования.

Возможность описания трещины в виде физического разреза дается в работе Ф. Маклинтока[30,59], однако, выбор величины характерного размера физического разреза, а также постановок соответствующих задач механики сплошной среды  в данной работе и последующих работах автора не приведено.

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

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

В работах В.В. Глаголева и А.А. Маркина [7,8,11-13,53] была предложена модель трещины в виде физического разреза и материального слоя на его продолжении. Напряженное состояние слоя описывается средними и граничными напряжениями, связанными условиями равновесия. Использование средних по толщине слоя напряжений позволяет отказаться от рассмотрения формы окончания физического разреза. В статье [8] на основе предложенной модели была предложена постановка и решена задача о развитии тонкой пластической области в окрестности трещины при нагружении нормальным отрывом.

Таким образом, является актуальным является развитие моделей трещины в которых отсутствует сингулярность напряжений [6,13],  что позволяет в рамках естественных критериев механики сплошной среды определить момент перехода в состояние пластичности. В данной работе на основе модели [13] рассмотрена вариационная постановка задачи об определении напряженно-деформированного состояния тела конечных размеров с трещиной. Методом конечного элемента представлены решения для частных случаев нагружения. Найдены области, согласно критерию Треска - Сен-Венана, переходящие в состояние пластичности.

1.Постановка задачи.

На Рис. 1  представлено тело с трещиной в виде физического разреза с толщиной  и неопределенной границей его окончания. Процесс нагружения предполагаем квазистатическим и изотермическим. В этом случае условие равновесия запишем в вариационной форме:

,                                                  (1.1)

где под номером 3 определена материальная область, лежащая на продолжении физического разреза; 1,2 - смежные с 3 области;  - внешняя нагрузка на контуре ;  - тензор напряжений;  - тензор напряжений;  - поле перемещений.

Воспользуемся следующими обозначениями для напряжений на границах слоя: , , , .

Рис.1.1. Тело с трещиной.

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

В слое 3 средние напряжения, деформации и перемещения определяем через их граничные значения следующим образом [13]:

,                     (1.2)    ,                     (1.3)

,                                           (1.4)

,   ,            (1.5)

  ,            (1.6)

 ,                 (1.7)

где , - вектора перемещения верхней и нижней границы области 3.

Из выражений (1.6) приходим к представлению средней сдвиговой деформации вдоль слоя:

 (1.8)

Напряжения по границе области 3: ,,, связаны со средними напряжениями условиями равновесия:

,                                                (1.9)

                                               (1.10)

при условии симметрии касательных напряжений: ,.

К внешней нагрузке, наряду с силой  относим и нагрузку со стороны слоя 3 на тело 1:

,               (1.11)

и тело 2:

.                (1.12)

Выразим граничные со слоем напряжения для тела 1:

,                              (1.13)

,                             (1.14)

и для тела 2:

,                              (1.15)

.                            (1.16)

Запишем выражения (1.11) и (1.12) с учетом (1.13)-(1.16):

           (1.17)

           (1.18)

Преобразуем слагаемые, содержащие производную от среднего напряжения.

                           (1.19)

                           (1.20)

                           (1.21)

                           (1.22)

При отсутствии нагрузок по торцам слоя получаем:

              (1.23)

              (1.24)

Подставим (1.23),(1.24) в (1.1), тогда решение задачи о равновесии тела с трещиной сводится к совместному решению двух вариационных уравнений [53] для тела 1:

                        (1.25)

и тела 2:

                      (1.26)

где ,- верхняя и нижняя границы тела 3; , - контуры приложения внешней нагрузки для тела 1 и 2.

Для решения уравнений (1.25) и (1.26) необходимо замкнуть модель конкретными определяющими соотношениями, связывающими напряжения с деформациями. В рамках данной работы связь между напряжениями и деформациями представим в виде закона Гука:

,                                                 (1.27)

.                                                  (1.28)

Здесь - модуль сдвига;  -гидростатическое напряжение;  - относительное изменение объема;  коэффициент объемного расширения,  - компоненты девиаторов напряжения и деформации.

В результате решения (1.25)-(1.28) определяется поле перемещений в телах 1 и 2, в том числе и по границе со слоем взаимодействия, что позволяет с учетом (1.5), (1.8) найти средние деформации слоя.

Для решения вариационных уравнений (1.25), (1.26) совместно с определяющими соотношениями  (1.27), (1.28) предлагается использовать метод конечного элемента с линейной аппроксимацией поля перемещений. Предполагается, что размер грани конечного элемента соизмеримы с введенным параметром .

В предлагаемой модели трещины основным параметром для конкретного материала является масштабный параметр . Решение задачи о нахождении НДС в теле с трещиной может быть получено при произвольной внешней нагрузке.

В работе [53] показано, что введенный линейный размер может быть получен из эксперимента по податливости образца с трещиной в режиме упругого деформирования. В статьях [6],[53] даны оценки толщины слоя через известные механические характеристики материала.

Рассмотрим образец с трещиной, показанный на Рис.1.2, со следующими геометрическими характеристиками: м, м, м,  м.

Рис. 1.2.  Конфигурация образца

На примере образца Рис.2 смоделируем известные в механике виды деформирования типа разрушения по моде I и моде II для материала со следующими механическими характеристиками: Па, Па. После нахождения НДС в элементах слоя и сопряженных с ним тел определяются элементы с максимальным касательным напряжением [21]:

,

где - главные напряжения тензора напряжений.

Максимальное касательное напряжение согласно критерию Треска-Сен-Венана будет определять тенденцию к началу развития пластической зоны. Для решения данной задачи используем метод конечных элементов.

2.Решение задачи.

  1.  . История развития метода конечных элементов.

Возникновение метода конечных элементов связано с решением задач космических исследований в 1950-х годах (идея МКЭ была разработана советскими учёными ещё в 1936 году, но из-за неразвитости вычислительной техники метод не получил развития). Этот метод возник из строительной механики и теории упругости, а уже затем было получено его математическое обоснование. Существенный толчок в своём развитии МКЭ получил в 1963 году после того, как было доказано то, что его можно рассматривать как один из вариантов распространённого в строительной механике метода Рэлея — Ритца, который путём минимизации потенциальной энергии сводит задачу к системе линейных уравнений равновесия. После того, как была установлена связь МКЭ с процедурой минимизации, он стал применяться к задачам, описываемым уравнениями Лапласа или Пуассона. Область применения МКЭ значительно расширилась, когда было установлено (в 1968 году), что уравнения, определяющие элементы в задачах, могут быть легко получены с помощью вариантов метода взвешенных невязок, таких как метод Галёркина или метод наименьших квадратов. Это сыграло важную роль в теоретическом обосновании МКЭ, так как позволило применять его при решении многих типов дифференциальных уравнений. Таким образом, метод конечных элементов превратился в общий метод численного решения дифференциальных уравнений или систем дифференциальных уравнений.

К семидесятым годам относится появление математической теории конечных элементов. Здесь можно выделить труды И. Бабушки, Р. Галлагера, Ж. Дек-лу, Дж. Одена, Г. Стренга, Дж. Фикса. Значительный вклад в разработку теоретических основ МКЭ внесли и российские ученые. В. Г. Корнеев указал на совпадение математической сущности МКЭ и ВРМ. Сопоставление МКЭ с рядом вариационных методов приведено в трудах Л. А. Розина. Под руководством А. С. Сахарова разработана моментная схема конечных элементов.

С развитием вычислительных средств возможности метода постоянно расширяются, также расширяется и класс решаемых задач. Практически все современные расчёты на прочность проводят, используя метод конечных элементов.

2.2. Метод конечных элементов. Введение локальной и глобальной нумерации узлов.

Метод конечных элементов есть способ аппроксимации непрерывной функции дискретной моделью, представляющей собой множество значений заданной функции в конечном числе точек области ее определения в совокупности с кусочными аппроксимациями этой функции этой функции на конечном числе подобластей.  Данные области называются конечными элементами. Локальная аппроксимация на каждом конечном элементе единственным образом определяется значением этой функции в конечном числе предварительно выбранных точек области ее определения. Таким образом, при построении конечноэлементной модели заданной функции поступают следующим образом:

  1.  В области определения функции фиксируют конечное число точек и определяют значения функции в этих точках. Выбранные точки носят название узловых точек или узлов.
  2.  Область определения функции приближенно представляют в виде совокупности конечного числа непересекающихся подобластей, называемых конечными элементами. Таким образом, областью определения функции является совокупность конечных элементов, связанных между собой в узлах.
  3.  Заданную функцию локально аппроксимируют на каждом конечном элементе непрерывными функциями, однозначно определяемыми значениями функциями в узловых точках, принадлежащих этим элементам.

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

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

Рис. 2.1. Расположение области на плоскости.

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

Рассмотрим область АВСDF в виде совокупности отдельных элементов. Соответствующее представление показано на Рис. 2.2. В этом случае для каждого отдельного элемента можно ввести свою локальную нумерацию. Полагаем, что направление обхода при локальной нумерации одинаково для всех элементов модели и при дальнейшем изложении будем вести ее против часовой стрелки.

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

Рис. 2.2. Элементы области.

Таким образом, для однозначного определения положения узла необходимо указать либо его глобальный номер, либо локальный номер в паре с номером рассматриваемого элемента. Так узловая точка F может быть идентифицирована как глобальный узел 2 или локальный узел 3 третьего конечного элемента. В дальнейшем при определении координат узлов мы будем использовать как глобальную, так и локальную нумерации. В этом случае верхние индексы координаты будут определять номер элемента и номер его локального узла, а нижние индексы соответственно номер координатной оси и номер глобального узла. Обозначение рассматривает k-ый конечный элемент, его m-ый локальный узел, i-ую координатную ось, l -ый глобальный узел. Таким образом, представления , ,,  определяют координату в направлении первой координатной оси узла F, используя локальную и глобальную нумерацию узлов. Введенного обозначения будем придерживаться и для идентификации перемещений узловых точек.

Основным неизвестным при решении задач теории упругости является вектор перемещения .  Следуя идеологии метода конечного элемента, для каждого дискретного элемента вводится локальная аппроксимация определённой на области функции. В этом случае определяются локальные функции формы элемента , где n - номер элемента, i – номер локального узла элемента. Их количество соответствует числу узлов элемента. Функции формы элемента обладают одним важным свойством: ,  т.е. равна единице в узле, номер которого совпадает с номером функции формы и равна нулю в остальных узлах элемента. Таким образом, аппроксимация поля перемещения на элементе через функции формы и узловые значения перемещений определяется в виде [37,45]:

,                                              (2.1)

где по одинаковому индексу, расположенному вверху и внизу двух стоящих рядом переменных подразумевается суммирование.

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

,                                              (2.2)

где - постоянные i-ой функции формы для n-ого элемента.

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

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

                                            (2.3)

Решая систему (2.3) приходим к выражениям:

, , ,                 (2.4)

где .

Для остальных функций форм получаем следующие значения постоянных:

, , ,                 (2.5)

, , .                 (2.6)

Используя представления (2.1) и (2.2) найдем выражения деформаций для n-ого элемента:

,                                     (2.7)

,                                   (2.8)

,                          (2.9)

где по индексу i ведется суммирование .

Отметим, что в пределах симплекс-элемента деформации (2.7)-(2.9) принимают постоянные значения.

2.3. Равновесие конечного элемента

Предполагаем, что процесс деформирования квазистатический и изотермический [32,33]. Рассмотрим условие равновесия n-ого изолированного конечного элемента в виде вариационного принципа Лагранжа:

,                                                 (2.10)

где  - тензор напряжений;  - набла-оператор;  - площадь элемента; - контур, ограничивающий элемент;  - вектор внешней нагрузки. Левая часть соотношения (2.10) определяет работу внутренних напряжения элемента на возможных перемещениях, а правая часть (2.10) – работу внешней нагрузки по контуру элемента.

Рассмотрим правую часть соотношения (2.10):

,                       (2.11)

где по повторяющимся индексам подразумевается суммирование.

Введем следующее обозначение:

.                                                         (2.12)

Величина обозначает -ю проекцию узловой силы в -ом узле конечного элемента. Соотношение (2.12) позволяет при заданной внешней нагрузки и геометрии конечного элемента свести распределенную внешнюю нагрузку к узловым точкам элемента. С учетом определения (2.12) запишем соотношение (2.11) в следующей форме:

,                                                     (2.13)

где ; ; по соответствующим повторяющимся индексам производится суммирование; n – номер изолированного элемента.

Преобразуем работу внутренних напряжений конечного элемента в левой части соотношения (2.10):

.          (2.14)

Из соотношения (2.10) с учетом (2.13), (2.14) приравнивая члены при равных вариациях, приходим к уравнениям равновесия конечного элемента:

.                                           (2.15)

Запишем (2.15) в виде соответствующей системы из шести уравнений:

                                     (2.16)

Представим связь между напряжениями и деформациями  в виде закона Гука:

,                                 (2.17)

,                                     (2.18)

,                                                       (2.19)

где , ,  для случая плоского деформирования; , ,  для случая плоского напряженного состояния.

Подставим в (2.17)-(2.19) выражения для деформаций в конечном элементе (2.7)-(2.9):

,                                 (2.20)

,                                     (2.21)

.                                                       (2.22)

Система (2.16) с учетом (2.20)-(2.22) предстанет в виде:

          (2.23)

Отметим, что подынтегральные выражения являются постоянными и могут быть вынесены за знак интеграла. Площадь конечного элемента обозначим через .  Перепишем первое уравнение системы (2.23) сгруппировав слагаемые перед соответствующими узловыми перемещениями:

        (2.24)

Из аналитической геометрии известно, что площадь треугольника на плоскости может быть вычислена  как определитель матрицы, составленный из координат вершин треугольника:

,

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

По аналогии с представлением (2.24) перепишем систему (2.23):

 (2.25)

Запишем связь между узловыми перемещениями и узловыми силами для n-ого конечного элемента в матричном виде:

,             (2.26)

где ; ;

      ; ;

      ; ;

           ; ;

       ; ;

       ; ;

           ; ;

       ; ;

       ; ;

           ; ;

       ; ;

       ; ;

           ; ;

       ; ;

       ; ;

         ; ;

      ; ;       

По сложившейся терминологии матрица  носит название матрицы жесткости элемента [37,45]. При заданном векторе узловых сил из уравнения (2.26) может быть найден вектор узловых перемещений. Отметим, что в узловой точке могут быть заданы  узловые перемещения, а так же узловая сила в одном направлении и узловое перемещение в направлении ему ортогональном. В этих случаях необходимо преобразовать систему (2.26), перенеся неизвестные компоненты узловых сил в левую часть системы, а вычисляемые компоненты левой части в столбец свободных членов.  При задании узловых компонент для конечного элемента необходимо помнить,  что граничные условия должны обеспечивать равновесие конечного элемента.

2.4. Равновесие системы конечных элементов

Как правило, реальные конструкции состоят более чем из одного элемента. В этом случае условие равновесия конечного элемента (2.15) переносится на каждый элемент из ансамбля конечных элементов, образующих деформируемое тело. Таким образом, равновесие тела определяется следующим уравнением:

,                              (2.28)

где  - количество элементов, образующих тело;- площадь -ого конечного элемента.

Однако, напрямую применение (2.28) нецелесообразно, т.к. для каждого элемента имеем шесть уравнений, и необходимо систему (2.28) дополнить связями на локальные узловые перемещения и узловые силы. Поэтому применяют переход от локальной нумерации перемещений и узловых сил к глобальной. Рассмотрим соответствующее преобразование для правой части системы (2.28):

,                                                   (2.29)

где - количество узлов в глобальной модели; .

В (2.29) под проекцией глобальной узловой силы   подразумевается сумма всех локальных проекций узловых сил принадлежащих данному узлу. В качестве примера рассмотрим конечноэлементную модель, показанную на рисунках 1 и 2. Глобальный узел 2 является общим для трех элементов. Для первого элемента он совпадает с первым локальным узлом, для второго – с его первым локальным узлом, а для третьего элемента – с третьим локальным узлом. Следовательно,  . Для остальных глобальных узлов рассмотренной модели по аналогии запишем: ; ; ; . Отметим, что если глобальный узел является внутренним, то .

Рассмотрим левую часть равенства (2.28). В силу представления (2.26) запишем:

,                                   (2.30)

где ,,,,, представляют собой локальные перемещения -ого конечного элемента.  

В силу непрерывности вектора перемещений в точке связи конечных элементов запишем правую часть соотношения (2.30) через глобальные узловые перемещения:

,                  (2.31)

где  компоненты вектора глобальных перемещений; связаны к с глобальными узловыми перемещениями следующим образом:  ,  - номер координатной оси, - глобальный номер узла;  - локальная матрица жесткости элемента в котором вершины имеют глобальную нумерацию.

Локальная матрица жесткости  имеет размерность . Компоненты матрицы равны нулю за исключением тридцати шести компонент соответствующих матрице , образуемой при использовании локальной нумерации узлов (26).

Из представления (31) получаем связь между матрицами жесткости элементов и глобальной матрицей жесткости:

.                                                       (2.33)

Таким образом, приходим к следующему алгоритму сборки матрицы жесткости:

  1.  Триангулировать заданную область. Ввести локальную и глобальную нумерацию узлов элементов.
  2.  Для каждого элемента определить локальные функции формы (2.2)
  3.  Определить коэффициенты локальной матрицы жесткости  для каждого элемента.
  4.  Зная количество глобальных узлов ,  сформировать  матрицу  размером .
  5.  Преобразовать локальные матрицы жесткости  в матрицы жесткости  размерности .
  6.  Провести суммирование матриц жесткости элементов согласно формуле (2.33).

После формирования глобальной матрицы жесткости необходимо определить вектор узловых сил с компонентами  размерностью . При этом нужно учесть, что для некоторых узловых точек, в силу граничных условий, компоненты узловых сил могут являться искомыми величинами. Таким образом, согласно (2.28), (2.29) и (2.31) приходим к следующей системе линейных уравнений:

.                                                           (2.34)

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

  1.  Известная компонента вектора узлового перемещения, например стоящая в -ой строке системы (2.34), умножается на коэффициент матрицы жестокости  и переносится с обратным знаком в столбец свободных членов для соответствующих уравнений. Получаемый таким способом столбец свободных членов обозначим через .
  2.  Неизвестная компонента вектора узловой силы в данном случае занимает -ую позицию в векторе столбце . В матрице жесткости  -ый столбец обнуляется за исключением -ого элемента, где ставится (-1). Модифицированную матрицу жесткости назовем .

В итоге решение задачи сводится к решению следующей системы уравнений:

.                                                   (2.35)

После нахождения поля узловых перемещений через функции формы определяются компоненты линейного тензора деформаций согласно (2.7)-(2.9). Затем, используя определяющие соотношения (2.20)-(2.22), строится напряженное состояние элементов конструкции.

2.5. Ошибки метода конечных элементов

Наряду с обычными ошибками округления и погрешностью приближенных методов линейной алгебры, применяемых в МКЭ, есть и ошибки, имеющие непосредственное отношение к методу конечных элементов:

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

 ошибки аппроксимации, обусловленные разностью между действительным распределением искомых функций в пределах КЭ и их представлением с помощью аппроксимирующих функций.

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

1) непрерывность искомой функции и ее производных в области КЭ до степени m–1 включительно (m – наибольший порядок производных искомой функции, используемых в качестве основных неизвестных в эрмитовых элементах);

2) выполнение условий полноты, т. е. при уменьшении размеров КЭ аппроксимирующие функции должны обеспечить стремление значений искомой функции, а также ее производных к постоянным значениям;

3) выполнение условий совместности искомой функции и частично ее производных на границе между смежными элементами;

4) приближенное удовлетворение условий совместности не основных переменных (например, напряжений, если основные неизвестные – перемещения) на границах КЭ, а также граничных условий в рассматриваемой области;

5) исключение концентрации напряжений в КЭ, если в рассматриваемой области такие концентрации заведомо отсутствуют;

6) при перемещениях КЭ как жесткого целого в нем не должны возникать деформации.

Требование полноты аппроксимирующих функций необходимо для учета смещения КЭ как жесткого целого и обеспечения состояния постоянных деформаций в элементе. Механический смысл совместности заключается в непрерывности основных неизвестных на смежных границах соседних КЭ.

  1.  Результаты расчетов
    1.  Использование метода конечных элементов для решения задачи.

Рассмотрим образец на Рис.1.2. Используя представления (2.17)-(2.19),(1.5),(1.8) запишем (1.25) и (1.26) в следующем виде:

   (3.1)

 (3.2)

Рассмотрим первое слагаемое левой части (3.1). С учетом (2.14), (2.28) и (2.31) имеем:

.        (3.3)

Пусть при конечно-элементном разбиении граница со слоем FT представлена совокупностью k1 граней k1 конечных элементов:

,                                                   (3.4)

где  - длина грани n-го конечного элемента, лежащей на сегменте  FT с координатами  и . Для грани  аналогично:

.

Второе слагаемое левой части (3.1) на  n-ом элементе запишем следующим образом:

                       (3.5)

,           (3.6)

.              (3.7)

Обозначим номера локальных узлов, лежащих на грани FT n-го элемента через d и c, а соответствующую им глобальную нумерацию  - значениями p, k. Для локальных узлов, лежащих на грани  n1-го элемента через d1 и c1, а соответствующую им глобальную нумерацию  - значениями p1, k1. Тогда из выражений (3.6) и (3.7) получим:

    (3.8)

    (3.9)

В правых частях (3.8) и (3.9) по переменным d,d1 и c,c1 суммирование не ведется. Используем в (3.8) и (3.9) глобальную нумерацию узлов:

 (3.10)

 (3.11)

где матрицы  и  размера  имеют все нулевые элементы за исключением:

; ;

; ;

;;

;.

Введем следующее обозначение:

,                                                 (3.12)

где суммирование ведется по k1 элементам, грани которых образуют границу FT.

Запишем третье слагаемое левой части (3.1) на  n-ом элементе:

,            

 

Используя глобальную нумерацию узлов:

 

По аналогии с матрицами  и  для третьего слагаемого левой части (3.1) вводим матрицу и , для которых запишем элементы, отличные от нуля:

; ;

; ;

; ;

; ,

где под верхним штрихом понимаем первую производную по переменной .

Определяем матрицу :

.                                                (3.13)

Для пятого слагаемого в левой части (3.1) на  n-ом элементе:

.

,

.

Используя глобальную нумерацию узлов:

Вводим в рассмотрение матрицы  и , со следующими ненулевыми компонентами:

; ;

; ;

; ;

; ,

и, соответственно,  матрицу :

.                                                (3.14)

Для седьмого слагаемого левой части (3.1) определяем матрицы  и , с ненулевыми компонентами:

; ;

; ;

; ;

; ,

и матрицу :

.                                (3.15)

Для второго, четвертого, шестого и восьмого слагаемых будем иметь:

; ;

; ;

;;

;.

;

; ;

; ;

; ;

; ,

;

; ;

; ;

; ;

; ,

; ;

; ;

; ;

; ,

Таким образом, матрица жесткости для уравнения (3.1) определяется в виде:

,                                                    (3.16)

где .

Рассмотрим уравнение (3.2). Первый элемент уравнения (3.2) будет находиться аналогично первому элементу уравнения (3.1).

Второе слагаемое левой части (3.2) на n1-ом элементе запишем следующим образом:

                       (3.17)

,           (3.18)

.              (3.19)

Из выражений (3.18) и (3.19) получим:

    (3.20)

    (3.21)

Используем в (3.20) и (3.21) глобальную нумерацию узлов:

 (3.22)

 (3.23)

где матрицы  и  размера  имеют все нулевые элементы за исключением:

; ;

; ;

;;

;.

Введем следующее обозначение:

,                                                 (3.24)

где суммирование ведется по k1 элементам, грани которых образуют границу .

Запишем третье слагаемое левой части (3.2) на  n1-ом элементе:

,            

 

Используя глобальную нумерацию узлов:

 

По аналогии с матрицами  и  для третьего слагаемого левой части (3.2) вводим матрицу и , для которых запишем элементы, отличные от нуля:

; ;

; ;

; ;

; ,

где под верхним штрихом понимаем первую производную по переменной .

Определяем матрицу :

.                                                (3.25)

Для пятого слагаемого в левой части (3.2) на  n1-ом элементе:

.

,

.

Используя глобальную нумерацию узлов:

Вводим в рассмотрение матрицы  и , со следующими ненулевыми компонентами:

; ;

; ;

; ;

; ,

и, соответственно,  матрицу :

.                                 (3.26)

Для седьмого слагаемого левой части (3.2) определяем матрицы  и , с ненулевыми компонентами:

; ;

; ;

; ;

; ,

и матрицу :

.                                (3.27)

Для второго, четвертого, шестого и восьмого слагаемых будем иметь:

; ;

; ;

;;

;.

;

; ;

; ;

; ;

; ,

;

; ;

; ;

; ;

; ,

; ;

; ;

; ;

; ,

Таким образом, матрица жесткости для уравнения (3.2) определяется в виде:

,                                                    (3.28)

где .

Глобальная матрица жесткости будет равна:

,                                                    (3.29)

а вариационные уравнения (3.1) и  (3.2)  сводятся к системе линейных уравнений:

.                                                       (3.30)

Рассмотрим функции формы по элементам, грани которых образуют границу FT и :

,                                         (3.31)

,                                         (3.32)

,                                         (3.33)

.                                         (3.34)

С учетом (3.31)-(3.34) определим соответствующие коэффициенты:

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

;

.

Коэффициенты для формулы (3.2) получаем аналогичным образом.

  1.  Моделирование трещины нормального отрыва.

В этом случае в точках  и  приложим равные по модулю и противоположные по направлению силы.  Точки  и закрепим от перемещений. Таким образом, для образца на Рис.1.2 имеем следующие граничные условия:

для точки ;

для точки ;

для точки ;

для точки ;

вся остальная поверхность тела свободна от напряжений.

В силу симметрии задачи достаточно рассмотреть половинку тела 1. В этом случае имеют место следующие условия для границ слоя:

,                                                             (3.35)

,                                                           (3.36)

С учетом (3.35) и (3.36) из (1.8) в слое .  Следовательно из (1.25) получаем:

       (3.37)

Рис.3.1. Распределение напряжений в слое в состоянии плоской деформации

На Рис.3.1 приведены значения средних напряжений по длине элемента слоя взаимодействия. Напряжения на графиках отнесены к максимальному касательному напряжению .  Координата  определяет расстояние от вершины трещины, график 1 - напряжение , график 2 - напряжение , график 3 - напряжение . Из приведенной зависимости видно, что в окрестности трещины нормального отрыва при плоской деформации  реализуется высокое гидростатическое растяжение.

Рис.3.2. Форма пластической зоны.

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

Рис.3.3. Распределение напряжений в слое в плоском напряженном состоянии.

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

  1.  Моделирование трещины типа моды II.

Для рассматриваемого модельного представления трещины нагружение по моде II должно приводить к следующему напряженному состоянию в образце Рис.1.2: средние напряжения в слое взаимодействия должны удовлетворять условиям: ; ; . В областях 1 и 2 вне слоя взаимодействия должна наблюдаться симметрия по напряжениям и антисимметрия по напряжениям ,. Добиться соответствующего напряженного состояния возможно приложив к точкам  и  равные по модулю и противоположно направленные горизонтальные силы.  При этом грани  и  должны быть закреплены от вертикальных перемещений. Для образца на Рис.1.2 имеем следующие граничные условия:

для точки ;

для точки ;

для грани ;

для грани ;

вся остальная поверхность тела свободна от напряжений.

В данном случае в силу антисимметрии  нагружения:

,                                                             (3.38)

,                                                             (3.39)

      (3.40)

На Рис.3.4 приведены значения среднего напряжения  для слоя  в концевой зоне трещины.

Рис.3.4. Распределение напряжений в слое.

Рис.3.5. Форма пластической зоны для плоского деформированного состояния.

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

Заключение

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

Список литературы

  1.  Аэро Э.Л. Существенно нелинейная микромеханика среды с изменяемой периодической структурой // Успехи мех. 2002. - №3. - С. 13-176.
  2.  Аэро Э.Л., Кувшинский Е.В. Континуальная теория асимметрической упругости. Равновесие изотропного тела // ФТТ. 1964. - Т. 6, вып.9. -С. 2689-2699.
  3.  Баренблатт Г.И., Христианович С.А. О модуле сцепления в теории трещин // Инженерный журнал. Механика твердого тела. 1968. № 2. С. 69-75.
  4.  Белов, П. А., С. А. Лурье. "Общая теория дефектов сплошных сред." Механика композиционных материалов и конструкций. 9.4 2003 C: 210-222.
  5.  Буханько А. А., Хромов А. И.  Пластическое течение в окрестности вершины трещины. Энергетический критерий разрушения и его связь с J-интегралом // Прикладная механика и техническая физика. 2012 №6 С. 112-120
  6.  Гаврилкина М.В., Глаголев В.В., Маркин А.А. К решению одной задачи механики разрушения // Прикладная механика и техническая физика. – Т. 48.– №4. – 2007. – С. 121-127.
  7.  Глаголев В.В., Кузнецов К.А., Маркин А.А. Модель процесса разделения деформируемого тела // Изв. РАН. МТТ. 2003. № 6. С.61-68.
  8.  Глаголев В.В., Маркин А.А. Об одной постановке задачи упругопластического разделения //  Прикладная механика и техническая физика. – 2009. – Т. 50 – №4. – С. 187-195.
  9.  Глаголев В.В., Маркин А.А. Модель установившегося разделения материального слоя // Изв. РАН. МТТ. 2004. № 5. С. 121-129.
  10.  Глаголев В.В., Маркин А.А. Об одном способе определения связей между критическими значениями характеристик процесса установившегося разделения материала // Проблемы прочности. 2006. №2. С. 47-58.
  11.  Глаголев В.В., Маркин А.А. Определение термомеханических характеристик процесса разделения // Известия РАН. Механика твердого тела. - №6. - 2007. – С. 101-112.
  12.  Глаголев В.В., Глаголев Л.В., Кунашов Н.Д. Продольный сдвиг в рамках дискретного подхода к разрушению// Вестник Чувашского государственного педагогического университета. Серия: Механика предельного состояния. – №4(14) – 2012. – С. 17-25.
  13.  Глаголев В.В., Маркин А.А. Нахождение предела упругого деформирования в концевой области физического разреза при произвольном нагружении его берегов // Прикладная механика и техническая физика. – 2012. – Т. 53 – №5. – С. 174-183.
  14.  Гольдштейн Р.В., Осипенко Н.М. Структуры разрушения в окрестности макроразрыва продольного сдвига  // Известия РАН. Механика твердого тела 2012 № 5 Страницы 22-34
  15.  Гольдштейн Р.В., Осипенко Н.М. Инициирование разрушения на контакте при сдвиге // Известия РАН. Механика твердого тела. 2013. № 4. С. 72-79.
  16.  Гольдштейн Р.В., Перельмутер М.Н. Рост трещин по границе соединения материалов // В кн.: Проблемы механики. Сб. статей. – М.: ФИЗМАТЛИТ. –  2003. –  С. 221-239.
  17.  Гольдштейн Р.В., Сарычев М.Е. Влияние дислокаций на критерий роста трещин по границе соединения деформируемых материалов // Известия РАН. Механика твердого тела. 2006. № 1. С. 125-135.
  18.  Гольдштейн Р.В., Шаталов Г.А. Моделирование процессов разрушения в рамках обобщенной модели атомистической трещины нормального отрыва // Известия РАН. Механика твердого тела. 2006. № 4. С. 151-164.
  19.  Гольдштейн Р.В., Шифрин Е.И. Овозможной искривления трещины нормального отрыва в анизотропной плоскости // Известия РАН. Механика твердого тела. 2006. № 6. С. 173-182.
  20.  Гольдштейн Р.В., Шифрин Е.И. О возможной неустойчивости прямолинейного пути трещины в ортотропной плоскости в условиях одноосного нормального растяжения // Известия РАН. Механика твердого тела. 2007. № 3. С. 33-45.
  21.  Качанов Л.М. Основы теории пластичности. – М Гостехиздат, 1956. 324 c.
  22.  Клевцов Г.В., Ботвина Л.Р. Микро- и макрозона пластической деформации как критерии предельного состояния материала при разрушении // Проблемы прочности. – 1984. – №4. – С. 24-28.
  23.  Корнев В.М. Обобщенный достаточный критерий прочности. Описание зоны предразрушения // ПМТФ. 2002. Т. 43, № 5. С. 153-161.
  24.  Корнев В.М. Распределение напряжений и раскрытие трещин в зоне предразрушения (подход Нейбера-Новожилова) // Физическая мезомеханика. 2004. Т. 7. № 3. С. 53-62.
  25.  Корнев В.М., Кургузов В.Д. Многопараметрический достаточный критерий квазихрупкой прочности для сложного напряженного состояния // Физическая мезомеханика. 2006. Т. 9. № 5. С. 43-52.
  26.  Лавит И.М. Об устойчивом росте трещины в упругопластическом материале // Проблемы прочности.-1988.-№7.-С.18-23.
  27.  Лавит И.М., Толоконников Л.А. Силы сцепления и J-интеграл // Изв. Сев.-Кавказского научного центра высш. Школы. Естественные нуки.-1985.-№1.-С.28-30.
  28.  Лавит И.М.  Рост трещины в условиях квазихрупкого разрушения при монотонно возрастающей и циклической нагрузках // Известия РАН. Механика твердого тела. – 2001 - № 2 – С. 109-120.
  29.  Лебедев А.А., Чаусов Н.Г. Феноменологические основы оценки трещиностойкости материалов по параметрам спадающих участков диаграмм деформаций//Пробл. прочности. – 1983. №2. С. 6-10.
  30.  Макклинток Ф. Пластические аспекты разрушения // Разрушение. Т.3. М.: Мир, 1975. С. 67-262.
  31.  Макклинток Ф.А., Ирвин Дж. Р. Вопросы пластичности в механике разрушения. В. кн.: Прикладные вопросы вязкости разрушения. 1968. С.143-186.
  32.  Маркин А.А. Термомеханика сплошных сред Учеб. пособие. Тула: Изд-во ТулГУ, 2009. - 106 с.
  33.  Маркин А.А, Соколова М.Ю., Христич Д.В. Процессы упругопластического конечного деформирования. Монография. Тула: Изд-во ТулГУ, 2011. - 374 с
  34.  Махутов Н.А. Деформационные критерии разрушения и расчет элементов конствукций на прочность.-М.:Машиностроение, 1981.-270с.
  35.  Механика - от дискретного к сплошному/ [А.Н. Андреев и др.] - Новосибирск: Изд-во СО РАН, 2008. - 344 с.
  36.  Мир-Салим-заде М.В. Зарождение трещины в подкрепленной пластине // ПМТФ. 2007. Т. 48, № 4. С. 111-120.
  37.  Морозов Е.М., Никишков Г.П. Метод конечных элементов в механике разрушения. - М.: Наука. Гл. ред. физ.-мат. лит-ры, 1980. - 256 с.
  38.  Морозов Н.Ф., Петров Ю.В. Проблемы механики разрушения твердых тел – Спб.: Изд-во  С-Петербурского ун-та, 1997.- 132 с.
  39.  Назаров С.А., Паукшто М.В. Дискретные модели и осреднение в задачах теории упругости. – Л. , 1984. -  93 с.
  40.  Назаров С.А., Шпековиус-Нойгебауер М. Применение энергетического критерия разрушения для определения формы слабоискривленной трещины // ПМТФ. 2006. Т. 47, № 5. С. 119-130.
  41.  НейберГ.  Концентрация напряжений. ОГИЗ, Гостехиздат, 1947. – 204 с.
  42.  Новожилов В.В. О необходимом и достаточном критерии хрупкой прочности // ПММ. 1969. № 2. С. 212-222.
  43.  Новожилов В.В., Рыбакина О.В. О перспективах построения критерия прочности при сложном нагружении // Прочность при малом числе циклов нагружения / М.: Наука. – 1969. –  С. 71-80.
  44.  НоттДж.Ф. Основы механики разрушения. – М: Металлургия, 1978. – 256 с.
  45.  Дж. Конечные элементы в нелинейной механике сплошных сред. – М.: Мир. – 1976. – 464 с.
  46.  Панасюк В.В. Механика квазихрупкого разрушения материалов. – Киев: Наук. думка, 1991. – 416 с.
  47.  Петров Ю.В. О «квантовой» природе разрушения хрупких сред // Докл. АН. 1991. Т. 321. № 1. С. 66-68.
  48.  Черепанов Г.П. Механика хрупкого разрушения. М.: Наука, 1974. 640 с.
  49.  Barenblatt, G.I. The formation of equilibrium cracks during brittle fracture. General ideas and hypotheses. Axially-symmetric cracks, J. Applied Mathematics and Mechanics 23 (1959), 622-636.
  50.  Barenblatt G.I. On a model of small fatigue cracks // Eng. Fract. Mech. – 1987. – V.28. - №5/6. – P. 623-626.
  51.  Bolotin V.V., Lebedev V.I. Analytical model of fatigue crack growth retardation due to overloading // International Journal of Solids and Structures. – 1996. –  №9. – P. 1229-1242.
  52.  Goldstein R.V., Perelmuter M.N. Modeling of bonding at the interface crack // Internal J. of Fracture. – 1999. – V. 99. – №1-2. – P. 53-79.
  53.  GlagolevV.V., MarkinA.A. Stress-Strain State in Elastic Body with Physical Cut// World Journal of Mechanics - Vol. 3 - No. 7- 2013. - P. 299-306.
  54.  Griffith A.A. The phenomenon or rupture and flow in solids // Phil. Trans. Roy. Soc., Ser. A. – 1920. – V. 221. – P.163-198.
  55.  Griffith A.A. The theory of rupture // In: Proc. 1st Int. Congr. Appl. Mech.- Delft. – 1924. – P. 55-63.
  56.  Irwin G.R. Relation or stresses near a crack to the crack extension force // Proc. 9th Int. Congr. Appl. Mech.- Brussels. – 1957. – V. 8. – P. 245-251.
  57.  Irwin G.R. Analysis of stresses and stain near the end of a crack traversing a plate // J. Appl. Mech. 1958. – V. 24. –    № 3. – P. 361-364. (Discussion // J. Appl. Mech. 1958. –     V. 25. – № 2. – P. 299-303 ).
  58.  Irwin G.R. Plastic zone near a crack and fracture toughness. – 7th SamagoreArdance Materials Research Conference. – Syracuse: Syracuse Univ. Press, 1960.  
  59.  McClintock F.A. Ductile fracture instability in shear // J. Appl. Mech. – 1958. – V. 25. – P. 581-588.
  60.  Prandtl L. EinGedankenmodellfür den ZerreibvorgandspröderKörper // ZAMM Bd. 13. 1933. P. 129-133.
  61.  Rice J.R. The elastic-plastic mechanics of crack extension // Int. J. Fracture Mech. – 1968. – V. 4. – № 1. – P. 41-47.
  62.  Rice J.R. Some mechanics research topics related to the hydrogen embrittlement of metals // Corrosion. – 1976. –     V. 32. – № 1. – P. 22-26.
  63.  Rice J.R., Johnson M.A. The role of large crack tip geometry changes in plane strain fracture // Inelastic Behaviour in Solids. – New York: McGraw-Hill. – 1970. – P. 641-672.


X2

X1

0

EMBED Equation.3

1

2

3

EMBED Equation.3

EMBED Equation.3

1

2

3

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

F(2)

C(4)

A(1)

D(5)

1

2

3

B(3)

X2

X1

0

1(3)

2(1)

3(2)

3(3)

4(2)

2(1)

4(2)

5(1)

2(3)

1

2

3

EMBED Equation.3

EMBED Equation.3

1.5

1

0.5

0.0

25

50

75

100

0

1

2

3

EMBED Equation.3

1.5

1

0.5

0.0

25

50

75

100

0

1

2

EMBED Equation.3

EMBED Equation.3

EMBED Equation.3

 0.4

0.2

0.0

25

50

75

100

0

1 Работа выполнена при поддержке РФФИ (гранты №13-08-00134 А и №14-01-31138 мол_а)


 

А также другие работы, которые могут Вас заинтересовать

29058. Гражданские правоотношения: понятие, виды гражданских правоотношений. Субъекты и объекты гражданских правоотношений. Основания их возникновения 89.5 KB
  Иначе говоря специфические черты и признаки гражданских правоотношений предопределены особенностями самого гражданского права. Содержание гражданского правоотношения Содержание гражданского правоотношения составляют субъективные права и обязанности его участников. Юридические возможности как составные части содержания субъективного гражданского права называются правомочиями. При весьма большом разнообразии содержания субъективных гражданских прав можно обнаружить что оно является результатом разновариантных комбинаций трех правомочий: 1...
29059. Правосубъектность гражданина 49.5 KB
  Гражданская дееспособность определяется как способность гражданина своими действиями приобретать и осуществлять гражданские права создавать для себя гражданские обязанности и исполнять их ст. 21 ГК РФ можно определить содержание дееспособности которое понимается как предоставленная гражданину возможность реализации своей правоспособности собственными действиями и включает способность гражданина своими действиями: приобретать гражданские права; осуществлять гражданские права; создавать и исполнять гражданские обязанности; нести...
29060. Юридические лица: понятие, виды, создание и прекращение юридических лиц 77 KB
  Понятие юридического лица в ст. Законодатель убрал уточнение по имущественной обособленности на праве собственности праве хозяйственного ведения или оперативного управления Суть от этого особо не изменилась поскольку в обоих определениях указываются одни и те же признаки юридического лица. В данном определении указываются признаки юридического лица это такие внутренние присущие ему свойства каждое из которых необходимо а все вместе достаточны для того чтобы организация могла признаваться субъектом гражданского права. Правовая...
29061. Объекты гражданских прав 169.5 KB
  К объектам гражданских прав закон относит вещи включая деньги и ценные бумаги иное имущество в том числе имущественные права; работы и услуги; результаты интеллектуальной деятельности в том числе исключительные права на них интеллектуальная собственность; нематериальные блага ст. Таким образомобъекты гражданских прав подразделяются на материальные и нематериальные идеальные. К первой группе относят: вещи; работы и услуги а также их результаты имеющие овеществленный либо иной стоимостной эффект например ремонтные работы услуги по...
29062. Понятие и виды сделок. Требования, предъявляемые к форме сделок 245.5 KB
  Но поскольку субъектами сделок являются граждане и юридические лица муниципальные образования субъекты РФ государство иностранные граждане и иностранные юридические лица лица без гражданства а субъективные гражданские права и обязанности являются содержанием гражданского правоотношения то понятие сделки можно определить следующим образом: Сделка – это осознанное волевое и юридически значимое действие субъектов гражданского права прямо направленное на возникновение изменение и прекращении гражданских правоотношний. С помощью сделки...
29063. Недействительные сделки: понятие, виды, последствия недействительности 206.5 KB
  Недействительными являются сделки не создающие правового результата прав и обязанностей к которому стремились стороны. Эти основания связаны с теми или иными нарушениями условий действительности сделок а именно: законность содержания сделки соответствие воли и волеизъявления соблюдение требуемой по закону формы сделки совершение сделки лицом обладающим необходимой дееспособностью. Основания недействительности должны иметь место одновременно с совершением сделки как юридического факта.
29064. Исковая давность в гражданском праве 100 KB
  Институт исковой давности имеет большое значение для защиты имущественных прав граждан и юридических лиц. С его помощью во многом обеспечивается стабильность гражданского оборота: участники гражданских отношений знают что по истечении исковой давности никто не может предъявить к ним требований и провести ревизию принадлежащих им прав. С истечением срока исковой давности право на принудительную судебную защиту не прекращается а приобретает условный характер: оно сохраняется если сторона в споре не заявит о применении исковой давности. При...
29065. Право собственности: понятие, содержание, основания возникновения и прекращения 51 KB
  Нормы образующие институт права собственности находятся в постоянном контакте и взаимодействии с нормами других правовых институтов например с обязательственным правом. Субъективное право собственности это закрепленная за собственником юридически обеспеченная возможность владеть пользоваться и распоряжаться принадлежащим ему имуществом по своему усмотрению и в своих интересах путем совершения в отношении этого имущества любых действий которые не противоречат закону и иным правовым актам и не нарушают права и охраняемые законом интересы...
29066. Вещно – правовые средства защиты права собственности и других вещных прав 37 KB
  виндикационный об истребовании имущества из чужого незаконного владенияи негаторный об устранении препятствий в пользовании имуществом не связанных с лишением владения вещью. Виндикационный иск Собственник который оказался лишенным принадлежащего ему имущества без необходимого правового основания может истребовать свое имущество из чужого незаконного владения. Виндикационные иски предъявляются в отношении индивидуальноопределенного имущества причем сохранившегося в натуре; в иных ситуациях надлежит требовать возмещения причиненных...