Расчет по циклической долговечности

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

КАЗАНСКИЙ (ПРИВОЛЖСКИЙ) ФЕДЕРАЛЬНЫЙ УНИВЕРСИТЕТ

ИНСТИТУТ МАТЕМАТИКИ И МЕХАНИКИ ИМ. Н.И. ЛОБАЧЕВСКОГО

Направление: Механика и математическое моделирование

Индивидуальная работа № 2.2 по дисциплине

“Физико-механический практикум и вычислительный эксперимент”

Выполнил: студент 4 курса группы 05-301

Герасимов О.В.

Проверил: доцент кафедры теоретической механики

Великанов П.Г

Казань 2016

Решение задачи расчета диска на прочностную надежность

Задание:

1. Построить плоскую модель сечения диска.

2. Задать:

1. нелинейное распределение параметров в табличном виде.

2. граничные условия на линии соприкосновения с валом диска: фиксирование перемещений вдоль оси OY.

3. распределенную нагрузку на большем радиусе.

4. угловую скорость и ускорение свободного падения.

3. Вывести:

1. эпюру распределения температуры.

2. деформированное и недеформированное состояние.

3. напряжения по критерию Вон-Мизеса.

4. эпюру перемещений по оси OX.

5. эпюру распределения напряжений.

Решение

Расчет по циклической долговечности - student2.ru

Рис.1 Плоская модель: область сечения.

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

Расчет по циклической долговечности - student2.ru

Рис.2 Напряжения и температуры.

Создадим сетку:

Расчет по циклической долговечности - student2.ru

Рис.3 Разбиение сетки.

Зададим граничные условия, распределенную нагрузку, угловую скорость и ускорения свободного падения:

Расчет по циклической долговечности - student2.ru

Рис.4 Задание необходимых условий.

Осуществим решение. Выведем все необходимые эпюры.

ANSYS выводит на экран графический анализ итерационного процесса решения задачи:

Расчет по циклической долговечности - student2.ru

Рис.5 Графический анализ итерационного процесса.

Расчет по циклической долговечности - student2.ru

Рис.6 Температурное распределение.

Расчет по циклической долговечности - student2.ru

Рис.7 Деформированное и недеформированное состояние.

Расчет по циклической долговечности - student2.ru

Рис.8 Напряжения по критерию Вон-Мизеса.

Расчет по циклической долговечности - student2.ru

Рис.9 Распределение перемещений по оси OX.

Расчет по циклической долговечности - student2.ru

Рис.10 Распределение напряжений по оси OX.

Расчет по циклической долговечности - student2.ru

Рис.11 Распределение напряжений по оси OY.

Расчет по циклической долговечности - student2.ru

Рис.12 Распределение перемещений по оси OY.

Исследование на собственные частоты.

Расчет по циклической долговечности - student2.ru

Рис.13.1 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис.13.2 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис.13.3 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис.13.4 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис.13.5 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис.13.6 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис.13.7 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Табл.1 Таблица значений критических частот.

LOG-FILE

finish

/CLEAR,START

/units,si

/COLOR,PBAK,OFF

/REPLOT

/prep7

!* Построение точек.

K,1,0.0767,0.071/2,

K,2,0.0805,0.0845/2,

K,3,0.0875,0.0888/2,

K,4,0.0945,0.0888/2,

K,5,0.109,0.084/2,

K,6,0.12,0.0638/2,

K,7,0.138,0.0382/2,

K,8,0.152,0.0251/2,

K,9,0.179,0.0231/2,

K,10,0.194,0.0334/2,

K,11,0.196,0.0334/2,

K,12,0.21,0.0334/2,

K,13,0.2125,0.03/2,

K,14,0.22,0.022/2,

K,15,0.2627,0.0185/2,

K,16,0.305,0.016/2,

K,17,0.348,0.0134/2,

K,18,0.349,0.0122/2,

K,19,0.357,0.015/2,

K,20,0.367,0.0228/2,

K,21,0.374,0.0235/2,

K,22,0.374,,

K,23,0.0767,,

!* Соединение точек.

L, 1, 2

L, 2, 3

L, 3, 4

L, 4, 5

L, 5, 6

L, 6, 7

L, 7, 8

L, 8, 9

L, 9, 10

L, 10, 11

L, 11, 12

L, 12, 13

L, 13, 14

L, 14, 15

L, 15, 16

L, 16, 17

L, 17, 18

L, 18, 19

L, 19, 20

L, 20, 21

L, 21, 22

L, 22, 23

L, 23, 1

!* Создание области.

FLST,2,23,4

FITEM,2,7

FITEM,2,6

FITEM,2,8

FITEM,2,5

FITEM,2,4

FITEM,2,3

FITEM,2,2

FITEM,2,9

FITEM,2,1

FITEM,2,10

FITEM,2,11

FITEM,2,23

FITEM,2,22

FITEM,2,12

FITEM,2,13

FITEM,2,14

FITEM,2,15

FITEM,2,16

FITEM,2,17

FITEM,2,18

FITEM,2,19

FITEM,2,20

FITEM,2,21

AL,P51X

