Моделирование процесса склеивания древесины клеем, обработанным в электрическом поле

УДК 674.028

МОДЕЛИРОВАНИЕ ПРОЦЕССА СКЛЕИВАНИЯ ДРЕВЕСИНЫ КЛЕЕМ, ОБРАБОТАННЫМ В ЭЛЕКТРИЧЕСКОМ ПОЛЕ

UDC 674.028

MODELING OF PROCESS OF WOOD GLUING WITH A GLUE, PROCESSED IN ELECTRIC FIELD

Мозговой Николай Васильевич д.т.н., профессор

Воронежская государственная лесотехническая академия, Воронеж, Россия

Предложена математическая модель процесса склеивания древесины клеями, подвергнутыми обработке электрическим полем. Модель позволяет определять зависимость прочности клеевых соединений от параметров технологического процесса

Ключевые слова: ЭЛЕКТРИЧЕСКОЕ ПОЛЕ, КЛЕЕВОЕ СОЕДИНЕНИЕ, НАПРЯЖЕННОСТЬ, ПРОЧНОСТЬ, МОДЕЛЬ

Mozgovoy Nokolay Vasilyevich Dr.Sci.Tech., professor

Voronezh State Academy of Forestry and Technologies, Voronezh, Russia

A mathematical model of the process of adhesive bonding of wood with glues, treated by an electric field, is suggested. The model allows determining the dependence of the strength of adhesive joints on technological process parameters

Keywords: ELECTRIC FIELD, ADHESIVE JOINT, INTENSITY, STRENGTH, MODEL

Изделия из клееной древесины в последнее время составляют основной объем продукции, которую выпускают современные деревообрабатывающие предприятия. При склеивании изделий из массивной древесины основным требованием к клеевым соединениям является их высокая прочность. Для повышения прочности клееной древесины проводятся работы по совершенствованию технологии склеивания, разрабатываются новые марки клеев. Однако, эти мероприятия уже не в состоянии удовлетворить растущие требования производственников и потребителей готовой продукции. Более перспективными и дающими большой эффект повышения прочности клеевых соединений выглядят способы воздействия на клей физическими полями, приводящие к перестройке структуры полимерной основы клея. Достаточно эффективным показал себя способ обработки клеев постоянным магнитным полем [1, 2]. Аналогичное воздействие на структуру клея следует ожидать при его обработке электрическим полем [3].

Для реализации экспериментальной программы по изучению процесса склеивания на клее, подвергнутом воздействию электрическим полем были проведены теоретические исследования этого процесса. На этом эта-

пе была поставлена задача создать математическую модель процесса склеивания древесины наиболее часто применяемыми на деревообрабатывающих предприятиях клеями марок КФ-МТ-15 и КФ-Ж, которые предварительно были подвергнуты воздействию постоянным электрическим полем. Модель позволила бы изучить зависимость прочности клеевого соединения от параметров технологического процесса и на основе этого оптимизировать технологический процесс.

Для исследования сложных систем, к которым относятся полимерные клеи, на микроскопическом уровне в настоящее время наибольшее распространение получил метод молекулярной динамики [4, 5]. С его помощью моделируется структура и физические свойства твердых тел, жидкостей и газов, а также процессы плавления, кристаллизации, испарения и конденсации. В методе молекулярной динамики каждая микроскопическая частица (атом, ион, молекула) рассматривается как геометрическое тело простой формы (шар, диск), движущееся в соответствии с законами классической динамики Ньютона и взаимодействующее с другими частицами по определенному закону.

Основные особенности поведения реальной системы могут быть воспроизведены в двух измерениях, поэтому для большей наглядности получаемых результатов была построена двумерная модель. Исключения третьего измерения дает возможность при заданном числе элементов (в расчетах использовалось до 10000 элементов) увеличить линейные размеры образца материала в направлениях X и У. Метилольные группы (мономеры) в модели представляются отдельными круглыми элементами (в дальнейшем называемыми «молекулами»), имеющими три степени свободы, а именно: движение вдоль осей ОХ и ОУ и вращение вокруг собственной оси.

Каждая отдельная молекула в модели представляет собой электрический диполь с зарядами +q и -д, находящимися на расстоянии ё друг от

друга. Состояние і-й молекулы описывается тремя переменными: декартовыми координатами ее центра (х, у) и углом поворота фі і-й молекулы относительно оси ОХ. Для формирования начальной конфигурации 10000 кругов-молекул случайным образом помещаются в ячейку в форме прямоугольника. На всех границах прямоугольника используются зеркальные граничные условия: молекула, которая соприкасается со стенкой, испытывает абсолютно упругое соударение и возвращается в расчетную ячейку. Силы межмолекулярного взаимодействия складываются из взаимодействия, обусловленного отталкиванием электронных оболочек молекул на малых расстояниях, а также электростатического диполь-дипольного взаимодействия между молекулами. Отсюда сила, действующая на і-ю молекулу со стороны других молекул системы, равна

N ґ   Л

Fi=і к'+г”),

1=1 V У

У  эл

где Г, - сила упругого взаимодействия молекул і и у; Г, - сила электростатического взаимодействия і и у; N - число молекул в модели.

