Совершенствование вязальных механизмов кругловязальных машин тема диссертации и автореферата по ВАК РФ 05.02.13, кандидат технических наук Анашкина, Елена Владимировна

• Анашкина, Елена Владимировна
• кандидат технических наук
• 1998

Оглавление диссертации кандидат технических наук Анашкина, Елена Владимировна

СОДЕРЖАНИЕ

ВВЕДЕНИЕ

1. ОБЗОР ЛИТЕРАТУРНЫХ ИСТОЧНИКОВ И ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЯ

1.1. Обзор конструкций вязального механизма кругловязальных машин

1.1.1. Наклонное перемещение платины

1.1.2. Две платины

1.1.3. Сложное движение платины

1.2. Обзор литературы по исследованиям динамики петлеобразующих органов

1.2.1. Работы по снижению динамических нагрузок путем изменения конструкции игл

1.2.2. Работы по снижению динамических нагрузок путем модернизации клиньев

1.2.3. Работы по исследованию условий взаимодействия игл с пазами игольницы и клиновой системой, с учетом технологических зазоров

1.3. Обзор литературы по синтезу кулачковых механизмов 33 ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЯ

2. ТЕХНОЛОГИЧЕСКИЕ РАСЧЕТЫ ВЯЗАЛЬНОГО МЕХАНИЗМА

2.1. Сравнительный анализ процессов петлеобразования однофонтурных кругловязальных машин с постоянной и переменной отбойными плоскостями

2.2. Расчет усилия полезного сопротивления в зависимости от

угла кулирования

2.3. Движение и изменение натяжения нити в процессе вязания

за один игольный шаг

ВЫВОДЫ

3. ДИНАМИЧЕСКИЕ И МАТЕМАТИЧЕСКИЕ МОДЕЛИ ДВИЖЕНИЯ ИГЛИПЛАТИН В ЗАМКАХ ВЯЗАЛЬНОГО МЕХАНИЗМА

3.1. Динамическая модель движения иглы в замке вязального механизма

3.2. Математическая модель движения иглы в замке

вязального механизма

3.3. Динамическая модель движения платины в замке вязального механизма

3.4. Математическая модель движения платины в замке вязального механизма

3.5. Выбор закона движения петлеобразующих органов в замках вязального механизма

3.5.1. Выбор закона движения иглы в процессе вязания

3.5.2. Выбор закона движения платины в процессе вязания

3.5.2.1.Выбор закона поступательного перемещения платины

3.5.2.2. Выбор закона качательного движения платины 117 ВЫВОДЫ

4. МОДЕЛИРОВАНИЕ ДИНАМИКИ ПЕТЛЕОБРАЗУЮЩИХ ОРГАНОВ В ЗАМКАХ ВЯЗАЛЬНОГО МЕХАНИЗМА И ОПТИМИЗАЦИЯ ПРОФИЛЕЙ КЛИНЬЕВ ИГОЛЬНОГО И ПЛАТИННОГО ЗАМКОВ 125 4.1. Алгоритм моделирования движения петлеобразующих

органов в замках вязального механизма

4.1.1. Алгоритм моделирования движения иглы

4.1.2. Алгоритм моделирования движения платины в замке вязального механизма с переменной отбойной плоскостью

4.2. Моделирование динамики иглы в замке вязального механизма

4.3. Регулирование глубины кулирования

4.4. Моделирование динамики платины в замке вязального механизма с переменной отбойной плоскостью

4.4.1. Моделирование вертикального движения платины

4.4.2. Моделирование качательного движения платины 158 ВЫВОДЫ 161 ЗАКЛЮЧЕНИЕ 164 СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ 167 ПРИЛОЖЕНИЯ

Введение диссертации (часть автореферата) на тему «Совершенствование вязальных механизмов кругловязальных машин»

ВВЕДЕНИЕ

Актуальность темы. Трикотажные изделия занимают одно из ведущих мест среди товаров широкого спроса, они отличаются разнообразием, удобством в носке, гигиеничностью и красивым внешним видом. Ассортимент трикотажных изделий очень широк - это белье и одежда, полотна медицинского и технического назначения, обивочные и др. Всегда актуальными являются вопросы расширения ассортимента трикотажных изделий, повышения их качества и объемов выпуска. В настоящее время доля трикотажных изделий в общем объеме продукции данного назначения составляет (40-50)%, в том числе 20 % приходится на кругловязальные изделия1.

Одной из проблем, стоящих перед производителями трикотажа является неудовлетворительное качество сырья. При переработке на кругловязаль-ных машинах нитей низкого и среднего качества существенно повышается их обрывность. Для решения указанных задач на однофонтурных кругловязаль-ных машинах применяется техника петлеобразования с переменной отбойной плоскостью, которая предполагает встречное движение игл и платин. Техника вязания с переменной отбойной плоскостью известна достаточно давно, но в последние 20 лет появились новые конструктивные решения. Тем не менее в отечественной и зарубежной литературе машины с переменной отбойной плоскостью изучены недостаточно. Поэтому актуален анализ техники вязания с переменной отбойной плоскостью. Техника вязания с переменной отбойной плоскостью требует усложнения конструкции вязального механизма, так как платина совершает сложное движение, замки имеют моноблочную структуру. Поэтому актуальными вопросами являются исследование движения игл и платин в замках с точки зрения динамики, а также синтез законов их движения. Данная задача актуальна и для вязальных механизмов

1 по данным Заседания специалистов трикотажной промышленности Германии. IFWS - Frühjahrstagung, Melliand Textilber. - 1997. - 78, № 7-8. - с. 505

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

Цели и задачи исследования. Целью диссертационной работы является выдача инженерных рекомендаций по выбору: рациональной циклограммы движения петлеобразующих органов в вязальном механизме с переменной отбойной плоскостью, профилей клиньев игольного и платинного замков на основе математического моделирования динамики игл и платин; исследование натяжения, движения и продольных колебаний нити в процессе вязания; разработка алгоритмического и программного обеспечения для указанных целей.

Методы исследования. При выполнении работы применялись методы математического моделирования. Теоретические исследования базировались на основных законах математики и механики.

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

Разработаны динамические и математические модели движения игл и платин в замках вязального механизма, учитывающие: неудерживающий характер связи между иглой (платиной) и клином, податливость иглы (платины) в зоне контакта с клином, зазор в паре игла (платина) - клин, а также нелинейность сил, действующих на иглу (платину) со стороны нити и паза игольницы.

Разработано алгоритмическое и программное обеспечение на языке СИ для : определения натяжения нити и силы полезного сопротивления, преодолеваемого иглой при кулировании; моделирования движения иглы в замке

вязального механизма; моделирования движения платины в замке вязального механизма с переменной отбойной плоскостью;

Выполнено моделирование движения петлеобразующих органов в замках вязального механизма. По результатам математического моделирования на ЭВМ выбран компромиссный вариант профиля игольного замка, путем сравнения различных вариантов профилей по минимальным значениям: максимальных ускорений игл, протяженности вязальной системы. Даны инженерные рекомендации для проектирования клиньев замков вязального механизма.

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

