On the specifics of behavior of the sandwich plate composite facing layers under local loading

Abstract


The problem of a four-point bending of sandwich plates with external layers of a fiber reinforced plastic is considered, the results of numerical and experimental studies are presented. It is shown that in the vicinity of the loading roller, which exerts a local effect on the external layer, there is a strong decrease in the transverse shear secant modulus of fiber reinforced plastic with an increase in the transverse shear strains. The numerical solution of the problem of the plate bending is carried out in a physically and geometrically nonlinear formulation using various relations of the finite element method, two variants of geometrically nonlinear kinematic relations of the equations of the elasticity theory and different variants of the loading process parameter. Along with the classical nonlinear relations, the problem solutions are also constructed on the basis of consistent relations between strains and displacements, the use of which allows one to avoid the appearance occurrence of false bifurcation points. The results of the numerical calculations obtained with different methods and use of different ratios are presented, the analysis of which showed their small difference. It was revealed that the stability loss and failure of the external layers of a sandwich plate occurs due to the stability loss in the nonclassical transverse-shear mode. To determine the ultimate load, which is accompanied with a loss of strength of the external loaded layer, the Tsai-Wu criterion was used. A comparative analysis of the behavior of the plate external layers at different thicknesses and different diameters of the loading roller is carried out. It is shown that the ultimate load is practically not affected by the roller diameter, while the load at which the external layer loses its stability, is very sensitive to a change in its value.

Full Text