!* Отображаем зеркально область относительно вертикальной оси.

FLST,3,1,5,ORDE,1

FITEM,3,1

ARSYM,Y,P51X, , , ,0,0

/PNUM,ELEM,0

/REPLOT

!* Обединяем отображение и оригинал.

FLST,2,2,5,ORDE,2

FITEM,2,1

FITEM,2,-2

AADD,P51X

!* Задаем тип материала.

ET,1,PLANE82

!* Определяемм параметры материала.

KEYOPT,1,1,0

KEYOPT,1,3,1

KEYOPT,1,6,1

!* Задаем табличные значения.

TB,MELA,1,5,15,

TBTEMP,473

TBPT,,0,0

TBPT,,0.1/100,1880e5

TBPT,,0.2/100,3800e5

TBPT,,0.3/100,5720e5

TBPT,,0.35/100,6420e5

TBPT,,0.37/100,6550e5

TBPT,,0.4/100,6720e5

TBPT,,0.42/100,6810e5

TBPT,,0.46/100,6970e5

TBPT,,0.5/100,7100e5

TBPT,,0.55/100,7200e5

TBPT,,0.6/100,7300e5

TBPT,,0.7/100,7390e5

TBPT,,0.8/100,7460e5

TBPT,,0.86/100,7500e5

TBTEMP,573

TBPT,,0,0

TBPT,,0.1/100,1830e5

TBPT,,0.2/100,3700e5

TBPT,,0.3/100,5580e5

TBPT,,0.35/100,6310e5

TBPT,,0.37/100,6450e5

TBPT,,0.4/100,6630e5

TBPT,,0.42/100,6720e5

TBPT,,0.46/100,6880e5

TBPT,,0.5/100,7000e5

TBPT,,0.55/100,7100e5

TBPT,,0.6/100,7190e5

TBPT,,0.7/100,7260e5

TBPT,,0.8/100,7330e5

TBPT,,0.86/100,7380e5

TBTEMP,773

TBPT,,0,0

TBPT,,0.1/100,1700e5

TBPT,,0.2/100,3420e5

TBPT,,0.3/100,5100e5

TBPT,,0.35/100,5800e5

TBPT,,0.37/100,6000e5

TBPT,,0.4/100,6230e5

TBPT,,0.42/100,6330e5

TBPT,,0.46/100,6490e5

TBPT,,0.5/100,6620e5

TBPT,,0.55/100,6720e5

TBPT,,0.6/100,6760e5

TBPT,,0.7/100,6830e5

TBPT,,0.8/100,6900e5

TBPT,,0.86/100,6930e5

TBTEMP,873

TBPT,,0,0

TBPT,,0.1/100,1600e5

TBPT,,0.2/100,3200e5

TBPT,,0.3/100,4830e5

TBPT,,0.35/100,5600e5

TBPT,,0.37/100,5800e5

TBPT,,0.4/100,6000e5

TBPT,,0.42/100,6110e5

TBPT,,0.46/100,6250e5

TBPT,,0.5/100,6370e5

TBPT,,0.55/100,6470e5

TBPT,,0.6/100,6500e5

TBPT,,0.7/100,6560e5

TBPT,,0.8/100,6600e5

TBPT,,0.86/100,6630e5

TBTEMP,973

TBPT,,0,0

TBPT,,0.1/100,2470e5

TBPT,,0.2/100,3000e5

TBPT,,0.3/100,4570e5

TBPT,,0.35/100,5200e5

TBPT,,0.37/100,5130e5

TBPT,,0.4/100,5500e5

TBPT,,0.42/100,5620e5

TBPT,,0.46/100,5850e5

TBPT,,0.5/100,6000e5

TBPT,,0.55/100,6200e5

TBPT,,0.6/100,6300e5

TBPT,,0.7/100,6350e5

TBPT,,0.8/100,6400e5

TBPT,,0.86/100,6420e5

GPLOT

!*

MPTEMP,,,,,,,,

MPTEMP,1,475

MPTEMP,2,477.88

MPTEMP,3,483.17

MPTEMP,4,488.47

MPTEMP,5,499.45

MPTEMP,6,507.77

MPTEMP,7,521.39

MPTEMP,8,531.99

MPTEMP,9,552.42

MPTEMP,10,563.77

MPTEMP,11,565.29

MPTEMP,12,575.88

MPTEMP,13,577.77

MPTEMP,14,583.45

MPTEMP,15,615.77

MPTEMP,16,647.78

MPTEMP,17,680.32

MPTEMP,18,681.08

MPTEMP,19,687.13

MPTEMP,20,694.70

MPTEMP,21,700

UIMP,1,REFT,,,

MPDATA,ALPX,1,,12.92e-7

MPDATA,ALPX,1,,13.04e-7

MPDATA,ALPX,1,,13.25e-7

MPDATA,ALPX,1,,13.46e-7

MPDATA,ALPX,1,,13.9e-7

MPDATA,ALPX,1,,14.23e-7

MPDATA,ALPX,1,,14.78e-7

MPDATA,ALPX,1,,15.2e-7

MPDATA,ALPX,1,,16.02e-7

