CALCULATION OF STRESS-DEFORMED STATE OF ELASTIC CYLINDERS UNDER TEMPERATURE AND HUMIDITY IMPACT
Abstract and keywords
Abstract (English):
The actual scientific and practical problem of plane deformation of an elastic cylinder under conditions of temperature and humidity is considered. Various pipes, shafts, plain bearings, bushings from natural and modified wood, etc. have a cylindrical shape. An exact analytical solution of this problem is given for an isotropic cylinder in the case of stationary temperature and humidity exposure. All formulas for stresses and strains and stresses contain mechanical and thermophysical constants, which corresponds to the physical meaning. It is proved that the use of classical methods for solving this problem through the potentials of stresses and displacements leads to an undesirable result in which stresses or deformations do not depend on the thermophysical constant C_2, which contradicts the physical meaning.

Keywords:
deformation of an elastic cylinder, temperature and humidity effect on an isotropic cylinder, Hooke's law
Text
Publication text (PDF): Read Download

 

ВВЕДЕНИЕ

Проектирование и расчет напряженно-деформированного состояния упругих и упруго-вязко-пластических цилиндров из сжимаемого материала имеет важное научно-практическое значение [1, 10]. Цилиндр является основной деталью при изготовлении труб, валов, подшипников скольжения, различных втулок из естественной и модифицированной древесины и т.д. Однако до настоящего времени использующиеся известные методы решения таких задач нередко приводят к нежелательному результату – независимости полученных формул от механических и теплофизических констант. Поэтому получение общих решений, которые позволяют изучить зависимость напряжений и деформаций от механических и теплофизических постоянных имеет большое значение. Такие формулы могут успешно быть внедрены как в научно-исследовательских отчетах, так и в учебном процессе в курсах сопротивления материалов и древесиноведения.

 

1 ПЛОСКАЯ ДЕФОРМАЦИЯ УПРУГОГО ЦИЛИНДРА

ПРИ ТЕМПЕРАТУРНО-ВЛАЖНОСТНОМ ВОЗДЕЙСТВИИ

Рассмотрим плоскую деформацию упругого изотропного цилиндра при полярно-симметричном деформировании в условиях температурно-влажностного воздействия. Основная система уравнений имеет следующий вид [1, 2, 5, 8].

Уравнение равновесия:

rdσrdr+σr-σθ=0.                                                (1)

Геометрические соотношения Коши:

εr=dudr;     εθ=ur.                                                 (2)

Уравнение совместности деформаций:

ddrr2dεθdr-rdεrdr=0.                                            (3)

Обобщенный закон Гука в случае стационарного теплового воздействия:

εr=1Eσr-μσθ+σz+αTr;                                   (4)

εθ=1Eσθ-μσr+σz+αTr;                                   (5)

εz=1Eσz-μσθ+σr+αTr.                                   (6)

В случае плоской деформации:

εz=0;                                                              (7)

σz=μσθ+σr-EαTr.                                         (8)

Подставим формулу (8) в закон Гука (4)-(5):

εr=1Eσr-μσθ-μμσθ+σr-EαT+αT==1E1-μ2σr-μ1+μσθ+α1+μT;                  (9)

εθ=1E1-μ2σθ-μ1+μσr+α1+μT.                   (10)

Уравнение теплопроводности в случае полярной симметрии имеет вид [1, 4-8]:

ddrrdTdr=0.                                                 (11)

Решение этого уравнения имеет вид:

Tr=C1lnr+C2.                                            (12)

Граничные условия имеют вид:

При r=r1: T=T1; при r=r2: T=T2.                                                      (13)

С учетом граничных условий, получим:

C1=T2-T1lnr2r1;     C2=T1lnr2-T2lnr1lnr2r1.                   (14)

Обобщенный закон Гука (9)-(10):

εr=1E1-μ2σr-μ1+μσθ+α1+μC1lnr+α1+μC2;   (15)

εθ=1E1-μ2σθ-μ1+μσr+α1+μC1lnr+α1+μC2;   (16)

εr=1+μE1-μσr-μσθ+α1+μC1lnr+α1+μC2;     (17)

εθ=1+μE1-μσθ-μσr+α1+μC1lnr+α1+μC2.     (18)

Будем искать решение данной задачи в следующем виде:

σr=A1r2+A21+2lnr+2A3;                             (19)

σθ=-A1r2+A23+2lnr+2A3;                            (20)

εθ=B1r2+B21+2lnr+2B3;                             (21)

εr=-B1r2+B23+2lnr+2B3.                            (22)