Материалы диссертации будут использованы АО "Вулкан" при модернизации вязального механизма машин типа ЮЖ. Разработанное программное обеспечение используется в учебном процессе кафедры проектирования машин текстильной и легкой промышленности СПГУТД при подготовке инженеров - механиков спец. 170700 и бакалавров направления 551800 по дисциплине "Математические модели узлов машин и механизмов", в процессе курсового и дипломного проектирования.

Апробация работы. Основные положения диссертационной работы были доложены на конференции по машинам и аппаратам текстильной и легкой промышленности, посвященной 60-летию механического факультета

СПГУТД (1998 г.), научных семинарах кафедры "Проектирование машин текстильной и легкой промышленности", научном семинаре АО "Вулкан".

Публикации. По теме диссертационной работы опубликовано 4 статьи и тезисы доклада.

Структура и объем работы. Диссертация состоит из введения, четырех глав и заключения. Содержание изложено на 190 страницах машинописного текста, включающего 40 рисунков, 16 таблиц, 5 приложений, список использованных источников.

1. ОБЗОР ЛИТЕРАТУРНЫХ ИСТОЧНИКОВ И ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЯ

1Л Обзор конструкций вязального механизма кругловязальных

машин

В последнее время дальнейшее развитие трикотажного оборудования идет за счет совершенствования техники петлеобразования. Уже появились машины, где процесс петлеобразования обеспечивается встречным движением игл и платан, то есть применением так называемой техники переменной отбойной плоскости. Принцип техники вязания с переменной отбойной плоскостью получил основное свое развитие в последние двадцать лет за счет новых конструктивных решений. Это связано с тем, что техника переменной отбойной плоскости предъявляет высокие требования к надежности функционирования вязальной системы, используемым конструкционным материалам и методам их обработки. Реализация в практику соответствующих конструктивных решений стала возможной с достижением необходимой точности в машиностроении, появлением новых методов упрочнения и обработки металлов, новых технологий и т.д.

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

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

ным. Большой угол наклона кулирной части замка вызывает значительный износ игл и клиньев, ограничивает возможности увеличения скорости работы машины и применения пряжи с низкой разрывной прочностью.

Во избежание этих недостатков пытаются использовать конструкции вязального механизма, в которых платины перемещаются не только горизонтально, но и навстречу иглам. Термин "переменная отбойная плоскость" впервые использовал И.С. Мильченко в монографии [2]. В зарубежной литературе [3, 4] такая техника петлеобразования получила название Ке1апк.

Для анализа существующих конструкций механизмов с встречным движением игл и платан, были рассмотрены патентные материалы ведущих стран мира: США, ФРГ, Франции, Японии, Великобритании, СССР, Российской Федерации и стран Восточной Европы. При этом была использована реферативная информация из журналов: "Изобретения стран мира "Изобретения в СССР и за рубежом ", описания изобретений к авторским свидетельствам и патентам, реферативный журнал "Легкая промышленность", бюллетень "Открытия. Изобретения" и т.д. Поиск был произведен за двадцать лет, с 1975 по 1995 года.

Анализируя приведенные в патентах способы конструктивных решений процесса петлеобразования с применением техники переменной отбойной плоскости, можно выделить следующие типы конструкций:

1) платина перемещается относительно иглы по наклонной направ-

_ о

ляющеи;

2) две платины, одна из которых совершает вертикальное движение, другая движется горизонтально;

3) платина совершает сложное (поступательное и качательное) движение относительно иглы;

Рассмотрим указанные способы более подробно на базе конкретных конструктивных решений.

1.1.1. Наклонное перемещение платины

Конструкция, представленная в патенте №3986371 (США) [5] содержит платины 1 (см. рис. 1.1. а), закрепленные на вращающемся кольце внешнего цилиндра, и иглы 2, движущиеся по вертикальному пазу 3 игольного цилиндра 4. На одном конце платан 1 имеются горизонтальные кромки 5, образующие отбойную плоскость и расположенные перпендикулярно к движению игл 2. На другом конце платины 1 расположена пятка 6. С помощью кольцевой пружины 7 поддерживается контакт между пяткой 6 и поверхностью 8 крышки 9. Внутренняя поверхность 9 имеет определенный профиль, по которому скользит пятка 6, таким образом чтобы платина 1 двигалась поступательно под углом относительно иглы 2.

В патенте №1393352 США [6], пятка платины 1 и толкателя 2 (направляющей платины) двигаются в горизонтальном пазу (рис. 1.1.6), за счет относительного движения платины 1 и переносного движения толкателя 2 отбойная плоскость, которую образуют бородки платин 1, перемещается вверх и вниз.

В заявке №0351935 Европейского патентного ведомства [7] иглы 1 (рис. 1.1. в) машины, опираются на вращающийся игольный цилиндр 2 и движутся по вертикали относительно оси вращения цилиндра. Платина 3 находится на вращающемся кольце 4, которое расположено снаружи игольного цилиндра 2, и движется поступательно под наклоном вниз в радиальном направлении. Таким образом отбойная плоскость перемещается вверх при опускании иглы 1, и вниз при ее подъеме.

В качестве преимуществ данного типа конструкций можно назвать следующие: отсутствие наложения сил, т.к. платина имеет только одно направление движения; благоприятные условия трения, т.к. игла и платина движутся между жесткими штегами;

Конструкции вязального механизма с наклонным перемещением

платан

Заключение диссертации по теме «Машины, агрегаты и процессы (по отраслям)», Анашкина, Елена Владимировна

ВЫВОДЫ

1. Разработаны алгоритм и программа моделирования движения иглы в замке вязального механизма.

2. Разработаны алгоритм и программа моделирования движения платины в замке вязального механизма с переменной отбойной плоскостью.

3. По результатам математического моделирования движения игл вязального механизма на ЭВМ выполнен сравнительный анализ различных профилей клиньев игольного замка, синтезированных в п. 3.5. В качестве наилучшего профиля выбран вариант с переходными кривыми в виде парабол. Критериями оценки являлись наименьшая величина динамических нагрузок и протяженность вязальной системы. Наименьшую протяженность вязальной системы дает профиль с криволинейными участками в виде отрезков косинусоиды, но ускорения иглы в этом случае в среднем на 65 % больше чем при движении иглы по клину с переходными кривыми в виде отрезков парабол и синусоид. Сравнение двух последних профилей между собой (см. табл. 4.5) показывает, что в качестве компромиссного варианта можно выбрать профиль с криволинейными участками в виде отрезков парабол.

4. Как показали результаты моделирования, отрыв иглы от профиля (при идеальной его поверхности) происходит лишь при переходе с заключающего клина на кулирный и с кулирного на заключающий, поэтому в месте стыковки клиньев требуется особо высокое качество исполнения.

5. Для выбранного профиля клина выполнена оценка влияния на динамику иглы приведенной жесткости материала иглы, силы трения в пазу игольницы и коэффициента рассеяния при взаимодействии иглы с каналом игольного замка. Увеличение приведенной жесткости особенно влияет на максимальное ускорение иглы в месте стыковки заключающего и кулирного клиньев - оно увеличивается « на 25 %. Изменение силы трения пропорционально влияет на максимальные ускорения иглы. Увеличение коэффициента демпфирования способствует затуханию колебаний и уменьшению максимальных ускорений.