При расчете сил для каждой пары молекул предварительно вычисляется расстояние Ту между их центрами 8І (хі,уі) и S] (х,,у,):

Ту = .

Первый вид межмолекулярного взаимодействия описывается в рамках представления об упругом взаимодействии. Зависимость силы взаимодействия двух молекул і и у от расстояния между ними (рисунок 4.4) описывается следующими формулами:

С(т0 - Ту)(X - ху)/ту> если Ту < Т0 + ^;

0, если т.. > т0 + d0

ГУ =

гух

ГУ =

Гу¥

С(Т„ -Т,)(Уі -у,)/ті,; если Ту < г + d0;

0, если г, > г0 + d0,

где гух и ^,т - проекции силы Ё У на координатные оси ОХ и ОУ соответственно;

С - жесткость упругого взаимодействия молекул; г0 - диаметр молекулы;

г0 + й0 - расстояние, на котором молекулы упруго взаимодействуют между собой (^ = 0,2г0).

Рисунок 1. Зависимость силы упругого взаимодействия двух молекул I и] от расстояния между ними.

Моделирование производится при участии большого числа молекул и требует значительных затрат времени, поэтому целесообразно оптимизировать алгоритм с целью максимального уменьшения времени счета. Для этого используется следующий прием. Поскольку основной объем вычислений связан с нахождением сил межмолекулярного взаимодействия (3), то в процессе расчета сил упругого взаимодействия проверяется условие

г2 < Г2вз,

где гвз - радиус области, в пределах которой учитывается взаимодействие.

Вычисление функций Fу производится только при выполнении условия (4), в противном случае указанные величины считаются равными нулю. Так как на каждом шаге молекулы смещаются на расстояния, значительно меньшие среднего расстояния между ними, то расчет сил можно значительно ускорить с помощью следующей процедуры. Для каждой молекулы составляются массивы с номерами частиц, находящихся на расстояниях, не превышающих гвз от данной. Для i-й молекулы рассчитывается взаимодействие только с частицами, попавшими в соответствующий массив.

Электростатическое взаимодействие между электрическими зарядами диполей осуществляется по закону Кулона:

F3Jl = qqj r

4^rj ^

где q и qj - электрические заряды в узлах i и j, равные ± e ; e0 - электрическая постоянная; rj - радиус-вектор, направленный от заряда i к заряду j.

При приложении внешнего электрического поля с напряженностью E на каждый электрический заряд будет действовать дополнительная сила, направленная вдоль поля и вычисляемая по формуле:

Fвн = qfi.

Таким образом, на каждый электрический заряд действует сила электростатического взаимодействия:

N

—* ш щ -* -*

F = ЕF* + FT .

j=1

На каждую молекулу действует момент силы, вызывающий ее вращение относительно центра:

M =

где Кл и - проекции силы Кэ на ось ОХ, действующей на заряды q1 и q2 соответственно; и - проекции силы Рхэл на ось ОУ, дейст-

вующей на заряды д1 и q2 соответственно; d - плечо диполя.

В методе молекулярной динамики численно решаются уравнения движения в форме Ньютона. Для поступательного движения N молекул система уравнений движения имеет вид:

где ш,, V , Гі - масса, скорость и радиус-вектор і-й молекулы соответственно; Кі (т ) - сила, действующая на і-ю молекулу; ґ - время, і - номер молекулы (і = 1, 2, ..., Щ.

Вращательное движение молекул описывается системой уравнений:

dVi

ші 1

1 dt

dTi

1 = V

к dt 1

I ^=м,

о dt

ю.

где I, Ц, ф - момент инерции, угловая скорость и угол поворота относительно оси ОХ 1-й молекулы; М1 (ф;) - момент силы, действующий на /-ю молекулу.

Для численного интегрирования системы уравнений (9) применялся модифицированный метод Эйлера [6]:

г"+1 = г" + А ґї" +

V 2ші у

К"

у"+‘ = V" +

2ші у

где Ж - шаг интегрирования по времени; п - номер шага интегрирования.

Во всех расчетах шаг интегрирования составлял Д1 = 20-10 15 с. Система уравнений (10) интегрировалась аналогично:

ф"+‘ =ф" +д I Ю +

V 21 У

МП

'Д !Л

V 21 У

Данный численный алгоритм отличается простотой, надежностью и устойчивостью, и связанная с ним глобальная погрешность имеет третий порядок для координаты и второй порядок для скорости. Он является са-мостартующим и не приводит к накоплению погрешности округления [7].

Для того, чтобы в модели корректно учесть температуру клея, начальное распределение молекул по скоростям задавалось по максвелловскому закону для некоторой температуры Т^ Для достижения системой равновесия при заданной температуре необходимо производить нормировку скоростей на каждом временном шаге по формуле

Г1  в тТ1,

где Ь - нормировочный множитель:

ь=

1

ЗМкГ,

/

N

1=1

где к - постоянная Больцмана.

На каждом шаге по времени осуществлялся контроль над потенциальной энергией системы Ж:

*=- Ыг,),

2*.]=1

где и (ту) - потенциальная энергия взаимодействия молекул I и у.