Подставим формулы (19)-(22) в обобщенный закон Гука (17)-(18):

-B1r2+B23+2lnr+2B3==1+μE1-μA1r2+A21+2lnr+2A3--μ-A1r2+A23+2lnr+2A3+α1+μC1lnr+α1+μC2;  (23)

B1r2+B21+2lnr+2B3==1+μE1-μ-A1r2+A23+2lnr+2A3--μA1r2+A21+2lnr+2A3+α1+μC1lnr+α1+μC2;  (24)

-B1r2+2B2lnr+3B2+2B3==1+μEA1r2+21-2μA2lnr+1-3μA2++21-2μA3+α1+μC1lnr+α1+μC2;            (25)

B1r2+2B2lnr+B2+2B3==1+μE-A1r2+21-2μA2lnr+3-4μA2++21-2μA3+α1+μC1lnr+α1+μC2;            (26)

B1=-1+μEA1;                                             (27)

2B2=21+μ1-2μA2+α1+μC1;                      (28)

B2=1+μ1-2μA2E+α21+μC1;                       (29)

3B2+2B3=1+μE1-4μA2+21+μ1-2μEA3+α1+μC2;(30)

B2+2B3=1+μE3-4μA2+21+μ1-2μEA3+α1+μC2;(31)

3E1+μ1-2μA2+α21+μC1==1+μE3-4μA2+21+μ1-2μEA3+α1+μC2-2B3;  (32)

1+μ1-2μA2+α21+μC1==1+μ3-4μEA2+21+μ1-2μEA3+α1+μC2-2B3;  (33)

31+μ1-2μA2E+32α1+μC1==1+μ3-4μA2E+21+μ1-2μA3E+α1+μC2-2B3;    (34)

2B2=21+μ1-2μA2E+α1+μC1;                         (35)

B2=1+μ1-2μA2E+α1+μ2C1;                           (36)

3B2+2B3=1+μ1-4μEA2+21+μ1-2μEA3+α1+μC2;  (37)

B2+2B3=1+μ3-4μEA2+21+μ1-2μEA3+α1+μC2;  (38)

31+μ1-2μA2E+α1+μ2C1+2B3==1+μ1-4μEA2+21+μ1-2μEA3+α1+μC2;       (39)

1+μ1-2μEA2+21+μ1-2μEA3+α1+μC2+2B3==1+μ3-4μEA2+21+μ1-2μEA3+α1+μC2;   (40)

31+μ1-2μEA2-1+μ1-4μEA2==21+μ1-2μEA3+α1+μC2-32α1+μC1-2B3;       (41)

1+μ1-2μEA2-1+μ3-4μEA2==21+μ1-2μEA3+α1+μC2-α21+μC1-2B3;        (42)

A21+μE31-2μ-1-4μ=21+μ1-2μEA3++α1+μC2-32α1+μC1-2B3;        (43)

A21+μE1-2μ-3-4μ=21+μ1-2μEA3++α1+μC2-32α1+μC1-2B3;       (44)

21+μA21-μE=21+μ1-2μEA3++α1+μC2-32α1+μC1-2B3;                    (45)

21+μA2μ-1E=21+μ1-2μEA3++α1+μC2-α21+μC1-2B3.                    (46)

Приравняем соотношения (45) и (46):

21+μ1-2μEA3+α1+μC2-32α1+μC1-2B3==-21+μ1-2μEA3-α1+μC2+α21+μC1+2B3;      (47)

41+μ1-2μEA3=-2α1+μC2+2α1+μC1+4B3;      (48)

A3=E2α1+μC1-2α1+μC2+4B341+μ1-2μ.                      (49)

Используем соотношение (45):

2A21-μ2E=21+μ1-2μEA3++α1+μC2-32α1+μC1-2B3;                            (50)

A2=E21-μ221+μ1-2μEA3++α1+μC2-32α1+μC1-2B3.                           (51)

Проверим выполнение уравнения (46):

21+μ1-2μEA3+α1+μC2-32α1+μC1-2B3==-21+μ1-2μEA3+α1+μC2-α21+μC1-2B3;       (52)

41+μ1-2μEA3=-2α1+μC2+2α1+μC1+4B3.          (53)

Соотношение (53) совпадает с соотношением (48).

Подсчитаем величину A2:

A2=A31-μ+αEC221-μ-34αEC11-μ-EB31-μ2.                 (54)

Подсчитаем коэффициент B2 по  формуле (36):

