Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

THE CHOICE OF THE DETERMINISTIC MODEL WHEN THE CALCULATION DETAILS OF THE FINITE ELEMENT METHOD

Gots A.N. 1
1 Vladimir State University named after Alexander and Nikolay Stoletovs
In the calculation of machine parts using finite element method there is a difficulty in the choice of a deterministic model in determining the factors of safety. All software systems that use finite element method, as a deterministic model for a known stress tensor in the vicinity of some nodal points for the determination of amplitude and mean stress is used stress intensity defined by the maximum and minimum loading. However, this does not take into account the fact that areas where stress intensity is determined, under the action of the maximum and minimum loads will be different, and therefore the calculation is not quite correct. To determine if the intensity of the amplitudes of the alternating stresses calculated in the same sites, but as a medium stress to choose the constant stress cycle, the results of calculations of the strength margins will differ from the calculation of stress intensity.
stress tensor
a deterministic model
amplitude stress
mean stresses
stress intensity
factor of safety
1. Belov E.S. Raschet sostavnogo blok-kartera odnocilindrovogo dizelya /E.S. Belov, B.Yu. Golev, A.N. Goc i dr. // Fundamentalnye i prikladnye problemy sovershenstvovaniya porshnevyx dvigatelej: Materialy XII Mezhdun. nauchno-prakt. konferencii; pod red. A.N. Goca. VlGU. Vladimir. 2010. рр. 114–118.
2. Goc A. N. Raschet shheki kolenchatogo vala DVS / A.N. Goc, I. N. Kudelya // Traktory i selskoxozyajstvennye mashiny. 1998. no. 9. рр. 26–31.
3. Goc A.N. Determinirovannye modeli ustalostnoj dolgovechnosti detalej DVS / Traktory i selskoxozyajstvennye mashiny. 1997. no. 12. рр. 19–21.
4. Goc A.N. Raschety na prochnost detalej DVS pri napryazheniyax, peremennyx vo vremeni. 3-e izd., isp. i dop. M.: FORUM: infra-m, 2013. 208 р.
5. Goc A.N. Vybor determinirovannoj modeli ustalostnoj dolgovechnosti / A.N. Goc, A.B. Ivanchenko // Aktualnye problemy upravleniya kachestvom proizvodstva i ekspluatacii transportnyx sredstv: materialy IX mezhdunar. nauch.-prakt. konf. / Vladim. gos. un-t, Vladimir, VlGU, 2002. рр. 311–314.
6. Goc A.N. Vybor determinirovannoj modeli ustalostnoj dolgovechnosti pri raschete detalej porshnevyx dvigatelej na peremennye nagruzki // Sb. nauchnyx trudov Mezhdunarodnoj konferencii Dvigatel-2010, posvyashhennoj 180-letiyu MGTU im. N.E. Baumana // Pod red. N.A. Ivashhenko, V.A. Vagnera, L.V. Grexova. M.: MGTU im. N.E. Baumana, 2010. рр. 44–49.
7. Goc A.N. Raschety na prochnost detalej DVS pri napryazheniyax, peremennyx vo vremeni: ucheb. posobie. 2-e izd., ispr. i dop.; Vladim. gos. un-t, Vladimir: Izd-vo Vladim. gos. un-ta, 2011. 140 р.
8. Goc A.N. Opredelenie zapasa prochnosti shheki kolenchatogo vala s polozhitelnym perekrytiem / A.N. Goc, I.N. Kudelya // Dvigatelestroenie. 1999. no. 4. рр. 15–18.
9. Determinirovannye i statisticheskie modeli dolgovechnosti / I.A. Birger // Problemy prochnosti letatelnyx apparatov: sb. st.; pod red. I.F. Obrazcova, A.S. Volmira. M.: Mashinostroenie. 1985. 280 р.
10. Ponomarev S.D. Raschety na prochnost v mashinostroenii / S.D. Ponomarev [i dr.]. V 3 t. T. 3 / pod red. S.D. Ponomareva. M.: Mashgiz. 1959. 1118 р.

При расчете деталей поршневых двигателей методом конечных элементов (МКЭ) возникает вопрос о выборе детерминированных моделей усталостной долговечности. Как правило, при использовании наиболее распространенных программных комплексов расчета МКЭ в качестве гипотезы разрушения при сложном напряженном состоянии используется гипотеза удельной энергии формоизменения (гипотеза Губера – Мизеса – Генки), в которой эквивалентное напряжение определяется через интенсивность напряжений (или октаэдрическое касательное напряжение) по формуле [4]

gots01.wmf (1)

или через главные напряжения

gots02.wmf (1, а)

где σx, σy, σz, τxy, τyz, τzx – компоненты тензора напряжений в окрестности узловой точки конечно-элементной модели:

gots03.wmf (2)

σ1, σ2 и σ3 – главные напряжения, определяемые при том же тензоре напряжений.