MPDATA,ALPX,1,,16.47e-7

MPDATA,ALPX,1,,16.53e-7

MPDATA,ALPX,1,,16.96e-7

MPDATA,ALPX,1,,17.03e-7

MPDATA,ALPX,1,,17.26e-7

MPDATA,ALPX,1,,18.55e-7

MPDATA,ALPX,1,,19.83e-7

MPDATA,ALPX,1,,21.13e-7

MPDATA,ALPX,1,,21.16e-7

MPDATA,ALPX,1,,21.41e-7

MPDATA,ALPX,1,,21.71e-7

MPDATA,ALPX,1,,21.92e-7

!* Задание параметров разбиения.

FLST,2,1,5,ORDE,1

FITEM,2,3

AESIZE,P51X,0.005,

MSHAPE,0,2D

MSHKEY,0

!* Разбиение на КЭ.

CM,_Y,AREA

ASEL, , , , 3

CM,_Y1,AREA

CHKMSH,'AREA'

CMSEL,S,_Y

!*

AMESH,_Y1

!*

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

!* Закрепление по оси OY линии, каксательной к плоскости вала.

FINISH

/SOL

FLST,2,2,4,ORDE,2

FITEM,2,23

FITEM,2,48

!*

/GO

DL,P51X, ,UY,

!* Создание шарнира (закрепление по оси OY).

!FLST,2,1,1,ORDE,1

!FITEM,2,20

!*

!/GO

!D,P51X, , , , , ,UY, , , , ,

!* Задание нагрузки.

FLST,2,1,1,ORDE,1

FITEM,2,21

FLST,2,2,4,ORDE,2

FITEM,2,21

FITEM,2,47

/GO

!*

SFL,P51X,PRES,-2060000,

!* Задание температурного поля.

*get, n_count, node, 0, count

*do, i, 1, n_count

BF, i, temp, 756.811*nx(i)+416.953

*enddo

!* Задание угловой скорости и ускорения свободного падения.

OMEGA,0,-764,0,

ACEL,0,-9.81,0,

!*

FINISH

/PREP7

!* Задание параметров материала.

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,DENS,1,,8000

MPTEMP,,,,,,,,

MPTEMP,1,0

UIMP,1,REFT,,,

MPDATA,ALPX,1,,1.2e-6

!*

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,EX,1,,2.5e11

MPDATA,PRXY,1,,0.3

MSHAPE,1,3D

MSHKEY,0

FINISH

/SOL

!* Запуск на решение.

/STATUS,SOLU

SOLVE

!* Деформированное состояние.

/POST1

PLDISP,1

!* Распределение температуры.

/PSF,PRES,NORM,1,0,1

/PBF,TEMP, ,1

/PIC,DEFA, ,1

/PSYMB,CS,0

/PSYMB,NDIR,0

/PSYMB,ESYS,0

/PSYMB,LDIV,0

/PSYMB,LDIR,0

/PSYMB,ADIR,0

/PSYMB,ECON,0

/PSYMB,XNODE,0

/PSYMB,DOT,1

/PSYMB,PCONV,

/PSYMB,LAYR,0

/PSYMB,FBCS,0

!*

/PBC,ALL, ,1

/REP

FINISH

/POST1

!* Напряжения по критерию Вон-Мизеса.

/EFACET,1

PLNSOL, S,EQV, 0,1.0

!* Распределение перемещений по оси UX

PLNSOL,U,X,0,1

!* Распределение напряжений по оси UX

PLNSOL, S,X, 0,1.0

!* Распределение перемещений по оси UY

PLNSOL,U,Y,0,1

!* Распределение напряжений по оси UY

PLNSOL, S,Y, 0,1.0

!* Исследование на собственные частоты.

FINISH

/SOL

!Модальный анализ

pstres,on

solve

finish

/solu

antype,1

bucopt,subsp,7

mxpand,7,,,0

solve

finish

/POST1

SET,FIRST !Первая форма потери устойчивости.

PLDISP,1

SET,NEXT

PLDISP,1

SET,LAST

PLDISP,1

SET,LIST

FINISH

/SOL

Решим задачу в объемной постановке.

Построим объемную модель сектора диска (30°) с одним отверстием на биссектрисе под стяжной болт.

Граничные условия:

- распределенная нагрузка на внешней кромке;

- ограничение перемещений по оси OY внутренней кромки, граничащей с валом;

- условие симметрии на двух боковых поверхностях сектора.

Задаем:

- температурное поле, распределяющееся в радиальном направлении;

- угловую скорость;

- ускорение свободного падения в плоскости, перпендикулярной оси вала.

Построение модели было осуществлено в Siemens NX и импортировано для решения в ANSYS в формате Parasolid.

Расчет по циклической долговечности - student2.ru

Рис.14 Импортированная модель.

Расчет по циклической долговечности - student2.ru

Рис.15 Объемная модель с построенной сеткой.

Зададим граничные условия, угловую скорость, распределенную нагрузку, ускорение свободного падения, а также условия симметрии по боковым плоскостям (Рис.15):

Расчет по циклической долговечности - student2.ru

Рис.16 Модель с необходимыми условиями.

