Full Text

Введение Интенсивное развитие современной вычислительной техники и численных методов позволяют использовать для оценки прочности конструкций при малоцикловых термосиловых нагружениях подходы, основанные на непосредственном моделировании процессов деформирования и накопления повреждений в материале конструкций с использованием в качестве определяющих уравнений соотношений механики поврежденной среды (МПС) [1-11]. Необходимость применения таких подходов диктуется современными представлениями о явлении разрушения как процессе, определяемом многостадийным развитием повреждений в материале в процессе его деформирования. Согласно таким представлениям процесс разрушения происходит в результате реализации основных стадий: - первой стадии, характеризуемой зарождением и развитием микропор и микродефектов до образования макроскопического дефекта; - второй стадии, в пределах которой происходит устойчивое развитие образовавшихся макродефектов до критического значения, соответствующего нарушению сплошности материала в рассматриваемой точке; - третьей стадии, характеризуемой слиянием образовавшихся макродефектов в одну или несколько магистральных трещин, приводящих к разделению на части (фрагментации) рассматриваемого элемента конструкции. Соотношение между длительностью этих стадий изменяется в широких пределах в зависимости от физико-химической природы, структуры и условий нагружения [12]. Согласно сложившимся традициям исследование процессов на первых двух стадиях осуществляется на основе различных моделей термопластичности и накопления повреждений, позволяющих описать развитие процесса повреждений в начально-бездефектном материале, а на третьей стадии с использованием методов механики разрушения, предполагающих наличие начального дефекта в виде трещины [13, 14]. Однако трещина зарождается и развивается в объемах материала, подготовленных к разрушению на первой стадии. Поэтому моделирование реальных процессов разрушения конструкций должно обеспечивать возможность описания протекания этих процессов как на начальной, так и на заключительной стадиях. С этой целью могут быть использованы подходы, основанные на применении соотношений МПС, позволяющие прогнозировать развитие дефектов в материале от момента их зарождения до предельного раскрытия магистральной трещины. Реализация таких подходов при решении практических задач требует наличия надежных моделей, описывающих процессы упругопластического деформирования и накопления повреждений материалов на основе различных механизмов разрушения, эффективных методик решения нелинейных краевых задач, сформулированных на основе соотношений МПС, и разработки специальных алгоритмов, позволяющих моделировать процесс развития повреждений от исходного неповрежденного состояния до предельного раскрытия магистральной трещины. Для описания процессов пластического деформирования материалов могут быть использованы любые модели термопластичности, позволяющие описать основные закономерности пластического деформирования материалов при циклических нагружениях [9-11,15-20] а также кинетические уравнения разрушения, позволяющие учесть двухстадийный характер развития повреждений [9-11]. В настоящей статье рассмотрены алгоритмы, позволяющие в приемлемое для инженерных расчетов время проводить моделирование процессов деформирования и разрушения в элементах конструкций при малоцикловых термосиловых нагружениях в рамках соотношений механики поврежденной среды. Алгоритмы реализованы в составе программных средств вычислительного комплекса ВК УПАКС [21] решения на основе МКЭ задач деформирования и разрушения конструкций при квазистатических термосиловых нагружениях. 1. Методические основы моделирования процессов накопления повреждений при малоцикловых термосиловых нагружениях Исследование процессов упруговязкопластического деформирования и разрушения материалов конструкций осуществляется в ВК УПАКС на основе составной иерархической модели поврежденного материала в рамках соотношений МПС [9, 11]. В основу модели положена возможность представления сложного процесса развития взаимосвязанных эффектов деформирования и разрушения в виде последовательности формально независимых элементарных актов, описываемых соответствующими частными моделями пластичности и накопления повреждений. Учет взаимодействия и взаимного влияния таких элементарных актов при описании реальных процессов осуществляется в составной модели, обеспечивающей последовательную инициализацию частных моделей и коррекцию входящих в них параметров (напряжений, поврежденности, параметров, характеризующих историю упруговязкопластического деформирования материала). При этом описание взаимодействия различных видов поврежденности и влияния их на процесс деформирования строится на основе инвариантной по отношению к природе этих повреждений скалярной меры поврежденности ω впервые введенной Л.М. Качановым и Ю.Н. Работновым [1, 2]. Величина ω меняется от значения для неповрежденного материала до - для полностью разрушенного материала. Непосредственное влияние поврежденности на процесс деформирования учитывается в уравнениях равновесия путем введения зависимости упругих характеристик материала от текущего значения функции . В связи с этим при формулировке составной модели поврежденного материала в рассмотрение введены две характеристики напряжений: эффективные , действующие на поврежденных площадках, и приведенные , статически эквивалентные первым, но отнесенные к неповрежденным площадкам. Первые фигурируют во всех частных моделях, определяющих состояние материала в точке тела, вторые используются на уровне описания конструкции при формулировке уравнений равновесия и статических граничных условий. Общие соотношения модели поврежденного материала, устанавливают связь между изменениями приведенных напряжений и деформаций на элементарном шаге изменения внешних воздействий, а также параметрами, характеризующими текущее состояние материала на элементарном шаге изменения внешних воздействий. Они имеют вид уравнений упругости с дополнительными членами , обусловленными эффектами температуры, необратимого деформирования и деградации свойств материала, связанной с накоплением повреждений. (1) , где - модули объемной и сдвиговой деформации неповрежденного материала, отнесенные к уровню температур в исходном (в начале шага) и текущем (в конце шага) состояниях; - значения коэффициентов температурного расширения; - значения мер поврежденности в исходном и текущем состояниях; - значения девиаторных и шаровой составляющих тензора эффективных напряжений в исходном состоянии. Фигурирующие в общей модели величины изменения пластических деформаций и деформаций ползучести , описываемые соответствующими частными моделями, однозначно определяются уровнями эффективных напряжений и температур T в исходном и текущем состояниях, а также наборами скалярных и тензорных параметров, являющихся функционалами процесса и характеризующих истории упругопластического деформирования и ползучести соответственно. Определяющие соотношения в этих моделях записываются для неповрежденного материала и формально не включают в себя какой-либо зависимости от текущей поврежденности. Влияние последней проявляется лишь через уровень зависящих от эффективных напряжений . Исследование поведения конструкций на основе рассмотренного подхода строится путем пошагового интегрирования инкрементальных уравнений, описывающих их равновесное состояние в метрике текущей деформированной конфигурации с использованием специальной комбинированной шаговой схемы [9, 11]. Численное решение линеаризованных задач осуществляется МКЭ на основе изопараметрических квадратичных КЭ с сирендиповой аппроксимацией поля перемещений [11]. Вся физическая информация в элементе (температура, напряжения, деформации и т.д.) определяется на основе квадратичной интерполяции соответствующих параметров в «физических» узлах элемента. Для исследования процессов деформирования и разрушения конструкций от неповрежденного состояния до предельного раскрытия магистральной трещины, при котором происходит исчерпание конструкцией ее несущей способности, на основе рассмотренных методических положений может быть непосредственно использован описанный в [9, 11] общий алгоритм решения нелинейных задач деформирования и разрушения конструкций, заключающийся в пошаговом интегрировании определяющих уравнений для последовательности изменения внешних воздействий, соответствующих текущим этапам нагружения. При этом на каждом этапе решается нелинейная краевая задача для соответствующего изменения внешних воздействий с учетом всей предшествующей истории необратимого деформирования и накопленной поврежденности материала. Дальнейшее преобразование нелинейных уравнений на основе известных схем линеаризации приводит к необходимости решать на каждом шаге целую последовательность линейных краевых задач. Непосредственное применение описанного алгоритма для исследования поведения конструкций для достаточно протяженных траекторий нагружения (таких как малоцикловые нагружения) представляется весьма затруднительным, как из-за большой трудоемкости вычислений, связанной с необходимостью прослеживать решение задач для большого числа шагов по нагрузке, так и из-за возможности накопления численной погрешности. Кроме этого, применение его для моделирования третьей стадии процессов разрушения (анализа траекторий развития трещин) также требует учета ряда дополнительных обстоятельств и введения в общий алгоритм решения задачи некоторых изменений. Для моделирования процесса пластического деформирования материала в предлагаемом алгоритме использован реализованный в ВК УПАКС вариант модели термопластичности с комбинированным упрочнением, описывающий основные закономерности пластического деформирования при циклических нагружениях и вариант модели накопления повреждений, учитывающей двухстадийный характер развития этих процессов [9, 11, 21]. При описании процесса пластического деформирования вводится понятие поверхности текучести F, которая представляется в виде гиперсферы в пространстве девиаторов напряжений . Координаты центра поверхности определяются компонентами тензора остаточных микронапряжений . В модели предполагается, что в процессе развития пластических деформаций поверхность текучести меняет свои размеры (изотропное упрочнение, описываемое изменениями радиуса поверхности текучести ) и смещается в пространстве девиаторов напряжений (кинематическое упрочнение, характеризуемое изменением тензора остаточных микронапряжений ): , (2) где - тензор активных напряжений, . (3) Координаты центра поверхности текучести и ее радиус являются функционалами процесса пластического деформирования и зависят от температуры, пластических деформаций и истории их изменения. Для учета зависимости процесса пластического деформирования от циклиеских свойств материала в рассматриваемую модель также вводится дополнительная поверхность «памяти» для тензора остаточных микронапряжений радиусом , (4) и член, учитывающий стремление радиуса поверхности текучести при монотонном циклическом нагружении к некоторому стационарному значению , зависящему не от предшествующей истории пластического деформирования, а от параметров стабильной петли, характеризуемых значением . В результате изменения перечисленных величин , и на шаге нагружения представляются в виде , (5) , (6) где - стационарное значение радиуса поверхности текучести для данного материала; - материальные функции, получаемые на основе соответствующих экспериментальных данных [9]. При моделировании процесса накопления повреждений введено предположение, что разрушение материала в некоторой точке происходит при достижении критического значения энергии пластического разрыхления определяемой работой тензора остаточных микронапряжений на пластических деформациях: . (7) Предполагается также, что зависит от вида напряженно-деформированного состояния (НДС), температуры Т, но инвариантна по отношению к конкретным траекториям нагружения. Таким образом, разрушение представляется как процесс, характеризуемый изменением энергии разрушения W от значения для неповрежденного материала до - полного разрушения материала в рассматриваемой точке. Для идентификации процесса накопления повреждений в модель также вводится функция поврежденности , представляющая собой нормированный аналог энергии пластического разрыхления Изменение функции поврежденности связано с изменением энергии при фиксированном виде НДС, характеризуемым значением параметра (параметра вида НДС) соотношением , (8) где - изменение энергии пластического разрыхления на элементарном шаге изменения нагрузки. Для описания влияния текущего уровня поврежденности на характеристики процесса деформирования на второй стадии накопления повреждений в рассматриваемой модели вводится скалярная мера поврежденности w. Величина w рассматривается как величина, инвариантная по отношению к виду поврежденности, виду НДС и траекториям нагружения. Изменение меры поврежденности , в свою очередь, связано с накопленным значением , а также с изменением эффективной части функции поврежденности соответствующей второй стадии накопления повреждений: , (9) где - функция материала. Учет наличия фаз при описании накопления повреждений осуществляется по аналогии с моделями, описанными в [9], путем введения переменной, определяющей завершение первой фазы. В качестве такой переменной используется величина , определяемая значением упомянутой выше функции поврежденности к концу первой фазы [11], , (10) где - амплитудное значение энергии пластического разрыхления к концу первой фазы, являющееся функцией температуры Т и параметра , характеризующего размах амплитуды пластических деформаций текущего цикла нагружения [9]. При этом изменение эффективной части функции поврежденности представляется в виде [11] (11) Иначе говоря, мера поврежденности в точке тела становится отличной от нуля лишь тогда, когда функция поврежденности достигает некоторого амплитудного значения . 2. Алгоритм прогнозирования процессов накопления повреждений при малоцикловых термосиловых нагружениях В рамках настоящей статьи рассмотрены алгоритмы прогнозирования процессов деформирования и накопления повреждений в элементах конструкций при малоцикловых термосиловых нагружениях, сочетающие возможности детального описания особенностей процессов деформирования и накопления повреждений с минимизацией объема вычислений. В соответствии с рассмотренными выше предположениями о многостадийном характере развития повреждений в материале в процессе его деформирования в пределах первой стадии в процессе упруговязкопластического деформирования происходит преобразование структуры материала, приводящее к возникновению эффектов упрочнения и разупрочнения, а также зарождению рассеянных по объему материала повреждений в виде микропор, микротрещин, характеризуемых функцией поврежденности материала . Влияние накопленных повреждений на упругие характеристики материала на этой стадии, до достижения величиной некоторого амплитудного значения , незначительно, и мера повреждения в пределах этой стадии принимается равной нулю. При циклических нагружениях в пределах этой стадии может происходить стабилизация процесса циклического деформирования, проявляющаяся в стабилизации амплитудных значений напряжений, деформаций, пластических деформаций. На участке процесса от начала стабилизации до конца первой стадии накопления повреждений, соответствующего значению , параметры процесса пластического деформирования практически постоянны (амплитудные значения напряжений, деформаций, пластических деформаций, компонентов тензора ), а интегральные параметры (длина траектории пластической деформации), , характеризующие накопление повреждений, меняются по закону, близкому к линейному. Поэтому на этом участке можно с достаточной степенью достоверности прогнозировать значения перечисленных выше параметров путем их линейной экстраполяции на определенное число циклов вперед, исключив соответствующую часть трудоемкого процесса пошагового интегрирования уравнений исходной задачи. Стабилизация процесса деформирования при малоцикловых нагружениях обычно наступает по истечении первых 10-20 циклов либо не наступает вообще. В последнем случае применение каких-либо процедур прогнозирования процесса теряет смысл. Завершение периода стабилизации процесса циклического деформирования соответствует началу заметного влияния развивающейся поврежденности на характеристики процесса деформирования, т.е. началу второй фазы развития повреждений. При этом в отдельных точках конструкции будет происходить перераспределение действующих напряжений и деформаций и прогноз дальнейшего протекания процесса циклического деформирования может привести к неверным результатам. В связи с этим в качестве конечной точки, характеризующей состояние материала, до которого целесообразно экстраполировать решение, можно выбрать значение функции поврежденности , соответствующее концу первой стадии накопления повреждений. На второй стадии разрушения при в материале происходит коллективное взаимодействие развивающихся микропор и микротрещин, приводящее к значительному влиянию накопленной поврежденности на макроскопические характеристики материала (модуль упругости, скорость звука и т. д.) и соответственному росту меры повреждения . На этой стадии прогнозирование параметров процесса циклического деформирования с использованием схем линейной экстраполяции весьма проблематично, хотя и в этих случаях можно применить прогнозирование на основе нелинейных алгоритмов. Перечисленные обстоятельства послужили основой создания алгоритмов прогнозирования процессов деформирования и накопления повреждений в элементах конструкций при малоцикловых термосиловых нагружениях [22]. В основу таких алгоритмов положена возможность линейной экстраполяции по циклам нагружения параметров, характеризующих упругопластическое состояние и поврежденность материала в узлах конструкции, на основе установления закономерностей изменения этих параметров на различных участках процесса циклического деформирования. В качестве начальной точки экстраполяции выбирается число циклов, соответствующее началу стабилизации процесса циклического деформирования (точки выхода на линейный участок изменения функции , - номер цикла нагружения ). В конструкции такое условие должно выполняться одновременно для всех физических узлов области. При этом допуск на отклонение функции от линейного изменения может быть установлен опытным путем. В качестве конечной точки экстраполяции в рассмотренном алгоритме выбрано состояние, соответствующее наступлению второй стадии накопления повреждений, то есть когда функция достигнет амплитудного значения . При этом для получения предельного состояния материала исследуемого объекта, соответствующего завершению второй стадии процесса (достижению в наиболее нагруженной точке конструкции значения меры поврежденности ), потребуется продолжить процесс пошагового интегрирования исходных уравнений на основе обычной схемы либо применить алгоритм на основе рассмотренной ниже схемы нелинейного прогнозирования. Суть этой схемы заключается в экстраполяции параметров поврежденности материала и на число циклов определяемое заданным изменением меры поврежденности в наиболее нагруженном узле конструкции, с последующим уточнением ее равновесного состояния. Значения параметров для i-го узла конструкции можно получить из (9) и (11): (12) ; . Здесь - приращение функции поврежденности последнего просчитанного цикла. Число циклов экстраполяции определяется как минимальное из в узлах конструкции: . (13) После вычисления величины осуществляется экстраполяция параметров поврежденности в физических узлах конструкции: ; ; . (14) Экстраполяция второго уровня может производиться неоднократно, пока в некотором физическом узле конструкции мера поврежденности не достигнет критического значения, что будет свидетельствовать о появлении макроскопической трещины. При этом для грубой оценки предельного состояния вместо этого алгоритма можно принять и сразу получить приближенную оценку предельного числа циклов. Очевидно, что правомерность применения предложенных алгоритмов прогнозирования и области его применения существенно зависят от возможности реализации указанных выше условий стабилизации параметров циклического деформирования материалов при различных режимах циклических нагружений. Применительно к случаям малоциклового пластического деформирования элементов конструкций в условиях однородного НДС, которое обычно стремятся реализовать при проведении соответствующих экспериментов, можно выделить ряд типовых режимов с характерными для каждого из них закономерностями изменения пластических деформаций. В частности, известны жесткие нагружения с заданной амплитудой изменения деформаций и мягкие нагружения, в которых задаются амплитуды изменения компонент тензора напряжений, причем в пределах каждой из перечисленных групп могут быть реализованы различные варианты конкретных режимов, определяемых различными значениями параметров асимметрии цикла. В связи с этим была проведена серия численных экспериментов, позволившая установить основные закономерности кинетики процессов малоциклового деформирования стали 12Х18Н10Т на основе анализа результатов численного моделирования для некоторых наиболее важных типов нагружения в условиях однородного напряженного состояния. В результате проведенных исследований установлено, что в случаях жестких симметричных и несимметричных нагружений стабилизация петель пластического гистерезиса наступает достаточно быстро после нескольких первых циклов. При этом для симметричных режимов петли располагаются симметрично относительно начала координат, а для несимметричных - смещены вдоль оси деформаций на величину, соответствующую степени асимметрии цикла. Ширина петель не меняется и определяется амплитудными значениями торцевых перемещений. В случаях мягких симметричных нагружений по мере увеличения числа циклов координаты центров петель гистерезиса практически не меняются, при этом последовательно уменьшается ширина петель, асимптотически приближаясь к некоторому стационарному значению, зависящему от амплитудных значений торцевых нагрузок, которое для рассмотренных задач достигалось при 20 циклов. Для перечисленных (жестких и мягких симметричных) вариантов нагружений, с момента стабилизации процесса циклического деформирования, в пределах первой стадии накопления повреждений изменение функции поврежденности за цикл нагружения n практически постоянно. Также было установлено, что при мягких несимметричных нагружениях характерные особенности процесса малоциклового деформирования проявляются в смещении петель гистерезиса вдоль оси деформаций и уменьшении ширины петель с ростом числа циклов. При этом по мере увеличения значения параметра асимметрии цикла ширина петли и размах деформации за цикл монотонно убывают, а смещение петель меняется от нуля для симметричного цикла при ( , - максимальное напряжение в цикле; - минимальное напряжение в цикле) и достигает максимального значения в диапазоне ( - значение параметра асимметрии пульсирующего цикла). При циклическом деформировании реальных конструкций из-за неоднородности НДС в различных точках материала и возможности перераспределения напряжений и деформаций в процессе пластического деформирования редко реализуется какой-либо один тип нагружения, хотя особенности каждого из них могут превалировать в зависимости от режима нагружения самой конструкции. Полученные результаты позволили рекомендовать предложенный выше алгоритм прогнозирования для моделирования процессов малоциклового деформирования для всех рассмотренных типов нагружения в пределах первой стадии накопления повреждений. В пределах второй стадии накопления повреждений развивающаяся поврежденность начинает приводить к деградации исходных свойств материала и нарушению условий стационарности процесса циклического деформирования, поэтому для моделирования процесса на этой стадии приходится использовать традиционную схему пошагового интегрирования уравнений либо применять прогнозирование на основе нелинейных алгоритмов. 3. Моделирование процессов развития магистральных трещин в рамках соотношений механики поврежденной среды Завершение второй стадии накопления повреждений в наиболее нагруженной точке конструкции (предельного состояния материала объекта) еще не приводит к нарушению ее несущей способности (предельного состояния объекта), так как за ней следует третья стадия, характеризуемая слиянием образовавшихся макродефектов в одну или несколько магистральных трещин и продолжающимся сопротивлением конструкции приложенным воздействиям. При этом результаты моделирования поведения конструкции на первых двух стадиях могут служить отправной точкой для моделирования последующего развития трещин. В соответствии с рассмотренными выше модельными представлениями в процессе упругопластического деформирования материала конструкции в отдельных зонах могут зарождаться и развиваться зоны повреждений, интенсивность которых характеризуется мерой поврежденности . Увеличение меры в физическом узле конструкции приводит к снижению упругих характеристик материала (модулей G и K) в этом узле и тем самым к снижению его сопротивляемости. При этом в районе таких узлов происходит перераспределение напряжений по объему материала. К моменту завершения в узле второй стадии накопления повреждений (достижения предельного значения ) материал в районе такого узла полностью перестает сопротивляться дальнейшему деформированию, перераспределяя воспринимаемую ранее нагрузку на ближайшие соседние узлы. В результате в соседних узлах возрастает интенсивность роста поврежденности, и при достижении критического значения меры материал в районе таких узлов также выключается из работы. Вообще говоря, для дальнейшего моделирования процесса необходимо вводить в районе разрушения элементарную трещину путем введения вместо одного разрушенного узла двух несвязанных узлов и продолжить вычислительный процесс с измененной топологией расчетной области. Однако реализация такого подхода сопряжена со значительными трудностями, так как при этом для каждого случая локального нарушения прочности приходилось бы перестраивать сетку КЭ разбиения исследуемой области, а следовательно, менять структуру информационных массивов, осуществлять переинтерполяцию физических величин на новую сеточную область, не нарушая при этом условий равновесного состояния в локальной зоне повреждения и конструкции в целом. В предлагаемом алгоритме для моделирования процесса развития магистральной трещины предложено использовать способ «выключения» узлов при достижении в них критического значения меры . Напряжения в таком узле принимаются равными нулю, а влияние трещины на окружающую зону материала компенсируется дополнительными силами, определяемыми изменением полных деформаций в рассматриваемом узле и не нарушающими равновесного состояния в районе введенной несплошности [11]. Таким образом, процесс последовательного «выключения» соседних узлов в процессе нагружения конструкции позволяет моделировать процесс развития магистральной трещины без изменения исходной топологии дискретной области. При этом параметры, определяющие предельное состояния исследуемой конструкции, могут быть получены на основе реализованных в ВК УПАКС условий исчерпания конструкцией ее несущей способности [9, 11, 21]. В качестве иллюстрации работоспособности предложенного алгоритма для моделирования процессов зарождения и развития магистральных трещин ниже представлены результаты численного моделирования процесса упругопластического разрушения экспериментального образца с концентратором в условиях плоского изгиба, выполненные на основе программных средств ВК УПАКС [21]. Образец, представляющий собой прямоугольную призму длиной l = 65 мм и поперечным сечением b = 20 мм h = 10 мм с поперечным надрезом глубиной = 4 мм, постоянной по всей ширине образца, шарнирно оперт на две жесткие опоры, расположенные друг от друга но расстоянии = 50 мм. Образец выполнен из алюминиевого сплава АК-4 и нагружен вертикальной силой P, равномерно распределенной по его ширине вдоль линии в среднем поперечном сечении образца. В эксперименте нагружение образца осуществлялось за счет перемещения захватов испытательной машины УМЭ-10Т через индентор, установленный по ширине образца в его центральном сечении. При проведении численного исследования реальное нагружение моделировалось вертикальными перемещениями, приложенными через компенсирующую прокладку по узкой полосе на верхней грани образца. Аналогичные компенсирующие прокладки устанавливались также в районе опор для исключения возникновения нежелательных при проведении численных исследований пластических деформаций. Численное моделирование процессов деформирования и разрушения образца осуществлялось в пространственной постановке для симметричной ¼ его части, ограниченной плоскостями симметрии в продольном и поперечном направлениях, с использованием реализованных в ВК УПАКС двадцатиузловых изопараметрическх конечных элементов. Размер шага по нагрузке выбирался на основе проведения предварительных расчетов из условия удовлетворительного описания кривой P-V (V-вертикальное смещение верхней точки центрального сечения образца), характеризующей процесс разрушения образца. В процессе проведения численного моделирования была получена исчерпывающая информация о распределении всех параметров, характеризующих изменение НДС в материале образца в процессе его нагружения (интенсивности напряжений , главных напряжений (i = 1-3), параметра Одквиста , параметра вида НДС П, меры поврежденности ). На рис. 1 показаны картина распространения трещины в поперечном сечении образца для значения нагрузки t = 1,1 и фотография разрушенного в результате эксперимента образца. а б Рис. 1. Фотография разрушенного образца (а) и картина распространения трещины по результатам численного моделирования (б) Fig. 1. Picture of fractured sample (a) and crack growth according to results of numerical modeling (b) На рис. 2 показан график зависимости силы P от перемещения контролируемой точки V, построенный на основе результатов численного моделирования, в виде пунктирной линии и результатов эксперимента, обозначенных сплошной линией. Рис. 2. Зависимости силы P от перемещения контролируемой точки V (сплошная линия - эксперимент, пунктирная линия - численное моделирование) Fig. 2. Dependences between force P and displacement of controlled point V (solid line shows the experiment, dotted line shows numerical modelling) Приведенные результаты позволяют судить о том, что результаты численного моделирования процесса разрушения на основе рассмотренного выше алгоритма хорошо согласуются с экспериментальными данными. Максимальное значение силы P, отличающееся в расчете от экспериментального значения на величину 3 %, достигается в расчете и эксперименте при перемещениях, отличающихся не более чем на 3,6 %. 4. Оценка работоспособности и вычислительной эффективности предложенных алгоритмов моделирования процессов малоциклового разрушения С целью оценки работоспособности и вычислительной эффективности описанного выше алгоритма прогнозирования процессов развития повреждений при малоцикловых нагружениях рассмотрен ряд примеров численного моделирования процессов малоциклового деформирования элементов конструкций, выполненных с помощью программных средств ВК УПАКС. Численное решение задач осуществлялось в осесимметричной постановке с использованием изопараметрических КЭ осесимметричного тела с квадратичным законом распределения функций перемещений. В качестве одного из таких примеров ниже представлены результаты численного исследования процессов малоциклового деформирования и разрушения цилиндрического образца, выполненного из нержавеющей стали 12Х18Н10Т, с кольцевой выточкой, расположенной в среднем по длине образца. Длина рабочей части образца L = 30 мм, диаметр рабочей части Ø = 58 мм, радиус выточки R = 8,5 мм. Необходимые для расчетов материальные функции материала получены на основе экспериментальной информации, приведенной в [9]. Образец, предварительно нагретый до температуры Т = 350 °С, находился в условиях жесткого симметричного циклического нагружения с амплитудой перемещений на торцевой поверхности образца 0,095 мм. Численное моделирование осуществлялось для двух вариантов задачи. В первом варианте исследование процесса производилось до момента образования трещины без использования процедуры экстраполяции, а во втором принималось KCL = 20 (KCL-заданное в условиях задачи минимальное число циклов, начиная с которого допускается экстраполяция). В результате проведенных исследований было установлено, что на первом цикле нагружения при значении = 0,019 мм в наиболее нагруженной точке, расположенной на наружной поверхности центрального сечения образца (в дальнейшем точка B), появились пластические деформации, которые к концу первой четверти цикла при = +0,095 мм (в дальнейшем - максимальной точки цикла) охватили практически все сечение образца. При дальнейшем нагружении зона пластичности распространялась по объему образца и стабилизировалась при 15 циклов. В результате решения первого варианта задачи было установлено, что в точке B первая стадия накопления повреждений завершилась при 740 циклов, причем при дальнейшем увеличении числа циклов процесс повреждений в этой точке развивался наиболее интенсивно. При значении 1040 циклов мера поврежденности в этой точке достигла предельного значения 0,99. Для второго варианта задачи экстраполяция параметров процесса была произведена при n = 20 циклов на n = 750 циклов, после чего при n = 770 циклов в точке B завершилась первая стадия накопления повреждений. Прогнозируемое число циклов до разрушения материала в этой точке оказалось равным n = 1067. Параметры циклического деформирования для различных вариантов расчета Parameters of cyclic deformation for different computation variants i 1 2 NL 5000 20 N0 - 20 NE - 750 N1 740 770 Nr 1040 1067 (800) 0,700 0,673 (800) 30,46 29,79 Ef - 0,721 Er - 2,60 Результаты расчетов представлены в таблице, содержащей наиболее интересные для анализа задачи параметры, подсчитанные в точке B, для рассмотренных вариантов задачи: - номер цикла, начиная с которого произведена экстраполяция; = - число циклов экстраполяции; - номер цикла завершения первой стадии накопления повреждений (появление меры повреждений ); - общее число циклов до разрушения; - значение функции поврежденности для цикла с номером ; - значение параметра для цикла с номером ; - мера эффективности процедуры экстраполяции , где для варианта 1; - мера погрешности прогнозирования момента разрушения (при достижении мерой поврежденности значения = 0,99); Анализ представленных в таблице данных позволяют судить о приемлемой точности и достаточно высокой эффективности прогнозирования на основе предложенного алгоритма процессов накопления повреждений при жестких симметричных нагружениях в условиях неоднородного напряженного состояния. Для рассматриваемой задачи с помощью такого алгоритма удалось сэкономить более 70 % трудоемкости вычислительного процесса, необходимой для определения начала малоциклового разрушения образца. Исследование процесса дальнейшего разрушения образца осуществлялось на основе продолжения решения второго варианта задачи для третьей стадии разрушения с использованием рассмотренного выше алгоритма моделирования развития магистральных трещин. После завершения второй стадии разрушения материала образца в точке B, с увеличением числа циклов поврежденность развивалась в виде магистральной трещины вдоль поперечного сечения средней части образца. Конфигурация трещины для циклов n = 1122, 1181 и 1237, соответствующих различным уровням продвижения трещины, изображена на рис. 3. n = 1122 n = 1181 n = 1237 Рис. 3. Конфигурация трещины на различных циклах нагружения Fig. 3. Crack configuration at different loading cycles Предельное число циклов до полного разрушения образца (потери несущей способности) составило циклов. Изменение состояния материала образца в процессе его малоциклового деформирования сопровождалось перераспределением напряжений и деформаций в образце и, в частности, изменением осевой силы, действующей по его торцам. График изменения отношения осевой силы P, вычисленной по значению приведенных напряжений на торце образца в конце текущей четверти цикла, к ее значению P0 в конце первой четверти первого цикла в зависимости от числа циклов n показан на рис. 4 сплошной линией. На графике можно выделить три характерных участка, отражающих специфику поведения конструкции в условиях упругопластического малоциклового деформирования: - первый участок ( = 1-40) циклов, характеризуемый ростом торцевой силы за счет упрочнения материала в зонах пластичности; - второй участок ( = 40-1050), в пределах которого торцевая сила меняется незначительно; - третий участок ( 1050), характеризуемый заметной деградацией свойств материала, развитием магистральной трещины в центральном сечении образца и снижением торцевой силы. 1,00 1,05 1,10 0 100 1000 1100 1200 1300 P/P0 l/R n 0,2 0,4 0,6 Рис. 4. Отношения силы Р, вычисленной по значению приведенных напряжений на торце образца, в конце текущей четверти цикла, к ее значению P0 в конце первой четверти первого цикла (сплошная линия) и длины трещины l к радиусу выточки R (пунктирная линия) в зависимости от числа циклов n Fig. 4. Ratios between force calculated based on adjusted stresses at the sample end, at the end of the current quarter cycle to its value P0 at the end of the first quarter cycle (solid line) and the crack length l to the radius of the recess R (dashed line) depending on the number of cycles n На рис. 4 пунктирной линией изображен график изменения относительной длины трещины (отношения длины трещины l к радиусу выточки R) от числа циклов n. Заключение В статье представлено описание алгоритмов численного моделирования процессов зарождения и развития усталостных трещин при малоцикловых термосиловых нагружениях в рамках соотношений механики поврежденной среды, основанных на учете многостадийного характера развития повреждений в процессе разрушения материала. Рассмотрен алгоритм прогнозирования процессов деформирования и накопления повреждений в элементах конструкций при малоцикловых термосиловых нагружениях, сочетающий возможности детального описания особенностей процессов деформирования и накопления повреждений на первых стадиях разрушения с минимизацией объема вычислений при численном моделировании названных процессов на основе МКЭ. Для моделирования третьей стадии разрушения предложен алгоритм, позволяющий описать процесс развития трещины в элементах конструкций на основе результатов численного моделирования процесса на первых двух стадиях без изменения исходной топологии КЭ области исследуемых конструкции. Алгоритмы реализованы в составе программных средств ВК УПАКС. С целью обоснования эффективности и области применения предложенных алгоритмов в статье представлены результаты численного решения ряда задач квазистатического и малоциклового разрушения элементов конструкций.

About the authors

V A Gorokhov

Research Institute of Mechanics National Research Lobachevsky State University of Nizhny Novgorod

D A Kazakov

Research Institute of Mechanics National Research Lobachevsky State University of Nizhny Novgorod

S A Kapustin

Research Institute of Mechanics National Research Lobachevsky State University of Nizhny Novgorod

Yu A Churilov

Research Institute of Mechanics National Research Lobachevsky State University of Nizhny Novgorod

References

Statistics

Views

Abstract - 196

PDF (Russian) - 70