Введение Трехслойные элементы в виде стержней, пластин и оболочек имеют широкое применение в качестве конструктивных элементов в изделиях судостроительной и аэрокосмической отраслей. Как правило, такие элементы состоят из двух жестких несущих слоев, воспринимающих тангенциальные напряжения, и относительно маложесткого заполнителя, воспринимающего напряжения поперечного обжатия и поперечные касательные напряжения. Выбор материалов для несущих слоев и заполнителя зависит от специфики работы трехслойных конструкций. В изделиях аэрокосмической техники, как правило, применяется сотовый заполнитель на основе легких алюминиевых сплавов или полимерной бумаги. Такие заполнители обладают наиболее высокими удельными жесткостными и прочностными характеристиками при малом удельном весе. В трехслойных конструкциях, применяющихся в кораблестроении, используется пористый металлический заполнитель. От таких заполнителей требуется высокая удельная ударная вязкость и высокая способность поглощать энергию удара. В качестве материалов для несущих слоев используются как металлы, так и композитные материалы на основе стекло- и углепластиков. Разрушение трехслойных конструкций может быть обусловлено реализацией разных процессов деформирования их составных элементов в процессе их нагружения [1-4]: разрушение заполнителя по причине достижения предельных значений формирующимися в них поперечными касательными напряжениями [5, 6], усталостные разрушения заполнителя [7, 8]; отслоение несущих слоев от заполнителя [9, 10]; достижение предельных значений поперечными напряжениями обжатия в заполнителе [11]; потеря устойчивости несущих слоев по тем или иным формам [12, 13] или потеря их прочности. Из работ, посвященных действию локализованных нагрузок на пластины и оболочки, можно выделить монографию [14], в которой систематизирован обширный материал по концентрации напряжений от действия сосредоточенных нагрузок в тонкостенных конструкциях. При определении механических характеристик волокнистых композиционных материалов (ВКМ) применяются различные экспериментальные методы, для которых создано немало различных стандартов, зависящих от типа композитного материала (КМ). В отличие от металлических конструкционных материалов, для композитных материалов при этом требуется гораздо большее количество испытаний. Кроме того, как показывает анализ результатов экспериментальных исследований, в отличие от традиционных изотропных материалов, ВКМ проявляют специфические свойства при испытаниях на различные виды нагружения [15-17]. Поэтому при проектировании изделий из КМ коэффициенты запаса прочности принимаются значительно большими, чем в случае использования традиционных материалов. Одной из таких специфических особенностей поведения армированных КМ является потеря устойчивости фаз композита. В работах [18-21] установлено, что даже при простом сжатии стержня (волокна) может реализоваться одна из двух форм потери устойчивости: по неклассической поперечно-сдвиговой и классической изгибно-сдвиговой форме. Задача оценки устойчивости фаз композита усложняется тем, что эти фазы испытывают сложное напряженное состояние, причем вблизи закреплений и нагружающих элементов возникают сильные градиенты параметров напряженно-деформированного состояния (НДС). При анализе устойчивости фаз КМ важным является учет нелинейных свойств композита. Имеется много работ (см. например, в [22-29]), посвященных исследованиям процесса потери устойчивости конструкций при учете неупругих свойств композита. В них показано, что в ряде случаев потеря устойчивости может быть вызвана перераспределением напряжений, изменением геометрии тела и направления силового воздействия по причине наличия ползучести. При этом время потери устойчивости можно оценить обычным образом, анализируя возможность существования неединственного поля перемещений (бифуркации), принимая в этот момент упругий закон деформирования (ввиду практически мгновенной перестройки этого поля перемещений). В других случаях неустойчивость конструкций проявляется в том, что через некоторое время скорость деформации катастрофически быстро увеличивается. Таким образом, эти задачи становятся, с одной стороны, физически и геометрически нелинейными, а с другой - зависящими от времени. Однако для простоты нередко особенности поведения материала, обусловленные реологическими его свойствами, описываются моделями нелинейно-упругого деформирования. Различные виды потери устойчивости фаз композиционного материала являются факторами, приводящими или к мгновенным разрушениям, или к процессам приспособления, или к прогрессирующему накоплению микроповреждений и окончательному разрушению ВКМ, например, при циклических нагрузках [15-17]. В данной работе исследуется поведение трехслойной пластины с учетом возможной потери устойчивости или прочности внешних слоев на основе использования физически и геометрически нелинейных соотношений механики деформируемого твердого тела. Также считается, что для проведения качественного анализа поведения рассматриваемых конструкций можно вместо реологических соотношений использовать зависимости, описывающие нелинейно-упругое поведение материала внешних слоев как нелинейно-упругое тело. 1. Постановка задачи Основной целью работы является выявление механизма (причины) разрушения трехслойного тест-образца (пластины) при испытании на четырехточечный изгиб (рис. 1). Верхние и нижние несущие слои представляют собой ВКМ, армированные в направлении длинной стороны образца, внутренний слой является трансверсально-мягким заполнителем. Пластина опирается на верхние ролики, а через нижние ролики к ней прикладывается суммарная нагрузка 2P. Поскольку исследуются причины разрушения таких пластин, то в первом приближении можно заменить задачу деформирования с учетом деформаций ползучести на нелинейно-упругую задачу. Рис. 1. Схема нагружения трехслойного тест-образца на четырехточечный изгиб Fig. 1. Loading diagram of the sandwich specimen under four-point bending Зависимость между модулем сдвига и углом поперечного сдвига в несущих слоях закладывалась в виде физической нелинейности с учетом приведенной ниже диаграммы (рис. 2). Она была получена из экспериментов [17] на растяжение образцов с углами укладки волокон . Рис. 2. Экспериментальная зависимость секущего модуля от деформаций сдвига Fig. 2. Experimental dependence of the secant modulus on shear strain Полученная экспериментальная зависимость секущего модуля поперечного сдвига от деформации поперечного сдвига аппроксимировалась следующей функцией: (1) где - механические характеристики, получаемые минимизацией квадратичной невязки расчетных и экспериментальных значений секущего модуля, приведенного на рис. 2. Параметры, входящие в функцию (1), описывающую уменьшение секущего модуля поперечного сдвига, оказались следующими: ГПа, , . Длина контакта взаимодействия нагружающих и опорных роликов с несущими слоями тест-образца (рис. 3) считается постоянной, равной тому значению, при котором достигается разрушение нижнего несущего слоя или по причине потери прочности или потери устойчивости. При этом контактные давления аппроксимировались переменной по длине поверхностной нагрузкой p: (2) где - амплитуда нагрузки. Рис. 3. Схема контактного взаимодействия несущих слоев и нагружающих роликов Fig. 3. Diagram of the contact interaction of external layers and loading rollers 2. Основные соотношения и разрешающие уравнения Геометрические соотношения классической нелинейной теории упругости, построенные для случая малых деформаций и произвольных перемещений, имеют вид [33]: (3) (4) Как показано в работах [18-21], их использование при решении конкретных задач может привести к появлению «ложных» точек бифуркаций (в них дан детальный анализ соотношений (3), (4) и соответствующих им уравнений теории упругости и предложены непротиворечивые уравнения геометрически нелинейной теории упругости при малых деформациях). Поэтому, наряду с (3), (4), для решения рассматриваемой задачи были использованы полученные в [18-21] и геометрически корректные соотношения, в которых вместо (3) используются соотношения (5) а соотношения (4) остаются неизменными. Задача выявления возможных форм потери устойчивости внешних слоев из ВКМ в данной работе основана на использовании глобальной инкрементальной теории Лагранжа (Total Lagrange method) [34, 35]. Для дальнейшего его применения запишем основные соотношения теории деформирования в тензорном виде. Пусть текущая геометрия конструкции описывается уравнением , где X - радиус вектор, описывающий первоначальную геометрию тела (рис. 4). Рис. 4. Исходное и деформированное состояние сплошной среды Fig. 4. Initial and deformed state of the continuum Запишем тензор градиента деформации в виде , (6) компоненты которой будут: (7) Тогда тензор деформаций Грина можно представить в виде (8) где I - единичный тензор, а компоненты тензора деформаций Грина, соответствующие соотношениям (3)-(4), могут быть записаны в виде (9) Если вместо соотношений (3)-(4) в соответствии с [18-20] использовать геометрически корректные соотношения (5), (6), то тензор деформаций Грина (8) модифицируется следующим образом: (10) Здесь компоненты тензора I* имеют вид Вариационное уравнение Лагранжа представимо в виде (11) (12) Здесь - вариация работы внутренних сил; - вариация работы внешних сил; S - второй тензор напряжений Пиолы - Кирхгофа; f0 - объемные силы, - поверхностные силы; - вариация перемещений. Вариация тензора деформаций Грина (9) запишется в виде (13) Зависимость между тензорами напряжений Пиолы - Кирхгофа S и деформациями Грина E для физически нелинейного случая запишем в приращениях в виде где - тензор упругих характеристик материала (тензор касательных жесткостей). Вариационное уравнение (12) после линеаризации примет вид [35] (14) где (15) Для численного решения сформулированной задачи о четырехточечном изгибе трехслойной пластины с внешними слоями из ВКМ на основе вариационного уравнения (15) была разработана конечно-элементная модель с использованием трехмерных конечных элементов (КЭ) с учетом геометрической и физической нелинейности, в том числе и с возможностью учета соотношений (5). Рассматривались два типа КЭ: первый - восьмиузловой линейный изопараметрический КЭ [36], второй - модифицированный с дополнительным узлом, который «смягчает» элемент при изгибе [37]. Для двух типов КЭ были проведены численные эксперименты на тестовых задачах (растяжение и изгиб балки) при различном количестве элементов. Из анализа полученных результатов выявлено, что в задаче изгиба балки для получения решения, близкого к аналитическому, необходимо большое количество восьмиузловых КЭ первого типа. Это связано с тем, что такой элемент на изгиб довольно «жесткий». Аналогичный результат получился при решении задачи в геометрически нелинейной постановке. Поэтому для «смягчения» элемента использовались КЭ второго типа [37], в которых были введены дополнительные внутренние узлы. Это позволяло при небольшом количестве КЭ получить результаты, близкие к точным решениям как в задаче изгиба балки, так и при четырехточечном изгибе трехслойной пластины с внешними слоями из композиционного материала. 3. Методы расчета, результаты численных и физических экспериментов При использовании деформационного подхода под потерей устойчивости можно понимать два состояния. В первом случае под потерей устойчивости принимается состояние, при котором приращение перемещений или деформаций катастрофически быстро увеличивается при незначительном приросте нагрузки P. Во втором случае под потерей устойчивости понимается состояние, при котором на диаграмме «сила - перемещение» появляется экстремум. Нагрузку, при которой это происходит, в обоих случаях обычно называют критической. В соответствии с этими критериями при решении задачи о нахождении критической нагрузки для трехслойной пластины при четырехточечном изгибе были рассмотрены два подхода. Первый подход основан на глобальной инкрементальной теории Лагранжа [34, 35] (для краткости назовем его статическим), поскольку параметром процесса является величина нагрузки. При численном решении задачи потере устойчивости соответствует ситуация, в которой при приближении нагрузки к предельному значению решение задачи расходится. Во втором подходе использовалcя метод продолжения по другому параметру процесса нагружения, а именно, по длине дуги кривой «нагрузка - перемещение» (arc length method) [35, 38-40]. Ниже приведены результаты расчетов, полученные при использовании этих подходов. 3.1. Решение задачи статическим методом В качестве материала несущих слоев тест-образца принят однонаправленный ВКМ на основе углеленты ЭЛУР-П и эпоксидного связующего холодного отверждения. В численных расчетах использовались упругие характеристики, найденные в работе [41], и принятые равными: ГПа, ГПа, , а толщина несущих слоев см. Для трансверсально-мягкого заполнителя в виде сот из полимерной бумаги были приняты следующие характеристики: МПа, МПа, МПа, см. Расстояние между роликами принималось см при см, вся длина пластины равна 27 см, а ширина пластины b = 2,5 см. Ниже приводятся результаты расчетов трехслойного тест-образца с учетом геометрической и физической нелинейности (рис. 5-7). а б Рис. 5. Результаты расчета статическим методом: а - диаграмма «прогиб - нагрузка»; б - диаграмма «деформация сдвига нагрузка» Fig. 5. Calculation results by the static method: a is the deflection-load curve, b is the shear strain-load curve На рис. 5 графики построены для точки А, которая находится слева от приложенной нагрузки возле нагружающего ролика (см. рис. 3). Длина контактного взаимодействия нагружающих и опорных роликов с несущими слоями принята равной 0,3 cм. Картины напряженно-деформированного состояния на рис. 6-7 приведены для значения максимальной нагрузки, равной критической Pkr = 4,657 кН. Анализ полученной картины деформирования показал, что причиной неожиданного изменения процесса деформирования является скачкообразное изменение углов сдвига элементов над нагружающим роликом (см. рис 6-7). При этом зависимость перемещения ролика от нагрузки остается практически линейной, хотя зависимость угла сдвига от нагрузки в нижнем несущем слое, как видно из рис. 5, б, существенно нелинейна. Рис. 6. Картина деформирования и нормальные напряжения в правой половине тест-образца в момент потери устойчивости Fig. 6. Deformation pattern and normal stresses in the right half of the specimen at the buckling state Рис. 7. Картина деформирования и касательные напряжения τ13 в правой половине тест-образца в момент потери устойчивости Fig. 7. Deformation pattern and shear stresses τ13 in the right half of the specimen at the buckling state Из рис. 6-7 визуально обнаруживается, что потеря устойчивости происходит по сдвиговой форме в окрестности, близкой к нагружающему ролику. 3.2. Решение задачи методом дуги (arc length method) Решение задачи методом дуги основан на использовании метода продолжения решения по длине дуги [35, 38-40]. Основной идеей этого метода является то, что для получения решения на следующем шаге нагружения строится окружность (сфера) около найденного на данном шаге решения с заданным радиусом s (рис. 8, б). Уравнение окружности записывается в виде (16) Здесь - приращение перемещения на предыдущем шаге; δu - новое приращение перемещений; s - радиус окружности; - приращение параметра нагружения . а б Рис. 8. Схема метода продолжения по длине дуги Fig. 8. Scheme of the arc length method Точка пересечения окружности и диаграммы нагружения отыскивается путем решения системы уравнений (12), (16). Итерационно решение уточняется на каждом шаге параметра нагружения . Благодаря этому подходу нагрузка может увеличиваться или уменьшаться в зависимости от параметра нагружения , что позволяет обойти предельные точки А и В на рис. 8, а. С использованием описанного алгоритма были проведены численные расчеты при четырехточечном изгибе трехслойной пластины. Здесь, как и в статическом методе, бифуркация НДС приводит к резкому увеличению приращения угла сдвига (рис. 9). Хотя этот подход позволяет исследовать закритическое поведение конструкции (рис. 10, а), в рассматриваемой задаче не удается пройти далеко по дуге диаграммы от P. Оказалось, это вызвано тем, что в результате потери устойчивости внешнего слоя по поперечно-сдвиговой форме конечные элементы в области приложения нагрузки вырождаются, а именно: их площадь устремляется к нулю. Рис. 9. Картина деформирования и касательные напряжения τ13 в правой половине тест-образца в момент потери устойчивости (расчет методом дуги) Fig. 9. Deformation pattern and shear stresses τ13 in the right half of the specimen at the buckling state (the arc length method) а б Рис. 10. Результаты численных расчетов: а - зависимости деформаций сдвига от нагрузки (две методики); б - зависимости деформаций сдвига от нагрузки: 1 - стандартные геометрически нелинейные соотношения (3)-(4), 2 - непротиворечивые соотношения (5)-(6) Fig.10. Results of numerical calculations: a is dependence of shear strain on load (two methods); b is dependence of shear strain on load: 1 is classic geometric nonlinear relations of (3)-(4), 2 is consistent equations (5)-(6) Как видно из рис. 10, а, численные результаты расчета, представленные в виде зависимостей «нагрузка - сдвиг», по двум методикам близки в докритической области деформирования, метод дуги, в отличие от первого, позволяет проследить закритическое поведение конструкции. 3.3. Решение задачи на основе непротиворечивых кинематических соотношений Как показано в работах [18-20], использование соотношений (3)-(4) может в некоторых задачах привести к ложным бифуркационным решениям и неверным значениям определяемой критической нагрузки. Использование же непротиворечивых кинематических соотношений (5)-(6) позволяет в таких случаях получить реальные ее значения. Сравнение результатов применения обоих видов соотношений представлено на рис. 10, б, где приведены зависимости «нагрузка - сдвиговая деформация», полученные статическим методом. Видно, что в рассматриваемой задаче оба использованных подхода дают близкие результаты. 3.4. Исследование влияния геометрических параметров пластины на предельную нагрузку При проведении расчетов была принята упрощенная конечно-элементная схема, так как численная реализация задачи требует большого объема численных вычислений. Поэтому ширина тест-образца моделировалась одним слоем конечных элементов, задача решалась в постановке плоского напряженного состояния. Это упрощение может привести к тому, что не будет учтен кромочный эффект, а точное выявление сложного напряженного состояния возле нагружающих роликов требует большого объема вычислений. Такими особенностями можно пренебречь, если считать, что ширина несущих слоев тест-образца намного больше их толщины. Как было отмечено ранее, трехслойный тест-образец может разрушиться либо из-за потери прочности конструкции, либо из-за потери устойчивости несущего слоя в области нагружающего ролика. При этом от толщины несущего слоя трехслойной пластины зависит вид механизма, по которому происходит потеря несущей способности пластины. В связи с этим было проведено исследование влияния толщины несущего слоя, а также длины зоны контакта нагружающего ролика на механизм разрушения пластины. Для оценки прочности элемента несущего слоя пластины был принят критерий прочности Цая-Ву [42]. Для плоского напряженного состояния он запишется в виде (17) где Xt, Yt - пределы прочности на растяжение вдоль осей x, y соответственно; Xс, Yс - пределы прочности на сжатие вдоль осей x, y; S - предел прочности на сдвиг. Коэффициент F12 определяется при испытании на двухосное нагружение при (все остальные компоненты тензора напряжений равны нулю) и выражается формулой Однако из-за трудоемкости проведения испытания на двухосное растяжение коэффициент F12 обычно находят из испытаний на одноосное растяжение образцов из однонаправленных композитов с укладкой волокон [44]. Значение компонент будет составлять (остальные компоненты тензора напряжений равны нулю), где U - предел прочности образца при растяжении. Тогда При проведении расчетов коэффициент F12 был принят в более простом варианте, предложенном Хоффманом [44]: Кроме толщины несущего слоя на предельную нагрузку сильное влияние оказывают и геометрические размеры нагружающих элементов. Ролики могут быть разных диаметров, а область, на которую они действуют (она определяется отрезком , изображенным на рис. 3), меняется при нагружении. Оказалось, что в зависимости от величины отрезка несущий слой может или терять устойчивость, или достигать предела прочности. Пределы прочности для ВКМ несущих слоев были получены в работе [41]: Тогда На рис. 11 приведены результаты численных исследований, показывающие влияние отрезка и толщины несущего слоя t на предельную нагрузку. Было обнаружено (это видно и на рис. 11), что при малых толщинах несущих слоев (0,5-1,5 мм) пластина разрушается из-за потери прочности несущего слоя при сжатии, а при более толстых несущих слоях предел прочности не достигается, а происходит потеря устойчивости сжатого несущего слоя, причем, по чисто поперечно-сдвиговой форме. Можно отметить, что наибольший вклад в критерий прочности, как и ожидалось, вносят сжимающие напряжения. В момент достижения предела прочности в элементе над нагружающим роликом компоненты напряжений в сжатом несущем слое достигли следующих величин: При этом слагаемые в критерии (17) имеют значения: (18) t = 0,5 мм t = 1 мм t = 1,5 мм t = 2 мм t = 3 мм Рис. 11. Влияние толщины несущего слоя t и зоны контакта с0 на предельную нагрузку (синяя линия - предельная нагрузка, красная - критическая нагрузка) Fig. 11. Influence of the external layer thickness t and contact area c0 on the load limit (the blue line shows the ultimate load, the red line shows the critical load) Превышение функцией значения единицы вызвано, естественно, тем, что при численном решении задачи значение нагрузки меняется дискретно. Из (18) видно, что основной вклад в критерий прочности вносит сжимающее напряжение . Кроме этого, нужно отметить, что на предельную нагрузку, при которой функция f достигает значения 1, как это видно на рис. 11, практически не влияет длина контакта (т.е. значение отрезка , напрямую зависящего от диаметра ролика). Напротив, нагрузка, при которой происходит потеря устойчивости несущего слоя пластины, весьма чувствительна к изменению этой величины. Поэтому при испытаниях тест-образцов на четрырехточечный изгиб с целью определения пределов прочности на сжатие необходимо использовать либо тонкие несущие слои, либо ролики с большим диаметром. Иначе сжатые внешние слои тест-образцов могут потерять устойчивость раньше, чем будет достигнут предел прочности на сжатие. Это может привести в ряде случаев к неправильному определению его значения. Были проведены также численные эксперименты, в которых варьировалась величина a. Ожидаемо получилось, что предельная и критическая нагрузки повышались с увеличением a, что объясняется уменьшением изгибающего момента при одной и той же нагрузке. 3.5. Сравнение результатов физических и численных экспериментов Для испытаний были изготовлены тест-образцы с жестким заполнителем из древесины по всей длине образца (рис. 12). Общая длина тест-образца составила l = 270 мм, ширина b = 25 мм, толщина заполнителя h = 5 мм, номинальная толщина несущих слоев t = 1 мм, расстояние между опорными роликами составило 250 мм, а между нагружающими роликами - 70 мм. Рис. 12. Общий вид тест-образца в приспособлении Fig. 12. Set-up of the specimen in test fixture Диаграммы нагружения трех тест-образцов с жестким заполнителем приведены на рис. 13. w, мм Рис. 13. Диаграммы нагружения тест-образцов с жестким заполнителем («нагрузка - прогиб») Fig. 13. Loading curves of the specimens with a rigid core (load-deflection) Результаты экспериментов сведены в таблицу. В ней изгибные напряжения, определяемые по классической теории изгиба балок, обозначены через , а касательные напряжения, определяемые по формуле Журавского, - через . Результаты экспериментов на четыреточечный изигиб Four-point bending experimental results Номер образца , кН , МПа , МПа 1 -1,043 292,436 4,230 2 -0,943 269,012 3,824 3 -0,999 284,847 4,049 Среднее -0,995 282,098 4,034 Рис. 14. Уровень напряженного состояния правой половины трехслойного тест-образца по критерию Цая-Ву [42] Fig. 14. Level of the stress state of right half of the sandwich specimen using Tsai-Wu test [42] На рис. 14 приведена картина распределения функции f (показателя уровня напряженного состояния по критерию Цая-Ву) правой половины трехслойного тест-образца в предельном состоянии (толщина несущих слоев равна 1 мм). Следует отметить, что, в отличие от тест-образца с толщиной несущих слоев 3 мм (см. рис. 9), где разрушение происходит вследствие потери устойчивости несущего слоя по поперечно-сдвиговой форме, тест-образец с толщиной несущих слоев 1 мм разрушается вследствие потери прочности. В эксперименте среднее значение предельной нагрузки получилось равным 0,995 кН, а численный расчет предельной нагрузки приводит к значению 0,880 кН, при этом отличие от результата экспериментов составляет 11,6 %. Можно считать это неплохим согласованием расчетных и экспериментальных результатов. Заключение Проведено численное моделирование процесса четырехточечного изгиба трехслойного тест-образца. Зависимость поперечных касательных напряжений от деформаций поперечных сдвигов, формирующихся в несущих слоях, для простоты исследований принята нелинейно-упругой и получена в результате обработки данных испытаний на сдвиг материала несущих слоев. Численное решение задачи проведено в физически и геометрически нелинейной постановке с использованием метода конечных элементов. Сравнительный анализ использования различных кинематических соотношений и разных методов решения нелинейной задачи показал отсутствие существенных различий в результатах решений. Для расчетов применены два подхода, используемые при решении геометрически нелинейных задач (статический и метод продолжения по дуге). Установлено, что при некоторой нагрузке происходит потеря устойчивости несущего слоя по неклассической поперечно-сдвиговой форме в области приложения локализованной нагрузки, которая имитирует воздействие нагружающего ролика. Такой процесс удается обнаружить лишь визуализацией процесса деформирования. Использование статического метода и метода дуги приводит к тому, что итерационный алгоритм расходится при некоторых значениях параметров процесса, хотя зависимость перемещения ролика от нагрузки до этого момента остается практически линейной. Это объясняется тем, что в результате потери устойчивости несущего слоя по поперечно-сдвиговой форме конечные элементы в области приложения нагрузки вырождаются (их площадь устремляется к нулю). Показано, что классическая картина с предельной точкой, отражающая связь статических и кинематических характеристик конструкции, все же может быть обнаружена визуально при анализе зависимости деформаций поперечного сдвига внешнего слоя от нагрузки. Для определения предельной нагрузки, при которой происходит потеря прочности элементов несущего слоя, использован критерий Цая-Ву. Показано, что увеличение диаметра ролика после его некоторого значения мало влияет на предельную нагрузку, в то время как нагрузка, при которой происходит потеря устойчивости несущего слоя тест-образца, весьма чувствительна к изменению его величины. Проведен сравнительный анализ поведения несущего слоя при различных ее толщинах. Делается вывод о том, что при испытании тест-образцов на четырехточечный изгиб с целью определения пределов прочности сжатого несущего слоя необходимо использовать либо тонкие несущие слои, либо ролики с большим диаметром. Приведены результаты экспериментальных исследований трехслойных тест-образцов и дано сравнение экспериментальных значений предельных нагрузок с расчетными, показавшее их неплохое соответствие.

