Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Manual.pdf
Скачиваний:
32
Добавлен:
18.03.2015
Размер:
4 Мб
Скачать
☆

Задачи теории упругости

251

F = −λ

 

 

∂T

 

,

F = −λ

 

∂T

в плоском случае;

x ∂x

 

y ∂y

x

 

 

y

 

F = −λ

 

∂T

,

F = −λ

 

∂T

в осесимметричном случае;

z ∂z

r ∂r

z

 

 

r

 

Интегральные величины:

•Поток тепла через заданную поверхность

Φ= ∫F nds ,

где n - единичный вектор нормали к поверхности. Поверхность интегрирования задается контуром в плоскости модели, состоящим из отрезков и дуг окружностей.

Задачитеорииупругости

Комплекс ELCUT может решать задачи теории упругости в постановках плоских напряжений, плоских деформаций и осесимметричного напряженного состояния с изотропными или ортотропными свойствами материалов. Задача плоских напряжений подходит для анализа структур, тонких по глубине, которые нагружены в плоскости модели. Напряжение в направлении, нормальном к плоскости модели, предполагается отсутствующим. Задача плоских деформаций предполагает отсутствие деформаций вне плоскости модели. Эта задача подходит для моделирования объектов с весьма большой толщиной в направлении, нормальном к плоскости модели.

Перемещения, напряжения, деформации

Во всех постановках поле перемещений однозначно определяется двумя компонентами вектора перемещений δ в каждой точке:

δ

 

- в плоских задачах;

{δ} = x

 

δy

 

δz

 

- в осесимметричных задачах.

{δ} =

 

δr

 

 

252Глава 10 Теоретическое описание

Вобеих плоских постановках рассматривается только по три компоненты деформаций и напряжений. Деформация связана с перемещением соотношением:

 

 

 

 

 

∂δx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂x

 

 

 

ε

 

 

 

 

 

 

 

 

 

x

 

 

∂δy

 

 

 

{ε} =

εy

=

 

 

 

 

 

.

 

 

∂y

 

 

γ

 

 

∂δx

 

 

 

 

 

 

xy

 

 

+

∂δy

 

 

 

 

 

 

 

 

 

 

 

 

∂y

 

 

∂x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соответствующее ей напряжение выражается как

σx {σ} = σy .

τxy

В осесимметричном случае радиальное перемещение приводит к деформации εθ, в направлении, перпендикулярном плоскости модели, поэтому выражение для полной деформации имеет вид:

 

 

 

 

 

 

∂δz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂z

 

 

 

εz

 

 

 

 

 

 

 

 

 

 

 

 

 

∂δr

 

 

 

 

 

εr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂r

 

 

 

{ε} =

εθ

 

=

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

δr

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

γ

 

 

 

 

 

 

 

 

 

 

rz

 

 

∂δ

 

 

 

 

 

∂δ

 

 

 

 

 

z

 

 

+

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂r

 

 

 

 

 

∂z

 

 

 

 

 

 

 

 

 

 

 

Соответствующие компоненты напряжений:

σzσ {σ} = σr .θτrz

Задачи теории упругости

253

Уравнения статического равновесия для плоских задач имеют следующий вид:

 

∂σ

x

+

∂τ

xy

= − fx

 

 

 

 

 

 

∂y

 

∂x

 

 

,

 

∂τxy

 

 

∂σy

 

+

 

= − fy

 

∂x

 

∂y

 

 

 

 

а для осесимметричных задач:

1 ∂

(rσ

r

)

+

∂τ

rz = − fr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂r

 

 

 

∂z

 

,

r

 

 

 

∂σz

1 ∂

(rτrz )

+

= − fz

 

 

 

 

∂r

 

 

∂z

r

 

 

 

 

 

где ƒx, ƒy и ƒz, ƒr - компоненты вектора плотности объемной силы.

Соотношение между напряжениями и деформациями при упругом поведении материалов выражается зависимостью:

{σ} = [D]({ε}− {ε0 }),

где [D] - матрица упругости, и {ε0} - начальная термическая деформация, вызванная перепадом температур. Вид матрицы зависит от применяемой постановки.

Матрица упругости для плоского напряженного состояния (изотропный материал):

 

 

 

 

 

ν

0

 

 

 

 

E

1

 

 

[D] =

 

 

ν 1

0

 

;

 

− ν2

1

 

 

 

1 − ν

 

 

 

 

0

0

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

254 Глава 10 Теоретическое описание

плоские напряжения, ортотропный материал:

 

1

 

−

 

νyx

0

−1

 

 

 

 

 

 

 

 

 

 

 

 

Ex

 

Ey

 

 

 

 

 

 

 

 

 

 

 