6. Исследовано влияние на динамику иглы регулирования глубины ку-лирования путем изменения положения кулирного клина. Результаты моделирования говорят об увеличении зоны виброударного режима движения иглы в месте стыковки заключающего и кулирного клиньев и максимальной величины ускорений, имеющих место в указанной зоне.

7. Выполнено моделирование поступательного движения платин в вязальном механизме с переменной отбойной плоскостью, в котором переходные кривые клиньев имеют вид отрезков парабол. Показано, что сокращение протяженности вязальной системы в два раза (при тех же, что и для постоянной отбойной плоскости углах подъема и купирования), которое достигается за счет усложнения конструкции вязального механизма, отрицательно влияет на динамику платин - « 50 % времени платина движется в зазоре. В качестве компромиссного решения предлагается увеличить протяженность вязальной системы, исключив из профилей прямолинейные участки.

8. Оценено влияние на динамику вертикального движения платины приведенной жесткости материала платины и коэффициента демпфирования при движении платины в канале платинного замка. Увеличение приведенной жесткости материала платины в два раза вызывает рост ускорений платины на отдельных участках до 70 %. При увеличении коэффициента демпфирования колебания ускорений платины затухают и уменьшается количество отрывов от профиля платинного замка.

9. Выполнено моделирование качательного движения платин по клину с переходными кривыми в виде отрезков парабол.

ЗАКЛЮЧЕНИЕ

На основании проведенных в работе исследований получены следующие результаты.

1) Выполнен конструктивный анализ вязальных механизмов с переменной отбойной плоскостью, дана их классификация. Предложена методика построения циклограммы движения петлеобразующих органов из условия минимальности их вертикальных перемещений.

2) Сравнительный анализ зависимостей натяжения нити вязальных механизмов с постоянной и переменной отбойными плоскостями показал, что техника вязания с переменной отбойной плоскостью позволяет снизить натяжение нити более чем в семь раз, что достигается за счет уменьшения протяженности вязальной системы и количества купирующих, игл. Это дает возможность переработки пряжи с низкими прочностными характеристиками.

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

4) Разработаны динамическая и математическая модели движения игл (для многосистемной однофонтурной кругловязальной машины) и платин (за основу для случая переменной отбойной плоскости принята конструкция вязального механизма с вертикальным и качательным движениями платины, по патенту ФРГ №3311361) в замках вязального механизма. Полученные математические модели являются сугубо нелинейными. Нелинейность обусловлена наличием зазоров в замках и силами сухого трения, действующими в системе.

5) Исходя из технологических требований к процессу вязания, условий обеспечения плавности движения иглы, предложены законы движения игл и платин в замках вязального механизма, состоящие из шести участков, два из которых прямолинейные и сопрягаются друг с другом переходными кривыми. В качестве переходных кривых используются типовые законы движения: параболический, в виде кривых третьего - четвертого порядков, косинусои-дальный и синусоидальный. Получены аналитические выражения для указанных законов движения игл и платин.

6) Разработано алгоритмическое и программное обеспечения для моделирования на ЭВМ движения игл и платин по математическим моделям с использованием разложения решений в ряды Тейлора до членов шестого порядка малости. Программное обеспечение может использоваться для моделирования движения игл и платин в замках вязальных механизмов с постоянной и переменной отбойными плоскостями.

7) Выполнено моделирование движения иглы по клиньям различных профилей с целью выбора одного из предложенных законов движения, обеспечивающих минимальные протяженность вязальной системы и ускорения иглы. В качестве компромиссного решения предлагается профиль клиньев (закон движения) с криволинейными участками в виде отрезков парабол. Для данного профиля выполнена оценка влияния на динамику иглы силовых и упру-го-диссипативных характеристик вязального механизма. Увеличение приведенной жесткости особенно влияет на максимальное ускорение иглы в месте стыковки заключающего и кулирного клиньев - оно увеличивается » на 25 %.

8) Исследовано влияние на динамику иглы регулирования глубины кулиро-вания путем изменения положения кулирного клина. Результаты моделирования говорят об увеличении зоны виброударного режима движения иглы в месте стыковки заключающего и кулирного клиньев и максимальной величины ускорений, имеющих место в указанной зоне.

9) Выполнено моделирование вертикального и качательного движения платины в замке вязального механизма с переменной отбойной плоскостью. Показано, что сокращение протяженности вязальной системы в два раза (при тех же, что и для постоянной отбойной плоскости углах подъема и кулирования), отрицательно влияет на динамику платин - « 50 % времени платина движется в зазоре.

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

Список литературы диссертационного исследования кандидат технических наук Анашкина, Елена Владимировна, 1998 год

СПИСОК ИСПОЛЬЗОВАННЫХ источников

1. Гарбарук В.Н. Проектирование трикотажных машин.-Jl.: Машиностроение, Ленингр. отд-ние, 1980. - 472 с.

2. Мильченко И.С. Основы проектирования трикотажных машин.-М., Ростехиздат, 1962. - 226 с.

3. Stein W., Graband Н. Aspekten von Relativtechnick. - Wirkerei- und Strikerei - Technick, 1993, №4, s. 327-330.

4. Relanit und High - Tech - Elektronik in der Rundstrickerei. - Wirkerei- und Strikerei - Technick, 1992, №12, s. 1136-1139.

5. Патент, США, №3986371, Изобретения стран мира, № 6,1975.

6. Патент, США, №1393352, Изобретения стран мира, № 9, 1977.

7. Заявка, Европейское патентное ведомство, №0351935, Изобретения стран мира, №№ 9,10, 1990.

8. Заявка, ФРГ, №3735561, Изобретения стран мира, №11, 1988.

9. Патент, ФРГ, №3247767, 1985.

10. Патент, ГДР, №215594, 1984.

11. Масленников Е.С. Ликвидация ударного разрушения игл при повышении скорости чулочных автоматов,- Изв. вузов. Технология легкой промышленности, 1974, № 1, с. 113-120, №2, с. 146 - 150.

12. Баранов A.A., Масленников Е.С. . Исследование динамики механизма отбора одноцилиндровых чулочно-носочных автоматов,- Изв. вузов. Технология легкой промышленности. 1990, № 3, с. 110 -113.

13. Петров Е.И. Работа чулочного автомата на повышенных скоростях.-Легкая промышленность. 1958, № 12,

14. Петров Ю.И., Петров Е.И. Влияние ударных волн на разрушение крючков игл трикотажных машин. - Текстильная промышленность. 1960, №5, с. 39-43.

15. Петров Е.И. Новые формы игл для скоростных чулочных автоматов. -Текстильная промышленность. 1962, № 1, с. 69 -72 .

16. Гайдачук И.П. Исследование нагрузок в петлеобразующих системах однофонтурных кругловязальных машин. - Автореферат канд. диссертации, Киев, 1977, 22 с.