Расчет по циклической долговечности - student2.ru

Рис.17 Графический анализ итерационного процесса.

Расчет по циклической долговечности - student2.ru

Рис.18 Деформированное состояние.

Расчет по циклической долговечности - student2.ru

Рис.19 Температурное распределение.

Расчет по циклической долговечности - student2.ru

Рис.20 Напряжения по критерию Вон-Мизеса.

Расчет по циклической долговечности - student2.ru

Рис.21 Распределение перемещений по оси UX.

Расчет по циклической долговечности - student2.ru

Рис.22 Распределение напряжений по оси UX.

Расчет по циклической долговечности - student2.ru

Рис.23 Распределение напряжений по оси UY.

Расчет по циклической долговечности - student2.ru

Рис.24 Распределение напряжений по оси UY.

Исследование на собственные частоты.

Первая собственная форма (Рис.25):

Расчет по циклической долговечности - student2.ru

Рис.25.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.25.2 Вид спереди.

Расчет по циклической долговечности - student2.ru

Рис.25.3 Вид сбоку.

Вторая форма потери устойчивости (Рис.26):

Расчет по циклической долговечности - student2.ru

Рис.26.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.26.2 Вид спереди.

Расчет по циклической долговечности - student2.ru

Рис.26.3 Вид сбоку.

Третья форма потери устойчивости (Рис.27):

Расчет по циклической долговечности - student2.ru

Рис.27.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.27.2 Вид спереди.

Расчет по циклической долговечности - student2.ru

Рис.27.3 Вид сбоку.

Четвертая форма потери устойчивости (Рис.28):

Расчет по циклической долговечности - student2.ru

Рис.28.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.28.2 Вид спереди.

Расчет по циклической долговечности - student2.ru

Рис.28.3 Вид сбоку.

Пятая форма потери устойчивости (Рис.29):

Расчет по циклической долговечности - student2.ru

Рис.29.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.29.2 Вид спереди.

Расчет по циклической долговечности - student2.ru

Рис.29.3 Вид сбоку.

Шестая форма потери устойчивости (Рис.30):

Расчет по циклической долговечности - student2.ru

Рис.30.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.30.2 Вид спереди.

Расчет по циклической долговечности - student2.ru

Рис.30.3 Вид сбоку.

Седьмая форма потери устойчивости (Рис.31):

Расчет по циклической долговечности - student2.ru

Рис.31.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.31.2 Вид спереди.

Расчет по циклической долговечности - student2.ru

Рис.31.3 Вид сбоку.

Расчет по циклической долговечности - student2.ru

Табл.2 Таблица значений критических частот.

LOG-FILE

finish

/CLEAR,START

/units,si

/COLOR,PBAK,OFF

/REPLOT

/prep7

~PARAIN,'Sector','x_t','F:\ansys2\Models\',SOLIDS,0,0 !Импорт объемной модели.

/NOPR

/GO

!* Задаем тип КЭ.

ET,1,SOLID92

!* Задаем табличные значения.

TB,MELA,1,5,15,

TBTEMP,473

TBPT,,0,0

TBPT,,0.1/100,1880e5

TBPT,,0.2/100,3800e5

TBPT,,0.3/100,5720e5

TBPT,,0.35/100,6420e5

TBPT,,0.37/100,6550e5

TBPT,,0.4/100,6720e5

TBPT,,0.42/100,6810e5

TBPT,,0.46/100,6970e5

TBPT,,0.5/100,7100e5

TBPT,,0.55/100,7200e5

TBPT,,0.6/100,7300e5

TBPT,,0.7/100,7390e5

TBPT,,0.8/100,7460e5

TBPT,,0.86/100,7500e5

TBTEMP,573

TBPT,,0,0

TBPT,,0.1/100,1830e5

TBPT,,0.2/100,3700e5

TBPT,,0.3/100,5580e5

TBPT,,0.35/100,6310e5

TBPT,,0.37/100,6450e5

TBPT,,0.4/100,6630e5

TBPT,,0.42/100,6720e5

TBPT,,0.46/100,6880e5

TBPT,,0.5/100,7000e5

TBPT,,0.55/100,7100e5

TBPT,,0.6/100,7190e5

TBPT,,0.7/100,7260e5

TBPT,,0.8/100,7330e5

TBPT,,0.86/100,7380e5

TBTEMP,773

TBPT,,0,0

TBPT,,0.1/100,1700e5

TBPT,,0.2/100,3420e5

TBPT,,0.3/100,5100e5

TBPT,,0.35/100,5800e5

TBPT,,0.37/100,6000e5

TBPT,,0.4/100,6230e5

TBPT,,0.42/100,6330e5

TBPT,,0.46/100,6490e5

TBPT,,0.5/100,6620e5

TBPT,,0.55/100,6720e5

TBPT,,0.6/100,6760e5

TBPT,,0.7/100,6830e5

TBPT,,0.8/100,6900e5

TBPT,,0.86/100,6930e5

TBTEMP,873

TBPT,,0,0

TBPT,,0.1/100,1600e5

TBPT,,0.2/100,3200e5

TBPT,,0.3/100,4830e5