Известно [10], что результаты с использованием этой гипотезы достаточно хорошо согласуются с экспериментальными данными не только в области упругих деформаций, но и в пластической области в случае статического нагружения. Это дало повод использовать интенсивность напряжений при определении запасов прочности или долговечности деталей ДВС в случае нагружения их переменной нагрузкой. Другими словами, при переменном действии нагрузок для расчета интенсивности напряжений σi по формуле (1) в качестве исходного выбирают сначала максимальное нагружение, а затем минимальное, что позволяет найти σimax, а затем σimin, а далее оперировать ими, как обычно, как максимальными и минимальными напряжениями. Однако полученные при этом σimax и σimin получены для разных площадок, так как главные площадки могут быть во взаимно перпендикулярных площадках. Это было отмечено нами ранее в случае расчета даже для плоского напряженного состояния [4].

При действии сначала максимальных, а затем и минимальных нагрузок компоненты тензора напряжений (2) в окрестностях каждой точки конечно-элементной модели будут иметь разные значения. Значит, будут разными для каждого вида нагружения главные напряжения и положения площадок, в которых они действуют. Действительно, если при статическом нагружении определяется эквивалентное напряжение для какой-то определенной площадки, положение которой нас даже не интересует, то при использовании (1) в случае переменных нагрузок напряжения σimax и σimin будут определены для площадок, никак не связанных между собой (в некоторых случаях они могут быть взаимно перпендикулярными). Это значит, что запасы прочности при переменных напряжениях или запасы долговечности будут вычисляться, мягко говоря, не вполне корректно, поскольку амплитудные и средние напряжения определены с использованием напряжений, вычисление которых проводилось по формуле (1,а) в разных площадках.

Если графически интерпретировать прогнозирование разрушения по гипотезе удельной энергии формоизменения, то поверхность разрушения представляет круговой цилиндр, ось которого образует равные углы с тремя главными осями (рис. 1). Эти круговые цилиндры при разных видах нагружения могут пересекать друг друга, так как соотношение главных напряжений (и положение взаимно перпендикулярных площадок) σ1 > σ2 > σ3 будет различным при максимальной и минимальной нагрузке. Круговые цилиндры располагались бы по типу «матрешки» и стороны их были бы эквидистантными, если бы циклы переменных напряжений при максимальной и минимальной нагрузках были бы подобными (с одинаковыми коэффициентами асимметрии). Естественно, что вероятность такого нагружения весьма мала.

pic_4.tif

Рис. 1. Графическое представление гипотезы удельной энергии формоизменения для случая произвольного напряженного состояния

И.А. Биргером была предложена детерминированная модель усталостной долговечности, применимая при расчете любых конструкций в условиях переменного нагружения [9] и проверенная нами при расчетном исследовании широкого круга деталей ДВС [1–8]. Для расчета эквивалентного амплитудного напряжения σea при сложном напряженном состоянии (при объемном тензоре напряжений в случае расчета МКЭ) и действии переменных нагрузок предложены зависимости:

эквивалентное напряжение

gots04.wmf (3)

где Kσ – эффективный коэффициент концентрации напряжений; εσk – масштабный коэффициент (влияния абсолютных размеров детали); βσ – коэффициент поверхностной чувствительности; ψσ – коэффициент, учитывающий влияние средних напряжений на предел выносливости [4];

gots05.wmf (4)

