Помощь в написании студенческих работ
Антистрессовый сервис

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

Диссертация: Конструкционная прочность соединения с гарантированным натягом по прерывистой цилиндрической поверхности
ДиссертацияПомощь в написанииУзнать стоимостьмоей работы

В машиностроительной фирме «НИПРО» в г. Орле для получения компоновочного преимущества в трехступенчатом цилиндрическом зубчатом редукторе применена посадка зубчатого колеса тихоходного промежуточного вала на частично прорезанную (прерывистую) цилиндрическую поверхность вал-шестерни последней ступени. В разделе 1.2 будет приведен и другой пример конструкции, в которой посадка зубчатого колеса… Читать ещё >

Конструкционная прочность соединения с гарантированным натягом по прерывистой цилиндрической поверхности (реферат, курсовая, диплом, контрольная)

Содержание

  • 1. Состояние проблемы и постановка задач исследования
    • 1. 1. Обзор литературы
    • 1. 2. Конструкции, приводящие к задаче определения напряженно-деформированного состояния соединения с гарантированным натягом по прерывистой цилиндрической поверхности
    • 1. 3. Постановка задачи
    • 1. 4. Основные допущения, принятые в исследованиях
    • 1. 5. Выводы по главе
  • 2. Теоретическое исследование конструкционной прочности соединения
    • 2. 1. Соединение зубчатого вала со ступицей и условия на их границах в граничных интегральных уравнениях (ГИУ)
    • 2. 2. Метод решения — метод граничных конечных элементов (МГКЭ)
    • 2. 3. Геометрические параметры элементов зубчатого вала и ступицы
      • 2. 3. 1. Симметрия соединения
      • 2. 3. 2. Отображение исходного контура на контур зубчатого венца
      • 2. 3. 3. Зубчатый венец с эвольвентным зацеплением
      • 2. 3. 4. Зубчатый венец с зацеплением Новикова
      • 2. 3. 5. Отверстие полого вала
      • 2. 3. 6. Отверстие ступицы
      • 2. 3. 7. Внешний контур ступицы
      • 2. 3. 8. Координаты середины элементов, их длины и направляющие косинусы касательной и внешней нормали
      • 2. 3. 9. Преобразование координат
    • 2. 4. Условия на границах соединения в дискретной форме и важнейшие параметры напряженно-деформированного состояния
      • 2. 4. 1. Пакет машинных программ «PRESS»
    • 2. 5. Результаты теоретического решения
      • 2. 5. 1. Критерии конструкционной прочности соединения — коэффициенты прочности С0, CF, См
      • 2. 5. 2. Основные параметры, влияющие на конструкционную прочность соединения
      • 2. 5. 3. Исследование конструкционной прочности соединения цилиндрического вала со ступицей по прерывистой поверхности
      • 2. 5. 4. Исследование конструкционной прочности соединения зубчатого вала со ступицей с гладким цилиндрическим отверстием
    • 2. 6. Выводы по главе
  • 3. Оценка точности теоретического решения
    • 3. 1. Сравнение результатов расчета с тестовыми задачами
      • 3. 1. 1. Соединение цилиндрического вала со ступицей (задача Ляме)
      • 3. 1. 2. Задача о равновесии части соединения
    • 3. 2. Экспериментальное исследование конструкционной прочности соединения
      • 3. 2. 1. Модели зубчатого вала и ступицы
      • 3. 2. 2. Экспериментальные установки
      • 3. 2. 3. Результаты эксперимента, их анализ и сравнение с результатами теоретического решения
    • 3. 3. Выводы по главе
  • 4. Метод инженерного расчета соединения с гарантированным натягом по прерывистой цилиндрической поверхности
    • 4. 1. Разработка формулы для инженерного расчета коэффициента прочности С
    • 4. 2. Оценка точности формулы инженерного расчета для коэффициента прочности С
    • 4. 3. Пример расчета соединения с натягом по прерывистой цилиндрической поверхности

    4.4 Возможные направления применения метода граничных элементов в задачах расчета конструкционной прочности соединения с гарантированным натягом с контактом по неполной (прерывистой) рабочей поверхности.

    4.5. Выводы по главе.

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

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

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

Из термических методов более распространен нагрев деталь-ступицы, осуществляемый в печах или индукционным нагревом. Преимуществом метода является высокая прочность соединения — в 2−3 раза выше, чем при запрессовке [17, 30] и независимость от наличия пресса. Недостатком является опасность изменения структуры материала детали, ограничение в выборе антикоррозионных покрытий термостойкими покрытиями, возможность появления при охлаждении торцевых зазоров и практическая неразборность соединения.

Развитие холодильной техники расширило возможность осуществления посадок с помощью охлаждения жидким азотом (—196°С) или твердой углекислотой (-78°С) детали-вала. Метод рационален в первую очередь в случаях, когда деталь-ступица во много раз тяжелее (или имеет большие размеры), чем деталь-вал. Получаемое соединение на 10−15% прочнее [22, 23, 24], чем при использовании нагрева, свободно от опасений нарушения структуры материала деталей соединения и не имеет ограничений в выборе антикоррозионных или повышающих сцепление покрытий. Соединение сохраняет недостаток метода нагрева — практическую неразъемность.

Все большее распространение получает метод гидрораспора цилиндрических или конических поверхностей соединения при сборке и разборке. Нагнетаемое в зону соединения масло должно иметь давление выше контактного давления в соединении. При этом устраняется прямой контакт между соединяемыми деталями и осевое усилие, требуемое для их относительного перемещения снижается во много раз (при распрессовке конических соединений до нуля). Хотя прочность соединения при этом несколько ниже, чем при термических методах, зато оно лишено всех отмеченных для других способов недостатков. Однако и оно не лишено недостатка — необходимости конструктивного усложнения деталей соединения и наличия спецоборудования в местах производства и ремонта.

Актуальность работы. При расчете контактных давлений в области упругих деформаций в соединениях с круговыми цилиндрическими поверхностями обычно используют решение Ляме [80] 5 с1.

2 Е.