TBPT,,0.35/100,5600e5

TBPT,,0.37/100,5800e5

TBPT,,0.4/100,6000e5

TBPT,,0.42/100,6110e5

TBPT,,0.46/100,6250e5

TBPT,,0.5/100,6370e5

TBPT,,0.55/100,6470e5

TBPT,,0.6/100,6500e5

TBPT,,0.7/100,6560e5

TBPT,,0.8/100,6600e5

TBPT,,0.86/100,6630e5

TBTEMP,973

TBPT,,0,0

TBPT,,0.1/100,2470e5

TBPT,,0.2/100,3000e5

TBPT,,0.3/100,4570e5

TBPT,,0.35/100,5200e5

TBPT,,0.37/100,5130e5

TBPT,,0.4/100,5500e5

TBPT,,0.42/100,5620e5

TBPT,,0.46/100,5850e5

TBPT,,0.5/100,6000e5

TBPT,,0.55/100,6200e5

TBPT,,0.6/100,6300e5

TBPT,,0.7/100,6350e5

TBPT,,0.8/100,6400e5

TBPT,,0.86/100,6420e5

GPLOT

!*

MPTEMP,,,,,,,,

MPTEMP,1,475

MPTEMP,2,477.88

MPTEMP,3,483.17

MPTEMP,4,488.47

MPTEMP,5,499.45

MPTEMP,6,507.77

MPTEMP,7,521.39

MPTEMP,8,531.99

MPTEMP,9,552.42

MPTEMP,10,563.77

MPTEMP,11,565.29

MPTEMP,12,575.88

MPTEMP,13,577.77

MPTEMP,14,583.45

MPTEMP,15,615.77

MPTEMP,16,647.78

MPTEMP,17,680.32

MPTEMP,18,681.08

MPTEMP,19,687.13

MPTEMP,20,694.70

MPTEMP,21,700

UIMP,1,REFT,,,

MPDATA,ALPX,1,,12.92e-7

MPDATA,ALPX,1,,13.04e-7

MPDATA,ALPX,1,,13.25e-7

MPDATA,ALPX,1,,13.46e-7

MPDATA,ALPX,1,,13.9e-7

MPDATA,ALPX,1,,14.23e-7

MPDATA,ALPX,1,,14.78e-7

MPDATA,ALPX,1,,15.2e-7

MPDATA,ALPX,1,,16.02e-7

MPDATA,ALPX,1,,16.47e-7

MPDATA,ALPX,1,,16.53e-7

MPDATA,ALPX,1,,16.96e-7

MPDATA,ALPX,1,,17.03e-7

MPDATA,ALPX,1,,17.26e-7

MPDATA,ALPX,1,,18.55e-7

MPDATA,ALPX,1,,19.83e-7

MPDATA,ALPX,1,,21.13e-7

MPDATA,ALPX,1,,21.16e-7

MPDATA,ALPX,1,,21.41e-7

MPDATA,ALPX,1,,21.71e-7

MPDATA,ALPX,1,,21.92e-7

!* Задание параметров материала.

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,DENS,1,,8000

MPTEMP,,,,,,,,

MPTEMP,1,0

UIMP,1,REFT,,,

MPDATA,ALPX,1,,1.2*10e-6

!*

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,EX,1,,2.5e11

MPDATA,PRXY,1,,0.3

MSHAPE,1,3D

MSHKEY,0

!* Разбиение на КЭ.

ESIZE,0.005,0,

CM,_Y,VOLU

VSEL, , , , 1

CM,_Y1,VOLU

CHKMSH,'VOLU'

CMSEL,S,_Y

!*

VMESH,_Y1

!*

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

!* Условие закрепления по оси OY.

FLST,2,1,5,ORDE,1

FITEM,2,41

!*

/GO

DA,P51X,UY,

!* Условия симметрии.

FLST,2,2,5,ORDE,2

FITEM,2,42

FITEM,2,-43

DA,P51X,SYMM

!* Распределенная нагрузка.

FLST,2,1,5,ORDE,1

FITEM,2,22

/GO

!*

SFA,P51X,1,PRES,-2060000

!* Задание температурного поля.

*get,n_count,node,0,count

*do,i,1,n_count

Bf,i,temp,1580.46*sqrt(nx(i)*nx(i)+nz(i)*nz(i))+331.43

*enddo

!* Задание угловой скорости и ускорения свободного падения.

OMEGA,0,-764,0,

ACEL,0,-9.81,0,

!* Решение.

FINISH

/SOL

/STATUS,SOLU

SOLVE

!* Деформированное состояние.

/POST1

PLDISP,1

!* Распределение температуры.

/PSF,PRES,NORM,1,0,1

/PBF,TEMP, ,1

/PIC,DEFA, ,1

/PSYMB,CS,0

/PSYMB,NDIR,0

/PSYMB,ESYS,0

/PSYMB,LDIV,0

/PSYMB,LDIR,0

/PSYMB,ADIR,0

/PSYMB,ECON,0

/PSYMB,XNODE,0

/PSYMB,DOT,1

/PSYMB,PCONV,

/PSYMB,LAYR,0

/PSYMB,FBCS,0

!*

/PBC,ALL, ,1

/REP

FINISH

/POST1

!* Напряжения по критерию Вон-Мизеса.

/EFACET,1

PLNSOL, S,EQV, 0,1.0

!* Распределение перемещений по оси UX

PLNSOL,U,X,0,1

!* Распределение напряжений по оси UX

PLNSOL, S,X, 0,1.0

!* Распределение перемещений по оси UY

PLNSOL,U,Y,0,1

!* Распределение напряжений по оси UY

PLNSOL, S,Y, 0,1.0

!* Исследование на собственные частоты.

FINISH

/SOL

!Модальный анализ

pstres,on

solve

finish

/solu

antype,1

bucopt,subsp,7

mxpand,7,,,0

solve

finish

/POST1

SET,FIRST !Первая форма потери устойчивости.

PLDISP,1

SET,NEXT

PLDISP,1

SET,LAST

PLDISP,1

SET,LIST

FINISH

/SOL

Осуществим расчет лопатки без хвостовика.

Построим модель лопатки без хвостовика (Рис.31):

Расчет по циклической долговечности - student2.ru

Рис.32 Импортированная модель.

Расчет по циклической долговечности - student2.ru

Рис.33 Модель с наложенной сеткой и необходимыми условиями.

Расчет по циклической долговечности - student2.ru

Рис.34 Деформированное состояние.

Расчет по циклической долговечности - student2.ru

Рис.35 Температурное распределение.

Расчет по циклической долговечности - student2.ru

Рис.36 Напряжения по критерию Вон-Мизеса.

Расчет по циклической долговечности - student2.ru

Рис.37 Распределение перемещений по оси UX.

Расчет по циклической долговечности - student2.ru

Рис.38 Распределение напряжений по оси UX.

Расчет по циклической долговечности - student2.ru

Рис.39 Распределение напряжений по оси UY.

Расчет по циклической долговечности - student2.ru

Рис.40 Распределение напряжений по оси UY.

Исследование на собственные частоты.

Первая собственная форма (Рис.41):

Расчет по циклической долговечности - student2.ru

Рис.41.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.41.2 Вид спереди.

Расчет по циклической долговечности - student2.ru

Рис.41.3 Вид сбоку.

Вторая собственная форма (Рис.42):

Расчет по циклической долговечности - student2.ru

Рис.42.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.42.2 Вид спереди.

Расчет по циклической долговечности - student2.ru

Рис.42.3 Вид сбоку.

Третья собственная форма (Рис.43):

Расчет по циклической долговечности - student2.ru

Рис.43.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.43.2 Вид спереди.

+ Расчет по циклической долговечности - student2.ru

Рис.43.3 Вид сбоку.

Четвертая собственная форма (Рис.44):

Расчет по циклической долговечности - student2.ru

Рис.44.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.44.2 Вид спереди.

+ Расчет по циклической долговечности - student2.ru

Рис.44.3 Вид сбоку.

Пятая собственная форма (Рис.45):

Расчет по циклической долговечности - student2.ru

Рис.45.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.45.2 Вид спереди.

+ Расчет по циклической долговечности - student2.ru

Рис.45.3 Вид сбоку.

Шестая собственная форма (Рис.46):

Расчет по циклической долговечности - student2.ru

Рис.46.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.46.2 Вид спереди.

+ Расчет по циклической долговечности - student2.ru

Рис.46.3 Вид сбоку.

Седьмая собственная форма (Рис.47):

Расчет по циклической долговечности - student2.ru

Рис.47.1 Вид сверху.

Расчет по циклической долговечности - student2.ru

Рис.47.2 Вид спереди.

+ Расчет по циклической долговечности - student2.ru

Рис.47.3 Вид сбоку.

Расчет по циклической долговечности - student2.ru

Табл.3 Таблица значений критических частот.

LOG-FILE

finish

/CLEAR,START

/units,si

/COLOR,PBAK,OFF

/REPLOT

/prep7

~PARAIN,'Disc(sector)_net_cvyazej','x_t','F:\ansys2\Models\',SOLIDS,0,0 !* Экспорт модели.

/NOPR

/GO

FLST,2,5,6,ORDE,2

FITEM,2,2

FITEM,2,-6

VDELE,P51X

!*

ET,1,SOLID187

!* СОединяем узлы для построения совместной сетки.

nummrg,all,1e-4,1e-3

!* Разбиение на КЭ.

ESIZE,0.005,0,

MSHAPE,1,3D

MSHKEY,0

!*

CM,_Y,VOLU

VSEL, , , , 1

CM,_Y1,VOLU

CHKMSH,'VOLU'

CMSEL,S,_Y

!*

VMESH,_Y1

!*

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

!*

ESIZE,0.003,0,

FLST,5,5,6,ORDE,2

FITEM,5,7

FITEM,5,-11

CM,_Y,VOLU

VSEL, , , ,P51X

CM,_Y1,VOLU

CHKMSH,'VOLU'

CMSEL,S,_Y

!*

VMESH,_Y1

!*

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

!*

FINISH

/SOL

!* Закрепление плоскости по OY, касательной к валу.