About the authors

V. N Paimushin

Kazan National Research Technical University named after A.N.Tupolev; Kazan (Volga Region) Federal University

R. A Kayumov

Kazan National Research Technical University named after A.N.Tupolev; Kazan State University of Architecture and Engineering

F. R Shakirzyanov

Kazan National Research Technical University named after A.N.Tupolev; Kazan State University of Architecture and Engineering

S. A Kholmogorov

Kazan National Research Technical University named after A.N.Tupolev

References

  1. Petras A., Sutcliffe M.P.F. Failure mode maps for honeycomb sandwich panels // Composite Structures. - 2019. - Vol. 44. - Р. 237-252.
  2. Rupp P., Elsner P., Weidenmann Kay A. Failure mode maps for four-point-bending of hybrid sandwich structures with carbon fiber reinforced plastic face sheets and aluminum foam cores manufactured by a polyurethane spraying process // Journal of Sandwich Structures and Materials. - 2019. - Vol. 21, no. 8. - P. 2654-2679.
  3. Shi H., Liu W., Fang H. Damage characteristics analysis of GFRP-Balsa sandwich beams under four-point fatigue bending // Composites Part A: Applied Science and Manufacturing. - 2018. - Vol. 109. - P. 564-577.
  4. Experimental and analytical study of nonlinear bending response of sandwich beams / V.S. Sokolinsky, H. Shen, L. Vaikhanski, S.R. Nutt // Composite Structures. - 2003. - Vol. 60 - P. 219-229.
  5. Banghai J., Zhibin L., Fangyun L. Failure mechanisms of sandwich beams subjected to three-point bending // Composite Structures. - 2015. -Vol.133 - pp.739-745.
  6. An investigation of the flexural properties of balsa and polymer foam core sandwich structures: influence of core type and contour finishing options / A. Fathi, F. Woff-Fabris, V. Altstadt, R. Gatzi // Journal of Sandwich Structures and Materials. - 2013. - Vol. 15, no. 5. - P. 487-508.
  7. Crupi V., Epasto G., Guglielmino E. Comparison of aluminium sandwiches for lightweight ship structures: honeycomb vs. foam // Marine Structures. - 2013. - Vol. 30. - P. 74-96.
  8. Shi H., Liu W., Fang H. Damage characteristics analysis of GFRP-Balsa sandwich beams under four-point fatigue bending // Composites: Part A. - 2018. - Vol. 109. - P. 564-577.
  9. Viscoelastic behaviour investigation and new developed laboratory slamming test on foam core sandwich / F. Alila, J. Fajoui, R. Gerard, P. Casari, M. Kchaou, F. Jacquemin // Journal of Sandwich Structures and Materials. - 2020. - Vol. 22, no. 6. - P. 2049-2074. doi: 10.1177/1099636218792729.
  10. Piovar S., Kormanikova E. Sandwich beam in four-point bending test: experiment and numerical models // Advanced Materials Research. - 2014. - Vol. 969. - P. 316-319.
  11. Russo A., Zuccarello B. Experimental and numerical evaluation of the mechanical behaviour of GFRP sandwich panels // Composite Structures. - 2007. - Vol. 81. - P. 575-586.
  12. Geometrically nonlinear strain and buckling analysis of sandwich plates and shells reinforced on their edge / V.N. Paimushin, M.V. Makarov, I.B. Badriev, S.A. Kholmogorov // Shell Structures: Theory and Applications. Proceedings of the 11th international Conference. - Gdansk, 2018. - P. 267-270.
  13. Badriev I.B., Makarov M.V., Paimushin V.N. Solvability of physically and geomentrically nonlinear problem of the theory of sandwich plates with transversally-soft core // Russian Mathematics. - 2015. - Vol. 59, no. 10. - P. 57-60.
  14. Лукасевич С. Локальные нагрузки в пластинах и оболочках. Перевод с англ. и польского Б.Н. Ушакова. - М.: Мир, 1982. - 542 с.
  15. Thomsen O.T. Theoretical and experimental investigation local bending effects in sandwich plates // Composite Structures. - 1995. - Vol. 30 - P. 85-101.
  16. Ply drop-off effects in CFRP/honeycomb sandwich panels - experimental results / O.T. Thomsen, W. Rits, D.C.G. Eaton, O. Dupont, P. Queekers // Composites Science and Technology. - 1996. - Vol. 56. - P. 423-431.
  17. Mujikaa F., Pujanab J., Olave M. On the determination of out-of-plane elastic properties of honeycomb sandwich panels // Polymer Testing. - 2011. - Vol. 30. - P. 222-228.
  18. Паймушин В.Н., Шалашилин В.И. Непротиворечивый вариант теории деформаций сплошных сред в квадратичном приближении // Докл. РАН. - 2004. - Т. 396, № 4. - С. 492-495.
  19. Паймушин В.Н., Шалашилин В.И. О соотношениях теории деформаций в квадратичном приближении и проблемы построения уточненных вариантов геометрически нелинейной теории слоистых элементов конструкций // Прикладная математика и механика. - 2005. - Т. 69, № 5. - С. 864-882.
  20. Паймушин В.Н., Шалашилин В.И. О геометрически нелинейных уравнениях теории безмоментных оболочек с приложениями к задачам о неклассических формах потери устойчивости цилиндра // Прикладная математика и механика. - 2006. - Т. 70, № 1. - С. 100-110.
  21. Шклярчук Ф.Н. К расчету деформированного состояния и устойчивости геометрически нелинейных упругих систем // Изв. РАН, Мех. тв. тела. - 1998. - № 1. - С. 140-146.
  22. Работнов Ю.Н., Шестериков С.А. Устойчивость стержней и пластинок в состоянии ползучести // Приклад. мат-ка и механика. - 1957. - Т. 21, № 3. - С. 406-412.
  23. Шестериков С.А. О критерии устойчивости при ползучести // Прикладная математика и механика. - 1959. - Т. 23, № 6. - С. 1101-1106.
  24. Кузнецов А.П., Куршин Л.М. Решение некоторых задач устойчивости пластин и оболочек в условиях ползучести по теории упрочнения // Прикладная механика и техническая физика. - 1960. - № 4. - С. 84-89.
  25. Шестериков С.А. Выпучивание при ползучести с учетом мгновенных пластических деформаций // Прикладная механика и техническая физика. - 1963. - № 2. - С. 124-129.
  26. Терегулов И.Г. Устойчивость пластин и оболочек при неустановившейся ползучести // Исследования по теории пластин и оболочек. - Вып. 3. - Казань: Изд-во Казан. гос. ун-та, 1965. - С. 237-243.
  27. Pian T.H.H. Creep buckling of curved beam under lateral loading // Proc. of 3-rd U.S. Math. Congr. of Appl Mech. - N.Y., 1958. - Р. 58-63.
  28. Клюшников В.Д. Математическая теория пластичности. - М.: Изд-во Моск. ун-та, 1979. - 208 с.
  29. Федоренко А.Н., Федулов Б.Н., Ломакин Е.В. Задача потери устойчивости тонкостенных конструкций из композиционных материалов, свойства которых зависят от типа нагружения // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2019. - № 3. - С. 104-111.
  30. Образцов И.Ф., Васильев В.В. Нелинейные феноменологические модели деформирования волокнистых композитных материалов // Механика композитных материалов. - 1982. - № 3. - С. 390-393.
  31. Dumansky A.M., Liu Hao Analysis of anisotropy of time-dependent and nonlinear properties of unidirectional CFRP // Advances in Composite Science and Technology. IOP Conference Series: Materials Science and Engineering. - 2019. - Vol. 683. Р. 012093. doi: 10.1088/1757-899X/683/1/012093
  32. Паймушин В.Н., Холмогоров С.А., Каюмов Р.А. Экспериментальные исследования механизмов формирования остаточных деформаций волокнистых композитов слоистой структуры при циклическом нагружении // Ученые записки Казан. ун-та. Сер. Физ.-матем. науки. - 2017. - Т. 159, № 4. - С. 395-428.
  33. Новожилов В.В. Теория упругости. - Л.: СУДПРОМГИЗ, 1958. - 411 с.
  34. Васидзу К. Вариационные методы в теории упругости и пластичности. - М.: Мир, 1987. - 542 с.
  35. Bonet J., Wood D. Nonlinear Continuum Mechanics for Finite Element Analysis. - Cambridge University Press, 2008. - 318 p.
  36. Zienkiewicz O.C., Taylor R.L. The Finite Element Method. Fifth edition. Vol. 2: Solid Mechanics 2000. - 316 p.
  37. Голованов А.И., Бережной Д.В. Метод конечных элементов в механике деформируемых твердых тел. - Казань: ДАС, 2001. - 301 с.
  38. Riks E. An incremental approach to the solution of snapping and buckling problems // Solid Structures. - 1979. - Vol. 15. - Р. 529-551.
  39. Crisfield M.A. Non-linear finite element analysis of solids and structures. - Wiley, 1991. - Vol. 1. - 362 p.
  40. Crisfield M.A. A fast incremental/iterative solution procedure that handles “snap-through” // Computers & Structures. - 1981. - Vol. 13. - Р. 55-62.
  41. Паймушин В.Н., Холмогоров С.А. Физико-механические характеристики волокнистого композита на основе углеленты ЭЛУР-П и связующего ХТ-118 // Механика композиционных материалов. - 2018. - Т. 54. - С. 5-20.
  42. Гриневич А.В., Яковлев Н.О., Славин А.В. Критерии разрушения полимерных композиционных материалов (обзор) // Электронный журнал «Труды ВИАМ». - 2019. - № 7. - URL: http://viam-works.ru/ru/articles?art_id=1442.
  43. Тарнопольский Ю.М., Кинцис Т.Я. Методы статических испытаний армированных пластиков. - 2-е изд., перераб. и доп. - М.: Химия, 1975. - 262 с.
  44. Narayanaswami R., Adelman H. M. Evaluation of the Tensor Polynomial and Hoffman Strength Theories for Composite Materials // Journal of Composite Materials. - 1977. - Vol. 11. - Р. 366-377.

Statistics

Views

Abstract - 319

PDF (Russian) - 216

Cited-By


PlumX


Copyright (c) 2020 Paimushin V.N., Kayumov R.A., Shakirzyanov F.R., Kholmogorov S.A.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies