ВВЕДЕНИЕ
Проектирование и расчет напряженно-деформированного состояния упругих и упруго-вязко-пластических цилиндров из сжимаемого материала имеет важное научно-практическое значение [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=EαC1-C221-2μ+B31+μ1-2μ; (96)
B31+μ1-2μ=A3-EαC1-C221-2μ; (97)
B3=1+μ1-2μA3-Eα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; σθ=dφdr. (101)
Потенциал напряжений удовлетворяет уравнению равновесия (1).
Деформации примут следующий вид:
εr=1+μE1-μφr-μdφdr+α1+μC1lnr+α1+μC2; (102)
εθ=1+μE1-μdφdr-μφr+α1+μC1lnr+α1+μC2. (103)
Уравнение совместности деформаций:
rdεθdr+εθ-εr=0. (104)
Подставим деформации (99) и (100) в уравнение совместности деформаций (104):
1+μErddr1-μdφdr-μφr+1+μEdφdr-φr+α1+μC1r; (105)
r1-μd2φdr2+μrφr2-μdφdr+dφdr-φr=-αEC1r; (106)
1-μd2φdr2+1-μrdφdr+μ-1r2φ=-αEC1r2; (107)
d2φdr2+1rdφdr-φ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)
σθ=dφ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+μ