FLST,2,1,5,ORDE,1

FITEM,2,45

!*

/GO

DA,P51X,UY,

!* Условия симметрии.

FLST,2,2,5,ORDE,2

FITEM,2,46

FITEM,2,-47

DA,P51X,SYMM

!* Задание распределенной нагрузки на лопатки.

FLST,2,5,5,ORDE,5

FITEM,2,85

FITEM,2,96

FITEM,2,107

FITEM,2,118

FITEM,2,129

/GO

!*

SFA,P51X,1,PRES,5

!* Задание параметров материала.

FINISH

/PREP7

!*

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,DENS,1,,8000

MPTEMP,,,,,,,,

MPTEMP,1,0

UIMP,1,REFT,,,

MPDATA,ALPX,1,,1.2*10e-6

!*

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,EX,1,,2.5e11

MPDATA,PRXY,1,,0.3

MSHAPE,1,3D

MSHKEY,0

!* Температурное распределение.

*get,n_count,node,0,count

*do,i,1,n_count

Bf,i,temp,1580.46*sqrt(nx(i)*nx(i)+nz(i)*nz(i))+331.43

*enddo

!* Задание угловой скорости и ускорения свободного падения.

OMEGA,0,-764,0,

ACEL,0,-9.81,0,

!* Решение.

FINISH

/SOL

/STATUS,SOLU

SOLVE

!* Деформированное состояние.

/POST1

PLDISP,1

!* Распределение температуры.

/PSF,PRES,NORM,1,0,1

/PBF,TEMP, ,1

/PIC,DEFA, ,1

/PSYMB,CS,0

/PSYMB,NDIR,0

/PSYMB,ESYS,0

/PSYMB,LDIV,0

/PSYMB,LDIR,0

/PSYMB,ADIR,0

/PSYMB,ECON,0

/PSYMB,XNODE,0

/PSYMB,DOT,1

/PSYMB,PCONV,

/PSYMB,LAYR,0

/PSYMB,FBCS,0

!*

/PBC,ALL, ,1

/REP

FINISH

/POST1

!* Напряжения по критерию Вон-Мизеса.

/EFACET,1

PLNSOL, S,EQV, 0,1.0

!* Распределение перемещений по оси UX.

PLNSOL,U,X,0,1

!* Распределение напряжений по оси UX.

PLNSOL, S,X, 0,1.0

!* Распределение перемещений по оси UY.

PLNSOL,U,Y,0,1

!* Распределение напряжений по оси UY.

PLNSOL, S,Y, 0,1.0

!* Исследование на собственные частоты.

FINISH

/SOL

!Модальный анализ

pstres,on

solve

finish

/solu

antype,1

bucopt,subsp,7

mxpand,7,,,0

solve

finish

/POST1

SET,FIRST !Первая форма потери устойчивости.

PLDISP,1

SET,NEXT

PLDISP,1

SET,LAST

PLDISP,1

SET,LIST

FINISH

/SOL

Осуществим решение для одной из лопаток с хвостовиком отдельно.

Решение

Расчет по циклической долговечности - student2.ru

Табл.4 Механические свойства сплава.

Расчет по циклической долговечности - student2.ru

Рис.48 Объемная модель лопатки с хвостовиком.

Разобьем на КЭ:

Расчет по циклической долговечности - student2.ru

Рис.49 Разбиение на КЭ.

Зададим граничные условия, распределенную нагрузку на внутреннюю часть лопатки, угловую скорость и ускорения свободного падения:

Расчет по циклической долговечности - student2.ru

Рис.50 Задание необходимых условий.

Осуществим решение. Выведем все необходимые эпюры.

Расчет по циклической долговечности - student2.ru

Рис.51 Температурное распределение.

Расчет по циклической долговечности - student2.ru

Рис.52 Деформированное и недеформированное состояние.

Расчет по циклической долговечности - student2.ru

Рис.53 Напряжения по критерию Вон-Мизеса.

Расчет по циклической долговечности - student2.ru

Рис.54 Распределение перемещений по оси OX.

Расчет по циклической долговечности - student2.ru

Рис.55 Распределение напряжений по оси OX.

Расчет по циклической долговечности - student2.ru

Рис.56 Распределение напряжений по оси OY.

Расчет по циклической долговечности - student2.ru

Рис.57 Распределение перемещений по оси OY.

Исследование на собственные частоты.

Расчет по циклической долговечности - student2.ru

Рис.58.1 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис.58.2 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 58.3 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 58.4 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 58.5 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 58.6 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 58.7 Первая форма потери устойчивости.

Расчет по циклической долговечности - student2.ru

Табл.4 Таблица значений критических частот.

LOG-FILE

finish

/CLEAR,START

/units,si

/COLOR,PBAK,OFF

/REPLOT

/prep7

~PARAIN,'Lopatka_by_kek','x_t','F:\ansys2\Models\',SOLIDS,0,0 !Экспорт объемной модели.

/NOPR

/GO

!*

ET,1,SOLID92

!*

ESIZE,0.002,0,

VCLEAR, 1

CM,_Y,VOLU

VSEL, , , , 1

