CC BY
42
УДК 539.3
Вестник СПбГУ. Сер. 1. Т. 3(61). 2016. Вып. 1
О ПОСТРОЕНИИ МАТРИЦЫ ЖЕСТКОСТЕЙ ЕСТЕСТВЕННО-ЗАКРУЧЕННОГО СТЕРЖНЯ*
Н. В. Курбатова, Н. М. Романова, Ю. А. Устинов
Южный федеральный университет,
Российская Федерация, 344090, Ростов-на-Дону, ул. Мильчакова, 8а
В статье рассматриваются задачи Сен-Венана для естественно-закрученного стержня с прямоугольным поперечным сечением. Исследования задач проводятся на основе метода однородных решений в совокупности с методом конечных элементов. Зависимости элементов матрицы жесткостей от крутки и напряженно-деформированное состояние сечений представлены графически. Библиогр. 17 назв. Ил. 3.
Ключевые слова: естественно-закрученный стержень, задача Сен-Венана, матрица жест-костей.
Первые подходы построения решений задач Сен-Венана о растяжении, кручении и изгибе естественно-закрученных стержней (ЕЗС) были сделаны в работах [1, 2]. Попытки исследования задач изгиба без привлечения дополнительных гипотез не привели к существенным результатам. Однако потребность в решении конкретных практических задач обусловила развитие прикладных теорий [3-7].
В работах [8, 9] методом однородных решений и спектральной теории операторов решение задачи Сен-Венана представлено в виде бесконечного ряда по элементарным решениям (ЭР), в котором выделена группа из двенадцати ЭР. Шесть из них описывают смещение ЕЗС как твердого тела, а напряженное состояние в поперечном сечении остальных шести ЭР таково, что их главный вектор и главный момент отличны от нуля и могут быть уравновешены усилиями, приложенными к торцам стержня. Совокупность таких элементарных решений названа «решением Сен-Венана». Это понятие включает в себя шесть элементарных решений: решения Сен-Венана задачи растяжения-сжатия, задачи кручения, задачи чистого изгиба и изгиба поперечными силами. Построение решений каждой из перечисленных задач сведено к двумерным краевым задачам на поперечном сечении ЕЗС.
Заметим, что для построения аналогичных решений для призмы Сен-Венан использовал полуобратный метод [10, 11]. В [9] для цилиндра строго доказано, что главный вектор и главный момент напряжений, отвечающие остальным решениям бесконечного ряда (не входящим в выделенную группу ЭР), равны нулю, и эта часть общего решения локализуется у торцов стержня, а обоснование принципа Сен-Венана для ЕЗС приводится в [12].
Аналитические решения всех шести задач Сен-Венана построены в виде двух членов разложений по малому параметру безразмерной «крутки» то = кт, т — относительный угол закручивания; к — характерный линейный размер сечения. В первом приближении при некоторых дополнительных упрощениях эти решения совпадают с решениями теории Кирхгофа—Клебша.
* Работа выполнена при частичной поддержке проекта «Математическое моделирование неоднородных и многофазных структур» (в рамках Программы фундаментальных исследований по стратегическим направлениям развития науки Президиума РАН №1 «Фундаментальные проблемы математического моделирования»).
(¡5 Санкт-Петербургский государственный университет, 2016
Система уравнений пространственной теории упругости сведена к системам двумерных уравнений на сечении ЕЗС [13], и для малых круток в случае эллиптического сечения сделаны оценки эффективных жесткостей.
В работах [14, 15] методом конечных элементов (МКЭ) исследована задача Сен-Венана растяжения-кручения ЕЗС, в [16]—задача о чистом изгибе, а в работе [17] строится решение краевой задачи, отвечающее задаче изгиба поперечной силой.
В настоящей работе на основе решения Сен-Венана для ЕЗС с прямоугольным поперечным сечением строится матрица жесткостей в широком диапазоне изменения величины крутки т, 0 < т < ж.
1. Пусть xi,x2,хз —декартова система координат с ортами ii, Í2, Í3. Объем ЕЗС получается винтовым движением плоской фигуры S параллельно плоскости Х1Х2 вдоль оси Х3. Угол поворота сечения, находящегося на расстоянии x по оси Х3 от начала координат, равен р = тх. В сопутствующей системе координат £i,£2, £з = £ оси £i и £2 жестко связаны с S в процессе ее движения и направлены по главным осям инерции. Связь между координатами обеих систем определяется соотношениями
£i = xi cos р + Х2 sin р, £2 = —xi sin р + Х2 cos р, £ = хз.
Орты сопутствующей системы координат обозначим ea, а = 1, 2. Таким образом,
R = £a.ea + £Í3 (1)
— радиус-вектор произвольной точки объема ЕЗС. Уравнение контура dS области S £а0 = £а0 (s), £а0 — координаты точки контура и s — расстояние от некоторой фиксированной точки по дуге контура. Тогда касательный t и нормальный n орты контура dS выражаются следующим образом:
t = £'а 0 ea, n = Па ea = £2 0 e1 — £1 oe2, (2)
где £í0,£2o — производные £10,£20 по s. Проекции внешней нормали N к боковой поверхности Г ЕЗС на оси сопутствующей системы координат таковы:
Ni = ni = £2 0, N2 = П2 = — £i 0, N3 = тЪп d
= п2а - га1£2 = -
Будем считать, что 0 ^ х ^ Ь, где Ь — длина ЕЗС.
Для постановки и исследования задачи Сен-Венана приведем основные соотношения теории упругости в сопутствующей системе координат и, опираясь на соотношения (1), (2), введем основной и взаимный базисы:
дК Ж .
В-а = дт- = еа, 11з = -т— = 13 - т£2е 1 + т£1в2, д£а ох
К1 = е1 + т£21з, И2 = е2 - т£х 1э, И3 = 1э = е3.
Учитывая, что метрический тензор g в смешанном базисе совпадает с единичным (д3 = К • И.3 = ), и используя следующее выражение:
V = И3 83 = еада + ез(д + тП), (3)
д д
уравнение равновесия представляется в векторном виде
дааа + (д + тВ)а = 0, (4)
а а = ааев£, а = аз£в£, I = 1, 2, 3, а = 1, 2.
Проектируя (4) на оси сопутствующей системы координат, получаем
дхаи + д2а 12 + дст!з + т(Вахз - а2з) = 0,
д\о 12 + д2 022 + да2з + т (Ба2з + ахз) = 0, (5)
дхахз + д2 023 + дазз + тВазз = 0. Для компонент тензоров деформаций и напряжений имеем
£ц = <91^1, £12 = + 8^2), £22 = д2и2,
£13 = ^[дщз + (а + тВ)ьл - ти2], (6)
£23 = + (д + тВ)и2 + тих], £зз = (9 + тВ)и3.
он = 2^[(1 + к)дх их + кд2 «2 + к(д + тВ)пз], а 12 = ¡¡(дп + д2«х),
а 22 = 2^[(1 + к)д2 «2 + кдх пх + к(д + тВ)пз], ахз = ¡¡(дхпз - тп2 + (д + тВ)пх), а 2з = ¡¡(д2«з + тпх + (д + тВ)п2), азз = пх + кд2 «2 + (1 + к)(д + тВ)пз].
Будем считать, что боковая поверхность ЕЗС свободна от напряжений, т. е.
N • а
(7)
г=0, (8)
а на торцах при г = 0 и при г = Ь заданы следующие граничные условия:
и(0) = 0, а(Ь) = р.
2. Уравнение равновесия (5) и граничное условие (8) могут быть представлены в операторном виде:
Ьхт(д)и = [ьт(д)и, Ит(д)и} = 0, (9)
где операторы Ьт , Мт определены в [9]. Отыскивая решение в виде
и(£) = е^а, а = а(& ,6) (10)
и подставляя (10) в (9), получаем спектральную задачу на сечении:
Ьт х(7)а ={Ьт (7 )а, Мт (7)а> = 0. (11)
Вестник СПбГУ. Сер. 1. Математика. Механика. Астрономия. Т. 3 (61). 2016. Вып. 1 149
В [9] показано, что задача (11) имеет три четырехкратных собственных значения 70 = 0, 7±1 = ±гт, им соответствуют 12 элементарных решений Сен-Венана:
U1 (О = (О = ai, U3 (£) = u4 (£) = (£ai + аз), us(£)=a5, u6(£)=a6, u7(£) = eiT? i —ai + ya3 + £a9+a7
U8(0=U7(0, u9(e) = ui0(e) = e-«^a1+ea3 + a9^, (12)
uii(£) = £a5 + an, U12 (£) = £ae + ai2,
ai = a2 = {1,i,0}, a3 = a4 = {0, 0, -C}, a5 = {0,0,1}, ae = {-£2,£i, 0}, ag = a7, aio = a9,
()* —комплексное сопряжение. Решение Сен-Венана представим в виде
e i2
us(£) = 53 Qu(£) Ceue(£ - L).
Здесь первая сумма содержит шесть элементарных решений (12), определяющих смещение ЕЗС как твердого тела, а вторая — определяет напряженно-деформированное состояние типа Сен-Венана. В [9] сформулированы краевые задачи для определения ад, а7, ац и а12 и соответствующие им вариационные постановки.
Для координат векторов а; и компонент тензоров напряжений, соответствующих элементарным решениям с индексом I, используем обозначения , атп,1. Тогда
*тп,1 =0, I = 1,..., 6,
атп,7 = (£Ътп,9 + Ътп,7), атп,8 = а тп,7, , ,
_ _ *
атп,9 bmn,д, атп, 10 а mn,д,
атп,11 ^Ъmn,11, атп, 12 Ътп,12.
Подстановка (13) в соотношения обобщенного закона Гука (7) приводит к выражениям для следов напряжений ЪтпI = 7,..., 12 [9]. Компоненты главного вектора Р = и главного момента М = {М^1,М^2,М3} напряжений = {а 13, о"23,*33}, действующих в сечении Б, определяются следующим образом:
д* = / ^ ^ = / Ь*33^
Б Б
М52 = - J £lаззdS, М3 = J(£1*23 - £2Vlз)dS,
(14)
где — проекции векторов напряжений в сечении £ = const на оси сопутствующей системы координат.
Используя (14) и соотношения обобщенной ортогональности, а также выражения для следов напряжений [9], получаем ненулевые элементы матрицы жесткостей
dii — ¿22 — —dss — —¿44 — —d — j Z*bss,QdS, d^i — dl2 — — J Z*bss,rdS,
S S
d55 — J b33,iidS, d56 — d65 — J b33,udS, dee — J(£ib23,i2 — &bi3,i2)dS.
S S S
(15)
Постоянные Cg, i — 7,..., 12, определяются из следующей алгебраической системы:
—eiTn^dC7 — Qa — iQp,
eiTn M[(d3i + nd)C7 + dCg] — M^ + iM^,
l^(d55Cii + d56 Ci2 ) — Q3,
n(d56Cii + deeCi2) — M3, ц — £ — L.
3. На основе вариационной постановки задач, соответствующих задачам растяжения-кручения, чистого изгиба и изгиба поперечной силой, с помощью МКЭ в случае прямоугольного поперечного сечения ЕЗС была создана программа для определения координат векторов a7, ag, ац и ai2.
Особенности реализации МКЭ были отражены в [15-17].
На рис. 1 изображены графики жесткостей ЕЗС при различных значениях т. В задаче на растяжение-сжатие для d5j5 характерно увеличение жесткости с ростом т. Жесткость d56 является характеристикой взаимодействия кручения и растяжения-сжатия; с ростом т меняется знак d56, что отражает эффекты «раскручивания-закручивания» ЕЗС.
8 6 4 2 О -2 -4 -6
-Чб
/ \ 1 \
1 1 \ \
\ \ \
ч N
Ч
02468 10 12 0246
Рис. 1. Зависимости жесткостей от крутки в задаче растяжения-кручения
8 10.
На рис. 2 приведены графики жесткостей dзl, dзз при Н2/Н1 = 1 (кривые Ь и d соответственно) и при Н2/Н1 = 2 (кривые а и с соответственно), Ь = 1 м, Н = 0.1 м.
На рис. 3 представлено распределение нормированных напряжений <3з = <733/^ в сечении хз = 0.5235 м при изгибе поперечной силой ^2 = 100 Н с максимальным значением <3з, равным 2.7397 • 10~9. Расчеты показали, что с ростом то, при то > 4 происходит формирование круга, вписанного в сечение, для которого, как и для малых то, свойственна линейная зависимость величины нормального напряжения, пропорциональная расстоянию от нейтральной линии и имеющая противоположные знаки в точках, симметричных относительно нее. Это, на наш взгляд, подтверждает
1
0.9 0.8 0.7 0.6 D 0.5 0.4 0.3 0.2 0.1 0
* * ч \ * * (а)
V \ ......(Ъ)
! \ \ —■(с)
■ ! у \ 5 \ \ / \ Ч. » - ч (d)
............'
0 2 4 6 8 10 12 14 16 18 20 22 24 г
Рис. 2. Зависимости жесткостей стержня от крутки задачи изгиба
"33
1,: ol -11 -0.5
0.5 0.5
Рис. 3. Распределение по сечению нормальных напряжений 033: слева т = 0.02; справа т =
гипотезу автора [13] о том, что с ростом т работает ядро сечения, которое стягивается ко вписанному в сечение кругу, что является следствием быстрой осцилляции боковой поверхности ЕЗС.
Литература
1. Риз П. М. Деформация естественно закрученных стержней // Докл. АН СССР. 1939. Т. 23, №1. С. 18-21.
2. Джанелидзе Г.Ю., Лурье А. И. Задачи Сен-Венана для естественно скрученных стержней // Докл. АН СССР. 1939. Т. 24, №4. С.325-326.
3. Уздалев А. И., Иноземцев Г. Г., Зубков А. В., Алахазова О. В. Напряженное состояние естественно закрученного стержня // Прикладная механика. 1988. Т. 24, №14. С. 103-108.
4. Елисеев В. В. Задача Сен-Венана и упругие модули для стержней с кривизной и кручением // Изв. АН СССР, МТТ. 1991. №2. С. 167-176.
5. Заметалина Н.П., Прокопов В. К. Напряженное состояние естественно скрученных стержней типа спиральных сверл // Изв. АН Арм. ССР. 1974. Т. 27, №3. С. 3-9.
6. Корольков В. И. К решению задачи о растяжении естественно закрученного стержня произвольного поперечного сечения в трехмерной постановке // Прикладная механика. 1988. Т. 4, №12. С. 113-115.
7. Шорр Б. Ф. К теории закрученных неравномерно нагретых стержней // Изв. АН СССР, ОТН. Механика и машиностроение. 1960. 1. С. 141-151.
8. Друзь А. Н., Поляков Н. А., Устинов Ю. А. Однородные решения и задачи Сен-Венана для естественно закрученного стержня // ПММ. 1996. Т. 60, №4. С. 657-664.
9. Устинов Ю. А. Задачи Сен-Венана для псевдоцилиндров. М.: Наука, 2003.
10. Saint-Venant. Memoire sur la torsion des prisms, avec des conciderations sur leur flexion, ainsi que sur l'equilibre interieur des solides elastiques en general, et des formules pratiques pour le calcul de leur resistance a divers efforts s'exercant simultanement // Mem. Savants Etrang. 1856. XIV. P. 233-560.
11. Saint-Venant. Memoire sur la flexion des prismes // Liouville, J. Math. I. P. 89-189.
12. Устинов Ю. А. Обоснование принципа Сен-Венана для естественно-закрученного стержня // ВМЖ. 2010. Т. 12. Вып. 1. С. 58-71.
13. Бердичевский В. Л., Старосельский Л. А. Изгиб, растяжение и кручение естественно-закрученных стержней // ПММ. 1985. Т. 49, №6. С. 978-991.
14. Устинов Ю. А., Курбатова Н. В. Задачи Сен-Венана для стержней с физической и геометрической анизотропией // Изв. вузов. Сев.-Кавказ. регион. Естеств. науки. Мат. модел. 2001. Спецвыпуск. С. 154-157.
15. Kurbatova N. V. On a stretching-torsion of a naturally twisted rod // Proc. of XXXIII Summer School — Conference "Advanced Problems in Mechanics". St. Petersburg. Russia. 2005. P. 59-60.
16. Курбатова Н.В., Романова Н.М. Конечно-элементное решение задачи изгиба для естественно-закрученного стержня // Тр. IX конф. «Совр. пробл. механики сплош. среды». Ростов-на-Дону: ЦВВР. 2005. Т. 1. С. 123-126.
17. Курбатова Н.В., Устинов Ю.А., Чумакова Е. С. Анализ напряженно-деформированного состояния естественно закрученного стержня при изгибе поперечной силой на основе метода конечных элементов // Владикавказ. мат. журн. 2013. Т. 15, вып. 3. С. 45-53.
Статья поступила в редакцию 22 октября 2015 г. Сведения об авторах
Курбатова Наталья Викторовна — кандидат физико-математических наук, доцент; nvk-ru@yandex.ru, nvk@math.sfedu.ru
Романова Наталья Михайловна — ассистент; sunshinemail@yandex.ru
Устинов Юрий Анатольевич — доктор физико-математических наук, профессор;
ustinov@math.sfedu.ru
ON STIFFNESS MATRIX CONSTRUCTION FOR NATURALLY TWISTED ROD
Natalia V. Kurbatova, Natalia M. Romanova, Yuriy A. Ustinov
Southern Federal University, ul. Milchakova, 8a, Rostov-on-Don, 344090, Russian Federation; nvk-ru@yandex.ru, nvk@math.sfedu.ru, sunshinemail@yandex.ru, ustinov@math.sfedu.ru
Here we considered the Saint Venant problem for a naturally twisted rod with a rectangular cross-section. Research was conducted on the basis of homogeneous solutions in conjunction with the finite element method. The general solution is constructed as a linear combination of elementary solutions corresponding to three four-roots eigenvalues of spectral problem on the cross section. Elementary solutions determining the stress-strain state of Saint Venant type contain unknown eigenvectors and associated vectors. To determine of unknown solutions the authors formulated boundary value problems and their variational formulations previously. They are correspond to problems of stretching-torsion, pure bending and bending of the lateral force. Variational problems were solved using FEM. The stress-strain state of the rod was studied numerically and non-zero elements of the stiffness matrix were found in the case of square and rectangular cross-sections of the rod for the different values of twist. Graphically shows the numerical results for a wide range of change of the twist parameter t. Calculations showed that the identified patterns are consistent with the corresponding behavior untwisted rods (for small twist parameter), and with the growth of the twist, new effects, which confirm the hypothesis proposed earlier. Refs 17. Figs 3. Keywords: naturally twisted rod, Saint-Venaunt problem, stiffness matrix.
References
1. Reese P. M., "The deformation of naturally twisted rods", Dokl. AN USSR 23(1), 18-21 (1939) [in Russian].
2. Dzhanilidze G.Yu., Lurie A. I., "Saint-Venaunt problems for naturally twisted rods", Dokl. AN USSR 24(4), 325-326 (1939) [in Russian].
3. Uzdalev A. I., Inozemtsev G.G., Zubkov A. V., Alakhazova O.V., "The stress state of a naturally twisted rod", Int. Appl. Mech. 24(14), 103-108 (1988) [in Russian].
4. Yeliseyev V. V., "Saint-Venaunt problem and elastic modulus for rods with curvature and torsion", Izv. AN ASSSR, Mech. Solids (2), 167-176 (1991) [in Russian].
5. Zametalina N.P., Prokopov V. K., "The stress state of a naturally twisted spiral rod", Izv. AN Arm. SSR. 27(3), 3-9 (1974) [in Russian].
6. Korolkov V. I., "On stretching problem solution of a naturally twisted rod for arbitrary cross section in 3-D case", Int. Appl. Mech. 4(l2), 113-115 (1988) [in Russian].
7. Shorr B.F., "For the theory of twisted irregularly heated rods", Izv. AN USSR, J. Mech. and Machin. 1, 141-151 (1960) [in Russian].
8. Druz A. N., Polyakov N. A., Ustinov Yu.A., "Homogeneous solutions and Saint-Venant problems for a naturally twisted rod", J. Appl. Math. Mech. 60(4), 657-664 (1996) [in Russian].
9. Ustinov Yu.A., Saint-Venant Problem for Pseudo-Cylinder (Nauka, Moscow, 2003) [in Russian].
10. Saint-Venant, "Memoire sur la torsion des prisms, avec des conciderations sur leur flexion, ainsi que sur l'equilibre interieur des solides elastiques en general, et des formules pratiques pour le calcul de leur resistance a divers efforts s'exercant simultanement", Mem. Savants Etrang. XIV, 233-560 (1856) [in France].
11. Saint-Venant, "Memoire sur la flexion des prismes", Liouville, J. Math. I, 89-189 [in France].
12. Ustinov Yu.A., "Proof of Saint-Venant principle for a naturally twisted rod", VMJ 12(1), 58-71 (2010) [in Russian].
13. Berdichevskii V. L., Starosel'skii L. A., "Bending, extension, and torsion of naturally twisted rods", J. Appl. Math. Mech. 49(6), 746-755 (1985).
14. Ustinov Y. A., Kurbatova N. V., "The Saint-Venant problems for rods with physical and geometrical anisotropy", Izvestia of Universities. North-Caucas. Region. Math. Model. Natural Science. Special Issue 154-157 (2001) [in Russian].
15. Kurbatova N. V., "On a stretching-torsion of a naturally twisted rod", Proc. of XXXIII Summer School — Conference "Advanced Problems in Mechanics" (St. Petersburg, Russia, 59-60, 2005).
16. Kurbatova N. V., Romanova N. M., "Finite element solution of the bending problem for a naturally twisted rod", Proc. of IX intern. conference "Modern Problems of Continuum Mechanics" 1, 123-126 (Rostov-on-Don, 2005) [in Russian].
17. Kurbatova N. V., Ustinov Yu. A., Chumakova E. S., "Analysis of stress-strain state of the naturally twisted rod bending by transverse force on the basis of the finite element method", Vladikavkazian Math. J. 15(3), 45-53 (2013) [in Russian].