σia – интенсивность амплитуд переменных напряжений цикла σxa, σya, ..., τxya, τyza, ... (вычисленных по известному тензору напряжений [10]); σ1m – наибольшее постоянное напряжение цикла (первое главное напряжение.

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

Проведем сравнение расчетных данных с использованием (1) и (3) на примере расчета запасов прочности составного блок-картера одноцилиндрового дизеля [1]. На рис. 2 показана расчетная конечно-элементная модель блок-картера, которая состояла из 247028 элементов и 457782 узлов [4, 9]. Для расчета использовался программный комплекс Ansys 11.01 Наибольшие напряжения возникали в сечениях вокруг бобышки (см. рис. 2 – показаны стрелками), которая является коренной опорой. Максимальная нагрузка на блок-картер возникала при действии газовой силы и сил инерции деталей, совершающих возвратно-поступательное движение, а минимальная – только от действия сил инерции на режиме максимальной частоты вращения холостого хода [1].

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

В последних столбцах табл. 1 и 2 приведены значения эквивалентных напряжений, вычисленных по формуле (1), а на рис. 3 показаны точки конечно-элементной модели, где эквивалентные напряжения σi в соответствии с гипотезой потенциальной энергии формоизменения достигают максимального значения.

pic_5.tif

Рис. 2. Расчетная КЭ модель сборки блок-картера: узлов – 457782, элементов – 247028

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

Таблица 1

Компоненты тензоров напряжений при действии газовых и инерционных сил

№ п/п

Нормальные напряжения

Касательные напряжения

σe, МПа

σx, МПа

σy, МПа

σz, МПа

τxy, МПа

τyz, МПа

τxz, МПа

1

16,1

89,2

16,6

33,6

0,34

–4,7

93,6

2

18,8

89,6

16,3

37,0

8,0

9,4

98,8

3

19,6

64,4

11,4

–30,3

6,1

–5,6

73,5

4

39,8

103,6

43,1

27,8

–54,0

7,4

122,9

5

–21,6

66,1

26,8

6,7

–53,7

–25,7

128,7

6

–22,3

–112,2

–16,8

17,3

33

5,1

113,4

7

–11,1

–97,4

–12,9

–27,2

–7,4

–1,2

98,4

8

5,2

91,4

15,3

19,3

11

6

90,8

9

9,9

78,7

16,9

12,9

–24,9

1,1

81,6

Таблица 2

Компоненты тензоров напряжений при действии сил инерции

№ п/п

Нормальные напряжения

Касательные напряжения

σi, МПа

σx, МПа

σy, МПа

σz, МПа

τxy, МПа

τyz, МПа

τxz, МПа

1

0,93

43,6

0,93

1,8

–0,23

–0,36

42,8

2

1,2

5,8

0,9

2,4

0,5

0,5

6,4

3

4,2

16,5

2,4

–6,9

–1,0

–0,1

17,96

4

16,4

35,9

15,0

10,7

–17,9

3,6

41,9

5

0,85

–2,2

26,8

0,5

22,9

0,4

48,3

6

1,1

3,9

0,45

–1,03

–0,89

–0,17

3,97

7

0,36

0,28

–0,2

0,85

–0,2

–0,6

1,91

8

0,4

5,6

1,2

1,3

0,8

0,5

5,59

9

1

9,5

2,9

–1

–4,9

–0,1

11,6

pic_6.tif

Рис. 3. Значения эквивалентных напряжений

Все программные комплексы для расчета деталей ДВС методом МКЭ приводят в результате расчета значения интенсивности напряжений σi, которые и представляют собой значение эквивалентных напряжений, вычисленное по формуле (1). На рис. 3 приведена конечно-элементная модель блок-картера, а стрелками показаны эквивалентные напряжения при действии газовых и инерционных сил. В этих точках в эксплуатационных условиях появлялись трещины.

Мы уже отмечали, что эти напряжения являются эквивалентными в случае статического нагружения, но их использовать для расчета запасов прочности или долговечности при действии переменных напряжений (при максимальной нагрузке – как максимальные, а минимальной – минимальные) не вполне корректно, а кроме того, приводит к неверным результатам.

Для использования модели (3) проведем расчет амплитудных нормальных и касательных напряжений в соответствии с данными табл. 1 и 2. После этого определим амплитудные значения нормальных, касательных, а также интенсивность амплитуд напряжений переменных циклов в соответствии с зависимостью (3). Запасы прочности определим по зависимости [4]:

gots06.wmf (5)

Здесь σ–1 – предел выносливости при изгибе при симметричном цикле. Смысл значений Kσ, βσ, Ψσ, σ1m приведен выше.

Расчеты проводились при следующих исходных данных: K σ = 1,4; εσk = 0,9; βσ = 0,9; Ψσ = 0,2 [27];σ–1 = 120 МПа. Результаты расчетов приведены в табл. 10. Расчеты проводились с учетом того, что сначала определялись эквивалентные амплитудные напряжения, а второй расчет, когда амплитудные и средние напряжения определялись с учетом зависимости.

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

Таблица 3

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

№ п/п

Амплитуды нормальных напряжений

Амплитуды касательных напряжений

σae, МПа

σ1m, МПа

nσ

nσ

σxa, МПа

σya, МПа

σza, МПа

τxya, МПа

τyza, МПа

τxza, МПа

1

7,6

22,8

7,8

15,9

0,3

–2,6

31,7

104

1,59

2,08

2

8,8

41,9

7,7

17,3

3,7

4,8

46,3

106,8

1,18

1,33

3

7,7

24,0

4,5

–11,7

3,6

–4,6

29,0

81,4

1,81

2,1

4

11,7

33,9

14,1

8,6

–18,1

3,5

41,0

139,8

1,22

1,38

5

–11,2

34,2

0,0

3,1

–38,3

–12,8

81,2

107,8

0,74

1,38

6

–11,7

–58,1

–8,6

9,2

16,9

2,9

58,6

–1,2

1,19

1,13

7

–5,7

–48,8

–6,4

–14,0

–3,6

–0,9

49,6

–3,3

1,41

1,28

8

2,4

42,9

7,1

9,0

5,1

3,1

42,7

97,3

1,29

1,44

9

4,5

34,6

7,0

7,0

–10,0

0,6

35,8

89,4

1,50

1,72

Таким образом, для расчета на действия переменных напряжений необходимо внести в известные расчетные комплексы дополнительно расчет амплитудных межцикловых напряжений и амплитудного эквивалентного напряжения в соответствии с детерминированной моделью И.А. Биргера.

Рецензенты:

Кобзев А.А., д.т.н., профессор, ФГБОУ ВПО «Владимирский государственный университет имени Александра Григорьевича и Николая Григорьевича Столетовых», г. Владимир;

Кульчицкий А.Р., д.т.н., профессор, главный специалист ООО «Завод инновационных технологий» КТЗ, г. Владимир.

[1] Расчеты проводились Б.Ю. Голевым и Е.С. Беловым.