Заглавная страница Избранные статьи Случайная статья Познавательные статьи Новые добавления Обратная связь КАТЕГОРИИ: АрхеологияБиология Генетика География Информатика История Логика Маркетинг Математика Менеджмент Механика Педагогика Религия Социология Технологии Физика Философия Финансы Химия Экология ТОП 10 на сайте Приготовление дезинфицирующих растворов различной концентрацииТехника нижней прямой подачи мяча. Франко-прусская война (причины и последствия) Организация работы процедурного кабинета Смысловое и механическое запоминание, их место и роль в усвоении знаний Коммуникативные барьеры и пути их преодоления Обработка изделий медицинского назначения многократного применения Образцы текста публицистического стиля Четыре типа изменения баланса Задачи с ответами для Всероссийской олимпиады по праву Мы поможем в написании ваших работ! ЗНАЕТЕ ЛИ ВЫ?
Влияние общества на человека
Приготовление дезинфицирующих растворов различной концентрации Практические работы по географии для 6 класса Организация работы процедурного кабинета Изменения в неживой природе осенью Уборка процедурного кабинета Сольфеджио. Все правила по сольфеджио Балочные системы. Определение реакций опор и моментов защемления |
Свойства резин в сложном НДС
Возвращаясь к уравнению (2.6.8), воспользуемся формулами (П.2.7) (Приложение 2) для выражения главных условных напряжений f1 и f2 через производные упругого потенциала. В результате из (2.6.8) получится уравнение, которое будет строго выполняться, если аналитическое выражение для U идеально хорошо описывает эксперимент. Практически задача сводится, как и прежде, к минимизации функционала (2.2.4) или (2.2.4-1). Ее решение осуществляли численно методом многопараметрической оптимизации, описанным выше. Испытания проводили для каждого резинового образца девять раз, по числу сменных боковых планок рамки. Каждое следующее испытание проводили через время, достаточное для релаксации напряжения после предыдущего испытания. В результате каждого испытания записывали от 800 до 1200 пар значений l12 и F. Кроме испытаний на одноосное растяжение и с помощью рамки, для получения более полной информации о виде потенциальной поверхности проводили испытания образцов в виде пробки на сжатие. В таблице 2.7.1 приведены полные результаты аппроксимации всех полученных экспериментальных данных теми же аналитическими выражениями, которые использованы для одноосного растяжения (табл. 2.2.1). Предложенные в данной работе потенциалы использовались в инвариантном виде (2.2.12).
Таблица 2.7.1 Результаты аппроксимации экспериментальных данных, полученных в сложном НДС, потенциалами, приведенными в разделе 2.2 и потенциалами (2.2.12). Минимизируемый функционал (2.2.4).
Общее замечание ко всей таблице следующее. В графе «число шагов счета» стоят существенно различающиеся числа. Это связано с тем, что для разных вариантов использовался разный алгоритм. Малое число шагов соответствует «овражному» варианту. Результаты показывают, что оба варианта работают, но с разной скоростью. Примеры – позиции 11 и 12; 28 и 29. Видно, что из однопараметрических потенциалов лучше работает потенциал Хазановича. Такой же вывод был сделан по результатам аппроксимации одноосного растяжения, но там дисперсия была в три раза больше. Можно заключить, что однопараметрические потенциалы лучше работают во всей области деформирования, чем только при одноосном растяжении. Двухпараметрический потенциал Муни- Ривлина совсем ненамного лучше однопараметрического Хазановича. Здесь ситуация, похожая на ту, что наблюдалась для ненаполненных резин [42, 43, 44], где был сделан вывод о предпочтении потенциала Хазановича неогуковскому. Трехпараметрические потенциалы, как и по результатам обработки одноосного растяжения, существенно различаются по точности описания эксперимента. Потенциал Исихары и др. практически не лучше Муни- Ривлина. Причем его константы и точность аппроксимации не зависят ни от точности задания шага поиска (графа «точность определения констант»), ни от исходных значений констант. Среднеквадратическое отклонение порядка 28% слишком велико. О потенциале Муни- Ривлин + Хазанович можно сказать те же слова. Потенциал Черных относится к той же категории, но здесь налицо существенное влияние начального приближения на конечные значения констант, хотя точность на это реагирует слабо, оставаясь в пределах тех же 27% среднеквадратического отклонения. Существенно лучше других потенциал Харт-Смита. Значение S = 16.6% говорит само за себя. Этот потенциал превосходит, как мы увидим ниже, все известные четырехпараметрические потенциалы, кроме одного из них. Кроме того, этот потенциал вполне можно перевести в разряд двухпараметрических, т.к. его третий параметр, являющийся множителем перед слагаемым со вторым инвариантом тензора деформации, во всех случаях близок к нулю. Из известных четырехпараметрических потенциалы Александера и Бидермана по дисперсии заметно лучше потенциала Огдена, хотя по среднеквадратическому отклонению это не так бросается в глаза. Обращает на себя внимание одинаковая точность потенциала Огдена и большинства описанных выше трехпараметрических потенциалов. Но заметно выделяется потенциал Блатца, Шарды, Чоэгла, для которого S = 11%.
Прежде чем перейти к анализу предложенных нами выражений, обратим внимание на следующее. Потенциалы (2.2.12) содержат три параметра. Однако из них можно сделать четырехпараметрические, рассматривая в качестве четвертого параметра показатель n инварианта (2.2.10). Из содержания таблицы 2.7.1 понятно, когда рассматриваются 3, а когда 4 параметра. В случае трех параметров в скобках добавлено число, указывающее значение n. Для этих потенциалов налицо существенное (в 4 раза и более) уменьшение дисперсии в оптимальных вариантах всех четырех выражений (2.2.12) по сравнению с потенциалами Александера и Бидермана. Среднеквадратическое отклонение расчета от эксперимента S для них составляет 11%, т.е. такое же, как и для потенциала Блатца и др. Итак, из известных ранее потенциалов для описания свойств наполненных шинных резин в области малых и средних деформаций в условиях сложного НДС хорошо себя проявили два выражения – Харт-Смита (3 параметра) и Блатца, Шарды, Чоэгла (4 параметра). Второй из них также хорошо сработал и в случае одноосного растяжения, а Харт-Смит там не выделился. Все из предложенных нами выражений показали хорошие результаты в этих двух случаях. На примере логарифмического 1 видно, что для значения n = 6 точность описания эксперимента существенно зависит от точности определения констант (позиции 27-29). Интересен следующий результат. Для всех четырех потенциалов оптимальное значение степени инварианта n близко к 1.5. (Лишь для логарифмического 1 оно составляет 1.7). Можно сделать практический вывод о том, что достаточно использовать потенциалы с тремя параметрами, задав четвертый равным 1.5. На примере экспоненциального потенциала (позиции 43-44 и 46-47) видно, что и здесь точность определения констант существенно влияет на точность результата. Вполне достаточной является величина 10-6. Важная роль константы, ответственной за скорость изменения модуля в области малых деформаций, демонстрируется на примере усеченного логарифмического потенциала (выражение (2.2.5)), лишенного этой константы (позиции 24 - 26) Хотя этот потенциал и ведет себя лучше многих известных ранее, даже четырехпараметрических, он заметно хуже, чем потенциалы (2.2.12). Попутно отметим, что существенно разные начальные приближения дают один и тот же результат.
Теперь попробуем ответить на вопрос: а можно ли строить упругий потенциал только по данным одноосного растяжения? Насколько мы ошибемся, если будем этот потенциал использовать для описания произвольного сложного НДС? Может быть, правы авторы работ [63, 65, 66], которые так и рекомендуют поступать, и незачем городить весь этот огород со сложным НДС? Ответ на этот вопрос дают результаты, приведенные в таблице 2.7.2. В графе «дисперсия» приведены данные для одноосного растяжения и в скобках – точность аппроксимации остального эксперимента (использованы те же экспериментальные данные, что и при построении предыдущей таблицы). Внимательный читатель может нас заподозрить в повторе. Действительно, в таблице 2.2.1 уже были приведены данные по одноосному растяжению. И в таблице 2.7.2 для всех потенциалов, кроме предлагаемых нами, данные те же, что и в упомянутой таблице 2.2.1. Наши же потенциалы другие. Прежде они соответствовали выражениям (2.2.3), теперь – (2.2.12). Теперь они в инвариантном виде.
Таблица 2.7.2. Результаты аппроксимации экспериментальных данных в условиях одноосного растяжения для потенциалов, использованных в табл. 2.2.1. В столбце «дисперсия» в скобках приведены значения дисперсии для сложного НДС. Минимизируемый функционал (2.2.4).
Анализ полученных результатов дает следующее. Выражения (2.2.12) описывают одноосное растяжение с той же точностью, что и (2.2.3). Следовательно, в строгом смысле произвольный способ перехода от одноосного растяжения к инвариантному виду оказался удачным, что мы и предполагали. Это первый существенный вывод. Второй заключается в том, что если прогнозировать поведение резины в сложном НДС на основе данных, полученных в условиях одноосного растяжения, то можно ошибиться в 1.5 - 4 раза по среднеквадратическому отклонению. В следующей таблице 2.7.3 приведены аналогичные результаты, но полученные минимизацией функционала (2.2.4-1).
Таблица 2.7.3. Результаты аппроксимации зависимостей одноосного растяжения выражениями, построенными в инвариантном виде. Минимизируемый функционал (2.2.4-1).
Для функционала (2.2.4-1) лучшие значения дисперсии находятся на том же уровне, что и для функционала (2.2.4), табл. 2.7.2. Однако для малопараметрических потенциалов описание эксперимента более точное. Что касается прогноза свойств по результатам одноосного растяжения, то ситуация, как и прежде, неудовлетворительная. Ошибка в сложном НДС может быть непредсказуемо велика (Исихара и др., позиция 4; Огден, позиция 8). Лучшее предсказание по одноосному растяжению не гарантирует лучшего предсказания в сложном НДС (Блатц и др., позиции 13 и 14). Тот результат, что прогнозировать свойства в сложном НДС рискованно по данным одноосного растяжения, подтверждает выводы из работы [34]. Возникает вопрос: может быть, достаточно добавить к одноосному растяжению одноосное сжатие, и этого будет достаточно? В следующей таблице дана проверка этого предположения.
Таблица 2.7.4 Результаты аппроксимации эксперимента по одноосному растяжению-сжатию. Минимизируется функционал (2.24-1).
Из результатов этой таблицы следует, что, как и прежде, потенциалы (2.2.12) заметно лучше других, но ошибки в сложном НДС слишком велики. Они несколько меньше, чем для случая использования только одноосного растяжения, но все-таки велики. Итак, для повышения точности расчета констант потенциалов требуется проводить испытания в сложном НДС. С помощью предложенной нами рамки это не слишком утомительно. То, что по результатам растяжения-сжатия плохо прогнозировать поведение резины в сложном НДС, было известно и раньше. Это, например, убедительно аргументировано в монографии Уорда [34]. Правда, там вывод сделан для ненаполненных резин. Теперь мы можем распространить его на наполненные резины. Для полноты картины приведем результаты аппроксимации полного эксперимента в сложном НДС для функционала (2.2.4-1).
Таблица 2.7.5. Результаты аппроксимации эксперимента в сложном НДС. Минимизируется функционал (2.24-1).
Эксперимент, по которому строилась таблица 2.7.5, несколько скорректирован. Коррекция носит технический характер и заключается в стыковке данных растяжения и сжатия. Желательность стыковки вызвана наличием некоторого малого значения начального натяжения образца, необходимого для его выпрямления перед испытаниями. Эта процедура тонкая, а так как испытания претендуют на прецизионность, следует оценить влияние этого предварительного натяжения на результаты. Из сравнения результатов таблиц 2.7.1 и 2.6.5 видно, что погрешности лучших потенциалов остались на уровне того, что было в таблице 2.7.1. Изменилось то, что значительно сблизились дисперсии всех потенциалов. Особенно это заметно на однопараметрическом потенциале Хазановича. Но ранжирование потенциалов осталось прежним. Возникает вопрос: какому из функционалов - (2.2.4) или (2.2.4-1) – отдать предпочтение? Какой дает более правильные результаты? Ответ следующий. Те соображения, которые были высказаны в разделе 2.2, подтвердились. Более чувствительным оказался функционал (2.2.4). Этот вывод остается в силе и для скорректированного эксперимента (см. предыдущий абзац). Следующее соображение имеет общий характер. Существует проблема правильного задания исходных значений констант, о чем уже говорилось выше. Часто точность решения существенно зависит от выбора начального приближения. Повторим, что эта проблема фатальная. Не существует строгих теорем о числе решений и месте нахождения оптимального решения для минимизации функционалов вида (2.2.4) или (2.2.4-1). Всегда остается сомнение, а не проглядели мы более глубокий минимум? Указанная проблема особенно актуальна для многопараметрических потенциалов. Для них возможны такие ситуации, когда небольшие колебания значений констант существенно влияют на глубину минимума. Выходов из этой ситуации можно назвать четыре. Первый – считать задачу для большого числа начальных приближений. Второй – не использовать потенциалы с большим числом констант. Третий – строить потенциалы на основе моделей, позволяющих придавать константам определенный физический или геометрический смысл, тогда легче будет оценивать область их наиболее вероятных значений. Четвертый – следить за тем, чтобы при добавлении к выражению для упругого потенциала нового члена дисперсия уменьшалась. Если она увеличивается, то это явное свидетельство того, что глобальный минимум не найден. Примером последнего утверждения может служить последовательность следующих потенциалов: Неогуковский, Муни-Ривлин, Исихара и др., Бидерман. В этом ряду по данным всех таблиц дисперсия уменьшается. Есть еще один, пятый, способ проверки глобальности найденного минимума. Но он может быть применен не ко всем, а только к тем потенциалам, в которых константы входят только в виде множителей. Этот способ состоит в применении к таким потенциалам классического метода наименьших квадратов [215, с. 268], который дает единственное решение. Кроме поиска глобального минимума, этот метод может служить критерием правильности работы всего алгоритма многопараметрической оптимизации (решения этими двумя методами обязаны совпадать). Оптимальное число констант для нашей задачи – 3 или 4, что следует из приведенных данных. Это обеспечивает точность по среднеквадратическому отклонению около ±10%. Здесь следует уточнить, какие отклонения расчета от эксперимента мы считаем приемлемыми, а какие нет, и почему. Почему, например, ±10% по среднеквадратическому отклонению допустимо, а ±20% уже плохо? Все дело в том, для каких целей будет применяться в дальнейшем найденный упругий потенциал. Если на его основе рассчитывать равновесную конфигурацию шины с использованием различных численных схем, то, видимо, точность ±20%, даже ±40% вполне удовлетворит расчетчика. Практически всю нагрузку в шине держит корд, и это позволяет снижать требования к точности задания свойств резины. Если же использовать потенциал для определения наиболее нагруженных (напряженных, максимально деформированных, имеющих максимальную плотность энергии деформации, или нечто другое) мест, то здесь требования к точности расчета НДС существенно возрастают. Почему? Да потому, что нас нагруженность волнует, в первую очередь, с точки зрения числа циклов, которое шина выдержит до разрушения (что соответствует числу пройденных километров). А зависимость числа циклов от величины нагрузки (напряжения, деформации, энергии цикла или чего-то другого) растет как степенная функция с показателем около 7-10. Отсюда ясно, что погрешность в среднеквадратическом отклонении на ±10% или на ±20% даст прогноз по ходимости, который может отличаться на порядок. (Об этом уже говорилось в литературном обзоре). Это и есть та основная причина, которая вынуждает нас среди всех упругих потенциалов тщательно выискивать наилучший.
2.8 Порообразование. Задача Ламе для больших деформаций резины
Явление образования пор в резине в процессе вулканизации хорошо известно технологам-резинщикам [[265]]. Причиной его является наличие в составе резиновой смеси жидкой или газовой фазы, которая находится в растворенном состоянии при высоких давлениях. При определенных соотношениях между температурой, давлением и степенью структурирования вулканизата в резине образуются микропоры, которые при разрастании существенно ухудшают комплекс эксплуатационных показателей готового изделия. С точки зрения физики и механики это явление можно описать следующим образом. Исходя из характеристик газовой фазы, растворенной в резине (чаще всего это вода), можно определить ее растворимость, зависящую от температуры и давления. Зная кинетику вулканизации данной резиновой смеси (зависимость механических характеристик от температурно-временного режима вулканизации при избытке давления), можно рассчитать максимальные деформации или напряжения от воздействия растворенного газа, которые не вызовут возникновения очагов порообразования и их дальнейшего разрастания. Актуальность расчетной оценки критических параметров явления порообразования связана с тем, что визуально не всегда можно определить момент образования микропор, а имеющиеся методические подходы [[266]] или трудоемки, или не всегда позволяют гарантировать их отсутствие. Кроме явления образования и разрастания пор при вулканизации причиной возникновения пор могут быть процессы смешения и переработки резиновых смесей. В этом случае важно, чтобы соотношение давления, температуры и времени вулканизации способствовало зарастанию пор и отсутствию последних в готовом изделии. Задача о раздувании толстостенной сферической или цилиндрической оболочки в механике твердого тела, подчиняющегося закону Гука, носит название задачи Ламе [[267]]. Ее решение для малых деформаций осуществляется стандартным методом из уравнений равновесия. Рассмотрим задачу о раздувании толстостенной оболочки из высокоэластичного несжимаемого материала для случая больших деформаций с учетом его существенной нелинейности. В отличие от классической задачи Ламе решение будем искать не из уравнений равновесия, а с использованием упругого потенциала – зависимости плотности энергии деформации U от инвариантов тензора больших деформаций. Решение будем искать в общем виде для произвольного вида упругого потенциала и произвольных инвариантов. Для сравнения с экспериментом и с результатами классической задачи Ламе, конечно, потребуется явный вид и потенциала, и инвариантов, и их констант. Решим вначале задачу о раздувании тонкой сферической оболочки (рис. 2.8.1) с исходными внутренним диаметром D=2R и толщиной d, причем D>>d. Объем оболочки Vобол определяется выражением (2.8.1) Раздувание приводит к однородному двуосному растяжению оболочки или, другими словами, к ее одноосному сжатию в l1 раз. При раздувании D ®aD, d ® l1d (a > 1). Из равенства объема шаровой оболочки до и после раздувания (2.8.2) найдем l1: (2.8.3)
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Последнее изменение этой страницы: 2017-02-08; просмотров: 579; Нарушение авторского права страницы; Мы поможем в написании вашей работы! infopedia.su Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Обратная связь - 18.188.108.54 (0.054 с.) |