ВЫСОКОМОЛЕКУЛЯРНЫЕ СОЕДИНЕНИЯ, Серия А, 2004, том 46, № 11, с. 1875-1890
УДК 541.64:532.135
ВЯЗКОУПРУГОСТЬ АНИЗОТРОПНЫХ ПОЛИМЕРНЫХ СИСТЕМ1
© 2004 г. В. С. Волков, А. К. Терешин, В. Г. Куличихин
Институт нефтехимического синтеза им. A.B. Топчиева Российской академии наук 119991 Москва, Ленинский пр., 29 Поступила в редакцию 05.04.2004 г. Принята в печать 12.07.2004 г.
Представлены результаты общего анализа анизотропной вязкоупругости ЖК-полимеров и полимерных нанокомпозитов в жидком состоянии. Эти системы рассматриваются как анизотропные вязкоупругие жидкости с несимметричными напряжениями. В качестве иллюстрации релаксационной анизотропии теоретически исследовано осциллирующее сдвиговое течение полимерных нема-тиков для различных начальных ориентаций директора. Определено релаксационное уравнение ориентации. Впервые экспериментально измерены анизотропные коэффициенты продольной вязкости и продольных динамических модулей для ЖК-полимеров и полимерных суспензий, содержащих плоские силикатные наночастицы. Полученные результаты согласуются с предсказаниями теории.
ВВЕДЕНИЕ
Анизотропные полимерные системы получили в настоящее время широкое распространение. Прежде всего это касается ЖК-расплавов и растворов полимеров, их смесей с изотропными полимерными системами или друг с другом, сшитых ЖК-эластомеров и текучих предшественников (прекурсоров) полимерных нанокомпозитов в условиях эксфолиации. Несмотря на то, что они принадлежат к общему классу вязкоупругих материалов, специфические особенности их структуры влекут за собой анизотропию всех свойств, включая и реологические.
В этих случаях речь идет о термодинамически обусловленной анизотропии структуры и свойств. Однако существует большой круг гетерофазных полимерных систем, в которых при деформировании появляется анизотропия ориентации и морфологии. Прежде всего это относится к смесям несовместимых полимеров, расплавы которых представляют собой вязкоупругие эмульсии, и капли дисперсной фазы способны вытягиваться при течении вплоть до образования замкнутых колец в
1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (коды проектов 01-03-32111 и 04-03-32152), фонда INTAS (грант 00-278), фанта Президента РФ по поддержке ведущих научных школ (грант НШ-1598.2003.3).
E-mail: [email protected] (Куличихин Валерий Григорьевич).
ротационных приборах. Превращение таких текучих прекурсоров в твердые материалы сопровождается "наследованием" анизотропии структурных и морфологических образований, что приводит к анизотропии морфологии и соответственно механических характеристик, таких как модуль упругости, прочность, разрывное удлинение. Но в этих случаях следует иметь в виду, что речь идет о наведенной механическим полем анизотропии, исчезающей в равновесных условиях.
Аналогичная ситуация складывается и для полимеров вообще, которые вследствие анизодиа-метрии макромолекул могут рассматриваться как потенциально анизотропные системы. Однако только в случае ЖК-упорядочения анизотропия становится внутренним свойством полимерных систем, и реология именно таких объектов будет проанализирована в настоящей работе.
Разнообразие анизотропных полимерных систем и широкий спектр их использования диктуют необходимость разработки общего подхода к анализу их анизотропных вязкоупругих свойств. Простейший пример анизотропных полимеров -ЖК-нематики, которым присуще одно предпочтительное направление ориентации, характеризуемое директором. Самым существенным их отличием от обычных изотропных полимерных материалов является наличие анизотропных времен
1875
5*
релаксации. Мы должны представлять себе, что развивающиеся в материале напряжения релак-сируют с разной скоростью в зависимости от направления, а характер анизотропной релаксации, порождающей анизотропную вязкоупругость, определяется строением макромолекул и в конечном счете их жесткостью и межмолекулярным взаимодействием, т.е. теми факторами, которые влияют на структуру ЖК-фазы.
Таким образом, вязкоупругие характеристики анизотропных полимерных систем зависят от направления измерения и релаксационной анизотропии. Вследствие этого они являются тензорными и определяются (в отличие от изотропных систем) существенно большим числом скалярных базовых параметров. Экспериментальное измерение тензорных реологических свойств полимерных ЖК-систем - сложнейшая проблема. В целом, эта область реологии находится еще на начальной стадии развития. Не говоря уже о вязко-упругих характеристиках, более простое определение зависящих от ориентации среды коэффициентов вязкости Месовича или Лесли даже для низкомолекулярных жидких кристаллов не всегда выполнимо, хотя для таких систем практически идеальная одноосная ориентация может быть достигнута приложением магнитных полей разумной напряженности.
В дальнейшем мы рассмотрим наиболее интересные подходы к измерению вязкоупругих характеристик ЖК-полимеров, которые были развиты за последние 10 лет, однако следует принять как непреложный факт, что направляющую роль здесь пока играют теоретические исследования
[1-4].
Развитие молекулярных теорий вязкоупругос-ти было начато с лиотропных ЖК-полимеров [5-8]. Исследованию молекулярных механизмов анйзотропной вязкоупругости термотропных полимеров уделяется значительно меньшее внимание. Наибольший прогресс в описании анизотропной вязкоупругости термотропных ЖК-полимеров в настоящее время достигнут в рамках феноменологического подхода.
Континуальная теория анизотропной вязкоупругости ЖК-полимеров была инициирована работами [1,4]. В них на основе простых соображений симметрии сформулировано квазилинейное реологическое уравнение состояния с анизот-
ропными временем релаксации и вязкостью. Основная цель указанных работ - описание анизотропной вязкоупругости монодисперсных ЖК-полимеров в области медленной релаксации. Эта теория позволяет также анализировать анизотропную вязкоупругость и ориентационные свойства жидких полимерных нанокомпозитов -расплавов изотропных полимеров, наполненных анизодиаметричными плоскими частицами нано-размерной толщины. Такие системы можно условно рассматривать как лиотропные жидкие кристаллы, ибо при определенной степени наполнения анизотропные наночастицы формируют ЖК-фазу. Для них анизотропная вязкоупругость порождается совокупностью изотропных вязко-упругих свойств матрицы и анизотропной вязко-упругостью структур, образованных наночасти-цами, например слоистых силикатов.
Для анализа анизотропной вязкоупругости не-матических полимеров в разработанном подходе определяющую роль играет молекулярная гибкость ЖК-полимеров и связанные с ней упругие деформации, которые сопровождают их течение. При этом ориентационная упругость, обусловленная неоднородным распределением директора, играет второстепенную роль. Для низкомолекулярных жидких кристаллов, наоборот, развивающиеся при течении обычные обратимые деформации очень малы, поэтому основной тип упругости - ориентационная упругость. Хорошим примером систем, в которых молекулярная гибкость играет существенную роль, являются также сшитые ЖК-эластомеры. На этот факт впервые указал Эе Сеппев при анализе их анизотропной упругости в рамках феноменологической термодинамики [9].
В континуальной теории вязкоупругости ЖК-полимеров [1,4] учитывается эффект релаксационной анизотропии. Релаксационные свойства монодисперсных ЖК-полимеров в области медленной релаксации характеризуются тензором времен релаксации четвертого ранга. Термодинамический подход к реологии анизотропных полимерных систем [10, 11] привел к важному выводу о том, что релаксационная анизотропия ЖК-полимеров и эластомеров может быть естественно описана в рамках термодинамики неравновесных процессов. Таким образом, два разных континуальных подхода приводят к одинаковому результату.
Далее мы ограничимся детальным рассмотрением континуального подхода, который имеет основополагающее значение для исследования симметричной и несимметричной анизотропной вяз-коупругости не только ЖК-полимеров, но и текучих полимерных нанокомпозитов. В силу отсутствия удовлетворительной структурной теории этих полимерных систем в настоящее время феноменологический подход для описания их вяз-коупругих свойств более предпочтителен. Аналогичная ситуация наблюдается в теории низкомолекулярных жидких кристаллов. Наибольший успех в описании вязкой нематодинамики достигнут с помощью континуальных теорий [12-14].
ВЯЗКОУПРУГАЯ НЕМАТОДИНАМИКА
С макроскопической точки зрения ЖК-поли-меры представляют собой анизотропные вязко-упругие жидкости, реологические свойства которых различны в разных направлениях. В работе [1] предложено простейшее определяющее уравнение анизотропной вязкоупругой жидкости, учитывающее специфику релаксационных свойств монодисперсных ЖК-полимеров в области медленной релаксации:
о</ = - + <*</
Оаке - _ (1)
Х1¡ке + <5И ~ Л цк1екЬ
где р - изотропное давление, Су - тензор напряжений, а ву - тензор скоростей деформации. В этой модели тензор напряжений является симметричным. Динамические напряжения д/; описываются релаксационным уравнением первого порядка с тензорными вязкостью и временем релаксации Туь,. Релаксационное определяющее уравнение (1) основано на естественном предположении, что течение ЖК-полимеров, как и любых реальных полимерных жидкостей, сопровождается накоплением упругой деформации. Учет анизотропии упругих свойств ЖК-полимеров естественно приводит к анизотропии времен релаксации. Уравнение (1) представляет собой обобщение уравнения Максвелла на случай анизотропных жидкостей с произвольной симметрией. Оно может быть использовано для анализа анизотропной вязкоупругости ЖК-полимеров в области медленной релаксации. Для несжимаемых поли-
мерных нематиков тензор вязкости и тензор времен релаксации записываются как
_ 1 2 3
~ т1а17ке + + ХЗаре
Согласно выражениям (2), симметричная вяз-коупругость полимерных нематиков в общем случае может быть описана тремя базисными коэффициентами вязкости Г)1 и тремя базисными временами релаксации т,. Они представляют собой масштабные множители для трех независимых
базисных тензоров а^ке О = 1, 2, 3), которые являются нелинейными функциями директора п. Эти базисные тензоры определены в работе [1].
В частном случае, когда тензор времен релаксации Хуке равен нулю, реологические уравнения
(1), (2) сводятся к уравнению трансверсально-изо-тропной жидкости Эриксена [15]. Оно является обобщением реологического закона течения ньютоновских систем на случай одноосно анизотропных вязких жидкостей и может быть представлено в виде
= - + Чщекь (3)
где тензор вязкости Цуке определяется формулой
(2).
Хорошо известно, что с точки зрения характера проявления вязкоупругих свойств наблюдается существенное различие между монодисперсными и полидисперсными гибкоцепными полимерами. Особенность монодисперсных полимеров состоит в том, что переход от конечной зоны к зоне плато на частотной зависимости модуля накопления является резким. Таким образом, вязкоупру-гость линейных монодисперсных полимеров высокой молекулярной массы (с зацеплениями) в области медленной релаксации может быть охарактеризована практически одним временем релаксации. Поэтому естественно предположить, что конечная область анизотропного релаксационного спектра монодисперсных ЖК-полимеров определяется одним тензорным (анизотропным) временем релаксации. Это является главной характерной особенностью континуальной модели (1).
Выбор простейшей инвариантной (яуманнов-ской) производной DCy/Dt для скорости изменения напряжений во времени основывается на опыте исследования квазилинейной динамики монодисперсных гибкоцепных полимеров. Согласно работам Виноградова [16], Файтельсона и Якобсона [17], модель изотропной жидкости Максвелла с яуманновской производной достаточно хорошо описывает вязкоупругие свойства гибкоцепных полимеров с узким ММР в области медленной релаксации. Отметим, что этот результат согласуется с предсказаниями нерептаци-онной молекулярной модели нелинейной вязко-упругости гибкоцепных полимеров [18]. Она моделирует движение полимерных молекул движением линейной цепочки броуновских частиц в вязкоупругой (релаксирующей) среде, образованной макромолекулами окружения. Эта модель является естественным обобщением модели Каргина-Слонимского-Рауза на случай полимеров с зацеплениями.
АНИЗОТРОПНЫЕ ПОЛИМЕРЫ
И НЕСИММЕТРИЧНЫЕ НАПРЯЖЕНИЯ
В классической теории Навье и Стокса [19,20] элементарные частицы вязких жидкостей трактуются как материальные точки, для которых понятие внутреннего вращения не имеет смысла. Эта модель приводит к симметричным напряжениям. Напряжения, возникающие в жидкости при течении, могут быть связаны только со скоростью ее деформации. Если считать частицы жидкости действительно частицами, способными вращаться, можно получить гораздо более общее выражение для тензора напряжений, согласно которому напряжения в жидкости могут возникать и при отсутствии скольжения ее слоев относительно друг друга. Таким образом, в жидкостях с внутренними степенями свободы (имеется скрытый мбмент импульса) могут возникать напряжения, которые связаны с их вращением без деформации. Модель изотропной жидкости с внутренними вращениями предложил Born в 1920 г. [21]. Им впервые было введено понятие вращательной вязкости жидкости. Эта модель получила дальнейшее развитие в работе Сорокина [22].
Как известно, для ЖК-сред также необходимо учитывать вращательные степени свободы молекул. Введение внутренних степеней свободы приводит к возникновению несимметричных сило-
вых и моментных напряжений. При этом отсутствует синхронизация между скоростью вращения элемента среды и скоростью внутреннего вращения (вращение директора), что ведет к дополнительной диссипации энергии и дополнительным напряжениям в среде. Первым экспериментальным свидетельством существования несимметричных напряжений (внутренних вращений) в низкомолекулярных жидких кристаллах является эксперимент Цветкова [23]. Он показал, что вращающееся магнитное поле вызывает движение жидкого кристалла. Для объяснения данного эффекта Лесли ввел в реологическое уравнение Эриксена зависимость тензора напряжений от относительной скорости директора. Эта теория теперь называется теорией Лесли-Эриксена [24, 25]. Она описывает вязкую анизотропную жидкость с внутренними вращениями.
Однако теория Лесли-Эриксена пренебрегает спиновыми степенями свободы, т.е. взаимодействием между поворотами молекул вокруг собственных осей. Такой эффект становится существенным при течении нематиков между параллельными плоскостями в магнитном поле, когда директор перпендикулярен направлениям течения и градиента скорости. В данном случае поток не разворачивает директор, а вращает частицы среды вокруг директора. Модель вязкой немати-ческой жидкости с учетом спиновых эффектов рассматривалась в работах [26, 27].
В работах [28, 29], исходя из простых соображений симметрии, мы предложили модель анизотропной вязкоупругой жидкости с внутренними вращениями для ЖК-полимеров различной симметрии (нематиков, смектиков и холестериков). Введение спиновых степеней свободы и анизотропии релаксационных свойств среды связано с весьма радикальным видоизменением и обобщением теории Лесли-Эриксена. В этом случае полный тензор напряжений определяется соотношением
<5 и + рЬи = о,' + Оу, (4)
где р - изотропное давление; а* - тензор напряжений, связанный с ориентационной упругостью. Динамические несимметричные напряжения выражаются релаксационным уравнением первого порядка
Düke ~
Xijke jyt + Gij ~ 4ijkl£kl P)
(Л¡jke ~ тензор вязкости, %уке - тензор времен релаксации). Несимметричный тензор скоростей деформации еке записывается в виде
| = е + = cg - о/ (6)
Здесь g, оз - симметричная и антисимметричная
части градиента скорости. Внутреннее вращение
i i i характеризуется тензором 0)(; = £ljkoik, где со -
угловая скорость внутреннего вращения, eijk - антисимметричный единичный тензор. Антисимметричная часть тензора напряжений обусловлена отличной от нуля относительной скоростью внутреннего вращения.
Для нематических вязкоупругих жидкостей тензор вязкости и тензор времен релаксации являются трансверсально-изотропными относительно директора п. В случае несжимаемой жидкости нематический тензор вязкости определяется как сумма
Лре = &ijke + aijke О)
изотропного тензора
O-tke = Hl&ik&jc + ^i&U&jk и анизотропного тензора четвертого ранга
aijke = Лз ninjnkne + 1]4nink8ji; +
+ i\5njntSie + Лб щпед]к + л? преЬ1к
Нематический тензор вязкости характеризуется семью скалярными коэффициентами вязкости Г|а. Нематический тензор времен релаксации имеет аналогичную структуру [28].
В анизотропной вязкоупругой жидкости с несимметричными напряжениями, определяемой уравнениями (5), (6), могут возникать напряжения при ее вращении без сдвиговой деформации. При этом видимое вращение жидкости переходит в скрытую форму внутреннего вращения или, наоборот, вихревое движение жидкости создается за счет внутреннего вращения. Для специального класса анизотропных вязкоупругих жидкостей с
симметричными напряжениями реологическое уравнение (5) сводится к более простому определяющему уравнению (1).
Тензор напряжений для несимметричных жидкостей можно представить в виде суммы симметричной d¿j и антисимметричной частей óy. Реологическое уравнение (5) ведет к взаимосвязанной системе релаксационных уравнений для симметричных и антисимметричных напряжений
Ddke Daske ~а
+ XUj](ke)~pf + <*ij
т D®ke _i_ т , -i
Шке)~Ш~ Шке]~оГ 1
Здесь по индексам в круглых скобках проводится симметризация. Квадратные скобки обозначают антисимметризацию по соответствующим индексам.
Рассмотрим частные случаи, которые существенно упрощают рассмотренную модель вязко-упругой жидкости с внутренними вращениями. Если среда не имеет предпочтительного направления и ее релаксационными свойствами можно пренебречь, то уравнение (5) сводится к определяющему уравнению жидкости Борна [21]
0¡; + pS¿i = тк7 + ЛА, (9)
(t|r - вращательная вязкость), которое можно обобщить на вязкоупругие среды. В случае изотропного тензора вязкости и тензора времен релаксации из выражения (5) получаем следующее релаксационные уравнения для вязкоупругой жидкости Борна:
Gij = -Pbij + Gij Daи Doji ~ (Ю)
Xl~Dt + t2_Dí" + = Лву +
Несжимаемая изотропная вязкоупругая жидкость с внутренними вращениями характеризуется двумя коэффициентами вязкости Т), Г)г и двумя временами релаксации tt, т2. Если пренебречь эффектами вращения вокруг директора и релаксационными (полимерными) свойствами среды, определяющее уравнение (5) сводится к известному
~ r\¡ij)ke^ke
(8)
= 4(ij)ke£-ke
уравнению Лесли-Эриксена для низкомолекулярных жидких кристаллов [24]
= ампрмеь + агп1Ы] + а3АГ^ + + аАец + а 5п(пкек] + а6е1кпкп,,
гдеN¡=11-(о1епе = -Щепе. Материальные параметры а, обычно называются коэффициентами вязкости Лесли. Более простая теория анизотропной жидкости была предложена Эриксеном [15]. Она является специальным случаем жидкости Лесли-
Эриксена, когда Оу = 0, о^ = сх3 = 0 и а5 = а6. Анизотропная жидкость Эриксена характеризуется симметричными напряжениями и тремя независимыми коэффициентами вязкости.
ОРИЕНТАЦИОННАЯ ДИНАМИКА
Динамика директора п, входящего в реологические уравнения полимерных нематических фаз, описывается с помощью дополнительного уравнения, которое характеризует ориентационные свойства среды. Уравнение движения директора анизотропной вязкоупругой жидкости можно вывести непосредственно из уравнения для изменения внутреннего момента импульса [4, 30, 31]
= о, + е + т„ о, = е^Од (12)
(ц,^ - несимметричный тензор моментных напряжений). Внутренний момент импульса определяется соотношением = /1(,со', где /у = /±5,; + (/ц -- - тензор вращательной инерции; /,-, явля-
ется дополнительной характеристикой среды и зависит от природы молекул и их взаимной ориентации. Угловая скорость внутреннего враще-
/ I I
ния со = «»и + (й± складывается из скорости вращения (переориентации) директора со[, = г^пк и
/ I
скорости вращения вокруг директора С0ц, = соепеп1.
Уравнение (12) учитывает три вращательные степени свободы среды и показывает, что изменение внутреннего момента импульса происходит по трем причинам: в результате действия момента внешних сил т, вследствие действия поверхностных моментных сил е и как результат процесса преобразования внешнего момента импуль-
са во внутренний, описываемого вектором о,. Если этот вектор отличен от нуля, то тензор напряжений Су несимметричен.
Анализ влияния магнитного поля на жидкие кристаллы обычно основывается на предположении о том, что оно взаимодействует с жидкостью только через момент ш, который определяется в виде
ш = МН, (13а)
где Н - напряженность магнитного поля, М - вектор намагниченности. При этом для намагниченности в основном используется линейное соотношение
М = 5С±Н + Ха(пН)п (136)
Здесь Ха = Х\\ ~ Х±'а 1\\и Хх - главные значения тензора динамической восприимчивости вдоль и поперек директора соответственно. Магнитный момент ш = пН"1 можно рассматривать как результат действия на директор эффективного поля Ьт = = Ха(п Н)Н.
Используя уравнение вращательного движения (12), находим [29] ориентационное уравнение
1Пепе - /„«¡со!, = /г«1 + 2педае, (14)
и спиновое уравнение
1\М\ = (15)
которое описывает вращение среды вокруг директора. Здесь = - пп) - поперечная компонента тензора вращательной инерции, а Ь1 = = Ь - пп • 11 - поперечная компонента молекулярного поля. Внутренние упругие силы ЖК-немати-ка и внешние силы, действующие на директор со стороны магнитного поля, могут быть определены с помощью единого молекулярного поля Н = = Ье + Ьт, где Ьг играет роль молекулярного поля, стремящегося установить одинаковое направление директора во всем объеме нематика. Согласно работе [13], оно определяется как /г, = Ке1 е - Э/^/Эи,, где яу = Э/^/Эпд,, - свободная энергия Франка, связанная с искажением поля директора.
При выводе уравнений (14) и (15) учтен только равновесный вклад в моментные напряжения. Он имеет вид [12,32]
М-е/ -
(16)
Здесь не рассматриваются эффекты, связанные с неоднородными внутренними вращениеми, поэтому динамический вклад в моментные напряжения равен нулю.
Согласно уравнениям (14) и (15), ориентацион-ная и спиновая динамика определяется антисимметричными напряжениями, которые в свою очередь описываются релаксационным уравнением (8). Их анализ показывает, что инерционные и вязкоупругие эффекты могут быть связаны в ориентационной динамике вязкоупругих анизотропных сред при больших временах релаксации. Из моментного уравнения (12) видно, что если можно пренебречь инерцией внутреннего вращения и эффектами ориентационной упругости, то тензор напряжений является симметричным в отсутствие внешних полей. При этом получаем уравнение ориентации директора для полимерных нематиков
Dnl
= Р,
(17)
с релаксирующеи ориентационнои силои
где хи = х(51е - «,ие), йщ/£>/ = и, - со¡епе - относительная скорость директора, х - ориентационное время релаксации. В случае вязкой нематодина-мики оно сводится к известному уравнению ориентации Эриксена
Рп,
= Че^п.-щп^е^)
(18)
(к - безразмерный материальный параметр, определяющий характер ориентационного процесса).
Таким образом, ориентационная динамика полимерных (вязкоупругих) нематиков зависит от ускорения директора. Данным эффектом можно пренебречь для процессов переориентации директора, происходящих с малыми ускорениями.
градиент скорости
направление течения
Рис. 1. Касательные напряжения при простом сдвиговом течении.
Тогда для вязкоупругих нематиков можно использовать уравнение ориентации Эриксена, как это было сделано в работе [1]. В общем случае анализ ориентационной динамики высокомолекулярных жидких кристаллов на основе уравнений ориентации Лесли-Эриксена не является корректным.
В работах [33-35] начаты экспериментальные исследования ориентационной релаксации ЖК-полимеров и обнаружен существенный рост времени релаксации директора нематических гребнеобразных полимеров по сравнению с низкомолекулярными системами.
ПРОДОЛЬНАЯ И ПОПЕРЕЧНАЯ ВЯЗКОСТИ
Для общего описания напряженного состояния любой сплошной среды требуется девять величин -одно нормальное и два касательных напряжения для каждой из трех координатных плоскостей. Общая совокупность напряжений сводится к девяти компонентам: три нормальных напряжения и три пары касательных напряжений. Для сред с симметричными напряжениями из всех девяти компонентов только шесть отличны друг от друга. Это обстоятельство является чрезвычайно важным, так как число неизвестных в общих уравнениях движения среды уменьшается на три.
Рассмотрим простое сдвиговое течение. Для жидкостей с несимметричными напряжениями оно характеризуется только одной парой [4] касательных напряжений о21 и <у12 (рис. 1). Первый индекс в них указывает плоскость, на которую действует напряжение, второй - направление его действия. Компоненты касательных напряжений, изображенные на рис. 1, положительны. Таким образом, анизотропные жидкости с несимметрич-
ными напряжениями, к которым относятся и жидкие кристаллы, при простом сдвиговом течении должны характеризоваться двумя вязкостями сдвига
где у - скорость сдвига. Продольная (относительно направления течения) вязкость сдвига Г|1 связана с касательным напряжением, действующим в плоскости, параллельной направлению течения; поперечная вязкость сдвига цг обусловлена поперечным касательным напряжением, действующим в плоскости, перпендикулярной направлению сдвига. Для обычных полимерных жидкостей (расплавов и растворов гибкоцепных полимеров) напряжения о21 и а12 совпадают и поэтому достаточно исследовать только одну вязкость сдвига. Традиционные ротационные и капиллярные реологические приборы позволяют измерять продольную вязкость только в одном направлении. Изучать ее анизотропию можно только на специально созданном оборудовании (см. ниже). Измерение поперечной (относительно направления течения) вязкости сдвига анизотропных полимерных жидкостей составляет предмет нетрадиционной вискозиметрии. До сих пор не существует подходов к ее измерению. Это дело будущего.
Важной особенностью нематических жидких кристаллов является то, что сдвиговые вязкости весьма чувствительны к ориентации директора. В линейной области они зависят только от начальной ориентации. Если фиксировать начальную ориентацию директора с помощью внешнего поля параллельно направлению течения (а), параллельно градиенту скорости (Ь) и перпендикулярно к этим двум направлениям (с), то можно определить три вязкости сдвига Г|", и Т)^. Их обычно называют вязкостями Месовича [36]. Для низкомолекулярных жидких кристаллов они могут быть вычислены на основе теории Лесли-Эрик-сена (см., например, работу [13]) и определяются как линейные комбинации вязкостей Лесли.
Согласно изложенному, для жидких кристаллов при простом сдвиговом течении в общем случае необходимо определять дополнительно три
сдвиговые (поперечные) вязкости Г|2, Л 2 и Л 2 •
Наш расчет показывает, что они связаны с вязкостями Лесли следующими линейными соотношениями:
Л 2 = |(а4 + а5 + а2)
12 = 5(04 + 06-03) (20)
с 1
Л2 = 2а4
Необходимо отметить, что простое сдвиговое течение одноосных (нематических) жидкостей, когда отлична от нуля только одна пара касательных напряжений, наблюдается не для всех ориен-таций директора. Можно показать, что оно существует тогда и только тогда, когда директор находится в плоскости сдвига, определяемой направлениями течения и градиента скорости, или перпендикулярен этой плоскости. В противном случае простой сдвиг сопровождается появлением при течении еще двух пар касательных напряжений о13, о31 и о2з, О32.
Для анизотропных полимерных систем ситуация существенно усложняется наличием анизотропных релаксационных процессов. В данном случае рассмотренные выше продольная и поперечная вязкости зависят как от ориентации директора, так и от характера релаксационной анизотропии. Проиллюстрируем этот эффект на примере вязкоупругих характеристик, определяемых при малоамплитудном осциллирующем сдвиговом течении для трех базовых начальных ори-ентаций, создаваемых магнитным полем. В области линейной вязкоупругости продольная и поперечная (относительно направления течения) комплексные вязкости сдвига зависят только от начальной ориентации директора.
"Асимметрию" вязкости сдвига можно оценить на основе простого релаксационного уравнения для динамических напряжений с одним анизотропным временем релаксации [4, 30, 31]
Х']ке~5^ + = Лчь^ь + (21)
где N - вектор, характеризующий относительную скорость вращения (переориентации) директора. Оно является частным случаем релаксационного реологического уравнения (5) и учитывает толь-
ко две вращательные степени свободы, пренебрегая эффектами, связанными с вращениями вокруг директора. В случае нулевого тензора времен релаксации Ту!«, - 0 релаксационное уравнение (21) сводится к реологическому уравнению теории Лесли-Эриксена (11). Для несжимаемых немати-ческих (одноосных) жидкостей тензоры вязкости Цу/се, рр и тензор времен релаксации Хуке являются функциями директора и определяются в виде
= Л5«.&д + Лб"А* Л Цке = Л11щке) + Л гПцке + Л з"А'(*"«) + Л^Ах*"*)
= Ъу11](ке) + Х2П1]ке + ЬП&еПк + Т4"АЛ + ^ + х5п&кпе + т6^8;кпе
Здесь 1ике = 8,А» - единичный тензор четвертого ранга, пуке = щп^,,, п^(кпе} = | (п£]кпе + п^епк).
В силу слабой намагниченности рассматриваемых анизотропных полимерных жидкостей влиянием магнитного поля Н на материальные функции можно пренебречь. Это значительно упрощает их реологические уравнения. Входящие в выражения (22) материальные параметры Г|, и x¡ могут быть найдены из эксперимента или в рамках молекулярных теорий. Примеры вычислений коэффициентов вязкости Лесли (анизотропной вязкости) из простых молекулярных соображений хорошо известны [5-8].
Частотные зависимости продольной и поперечной вязкости были оценены [4], исходя из определяющих уравнений (21), (22) для простого случая независимой релаксации симметричных и антисимметричных напряжений. На рис. 2 и 3 представлены характерные безразмерные частотные зависимости продольной и поперечной комплексных вязкостей для трех базовых начальных ориентаций директора а, Ь, с. Анализ этих кривых показывает, что динамические свойства полимерных нематиков существенным образом зависят как от ориентации директора, так и от анизотропии их релаксационных свойств. При этом частотные зависимости поперечной вязкости обнаруживают особенности, не присущие изотропным полимерным жидкостям. Наибольшее отклонение от изотропного поведения жидкости Максвелла (штриховые кривые на рис. 2) при простом сдвиге происходит тогда, когда директор ориентирован параллельно градиенту скорости.
П ^/л !а(0) 1.0
-4
Л1С(*)/Л1С(0) 1.0
Рис. 2. Зависимости действительной (7) и мнимой (2) частей продольной комплексной вязкости л !(*)/?! х(0) от безразмерной частоты х для трех ориентаций директора: а = а (а), Ь (б) и с (в). Штриховые кривые характеризуют вязкоупру-гость жидкости Максвелла.
Следует отметить, что для ЖК-систем вязкость, соответствующая ориентации директора, параллельной направлению градиента скорости, также максимально отличается от изотропного случая. Как правило, она наибольшая [37-40].
ИЗМЕРЕНИЕ ПРОДОЛЬНОЙ ВЯЗКОСТИ И АНИЗОТРОПНОЙ ВЯЗКОУПРУГОСТИ
Впервые продольная вязкость, определяемая касательным напряжением с21> Для низкомолекулярных жидких кристаллов была измерена
Т12в(*)/тъ°(0)
Ц^хЩ^О)
штт
Рис. 3. Безразмерные частотные зависимости действительной (]) и мнимой (2) компонент поперечной комплексной вязкости Т| 2 (дс)/Г12 (0) для трех начальных ориентаций директора а = а (а), Ь (б) и с (в).
ЭсЬепск [41]. На температурной зависимости вязкости был обнаружен максимум в точке перехода в ЖК-состояние. Более поздние работы Месови-ча [36], Цветкова и Михайлова [42] и других авторов показали наличие зависимости продольной вязкости низкомолекулярных жидких кристаллов от ориентации директора относительно течения. Современное состояние этой проблемы отражено в обзоре [43]. В нем описаны различные
методы измерения продольной вязкости немати-ческих жидких кристаллов, наиболее широко используемых в практике.
При исследовании анизотропии вязкости ЖК-систем важное значение имеет возможность создания различных (хорошо определенных) ориентаций среды относительно направления течения. Для низкомолекулярных жидких кристаллов это достигается приложением сильных магнитных полей. Для ЖК-полимеров (особенно линейного строения) определенная ориентация в магнитном или электрическом поле разумной напряженности труднодостижима из-за высокой вязкости высокомолекулярных соединений.
Не обсуждая сейчас проблему ориентации полимерных жидких кристаллов с помощью электрического и магнитного поля для реологических исследований (хотя и известны измерения вращательной вязкости гребнеобразных ЖК-полиме-ров в магнитном поле ЯМР-спектрометра), укажем некоторые подходы к определению продольных коэффициентов вязкости и динамических модулей упругости при различных ориентациях среды относительно направления течения. В качестве ориентирующего средства использовали механическое поле или эффект взаимодействия анизотропной полимерной системы с подложкой. Мы не можем претендовать на количественное сопоставление экспериментальных данных с предсказаниями теории, однако о качественном сопоставлении можно говорить уже сейчас.
Первые попытки измерения анизотропной продольной вязкости были сделаны в 70-х годах XX столетия на примере анизотропных растворов поли-п-бензамида в ДМ А А с ЫС1 [44, 45]. Фактически речь шла об оценке трех коэффициентов вязкости Месовича для полимерных жидких кристаллов. Что касается первых двух коэффициентов вязкости т]? и г], , соответствующих параллельной и перпендикулярной (относительно направления течения) ориентации среды, то для их оценки впервые была использована простая и оригинальная схема измерения на ротационной паре цилиндр-цилиндр. Одновременно с регистрацией крутящего момента на валу внутреннего (подвижного) цилиндра измеряли и продолжительность падения стального шарика в зазоре при течении анизотропного раствора (рис. 4а). По смыслу рассчитанная по величине крутящего
(б)
[отн. ед.] 0.51-
Ч^1
со.
Рис. 4. Схемы измерительных узлов для оценки анизотропии продольных вязкости (а) и комплексного модуля упругости (б).
момента вязкость Г|ц должна быть близка к , а рассчитанная по Стоксу вязкость Г)х - к вязкости
т^ 1. Зависимости соответствующих вязкостен от скорости сдвига приведены на рис. 5. Как и ожидалось, вязкость ЛиСЛ?) снижается, а вязкость
Т|^(Г|1) возрастает с увеличением скорости деформации. При определенной скорости обе вязкости достигают постоянства, которое можно интерпретировать как достижение максимально возможной в данных условиях молекулярной ориентации (аномалия вязкости для этих растворов несущественна). Впоследствии ситуация с движением шарика в анизотропной среде была проанализирована теоретически [2,46].
Для измерения вязкости был использован другой оригинальный подход [45], основанный на определенной обработке поверхности стеклянной ячейки (плоского капилляра), приводящей к формированию гомеотропной текстуры ЖК-рас-твора (длинные оси арамидных макромолекул ориентированы перпендикулярно поверхности стекол). Наклон такой ячейки приводит к деформированию препарата, хотя в этом случае отме-
Рис. 5. Зависимость коэффициентов продольной вязкости г|± (7) и г|ц (2) от продольного градиента скорости Уц для 10%-ного раствора по-ли-и-бензамида в ДМАА.
1§<[Пас]
мин
Рис. 6. Кинетика изменения напряжения сдвига х (1) и вязкости т^ (2) гомеотропной текстуры ЖК-раствора поли-и-бензамида в наклонной плоской ячейке.
чается наличие предела текучести (рис. 6). После преодоления предела текучести начинается течение практически неразрушенной структуры, скорость которого измеряется микроскопом по скорости перемещения выбранного дефекта текстуры. Такое течение отвечает гомеотропной
ориентации ЖК-препарата. Наступающее вскоре лавинообразное снижение вязкости, сопровождающееся просветлением поля зрения между скрещенными поляроидами, свидетельствует о частичном переходе при сдвиге гомеотропной текстуры в планарную. Замедление потока путем уменьшения угла наклона ячейки с определенным запаздыванием вызывает повышение вязкости, сопровождающееся потемнением поля зрения поляризационного микроскопа (формирование исходной гомеотропной текстуры).
В результате впервые для полимерных жидких кристаллов при относительно низких скоростях сдвига были оценены три коэффициента вязкости, близкие по смыслу к коэффициентам вязкости Месовича: т^ = 6.5 Па с, т^ =35 Па с и
Т| 1 = 12 Па с. Порядок величин коэффициентов вязкости такой же, как и для низкомолекулярных жидких кристаллов, например метоксибензили-ден-и-бутиланилина (ММБА) [43]: минимален коэффициент вязкости Т)" и максимален коэффи-ь
циент вязкости %.
Следующим шагом была разработка подхода к экспериментальному определению анизотропных вязкоупругих свойств; Движущей силой экспериментов такого рода является необходимость использования для анизотропных сред соответствующего "анизотропного реометра", позволяющего проводить измерения динамических модулей в двух ортогональных направлениях, а не традиционных приборов для определения изотропных вязкоупругих характеристик. Прибор ВР-04 был сконструирован и изготовлен в СКВ ИНХС РАН в середине последнего десятилетия XX века [47].
При использовании рабочего узла в виде двух коаксиальных цилиндров (рис. 46) возникает возможность тестирования анизотропното образца (ЖК-систем, наночастиц Иа-монтмориллонита, распределенных в полимерной матрице, анизоди-аметричных капель дисперсной фазы в расплавах смесей полимеров и т.п.) вдоль и поперек направления его преимущественной ориентации. Главной особенностью прибора является возможность предварительной ориентации испытуемой среды непрерывным сдвигом с последующим воздействием на ориентированный образец малоамплитудных частотных деформаций в угловом
(вдоль ориентации) и осевом (поперек ориентации) направлениях. Более того, в приборе предусмотрен режим наложения синусоидальных колебаний на непрерывный сдвиг, что с учетом возможности работы при температурах до 300°С делает его единственным в мире аппаратом такого рода. В работах бельгийской группы реологов аналогичные возможности достигаются использованием двух различных приборов, изготовленных и испытанных позже виброреометра ВР-04 [48].
Таким образом, прибор позволяет измерять продольную вязкость в установившемся или динамическом (частотном) режиме, а также продольный комплексный модуль вдоль и поперек направления преимущественной ориентации. К сожалению, прибор не позволяет строго задавать и контролировать начальную ориентацию образца, поэтому в начале опыта образец подвергается интенсивному сдвигу в течение определенного времени для создания ориентации, после чего проводятся частотные измерения. Как правило, релаксация созданной непрерывным сдвигом ориентации происходит в течение продолжительного времени, поэтому можно считать, что она не изменяется в процессе определения комплексного динамического модуля. В пользу этого соображения свидетельствует то обстоятельство, что измерение комплексного модуля при первой, самой низкой частоте после снятия частотной зависимости приводит к тому же результату, что и в начале опыта.
Итак, чего мы можем ожидать от эксперимента для подтверждения предсказаний теории? Прежде всего существенного различия продольных реологических характеристик вдоль и поперек директора. В качестве примера функций такого рода приведем частотные зависимости продольного модуля накопления G' и продольной динамической вязкости r|' = G"/со для водного ЖК-раствора гидроксипропилцеллюлозы Klucel (фирма "Hercules", M = 1.3 х 105), измеренные в различных направлениях в чисто колебательном режиме и в режиме суперпозиции (рис. 7а и 76 из работы [49]).
В принципе получен ожидаемый результат, а именно: при наложении угловых колебаний на установившееся сдвиговое течение происходит снижение, а в случае осевых колебаний - повышение
[Па], 1ёту [Пас] (а)
18СГ [Па]
[с-1]
Рис. 7. Частотные зависимости модуля накопления (1,3) и динамической вязкости (2,4) для водных ЖК-растворов гидроксипропилцеллюлозы в режиме угловых колебаний и суперпозиции (а), а также в режиме осевых колебаний и суперпозиции (б), у = 0 (7,2) и 0.11 с"1 (3,4).
модуля упругости и динамическои вязкости по сравнению с данными для чисто колебательных режимов. Обращает на себя внимание асимметрия эффекта, проявляющаяся в том, что динамические характеристики при угловой (параллельной) суперпозиции снижаются в меньшей степени, чем возрастают при осевой (перпендикулярной) суперпозиции. Не исключено, что это может быть обусловлено выходом директора из плоскости сдвига и нарушением условий простого сдвига.
Другой пример относится к жидким прекурсорам нанокомпозитов, т.е. расплавам полипропилена, наполненным 10% слоевого силиката - Иа-монтмориллонита. Поверхность элементарных
3.5
2.5
1.5
-1
1
[с-1]
Рис. 8. Частотные зависимости модуля потерь для расплава полипропилена, наполненного 10% модифицированного ^-монтмориллонита, фазу после загрузки (7) и при угловых (2) и севых (3) малоамплитудных колебаниях после 30-минут-
ного ориентирующего сдвига при у = 180°С.
= 1 с"1. Т =
пластин глины модифицирована диблок-сополи-мером с блоками ПЭО и сополимера этилена с изобутиленом. В этом случае для ориентации системы использовали предварительный сдвиг со скоростью 1 с-1 в течение 30 мин, после чего определяли частотные зависимости модулей в одном из взаимно перпендикулярных направлений. Затем снова включали ориентирующий непрерывный сдвиг с той же скоростью и на то же время, что и в первом случае, и измеряли частотные зависимости модулей в ортогональном направлении.
Экспериментальные результаты для модуля потерь показаны на рис. 8. И для этой системы наблюдается существенная анизотропия продольного модуля потерь, а следовательно, и динамической вязкости. Как и в предыдущем случае, в направлении ориентации модуль ниже, чем в направлении, перпендикулярном ей. Здесь также имеет место небольшая асимметрия эффекта: прирост модуля в ортогональном направлении ниже, чем его снижение в направлении, параллельном ориентации. По-видимому, специфическая ориентация директора с выходом элементов структуры из плоскости сдвига является скорее правилом, чем исключением.
Таким образом, во всех случаях продольные реологические характеристики, измеренные в направлении ориентации директора, оказываются значительно ниже соответствующих характеристик в трансверсальном направлении.
ЗАКЛЮЧЕНИЕ
В данной работе представлена общая методология анализа анизотропной вязкоупругости ЖК-полимеров и жидких полимерных нанокомпози-тов, которые характеризуются анизотропными релаксационными свойствами и несимметричными напряжениями. Особенность этих систем состоит в том, что они сочетают в себе анизотропные свойства жидких кристаллов и релаксационные (вязкоупругие) свойства полимеров. ЖК-полимеры и полимерные нанокомпозиты рассматриваются с позиций анизотропных вязкоупругих жидкостей с внутренними вращениями. Предложенные из простых соображений симметрии определяющие уравнения описывают наряду с анизотропией вязкости анизотропию динамических характеристик и ориентационной релаксации. Представленный теоретический анализ ограничивается областью квазилинейной динамики.
Прямая экспериментальная проверка теории несимметричной вязкоупругости в отсутствие внешних полей в настоящее время невозможна, поскольку не существует рациональных подходов к измерению поперечной вязкости сдвига. Что же касается анизотропии продольной вязкости, то она была экспериментально исследована на примере ЖК-растворов ароматических полиамидов. Анизотропия продольных вязкоупругих характеристик была измерена для ЖК-растворов окси-пропилцеллюлозы, ЖК-расплавов алкилен-аро-матического сополиэфира, их смесей с изотропными растворами и расплавами соответственно [49], а также для расплавов изотактического полипропилена, наполненного частицами Ка-монт-мориллонита. Последняя система представляет собой жидкий предшественник нанокомпозита на основе полимерной матрицы и частиц слоевых силикатов, способных формировать ЖК-фазу.
Основным результатом проведенных экспериментов является обнаружение релаксационной анизотропии в ЖК-полимерах и нанокомпозитах. Этот факт поддерживает основное положение те-
ории о существовании в полимерных системах анизотропных релаксационных процессов.
Теоретическое рассмотрение сильно нелинейных вязкоупругих эффектов в ЖК-полимерах достаточно сложно. Заметного продвижения в этом направлении следует ожидать в рамках термодинамического подхода. В работах [10, 11, 50] начаты исследования анизотропной нелинейной вязкоупругости разных полимерных систем с различной степенью жесткости в рамках термодинамики необратимых процессов. Весьма существенным здесь является то, что физика таких сплошных сред не может быть описана в общей форме ни механикой, ни электро(магнито)динамикой в отдельности. Только термодинамика необратимых процессов дает возможность получить общее описание макроскопических свойств анизотропных полимерных систем. При этом в рамках единого подхода можно исследовать анизотропную вязкоупругость полимерных жидкостей и твердых или высокоэластических тел, например ЖК-эластомеров [51,52]. Неравновесная термодинамика дает общую основу для описания макроскопических свойств различных анизотропных полимерных систем в равновесной и неравновесной ситуациях.
СПИСОК ЛИТЕРАТУРЫ
1. Volkov VS., Kulichikhin V.G. // J. Rheol. 1990. V. 34. P. 281.
2. Kulichikhin V.G., Volkov VS., Plate NA. Comprehensive Polymer Science, First Supplement / Ed. by Aggar-wal S., Russo S. New York: Pergamon Press, 1992. P. 385.
3. Волков B.C., Куличихин В.Г. // Высокомолек. co-ед. A. 1994. Т. 36. № 11. С. 1776.
4. Volkov VS., Kulichikhin V. G. I I Rheol. Acta. 2000. V. 39. P. 360.
5. Doi M., Edwards S.F. The Theory of Polymer Dynamics. New York: Oxford Press, 1986.
6. Семенов A.H., Хохлов A.P. // Успехи физ. наук. 1988. Т. 156. № 3. С. 427.
7. Marrucci G., Greco F. // Adv. Chem. Phys. 1993. V. 86. P. 331.
8. Larson R.G. The Structure and Rheology of Complex Fluids. Oxford: Oxford Univ. Press, 1998.
9. De Gennes P.G. Liquid Crystals in One- and Two Dimensional Order / Ed. by Helfrich W„ Kleppke G. Berlin: Springer, 1980. P. 231.
10. Leonov A J., Volkov V.S. // Los Alamos e-print arxive: cond mat. ArXiv.org. 2002. N. 0202275.
11. Leonov A.I., Volkov V.S. // Инж.-физ. журн. 2003. Т. 76. № 3. С. 24.
12. Аэро А.Э., Булыгин А.Н. // Гидромеханика. Итоги науки и техники. 1973. Т. 7. С. 106.
13. De Gennes P.G., Prost G. The Physics of Liquid Crystals. Oxford: Clarendon Press, 1993.
14. Кац Е.И., Лебедев B.B. Динамика жидких кристаллов. М.: Наука, 1988.
15. EricksenJL.U KolloidZ. 1960. V. 173. Р.117.
16. Виноградов Г. В. // Механика полимеров. 1977. №6. С. 1062.
17. Якобсон Э.Э., Файтельсон Л.А. // Механика композит. материалов. 1985. № 2. С. 328.
18. Volkov V.S., Vinogradov G.V. // Rheol. Acta. 1988. V. 26. P. 96.
19. Navier M. // Mem. de 1'Acad, des Sci. 1826. № 6. P. 389.
20. Stokes G.G. I I Math, and Phys. Pap. 1847. № 1. P. 75.
21. Born M. // Z. Physik. 1920. V. 1. P. 221.
22. Сорокин B.C. // Журн. эксперим. и теорет. физики. 1943. Т. 13. № 7-8. С. 306.
23. Цветков В.Н. // Журн. эксперим. и теорет. физики. 1939. Т. 9. № 5. С. 603.
24. Leslie F.M. // Arch. Ration. Mech. Anal. 1968. V. 28. P. 265.
25. Carlson Т., Leslie F.M. // Liq. Cryst. 1999. V. 26. № 9. P. 1267.
26. Allen SJ„ de Silva C.N. // J. Fluid Mech. 1966. V. 24. P. 801.
27. Аэро ЭЛ., Булыгин А.Н. // Гидромеханика. 1973. Т. 7. № 1. С. 106.
28. Volkov V.S. И Macromol. Symp. 2000. V. 160. P. 261.
29. Volkov V.S., Kulichikhin V.G. // Proc. 6 Eur. Conf. on Rheology / Ed. by Muenstedt H., Kaschta J., Merten A. Erlangen, Germany, 2002. P. 407.
30. Волков B.C. // Высокомолек. соед. A. 1998. Т. 40. №7. С. 1150.
31. Volkov V.S. Proc. 4 Eur. Conf. on Rheology "Progress and Trends in Rheology" / Ed. by Emri I. Portoroz. Slovenia, 1998. P. 240.
32. Leslie F.M. // Continuum Mech. Thermodyn. 1992. V. 4. P. 167.
33. Götz S., Stille W., Stroble G., Scheuermann H. // Macro-molecules. 1993. V. 26. № 7. P. 1520.
34. Van der Putten D„ Schwenk N.. Spiess H.W. // Liq. Cryst. 1989 V. 4. № 3. P. 341.
35. Bock FJ., Kneppe H„ Schneider F. // J. Phys. Chem. 1989. V. 93. P. 3848.
36. Miesowicz M. // Nature. 1946. V. 158. № 4001. P. 27.
37. Porter R. S., Johnson J. F. Ц Rheology / Ed. by Erich F.R. New York: Acad. Press, 1967. V. 4.
38. Gahwiller Ch. // Mol. Cryst. Liq. Cryst. 1973. V. 20. P. 301.
39. Kneppe H., Schneider F. // Mol. Cryst. Liq. Cryst. 1981. V. 65. P. 23.
40. Benicewicz B.C., Johnson J.F. // Mol. Cryst. Liq. Cryst.
1981. V. 65. P. 111.
41. SchenckR. //Phys. Chem. 1898. V. 27. P. 167.
42. Tsvetkov V.N., Mikhailov G.M. I I Acta Physicochim. USSR. 1938. V. 8. P. 77.
43. Беляев B.B. // Успехи физ. наук. 2001. Т. 171. № 3. С. 267.
44. Малкин А.Я., Васильева Н.В., Белоусова ТА., Ку-личихин В.Г. // Коллоид, журн. 1979. Т. 61. № 1. С. 200.
45. Kulichikhin V.G. // Mol. Cryst. Liq. Cryst. 1989. V. 169. Ms 1. P. 51.
46. Волков B.C. // Прикл. математика, и механика.
1982. Т. 46. № 2. С. 248.
47. Kulichikhin V., Volkov V., Vasil'eva О. / Proc. Polymer Processing Society Meeting. Stuttgart, Germany, 1995. P. 4.2.
48. Vermont J., Walker L., Moldenaers P., Mewis J. // J. Non-Newton. Fluid. Mech. 1998. V. 79. № 1. P. 173.
49. Куличихин В.Г., Плотникова Е.П., Терешин А.К., Субботин A.B., Цамалашвили Л.А. // Высокомолек. соед. С. 2000. Т. 42. № 12. С. 2235.
50. Leonov A.I., Volkov V. S. I I Proc. 6 Eur. Conf. on Rheology / Ed. by Muenstedt H., Kaschta J., Merten A. Erlangen, Germany, 2002. P. 315.
51. Leonov A J., Volkov V.S. // Los Alamos e-print arxive: cond mat. ArXiv. org. 2003. Ms 0302128.
52. Леонов A.M., Волков B.C. // Инж.-физ. журн. 2004. Т. 77. №4. С. 36.
Viscoelasticity of Anisotropic Polymer Systems
V. S. Volkov, A. K. Tereshin, and V. G. Kulichikhin
Topchiev Institute of Petrochemical Synthesis, Russian Academy of Sciences, Leninskii pr. 29, Moscow, 119991 Russia
Abstract—Results of the general analysis of anisotropic viscoelasticity of LC polymers and polymer nanocom-posites in a liquid state are reported. These systems are considered as anisotropic viscoelastic liquids with asymmetric stresses. The oscillating shear flow of polymer nematics is theoretically analyzed at various initial director orientations as an illustration of relaxation anisotropy. The relaxation equation of orientation is derived. The anisotropy of longitudinal viscosity and longitudinal dynamic moduli for LC polymers and polymer suspensions containing plain silicate nanoparticles are experimentally measured for the first time. The results obtained agree with theoretical predictions.