B2=1+μ1-2μEA31-μ+αEC221-μ-34αEC11-μ-EB31-μ2++α21+μC1;  (55)

B2=1+μ1-2μE1-μA3+1+μ1-2μEαEC221-μ--34αEC11-μ-EB31-μ2;           (56)

B2=1+μ1-2μE1-μA3+1+μ1-2μαC221-μ-3αC14+B31+μ;(57)

B2=1+μ1-2μE1-μE2α1+μC1-2α1+μC241+μ1-2μ++1+μ1-2μE1-μ4EB3+1+μ1-2μ1-μαC22-3αC14+B31+μ;   (58)

B2=1+μ1-2μ2α1+μ41-μ1+μ1-2μC1-C2++1+μ1-2μ41-μB3+1+μ1-2μ1-μαC22-3αC14+B31+μ;   (59)

B2=α1+μ2C1-C2+1+μ1-2μ1-μB3+1+μ1-2μ1-μαC22-3αC14+B31+μ.                  (60)

Запишем выражение для тангенциальной и радиальной деформаций:

εθ=B1r2+B21+2lnr+2B3;                                 (61)

εr=-B1r2+B23+2lnr+2B3.                                (62)

Однозначное радиальное перемещение будем находить по следующей формуле [2, 9, 11-14]:

ur=12εθr+εrdr.                                     (63)

Подставим формулы (61) и (62) в соотношение (63):

ur=12rB1r2+B2+2B2lnr+2B3+

+-B1r2+3B2+2B2lnr+2B3dr=

=12B1r+B2+2B3r+2B2rlnr+

+B1r+3B2+2B3r+2B2rlnr-r=

=122B1r+4B2+4B3r-2B2r+4B2lnr=

=122B1r+2B2+2B3r+4B2lnr;                            (64)

ur=B1r+B2+2B3r+2B2lnr.                              (65)

Пусть заданы граничные условия в перемещениях:

ur=r1=u1;     ur=r2=u2;                               (66)

B1r1+B2+2B3r1+2B2lnr1=u1;B1r2+B2+2B3r2+2B2lnr2=u2.                               (67)

Используем формулу (60):

B2=aB3+b;                                                  (68)

a=1+μ1-2μ1-μ;                                           (69)

b=α1+μC1-C22+1+μ1-2μ1-μαC22-3αC14+B31+μ;    (70)

B1r1+aB3+br1+3B3r1+2aB3+blnr1=u1;              (71)

B1r2+aB3+br2+2B3r2+2aB3+blnr2=u2;              (72)

B1r1+ar1+2r1+2alnr1B3=u1-br1-2ablnr1;          (73)

B1r2+ar2+2r2+2alnr2B3=u2-br2-2ablnr2;         74

B1+ar12+2r12+2ar1lnr1B3=u1r1-br12-2abr1lnr1;    (75)

B1+ar22+2r22+2ar2lnr2B3=u2r2-br22-2abr2lnr2;   (76)

B3ar22-r12+2r22-r12+2ar2lnr2-r1lnr1==u2r2+br12+2abr1lnr1-u1r1-br22-2abr2lnr2;        (77)

B3=ar22-r12+2r22-r12+2ar2lnr2-r1lnr1u2r2+br12+2abr1lnr1-u1r1-br22-2abr2lnr2;     (78)

B1=r1u1-br1-2ablnr1-r1ar1+2r1+2alnr1;         (79)

A1=-E1+μB1.                                           (80)

Коэффициенты A2 и A3 находятся по формулам (49) и (53).

Пусть граничные условия заданы в напряжениях:

σrr=r1=-p;     σrr=r2=-q.                            (81)

Имеем выражения для напряжений:

σr=A1r2+A21+2lnr+2A3;                              (82)

σθ=-A1r2+A23+2lnr+2A3.                             (83)

Подсчитаем величину A2 по формуле (51):

A2=1-2μ1-μA3+αEC221-μ-34αC1E1-μ-EB31-μ2;            (84)

A2=1-2μ1-μA3+E1-μαC22-34αC1-B31+μ.            (85)

Подставим в формулу (82) соотношение (49):

A2=1-2μ1-μ2αE1+μC1-C241+μ1-2μ+1-2μ1-μEB31+μ1-2μ++E1-μαC22-34αC1-B31+μ;        

A2=αEC1-C221-μ+E1-μαC22-34αC1-B31+μ+EB31-μ2=

=αEC1-C221-μ+αE21-μC2-32C1=αE21-μC1-C2+C2-32C1=