CM,_Y1,VOLU

CHKMSH,'VOLU'

CMSEL,S,_Y

!*

VMESH,_Y1

!*

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

!* Задаем параметры материала.

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,DENS,1,,8400

MPTEMP,1,500

MPTEMP,2,600

MPDE,EX,1

MPDE,PRXY,1

MPDATA,EX,1,,1.7E+011

MPDATA,EX,1,,1.65E+011

MPDATA,PRXY,1,,0.3

MPDATA,PRXY,1,,0.3

MPTEMP,,,,,,,,

MPTEMP,1,500

MPTEMP,2,600

UIMP,1,REFT,,,

MPDE,ALPX,1

MPDATA,ALPX,1,,1.48E-006

MPDATA,ALPX,1,,1.58E-006

FINISH

/SOL

!* Задание температурного поля.

*get, n_count, node, 0, count

*do, i, 1, n_count

BF, i, temp,1580*ny(i)+780

*end do

!* Задание угловой скорости и ускорения свободного падения.

OMEGA,0,0,764,

ACEL,0,0,9.81,

!*

FINISH

/PREP7

!*

FLST,2,2,5,ORDE,2

FITEM,2,1

FITEM,2,5

!*

/GO

DA,P51X,ALL,

FLST,2,1,5,ORDE,1

FITEM,2,9

/GO

!*

SFA,P51X,1,PRES,5

!* Запуск на решение.

FINISH

/SOL

/STATUS,SOLU

SOLVE

!* Деформированное состояние.

/POST1

PLDISP,1

!* Распределение температуры.

/PSF,PRES,NORM,1,0,1

/PBF,TEMP, ,1

/PIC,DEFA, ,1

/PSYMB,CS,0

/PSYMB,NDIR,0

/PSYMB,ESYS,0

/PSYMB,LDIV,0

/PSYMB,LDIR,0

/PSYMB,ADIR,0

/PSYMB,ECON,0

/PSYMB,XNODE,0

/PSYMB,DOT,1

/PSYMB,PCONV,

/PSYMB,LAYR,0

/PSYMB,FBCS,0

!*

/PBC,ALL, ,1

/REP

FINISH

/POST1

!* Напряжения по критерию Вон-Мизеса.

/EFACET,1

PLNSOL, S,EQV, 0,1.0

!* Распределение перемещений по оси UX

PLNSOL,U,X,0,1

!* Распределение напряжений по оси UX

PLNSOL, S,X, 0,1.0

!* Распределение перемещений по оси UY

PLNSOL,U,Y,0,1

!* Распределение напряжений по оси UY

PLNSOL, S,Y, 0,1.0

!* Исследование на собственные частоты.

FINISH

/SOL

!Модальный анализ

pstres,on

solve

finish

/solu

antype,1

bucopt,subsp,7

mxpand,7,,,0

solve

finish

/POST1

SET,FIRST !Первая форма потери устойчивости.

PLDISP,1

SET,NEXT

PLDISP,1

SET,LAST

PLDISP,1

SET,LIST

FINISH

/SOL

Расчет по циклической долговечности

Зоны диска r, м T, K σr, МПа σt, МПа σнэкв, Мпа αT σ, Мпа σ=αT∙σнэкв
Центральное отверстие 7,67 475,00 0,00 8600,02 8600,02 1,68 14448,03
Отверстие в полотне (точка D1) 15,70 542,20 3534,42 5079,71 4516,12 3,05 13774,18
Дно паза замка 37,40 700,00 2649,81 -531,64 2951,76 7,98 23567,36

Табл.5 Исходные данные.

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru , в силу симметричности цикла.

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru - число циклов до образования трещины,

Расчет по циклической долговечности - student2.ru - число циклов нагружения,

Расчет по циклической долговечности - student2.ru - ресурс,

Расчет по циклической долговечности - student2.ru - средняя продолжительность работы,

Расчет по циклической долговечности - student2.ru - предел длительной прочности,

Расчет по циклической долговечности - student2.ru - среднее напряжение,

Расчет по циклической долговечности - student2.ru - размах упруго - пластической деформации.

Решим уравнение Мэнсона в программе Mathematica (см. рис. 98).

Расчет по циклической долговечности - student2.ru

Рис.59 Решение уравнения Мэнсона.

Расчет по циклической долговечности - student2.ru

Расчет по циклической долговечности - student2.ru , Расчет по циклической долговечности - student2.ru .

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

Расчет по циклической долговечности - student2.ru

Рис.60.1 Эпюра перемещений вдоль оси OY для первой формы потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 60.2 Эпюра перемещений вдоль оси OY для второй формы потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 60.3 Эпюра перемещений вдоль оси OY для третьей формы потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 60.4 Эпюра перемещений вдоль оси OY для четвертой формы потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 60.5 Эпюра перемещений вдоль оси OY для пятой формы потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 60.6 Эпюра перемещений вдоль оси OY для шестой формы потери устойчивости.

Расчет по циклической долговечности - student2.ru

Рис. 60.7 Эпюра перемещений вдоль оси OY для седьмой формы потери устойчивости.

Наши рекомендации