17. Wray G.R., Burns N.D. Cam - to - neede impact forces in weft-knitting. Part V - Part VIII. Journal the Textil Institute, 1978, № 8.

18. Wray G.R., Burns N.D. Cam - to - neede impact forces in weft-knitting. Part IX - Part XI. Journal the Textil Institute, 1978, № 10.

19. Малышев А.П., Галушко И.В. Влияние формы антиударных игл на распространение ударных импульсов. - Изв. вузов. Технология легкой промышленности. 1992, № 1, с. 103 - 105.

20. Волощенко В.П., Пипа Б.Ф., Шипуков С.Т. Эксплуатационная надежность машин трикотажного производства. - Киев, "Техн1ка", 1977, 136 с.

21. Кот А.Н., Щербань В.Ю., Хомяк О.Н. Эффективность использования в петлеобразующих системах вязальных машин клиньев с подвижными элементами. - Изв. вузов. Технология легкой промышленности. 1989, № 3, с. 110-116.

22. Соловьев Л.И., Гарбарук В.Н. Клинья трикотажных машин, спроектированные по составным законам движения. - Изв. вузов. Технология легкой промышленности. 1967, № 5, с. 111 -117.

23. Домитеев Ф.А. , Полухин В.П. Анализ динамики пары "игла - клин" трикотажных машин. - Тез. докл. Всерос. научно-техн. конф. "Соврем.

технол. текстил. пром-сти", (Текстиль-95). Москва, 28-29 нояб. 1995. - М., 1995,- с. 67 - 68.

24. Кот А.Н., Хомяк О.Н., Щербань В.Ю. Оптимизация суммарного бокового зазора между иголой (игловодом) и стенками игольного паза перчаточного автомата. - Изв. вузов. Технология легкой промышленности. 1992, № 2, с. 82 - 87, № 3, с. 86 - 90.

25. Молчанов К.И., Соломонов В.И. Ударное взаимодействие иглы в пазу игольницы. - Изв. вузов. Технология легкой промышленности. 1985, №2, с. 118-120.

26. Соломонов В.И. Исследование и усовершенствование эксплуатационных параметров вязальной системы скоростных кругловязальных машин. -Канд. дисс., ЛИТЛП, 1985,118 с.

27. Вульфсон И.И. К проблеме динамической устойчивости при изгибных колебаниях игл кругловязальных машин. - В кн.: Машины и технология текстильного и трикотажного производства. Вып. 5. Рига, РПИ, 1974, с. 116-122.

28. Преображенская М.В., Вульфсон И.И. К методике определения собственной частоты иглы кругловязальной машины при изгибных колебаниях. - В кн.: Машины и технология текстильного и трикотажного производства. Вып. 5. Рига, РПИ, 1974, с. 109 -115.

29. Кот А.Н., Хомяк О.Н. Оптимизация зазора между клином и поверхностью игольницы трикотажной машины. - Изв. вузов. Технология легкой промышленности. 1989, №5,с.115-119.

30. Кот А.Н., Хомяк О.Н. Исследование влияния зазоров на нагруженность элементов механизма вязания вязальной машины. - Изв. вузов. Технология легкой промышленности. 1992, № 5-6, с. 74 -78.

31. Полухин В.П., Домитеев Ф.А., Масленников Е.С. Анализ взаимодействия клина с иглой в чулочно - носочных автоматах . - Изв. вузов. Технология текстильной промышленности. 1995, № 6, с. 82 - 84.

32. Левитский Н.И. Кулачковые механизмы. М., "Машиностроение", 1964, 265 с.

33. Артоболевский И.И. Теория механизмов и машин. М., 1988, 640 с.

34. Колчин Н.И. Механика машин. Т. 1 . М.-Л., Машгиз, 1971, 560 с.

35. РотбартГ.А. Кулачковые механизмы. Л., Судпромгиз, 1960, 336 с.

36. Тир К.В. Комплексный расчет кулачковых механизмов. М., Машгиз, 1958.

37. Юдин В.А. , Петрокас Л.В. Теория механизмов и машин. М., "Высшая школа", 1967, 528 с.

38. Вейц В.Л. Динамика машинных агрегатов. Л., "Машиностроение", 1969, 370 с.

39. Вульфсон И.И. Динамические расчеты цикловых механизмов. Л., "Машиностроение", 1976, 328 с.

40. Вульфсон И.И., Коловский М.З. Нелинейные задачи динамики машин. Л., "Машиностроение", 1968, 281 с.

41. Кобринский А.Е. Механизмы с упругими связями. М.,"Наука", 1964, 390 с.

42. Кожевников С.Н. Динамика машин с упругими связями. Киев, АН УССР,1961, 160 с.

43. Рагульскис K.M. Механизмы на вибрирующем основании ( Вопросы динамики и устойчивости). Каунас, АН Лит. ССР, 1963, 232 с.

44. Вульфсон И.И. Колебания в машинах. - Спб.: СПГУТД, 1996. - 185 с.

45. Корн Г., Корн Т. Справочник по математике для научных работников и иженеров. М.,"Наука", 1968, 720 с.

46. Бабаков И.М. Теория колебаний. М., "Наука", 1968, 560 с.

47. Лурье А.И. Аналитическая механика. М., Физматгиз, 1961, 824 с.

48. Вибрации в технике: Справочник. В 6-ти т./Ред. совет: В.Н. Чело-мей (пред.). - Т.1. Колебания линейных систем/ Под. ред. В.В. Болотина. -М., 1978, 352 с.

49. Деч Г. Руководство к практическому применению преобразования Лапласа и Z- преобразования. М, 1971

50. Бахматов К.И, Анашкина Е. В. Влияние применения техники вязания с переменной отбойной плоскостью на число вязальных систем. Вестник СПГУТД. СПб: Изд. СПГУТД . 1996 , №1.

51. Мазин Л.С., Бахматов К.И, Анашкина Е. В. Математическая модель движения иглы по заключающему клину. Изв. вузов. Технология текстильной промышленности,- 1996, №6, с. 63-67.

52. Анашкина Е. В., Бахматов К.И, Мазин Л.С. Моделирование движения иглы по заключающему клину кругловязальной машины. Изв. вузов. Технология текстильной промышленности. - 1997, №4, с. 69 -72.

53. Анашкина Е. В. О постановке задач проектирования трикотажного оборудования для производства искусственного меха . Исследование и разработка механизмов машин и средств автоматизации в текстильной и легкой промышленности. Меж-вуз. сб. науч.тр./ СПб: Изд. СПГУТД, 1998, с.82-85.

54.Анашкина Е. В., Бахматов К.И, Мазин Л.С. Определение предельных значений натяжения и проскальзывания перерабатываемой нити с учетом ее деформации. Тезисы докладов конференции по машинам и аппаратам

текстильной и легкой промышленности, посвященной 60-летию механического факультета. СПб: Изд. СПГУТД, 1998, с.45

ПРОГРАММА ОПРЕДЕЛЕНИЯ НАТЯЖЕНИЯ НИТИ И СИЛЫ ПОЛЕЗНОГО СОПРОТИВЛЕНИЯ, ПРЕОДОЛЕВАЕМОГО ИГЛОЙ ПРИ