Процедура интегрирования уравнений движения по формулам (10) и (11) продолжалась до тех пор, пока значение потенциальной энергии сис-

темы уменьшалось. Для достижения системой равновесия требовалось около 20000 временных шагов, после чего величина W начинала осциллировать относительно среднего равновесного значения.

После релаксации системы в течение заданного времени рассчитывалась сила адгезии FA между древесиной и клеем:

FA = Fk +F ,

где F3n - сила электростатического взаимодействия между молекулами клея и древесины, обусловленная поляризацией молекул; Fk - обычная сила сцепления между клеем и древесиной, обусловленная прочими факторами.

Предел прочности на скалывание t пропорционален силе адгезии склеенных поверхностей, поэтому

t = Ft = к • Fa,

где S - площадь склеиваемых поверхностей; Ft - минимальная сила, которую нужно приложить к клеевому соединению для его разрушения; к - коэффициент пропорциональности.

Разработанная математическая модель представляет собой систему дифференциальных и алгебраических уравнений. Для удобства исследования системы уравнений составлена компьютерная программа "Программа для моделирования процесса склеивания древесины клеями КФ, обработанными в электрическом поле" на языке Object Pascal в интегрированной среде программирования Borland Delphi 7.0.

При запуске программы появляется интерфейсная форма с большим количеством окон ввода, в которых можно задать 23 основных параметра модели.

Созданное в объеме клея электрическое поле способствует упорядочению диполей. Уже существующие в клее маленькие области сонаправ-

ленных диполей начинают укрупняться и выстраиваться вдоль поля. Для того, чтобы изучить, как напряженность электрического поля влияет на доменную структуру и прочность получаемого в итоге клеевого соединения проведена серия компьютерных экспериментов с моделью. В рамках серии изменяли значение Е от 0 до 1200 В/см с шагом 150 В/см. Полученные данные представлены в виде графика на рис. 2. На графике видно, что с увеличением напряженности электрического поля происходит увеличение сил адгезии.

4,5x1 0',п 4,0x1 0''°

3,5x1 О'1"

3,0x1 0'ш 2,5x1 О'1"

^ 2,0x1 0',п 1 ,5x1 О'1"

1 ,0x1 О'1"

5,0x1 О'11 0,0

0 200 400 600 800 1000 1 200

£, В/см

Рисунок 2. Зависимость силы адгезии РА от напряженности электрического поля Е.

От времени выдержки клея в поле зависит то, насколько прочно зафиксируется созданная полем дипольная структура. Для всех экспериментов напряженность поля составляла 1200 В/см. Полученные данные представлены на рис. 3 из анализа которых следует, что время выдержки имеет параболическую зависимость, максимальные значения силы адгезии под действием электрического поля достигают на 700 шаге, что равняется примерно 12 мин.

4,5x1 0'1П 4,0x1 0'1П 3,5x1 0''°

^ 3,0x10''°

2,5x1 0''°

2,0x1 0''°

О 200 400 ВОО 800 1000 1200

и

Рисунок 3. Зависимость силы адгезии РА от времени ? выдержки клея

в электрическом поле.

От температуры отверждения клея зависит то, насколько быстро будет протекать реакция поликонденсации, от которой зависит величина внутренних напряжений в клеевом шве, а следовательно и предел прочности клеевого шва. Полученные данные представлены на рис. 4. Они показывают, что увеличение температуры до 30°С приводит к увеличению предела прочности клеевого шва до максимума, а затем при дальнейшем увеличении температуры выше 30°С происходит резкий спад прочности.

4.50Е-010 4.00Е-010 3.50Е-010 _ З.ООЕ-ОЮ

и_т

2.50Е-010 2.00Е-010 1.50Е-010

Т,'С

Рисунок 4. Зависимость силы адгезии РА от температуры клея Т.

Список использованной литературы

1. Попов В.М., Иванов А.В. Интенсивная технология получения клееной древесины повышенной прочности // Вестник МГУЛ «Лесной вестник». 2007. № 4. С. 89-91.

2. Попов В.М., Иванов А.В., Шендриков М.А. Влияние магнитного поля на прочность клееной древесины // Вестник МГУЛ «Лесной вестник». 2008. № 6. С. 80-81.

3. Воронежцев Ю.И., Гольдаде В. А., Пинчук Л.С. Электрические и магнитные поля в технологии полимерных композитов. Мн.: Навука 1 тэхшка, 1990. 263 с.

4. Хеерман Д.В. Метод компьютерного эксперимента в теоретической физике. М.: Наука, 1990. 176 с.

5. Лагарьков А.Н., Сергеев В.М. Метод молекулярной динамики в статической физике // УФН. 1987. Т. 125. № 7. С. 409-448.

6. Скрышевский А.Ф. Структурный анализ жидкостей и аморфных тел. М.: Высшая школа, 1980. 328 с.

7. Самойлов В.Г., Батаронов И.Л., Рощупкин С.А. Использование метода сглаживания экспериментальных данных для решения обратных задач матеметической физики // Изв. РАН. Сер. физ. 1995. Т. 59. № 10. С. 103-107.