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

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

Модуль перемещения (при z = 0.5l ) u = 2.9139E-005 м.

Таблица 2. Результат «Угол закручивания»

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

 

 

 

ϕ

ψ

 

 

 

 

Угол закручивания

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

Ошибка δ =

 

 

 

 

×100%

 

 

 

u

 

Угол закручивания ϕ , рад

 

 

 

 

 

 

 

рад

 

 

ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ = arcsin

,

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.4570E-003

 

1.4551E-003

1.2608E-001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Стержень под действием собственного веса

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

Примем следующие исходные данные: длина стержня l равна 1м, радиус поперечного сечения стержняR равен 0.02м.

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

Полное удлинение стержня под действием собственного веса можно определить по формуле [7]:

l = γ l 2 ,

2E

где γ – удельный вес материала стержня, то есть γ = ρ g , g = 9.80665 см2 .

Напряжение в сечении стержня, находящемся на расстоянии x от нижнего (незакрёпленного) конца, можно вычислить по формуле [7]:

σ = γ x .

Таким образом, l =1.8212×107 м; σ = 3.8246×104

Н

 

при

x =

l

.

 

 

 

 

 

 

 

м2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

Таблица 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

 

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

 

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

узлов

 

 

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

 

 

элементов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

2482

 

 

 

 

 

16089

 

 

9518

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

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

 

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

 

Ошибка δ =

w w*

 

×100%

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

*

, м

 

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

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.8177E-007

 

 

 

1.8212E-007

 

 

 

 

 

1.9410E-001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

188

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

Таблица 3. Результат «Напряжение»

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

 

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

 

 

 

σ σ*

 

 

 

 

 

 

 

 

Напряжение σ

*

,

Н

 

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

Н

 

Ошибка δ =

 

 

 

 

 

 

 

×100%

 

 

 

 

 

 

 

 

 

 

 

σ

 

 

 

 

 

 

 

 

 

 

м2

 

м2

 

 

 

 

 

 

 

 

 

 

 

8.0139E-003

 

3.8249E+004

 

 

3.8246E+004

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчёт вращающегося сплошного диска постоянной толщины

Рассмотрим диск радиуса R и толщиной h (см. рис.).

В расчёте рассматривается 14 часть диска с наложением условий симметрии на соответствующих

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

Диск находится под действием центробежной силы F = ρ ω2 R , где ρ - масса единицы объёма материала диска, ω - угловая скорость вращения.b

Примем следующие исходные данные: радиус диска R равен 0.457м , толщина диска h равна

0.01м, величина угловой скорости вращения ω равна 300 радс .

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

189

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

Перемещения u в таком случае можно определить по формуле [8]:

 

 

1

 

 

 

1 (

1ν

2

)

2

3

 

 

 

 

 

 

(1ν )C1r (1ν )C2 r

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

u = E

8

 

 

ρω r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где постоянные C =

3 +ν

ρω2 R2

, C

2

= 0 определяются из граничных условий.

 

1

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Максимальное перемещение umax

ожидается при r = R , то есть umax =

(1ν )ρω2 R3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

4E

Компоненты напряжений σr , σθ находятся в виде [8]:

σr = 3 +8ν ρω2 (R2 r 2 ),

σθ = 3 +8ν ρω2 R2 1+83ν ρω2r2 .

Эти напряжения принимают максимальное значение в центре диска, где

 

 

 

 

 

 

 

σr

 

=σθ

=

3 +ν

ρω2 R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, umax = 5.7430×105 м, σ

 

 

= σθ

=σ = 6.0111×107

 

Н

.

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м2

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

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

 

 

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

 

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

 

 

 

узлов

 

 

 

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

 

 

 

 

 

элементов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

1553

 

 

 

 

9010

 

 

 

 

4372

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

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

 

 

Ошибка δ =

 

 

 

w w*

 

 

×100%

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

*

, м

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.8033E-005

 

 

 

 

 

 

 

5.7430E-005

 

 

 

 

 

 

1.0390E+000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Результат «Напряжение»

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

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

 

 

 

 

 

 

 

 

σ σ*

 

 

 

 

 

 

 

 

 

 

 

 

 

Напряжение σ

*

,

 

Н

 

 

 

 

 

Напряжение

σ ,

Н

 

 

 

Ошибка δ =

 

 

 

 

 

 

 

 

×100%

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м2

 

 

 

 

м2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.4740E+007

 

 

 

 

 

 

 

6.0111E+007

 

 

 

 

 

 

7.7009E+00

 

190

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

Свободно опертая прямоугольная пластинка под синусоидальной нагрузкой

Рассмотрим прямоугольную пластинку со сторонами a , b и толщиной h (см. рис.).

Толщина пластины h существенно меньше длин её сторон a , b .

Пластина находится под действием нагрузки, распределённой по поверхности пластины согласно

закону: q = q

sin

 

π x sin

 

π y

, где q0 представляет собой интенсивность нагрузки в центре

0

 

 

a

 

 

b

 

 

 

 

 

 

 

 

 

пластинки.

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

Примем следующие исходные данные: длина стороны a пластины равна 0.5м, длина стороны b пластины равна 0.4м, толщина пластины h = 0.003м , интенсивность нагрузки в центре пластинки

q0 =100 мН2 .

Характеристики материала: E = 2.1×1011 Па, ν = 0.28. Аналитическое решение этой задачи имеет вид [5]:

 

 

 

w =

 

 

 

 

q0

 

 

 

 

 

 

 

 

 

π x

 

 

 

π y

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

1

 

1

 

2

a

 

 

 

 

 

 

 

 

π

4

D

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

b

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где D =

Eh3

 

– цилиндрическая жесткость пластины.

 

 

 

 

 

 

 

 

 

 

12(1ν 2 )

 

 

 

 

 

 

 

 

 

 

Максимальный прогиб получается в центре пластинки.

 

Подставив

x =

a

и y =

b

в уравнение,

 

 

 

получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wmax

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

1

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

D

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

b

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

w =1.9059×105

м.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Врасчёте рассматривается 14 часть пластины с наложением условий симметрии на

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

191

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

 

 

 

Модель пластины с нагрузками и закреплениями

 

 

 

 

 

 

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

 

 

 

Таблица 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Номер

 

 

 

 

 

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

 

 

 

Число

 

 

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

 

 

Число узлов

 

 

 

конечных

 

 

сетки

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

элементов

 

 

1

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

 

5836

 

 

34698

 

 

 

 

17318

 

 

2

 

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

 

2747

 

 

2747

 

 

 

 

5304

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

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

2747

 

 

10797

 

 

 

 

5304

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

Номер

 

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

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

 

 

 

w w*

 

 

 

сетки

 

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

 

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

 

Ошибка δ =

 

 

 

 

 

×100%

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1.8741E-005

 

 

1.9059E-005

 

 

1.6660E+000

 

 

 

 

 

 

 

 

 

 

 

2

 

1.8798E-005

 

 

1.9059E-005

 

 

1.3693E+000

 

 

 

 

 

 

 

 

 

 

 

3

 

1.8752E-005

 

 

1.9059E-005

 

 

1.6098E+000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

192