КУЛИРОВАНИИ

#include <math.h> #include <stdio.h> #define PI 3.1415926 float b,zl,ti;

main() {

FILE *st; int n,i,ii;

float v,qmax,qO,mj;

float betk,ll ,12,l,h,al ,a2,a3,alfl ,alf2,alf3,P; double fimction();

/*Ввод исходных данных*/

zl=2.23e-3; /*Глубина кулирования*/

ti= 1. 2е-3; /*Игольный шаг*/

q0=10e-2; /*Начальное натяжение нити*/

mj=0.2; /*Коэффициент трения между п/о органами*/ betk=55*PI/180; /*Угол кулирования*/

y=0.7; /*Скорость вращения цилиндра*/

b=tan(PI-betk); /*Угловой коэффициент прямолинейного уч-ка*/ 11=0;

12=-2*zl/b;

п=30; /*Количество точек счета*/ h=(12-ll)/n;

st=fopen("u55.dat","w");

for(l=ll;l<=12;l=l+h) {

if(l>=ll && l<(ll+ti)) {

al=l-12;

alfl =fimction(al); qmax=q0*exp(2*mj*alfl);

}

else if(l>=(H+ti) && l<(ll+2*ti)) {

al=l-12;

alfl =function(al);

a2=l-ti-12;

alf2=function(a2);

qmax=qO*exp(mj *(3 *alfl+3 *alf2));

}

else if(l>=(ll+2*ti) && 1<=12) {

al=l-12;

alfl=function(al);

a2=l-ti-12;

alG=function(a2);

a3=l-2*ti-12;

alf3=function(a3);

qmax=qO*exp(mj *(3 *alfl+4*alf2+3 *alf3));

}

P=qmax*sin(alfl )*(exp(2*mj *alfl)+1 )/exp(2*mj *alfl); fjprintf(st," %e\t%e\t%e\n" ,1, qmax, P);

}

fclose (st); }

double functional)

float 11; {

float ksi,hk,alfa; ksi=(b*b/(4*zl))*ll*ll; hk=zl-ksi; alfa=atan(2 *hk/ti); return(alfa);

}

РЕЗУЛЬТАТЫ СЧЕТА ПО ПРОГРАММЕ ОПЕДЕЛЕНИЯ НАТЯЖЕНИЯ НИТИ И УСИЛИЯ ПОЛЕЗНОГО СОПРОТИВЛЕНИЯ, ПРЕОДОЛЕВАЕМОГО ИГЛОЙ ПРИ КУЛИРОВАНИИ

Постоянная отбойная плоскость

1, м Р,Н

0.00000е+00 1.00000е-01 0.00000е+00

1.04098е-04 1.10032е-01 4.97192е-02

2.08195е-04 1.19566е-01 9.48579е-02

3.12293е-04 1.27882е-01 1.31451е-01

4.16390е-04 1.34796е-01 1.59440е-01

5.20488е-04 1.40429е-01 1.80444е-01

6.24585е-04 1.45004е-01 1.96249е-01

7.28683е-04 1.48737е-01 2.08294е-01

8.32780е-04 1.51810е-01 2.17625е-01

9.36878е-04 1.54363е-01 2.24974е-01

1.04098е-03 1.56504е-01 2.30855е-01

1.14507е-03 1.58316е-01 2.35629е-01

1.24917е-03 2.16637е-01 3.24633е-01

1.35327е-03 2.51233е-01 3.78460е-01

1.45737е-03 2.84432е-01 4.30241е-01

1.56146е-03 3.14089е-01 4.76660е-01

1.66556е-03 3.39570е-01 5.16685е-01

1.76966е-03 3.61102е-01 5.50617е-01

1.87376е-03 3.79222е-01 5.79250е-01

1.97785е-03 3.94503е-01 6.03448е-01

2.08195е-03 4.07448е-01 6.23977е-01

2.18605е-03 4.18470е-01 6.41473е-01

2.29015е-03 4.27898е-01 6.56445е-01

2.39424е-03 4.35993е-01 6.69299е-01

2.49834е-03 6.42869е-01 9.87397е-01

2.60244е-03 7.41842е-01 1.13989е+00

2.70653е-03 8.33847е-01 1.28169е+00

2.81063е-03 9.14784е-01 1.40645е+00

2.91473е-03 9.84010е-01 1.51314е+00

3.01883е-03 1.04252е+00 1.60328е+00

3.12292е-03 1.09181 е+00 1.67913е+00

Переменная отбойная плоскость

1, м

0.00000Е+00 5.20488Е-05 1.56146Е-04 2.08195Е-04 3.12293Е-04 3.64341Е-04 4.68439Е-04 5.20488Е-04 6.24585Е-04 6.76634Е-04 7.80731Е-04 8.32780Е-04 9.36878Е-04 9.88927Е-04 1.09302Е-03 1.14507Е-03 1.24917Е-03 1.30122Е-03 1.40532Е-03 1.45736Е-03 1.56146Е-03 1.61351Е-03 1.71761Е-03 1.76966Е-03 1.87375Е-03 1.92580Е-03 1.97785Е-03 2.02990Е-03 2.08195Е-03 2.13400Е-03

0,н

1.06451Е-01 1.06764Е-01 1.07421Е-01 1.07761Е-01 1.08455Е-01 1.08802Е-01 1.09478Е-01 1.09798Е-01 1.10379Е-01 1.10629Е-01 1.11024Е-01 1.11162Е-01 1.11301Е-01 1.11300Е-01 1.11155Е-01 1.11016Е-01 1.54066Е-01 1.54623Е-01 1.55424Е-01 1.55672Е-01 1.55869Е-01 1.55820Е-01 1.55424Е-01 1.55075Е-01 1.54066Е-01 1.53403Е-01 1.52632Е-01 1.51755Е-01 1.50773Е-01 1.49690Е-01

Р,Н 2.04640Е-02 2.22010Е-02 2.60280Е-02 2.81098Е-02 3.25620Е-02 3.48998Е-02 3.96802Е-02 4.20604Е-02 4.66127Е-02 4.87011Е-02 5.22896Е-02 5.37108Е-02 5.56471Е-02 5.61201Е-02 5.60368Е-02 5.54930Е-02 7.45105Е-02 7.29985Е-02 6.89586Е-02 6.65348Е-02 6.11463Е-02 5.82830Е-02 5.24395Е-02 4.95324Е-02 4.38881Е-02 4.11915Е-02 3.85981Е-02 3.61176Е-02 3.37566Е-02 3.15187Е-02

ПРОГРАММА МОДЕЛИРОВАНИЯ ДВИЖЕНИЯ ИГЛЫ В ЗАМКЕ ВЯЗАЛЬНОГО МЕХАНИЗМА

#include<stdio ,h> #include<math.h> double regk(),regz(); double c,p,ft,g,m,d,t60;

main() {

double regk(double*, double*, double*, double*, double*, double*,double*,double*,int*);

double regz(double*,double*,double*,double*,double*,int*); double zl ,z2,z3,z4,z5,tig,b,bb,v,del,psi,z,11,12,13,14,15,16,

z0,dz0,t0,dz,d2z,t,tl0,t20,t30,t40,t50,dt,tl,ksi,pi,q,dksi,d2ksi,om;

int nj;

extern double c,p,ft,g,m,d,t60; FILE *kp, *kpl, *kp2, *kp3, *kp4, *kp5,*kp6; kp=fopen("pz.dat"V); kpl=fopen("pdz.dat","w"); kp2=fopen("pd2z.dat","w"); kp3=fopen("pksi.dat","w"); kp4=fopen("pdksi.dat","w"); kp5=fopen("pd2ksi.dat","w"); kp6=fopen("time","w"); z0=0.; dz0=0.; del=0.24e-3; zl=2.23e-3; z2=7.e-3; z3=14.04e-3;

z4=7.99e-3; z5=2.23e-3; pi=3.1415926; tig=1.2e-3; b=tan(50*pi/l 80); bb=tan(pi-55*pi/l 80);

v=.7; c=le4; p=4.e-2;

ft=24.e-2; g=6.987e-3; m=.713e-3;

psi=0.4; om=sqrt(c/m); d=c*psi/(2*pi*om); ll=2.*zl/b; tl0=ll/v;

12=(2.*zl+z2)/b; t20=12/v;

13=(2.*z3-z2)/b; t30=13/v;

14=13-2. *z4/bb; t40=14/v;

15=14+(-z3+z4+z5)/bb; t50=15/v; 16=15-2 *z5/bb; t60=16/v;

t0=t60; dt=.001*t0; z=z0; dz=dz0; d2z=0.;

ksi=0.; j=0; n=0;

/*Выбор вариантов режимов работы системы*/

for(t=0.; t<tO; t=t) {

if(t<tlO) {

ksi=b/2./ll*v*v*t*t; dksi=b*v*v/ll *t;

d2ksi=b*v*v/ll; }

else if(t<t20) {

ksi=b*(v*t-.5*ll); dksi=b*v; d2ksi=0.;

/

else if(t<t30) {

ksi=z3-b/2./(13-12)*(13-v*t)*(13-v*t); dksi=b*v/(13-12)*(13-v*t);

d2ksi=-b*v*v/(13-12); >

/

else if(t<t40) {

ksi=z3+bb/2./(14-13)*(v*t-13)*(v*t-13); dksi=bb*v/(14-13)*(v*t-13);

d2ksi=bb*v*v/(14-13); }

else if(t<t50) {

ksi=bb*(v*t-14)+z3-z4; dksi=bb*v; d2ksi=0.; }

else if(t<t60) {

ksi=bb/2./(15-16)*(v*t-16)*(v*t-16); dksi=bb*v/(15-16)*(v*t-16);

d2ksi=bb*v*v/(15-16); >

/

else if(t>=t60) {

ksi=ksi; dksi=0.; d2ksi=0.; }

fpnntf(kp," %f %f режим %d%d\n", t,z,j,n);

fprintf(kp1 ,"%f %f режим %d%d\n", t,dz,j,n); fprintf(kp2,"%f %f режим %d%d\n", t,d2z,j,n); iprintf(kp3,"%f %f режим %d%d\n", t,ksij,n); fprintf(kp4,"%f %f режим %d%d\n", t,dksi,j,n);

%f режим %d%d\n", t,d2ksi,j,n); tl=t; t=t+dt;

if(z-ksi>0 && z-ksi-del<0. && dz>0.) {

j=3; n=l; regz(&z,&dz,&t,&tl,&d2z,&n); continue; }

else if(z-ksi>0. && z-ksi-del<0. && dz<=0.) {

j=3; n=2; regz(&z,&dz,&t,&tl ,&d2z,&n); continue; }

else if(z-ksi<=0. && dz>0.) {

j=l; n=l; regk(&z,&dz,&t,&tl,&ksi,&dksi,&d2ksi,&d2z,&n);

continue; }

else if(z-ksi<=0. && dz<=0.) {

j=l; n=2; regk(&z,&dz,&t,&tl,&ksi,&dksi,&d2ksi,&d2z,&n);

continue; }

else if(z-ksi-del>=0. && dz>0.) {

j=2; n=l; ksi=ksi+del;

regk(&z,&dz,&t,&tl ,&ksi,&dksi,&d2ksi ,&d2z,&n); continue; }

else if(z-ksi-del>=0. && dz<=0.) {

j=2; n=2; ksi=ksi+del;

regk(&z,&dz ,&t ,&tl ,&ksi ,&dksi ,&d2ksi,&d2z,&n); continue; }

}

printf("nporpaMMa завершена\п"); priiltfj^"

return 0; }

/* Функция для анализа движения в контакте*/ double regk(y,dy,u,ul ,ks,dks,d2ks,d2y,nu)

double *y,*dy,*u,*ul,*ks,*dks,*d2ks,*d2y; int *nu; {

int n;

double z,dz,d2z,d3z,d4z,d5z,t,tl,dt,ksi, dksi,d2ksi,sf; extern double c,m,g,ft,p,d,t60; z=*y; dz=*dy; d2z=*d2y; t=*u; tl=*ul; ksi=*ks; dksi=*dks; d2ksi=*d2ks; n=*nu; if(n==l II t>=t60) sf=-g-ft-p; eise if(n==2)

sf=-g+ft+p; dt=t-tl;

d2z=l ,/m*(-c*z-d*dz+c*ksi+sf);

d3z=l Vm*(-c*dz-d*d2z+c*dksi);

d4z=l./m*(-c*d2z-d*d3z+c*d2ksi);

d5z=l ,/m*(-c*d3z-d*d4z);

z=z+dz*dt+.5 *d2z*dt*dt+1 ./6. *d3z*dt*dt*dt+

1 ,/24.*d4z*dt*dt*dt*dt+l ,/120.*d5z*dt*dt*dt*dt*dt;

dz=dz+d2z*dt+. 5 *d3z*dt*dt+1 ./6. *d4z*dt*dt*dt+

l./24.*d5z*dt*dt*dt*dt;

d2z=d2z+d3z*dt+.5*d4z*dt*dt+l./6.*d5z*dt*dt*dt; *y=z; *dy=dz; *d2y=d2z;

return 0; }

/* Функция для анализа движения в зазоре*/ double regz(y,dy,u,ul,d2y,nu)

double *y,*dy,*u,*ul,*d2y; int *nu; {

int n;

double z,dz,d2z,d3z,d4z,d5z,t,tl ,dt,sf; extern double m,g,ft,p,t60; z=*y; dz=*dy; d2z=*d2y; t=*u; tl=*ul;n=*nu; if(n==l || t>=t60) sf=-g-ft-p; else if(n==2)

s^-g+ft+p; dt=t-tl; d2z=l./m*sf; z=z+dz*dt+.5 *d2z*dt*dt; dz=dz+d2z*dt; d2z=d2z;

*y=z; *dy=dz; *d2y=d2z; return 0; }

ОБОЗНАЧЕНИЯ ПЕРЕМЕННЫХ В ПРОГРАММЕ МОДЕЛИРОВАНИЯ

ДВИЖЕНИЯ ИГЛЫ

Таблица 4.1.

Наименование Обоз- Машинные Численные

величин начение переменные значения

1 2 3 4

Исходные данные:

Начальные условия: 2(0)

перемещение иглы, м скорость иглы, м/с гО сЫ) 0 0

Максимальная величина зазора

в паре "игла-клин", м М 0,24 • Ю-3

Перемещение иглы из нижнего положения до уровня отбойной г! 2,23- Ю-3

плоскости (о.п.), м

Величина подъема иглы от о.п. до

момента схода петли с открытого 22 т1 7,00-10~3

язычка, м

Максимальное

перемещение иглы, м 23 гЗ 14,04 • 10~3

Опускание иглы до момента

прессования, м 24 г4 7,99-10"3

Глубина купирования, м ъЪ 2,23-10~3

Угол подъема, град. рб1р 50

Угол купирования, град. Рк 55

Линейная скорость цилиндра машины, м/с v v 0.7

Приведенная жесткость материала иглы , Н/м с с 104

Масса иглы, кг т т 0.713-Ю-3

Сила тяжести, Н О 0,698-10"2

Амплитудное значение силы по-

лезного сопротивления, Н ро Р 4,0 • 10~2

Амплитудное значение силы трения, Н а" йг 24,0-10-2

1 2 3 4

Остальные переменные:

Собственная частота, с"1 со от

Коэффициент рассеяния V <1

Коэффициенты в уравнении прямолинейного участка: для подъемного клина для опускающего клина Ъ1 ь2 Ь ЬЬ

Протяженности каждого из участков , м : подъемного клина опускающего клина и Лб 11,12,13 14,15,16

Время конца каждого из участков игольного замка, с П 0,120,130 140,150,160

Шаг счета, с Л

Обобщенная координата и ее производные, м , м/с, м/с2 х, йг, 62ъ

Функция положения клина и ее производные, м , м/с, м/с2 кв^сНодсШсз!

Переменные, определяющие режим работы иглы

ПРОГРАММА МОДЕЛИРОВАНИЯ КАЧАТЕЛЬНОГО ДВИЖЕНИЯ ПЛАТИНЫ В ЗАМКЕ ВЯЗАЛЬНОГО МЕХАНИЗМА С ПЕРЕМЕННОЙ

ОТБОЙНОЙ ПЛОСКОСТЬЮ

#include<stdio ,h> #include<math.h> #defme pi 3.1415926 double kont_bc(),zazor(); double c,lb,lc,im,mtr,g,d,t40;

main() {

double kont_bc(double*,double*,double*,double*,double*, double*,double*,double*,int*);

double zazor(double*,double*,double*,double *,double * ,int*);

extern double c,lb,lc,im,mtr,g,d,t40;

double zl,z2,z3,z4,z5,al,a2,v,delh,m,L,I,psi,om,bl,b2,b3,b4;

double tl0,t20,t30;

double zO,dzO,d2zO,z,dz,d2z,dt,t,tl;

double ksil ,dksil ,d2ksil ,ksi2,dksi2,d2ksi2;

double teta,fi,sl,s2,kb,sl2,s22,kc,s;

double zpk,hik,hpl ,zp 1 ,zp2,zp3,zp4,zp5;

double 111,112,113,114,115,116,tl 1 ,t21 ,t31 ,t41 ,t51;

double k,p,gl ,g2,g3,g4,kl ,pl ,bl2,b22,hpl,wpl;

intn,j;

FILE *kp, *kpl, *kp2, *kp3, *kp4, *kp5,*kp6,*kp7;

kp=fopen("z.dat","w");

kpl=fopen("dz.dat","w");

kp2=fopen("d2z.dat,,,"w");

kp3=fopen("ksi.dat","w");

kp4=fopen("dksi.dat",,,w");

kp5=fopen("d2ksi.dat","w");

kp6=fopen("time'V'w");

kp7=fopen("rasp.dat","w");

/*Ввод исходных данных*/ zl=2.23e-3; /*Глубина кулирования ,м*/ z2=7.00e-3; /*Перемещение до чаши ,м*/ z3=14.04e-3; /*Общее перемещение иглы ,м*/

al=l 1.82е-3; /*Перемещение иглы на заключение, м*/ а2=3.8 5е-3; /*Перемещение на прессование, м */ /*Радиальные перемещения платин, м*/

з1=3.528е-3; з2=0.209е-3; кЬ=2.9е-3*24.5/37.5; 812=-2.102е-3; 822=-0.135е-3; кс=-1.129е-3; 1е1а=55 *(р1/180); /*Угол подъема платины*/

й=50*(р1/180); /*Угол опускания платины*/ у=0.7; /*Линейная скорость цилиндра,м/с*/

с=1е+3; /*Приведенная жесткость системы платина-клин В, Н/м*/ р81=0.2; /*Коэффициент рассеяния*/ ёеШ=-0.1 е-3; /*Зазор при геометрическом замыкании,м*/ ш=0.913 е-3; /*Масса платины,кг*/ g=9.8*m; /*Сила тяжести ,Н*/ пйг=0.24е-3; /*Момент трения, Нм*/

/*Плечи сил действующих на платину при качательном движении*/ 1Ь=24.5е-3; 1с=-14.6е-3; Ь=8.25е-3;

/*Геометрические размеры платины - ширина и высота*/ лур1=7.е-3; Ьр1=55.е-3;

/* Предварительный технологический расчет: */ лрк=ъ\12\, /* Перемещение платины в процессе кулирования,м*/ Ык=г1/2; /*Глубина кулирования для переменной отбойной плоскости,м*/ Ьр1=а1/2; /*Перемещение платины в процессе заключения^*/

/*Расчет траектории движения платин*/ /*Углы наклона прямолинейных участков клина А*/ Ь3=1ап(2*р1-й); Ь4=1ап(1е1а);

/*Параметры траектории вертикального движения платин*/

гр1=:грк;

7ф2=ъ212\

2рЗ=Ьр1+2рк;

2р4=Ир1-Ык-а2/2;

гр5=грк;

/*Протяженность отдельных участков клина

вертикального движения платин,м*/ 111=(-2*гр1)/ЬЗ; 112=111 -гр2/ЬЗ; 113=112-2 -гр2)/ЬЗ;

114=113+2 * гр4/Ь4; 115=114+(грЗ-2р4-гр5)/Ь4;

116=115+2 *zp5/b4;

/*Параметры траектории клина В*/ s=sl+s2; k=s2+kb; p=kb;

/^Расчет протяженности и углов наклона клина В*/ g3=(s-p) * (s *113 -k*116)/(s *(p-k)) ; g4=116-113-g3; gl=:k*g4/(s-k); g2=p*g3/(s-p); bl=2*k/gl; b2=-2*p/g2; tlO=gl/v; t20=(gl+g2)/v; t30=(gl+g2+g3)/v; t40=(gl+g2+g3+g4)/v; /*Параметры траектории клина С*/ kl=s22+kc; pl=kc; bl2=2*kl/gl; b22=-2*pl/g2;

I=m* (wpl* wpl+hpl*hpl)/12 ; im=I+m*L*L; om=sqrt(c/I); d=c*psi/(2*pi*om);

/* Начальные условия */ z0=(-0. 5 *b 1 *gl+s2)/lb; dzO=bl*v/lb; d2z0=(-bl*v*v/gl)/lb; z=zO; dz=dzO; d2z=d2z0;

j=0; n=0; dt=.001*t40;

/*Выбор вариантов режимов работы системы*/

for(t=0.; t<=t40; t=t) {

/*ksil - профиль клина В; ksi2- профиль клина С*/

if(t<tlO) {

ksil=-0.5*bl/gl*(v*t-gl)*(v*t-gl)+s2;

dksil=-bl/gl*(v*t-gl)*v;

d2ksi 1=-b 1 * v* v/g 1 ;

ksi2=-0.5*bl2/gl*(v*t-gl)*(v*t-gl)+s22; dksi2=-bl2/gl *(v*t-gl)*v;

d2ksi2=-b 12 * v* v/g 1 ; }

else if(t<t20) {

ksil=0.5*b2/g2*(v*t-gl)*(v*t-gl)+s2;

dksil=b2/g2*(v*t-gl)*v;

d2ksi 1 =Ь2 * v* v/g2;

ksi2=0.5 *b22/g2 * (v*t-g 1) * (v*t-g 1 )+s22; dksi2=b22/g2*v*(v*t-gl);

d2ksi2=b22*v*v/g2; }

else if(t<t30) {

ksil=-(0.5*b2/g3)*(v*t-(gl+g2+g3))*(v*t-(gl+g2+g3))-sl;

dksil=-b2/g3*v*(v*t-(gl+g2+g3));

d2ksil=-b2*v*v/g3;

ksi2=-(0.5*b22/g3)*(v*t-(gl+g2+g3))*(v*t-(gl+g2+g3))-sl2; dksi2=-b22/g3*(v*t-(gl+g2+g3));

d2ksi2=-b22*v*v/g3; }

else if(t<-t40) {

ksil=(0.5*bl/g4)*(v*t-(gl+g2+g3))*(v*t-(gl+g2+g3))-sl;

dksil=bl/g4*v*(v*t-(gl+g2+g3));

d2ksil=bl *v*v/g4;

ksi2=(0.5*bl2/g4)*(v*t-(gl+g2+g3))*(v*t-(gl+g2+g3))-sl2; dksi2=b 12/g4 * v*(v*t-(g 1+g2+g3));

d2ksi2=b 12 * v* v/g4; }

fprmtf(kp,'' %f %f режим %d%d\n", t,z*lb,j,n); fprmtf(kpl ,"%f %f режим %d%d\n", t,dz*lb,j,n); fprintf(kp2,M%f %f режим %d%d\n", t,d2z*lbj,n); ¿rintf(kp3,"%f %f %f режим %d%d\n", t,ksil,ksi2,j,n); iprintf(kp4,"%f %f режим %d%d\n", t,dksilj,n); fprintf(kp5,"%f %f режим %d%d\n", t,d2ksil,j,n); tl=t; t=t+dt;

if(-lb*z+ksil<=0 && -lb*z+ksil>=delh && dz>0.) {

j=3; n=l; zazor(&z,&dz,&t,&tl,&d2z,&n); continue; }

else if(-lb*z+ksil<=0 && -lb*z+ksil>=delh && dz<=0.) {

j=3; n=2; zazor(&z,&dz,&t,&tl ,&d2z,&n); continue; }

else if(-lb*z+ksi 1 <delh && dz>0.)

{

j=l; n=l; ksil=ksil+delh;

kont_bc(&z ,&dz ,&t ,&tl ,&ksil ,&dksil ,&d2ksil ,&d2z,&n); }

else if(-lb*z+ksi 1 <delh && dz<=0.) {

j=l; n=2; ksil=ksil+delh;

kont_bc(&z,&dz ,&t,&tl ,&ksil ,&dksil ,&d2ksil ,&d2z,&n); }

else if(-lb*z+ksil>0. && dz>0.) {

j=2; n-l;

kont_bc(&z ,&dz ,&t,&tl ,&ksil ,&dksil,&d2ksil ,&d2z,&n); }

else if(-lb*z+ksil>0. && dz<=0.) {

j=2; n=2;

kont_bc(&z,&dz,&t,&tl ,&ksil ,&dksil ,&d2ksil ,&d2z,&n); }

}

printf("nporpaMMa завершена\п");

return 0; }

/* Функция для анализа движения в контакте*/ double kont_bc(y,dy,u,ul ,ks,dks,d2ks,d2y,nu)

double *y,*dy,*u,*ul,*ks,*dks,*d2ks,*d2y; int *nu; {

intn;

extern double c,lb,lc,im,mtr,g,d,t40; double z,dz,d2z,d3z,d4z,d5z,t,tl,dt,ksi, dksi,d2ksi,mf; z=*y; dz=*dy; d2z=*d2y; t=*u; tl=*ul; ksi=*ks; dksi=:*dks; d2ksi=*d2ks; n=*nu; if(n=l)

mf=mtr; else if(n==2)

mf=-mtr; dt=t-t 1;

d2z=l./im*(c*lb*(-lb*z+ksi)-d*dz+mf); d3z=l./im*(c*lb*(-lb*dz+dksi)-d*d2z); d4z=l./im*(c*lb*(-lb*d2z+d2ksi)-d*d3z);

continue;

continue;

continue;

continue;

d5z=l ,/im*(-c*lb*lb*d3z-d*d4z); z=z+dz*dt+.5*d2z*dt*dt+l ,/6.*d3z*dt*dt*dt+ 1 ,/24.*d4z*dt*dt*dt*dt+l ,/120.*d5z*dt*dt*dt*dt*dt; dz=dz+d2z*dt+.5*d3z*dt*dt+l./6.*d4z*dt*dt*dt+ l./24.*d5z*dt*dt*dt*dt;

d2z=d2z+d3z*dt+.5*d4z*dt*dt+l./6.*d5z*dt*dt*dt; *y=z; *dy=dz; *d2y=d2z;

return 0; }

/* Функция для анализа движения в зазоре*/ double zazor(y,dy,u,ul,d2y,nu)

double *y,*dy,*u,*ul,*d2y; int *nu; {

intn;

extern double c,lb,lc,im,mtr,g,m,d,t40; double z,dz,d2z,d3z,d4z,d5z,t,tl ,dt,mf; z=*y; dz=*dy; d2z=*d2y; t=*u; tl=*ul; n=*nu;

if(n==l)

mf=mtr;

else

if(n==2)

mf=-mtr;

dt=t-tl; d2z=l./im*mf; z=z+dz*dt+.5*d2z*dt*dt; dz=dz+d2z*dt; d2z=d2z;

*y=z; *dy=dz; *d2y=d2z;

return 0; }