νyx

 

 

1

 

 

 

 

 

[D]= −

 

 

 

 

 

 

 

 

0

 

;

 

Ey

 

 

Ey

 

 

 

 

1

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Gxy

 

плоские деформации в изотропном материале:

 

 

 

 

1

 

 

ν

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 − ν

 

 

 

E(1 − ν)

 

ν

 

 

 

 

 

[D] =

 

 

 

1

0

 

;

(1

 

 

ν

 

+ ν)(1 − 2ν) 1 −

 

 

 

1 − 2ν

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2(1 − ν)

 

 

 

 

 

 

 

 

 

 

 

 

 

плоские деформации в ортотропном материале:

 

1

 

−

νzx2

−

νyx

 

−

νzxνzy

0

−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ex

Ez

 

Ey

 

Ez

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

νyx

 

νzxνzy

 

1

 

 

 

νzy2

 

 

 

 

 

[D]= −

 

 

 

−

 

 

 

 

 

−

 

 

 

0

 

;

 

Ey

Ez

 

 

Ey

 

 

Ez

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xy

 

осесимметричная задача, изотропный материал:

 

 

 

 

 

1

 

 

ν

 

ν

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 − ν

1 −ν

ν

 

 

 

 

 

 

 

ν

 

 

 

 

 

 

 

 

E(1 − ν)

 

 

 

 

 

1

 

 

 

0

 

 

 

 

 

 

 

−

ν

1 −

ν

 

 

[D] =

 

 

1

 

 

 

 

 

 

;

 

 

 

 

ν

 

 

ν

 

 

 

 

 

 

 

 

(1

+ ν)(1 − 2

ν)

 

 

 

 

 

1

 

0

 

 

 

 

 

− ν

1 − ν

 

 

 

 

 

 

 

1

 

 

 

 

1 − 2

ν

 

 

 

 

 

 

 

0

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

2(1 −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ν)

 

Eθ
0
Er
νθz
νrz

Задачи теории упругости

255

осесимметричная задача, ортотропный материал:

1

Ez−

[D]=

−

−

 

ν

rz

 

−

 

νθz

0

−1

 

 

 

 

 

 

 

 

Er

 

Eθ

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

−

 

νθr

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Er

 

Eθ

 

 

 

 

 

 

 

 

 

.

 

 

 

νθr

 

 

 

1

 

 

 

 

−

 

 

 

 

 

0

 

 

Eθ

 

 

 

Eθ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zr

В этих формулах E - модуль Юнга в изотропном случае; Ex, Ey, Ez, Er, и Eθ - модули Юнга в направлении соответствующих осей в материале с ортотропными свойствами; ν - коэффициент Пуассона изотропного материала; νyx, νzx, νzy, νrz, νθz, νθr - коэффициенты Пуассона для указанных пар осей в ортотропном случае; Gxy или Gzr - модуль сдвига.

Температурные деформации

Температурная (начальная) деформация материала определяется коэффициентами линейного расширения и изменением температуры относительно температуры недеформированного состояния. Составляющие начальной деформации для плоского напряженного состояния и изотропного материала определяются соотношением:

α

{ε0 } = α T ;0

плоское напряженное состояние, ортотропный материал:

αx {ε0 } = αy T ;

0

плоские деформации, изотропный материал:

α {ε0 } = (1 + ν) α T ;

0

256 Глава 10 Теоретическое описание

плоские деформации, ортотропный материал:

 

αx

+νzx

αz

 

{ε0

 

+νzy

 

T ;

}= αy

αz

 

 

0

 

 

 

 

 

 

осесимметричная задача для изотропного материала:

αα

{ε0 } = α T ;0

осесимметричная задача для ортотропного материала:

αzα

{ε0 } = αr T ,θ0

где α - коэффициент линейного расширения в изотропном материале; αx, αy, αz, αr, αθ - коэффициенты линейного расширения в направлении соответствующих осей в ортотропном материале; T - перепад температуры между деформированным и недеформированным состоянием.

Внешние силы

ELCUT позволяет учитывать сосредоточенные, поверхностные и объемные силы, действующие на модель. Сосредоточенные нагрузки задаются в вершинах в форме двух компонент вектора. Поверхностные силы задаются на сторонах модели либо компонентами вектора, либо величиной нормального давления. Объемные силы задаются вектором в пределах блоков. Каждая компонента объемной силы в пределах подобласти может быть задана функцией координат. Это применяется, например, для учета центробежных сил в осесимметричных задачах. Давление также может быть задано функцией, зависящей от координат, например, для моделирования гидростатического давления.

Задачи теории упругости

257

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

