Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
T-FLEX CAD - Руководство пользователя / T-FLEX Анализ. Пособие по работе с системой.pdf
Скачиваний:
243
Добавлен:
16.05.2015
Размер:
7.83 Mб
Скачать

Верификационные примеры

Эквивалентные напряжения определяются по формуле:

σэкв(ρ) = (σρ (ρ))2 + (σt (ρ))2 2 σρ (ρ) σt (ρ) .

Эквивалентные напряжения на внутренней поверхности сферы равны σэкв(r) =1148 МПа. Выполнив расчет при помощи T-FLEX Анализ, получаем следующие результаты:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1.

 

 

 

 

 

 

 

 

 

 

 

Параметры конечно-элементной сетки

 

 

 

 

 

 

 

Число конечных элементов

 

 

Тип конечных элементов

 

Число узлов

 

 

 

 

 

 

 

 

 

 

5338

 

 

 

 

 

 

 

 

 

квадратичный тетраэдр (10 узлов)

 

1846

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.

 

 

 

 

 

 

 

 

 

 

 

 

Результат «Перемещение»

 

 

 

 

 

 

 

 

 

 

w w*

 

 

 

 

 

 

 

Численное решение

 

Аналитическое решение

 

Ошибка δ =

 

 

 

 

 

 

 

×100%

 

 

Перемещение w

*

, м

 

 

Перемещение w , м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.03

 

 

 

 

 

 

 

 

 

1.4067E-003

 

 

 

 

 

1.4063E-003

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3.

 

 

 

 

 

 

 

 

 

 

Результат «Эквивалентные напряжения»

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Численное решение

 

Аналитическое решение

 

 

 

 

 

σ σ*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Напряжение σ

*

,

 

Н

 

 

Напряжение σ ,

Н

 

 

Ошибка δ =

 

 

 

 

 

 

 

 

 

×100%

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м2

 

м2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.17

 

 

 

 

 

 

 

 

 

1.1499E+009

 

 

 

 

 

1.1480E+009

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Рассматривается жестко закреплённая квадратная пластина под сосредоточенной силой в центре [4].

Задача статического анализа решается с использованием пластинчатых конечных элементов. В расчёте рассматривается 14 часть пластины с

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

183

Руководство пользователя T-FLEX Анализ

Конечно-элементная модель пластины с нагрузками и закреплениями

Решение этой задачи для прогиба под силой даёт формула: w = 0.0224

Pa

2

Eh3

 

 

, D =

 

.

D

12(1ν 2 )

 

 

 

Примем следующие исходные данные: длина и ширина пластины a = 500мм, толщина пластины h = 3мм, приложенная сосредоточенная сила P = 50кгс (или P = 490.3325Н ).

Характеристики материала примем E = 2.1×1011 Па, ν = 0.28. Аналитическое решение имеет вид: w = 5.3557×10-3 м.

Выполнив расчет при помощи T-FLEX Анализ, получили следующие результаты:

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1.

 

 

 

 

Параметры конечно-элементных сеток

 

 

 

 

 

 

 

 

 

 

 

 

 

Номер

Тип конечных элементов

 

Число узлов

 

Число конечных

сетки

 

 

элементов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

квадратичный тетраэдр (10 узлов)

 

8083

 

 

24222

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

линейный треугольник (6 узлов)

 

3818

 

 

7426

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

квадратичный треугольник(6 узлов)

 

3818

 

 

7426

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.

 

 

 

 

 

 

Результат «Перемещение»

 

 

 

 

 

 

 

 

 

Номер

Численное решение

Аналитическое решение

Ошибка δ =

 

w w*

 

×100%

сетки

Перемещение w* , м

Перемещение w , м

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

5.3700E-003

 

5.3557E-003

 

 

0.27

 

 

 

 

2

5.4078E-003

 

5.3557E-003

 

 

0.97

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

5.3763E-003

 

5.3557E-003

 

 

0.38

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Цилиндрический резервуар со стенками постоянной толщины

Резервуар подвергается действию давления жидкости, как показано на рисунке. Днище резервуара защемлёно в абсолютно жёсткий фундамент.

184

Верификационные примеры

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

Конечно-элементная модель конструкции с нагрузками и закреплениями

Аналитическое решение этой задачи имеет вид [5]:

w = eβ zˆ (C1 cos(βzˆ)+C2 sin(βzˆ))ρg(d zˆ)R2 , Eh

