ITERATIVE METHOD OF CALCULATION FOR THIN-WALLED ANISOTROPIC CORES OF OPEN PROFILE ON TORSION

Research article
DOI:
https://doi.org/10.23670/IRJ.2017.55.091
Issue: № 1 (55), 2017
Published:
2017/01/25
PDF

Аллахвердов Б.М.1, Полинкевич К.Ю.2

1Кандидат технических наук, доцент,

Петербургский государственный университет путей сообщения императора Александра I г. Санкт-Петербург, Россия

2Аспирант,

Петербургский государственный университет путей сообщения императора Александра I г. Санкт-Петербург, Россия

ИТЕРАЦИОННЫЙ СПОСОБ РАСЧЕТА ТОНКОСТЕННЫХ АНИЗОТРОПНЫХ СТЕРЖНЕЙ ОТКРЫТОГО ПРОФИЛЯ НА КРУЧЕНИЕ

Аннотация

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

Ключевые слова: теория упругости, анизотропия, метод итераций, напряжения, деформации.

Allakhverdov B.M.1, Polinkevich K.Y.2

1 PhD in Engineering, lecturer,

Petersburg State University of Railways of Emperor Alexander I St. Petersburg, Russia

2 Postgraduate student,

Petersburg State University of Railways of Emperor Alexander I St. Petersburg, Russia

ITERATIVE METHOD OF CALCULATION FOR THIN-WALLED ANISOTROPIC CORES OF OPEN PROFILE ON TORSION

Abstract

The article provides the description of the calculation technique which determines the strength of thin-walled anisotropic cores of an open profile based on an iterative method of successive accordance with the conditions of equilibrium and deformations compatibility. The first cycle of iterations allowed us to obtain V.Z.Vlasov solution. In subsequent cycles, we obtained the solution which takes into account the anisotropy, the deformation of the shift and normal transverse deformation. The question of the results convergence is studied in the article. The example based on modern carbon fiber illustrates how the material characteristics differ depending on the direction of the reinforcing fibers.

Keywords: elasticity theory, anisotropy, iterative method, stress, deformation.

Задача о кручении изотропного тонкостенного стержня открытого профиля без учета сдвигов решена В.З. Власовым [1]. Для уточнения этого решения используется метод итераций [2], [3] основанный на последовательном удовлетворении условиям равновесия и совместности деформаций. Аналогичный путь решения  распространен на стержни, выполненные из анизотропных материалов [4].

Предварительный этап

Будем считать, что поперечное сечение тонкостенного стержня открытого профиля состоит из n анизотропных пластин толщиной ti, жестко соединенных между собой вдоль длинных ребер. (рис.1). В качестве основной принимается гипотеза недеформируемости контура поперечного сечения.

Стержень закручивается заданным крутящим моментом М. При этом функция углов закручивания 24-01-2017 16-43-30 поперечных сечений по длине стержня подлежит определению.

24-01-2017 16-52-52

Рис. 1.

а) Поперечное сечение тонкостенного стержня открытого профиля

(до и после поворота на угол φ).

б) Один из анизотропных элементов сечения (главные оси анизотропии под углом α к оси yi)

Для всего стержня выбирается произвольная глобальная система координат Х, Y. Координата Z направлена вдоль оси стержня (рис.1.а). В то же время для каждого i –го элемента вводится местная координатная система xi, yi, zi (рис.1.б). При кручении любое сечение стержня с координатой z поворачивается на угол 24-01-2017 16-34-19 относительно центра кручения Р, имеющего координаты Хр, Yр подлежащие определению. При таком повороте каждый элемент перемещается как в своей плоскости, так и из плоскости.  Перемещения Vi в плоскости yi, zi соответствуют изгибу элемента, перемещения же из плоскости вызваны кручением этой пластины.

Перемещения Vi отдельной пластинки можно описать через координаты центра кручения Xp, Yp и угол закручивания поперечных сечений относительно этого центра 24-01-2017 16-34-19:

24-01-2017 16-35-15                        (1)

где ri – нормаль к оси i –ой пластинки, проведенная из точки P.

 Таким образом кривизна i элемента в его плоскости становится функцией:

24-01-2017 16-35-27                                                  (2)

Считается, что пластины выполнены из различных анизотропных материалов, упругие свойства которых обычно задаются в главных осях анизотропии у, z (вдоль и поперек волокон материала). Обобщенный закон Гука для случая плоского напряженного состояния принимает следующий вид:

24-01-2017 16-35-44                                        (3)