Любая поверхностная сила, действующая по нормали к стороне в каждой ее точке, может быть описана как "давление". При этом на сторонах, составляющих внешнюю границу области, положительным считается давление, направленное извне вовнутрь области, а на внутренних сторонах - вправо при движении от начала к концу стороны. (В момент создания сторона получает естественную ориентацию, определяемую порядком указания вершин, являющихся ее концами. Для дуг положительным всегда считается направление против часовой стрелки).

Условия закрепления

Жесткое закрепление по любой из координат, или по обеим сразу, может быть задано в любой вершине или вдоль любой стороны области. Фиксированное перемещение вдоль линий закрепления описывается линейной функцией координат.

Условие пружинного подвеса описывает вершину, подвергающуюся действию силы, пропорциональной разнице между фактическим и начальным перемещением этой вершины. Оно характеризуется жесткостью пружины и начальным перемещением ее закрепленного конца.

Замечание. Для того чтобы задача была сформулирована корректно, необходимо задать жесткое или пружинное закрепление таким образом, чтобы исключить любое перемещение модели и ее частей, не изменяющее энергии состояния. В плоских задачах такое закрепление должно ограничивать сдвиг всей модели по осям x и y, а также поворот; в осесимметричных задачах достаточно ограничить сдвиг параллельно оси z.

258 Глава 10 Теоретическое описание

Вычисляемые физические величины

При анализе результатов решения задач теории упругости, ELCUT позволяет оперировать со следующими локальными и интегральными физическими величинами:

Локальные величины:

•Модуль вектора перемещения

δ= δ2x + δ2y , или δ = δ2z + δ2r ;

•Максимальное и минимальное главные напряжения в плоскости модели σ1 и

σ2;

•Нормальные и касательные напряжения по координатным осям σx, σy и τxy (σz, σr и τrz в осесимметричном случае);

•Нормальное напряжение в направлении, перпендикулярном к плоскости

модели (σz - для плоских задач, σθ - в осесимметричном случае). В задачах плоско-напряженного состояния эта величина равна нулю по определению;

•Критерий прочности фон Мизеса (потенциальная энергия формоизменения):

σ

 

=

1

(σ −σ

 

)2

+(σ

 

−σ

 

)2

+(σ

 

−σ )2

 

;

 

 

 

 

 

 

 

e

 

2

1

2

 

 

2

 

3

 

 

3

1

 

 

где σ1, σ2 и σ3 - главные напряжения, упорядоченные по убыванию.

•Критерий прочности Треска (максимальное сдвиговое напряжение):

σe = σ1 − σ3;

•Критерий прочности Мора-Кулона (модификация критерия Треска, учитывающая разное сопротивление материала растяжению и сжатию):

σe = σ1 − χσ3,

где

χ= [σ+ ] , [σ− ]

[σ+] и [σ−] - предельно допустимые напряжения растяжения и сжатия соответственно.

• Критерий прочности Друкера-Прагера:

Задачи теории упругости

259

σ

e

=

(

 

 

 

)

 

i

 

 

 

 

χ −χ

 

 

1

 

1

−

χ

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ χ

 

σ

 

−

 

 

 

 

σ+

 

 

 

 

 

 

 

 

σ ,

 

 

 

 

 

 

 

 

 

 

 

 

1

+ χ

 

 

[σ− ] 1

+ χ

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

=

 

1

(σ −σ

 

)2

+(σ

 

−σ

 

)2

+(σ

 

−σ )2

;

 

 

 

 

 

 

 

 

i

 

 

2

1

 

 

 

2

 

 

 

2

 

3

 

 

 

 

 

3

 

1

 

 

σ = σ1 +σ2 +σ3 . 3

• Коэффициент запаса прочности по Хиллу для ортотропных материалов:

 

σ2

 

σ σ

 

σ2

 

τ2

 

F.I. =

1

−

1 2

+

2

+

12

,

X12

X 22

S122

 

 

X12

 

 

 

где σ1, σ2 и τ12 - вычисляемые напряжения в материале и,

X1 = X1T если σ1>0; X1 = X1C если σ1<0;

X2 = X2T если σ2>0; X2 = X2С если σ2<0;

S12 = S12+ если τ12>0; S12 = S12- если τ12<0,

где X1T, X2T, X1C, X2C, S12+ и S12- - предельно допустимые напряжения растяжения, сжатия и сдвига.

Интегральные величины:

• Полная сила, действующая на заданный объем

F = С∫ (σn)ds ,

где σ - тензор напряжений.

Интеграл вычисляется по границе указанного объема, и n означает единичный вектор внешней нормали.

•Полный вращающий момент относительно начала координат, действующий на заданный объем

T = 12 С∫ (r ×(σn))ds ,

где r – радиус-вектор точки интегрирования.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]