185

Руководство пользователя T-FLEX Анализ

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3(1ν 2 )

 

 

 

 

 

 

 

 

 

 

z = z + d ,

C1 =

 

ρgR2d

,

 

ρgR2

 

1

 

,

β

4

=

,

 

 

 

 

 

 

 

 

 

 

C2 =

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

Eh

Eh

 

β

 

 

2

h

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

ρ - плотность жидкости,

g - ускорение свободного падения ( 9.8

м

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с2

 

 

 

 

 

 

 

 

 

 

Примем

следующие исходные

данные:

глубина

резервуара

d =1000мм, радиус

R = 200мм ,

толщина резервуара h = 3мм, плотность жидкости ρ =1000 кг

м

3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Характеристики материала примем: E = 2.1×1011 Па, ν = 0.28.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, w = 6.1366×10-7

 

м (при zˆ = 0.056 м).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выполнив расчет при помощи T-FLEX Анализ, получили следующие результаты:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Параметры конечно-элементных сеток

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Номер

 

 

Тип конечных элементов

 

 

 

Число узлов

 

 

 

Число конечных

сетки

 

 

 

 

 

 

 

 

элементов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

квадратичный тетраэдр (10 узлов)

 

 

 

 

23092

 

 

 

 

 

 

 

70280

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

линейный треугольник (6 узлов)

 

 

 

 

 

11676

 

 

 

 

 

 

 

23140

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

квадратичный треугольник(6 узлов)

 

 

 

 

11676

 

 

 

 

 

 

 

23140

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Результат «Перемещение»

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Номер

 

 

Численное решение

 

Аналитическое решение

 

 

Ошибка δ =

 

w w*

 

×100%

сетки

 

 

Перемещение w* , м

 

Перемещение w , м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

6.2345E-007

 

 

 

 

 

6.1366E-007

 

 

 

 

 

 

 

 

 

 

 

1.60

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

6.1938E-007

 

 

 

 

 

6.1366E-007

 

 

 

 

 

 

 

 

 

 

 

0.93

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

6.1743E-007

 

 

 

 

 

6.1366E-007

 

 

 

 

 

 

 

 

 

 

 

0.61

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кручение бруса с круглым поперечным сечением

Рассмотрим брус с круглым поперечным сечением радиуса R . Длина бруса l (см. рис.).

186

Верификационные примеры

Выберем систему координат с осью z вдоль оси бруса, при этом z = 0 на его левом конце.

Брус находится под действием внешнего скручивающего момента Mкр . Момент приложен на правом конце бруса, левый конец бруса жёстко закреплён.

Примем следующие исходные данные: длина l бруса равна 0.6м, радиус поперечного сечения R бруса равен 0.02м, величина приложенного крутящего момента M кр =100 Н м .

Характеристики материала: E = 2.1×1011 Па, ν = 0.28.

Для определения угла закручивания воспользуемся соотношением [6]:

 

 

z

M

êð

 

 

 

 

 

 

 

ϕ =

0

 

dz +ϕ0

,

 

 

 

 

 

 

 

 

 

 

 

GJ p

 

 

 

 

 

где ϕ0

– угол поворота сечения z = 0 , G =

 

E

 

 

- модуль сдвига,

J p =

πR4

- полярный момент

2(1+ν)

2

 

 

 

 

 

 

инерции круглого сечения.

Так как, по условию, левый конец бруса жёстко закреплён, то ϕ0 = 0 . Тогда на расстоянии z = 0.5l от жёстко закреплённого конца бруса угол закручивания ϕ определяется формулой:

 

 

ϕ0.5l =

0.5Mкрl

 

 

 

 

 

 

GJ p

 

 

 

 

 

 

 

 

 

 

 

Таким образом, ϕ0.5l =1.4551×103 рад.

 

 

 

 

 

 

Выполнив расчет при помощи T-FLEX Анализ, получаем следующие результаты.

Таблица 1.

 

 

 

 

 

 

 

 

 

 

 

 

Параметры конечно-элементной сетки

 

 

 

 

Число узлов при

 

 

Тип конечных элементов

 

Число главных

 

 

Число конечных

 

узлов

 

расчёте задачи

 

элементов

 

 

 

 

 

 

 

 

30188

 

 

квадратичный тетраэдр (10 узлов)

 

4526

 

 

18345

187