Е^ ^ с1 — с/) у 1 й- + (I1 /?2 у — <1~ где р — контактное давлениес/ - номинальный диаметр соединенияс/| и 4 — внутренний диаметр пустотелой охватываемой детали и наружный диаметр охватывающей детали соответственно- 5 — натяг соединения — разность диаметров контакта охватываемой и охватывающей деталейнЕ и Е2, и ц.2 — модули упругости и коэффициенты Пуассона охватываемой и охватывающей деталей соответственно. Однако, в последнее время наряду с обычными соединениями нашли применение и более сложные, например: эксцентриковые [91, 92], соединения по поверхности многовершинного овала [10] и другие, для которых приведенная выше формула является приближенной, так как не учитывает специфику каждого из них. Для более точного расчета таких соединений применяются, как правило, численные методы. Однако во всех перечисленных нестандартных соединениях с натягом контакт происходит по сплошной цилиндрической поверхности.

В специальную группу можно выделить соединения с гарантированным натягом по макропрерывистой цилиндрической поверхности. Под макропрерывистостью следует понимать не неполный контакт сопрягаемых поверхностей за счет их шероховатости (микропрерывистость) [36, 52, 53, 56], а неполный контакт за счет геометрического выпадения отдельных участков контакта. К таким соединениям можно отнести комбинированные соединения с натягом (передающие нагрузки одновременно за счет натяга и стопорящих элементов, например, шпонок) [34, 35, 72, 73], соединения сплошного вала со втулкой, контактная поверхность которой расчленена кольцевыми пазами или пазами параллельными оси вала, при этом контактная поверхность втулки становится прерывистой [39]. Расчет напряженно деформированного состояния (НДС) соединений этой группы разработан слабо, поэтому формула Ляме в этих случаях неприемлема.

В машиностроительной фирме «НИПРО» в г. Орле для получения компоновочного преимущества в трехступенчатом цилиндрическом зубчатом редукторе применена посадка зубчатого колеса тихоходного промежуточного вала на частично прорезанную (прерывистую) цилиндрическую поверхность вал-шестерни последней ступени. В разделе 1.2 будет приведен и другой пример конструкции, в которой посадка зубчатого колеса на частично или полностью прорезанную поверхность вал-шестерни увеличивает жесткость конструкции. Особенностью соединения в обоих этих случаях является контакт по прерывистой цилиндрической поверхности — вершинам зубьев и свободной от контакта остальной поверхностью сложной геометрической формы — боковые поверхности зубьев, впадины. Исследование конструкционной прочности и разработка методов расчета НДС такого соединения представляет значительный научный и практический интерес и составляет предмет настоящей работы.

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

Научная новизна.

Предложена математическая модель определения напряженно-деформированного состояния соединения с гарантированным натягом с прерывистым цилиндрическим контактом.

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

Разработан пакет программ «PRESS».

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

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

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

На защиту вынесены следующие положения:

1. Исследование напряженно-деформированного состояния нового вида соединения с натягом по прерывистой цилиндрической поверхности:

— выбор и обоснование расчетной модели соединения;

— выбор и обоснование метода теоретического решения;

— разработка теоретического решения- 9.

— разработка комплекса машинных программ «PRESS» для реализации теоретического решения;

— исследование влияния геометрических параметров соединения на его конструкционную прочность;

— обоснование достаточности нагрузочной способности соединения;

— оценка точности разработанного метода на примере решения задач-эталонов;

— экспериментальная оценка точности расчета конструкционной прочности соединения.

2. Рекомендации по применению соединения и методики расчета его на конструкционную прочность:

— точная (компьютерная), основанная на использовании разработанных в работе машинных программ;

— инженерная, основанная на аналитическом и графическом описании результатов теоретического исследования.

9. Результаты работы — пакет программ «PRESS» и инженерный метод расчета — внедрены в производство машиностроительной фирмой ООО «НИПРО» (г. Орел), поставившей в эксплуатацию на металлургические комбинаты России («ММК», «НЛМК», «НТМК» и «Северсталь») 20 тяжелых крановых редукторов ВКУ-965П. В редукторах произведена бесшпоночная посадка зубчатого колеса на центральную часть шевронной вал-шестерни, прорезанную впадинами зубьев. Результаты работы приняты к использованию на Майкопском и Ижевском редукторных заводах.

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

Показать весь текст

Список литературы

  1. И.В., Турыгин Ю. В. МКЭ для расчета НДС соединений с ав-тофретированными охватывающими деталями // Проблемы прочности. — Киев: Наукова думка, 1987. — № 3. — С. 105−108.
  2. И.В., Турыгин Ю. В., Клековкин B.C. Нагрузочная способность соединений с автофретированными деталями // Вестник машиностроения. 1983. — № 4. — С. 28−30.
  3. Ю.А. Об одной задаче упругого равновесия сопряженных посредствам посадки кусочно-однородных сред // Изв. АН СССР. Механизация и машиностроение. 1963. — № 5. — С. 179−186.
  4. Ю.А. К определению напряженного состояния в эксцентрическом кольце сопряженной посредством натяга круглой шайбой из другого материала//Изв. АН СССР. Механика. -1965. № 1. — С. 142−145.
  5. Ю.А. Упругое равновесие круглой пластинки с эллиптическим отверстием, в которое посредством натяга вставлена шайба из другого материала / Инж. журнал. МТТ. — 1966. — № 6. — С. 122−129.
  6. Ю.А. К проблеме напряженной посадки кусочно-однородных сред / В кн.: Контактные задачи и их инженерные приложения. -М.: Издательство НИИмаш, 1969 С. 111−122.
  7. Г. Я. Тепловая сборка колесных пар. Издательство Харьковский университет, 1965.
  8. П.К. Определение напряжений в круглой пластинке после запрессовки в нее диска квадратной формы // Известия вузов. Машиностроение. 1964. — № 7. — С.5−10.
  9. П.К. Напряженное состояние круглой пластинки с запрессованным в нее скругленным п угольным диском / Сб.: Расчеты на прочность. -М.: Машиностроение, 1969. — Вып.14. — С. 175−188.
  10. П.К. К задаче о прессовой посадке втулки на некруглый вал. Механическое оборудование металлообрабатывающих заводов / Сб. научных трудов ЖМИ. Высшая школа, 1971. -Вып.16. — С. 355−366.
  11. Ахмедов А.Д. О. Некоторые плоские контактные задачи для кусочно-однородных деталей различной конфигурации. Автореферат диссертации на соискание ученой степени к.т.н. АН Каз.ССР. Институт математики и механики. Алма-Ата, 1987. — 20 с.
  12. Л.Т. Усталость валов в соединениях. Киев: Техника, 1982.- 151 с.
  13. И.А. Об одной задаче определений поля напряжений в составной двусвязной пластинке / В кн.: Проблемы расчета пространственных конструкций. -М., 1980. С. 3−10.
  14. И.А. Некоторые задачи теории упругости для однородных и кусочно-однородных сред. Диссертация на соискание ученой степени д. ф-м.н. Баку, 1973.-381 с.
  15. И.А., Вольпе Г. С. К определению напряженного состояния в деталях с различными упругими константами возникающего при запрессовке // Изв. АН Азерб. ССР. 1969. — № 4. — С. 42−49.
  16. H.A., Караев А. Б. К вопросу определения поля напряжений в составных деталях, изготовленных из трех разных однородных материалов, соединенных между собой с натягом // Аз НИИНТИ. 1986. — № 727. — 9 с.
  17. Е.Ф. Расчет и выбор посадок с натягом из системы ИСО. / -Пенза: Издательство Саратовский университет, 1975. — 187с.
  18. Д.В. Дополнительные главы линейной алгебры. М.: Наука, 1983.-336 с.
  19. Н.М. Сопротивление материалов. -М.: ГИТТЛ, 1964. 583 с.
  20. П., Баттерфильд Р. Метод граничных элементов в прикладных науках. М.: Мир, 1984. — 494 с.
  21. E.H. Посадки с натягом в машиностроении. М.: Машиностроение, 1966. — 188 с.
  22. Г. А. Прочность посадок, осуществляемых с применением холода. М.: Машиностроение, 1971. — 95 с.
  23. Г. А. Сборка в машиностроении с применением холода. -М.: Машгиз, 1959.
  24. Г. А., Михайленко Л. Ф., Зенкин A.C., Беляев Н. С. Влияние технологических факторов на прочность соединений, осуществляемых с помощью глубокого холода / Технология и организация производства. -ИССН, 1974.-№ 5.-С. 16−18.
  25. К., Уокер С. Применение метода граничных элементов в технике. М.: Мир, 1982. — 246 с.
  26. И.Н. Приложение метода акад. Н. И. Мусхелишвили к решению граничных задач плоской теории упругости анизотропной среды / Сообщение Груз, филиала АН СССР. Т. 1.-№ 10. — 1940. — С. 10−38.
  27. И.Н., Мусхелишвили Н. И. Методы теории аналитических функций в теории упругости / Труды Всесоюзного съезда по теоретической и прикладной механике. М. -Л.: Издательство АН СССР, 1962. — С. 310−338.
  28. Ф., Миррей У., Райт М. Практическая оптимизация. М.: Мир, 1985.-510 с. ^
  29. И.С., Ильяшенко A.A. Соединения с натягом. Расчет, проектирование, изготовление. М.: Машиностроение, 1981. — 274 с.
  30. Е.С., Кречетов И. Э. Исследование прочности соединения оси с колесом тепловозных колесных пар / Труды ВНИТИ. ИССН, 1970. -Вып. 34.-С. 216−224.
  31. Н.Б. Контактирование шероховатых поверхностей. М.: Наука, 1970.
  32. Н.Б. Фактическая площадь касания твердых поверхностей. -М.: Издательство АН СССР, 1962.
  33. А., Лю Дж. Численное решение больших разряженных систем уравнений. -М.: Мир, 1984.-334 с.
  34. М.М. и др. К вопросу о жесткости соединения вала и втулки призматической шпонкой / В сб.: Научные труды Ленинградского института текстильной и легкой промышленности. 1971. -№ 11.
  35. М.М. Приближенный метод расчета деформаций штив-тового соединения вала и втулки //Известия вузов. Приборостроение. 1968. -№ 1.
  36. В.А. К вопросу о теоретическом обосновании закона Амонтана Кулона для трения несмазанных поверхностей / Ж.Т.Ф. — 1940. -Вып. 17. — № 10.
  37. А.Н., Медведева Р. В., Патрикевич Ф. В. Конические соединения. / Справочное пособие. М.: Машиностроение, 1968. — 142 с.
  38. А.Л. Исследование и разработка методов технологического обеспечения эксплуатационных свойств посадок с натягом. Диссертация на соискание ученой степени к.т.н. Киев, 1980. — 160 с.
  39. О. Метод конечных элементов в технике. М.: Мир, 1975.544 с.
  40. И.А. Напряженное состояние квадратной пластинки с запрессованной круглой шайбой из другого материала. Издательство ВУЗов. Машиностроение, 1969. -№ 9. — С.48−53.
  41. И.А. Напряженное состояние пластинки, ограниченной кривой частного вида, с несколькими запрессованными в нее круглыми шайбами из другого материала. Издательство АН СССР. МТТ. — 1971. — № 1. — С. 89.
  42. И.А. О напряжениях в кусочно-однородной плоскости с многогранным отверстием / В кн.: Методы решения задач упругости и пластичности. Горький, 1974. — Вып.8. — С. 43−51.
  43. И.А. Решение плоской задачи теории упругости для кусочно-однородных сред, сопряженных посредством посадки. Диссертация на соискание ученой степени к.т.н. Пермь, 1969. — 106 с.
  44. Г. Б., Лукащук Ю. В. Влияние некоторых конструктивных факторов на распределение напряжений в соединениях с натягом // Вестниц машиностроения. 1980. — № 4. — С. 38−42.
  45. A.M., Эфендиев Б. Р. Упругое равновесие квадратной пластинки с круговым отверстием, в которое посредством натяга вставлена шайба из другого материала / Уч. записки АГУ им. С. М. Кирова. Серия физ.-мат. наук. -1967.-№ 4.-С. 72−79.
  46. Караев А.Б. О. Методика расчета составных ослабленных пластинок различной конфигурации, соединенных между собой с натягом. Автореферат диссертации на соискание ученой степени к.т.н. Азерб. политехнический институт. Баку, 1988. — 13 с.
  47. B.C. Конструкторско-технологические основы управления нагрузочной способностью соединений с натягом. Автореферат диссертации на соискание ученой степени д.т.н. Ижевск, 1995.
  48. B.C., Абрамов И. В. Исследование НДС деталей в соединениях с автофретированной охватывающей деталью / Сборник трудов: Динамика и прочность механических систем. Пермь: ППИ, 1976. — № 186. — С. 112−116.
  49. B.C., Абрамов И. В. Исследование статической прочности цилиндрических соединений с автофретированными охватывающими деталями, // Вестник машиностроения. 1977. — № 8. — С. 21−24. ^
  50. М.М. Прочность прессовых соединений при повторно-переменных нагрузках. М.: Машгиз, 1954. — 256 с. ^
  51. И.В. Трение несмазанных поверхностей. Автореферат диссертации на соискание ученой степени д.т.н. ИМАШ, 1943.
  52. И.В. Трение покоя двух шероховатых поверхностей // Изв. АН СССР. 1940. — № 10.
  53. Г. В. Об одном приложении теории функций комплексного переменного к плоской задаче математической теории упругости. Юрьев, 1909.- 187 с.
  54. Г. В. Применение комплексной переменной в теории упругости. М: Гостехиздат, 1935. — 224 с.
  55. Н.Е. Исследование величин фактической площади контактаи ее влияние на качество соединений с натягом. Автореферат диссертации на V соискание ученой степени к.т.н. М.: МВТУ, 1976. — 16 с.
  56. С.Г. Теория упругости анизотропного тела / 2-е изд. М.: Наука, 1977.-415 с.
  57. Ляв А. Математическая теория упругости. -М.: ОНТИ, 1935.
  58. Л.В. О деформации цилиндрических колец при напряженной посадке с учетом малой некруглости сопрягаемых поверхностей. «Изв. Вузов. Машиностроение», 1973, № 7
  59. Т.М. Расчеты прессовых посадок шайб за пределами упругости при степенном упрочнении материала // Расчеты на прочность. Машгиз.-1960. — Вып. 6.-С.113−116.
  60. С.Г. Вариационные методы в математической физике. М., 1957.-476 с.
  61. С.Г. Плоская задача теории упругости для неоднородной среды. М.-Л.: Издательство АН СССР, 1935. — 215 с.
  62. A.A. К расчету соединений с гарантированным натягом по прерывистой цилиндрической поверхности / Сб. научных трудов ОрелГТУ. -Орел, 1998.-Т. 13.-С. 68−72.
  63. К.Я. Учет овальностей в расчетах соединений с натягом. -«Изв. АН Лат. ССР», 1954, № 11/88.
  64. С.Д. и др. Расчеты на прочность в машиностроении. М.: Машгиз, 1958. — Т.2.
  65. Расчеты машиностроительных конструкций методом конечных элементов. Справочник / Под редакцией Мячникова В. И. М.: Машиностроение, 1989.- 519 с.
  66. Д.Н. Детли машин. -М., Машгиз, 1975, 723с.
  67. А.Г. Конические прессовые посадки гребных винтов и муфт. М.: Судпромгиз, 1960.-91 с.
  68. Г. Н. Концентрация напряжений около отверстий. М.: Гос-техиздат, 1951. — 496 с.
  69. Г. Н. Распределение напряжений около отверстий. Киев: Наукова думка, 1968. — 888 с.
  70. Л. Применение метода конечных элементов. М.: Мир, 1979.-392 с.
  71. A.A. О натяге посадки деталей при шпоночном соединении / Детали машин и подъемно-транспортные машины // Респ. межвед. научн.-техн. сб. 1968. — Вып. 8.
  72. В.Н. Исследование прочности комбинированных соединений с натягом. Диссертация на соискание ученой степени к.т.н. Львов, 1979. -176 с.
  73. И.П., Ушаков Б. Н. Исследование деформаций и напряжений методом муаровых полос. -М.: Машиностроение, 1969. -208 с.
  74. Н.Д. Определение напряжений в круглой пластинке с несколькими запрессованными в нее круглыми шайбами. ДАН СССР. — 1948. -Т.63. -№ 1, — С. 15−18.
  75. Н.Д. Расчеты на прочность прессовых соединений. ДАН СССР. — 1949. — Т.67. — № 4. — С.615−980.
  76. Н.Д. Расчеты на прочность составных кольцевых соединений. ДАН СССР. — 1950. — Т.70. — № 6. — С. 977−980.
  77. Н.Д. Определение в некоторых деталях напряжений, возникающих от напряженной посадки / В кн.: Расчеты на прочность. 1958. -Вып.2. — С. 142−181.
  78. Н.Д. Новый прием решения задач прессовых соединений / В кн.: Расчеты на прочность. -М., 1969. Вып. 14. — С. 62−75.
  79. С.П. Сопротивление материалов. М., 1965. — Т.2.
  80. А.Б. Несущая способность многослойных прессовых соединений при действии циклических нагрузок. Автореферат диссертации на соискание ученой степени к.т.н. Ижевск, 1996. — 20 с.
  81. А.Г. Построение конформно отображающих функций при помощи электромоделирования и интерполяционных полиномов Лагранжа. -Киев: Наукова думка, 1966. 208 с.
  82. А.Г. Определение напряжений в пластинке, ограниченной улиткой Паскаля, с несколькими запрессованными в нее круглыми шайбами / Инж. сборник. 1953. — Т. 17. — С. 203−206.
  83. А.Г. Определение напряжений в пластинке при запрессовке в нее нескольких круглых шайб, нагруженных сосредоточенными силами и моментами / Труды ГИСИ им. В. П. Чкалова. 1953. — Вып. 24. — С. 146.
  84. А.Г. О расчетах посадочных напряжений в некоторых типах прессовых соединений / Труды ГИСИ им. В. П. Чкалова. 1956. — Вып. 25. -С.31−43.
  85. А.Г. Концентрация посадочных напряжений около отверстий / Сб.: Проблемы прочности в машиностроении. 1962. — Вып.9. — С. 5−14.
  86. Ф.Ф. Разработка методов повышения нагрузочной способности и расчета конических соединений. Автореферат диссертации на соискание ученой степени к.т.н. / Ротопринт. Устинов: ИМИ, 1985.
  87. В.Т. Теоретические основы расчета и обеспечения работо-у/ способности крупных деталей, соединенных натягом. Автореферат диссертации на соискание ученой степени д.т.н. МГТУ, 1990.
  88. В.Т., Гречушкин Г. М. Напряженное состояние и работоспо^ собность крупногабаритных прессовых соединений. // Проблемы прочности. -1990.-№ 3.-С. 77−82.
  89. Дж., Мол ер К. Численное решение систем линейных алгебраических уравнений. М.: Мир, 1969. — 280 с.
  90. Н.В. Расчет эксцентричного цилиндро-конического соединения на контактную прочность и жесткость // Изв. Куйбышевского с.-х. института. 1967. — Т. 22. — С. 245−256.
  91. Н.В. Силовой расчет эксцентричного соединения / Сб. научных трудов Куйбышевского филиала ВЗИИТ. 1967. — С. 120−129.
  92. Д.И. К определению напряжений в анизотропной, упругой и однородной среде, состоящих из нескольких сопряженных между собой посредством посадки тел. Прикладная математика и механика. 1945. — Т. 9. — С. 41−49.
  93. Д.И. О напряженном состоянии некоторых запрессованных^ деталей // Изв. АН СССР. ОТН. — 1948. — № 9. — С.1371−1388. Vх
  94. Д.И. Метод интегральных уравнений в плоских и пространственных задачах статической теории упругости / В кн.: Труды I Всесоюзного съезда по теоретической и прикладной механике. (Обзорный доклад). М.: Наука, 1962.-С. 405−467.
  95. И.В. Эксцентриковые соединения деталей машин типа вал-ступица. Диссертация на соискание ученой степени к.т.н. Харьков: ХПИ, 1984.- 18 с.
  96. А. с. № 1 805 379 СССР. Способ неразрушающего контроля соединений / Шуваев В. Г. // Бюллетень изобрет. 1993. — № 12,.
  97. Хейгеман JL, Ягн Ю. Прикладные итерационные методы. М.: Мир, 1986.-448 с.
  98. А.С. К определению напряжений изгиба в зубьях цилиндрических передач методом граничных конечных элементов. Машиноведение. — 1982.-№ 2.
  99. А.С., Моисеенко А. А. Регулировка локализации контакта зацепления Новикова смещением исходного контура // Известия вузов. Машиностроение. 1989. -№ 10.
  100. Aufgebaute Welle, Amboru Peter и др. / Uni Cardan A. G/ - N 3 844 699. 5: Заявлено 14. 09. 88. Опубликовано 28.06.90. ФРГ.
  101. Besuner P.M., Snow D.W. Application in two-dimentional integral equation method to engineering problems, in: Boundary Integral Equation Method: Computational Application in Applied Mechanics, ed. T.A. Cruse, F.J. Rizzo, ASME, 1975.
  102. Cruse T.A., Application of the boundary integral equation method to three dimensional stress analysis. — J. Computers and Structures. — 3. — 509—527 (1973).
  103. Cruse T.A., Swedlow J.L., Interactive program for analysis and design problems in advanced composites technology, Techn. Rep. AFML TR— 71 — 268, 1971.
  104. Dyson I., Hirst W. Proc. Poy. Soc.. V.67. -№ 412. — Sec. B. — 1954.
  105. Hisada Т., Nikagiri S. A note on stochastic finite element method (Part 2. Variation of stress and strain caused by shape fluctuation). SEISAN — KENKYU, 1980. — v. 32. — № 5. — P. 262−265.
  106. Improving the load carrying capacity of interference fits (Ramamoorthy В., Redhakrishnan V. // Proc. Inst. Mach. Eg. B. — 1989. — 203, № 2. — C. 83−90. -Англ.
  107. Jaswon M.A., Integral eguation methods in potential theory. I, Proc. Roy. Soc., A 275. C. 23−32 (1963).
  108. Mendelson A. Elastic plastic torsion problems for strain-hardening materials, NASATN D 4391, 1968.
  109. Niepage Peter, Grahn Bernd, Wirnich Karl Heinz // Draht, 1993−94. № 4. — C. 224−227. Влияние внутреннего напряженного состояния на вибропрочность тарельчатых пружин.129
  110. Nikagiri S., Hisada T. A note on stochastic finite element method (Part 1. Variation of stress and strain caused by shape fluctuation). SEISAN — KEN-KYU, 1980. — V. 32. — № 2. — P. 39−44.
  111. Riccardella P.C., An improved implementation of the boundary- integral technique for two dimensional elasticity problems, Carnegie-Mellon University Rep. SM-72−26, 1972.
  112. Rizzo F. J., An integral equation approach to boundary value problems of classical elastostatics, Quart. Appl. Math., 25, 83−95 (1967).
  113. Rizzo F. J., Shippy D. J., A method for stress determination in plane anisotropic elastic bodies, J. Composite Mater., 4, 36 61 (1970).
  114. Rizzo F.J. The boundary integral equation method: a modern computational procedure Cn applied mechanics, Appl. Mechanics Conference, New York, 1975.-P. 11−17.
  115. Si Karskic D.L., Cheatham J.B. Penetration problem in roch mechanics, in «RockMechanics Symposium», Vol.3, Appl. Mech. Div., ASME, 41−71, 1973.
  116. Snyder M.D., Cruse T.A. Boundary- integral equation analysis of cracked anisotropic plates, Intern, J. Fracture, 11, 315−328 (1975).
  117. Speek Lefotografische Deformations — und Spannungsanorlyse tor-sionsbeanspruchter Langspre Bverbindungen / Vogel jurgen U Techn Mech. — 1984. — 10. -№ l.-C.l, 17−25.
  118. Yannakakis M. Computing the minimum fill in is NP complete // SIAM J. Alg. Disc. Meth. — 1981. — № l.-P. 77−79.
  119. Г. Б. Детали машин. -M.Машиностроение, 1968.-368с.
  120. Мак-Кракен Д., Дорн У. Численные методы программирования на Фортране. М.: Мир, 1977. С. 331.
  121. ХеммингР. В. Численные методы. М.: Наука, 1968. С. 233.1.> (•1. Т'|ПРО"im,?.niiu ¦"< 10'i В1. Щ ц 4 AHK1. АКТвнедрения программ по расчету прочности соединений с гарантированным натягом по прерывистой цилиндрической поверхности
  122. SET SAFETY OFF && HE ВЫДАВАТЬ ПРЕДУПРЕЖДЕНИЯ ПО ZAP, INDEX ON,.1. SET PROC TO PRESS1. CLOS ALL1. CLEA ALL1. RELE ALLset date to dmy
  123. USE BAZA ORDE TAG KL @ 0,0 CLEAR1,30 SAY 'Программа PRESS' ETAP= 1* индексы для отладочной печати @ 0,0 CLEAR1. ЕТАР=5rest from press addl1. ETAP=51. ENDI1. ETAP<2* индексы для отладочной печати J0=0
  124. STOR 0 TO ALA, ALF, GAM, GF, HA, HF, LA, LF, R1, R2,RF, RG, Q1, Q2,WESS1,XA, XF, US1 US2
  125. DO PECHW WITH ' MU=' +TRAN (MU ,'@A')1. ENDIF
  126. AZ=INPUT ('AZ', IX0, IX, IZ, 0,0) A0=INPUT ('A0', IKN, IX0, I, K,0) AIZ=AIZ+A0* AZ1. ENDF
  127. OUTPUT ('AIZ', AIZ, IKN, IX, I, IZ, K) ENDF ENDF^Q ***
  128. WAIT WIND ' запись A ()=0, ждите' NOWA1. FOR J=1 TO 2* 1ST &&D0401. FOR 1=1 TO 2*IST
  129. OUTPUT (, A,0,I, J, D, 0,0) && запись A (I, J) = 0 =OUTPUT ('AU', 0,1,J, 0,0,0) && запись U (I, J) = 01. ENDF1. ENDF
  130. WAIT WIND ' запись S 10=0, ждите' NOWA1. FOR J=1 TO 1ST &&D0 411. FOR 1=1 TO 1ST
  131. OUTPUT ('S1', 0,1,U, 0,0) && запись S 1(1,1,J) = 0 =OUTPUT ('S 1 ', 0,2,1,J, 0,0) && запись S1(1,I, J) = 01. ENDF1. ENDF
  132. DELE FILE 'PRESS.MEM' save to press ENDIF&& ETAP<21. ETAP<3. OR. ETAP=51. ETAP<31. REST FROM PRESS1. ETAP=21. J=2*IST
  133. PUBL ARRA B (J), SJ (2,2), A00(2,2), XIZ0(2) && текущие XIZ, AO PUBL ARRA UJ (2,2), S (2,2), U (2,2), XIZKJ (2), AIZKJ (2,2)
  134. PUBL ARRA SIG1(J), SIG1(J), SIGZ (J)
  135. STOR 0 TO B, SJ, UJ, S, U, XIZKJ, AIZKJ, SIG1, SIG2,SIGZ
  136. PUBL A1, A2,A3,A4,A5 *** ***
  137. A1 =(3. -2. * MU)/(1. -MU)/(4 * PI) A2=(I.-2.*MU)/(I.4VIU)/(4*PI) A3=(1 .+2.*MU)/(1 .-MU)/(4*PI) A4=(3.-4.*MU)*(1.+MU)/(4*PI*E*(1 ,-MU)) A5=(1.+MU)/(4*PI*E*(1 ,-MU)) A6 = A2 + A3 && 9.03.98
  138. CDL = CDL * (2. CDS) && 9.03.981. J00=11. ELSE1. ETAP=51. ETAP=2 J00=J'Расчет продолжается с ETAP=2 J-, j wait 'нажмите Enter'1. ENDI0,0 CLEA
  139. WAIT WIND ' второй цикл J от '+str (j00,3)+' до '+str (ist, 3) NOWA REPL VALUE WITH 0 FOR
  140. NAME-S1 '. AND. I3≥J00). OR.NAME=A'. OR.NAME-AU'
  141. FOR J=J00 TO 1ST && DO 12 J= 1, ISTначало формирования A (I, J), B (J)
  142. SJ (1,1)=(1.47.40*MU-2*CDS*(1−2*MU)) * NU (J)/(2*PI*(1-MU))) /(2*DL (J)) SJ (2,2)=(1 .-(5−2*CDS) * (1−2*MU) * NU (J) / (2*PI*(1-MU))) / (2*DL (J))
  143. UJ (l, l)=A4*(l.-LOG (DL (J)/2.)) UJ (2,2)=UJ (1,1)-A5 IF J≤IOT 111=1 II2=IOT1. ELSE *** ***1.l=IOT+l II2=IST1. ENDI *** ***
  144. XIZ0(1)= INPUT ('XIZ', 1,1, IZ, K,0) && это есть XIZ (IX0,I, IZ, K) XIZ0(2)= INPUT ('XIZ', 2, I, IZ, K,0) FOR IKNJ=1 TO 2 && DO 15
  145. XIZKJ (IKNJ)=0. FOR 1X0=1 TO 2 && DO 15
  146. XIZKJ (IKNJ)=XIZKJ (IKNJ) + XIZO (IXO) * AOO (IKNJJXO)1. ENDF1. ENDF && IKNJ=1 TO 2j^ ***1. FOR IKN=1 TO 2 && DO 16
  147. S (1,1)=S (1,D+SJ (1,1) S (2,2)=S (2,2)+SJ (2,2)1. ELSE
  148. DX1 =XIKNJ (J, 1) XIZKJ (l) && INPUT ('XIZ', 1 .J, IZ, K,0)1. XIZKJ (l)1. DX11=0X1*0X1
  149. DX2=XIKNJ (J, 2) XIZKJ (2) && INPUT ('XIZ', 2, J, IZ, K,0)1. XIZKJ (2)
  150. DX22=DX2 * DX2 DX12=DX1*DX2 R=SQRT (DX1 RDX22) *** 27/02/1998 (ИЗМЕНЕНИЕ ПРОГРАММЫ) IF CDS=2.1. DLJ=DL (J)
  151. XX2=DX2 -D L J * NU (J)/4 XX22=XX2*XX2. XI l=DXl-DLJ/2 X12=DXl+DLJ/2 XI 12=X11*X11 X122=X12*X12 R1=SQRT (XX22+X112)1. RR2=SQRT (XX22+X 122)1. R12=R1*R11. R22=RR2*RR2
  152. SI1=(ATAN (X 12/xx2)-ATAN (Xl l/xx2))/DLJr02 = xx22 + xl 1 * xl2
  153. SIl=ATAN (xx2 * dlj / r02) / DLJ S I2=LOG (R22/R 12)1(2. * DLJ)
  154. S I3=XX2* (X12/R22-X11 /R12)/(2. * DL J)+S 11 /2
  155. S I4=XX22 * (1 ,/R 12−1 ,/R22)/(2. * DLJ)1. SI5=SI1-SI31. SI6=SI2-SI41. UI1=XX2*SI11. UI2=XX2*SI21. UI3=1.-UI1
  156. UI4=(X 12 * L0G (RR2)-X 11*L0G (R1))/DL J-UI31. SOI 11=A2*SI2+A6*SI61. S0221=A3*SI2-A6*SI61. S0121=A2*SI1+A6*SI51. S0222=A2*SI1+A6*SI31. SOI 12=A3*SI1-A6*SI31. S 0212=A2 * S12+A6 * S14
  157. SO тензор средних напряжений на DL (J) от единичной силынаправления 3-го индекса1. U011 =-A4*U 14-А5 *UI 11. U012=A5*UI21. U021=U0121. U022=-A4*UI4-A5*UI3
  158. UO вектор среднего перемещения на DL (J) от единичной силы 2-го индекса
  159. U11 =U011 * AIZKJ (1,1)+U012* AIZKJ (1,2) U21=U021*AIZKJ (1,1)+U022*AIZKJ (1,2) U12=U011 * AIZKJ (2,1)HJ012* AIZKJ (2,2) U22=U021*AIZKJ (2,1)+U022*AIZKJ (2,2)
  160. UIJ среднее перемещение направ.1 HaDL (J) от ед. нап. силы направления J на элементе I1. ELSE
  161. R2=R*SQRT (DX11+DX22-CDL*DL (J)**2) R4=R2*R*R
  162. W1=DX2*(A1 *DX11+A2*DX22) W2=DX1*(A3*DX22-A2*DX11) W3=A1 *DX22+A2*DX11
  163. SSS1=L/R4*(W1*AIZKJ (L1)+DX1*(A2*DX11+A1*DX22)*AIZKJ (1,2))
  164. SSS2=1 ./R4*(W2* AIZKJ (1,1)+DX2* W3 * AIZKJ (1,2)) SSS3=1 ./R4*(W1 *AIZKJ (2,1)+DX1 * W3*AIZKJ (2,2))
  165. SSS4=1./R4* (W2* AIZKJ (2,1)+DX2* W3 * AIZKJ (2,2)) ENDI
  166. S (U)=S (1,1)-SSS1 S (2,1)=S (2,1)-SSS2 S (1,2)=S (1,2)-SSS3 S (2,2)=S (2,2)-SSS41. ENDI30 ***1. MU)*DL (J))
  167. I=J.AND.IZ=1, AND. K=1 SSS1=0
  168. SSS2=MU*(l.+(5−6.*MU)*NU (J)/(2.!, iPI*MU))/(2.*(l .1. ELSE27/02/98 Начало изменений IF CDS=2. SSS1=-S111 SSS2=-S112 ELSE27/02/98 окончание изменений
  169. W1=DX1*(A1*DX11+A2*DX22) W2=DX2 * (A3 * DX11-A2*DX22) S S S1 =-1 ,/R4*(W 1 * AIZKJ (1,1)+W2 * AIZKJ (1,2)) SSS2=-1 ./R4*(W1 * AIZK J (2,1)+W2 * AIZK J (2,2)) ENDI1. ENDI1. ПЕЧАТЬ
  170. S11 =INPUT ('S 1', 1,1,J, 0, 0)+SSSl
  171. S 12=INPUT ('S 1 ', 2, I, J,0, 0)+SSS2 =OUTPUT ('S 1 ', S 11,1,1,J, 0,0) =OUTPUT ('S 1 ', S 12,2,I, J,0,0)28***
  172. SSS1=UJ (1,1) SSS2=UJ (2,2) U (1,1)=U (1,1)+SSS1 U (2,2)=U (2,2)+SSS21. ELSE27/02/98 Начало изменений IF CDS=2.
  173. SSS1=-U11 SSS2=-U21 SSS3=-U12 SSS4=-U221. ELSE27/02/98 Окончание изменений1. Wl=A4*R2*LOG®+A5*DX22
  174. SSS1=1VR2*(W1 *AIZKJ (1,1)-A5*DX12*AIZKJ (1,2))1. SSS2=1./R2*(
  175. A5 *DX 12* AIZKJ (1,1)+(A4*R2*LOG®+A5*DX 11)* AIZKJ (1,2))
  176. SSS3=1 ,/R2*(W 1 * AIZKJ (2,1)-A5 * DX12* AIZKJ (2,2)) SSS4=1./R2*(
  177. A5*DX12*AIZKJ (2,l)+(A4*R2*LOG®^5*DXl 1)*AIZKJ (2,2))1. ENDI
  178. U (1,1)=U (1,1)-SSS1 U (2,1)=U (2,1)-SSS2 U (1,2)=U (1,2)-SSS3 U (2,2)=U (2,2)-SSS41. ENDI
  179. W1 =A4*R20*LC)G (R0)+A5 *DX220
  180. SSS1=L/R20*(W1*AIZKJ (1,1)-A5*DX120*AIZKJ (1,2))1. SSS2=1./R20*(
  181. A5*DX120*AIZKJ (1,1)+(A4 * R20* LOG (RO)+A5 *DX110) * AIZKJ (1,2))
  182. SSS3=1 ./R20*(W1 *AIZKJ (2,1)-A5*DX120*AIZKJ (2,2))
  183. SSS4=1./R20*(-A5*DX120*AIZKJ (2,1)+(A4*R20*LOG (RO)+A5*DX110)*AIZKJ (2,2))
  184. U (1,1)=U (1,1)+SSS1 U (2,1)=U (2,1)+S S S 2 U (1,2)=U (1,2)+SSS3 U (2,2)=U (2,2)+S S S41. ENDI
  185. ENDF && K=1 TO 2 && DO 20 ENDF && IZ=1 TO IZS && DO 20* * 20 * * *1. FOR IKN=1 TO 2 && DO 211.=2*(I-1)+IKN
  186. FOR IKNJ=1 TO 2 && DO 21 JJ=2*(J-1)+IKNJ1. A (II, JJ)=S (IKNJ, IKN)
  187. OUTPUT ('A', S (IKNJ, IKN), II, JJ, 0,0,0)1. AU (II, JJ)=U (IKNJ, IKN)
  188. OUTPUT ('AU', U (IKNJ, IKN), II, JJ, 0,0,0)1. ENDF1. ENDF
  189. ENDF && 1=111 TO 112 && DO 22* 22 * * *
  190. ETAP<0 && прерван расчет ETAP=31. DELE FILE 'PRESS.MEM'save to press1. QUIT1. ENDI
  191. ENDF && J=j00 TO 1ST && DO 12 J=1,IST12 ***1. FOR J=1 TO 2*IST && DO 25
  192. SJJ=l+INT (J/2)*4−2*J IF J<=(2*IOT).OR.J>(2*(IOT+Il)) B (J)=0.1. ELSE1. B (J)=NAT* (S JJ+1)/4.1. ENDI
  193. OUTPUT ('B', B (J), J,0,0,0,0) && запись B (J) свободный член FOR 1=1 TO 2*IST && DO 25'вычисление i=', i, 'j=', j IF J>(2*I1) *** 24 *** запись матрицы1. ! (J<=(2^IOT).OR.J>(2^(IOT+Il)))
  194. A (I, J)=-AU (I, J)-AU (I, J-2*IOT)
  195. OUTPUT ('A',-INPUT ('AU', I, J, 0,0,0)-INPUT ('AU', I, J2*IOT, 0,0,0), I, J, 0,0,0)1. ENDI1. ELSE
  196. A (I, J)=A (I, J)-A (I, J+2*IOT)
  197. OUTPUT ('A', INPUT ('A', I, J, 0,0,0)-INPUT ('A', I, J+2*IOT, 0,0,0), 1., 0,0,0)1. ENDI1. ENDF1. ENDF E=E0j" ^ ^^ H4% %
  198. DELE FILE 'PRESS.MEM' save to press1. ENDIF && ETAP<3 RELE ALL1. DO PRESS2 && В НЕЙ ВЫХОД1. PROC GOT1. PRIV A00, X00,C1. C=ROT/RST1. FOR I=IS+1 TO IOT
  199. DL (I)=C*DL (I+ISTF-IS) NU (I)=-NU (I+ISTF-IS) FOR IT=1 TO 21. FOR K=1 TO 2
  200. X00=INPUT ('XO', IT, I+ISTF-IS, K,0,0) X00=X00*C
  201. OUTPUT ('X0', X00, IT, I, K, 0,0) FOR IU=1 TO 2
  202. A00=-INPUT ('A0', IU, IT, I+ISTF-IS, K,0) =OUTPUT ('A0', A00, IU, IT, I, K, 0)1. ENDF1. ENDF1. ENDF1. ENDF RETU1. PROC STUPif kstup = 1 && старый вариант1. Jl=IOT+l1. J2=IOT+Il1. FOR J=J1 TO J2 && DO 1
  203. J-IOT nu (j)=-nu (i) DL (J)=DL (I)
  204. XK (1,J)=XK (1,I) XK (2,J)=XK (2,I)1. ENDF* J ** *s=pi/izs-fk1
  205. S -угловая протяженность эл-тов отверстия ступицы без первых ее II эл-тов SA=S*RA/DL (I1)1. QQ=1.51.=ISTF-I0T-I11. FOR 1=1 ТО 10 && DO 2
  206. Q=(l .+SA*(QQ-1 .)/QQ)**(l ./ID) QQ=Q1. ENDF1. C ifc 5}C ^ ^ ^
  207. DL1=DL (I1) DFK=DL (I1)/RA FK=FK1 J=IOT+Il+l1. FOR I=J TO ISTF && DO 3
  208. DFK=DFK*Q DL1=DL1*Q DL (I)=DL1 FK=FK+DFK NU (I)=-DFK/2. XK (1,I)=RA* SIN (FK) XK (2,I)=RA*COS (FK)
  209. DFK и FK угловые протяженности1-го эл-та и его конца (координаты) ENDF1. Н4 Н4 Н4 ^ 4е 4s И41. XN1=0. XN2=RA
  210. FOR I=IOT+l TO ISTF && DO 4
  211. DX1 =XK (1,1)-XN 1 DX2=XK (2,I)-XN2 SDX=SQRT (DX 1 * * 2+DX2 * * 2) A01I=-DX1/SDX1. A02 I=DX2/SDX
  212. OUTPUT ('AO', AO11,1,1,1,1,0)
  213. OUTPUT (AO', AO2I, 2,1,1,1,0)
  214. OUTPUT (A0',-A02I, l,2,I, l,0)
  215. OUTPUT (A0', A0lI, 2,2,I, l,0)
  216. OUTPUT (A0',-A0lI, l, l, I,2,0)
  217. OUTPUT (A0',-A02I, 2, l, I,2,0)output (ao',-ao2i, 1,2,1,2,0)
  218. OUTPUT (AO', AO11,2,2,1,2,0)1. X01I=-RA*A02I1. XO2I=-RA*AO1I
  219. OUTPUT ('XO', XO11,1,1,1,0,0)
  220. OUTPUT ('XO', XO2I, 2,1,1,0,0)
  221. OUTPUTCXO',-XO11,1,1,2,0,0)
  222. OUTPUT ('X0', X02I, 2, I, 2,0,0)1. XN1=XK (1,I)1. XN2=XK (2,I)endf1. Ф ^ Ф ^ H®dfk=pi/izs/(ist-istf)else && kstup = 1 *** новый вариант STUP
  223. DFK=PI/IZS/(IST-ISTF) ISTF = IOT+IS
  224. FOR i=IOT+l TO ISTF && DO 1 DL (I)=DL (I-IOT) NU (i)=-NU (I-IOT) FOR 1X0=1 TO 21. FOR K=1 TO 2xoo=input ('xo', ixo, i-iot, k, o, o)
  225. OUTPUT ('X0', X00, IX0,I, K,0,0) FORIKN=1 TO 2aoo=-input ('ao', ikn, ixo, i-iot, k, o)
  226. OUTPUT (A0', A00, IKN, IX0, I, K,0)1. ENDF1. ENDF1. ENDF
  227. ENDF endi && kstup = 1 * * * общая часть
  228. FOR I=ISTF+1 TO 1ST && DO 5
  229. FI=DFK*(I-ISTF-0.5) AO1I=COS (FI) A02I=SIN (FI)
  230. OUTPUT ('AO', A011,1,1,1,1,0) =OUTPUT (AO', AO2I, 2,1,1,1,0) =OUTPUT ('AO',-AO2I, 1,2,1,1,0) =OUTPUT (A0', A0lI, 2,2,I, l,0) =OUTPUT (A0',-A0JI, l, l, I,2,0) =OUTPUT (A0',-A02I, l,2,I, 2,0) =OUTPUT (A0',-A02I, 2,1,1,2,0)
  231. XK (2,i)---------±------DL (i)------до 15 XK (2,i)---------±------DL (i)
  232. XK (2,i)---------±------DL (i)1. ПОДПРОГРАММА «GEW"0,0 clea WESS3=2.20,0 SAY 'Задайте WES S3 =' GET WES S3 PICT '99.99 999 999' READ
  233. DO PECHW WITH ' WES S3- +TRAN (WES S3 ,'@A')0,0 clea1. ALF = GR*ALF1. HAL = HA*(1 ,-CDL)1. AA1 = PI/4. HAL*TAN (ALF)1. CEL1 = WESS1*AA1
  234. CEL2 = HA*(2.-CDL) / COS (ALF)
  235. CEL3 = WESS3*RF* (PI/2. ALF)
  236. Начало вставки 13.01.99. IF CDS=21. RA=IZS/2.+(l ,-CDL)*HA+X1. ENDI1. Конец
  237. ALA = ACOS (COS (ALF) * IZS/(2.*RA))
  238. FK1 = (PI/2. + 2.*X*TAN (ALF)) / IZS + TAN (ALF) ALF — TAN (ALA) + ALA F2 =0.1. FOR 1=1 TO II1. I>(Il-5)1. QWl=Ql**(I-Il+5)
  239. F1 = FKl*(Il-5+2.*(l.-QWl) / (l.-Ql**5)) / (11−3)1. ELSE1. F1 = FK1 * 1/(11−3)1. ENDI
  240. FK1 угловая протяженность участка выступов
  241. NU (I) = (F1-F2) / 2. ХК (1 Д) = RA* SIN (F 1) XK (2,I) = RA*COS (Fl)31. F2 = F11. ENDF2. й участок
  242. RI = SQRT (IZS+>1:2/4. IZS*(HA-X) + ((HA-X) / SIN (ALF))**2)1. RI мин радиус
  243. ALI = ACOS (COS (ALA)*RA/RI)
  244. W = IZS*C0S (ALF)*(1./C0S (ALA)**2 l./COS (ALI)**2) / 4.
  245. W лин протяженность эвольв участка
  246. F2 = TAN (ALA) TAN (ALF) — (PI/2. +2.*X*TAN (ALF)) / IZS
  247. F2 угол на доп диаметре точки нормали начала 2-го участка
  248. FOR 1=11 + 1 ТО 12 IF I> (11+6)
  249. I = LEW*(I-I 1−4) / (12−11−4)1. ELSE1. QW1=Q2**(i-i1)
  250. I = LEW*2.*(QW1 -1.) / (Q2**6 1.) / (12−11−4)1. ENDI1. 1 линейная протяженность (общая) I≤6 эл-тов
  251. ALI = ACOS (SQRT (IZS*C0S (ALF)*C0S (ALA)**2 /- (IZS*COS (ALF) 4.*LE1*C0S (ALA)**2)))
  252. F1 центральный угол (влево от Х2),
  253. NU (I)=(F2-Fl)/2. FK=F 1+ALF
  254. XK (l, I)=-SIN (Fl)*IZS/2.+Ll*COS (FK) XK (2,I)=COS (Fl)*IZS/2.+Ll*SIN (FK)
  255. FK угол между нормалью и осью XI1. F2 = Fl1. ENDF3. й участок
  256. XRF=PI/4 ,+HF * TAN (ALF)+RF * TAN (PI/4. ALF/2.)1. A2=ALF1. HRF=HF-RF-X
  257. F2=2.*XRF/IZS+2.*HRF/(IZS*TAN (ALF))
  258. F2 угол нормали начала 3-го участка с осью XI FOR 1=12+1 ТО 13
  259. Al=ALF+(PI/2.-ALF)*(I-I2)/(I3-I2)
  260. F1 =2. *XRF/IZS+2. *HRF/(IZS * TAN (A 1))1. NU (I)=(A2-F2-A 1+F1)/2.1.=RF+HRF/SIN (A1)
  261. XK (l, I)=IZS*SIN (Fl)/2.-Ll*C0S (Al-Fl)
  262. XK (2,I)=IZS*COS (Fl)/2.-Ll*SIN (Al-Fl)1. A2=A11. F2=F11. ENDF && ************ 91. RFF=IZS/2.-HF+X1. F2=2.*XRF/IZS1. F20=F21. ALA=GAM ENDI END I1. Конец1. RF=R21. ALFF=60.*GR1. ENDI
  263. CEL1=(R1 *COS (ALA)-LA)*WESSl1. CEL2=R1 *(ALA-GAM)1. CEL3=RG*(GF-GAM)1. CEL4=R2*(ALFF-GF)
  264. CEL5=RF*(PI/2.-ALFF)*WESSl
  265. CELS=CEL1+CEL2+CEL3+CEL4+CEL51. PRIN > 2
  266. Начало вставки 13.01.99. RA=IZS/2.+H A+X IF CDS=21. RA=IZS/2.-fHA*(l .-CDL)+X1. ENDI
  267. Конец IF RF=R2 AA=PI/2 L=R1-XA+X1. ELSE1=1
  268. A A1 = AS IN ((H A+X A)/R 1) if- J * >fi *1. DO WHIL .T.
  269. R1-(XA-X)/SIN (AA1) AA=ASIN ((RA**2-(0.5*IZS)**2-L**2)/(IZS*L)) IF ABS (AA-AA1)<=0.00001 EXIT1. ENDI IF I>50
  270. WAIT WIND 'Нет сходимости' QUIT
  271. А1 =AA-(AA-GAM)*(I-I1 -4)/(I2−11 -4)1. ELSE1. QW1=Q2**(I-I1)
  272. A1=AA-2.*(AA-GAM)*(QW1−1.)/(Q2**6−1.)/(I2-I1 -4)
  273. Q2 коэф возрастающей геом прогрессии1. END I1. Вставка 22.09.99.1. END I1. Конецi ^ g ^? ^
  274. Fl=2.* (L A-(X A-X)/T AN (A1))/IZS1.R1-(XA-X)/SIN (A1)1. NU (I)=(F2-Fl+A2-Al)/2.
  275. XK (1,I)=-0.5*IZS*SIN (F1)+L*COS (F1+A1)
  276. XK (2,I)=0.5*IZS*COS (F 1)+L* SIN (F 1+A1)1. A2=A11. F2=F11. ENDF^ H5 ^ ^ ^1.=(Rl+RG)*COS (GAM)-LA XG=(R 1+RG) * SIN (G AM)-X A A2=GAM
  277. F2=2.*(LG-(XG+X)/TAN (GAM))/IZS FOR 1=12+1 TO 13 && DO 9
  278. A1 =GAM+(GF-GAM)*(I-I2)/(I3 -12) F1 =2. *(LG-(XG+X)/TAN (A1))/IZS L=RG-(XG+X)/SIN (A1) NU (I)=(A2-F2-Al+Fl)/2. XK (1,I)=0.5*IZS*SIN (F1)-L*COS (A1-F1) XK (2,I)=0.5*IZS*COS (F 1)-L* SIN (A 1 -F 1) A2=A1 F2=F11. ENDF^ ^ Q % ^1. A2=GF
  279. F2=2.*(PI/2.+LF-(XF+X)/TAN (GF))/IZS FOR 1=13+1 TO 14 &&D0 10
  280. A1 =GF+(ALFF-GF)* (I-I3)/(I4-I3) Fl=2.*(PI/2.+LF-(XF+X)/TAN (Al))/IZS L=R2-(XF+X)/SIN (A1) NU (I)=(A2-F2-A1 +F1)/2. XK (1,I)=0.5*IZS*SIN (F1)-L*COS (A1-F1) XK (2,I)=0.5*IZS*COS (F1)-L*SIN (A1-F1) A2=A1 F2=F11. ENDF Q * * =f=1. HFX=HF-RF-X A2=ALFF
  281. F2=2.*(PI/2.+HFX/TAN (ALFF))/IZS1. FOR 1=14+1 TO 15 && DO 11
  282. A1 =ALFF+(PI/2. ALFF) * (I-I4)/(I5 -14)
  283. F1 =2. * (PI/2.+HFX/T AN (A 1))/IZS1.RF+HFX/SIN (A1)1. NU (I)=(A2-F2-A 1+F1)/2.
  284. XK (1,I)=0.5*IZS*SIN (F1)-L*COS (A1-F1)
  285. XK (2,I)=0.5*IZS*COS (F1)-L*SIN (A1-F1)1. A2=A11. F2=F11. ENDF1. J i ***1. XI N=0. X2N=RA1. FOR 1=1 TO 15 &&D0 12
  286. DO PECHW WITH ' ПЕЧАТЬ ДЛЯ RST =101'
  287. DO PECHJW WITH 'CD ='+TRAN (CD,'@A')+' DO = '+TRAN (D0,'@A')+' DD = '+TRAN (DD,'@A')+' CM = '+TRAN (CM,'@A')+' X = '+TRAN (X,'@A')
  288. DO PECHJW WITH 'HA ='+TRAN (HA,'@A')+' HF = '+TRAN (HF,'@A')+' R1 = '+TRAN (R1,'@A')+' R2 = '+TRAN (R2,'@A')+' GF = '+TRAN (GF,'@A') FOR I = 12 TO INT (I3 + 0.3 * (14 -13))
  289. DO PECHW WITH 'I ='+TRAN (I,'999')+' XK (1,I) = '+TRAN (XK (1,1),'@A')+' XK (2,I) = '+TRAN (XK (2,I),'@A')+' DL (I) = '+TRAN (DL (I),'@A') ENDF1. ENDI RETU1. ПОДПРОГРАММА «GN50"0,0 С lea ISTF=IOT+IS IF IOT>IS
  290. T=ISTF+IOT-IS ELSE IST=IST END I I5=IS4,0 SAY 'KONTUR=50'6,0 SAY ' IZS='+STR (IZS, 3)+' IS='+STR (IS, 3)+' ISTF='+STR (ISTF, 3)+' IST='+STR (IST, 3)8,0 say1 11=' get ii pict '999' read @ 0,0 clea
  291. DIME F (2*IST), FK (2*IST) STOR 0 TO F, FKci = (Q2**I1 1) / (Q2**H 4- Q2**(IS-I1) — 2) fz = PI ()/IZS fa = fz*c1 fd = fz-fa fn = 0.
  292. DO PECHJW WITH 'CI =' + TRAN (C1,'@A') IF priN > 2do PECHW with ' ***** gn50 *****. do pechw with 'fz, fa, fd='±tran (fz,'@a')+tran (fa-@a')+tran (fd,'@a')endi1. II=IS1. FOR 1=1 TO IS &&D0 56
  293. F (I) = FZ * (1−0.5) / IS DL (I) = FZ / IS NU (I) = DL (I) * 0.5^ ^ 56 ^ ^ ^ ENDF1. ELSE 5:> ^ ^1. FOR 1=1 TO II &&D0 51
  294. FK (I) = FA * (1-Q1**I) / (I-Q1**I1) F (I) = (FK (I)+FN) * 0.5 NU (I) = FK (I) F (I) DL (I) = 2 * NU (I) FN = FK (I)endf* * ^ J * * *1. FOR 1=11+1 TO IS &&D0 52
  295. FK (I) =FA+ FD * (Q2**(I-I1)-1) / (Q2**(IS-I1)-1)
  296. F (I) = (FK (I)+FN) * 0.5 NU (I) = FK (I) F (I) DL (I) = 2 * NU (I) FN = FK (I)1. ENDF1. ENDI1. RA=1 D=2.1. DO 55
  297. DO PECHW WITH 'I, F (I)"DL, A0(1,1,I, 1), A0(2,1,1,1) ='±
  298. TRAN (I,'999')+TRAN (F (I),'9 999 999.9999)+TRAN (DL (I),'9 999 999.9999')+TElAN (A09999.9999')+TRAN (A02I,'9 999 999.9999')
  299. ENDI X01I=A02I XO 2I=A0 II
  300. DO PECHW WITH T, X0(1,I, 1), X0(2,I, 1)
  301. TRAN (I,'999')+TRAN (X0l I,'9 999 999.9999')+TRAN (X02I,'9 999 999.9999') * ENDI
  302. OUTPUT ('XO', XO11,1,1,1,0,0) =OUTPUT ('XO', XO2I, 2, I, 1,0,0) =OUTPUT ('XO',-XO11,1,1,2,0,0) =OUTPUT ('X0', X02I, 2, I, 2,0,0)
  303. NU (I)=-NU (I-IOT) DL (I)=DL (I-IOT) FOR IKN=1 TO 21. FOR K=T TO 2
  304. X0lI=INPUT ('X0', IKN, I-IOT, K,0,0) =OUTPUT ('X0', X0lI, IKN, I, K, 0,0) FOR 1X0=1 TO 2
  305. A0lI=-INPUT ('A0', IKN, IX0, I-IOT, K,0) && ??? знак + =OUTPUT ('A0', A0lI, IKN, IX0, I, K,0)1. ENDF1. MO
  306. FOR I=IOT+l TOISTF FOR I=IOT+l TO IOT+I11. DO 56 && DO 561. ENDF1. ENDF56 ***1. Вставка от 04.12.99.1. FN=FA
  307. FOR I=IOT+I1 +1 TO ISTF && DO 57
  308. FK (I)=FA+FD*(Q2* *(I-IOT-I1)-1 .)/(Q2* * (ISTF-IOT-11)-1.)1. F (I)=(FK (I)+FN)/2.1. FN=FK (I)1. NU (I)=F (I)-FK (T)1. DL (I)=-2.*NU (I)1. AO1I=-COS (F (I))1. A02I=-SIN (F (I))
  309. ОиТРиТ ('АО', АО11,1,1,1,1,0)
  310. OUTPUT (' AO', A02I, 2,1,1,1,0)
  311. OUTPUT ('AO',-AO21,1,2,1,1,0)
  312. OUTPUT ('A0', A0lI, 2,2,I, l,0)
  313. OUTPUT ('AO',-AO11,1,1,1,2,0)
  314. OUTPUT (, AO',-AO2I, 2,1,1,2,0)
  315. OUTPUT ('AO',-AO21.1,2,1,2,0)
  316. OUTPUT ('AO', A0 11,2,2,1,2,0)1. XO1I=-AO2I1. X02I=-A01I
  317. OUTPUT ('XO', XO11,1,1,1,0,0) =OUTPUT ('XO', XO2I, 2, I, 1,0,0) =OUTPUT ('XO',-XO11,1,1,2,0,0) =OUTPUT ('X0,X02I, 2, I, 2,0,0)1. ENDF
  318. FOR I=ISTF+1 TO 1ST && DO 57
  319. FOR I=IS+1 TO IOT J=ISTF+I-IS NU (I)=-NU (J) DL (I)=DL (J)*COST FOR K=1 TO 2 FOR 1X0= 1 TO 2
  320. W = INPUT (X0', IX0, J, iL0,0) W = W* COST
  321. OUTPUT ('X0,W, IX0, I, IC, 050) FOR IKN=1 TO 2
  322. W=INPUT ('A0', IKN, IX0, J, K, O) =OUTPUT ('A0',-W, IKN, IX0, I, K,0)1. ENDF ENDF ENDF ENDF1. PRin > 2
  323. FOR 1=1 TO 1ST DO PECHW WITH 'I, F, FK, DL, NU ='±
  324. TRAN (I,'99')+' '+TRAN (F (I)-'@A')+' 4TRAN (FK (I),'@A')+' '+TRAN (DL (I),'@A')+' '+TRAN (NU (I),'@A')1. ENDF1. ENDI1. RETU *1. ПОДПРОГРАММА «GAUSS»
  325. PARA КАК && .T. распечатать исходную и результат на экране PRIV B, C, J, T, S, L, YES SET TALK OFF SET RESO OFFрешение системы лин ур-ний методом Жордана-Гаусса1. N число уравнений
  326. A (N, N+1) расширенная матрица коэффициентов:
  327. А (1,1)*Х1+А (1,2)*Х2+. = А (1, N+1) ур-ние 1
  328. А (2,1)*Х1+А (2,2)*Х2+. =A (2,N+1) ур-ние 2
  329. A (N, 1)*X1+A (N, 2)*X2+. =A (N, N+1) ур-ние N textпример: DIME А (10Д0)1. N=3
  330. А (1,1)=4 А (1,2)=3 А (1,3)=1 А (1,4)=13
  331. А (2,1)=2 А (2,2)=-1 А (2,3)=-1 А (2,4)=-3
  332. А (3,1)=7 А (3,2)=1 А (3,3)=-3 А (3,4)=0решение 1,2,3 endt1. КАК’Расширенная матрица до решения'
  333. DO Р56 000 IN GAUSS && вывести расширенную матрицу на экран1. ENDI1. FOR S=1 ТО N1. YES=.T. && нет решения1. FOR T=S ТО N1. A (T, S) oO
  334. DO P53000 IN GAUSS && переставить ур-ние T на место S C=1/A (S, S) && умножить на этот коэф, чтобы A (S, S) стал =1
  335. DO Р54 000 IN GAUSS && разделить строку на С FOR L=1 TON IF L<>S1. C=A (L, S)
  336. DO P55000 IN GAUSS && вычесть из строки Lстроку S, умноженную на С1. ENDF1. YES=.F.exit1. ENDI1. ENDF IF YES
  337. WAIT WIND 'Нет единственного решения, ранг <'+STR (T, 4) SUSP1. ENDI1. ENDF IF КАК’Результат решения'
  338. T='Строка '+STR (J, 3)+':' FOR S=1 TO N+1
  339. T=T+TRAN (A (J, S),'999 999.999')1. ENDF? T1. ENDF RETU-----------------конец GAUSS1. ПОДПРОГРАММА «FINISH"расчет после решения системы set reso off
  340. SET TALK OFF SET DELE ON SET EXAC ON
  341. SET SAFETY OFF && HE ВЫДАВАТЬ ПРЕДУПРЕЖДЕНИЯ ПО ZAP, INDEX ON,. SET DECI TO 12 set blin off
  342. SET PROC TO PRESS CLOS ALL CLEA ALL RELE ALL @ 0,0 CLEA
  343. USE BAZA ORDE TAG KL REST FROM PRESS *prin = 3
  344. N=2*IST && ПОРЯДОК МАТРИЦЫ? '1ST, ЮТ, II, 12.13,14,15=', 1ST, ЮТД1,12,13,14,15
  345. DIME P (N), SIGST (N), SIG1(N), SIG2(N), SIGZ (N), SIGST (N), SUMST (N) STORO TO S, SIG1, SIG2,SIGZ, SIGST, SIGSTM, SUMST'ЧТЕНИЕ РЕШЕНИЙ ИЗ XI. OUT, ЖДИТЕ'
  346. HI= FOPEN ('XI.OUT', 2) IF HKO
  347. WAIT WIND 'Нет файла решения XI. OUT QUIT1. ENDIчтение решения в P ()
  348. FOR J=1 TO N && N порядок матрицы1.= FGETS (HI) && Store contents to memory каждое данное должно быть записано одной строкой P (J)=VAL (I)1. ENDF
  349. FCLOSE (hi)&& Close the file if prin > 2'ПРОВЕРКА РЕШЕНИЯ СИСТЕМЫ'
  350. DIME PP (N) && РР (1)=СУММА A (U)*P (J) STOR 0 TO PP
  351. FOR 1=1 TO N && ЗДЕСЬ ЭТО НОМЕР УРАВНЕНИЯ FOR J=1 TON
  352. FOR I=2*IOT+l TO 2*IST && DO 351. SS=INPUT ('A', I, J, 0,0,0)1. SIG2(J)=SIG2(J)-P (I)*SS1. ENDF1. ENDF *** 28 ***1. FOR J=1 TO IOT && DO 311. SIGZ (J)=0.1. FOR 1=1 TO 2*IOT && DO 311.=INT ((I+l)/2) IKN=2+INT (I/2)*2-I
  353. SIGZ (J)=SIGZ (J)+P (I) * INPUT ('S1', IKN, II, J,0,0)1. ENDF1. ENDF * * * ***
  354. FOR J=IOT+l TO 1ST && DO 321. SIGST (J)=0.
  355. FORI=2*IOT+l TO 2*IST && DO 321.=INT ((I+l)/2) IKN=2+INT (I/2)*2-I
  356. SIGST (J)=SIGST (J)+P (I) * INPUT ('S 1 ', IKN, II, J,0,0) SUMST (J)=SUMST (J)+INPUT ('Sr, IKN, II, J,0,0) * IIF (IKN=1, 0, 1)1. ПЕЧАТЬ СУММЫ SI1. DO PECH W IN PRESS WITH —
  357. FOR J=IOT+l TO 1ST && DO 32
  358. SIGSTM=MAX (SIGSTM, SIGST (J))1. ENDFstor 0 to SDL, SP2, nz, NOT, snz, snotdime natz (ist), natot (ist), snat (ist) FOR J=2 TO 2*11 STEP 2 && DO 33jj=j+2*iot
  359. SDL=SDL + DL (INT (J/2)) SP2=SP2 + SIG1(J) * DL (INT (J/2"storOto nz, NOT
  360. FOR 1=1 TO 2*ist if i≤2*iotnz = nz + INPUT ('A', I, jj, 0,0,0) * P (I)eise
  361. NOT=NOT+ INPUT ('A', I, jj, 0,0,0) * P (I)endi ENDF && 35 natz (j)=2*nznatot (j)=2*NOTsnat (j)=2 * (nz+NOT)snz=snz+2*nz*DL (INT (J/2)) snot=snot+2*NOT*DL (INT (J/2))1. ENDFnzm=snz/sdl notm=snot/sdl1. CZM = NZM / (NZM + NOTM)
  362. OSP2/SDL IF KONTUR = 50 RA = 1.1. ENDI1. CM = SP2 * IZS/ (PI * RA) vzm=nzm/cvotm=notm/c
  363. CI1=(ASIN (INPUTCX0', 1,11,1,0,0)/RA) + NU (I1)) * IZS /Р1() CSUM = C * CI1
  364. ROT=SQRT (X0(l, IOT, l,0,0)**2+X0(2,IOT, l,0,0)**2)
  365. CLAM=E/(2*RA*((RA**2+ROT**2)/(RA**2- ROT**2)+(RST**2+RA**2)/(RST**2ra**2)))*nat1. CFM=CM/CLAM
  366. SFN = 2 * SP2 * IZS CSTM = SIGSTM / SFN CNFN = NAT / SFN
  367. DO PECHW IN PRESS WITH ' '
  368. DO PECHW IN PRESS WITH ' Результаты расчета'
  369. DO PECH W IN PRESS WITH ' '1. DO PECHW IN PRESS WITH —
  370. С среднее значение по вершине конт давления'+ТКАМ (С,'@л') DO PECHJW IN PRESS WITH —
  371. CM средняя жесткость соединения '+TRAN (CM,'@A')1. DO PECHW IN PRESS WITH —
  372. CZM доля натяга вала '+TRAN (CZM,'@A')1. DO PECH W IN PRESS WITH —
  373. CLAM теоритеческая жесткость по Ляме '+TRAN (CLAM,'@A')1. PRIN > 1
  374. DO PECHW IN PRESS WITH — 'NZ, NOT нормальные прогибы на J ' + TRAN (NZ,'@A') + TRAN (NOT,'@A')
  375. DO PECHW IN PRESS WITH «DO PECHW IN PRESS WITH —
  376. Значения прогибов по эл-там J зуба и ступицы и их сумма’do pechjw in press with -1. J NATZ NATOT SNAT'
  377. FOR J=2 TO 2*11 STEP 2 && DO 33
  378. DO PECHW IN PRESS WITH — TRAN (J,'999')+' '+TRAN (NATZ (J),'999 999.99999')+' '+TRAN (NATOT (J),'999 999.99999')+' '+TRAN (SNAT (J),'999 999.99999') ENDFdo pechjw in press with «
  379. DO PECHJW IN PRESS WITH — 'J SIG1 SIG2 SIGZ / SIGST SIG1-SIG2' FOR J=1 TO ist && MAX (IOT, 2*IOT) IF J≤IOT && SIGZ
  380. DO PECFIW IN PRESS WITH — TRAN (J,'999')+' '+TRAN (SIG1 (J),'999 999.99999')+' '+TRAN (SIG2(J)/999 999.99999')+' '+TRAN (SIGZ (J),'999 999.99999')+' '+TRAN (SIG1(J)-SIG2(J),'999 999.99999') ELSE
  381. DO PECHW IN PRESS WITH — TRAN (J,'999')+' '+TRAN (SIG 1 (J),'999 999.99999')+' '+TRAN (SIG2(J),'999 999.99999')+' '+TRAN (SIGST (J),'999 999.99999')+' '+TRAN (SIG1(J)-SIG2(J),'999 999.99999') ENDIendf1. ENDI1. PRIN > 2do pech w in press with —
  382. SNZ площадь прогиба вершины 3y6a'+TRAN (SNZ,'@AI) DO PECH W IN PRESS WITH -snot площадь прогиба вершины cTynHn, bi'+tran (snot,'@a') do pech w in press with —
  383. NZM, NOTM средние натяги на вершинах вала и отверстия ступицыtr. an (nzm,'@a')+tran (notm-@a')do pechjw in press with —
  384. VZM, VOTM средним значения податливостей зуба и otb'+TRAN (VZM,'@a')+TRAN (YOTM,'@a')1671. DO PECH W IN PRESS WITH —
  385. SDL='+TRAN (SDL/@A')+' SP2='+TRAN (SP2,'@A')+' 0'+TRAN (C,'@A') DO PECHW IN PRESS WITH — 'CSUM='+TRAN (CSUM,'@A') DO PECHW IN PRESS WITH —
  386. CSTM-+TRAN (CSTM,'@A')+' CNFN='+TRAN (CNFN,'@A') DO PECHW IN PRESS WITH ' '
  387. DO PECHW IN PRESS WITH ' ' DO PECHW IN PRESS WITH — ' J P (J) J P (J) '1. FOR J-l TO ist1. DO PECHW IN PRESS WITH —
  388. R=SQRT (DX11+DX22) && 2*C0S (FI)
  389. R2=R*SQRT (DX11+DX22-CR2*0.68*NU*NU) && корректировка 9.8.95 R4=R2*R*R
  390. S=S-1./R4*(DX1*(A3*DX22-A2*DX11)*AI1+DX2*(A1 *DX22+A2*DX11)*AI2)
  391. СР2, UN = ', СР2, UN? 'CSU, CLAM = ', CSU, CLAM set de vi to sere set prin off set prin toclos all *
  392. ПРОГРАММА «Эталон-KONTROL»
  393. N=2* 1ST && ПОРЯДОК МАТРИЦЫ DIME P (N)
  394. DO PECHW in press WITH ' ПРОГРАММА PRESK '1. Mo PRINTG in press’КОНТРОЛЬ''ЧТЕНИЕ РЕШЕНИЙ ИЗ XI. OUT, ЖДИТЕ'
  395. HI= FOPEN ('XI.OUT', 2) IF HKO
  396. WAIT WIND 'Нет файла решения XI. OUT' QUIT1. ENDIчтение решения в P ()
  397. FOR 1=1 TO N && N порядок матрицы1.= FGETS (HI) && Store contents to memory каждое данное должно быть записано одной строкой P (J)=VAL (I)1. ENDF
  398. FCLOSE (hi)&& Close the file1. FOR N=1 TO 2 && F11. SN2(N)=0 IF N=11. Jl = l
  399. J2=IOT IF RA=1. XN2= 1. ELSE
  400. XN2=SQRT (XK (1,I, 1)**2+XK (2,I, 1)**2) ENDI && RA XN1=-XN2
  401. Jl=IOT+l J2=IST XN1=RA XN2=RST ENDI && N DLN=XN2 -XN1 FOR IZ=1 TO IZS SNZ2(IZ, N)=0 PIZ=PI*(IZ-1)*2./IZS SPZ1=0 SPZ2=0 FOR I=J1 TO J2 111=2*1−1 112=111+1 SNZ21=0 SNZ22=0 SA11=0 SA12=0 S A21 =0 SA22=01. FOR K=1 TO 2
  402. XX2=-INPUT ('XIZ2,I, IZ, K,0) X22=XX2*XX2
  403. SNZ2(IZ, N)=SNZ2(IZ, N)+SNZ21 * P (II 1)+SNZ22*P (II2)172
  404. SK2(N)=SN2(N)/(CM*RA*(3-N)) ENDF && N=1 TO 21. DO PECHW in press WITH —
  405. CK2(1)-+ TRAN (CK2(1),'@A')+' CK2(2)='+ TRAN (CK2(2),'@A')± 1 SN2(1)='+ TRAN (SN2(1),'@A')+' SN2(2)='+ TRAN (SN2(2),'@A')
  406. Pechatatb SN2(N), CK2(N) 4 znachenij**
  407. Конец PROG PRESKZ от 12 05 98quit1. ПРОГРАММА «KWADRAT"modi comm rezpr
  408. Ctrl + N вставить новую строку
  409. Ctrl + T удалить текущую строку
  410. Ctrl + W сохранить введеннуе строки1. Esc выход без сохраненияset reso off SET TALK OFF SET DELE ON SET EXAC ON SET DECI TO 12 set blin off
  411. SET SAFETY OFF && HE ВЫДАВАТЬ ПРЕДУПРЕЖДЕНИЯ ПО ZAP, INDEX ON,.1. SET PROC TO REZPR1. CLOS ALL1. CLEA ALL1. RELE ALLset date to dmy0,0 CLEAR1,0 SAY 'Программа обработки решений1
  412. PUBL ARRA A (5,6), AIL (5,5), BI (5), PII (5), QI (5), M (5), U (5), u (5), V (4) DELE FILE 'REZPR.out'
  413. DO PECHW WITH ' Программа обработки решений '+DTOC (DATE ())+' исходные данные: '1. SEle 1use rezpr orde tag nnjdo while .T.brow fiel j, z, 1, rs, rO, rf, h, d, rpackrein1. GO TOP1. ПРОВЕРКИ1. ERR = .F. IF ! EOF ()1. DO WHIL! EOF () IF Z 0 &&
  414. DO PECHJW WITH 'Размерность CHCTeMbi'+str (imax)
  415. DO PECHW WITH «DO PECHW WITH 'Индексы no i:'for i = 1 to imax
  416. DO PECHW WITH ' для индекса '+str (i) +' берется q с индексом '+str (u (i))next1. DO PECH^W WITH «
  417. DO PECHW WITH 'Для JO используется решение '+str (jjO) if v (3) > 0
  418. DO PECH W WITH 'Интервалы решений: с '+tran (v (l),'999') + ' no '+tran (v (2),'999') + ' и с ' +tran (v (3),'999') + ' по '+tran (v (4)/999') else
  419. DO WHILE! EOF () stor 0 to PII PII (1)=Z0/Z ZN = (1 LOJ IF ZN = 0
  420. WAIT WIND 'Знаменатель = 0 у P2 для j =' + str (j) quit1. ENDI
  421. PII (2) = ((1 L) / ZN) * L0 / L P1I (3) = (1 + R0**2 / RS0**2) / (1 + R**2 / RS ** 2) * РП (3) = 1/PII (3) if IMAX > 3
  422. ZN = (2 + H0 / LO J IF ZN = 0
  423. WAIT WIND 'Знаменатель = 0 у P5 для j =* + str (j) quit1. ENDI
  424. PII (4) = (2 + H / L) / ZN * PII (4) = 1/PII (4)1. PII (4) = (2+H) / (2±H0 Jendiif IMAX > 4
  425. ZN = (1 +(R0/RF)**2) IF ZN = 0
  426. WAIT WIND 'Знаменатель = 0 у P4 для j =' + str (j) quit1. ENDI
  427. BI (I) = BI (I) + Q0 * QI (U (I))1. ENDF1. FOR I = 1 TO IMAX1. FOR LL = 1 TO IMAX
  428. AIL (I, LL) = AIL (I, LL) + QI (U (I)) * QI (U (LL))1. ENDF1. ENDFj, 'p=', pii (l), pii (2), pii (3) — q= Q0, qi (l), qi (2), qi (3)
  429. REPL PI WITH PII (l), P2 WITH PII (2), P3 WITH PII (3), P4 WITH PII (4), P5 WITH
  430. REPL Ql WITH QI (1), Q2 WITH QI (2), Q3 WITH QI (3), Q4 WITH QI (4), Q5 WITH1. REPL QO WITH Q0 SKIP ENDD4=1. FOR I = 1 TO IMAXfor 11 = 1 to IMAX1. A (I, 11) = AIL (1,11)next1. NEXT1. FOR I = 1 TOIMAX1. A (I, IMAX + 1) = В 1(1)1. NEXT
  431. DO PECH^W WITH 'Расширенная матрица (A (l, l), A (l, 2),.):' FOR I = 1 TOIMAX ww = «for 11 = 1 to imax + 1ww = ww + TRAN (A (i, ll) ,'999.99 999')+' 'next1. DO PECHW WITH ww1. NEXT use
  432. СИСТЕМА СФОРМИРОВАНА publN
  433. N= IMAX && ЧИСЛО УРАВНЕНИЙ *suspdo GAUSS in GAUSS. prg with .F. && .Т. с выводом на экран * печать решения1. WW = «for 11 = 1 to IMAX
  434. WW = WW + TRAN (A (ll, imax + 1) ,'999.99 999')+' 'next
  435. DO PECHW WITH 'Найдено решение '+ww **** после решения1. USE REZ PR ORDE TAG NNJset flit to G ≥ v (l) .and. j <= v(2)) .or. 0>= v (3) .and. j ≤ v (4))sd = 0 && сумма Dj tc = 0 pp = 0 JMAX = 0jm = 0 && число решений TCKMAX = 0 go top1. DO PECHW WITH «
  436. DO PECHW WITH ' J PP D (D-PP)**2' do whil ! eof ()jm = jm + 1 sd = sd + (1. d) PII (l) =P1 PII (2) = P2
  437. PII (3) = P3 PII (4) = Р4 PII (5) = Р5ww = «www = «1. FOR I- 1 TO 3ww= ww + TRAN (PII (i) ,'999 999 999.99999')+' '1. NEXT1. FOR I = 4 TO 5www = www + TRAN (PII (i) ,'999 999 999.99999')+' '1. NEXT
  438. PP = D0 && произведен Pij ** К FOR I=1 TOIMAX
  439. PP = PP * (PII (U (I))**A (I, IMAX + 1))1. NEXT
  440. DO PECHW WITH tran (j,'999') + ' '+ tran (PP, '999.99 999') +' '+ tran (D, '999.99 999') + ' '+ tran ((D PP) ** 2, '999.9 999 999 999 999')
  441. ТС = TC + (D PP)**2 && сумма под большим корнемabc = ABS ((D PP) / (1. — D)) abc = abc*(l-D) REPL TCK WITH abc IF TCKMAX < abc1. TCKMAX = abc JMAX = J1. END I
  442. DO PECHW WITH 'j '+tran (j ,'9999') + ' p (i)=' +ww1. DO PECHW WITH ' ' +www
  443. DO PECHW WITH ' q (i)='+ tran (ql,'999 999 999.99999') + «+ tran (q2,'999 999 999.99999') + «+ tran (q3,'999 999 999.99999')
  444. DO PECH W WITH 'J max ='+ stiljmax) + * Tmax = '+ tran (TCKMAX, '9 999 999 999.99999') + ' TCK = '+ tran (tck, '9 999 999 999.99999')
  445. DO PECHW WITH 'Число использованных решений = '+ str (jm)go top browclose allquitretu
Заполнить форму текущей работой

Сессия? Спокойно!