ПЕРЕХОД К КОНТУРНОМУ J-ИНТЕГРАЛУ
 25 июля, 2014
 25 июля, 2014  Oleg Maloletnikov
 Oleg Maloletnikov В конструкциях дефекты могут иметь сложную форму. Кроме того, форма дефектных узлов конструкции может быть самой разнообразной. Поэтому для дефекта в конструкции действующее значение J-интеграла чаще всего находится на основе численных решений упругопластической задачи. В этом случае он вычисля - а ется как контурный интеграл по
| 
 Рис. 6.53 Контур интегрирования у вершины трещины и направления интегрирования по оси у | 
произвольному контуру S, включающему вершину трещины и связанную с ней пластическую зону, как показано на рис. 6.53а.
| 
 | 
| (V) | 
| 
 | 
| (S) | 
| (6.127) | 
| 
 | 
| Считается, что нагрузка (, приложенная к контуру, не меняется при росте трещины. Этот интеграл берется в неподвижной системе координат х, у, показанной на рис. 6.53а. Второй интеграл берется по объему, | 
В основе определения контурного интеграла лежит выражение (6.121), связывающее J с механическим потенциалом П. Перенося дифференцирование под знак интеграла, получим
заключенному внутри контура S, по которому берется первый интеграл.
Но поле напряжений и перемещений у вершины трещины описывается в подвижной системе координат X, Y, центр которой расположен в вершине трещины. Поле усилий q остается постоянным только в подвижной системе координат.
Связь подвижной системы координат с неподвижной системой дается формулами:
X = x - l; Y = y, откуда можно вычислить производные:
X = _ 0; dl, ;
it=+1,o;
dx
(6.128)
| d dl | 
| d ' dl | 
| d dl | 
| d dX | 
| d dl | 
| d dx' | 
| d dX dX dl | 
Используем оператор последней строки формул (6.128) для дифференцирования по длине трещины в формуле (6.127):
| 1 • | f q ■4й • ds - . dV 1 | -1. | fa •4u • ds - • dV | 
| t | J 4 4l j 4l (S) (V) | t | 1 x 4 x 4 q I | 
| . (6.129) | 
Докажем, что первая квадратная скобка в этом выражении равна нулю.
Если тело находится в равновесии, то согласно принципу возможных перемещений вариация механического потенциала 8П на возможных перемещениях равна нулю:
| (6.130) | 
8П = J 8w ■ dV - J q -8й ■ ds = 0,
| (S) | 
(V)
где 8w и 8й — вариации удельной энергии и вариации перемещений, которые могут быть заданы произвольно, но должны быть связаны друг с другом уравнениями теории упругости или теории пластичности (уравнения сплошности, равновесия, связи между деформациями и напряжениями, и т. п.). В частности, можно назначить:
8w =^W,§l; dl
8й = du -8l.
dl
Если эти значения вариаций подставить в уравнение (6.130) и сократить его на 81, то левая его часть не будет отличаться от первой скобки формулы (6.129). Следовательно, для тела, находящегося в равновесии:
| 
 | 
(6.131)
Остается преобразовать первый интеграл, взятый по объему, в контурный:
| 
 | 
(6.132)
| 
 | 
Последнее равенство в (6.132) ясно из рис. 6.53б. Точки, в которых вычисляются w(xj) и w(x2), лежат на контуре. При увеличении у от ушіп до ушах в процессе интегрирования по у, точка х2 будет перемещаться по правой стороне контура против часовой стрелки, а точка х1 — по левой по часовой стрелке. Знак «минус» в контурном интеграле появится в связи с изменением угла между направлениями dy и ds при переходе контурного интеграла из одной половины контура в другую.
| 
 | 
Остается только подставить последнее выражение (6.132) в (6.131). Тогда при единичной толщине получим:
(6.133)
Если получено распределение напряжений, деформаций и перемещений по сетке конечных элементов и намечен наиболее удобный для вычисления интеграла контур S вокруг вершины трещины, то используя формулу (6.133), нетрудно написать алгоритм для вычисления J-интеграла.
Условие прочности (отсутствие старта трещины) проверяют по формуле:
| J < Jc, | 
(6.134)
где Jc — вязкость разрушения этого материала, полученная экспериментально по ГОСТ 25.506-85.









 Опубликовано в
 Опубликовано в