MATHEMATICAL MODEL OF STIFFNESS PROPERTIES OF THIN-WALLED BARS OF CLOSED PROFILE OF THE SHIP STRUCTURES
Authors:
1.
Astrakhan State Technical University
(Doctor of Technical Sciences, Professor; Professor of the Department of Technique and Technology of Land Transport)
Astrakhan, Russian Federation
Astrakhan, Russian Federation
2.
Astrakhan Institute of Civil Engineering
Type:
Article
Pages:
from 41
to 52
Status:
Published
Received:
20.05.2016
Accepted:
25.05.2016
Published:
25.05.2016
Subject area:
UDK 539.4(076.5)
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
Language:
Russian
Keywords:
ship hull, thin-walled spatial bar of closed profile, stiffness matrix, mathematical model, finite element method
Abstract (English):
The calculation analysis of ship, crane and other engineering structures consisting of thin-walled bars of closed profile, based on the theory of thin-walled welded and hot rolled bars, still remains the subject of the research. The theoretical bases for building a mathematical model of thin-walled bar of closed profile in the local coordinate system at the spatial deformation taking into account the shift of the middle surface are given. The mathematical correlations for constructing the stiffness matrix of thin-walled spatial bar of closed profile, which can be used for static and dynamic analysis of ship structures calculated by finite element method, are obtained.
The calculation analysis of ship, crane and other engineering structures consisting of thin-walled bars of closed profile, based on the theory of thin-walled welded and hot rolled bars, still remains the subject of the research. The theoretical bases for building a mathematical model of thin-walled bar of closed profile in the local coordinate system at the spatial deformation taking into account the shift of the middle surface are given. The mathematical correlations for constructing the stiffness matrix of thin-walled spatial bar of closed profile, which can be used for static and dynamic analysis of ship structures calculated by finite element method, are obtained.
Keywords:
ship hull, thin-walled spatial bar of closed profile, stiffness matrix, mathematical model, finite element method
ship hull, thin-walled spatial bar of closed profile, stiffness matrix, mathematical model, finite element method
Введение Корпус судна проектировщики ранее считали абсолютно жестким, недеформируемым телом. В теории корабля интерес представляла только внешняя форма обводов корпуса, соответственно, изучалось влияние этой формы на отдельные мореходные качества. В настоящее время в строительной механике корабля [1] корпус плавающего судна рассматривается как тонкостенная балка переменного по длине коробчатого сечения (рис. 1, а). а б в г д Рис. 1. Тонкостенные элементы корабельных конструкций замкнутого профиля: тонкостенная балка по длине судна коробчатого сечения (а); днищевые перекрытия из коробчатых балок поперечной системы набора (б, в); продольной системы набора (г); коробчатое сечение корпуса танкера по бортовому отсеку (д) Под действием совокупности внешних и внутренних сил различной природы эта балка деформируется, в ее элементах возникают нормальные и касательные напряжения, следовательно, она должна обладать достаточной прочностью, которую можно определить как несущую способность корпуса воспринимать, не разрушаясь, нагрузки, возникающие в процессе эксплуатации судна. Корпус должен обладать и достаточной жесткостью, т. е. его деформации должны быть относительно невелики и не оказывать влияния на мореходные качества судна. Обеспечение прочности и жесткости корпуса при наименьшей затрате материала - одна из основных задач, решаемых при строительстве судна. Кроме корпуса судна, как тонкостенной балки коробчатого сечения, отдельные элементы корпуса - перекрытия - представляют собой близкие к плоскостным конструкции, составленные листами обшивки и балками набора. Участки обшивки, опирающиеся на балки, имеют коробчатую прямоугольную форму. Кроме того, поперечные и продольные системы набора днищевых перекрытий (рис. 1, б-г) представляют собой тонкостенные стержни замкнутого профиля. Особенно часто замкнутые профили применяются для поперечных сечений корпусов танкеров по бортовому отсеку (рис. 1, д). Следует отметить, что основным этапом расчета сложных пространственных машиностроительных конструкций произвольного вида, составленных из тонкостенных стержней замкнутого профиля методом конечных элементов в перемещениях, является формирование матриц жесткости и масс отдельного стержня, находящегося в условиях пространственного деформирования в местной системе координат (МСК). Математическая модель тонкостенного стержня замкнутого профиля Следуя методике, изложенной в [2], посвященной стержням открытого профиля, получим матрицу жесткости тонкостенного стержня замкнутого профиля. Положительные направления узловых перемещений и усилий в стержне указаны на рис. 2, согласно которому компоненты вектора перемещений (по узлам j и k) примут вид . (1) Рис. 2. Правило знаков для узловых перемещений тонкостенного стержня замкнутого профиля В (1) - линейные перемещения узла; - углы поворота; - производная от угла закручивания (депланация); подстрочные индексы x, y и z обозначают оси МСК; надстрочные индексы j и k указывают на начало (н) и конец (к) стержневого конечного элемента (КЭ); Т - индекс транспонирования матриц и векторов. Вектору перемещений (1) соответствует вектор внутренних усилий в дискретных граничных узлах j и k: , (2) где и - внутренние поперечные и продольная силы; и - изгибающие и крутящий момент; B - изгибно-крутящий бимомент. Правила знаков для внутренних усилий (2) аналогичны правилу знаков для узловых перемещений КЭ (1). Стержень считаем элементом с четырнадцатью степенями свободы, пространственное положение которого определяется вектором обобщенных координат (1). Значения элементов искомой матрицы жесткости при прочностном расчете зависят от жесткостных параметров КЭ и принятого закона изменения компонентов перемещений, в качестве которых примем аппроксимирующие функции Эрмита [3, 4]: (3) В выражениях (3) qs - узловое перемещение, соответствующее вектору (1); ys(z) - аппроксимирующие функции Эрмита, которые для тонкостенного стержня замкнутого профиля с двумя осями симметрии, жестко защемленного по концам, с учетом сдвига срединной поверхности принимают вид [5]: (4) . (5) При выводе выражений (4) использовались результаты работы [5], согласно которой функции , являются приближенными. Кроме того, поскольку функции в (4) получены для тонкостенного стержня открытого профиля [6], для тонкостенного стержня замкнутого профиля в (4) и (5) необходимо произвести замену . (6) Таким образом, матрица жесткости для пространственного тонкостенного стержня замкнутого профиля имеет размер . Чтобы построить элементы этой матрицы, воспользуемся вариационными принципами, для чего будем следовать методике, изложенной в работе [7]. Выпишем формулы для нормальных и касательных напряжений в стержне ; (7) , (8) причем в (8) полный крутящий момент Mz будем определять как сумму , (9) в которой момент чистого кручения H выражается через угол закручивания Θ формулой , (10) а бимомент В в (7) и момент стесненного кручения в (8), согласно полусдвиговой теории А. А. Уманского [8], связываются с мерой депланации β (11) соотношениями (12) Кроме того, в (11) , m - коэффициента депланации; - распределенная бимоментная нагрузка, зависящая от координаты z; - направленный момент инерции сечения (рис. 3), определяемый по формуле ; (13) - крутильный момент инерции сечения стержня однозамкнутого профиля, принимается, как в формуле (10): , (14) где - удвоенная площадь, охватываемая срединной линией поперечного сечения стержня (рис. 3). а б Рис. 3. Распределение касательных напряжений в тонкостенных стержнях замкнутого профиля при свободном кручении: общий вид обобщенного профиля (а); замкнутый силовой многоугольник погонного касательного усилия (б) Дополнительно укажем, что в (7), (8) и (12) - секториальный момент инерции сечения стержня замкнутого профиля, который имеет вид , где d, как и ранее, - толщина контура поперечного сечения тонкостенного стержня; определяется по формуле . (15) Введенная формулой (15) функция дуговой координаты называется обобщенной секториальной координатой и играет важную роль во всей теории тонкостенных стержней замкнутого профиля, роль которой такая же, как роль обычной секториальной координаты в теории стержней открытого профиля [2, 7]. Обратим внимание на одно существенное свойство обобщенной секториальной координаты из (15), а именно: обобщенная секториальная координата не зависит от выбора точки начала отсчета дуговой координаты s (рис. 3). Укажем, что величина в (15) характеризует распределение депланационных перемещений вдоль сечения замкнутого профиля при чистом кручении тонкостенного стержня. Именно поэтому именуют функцией депланации (11), в которой принимается согласно формуле (14), а - момент инерции чистого кручения сварного стержня, составленного из вертикалов и полок с параметрами : , где . Для получения матрицы жесткости тонкостенного стрежня замкнутого профиля вычислим потенциальную энергию внутренних сил через напряжения σ и τ из (7) и (8): . (16) Подстановка (7) и (8) в (16), с учетом (9) и (10), после преобразований дает [9] . (17) Коэффициенты формы сечения в (17) определяются по формулам ; ; ; ; (18) ; . В (17) и (18) и - это полярный момент инерции и полярный радиус инерции стержня соответственно; - геометрическая характеристика сечения стержня, определяемая по формуле [9]: . (19) Дополнительно укажем, что входящие в состав (18) и (19) функции, как и из (15), ; ; , (20) не зависят от выбора начала отсчета дуговой координаты s и являются однозначной характеристикой всякой точки, принадлежащей профилю поперечного сечения стержня (рис. 4). Рис. 4. Равновесие выделенного элемента тонкостенного стержня замкнутого профиля Как следует из рис. 4, в (19) и (20) представляют собой площадь и статические моменты отсеченной части сечения, расположенной между начальной точкой сечения О и текущей точкой М с координатой S. Величина называется секториальным статическим моментом отсеченной части сечения тонкостенного стержня замкнутого профиля (рис. 4). Очевидно, что . Пусть ξ, η, ζ - перемещения точек линии центров изгиба стержня в направлении осей x, y и z его МСК. Подставим (4) с заменой коэффициентов по (5) и (6) в (3) и учтем известные дифференциальные зависимости в которых - направленный момент инерции (13), а также формулы (12) и, подставляя (3) в (17), после преобразований получим потенциальную энергию внутренних сил [7] тонкостенного стержня замкнутого профиля с учетом сдвига срединной поверхности: (21) Далее представим выражение для потенциальной энергии деформации (21) квадратичной формой обобщенных перемещений , которое можно записать в виде двойной суммы . (22) Из (22) следует, что любой элемент матрицы жесткости можно получить как вторую производную от V из (21): . (23) Формирование матрицы жесткости тонкостенного стержня замкнутого профиля Определение элементов матрицы жесткости по формуле (23) позволило сформировать матрицу жесткости тонкостенного стержня замкнутого профиля в МСК при пространственном деформировании с учетом сдвига срединной поверхности: , (24) симметрично ненулевые элементы которой представлены формулами: где приняты следующие обозначения: . (25) Выводы В заключение отметим, что для целей практических расчетов корабельных, крановых и других машиностроительных конструкций при формировании матриц жесткости полной системы со многими степенями свободы матрицу (24) отдельного стержня следует представить в блочном виде, для случая, в расчетной модели конструкции, когда номер узла k не следует за номером узла j: , (26) в которой, как известно, . Кроме того, при использовании матрицы (26) следует учитывать возможность получения из формул (25) матрицы жесткости конечного элемента массивного стержня размером , широко представленную, в том числе в нормативной литературе [10], для чего в формулах (25) следует исключить влияние на ее компоненты как стесненного кручения, так и сдвига срединной поверхности, обусловленного касательными напряжениями (8).
1. Zhinkin V. B. Teoriya i ustroystvo korablya: uchebnik / V. B. Zhinkin. SPb.: Sudostroenie. 2002. 336 s.
2. Yuzikov V. P. Raschet tonkostennyh sterzhney otkrytogo profilya s uchetom sdviga sredinnoy poverhnosti / V. P. Yuzikov, O. B. Zav'yalova // Izv. vuzov. Stroitel'stvo. 2011. № 1. S. 108-115.
3. Klaf R. Dinamika sooruzheniy / R. Klaf, Dzh. Penzien. M.: Stroyizdat, 1979. 320 s.
4. Sinel'schikov A. V. Dinamika i seysmostoykost' mostovyh kranov: dis. … kand. tehn. nauk / A. V. Sinel'schikov. Astrahan': AGTU, 2000. 276 s.
5. Semenov P. I. Raschet prochnosti i deformativnosti anizotropnyh tonkostennyh sterzhney otkrytogo profilya / P. I. Semenov. Kiev: Vyscha shk., 1974. 184 s.
6. Panasenko N. N. Konechno-elementnaya model' prostranstvennyh konstrukciy iz tonkostennyh sterzhney otkrytogo profilya. V 2-h chastyah. Ch. 1 / N. N. Panasenko, V. P. Yuzikov, A. V. Sinel'schikov // Vestn. Astrahan. gos. tehn. un-ta. Ser.: Morskaya tehnika i tehnologiya. 2015. № 2. S. 89-100.
7. Yuzikov V. P. Stroitel'naya mehanika tonkostennyh sterzhney / V. P. Yuzikov, N. N. Panasenko; pod red. d-ra tehn. nauk N. N. Panasenko. Volgograd: Volgograd. nauch. izd-vo, 2013. 361 s.
8. Umanskiy A. A. Stroitel'naya mehanika samoleta / A. A. Umanskiy. M.: Oborongiz, 1961. 529 s.
9. Slivker V. I. Stroitel'naya mehanika. Variacionnye osnovy: ucheb. posobie / V. I. Slivker. M.: Izd-vo Associacii stroit. vuzov, 2005. 736 s.
10. RD 24-090-83-87. Normy rascheta prostranstvennyh metallokonstrukciy gruzopod'emnyh kranov atomnyh stanciy na ekspluatacionnye i seysmicheskie vozdeystviya