=-αEC141-μ;     A2=αEC141-μ.                                 (86)

Используем формулу (80):

σr=A1r2+A2+2A3+2A2lnr;                                 (87)

A1r12+A2+2A3+2A2lnr1=-p;                                (88)

A1r22+A2+2A3+2A2lnr2=-q;                               (89)

A11r22-1r12+2A2lnr2r1=p-q;                              (90)

A1=r12r22p-q2A2lnr2r1r12-r22;                                    (91)

2A3=-p-2A2lnr1-A2-A1r12;                             (92)

A3=-p2-A2lnr1-A22-A12r12;                             (93)

B1=-1+μA1;                                              (94)

B2=1+μ1-2μEA2+α1+μ2C1;                       (95)

A3=C1-C221-2μ+B31+μ1-2μ;                        (96)

B31+μ1-2μ=A3-C1-C221-2μ;                        (97)

B3=1+μ1-2μA3-1+μC1+C22.                (98)

Таким образом, все полученные выше формулы для напряжений, деформаций и радиального перемещения содержат в явной форме теплофизические константы C1 и C2.

 

2 РЕШЕНИЕ ЗАДАЧИ В НАПРЯЖЕНИЯХ

Рассмотрим алгоритм решения данной задачи в напряжениях.

Деформации имеют следующий вид:

εr=1+μE1-μσr-μσθ+α1+μC1lnr+α1+μC2;    (99)

εθ=1+μE1-μσθ-μσr+α1+μC1lnr+α1+μC2.    (100)

Введем в рассмотрение потенциал напряжений [3, 5]:

σr=φr;     σθ=dr.                                              (101)

Потенциал напряжений удовлетворяет уравнению равновесия (1).

Деформации примут следующий вид:

εr=1+μE1-μφr-μdr+α1+μC1lnr+α1+μC2;    (102)

εθ=1+μE1-μdr-μφr+α1+μC1lnr+α1+μC2.  (103)

Уравнение совместности деформаций:

rdεθdr+εθ-εr=0.                                           (104)

Подставим деформации (99) и (100) в уравнение совместности деформаций (104):

1+μErddr1-μdr-μφr+1+μEdr-φr+α1+μC1r;  (105)

r1-μd2φdr2+μrφr2-μdr+dr-φr=-αEC1r;               (106)

1-μd2φdr2+1-μrdr+μ-1r2φ=-αEC1r2;                 (107)

d2φdr2+1rdr-φr2=-αEC11-μr2.                            (108)

Решение уравнения (108) имеет вид:

φr=D1r+D2r+D3;                                       (109)

-D3=-αEC11-μ;     D3=αEC11-μ;                             (110)

φr=D1r+D2r+αEC11-μ;                                   (111)

σr=φr=D1+D2r2+αEC11-μr;                                (112)

σθ=dr=D1-D2r2.                                         (113)

Имеет следующие граничные условия:

σrr=r1=-p;     σrr=r2=-q;                          (114)

D1+D2r12+αEC11-μr1=-p;                                  (115)

D1+D2r22+αEC11-μr2;                                       (116)

D21r22-1r12=p-q-αEC11-μ1r22-1r12;                    (117)

D2r12-r22r12r22=p-q+αEC11-μ1r12-1r22;                     (118)

D2r12-r22r12r22=p-q+αEC1r22-r12r12r22;                        (119)

D2=p-qr12r22+αEC1r22-r12r12-r22;                         (120)

D1=-p-αEC11-μr1-D2r12.                                  (121)

Из формул (109), (112), (113), (120) и (121) следует, что напряжения σr и σθ не зависят от физической константы C2, что нежелательно.

 

3 РЕШЕНИЕ ЗАДАЧИ В ПЕРЕМЕЩЕНИЯХ

Обобщенный закон Гука при температурно-влажностном воздействии имеет следующий вид:

εr=1+μE1-μσr-μσθ+α1+μC1lnr+α1+μC2;     (122)

εθ=1+μE1-μσθ-μσr+α1+μC1lnr+α1+μC2.     (123)

Сложим и вычтем соотношения (122) и (123):

εr+εθ=1+μ1-2μEσθ+σr+2α1+μC1lnr+2α1+μC2; (124)

εr-εθ=1+μEσr-σθ;                                 (125)

σr-σθ=E1+μεr-εθ;                                 (126)

1+μ1-2μEσθ+σr=εr+εθ-2α1+μC1lnr-2α1+μC2;  (127)

σr+σθ=Eεr+εθ1+μ1-2μ-2αE1+μ