ВВЕДЕНИЕ
В современной технике возрастает значение проблем прочности. Это объясняется увеличением сложности технических изделий, необходимостью повышения качества, эффективности, надежности и долговечности. Проектирование и эксплуатация новой техники невозможно без всестороннего анализа прочности и надежности ее элементов. Для обеспечения заданного ресурса безопасной работы конструкции необходимо знать причины возникновения разрушения материала и конструкции в целом, а также характер развития процесса в зависимости от заданных условий внешнего воздействия, рабочей среды и структуры материала. Поэтому проблема разрушения является основной проблемой учения о сопротивлении материалов внешним воздействиям.
Направление науки о прочности, которое связано с исследованием несущей способности конструкции с учетом начального распределения повреждений и с изучением закономерностей зарождения и развития трещин, получило название механики разрушения. Создание механики разрушения явилось одним из величайших достижений ученых-механиков и материаловедов за последние 60 лет. Достигнутое понимание основных механизмов разрушения, а также внедрение контроля поврежденности в методологию технического проектирования оказало большое влияние на экономический аспект проблемы разрушения конструкций.
Появление трещин в конструкциях и их разрушение, которое происходило при средних расчетных напряжениях ниже предела текучести, показали необходимость дополнить классические методы расчета на прочность дополнительными условиями, которые учитывают развитие трещин, и новыми характеристиками материала, описывающими стадию разрушения.
Впервые задача о критическом состоянии тела с трещиной была решена Гриффитсом [229] с позиций общего энергетического баланса исследуемого объекта. Затем Вестергаард [260] и Снеддон [252] аналитически описали распределение напряжений у вершины трещины в упругом теле. Основываясь на этих результатах, Ирвин [236] предложил в качестве критических величин использовать коэффициенты интенсивности напряжений. Силовой критерий локального разрушения Ирвина состоит в сравнении рассчитанных значений коэффициентов интенсивности напряжений с их критическими значениями. Также Ирвин показал эквивалентность силового критерия разрушения и энергетического подхода Гриффитса в условиях упругой работы материала.
Крупным достижением механики разрушения явилась концепция квазихрупкого разрушения, впервые сформулированная Ирвином [235] и Оро-ваном [244]. Начиная с работ Гриффитса, Орована, Ирвина развитие теории прочности твердых тел пошло по пути изучения процесса разрушения — распространения трещин в твердых телах. При исследовании этого процесса классические подходы в рамках моделей механики сплошной среды оказались непригодными. Гриффите, Орован, Ирвин первыми сформулировали новые критерии прочности, вводя при этом такие характеристики материалов, которые являются определенными инвариантами в моделях механики сплошной среды. Эта теория получила название линейной механики разрушения. Линейная механика разрушения изучает состояние тел с трещинами в предположении, что материал сохраняет свойство линейной упругости вплоть до разрушения во всем объеме тела, за исключением, быть может, небольшой окрестности вершины трещины.
Линейная механика разрушения, созданная Ирвином и Орованом как обобщение теории Гриффитса на случай разрушения металлов, дала импульс для проведения огромного числа работ как теоретических, так и прикладного характера. Она была применена не только для исследования уже разрушившихся конструкций с целью выявить закономерности разрушения, но и интенсивно использовалась для определения скорости роста и усталостных трещин и устойчивого роста трещин коррозионного происхождения. Подробный анализ различных аспектов проблемы разрушения и широкий обзор литературы содержится в работах [156, 203, 209, 171, 181, 51, 196, 25].
Зарождение и рост трещин в элементах конструкций происходит в сложном неоднородном поле напряжений. Поэтому первой и главной задачей при исследовании механизма разрушения конструкции является детальное изучение поля напряжений в окрестности вершины трещины. Анализ напряженного состояния можно провести с помощью трех принципиально различных подходов: аналитическое решение, численный расчет и экспериментальное исследование. Каждое направление имеет известные достоинства и недостатки. Очевидно, что наиболее надежные результаты можно получить при совместном использовании этих методов.
Из приведенного обзора литературы и анализа состояния вопроса о концентрации напряжений в окрестности вершины трещины можно заключить, что к настоящему времени хорошо изученными в теоретической постановке и в смысле получения решений, удобных для практических инженерных расчетов, оказались лишь вопросы о концентрации напряжений в условиях упругости применительно к областям относительно простой фор-
мы. Для многосвязных областей, ограниченных линиями сложного очертания, особенно при наличии угловых точек и усложненных граничных условий, использование аналитических методов становится проблематичным. Даже численные методы при традиционных способах их применения оказываются малоэффективными [122].
Еще более трудную проблему представляет исследование поля напряжений в области пластических, упругопластических деформаций и в условиях ползучести. Точные теоретические решения получены для ограниченного числа случаев. Немногочисленны также и приближенные решения. Большинство этих решений выявляет лишь величину напряжений на контуре концентратора или по опасному сечению, оставляя открытым вопрос о напряженном состоянии по всей окрестности источника концентрации, что весьма существенно для оценки несущей способности элементов конструкций.
Для ряда случаев распределение напряжений в зоне концентрации применительно к реальным элементам конструкции в условиях упругих и упругопластических деформаций исследовано экспериментальным путем. Однако в связи со сложностью экспериментов, многообразием форм образцов и условий их нагружения не удается охватить необходимую для исследований область и сделать достаточно широкие обобщения. Поэтому в последние годы в практике научных исследований и инженерных расчетов в области прочности все чаще прибегают к использованию приближенных численных методов решения задач теории упругости и пластичности.
Численные методы, целесообразность использования которых была замечена уже давно, приобретают особенное значение в настоящее время. С одной стороны, это связано с повышением требований к надежности современных инженерных конструкций, расчет которых должен базироваться на новейших представлениях о поведении материала в различных условиях нагружения. При этом должна быть обеспечена высокая точность расчета. С другой стороны, появление мощных компьютеров привело к переоценке эффективности приближенных методов, связанных со значительными объемами вычислительных операций, что в прошлом ограничивало их применение.
Существенно, что при инженерных расчетах на прочность, как правило, не требуется абсолютная точность. Более того, она обычно лишена смысла, поскольку используемые в расчетах исходные данные о свойствах материала, спектре действующих нагрузок, режимах температуры и т. п. не являются точными величинами. Ясно, что в расчетах должны выдерживаться приближения, учитывающие степень точности исходных данных, и
требования, предъявляемые к конечной точности прочностных расчетов соответствующих конструкций. Все это говорит в пользу численных методов расчета.
Цель диссертации — численное моделирование явления разрушения материалов и конструкций. В диссертации сочетаются как континуальный (механика деформируемого твердого тела), так и дискретный (физика твердого тела) подходы к проблеме разрушения. Численный расчет основывается на аналитических решениях и дополняется экспериментальными исследованиями.
В первой главе рассматриваются численные методы механики твердого деформируемого тела, которые в последующих главах применяются для решения задач механики разрушения.
Подробный обзор и анализ различных подходов к решению задач динамической теории упругости и пластичности можно найти, например, в работах [13, 112,114,113, 147]. Существующие методы решения задач динамики деформируемых твердых тел достаточно условно можно представить в виде трех направлений [114]:
• метод конечных элементов;
• характеристические и сеточно-характеристические методы;
• сеточные или конечно-разностные методы.
Следует сказать, что все три подхода не противопоставляются и их взаимное проникновение все более заметно в последнее время.
Под методом конечных элементов понимают подход, основанный на дискретизации расчетной области и формировании конечных соотношений между искомыми величинами (действующими в узлах силами и их перемещениями) на основе законов механики в вариационной форме, минуя стадию формулировки краевых задач для систем дифференциальных уравнений. Такой подход дает определенные преимущества при описании процесса деформирования тел со сложной геометрией. Метод надежно зарекомендовал себя для решения статических задач и интенсивно используется при исследовании нестационарных процессов в деформируемых твердых телах. Среди отечественных работ этого направления отметим работы [6, 14, 21, 22, 24, 87, 106, 125, 184, 185, 207].
Как сочетание и обобщение методов конечных элементов и вариационно-разностных можно упомянуть дискретно-вариационный метод [109, 110], разработанный для исследования нестационарных процессов в слоистых и композиционных средах. Конечно-элементный подход активно развивается за рубежом. В качестве примера здесь можно указать работы [19, 57, 217, 218, 216, 219, 266].
Детальному изложению, подробному обзору и анализу характеристических и сеточно-характеристических методов посвящена монография [124]. Эти подходы основаны на записи системы дифференциальных уравнений в характеристической форме с последующей их конечно-разностной аппроксимацией. Различают прямой и обратный характеристический метод [123].
Прямой метод состоит в следующем. В начальный момент времени в среде выбирается некоторая сетка, в узлах сетки выстраиваются характеристические поверхности и с помощью соотношений на них определяется решение на некотором удалении от начального момента времени. В случае, когда характеристические поверхности существенно зависят от решения, реализация метода достаточно сложна (определенные трудности вызывает и неединственность характеристической формы системы в многомерном случае) — решение получается в точках, нерегулярно распределенных по пространству и на разном расстоянии от начального момента времени. Однако для решения задач динамики упругих материалов с малыми деформациями, когда скорости распространения возмущений в среде постоянны, метод представляет собой процедуру пересчета решения на один и тот же слой по времени и при этом сохраняет начальную регулярность выбранной сетки. При решении многомерных задач характеристический метод позволяет максимально сблизить области зависимости разностной задачи и исходной системы дифференциальных уравнений. В то же время даже для одномерных линейных задач в случае, когда из узла сетки выходят несколько поверхностей с постоянными, но различными скоростями, требуется интерполяция построенного решения, что расширяет область зависимости разностной задачи и в итоге снижает точность решения.
В обратном характеристическом методе решение для всей области вычисляется в точках, отвечающих одному и тому же шагу по времени. При этом используется конечно-разностная аппроксимация соотношений на характеристических плоскостях, касательных к характеристическим конусам, выходящим из этой точки на предыдущий (нижний) слой по времени. Метод требует интерполирования на предыдущем слое, при этом расширяется область зависимости разностной задачи. В некоторых работах [124, 123] метод называется сеточно-характеристическим.
Описанный подход допускает большое многообразие модификаций, основанных и на различной интерполяции, и на различном выборе шаблона. В результате могут получаться как схемы первого порядка аппроксимации, так и методы второго порядка [77, 111, 173, 182, 186, 207]. Иногда рассматриваются полухарактеристические схемы, которые получаются в результате записи в характеристической форме системы уравнений меньшей раз-
мерности, после предварительной конечно-разностной аппроксимации по одной из пространственных переменных.
Среди работ, посвященных применению сеточно-характеристических методов для решения динамических задач деформирования упругих и упру-гопластических тел, можно указать следующие [78, 80, 79, 172].
Сеточные методы решения нестационарных задач механики деформируемого твердого тела основаны на аппроксимации гиперболической системы дифференциальных уравнений, описывающей движение среды, краевых и начальных условий для нее. В настоящее время это наиболее разработанный и часто используемый способ численного интегрирования задачи. Как правило, алгоритм представляет собой пересчет известного решения с нижнего слоя по времени (начиная с известных начальных условий) на следующий — верхний слой. Однако известны и многослойные методы, когда в вычислении решения на некотором шаге участвуют несколько предыдущих слоев.
В зависимости от того, дает ли процедура такого вычисления непосредственно значения искомых величин на верхнем слое или же для их определения необходимо решать систему перевязанных между собой уравнений, различают явные и неявные схемы. Одним из основных вопросов при выборе разностной схемы решения задачи является вопрос о целесообразности использования явных и неявных схем с точки зрения их точности и экономичности.
В пользу применения неявных схем при решении динамических задач говорит тот факт, что в большинстве своем неявные схемы абсолютно устойчивы, что, в принципе, позволяет вести интегрирование с большим шагом по времени. Кроме того, им может быть отдано предпочтение при решении задач с сильной неоднородностью рассчитываемого процесса, так как использование в этом случае явных схем связано с большим различием величины шага интегрирования в разных точках области, что приводит к необходимости использования малого шага по времени. Несомненно, неявные схемы более экономичны и при вычислении гладких решений.
Однако, как отмечается в [114], при расчете волновых процессов с большими градиентами, на шаг по времени все равно возникают ограничения, вызванные соображениями не устойчивости, а точности, которые не позволяют выбирать его достаточно крупным. В то же время, использование неявных схем с шагами интегрирования, сравнимыми с допускаемыми явными схемами, менее экономично ввиду дополнительных сложностей реализации неявных схем.
Среди используемых неявных разностных схем наибольшее применение получили схемы расщепления [10, 55, 129, 175, 189, 190] и схемы, основанные на методе дробных шагов [45, 46, 48, 81, 82, 83, 84, 85, 130, 213]. Примеры решенных с использованием этих схем задач содержатся в работах [29, 27, 28]. Достаточно эффективными оказываются подходы, допускающие использование комбинированных схем, в частности схем явных в одном направлении и неявных в другом.
Для численного решения двумерных задач динамической теории упругости Г. В. Иванов предложил использовать несколько локальных аппроксимаций неизвестных функций полиномами Лежандра [59, 63]. В отличие от классических конечно-элементных подходов процедура нескольких аппроксимаций для каждой из искомых функций дает возможность сформировать достаточную для обеспечения монотонности численного решения искусственную диссипацию с одновременным расщеплением многомерной задачи на ряд независимых одномерных задач по пространственным направлениям. Содержащиеся в одномерных задачах параметры — константы диссипации — позволяют регулировать в получаемых схемах величину искусственной вязкости. При этом процедура решения каждой из этих одномерных задач является независимой и может быть как явной, так и неявной.
При частном выборе констант диссипации, для случая регулярных сеток, получающаяся явная схема полностью совпадает со схемой С. К. Годунова [42]. В то же время, одно из преимуществ подхода состоит в возможности построения схем, обладающих лучшими диссипативными свойствами по сравнению с методом распада разрыва. Существенно, что при этом не происходит увеличение числа арифметических операций, т. е. помимо улучшения диссипативных характеристик возрастает экономичность схемы.
В первых двух параграфах первой главы диссертации подход Г. В. Иванова, основанный на нескольких локальных аппроксимациях каждой из неизвестных функций, применяется к построению эффективных численных алгоритмов интегрирования двумерных задач динамической теории упругости и моделированию на их основе процессов разрушения твердых тел.
В § 1 на основе нескольких локальных аппроксимаций неизвестных функций линейными полиномами строится схема решения плоской двумерной динамической задачи теории упругости. Алгоритм основывается на расщеплении полной задачи на четыре независимые одномерные, причем каждая из этих одномерных задач содержит максимальный набор констант
10
диссипации. При этом алгоритм приобретает достаточно большой произвол, распорядиться которым можно так, что без дополнительных вычислительных затрат удается добиться более высокой точности при сохранении таких положительных качеств, как монотонность, диссипативность и т. д. Исследована диссипативность, а следовательно, и устойчивость в энергетической норме построенной схемы. Показано, что при лучших диссипа-тивных характеристиках и без увеличения вычислительных затрат, шаг интегрирования по времени в предложенном алгоритме может быть увеличен в л/2 раз по сравнению со схемой распада разрыва. Излагается процедура моделирования процессов хрупкого разрушения. При моделировании процесса хрупкого разрушения в рамках построенных алгоритмов предполагается, что разрушение происходит по границам элементов при достижении нормальными или касательными напряжениями предельных значений. Приводится решение задачи о хрупком разрушении однородного и слоистого колец.
В § 2 рассматриваются два универсальных способа построения неотражающих условий на искусственных границах при численном решении двумерных задач динамики по явным схемам: 1) моделирование, основанное на аналогии с точными условиями для одномерных задач; 2) использование пространственных и пространственно-временных экстраполяции. При решении задач для бесконечных, полубесконечных или протяженных областей возникает необходимость ограничить вычисления в конечной области, поэтому возникает вопрос формулировки граничных условий на границе этой области — так называемых неотражающих условий, обеспечивающих отсутствие всякого влияния на решение извне.
В § 3 подход, основанный на аппроксимации неизвестных функций несколькими полиномами, применяется к решению статических задач теории упругости. Процедура нескольких аппроксимаций напряжений и смещений позволяет строить четырехугольный конечный элемент, для которого условия сопряжения формулируются в виде условий непрерывности усилий и смещений на его гранях. В задачах с сингулярными особенностями в напряженных состояниях применение таких элементов оказывается намного эффективнее, чем использование традиционных конечных элементов. Эти элементы не содержат сингулярных точек, поскольку они содержат только величины, осредненные по граням, и поэтому являются естественными регуляризаторами в задачах с особенностями в напряжениях.
В § 4 рассматривается итерационное решение статических задач механики деформируемого твердого тела методом самоуравновешенных невязок. Каждая итерация начинается с некоторого приближения, которое не
11
удовлетворяет решаемой системе уравнений. Возникающие в уравнениях невязки в статических задачах можно интерпретировать как сосредоточенные силы и моменты, приложенные в узлах конечноэлементной сетки. Целью итерационного процесса является устранение этих сосредоточенных сил и моментов или сведение их в соответствии с некоторым критерием к минимальным значениям. В статических задачах механики деформируемого твердого тела каждая итерация релаксационного метода приводит, как правило, к уменьшению значения положительно определенного квадратичного функционала, что обеспечивает сходимость итераций.
Вторая глава диссертации посвящена численному моделированию процессов разрушения структурно-неоднородных сред.
Основная проблема механики деформируемого твердого тела состоит в установлении связи между внешним воздействием, изменением исходной структуры среды и возникающими вследствие этого механическими полями.
В процессе исследования этой общей проблемы предложены различные модели, описывающие отдельные аспекты поведения деформируемого твердого тела. Каждая конкретная модель характеризуется определенным способом задания состояния системы (упругость, пластичность, вязкоупру-гость и т. д.). Естественно, что выбор модели, отвечающей данному физическому явлению, может быть оправдан только сравнением с экспериментом.
Наиболее интересными с точки зрения механики твердого тела являются модели реальных сред с дефектами. В одной из конкретных моделей (континуальная механика дефектов) сделана попытка описания таких сред. Исходным состоянием данной модели предполагается идеальное упругое тело, задаваемое вектором поля смещений, называемая линейным кристаллом [194, 208, 69].
Присутствие дефектов (дислокаций) приводит к нарушению линейности и, как следствие, к разрушению исходного состояния. Теперь состояние системы характеризуется тензором дисторсии (шесть компонент упругой деформации е и три компоненты вектора поворота со). Следовательно, появление дислокаций неизбежно приводит к реализации вращательных степеней свободы. Разрывность вектора поворота вызывает появление коллективных или ротационных мод деформации. В континуальной механике дефектов рассматриваются три приведенных выше состояния, из которых каждое последующее является обобщением предыдущего [205, 52].
В работах школы академика В. Е. Панина развивается синергетиче-ский подход, рассматривающий деформируемое твердое тело как открытую сильнонеравновесную систему, а пластическое течение — как диссипа-
12
тивный процесс, снижающий уровень напряжений кристалла [170, 161, 169, 163, 164, 157, 158, 165, 162]. В соответствии с этим подходом пластическая деформация твердых тел может протекать только в условиях неоднородного напряженного состояния. Пластический сдвиг зарождается в зонах концентраторов напряжений как локальная потеря устойчивости атомной решетки и проявляется как локальное структурное превращение. Это структурное превращение может распространяться только в поле концентратора напряжений как сугубо релаксационный процесс. Основные концентраторы напряжений возникают на границах раздела и в различного рода зонах стесненной деформации. Их стохастическое распределение в твердом теле определяет так называемую распределенную систему. Взаимодействие релаксационных потоков в распределенной системе обусловливает волновой характер пластического течения твердого тела. При этом особенно существенно, что процесс пластического течения протекает одновременно или последовательно на различных структурных уровнях, масштаб которых определяется геометрией образца и размерами структурных неоднородно-стей, характерных для каждой стадии деформирования, а известные изменения дефектной структуры, сопровождающие пластическую деформацию, — не только следствие протекания процесса, но и причина перехода от одного уровня к другому.
Эти представления получили убедительное теоретическое и экспериментальное подтверждение [167, 168, 159], и в настоящее время волновой характер пластического течения твердых тел не вызывает сомнений.
В условиях кристаллографической направленности сдвига сохранение сплошности материала требует одновременного участия в скольжении не менее пяти систем плоскостей [174]. Однако результаты многочисленных экспериментальных исследований приводят к заключению, что в реальных условиях такая схема не реализуется: число действующих систем, как правило, меньше пяти, а во многих случаях скольжение происходит преимущественно по одной — трем системам скольжения, сопровождаясь эффектами поворота элементов структуры материала [183, 163, 166, 159]. Поворотные моды деформации в настоящее время являются объектами тщательного исследования.
Обнаружено, что повороты структурных элементов деформации наблюдаются с самого начала пластического течения [170]. Теоретическое и экспериментальное изучение этих эффектов привело к формулировке концепции структурных уровней деформации твердых тел [163]. Согласно этой концепции, при заданных граничных условиях любые пластические сдвиги в пределах определенного структурного элемента деформации вслед-
13
ствие их неизотропности обусловливают появление материального поворота, который вызывает действие на структурный элемент деформации со стороны окружения поля поворотных моментов. Как следствие, возникают поворотные моды деформации, вовлекающие в движение всю иерархию структурных уровней деформации материала. Это кардинально меняет механику пластического формоизменения и разрушения твердого тела. Учет данного обстоятельства позволил успешно описать кривые пластического течения в различных условиях нагружения.
Физика волнового характера пластического течения связана с особенностями вовлечения в деформацию множественного скольжения. Ведущим механизмом деформации является первичное скольжение под действием максимальных касательных напряжений, которое всегда порождает первичный материальный поворот. Все остальные механизмы деформации являются аккомодационными поворотными модами, обеспечивающими релаксацию поля поворотных моментов, действующих на структурный элемент деформации со стороны окружающего материала. Аккомодационные механизмы деформации осуществляются вторичными потоками дефектов и могут обусловливать как материальный поворот структурного элемента деформации (множественное скольжение), так и кристаллографический (зернограничное скольжение, миграция границ зерен, фрагментация и др.). Учитывая, что множественное скольжение зарождается только на границах раздела, включая боковую поверхность образца, корректное описание пластической деформации твердого тела может быть проведено только с учетом границ раздела и зависящих от времени релаксационных потоков деформационных дефектов.
Численное моделирование процессов разрушения структурно-неоднородных сред во второй главе диссертации основывается на представлении сплошной среды в виде жестких блоков, соединенных упругопластически-ми прослойками. Блоки могут смещаться и поворачиваться как твердое целое, прослойки же могут деформироваться упругопластически.
В § 5 формулируется модель плоской деформации прослоек, имеющая целью моделирование всего процесса упругопластического деформирования прослоек от момента возникновения пластических деформаций до достижения предельного состояния при комбинированном действии нагрузок растяжения (сжатия), сдвига и изгиба прослоек в общем случае, когда прослойка может быть криволинейной, а толщина ее переменной. Используется представление прослойки в виде слоя четырехугольных элементов. Строятся уравнения, определяющие зависимость усилий на гранях элементов от средних скоростей граней.
14
В § 6 излагаются результаты численного моделирования волн смещений и локализации деформаций при растяжении полосы из жестких недефор-мируемых блоков с упругопластическими прослойками.
В § 7 рассматривается плоская деформация тонких прослоек между жесткими блоками в случае, когда наряду с упругопластическим деформированием прослойки происходит ее вязкое деформирование (ползучесть). Проводится сопоставление построенной модели с моделью установившейся ползучести прослоек, предложенной Л. М. Качановым, и с моделью упру-гопластического деформирования прослоек без ползучести.
Исследование прочности материалов предполагает определение тех условий, при которых материал утрачивает способность противостоять внешнему нагружению и разрушается (разделяется на отдельные части). В рамках классического подхода к проблеме прочности принято, что разрушение происходит тогда, когда определенная комбинация, включающая в себя напряжения, деформации, температуру и некоторые другие параметры, характеризующие состояние материала и его конкретные свойства, достигает критического значения. Таким образом, в пространстве всех возможных значений указанных параметров существует замкнутая предельная поверхность, которая ограничивает область допустимых (с точки зрения прочности) состояний материала [171, 209].
Обычно окончательному разрушению материала, т. е. разделению его на части, предшествуют глубокие изменения в его структуре на макро-, мезо- и микроуровнях. Накопление этих необратимых изменений на различных стадиях деформирования можно отразить, если в уравнение предельной поверхности дополнительно ввести параметры, характеризующие поврежденность материала. В этом случае при помощи критерия прочности можно описать любое промежуточное состояние материала от начала появления пластических деформаций до полного разрушения и таким образом определить работоспособность конструкции.
Наличие в реальном теле остроконечных концентраторов напряжений, в частности дефектов типа трещин, принципиально усложняет его расчет на прочность. В таких случаях классические подходы механики сплошной среды приводят к некорректным результатам. Суть в том, что условие теоретической прочности недостижимо в практических задачах, а критерий Гриффитса-Ирвина успешно «работает» при расчетах трещин, но приводит к несообразностям в случае угловых вырезов и включений, удлиненных гладких дефектов и т. д. [150].
Необходимость корректировки классических критериев прочности особенно четко видна при одновременном растяжении и сдвиге упругой плос-
15
кости, ослабленной лункой, близкой к прямолинейной трещине [137]. В случае трещины характеристические числа для растяжения и сдвига совпадают, в случае лунки они не равны друг другу. В связи с этим, применяя критерий Ирвина, во втором случае получим несообразность.
Определенный прогресс в проблеме выбора критерия разрушения связан с подходом В. В. Новожилова, который в [148] предложил осреднять напряжения в пределах межатомных расстояний и сравнивать их с теоретической прочностью материала на разрыв. Проверено, что оценка Новожилова, взятая в качестве критерия разрушения, снимает ряд противоречий теории Грифитса-Ирвина и совпадает с ней в простейших случаях [136]. Недостаток критерия Новожилова — в неоднозначности интервала осреднения, поэтому наиболее успешно его можно применять в сравнительных оценках для разной геометрии вырезов в изделиях из одного материала.
В связи с широким внедрением расчетов на прочность и трещиностой-кость в современной технике особую актуальность приобрела проблема разработки достаточно простого критерия хрупкого разрушения. Введенный Новожиловым критерий решает поставленную проблему, однако ряд аспектов остаются невыясненными прежде всего из-за «белых пятен» на межатомном (дискретном) уровне, например, в задаче об угловом вырезе на границе двух сред.
Определенные надежды на решение указанных проблем исследователи связывают с дискретным подходом. Достаточно указать на работы М. Я. Леонова [120], В. В. Панасюка [155, 153], Н. Ф. Морозова [131, 132, 134, 135], Г. Зорского [58], Л. И. Слепяна [197], С. А. Назарова [141]. Дискретные модели, особенно в задачах разрушения, позволили обнаружить ряд эффектов, не улавливаемых континуальными моделями, кроме того, дискретные модели представляются привлекательными в силу моделирования реальной атомной структуры.
В третьей главе диссертации на основе объединения континуального и дискретного подходов формулируются необходимый и достаточный критерии хрупкой и квазихрупкой прочности, причем теоретическая прочность монокристалла определяется из решении задачи о деформировании атомных ячеек. Если рассматривать твердые тела со структурой от микроуровня, соответствующего характерному линейному размеру кристаллической решетки, до макроуровня — размер образца, то методы исследования, опирающиеся только на подходы механики сплошной среды обречены на неудачу, поэтому в механике разрушения необходимо применять такие методы исследования, которые не входили бы в противоречие с общепринятыми взглядами физики и химии твердого тела. По мнению автора наиболее
16
плодотворный подход к решению поставленных задач связан с идеями работ В. В. Новожилова [148, 150, 151], вероятно среди гибридных моделей механики разрушения критерии В. В. Новожилова наиболее перспективны.
В § 8 хорошо развитая техника численного решения нелинейных задач механики деформируемого твердого тела методом конечных элементов применяется к решению нелинейных задач по деформированию атомных решеток, характер деформирования которых близок к характеру деформирования стержневой конструкции. Предлагается конечный элемент атомной пары, матрицы и векторы которого близки к соответствующим матрицам и векторам стержневого элемента. При построении системы алгебраических уравнений для атомной решетки ошибки дискретизации не вводятся, в отличие от построения аналогичной системы в задачах механики сплошной среды. Кроме того, для атомных пар используются потенциальные законы взаимодействия атомов. Поэтому использование метода Ньютона приводит к итерационному процессу, сходящемуся (если этот процесс сходится) к точному решению нелинейной задачи по деформированию атомной решетки.
В § 9 рассматривается плотноупакованный слой атомов, имеющий макротрещину и вакансии. Моделируется поведение атомной структуры в окрестности вершины трещины. Изучается потеря устойчивости трехатомной ячейки в плотноупакованном слое атомов при сложном напряженном состоянии. Выявлена существенная зависимость критических нагрузок потери устойчивости трехатомной ячейки от закрепления одного из атомов этой ячейки (схема закрепления определяется наличием или отсутствием вакансий в окрестности вершины трещины). Кривая теоретической прочности для сложного напряженного состояния определена из решения задачи о деформировании и потери устойчивости трехатомной ячейки. Эта кривая для простейшего случая напоминает огибающую кругов Мора для разных напряженных состояний. Обнаружена угловая точка, соответствующая бифуркации решений, из которой возможны два пути закритическо-го деформирования трехатомной ячейки. Проведено сравнение полученных численных значений теоретических прочностей с существующими экспериментальными наблюдениями и с теоретическими оценками. Предложен дискретно-интегральный критерий прочности для сложного напряженно-деформированного состояния. В отличие от классических критериев предложенный критерий допускает предельный переход по параметру длина трещины и описывает как прочность трещиноватых тел, так и прочность тел без трещин. Предложенный критерий описывает хрупкое и квазихрупкое разрушения и частично описывает пластическое разрушение.
17

Получить работу