Здесь Ey, Ez, Gyz - соответственно модули упругости и сдвига;

εz, εy – относительные удлинения (линейные деформации) в направлении действия нормальных напряжений;

γyz - относительный сдвиг - величина изменения прямого угла между площадками, на которых действуют соответствующие касательные напряжения;

νzy, νуz  – коэффициенты  влияния  линейной деформации по направлению z  на линейную деформацию по направлению у и наоборот, при этом соблюдается отношение .

ησ,y, ησ,z , ηt,y, ηt,z – коэффициенты  влияния  линейной деформации на  сдвиговую деформацию   по направлениям z, у и наоборот; при этом соблюдаются отношения 24-01-2017 16-37-31 и 24-01-2017 16-37-40.

В случае, если при заданных для некоторой системы координат у, z упругих постоянных требуется найти упругие постоянные для новой системы y1, z1, повернутой по отношению к первой на угол α (рис.1 б), модули и коэффициенты для новых осей определяются по формулам, приведенным в [2]:

24-01-2017 16-38-01

На рис.2 показаны графики зависимости модуля упругости и модуля сдвига углепластика марки M60J/Epoxy [5] от угла поворота осей анизотропии.

image018

Рис.2 - Полярные диаграммы изменения модулей углепластика M60J/Epoxy при повороте координатных осей:

а) продольной упругости;

б) сдвига

Цикл 1.

  1. Первое приближение (первый шаг итераций).
1.1. Принимаем на первом этапе для каждой i –ой пластины анизотропный материал со следующими упругими характеристиками: Ezi (та же величина, что и в реальном материале); остальные модули стремятся к бесконечности Eyi= ∞; Gyzi=∞. В этом случае закон Гука имеет вид:

24-01-2017 16-39-29                                                      (5)

1.2. При таком виде закона Гука 24-01-2017 16-39-39,

что соответствует гипотезе о ненадавливании волокон.

1.3. Тогда 24-01-2017 16-39-49. 1.4. Также из (5) 24-01-2017 16-39-58. 1.5. Тогда 24-01-2017 16-40-08, 1.6. и  24-01-2017 16-40-17(6)

Получена известная в сопротивлении материалов зависимость - гипотеза плоских сечений (отдельно для каждого элемента), где 24-01-2017 16-40-27и 24-01-2017 16-40-36 неизвестные функции. Здесь 24-01-2017 16-40-27 можно трактовать как продольную деформацию оси стержня, а 24-01-2017 16-40-48- как кривизну его оси.

При этом на гранях пластин должны соблюдаться условия совместности деформаций: 24-01-2017 16-41-05 - так выглядят эти условия между двумя соседними пластинами (индекс к означает деформацию волокна 24-01-2017 16-41-12, индекс н - . деформацию волокна 24-01-2017 16-41-27).

24-01-2017 16-41-40- между второй и третьей и т.д.

1.7. Продольная деформация оси любой пластины выражается через продольную деформацию оси первой из них (похожая идея высказана в [3]):

24-01-2017 16-41-50                             (7)

Выражение (7) содержит четыре неизвестных величины: Xp, Yp, 24-01-2017 16-42-01  Для их определения следует составить глобальные уравнения равновесия.

1.8. Первое уравнение 24-01-2017 16-42-11 позволяет определить функцию 24-01-2017 16-42-18(z).

1.9. Два других уравнения 24-01-2017 16-42-30 служат для определения координат центра кручения Xp, Yp.

Здесь следует учесть, что изгибающий момент в плоскости каждой пластины записывается так:

24-01-2017 16-42-41                (8)

1.10. Возникающие в элементе поперечные силы представляются в следующем виде:

24-01-2017 16-42-51,                (9)

1.11. Поперечные силы создают часть крутящего момента, вызванного изгибом элементов (то, что у Власова [1] названо моментом стесненного кручения):

24-01-2017 16-42-59                  (10)

Здесь 24-01-2017 16-43-07  - плечо силы по (2).

1.12. Вторая составляющая крутящего момента связана с свободным кручением каждой пластины:

24-01-2017 16-43-14                                                          (11)

1.14.  Уравнение равновесия 24-01-2017 16-43-22, из которого можно найти функцию угла закручивания 24-01-2017 16-43-30, становится дифференциальным уравнением:

24-01-2017 16-43-41                               (12)

где 24-01-2017 16-43-49  - аналог секториальной жесткости по [1], 24-01-2017 16-43-59-   крутильная жесткость, что повторяет уравнение В.З. Власова [1].

Решив его при заданных граничных условиях, зависящих от способов закрепления стержня, получаем функцию углов закручивания 24-01-2017 16-43-30.

На этом заканчивается определение неизвестных первого этапа.

1.15. Далее на этом этапе окончательно определяется функция нормальных напряжений в каждом элементе 24-01-2017 16-44-10.

1.16. Из уравнения равновесия 24-01-2017 16-44-20 определяется функция касательных напряжений:

24-01-2017 16-44-30.                             (13)

Здесь, согласно теореме Бредта, поток касательных напряжений на границе элементов сохраняется т.е. 24-01-2017 16-44-39. Исходя из этого условия, и зная, что на свободных гранях отсутствуют касательные напряжения, определяются произвольные функции 24-01-2017 16-44-47

1.17. Из второго уравнения равновесия 24-01-2017 16-44-55 определяется функция поперечных нормальных напряжений:

24-01-2017 16-45-06.                               (14)

Теперь, если по теореме Бредта, поток касательных напряжений на границе элементов сохраняется, то также должен сохраняться поток поперечных нормальных напряжений: 24-01-2017 16-45-25.   Исходя из этого условия, и зная, что на свободных гранях отсутствуют эти напряжения, определяются произвольные функции 24-01-2017 17-05-20.

Следует обратить внимание на следующие факты – в этом цикле:

  • Продольные нормальные напряжения по высоте –линейные функции в пределах каждого элемента.
  • Продольные нормальные напряжения по длине –являются функциями, зависящими от 24-01-2017 16-45-33.
  • Касательные напряжения по высоте –квадратные параболы в пределах элемента.
  • Касательные напряжения по длине –являются функциями, зависящими от 24-01-2017 17-05-57.
  • Поперечные нормальные напряжения по высоте –кубические параболы в пределах элемента.
  • Поперечные нормальные напряжения по длине –являются функциями, зависящими от 24-01-2017 16-45-42.

На определении напряженного состояния заканчивается первый цикл итераций.

Цикл 2.

На втором цикле остается та же гипотеза неизменяемости контура и все предпосылки первого цикла. За неизвестные также принимаются:

  • Функция углов закручивания стержня 24-01-2017 16-47-30
  • Координаты центра кручения Хр, Yр.
  • Продольная деформация оси элементов 24-01-2017 16-46-01
  • Кривизна оси элементов 24-01-2017 16-46-08

Здесь и далее (2) –принадлежность второму циклу, (1) –принадлежность первому.

Матрица модулей принимает обычный вид по (5):

24-01-2017 16-46-18

В этом варианте учтено влияние сдвигов и поперечных нормальных напряжений на общее напряженное состояние стержня при кручении.

В начале второго цикла по закону Гука для анизотропного тела (5) находим для каждого элемента функцию деформаций вдоль оси x -  24-01-2017 16-46-26:

24-01-2017 16-46-35                      (15)

Зная функцию деформации 24-01-2017 16-46-26, можно найти соответствующую функцию перемещения 24-01-2017 16-46-44, пользуясь уравнением связи между деформациями и перемещениями:

24-01-2017 16-46-52                                                         (16)

с точностью до произвольной функции 24-01-2017 16-47-01.

Далее определяется функция деформации сдвига по обобщенному закону Гука (5):

24-01-2017 16-47-10          (17)

Затем вновь находим функцию продольных деформаций 24-01-2017 17-09-44 по процедуре, описанной в первом цикле:

  24-01-2017 16-47-22          (18)

Продольная деформация (17) состоит из двух частей: последние два слагаемых повторяют линейное распределение этой деформации по высоте элемента на первом цикле и зависят от производных функции 24-01-2017 16-47-30, первые же слагаемые создают более высокий порядок распределения по высоте, соответствующий напряжениям, определенным на первом этапе - они зависят от производных функции 24-01-2017 16-47-41.

Все дальнейшие вычисления ведутся по последовательности, изложенной выше для первого цикла, причем для величин, зависящих от 24-01-2017 16-47-30, решение сводится к тому же виду дифференциального уравнения (12). Однако величины, зависящие от уже определенной функции 24-01-2017 16-47-41, будучи перенесенными в правую часть этого уравнения, играют роль увеличения или уменьшения величины заданного крутящего момента. Таким образом уравнение (12) трансформируется в следующее:

24-01-2017 16-47-56                 (19)

с теми ж граничными условиями, как и на первом цикле.

Дальнейшие уточнения решения (следующие циклы) проводятся в той же последовательности.

Как следует из вышеизложенного, на каждом этапе при вычислении напряжений происходят операции дифференцирования по координате z и интегрирования по координате у. Такой подход к задаче определения напряжений приводит результату, представленному в табл.1.

Табл. 1

24-01-2017 16-32-11

На этих фактах базируется сходимость данного итерационного процесса -  при увеличении циклов до числа m напряженное состояние по длине стержня l зависит от возрастания степени производных функции 24-01-2017 16-34-19, т.е. приращение напряжений пропорционально величине 1/l 2m. С другой стороны напряженное состояние по высоте стержня h с доведением числа циклов до m характеризуется увеличением влияния высоты элементов в виде h2m. Следовательно на каждом цикле можно оценивать добавку к первоначальному решению пропорциональной величине (h/l)2m.

При отношении h/l=1/3 на втором цикле поправка к напряжениям 24-01-2017 16-48-18 пропорциональна величине (1/3)4 =0.012, а на третьем уже 0.15-4. Процесс вычисления напряжений имеет высокую скорость сходимости, обычно достаточно 2 – 3 циклов, чтобы получить необходимую точность решения.

Список литературы / References

  1. Власов В.З., Тонкостенные упругие стержни, Стройиздат, 1940 г.
  2. Аллахвердов Б.М. Итерационный метод расчета балок с изменяющимися по высоте характеристиками. Исследования по механике материалов и конструкций.(сб. научн.статей)/Вып.12/ Петерб. Гос. Универ. Путей сообщ. –СПб,2002.-С.30 –34, Деп. ВИНИТИ.№ 1400-В2002.
  3. Полинкевич К.Ю. Итерационный способ расчета слоистых балок на прочность. Известия Петербургского университета путей сообщения. /Вып 2 (35)/ Петерб. Гос. Универ. Путей сообщ. –СПб,2013.-С.148 –153
  4. Аллахвердов Б.М. Полинкевич К.Ю. Итерационный способ расчета анизотропной балки на прочность. Известия Петербургского университета путей сообщения. /Вып 2 (39)/ Петерб. Гос. Универ. Путей сообщ. –СПб,2014.-С.73 –79
  5. Попов А.А. Сопротивление материалов (теория и задачи). Государственное научно-техническое издательство машиностроительной литературы. Москва, 1956, стр. 476.
  6. Лехницкий С.Г. Анизотропные пластинки. Государственное издательство технико-теоретической литературы. Москва, 1957, стр. 463.

Список литературы на английском языке / References in English

  1. VZ Vlasov, Tonkostennie uprugie stergni [Thin walls elastic beams], Stroiizdat [construction publishing], 1940.
  2. BM Allakhverdov. Iteracionnii metod rascheta balok s izmenyaushimisya po visote harakteristikami [Iterative method of calculation of beams with variable height characteristics]. Issledovaniya po mechanike materialov I konstrukcii (sb. nauchn. statei)/Vip 12/ Peterb. Gos. Univer. Putei soobch. [Research on the mechanics of materials and structures.(collected scientific articles)/Edd.12/ Petersburg. The State University of Communications. St. Petersburg]. –SPb,2002.-P.30 –34, Dep. VINITI.№ 1400-В2002.
  3. KY Polinkevich. Iteracionnii sposob rascheta sloistich balok na prochnost’[An iterative method of calculation of layered beams for strength] Izvestiya Peterburgskogo universiteta putei soobcheniya [Proceedings of the St. Petersburg University of Railways]. / Issue 2 (35) / Petersburg State University of Railways. SPb, 2013.-p.148 -153
  4. BM Allakhverdov, KY Polinkevich. Iteracionnii sposob rascheta anizotropnoi balki na prochnost’[The iterative method of calculation of the anisotropic beam strength] Izvestiya Peterburgskogo universiteta putei soobcheniya [Proceedings of the St. Petersburg University of Railways]. / Issue 2 (39) / Petersburg State University of Railways. SPb, 2014.-p.73 -79
  5. AA Popov. Soprotivlenie materialov (teoriya I zadachi) [Strength of Materials (theory and tasks).]. Gosudarstvennoe nauchno-technicheskoe izdatel’stvo mashinostroitel’noi literature [State scientific and technical publishing engineering literature]. Moscow, 1956, P. 476.
  6. SG Lekhnickii. Anizotropnie plastinki [Anisotropic plates]. Gosudarstvennoe izdatel’stvo techniko-teoreticheskoi literature [State Publishing House technical and theoretical literature]. Moscow, 1957, P. 463.