Повышение эффективности механической лезвийной обработки на основе имитационного моделирования динамики технологической системы с учетом процесса стружкообразования
Диссертация
Появление возмущения в упругой системе приводит к изменению состояния деформированной зоны и к соответствующему изменению сил резания. Это изменение не может распространяться мгновенно на всю зону, что вызывает запаздывание в изменении сил резания. Наличие запаздывающих сил, возбуждающих замкнутую технологическую систему, может вызвать автоколебания в процессе резания. Потеря устойчивости… Читать ещё >
Список литературы
- H., Васильков Д. В., Вейц В. Л., Хитрик В. Э. Задачи динамики ГАП механообработки// Вибротехника. Межвуз.тематич. сб.науч.трудов. — Вильнюс, 1987. — 2(55). — С.73−83.
- Айзерман М.А. Теория автоматического регулирования.-М.:Наука, 1966.-452с.
- Амосов И.С. Осцилографические исследования вибраций причрезании металлов / Точность механической обработки и пути ее повышения//М.: Машгиз, 1951.
- Амосов И.С., Скраган В. А. Точность, вибрации и чистота поверхности при токарной обработке. М. — Л.: Машгиз, 1953. — 67 с.
- Армарего И.Дж.А. Браун Р. Х. Обработка металлов резанием/ Пер. с англ.- М.:Машиностроение, 1977. 325с.
- Бейлин И.LU., Вейц В. Л. Синтез параметров механической системы машинного агрегата / Зубчатые и червячные передачи: Некоторые вопросы кинематики, динамики, расчеты и производство.-Л.Машиностроение, 1974. С.267−285.
- Белламан Р. Введение в теорию матриц. М.:Наука, 1969. — 368 с.
- Блек У, Модель напряжения пластического течения при резании металла // Конструирование и технология машиностроения, 1979. № 4. — С.124 -139.
- Бродский А.Д., Кан.В. Л. Краткий справочник по математической обработке результатов измерений. М.:Стандартгиз, 1960. — 167с.
- Васильков Д.В. Оптимизация рабочих процессов на основе динамического моделирования технологической системы/ Машиностроение и автоматизация производства: Межвуз.Сборник. -Вып.2.-СПб.: СЗПИ.1996. С.93−103.
- Васильков Д.В. Теория и практика обеспечения стабильности и качества механической обработки маложестких заготовок/
- Машиностроение и автоматизация производства: Межвуз.Сборник. -Вып.З.-СПб.:СЗПИ, 1996.- С.54−76.
- Васильков Д.В. Теория и практика оптимизационного проектирования механической обработки маложестких заготовок: Дис.докт.техн.наук / СПб: ГТУ, 1997. 368с.
- Васильков Д.В., Вейц В. Л., Лонцих П. А. Динамика технологической системы при механической обработке маложестких заготовок.
- Иркутск: Изд-во Иркутск. ун-та. 1994.- 98с.
- Васильков Д.В., Вейц В. Л., Максаров В. В. К вопросу упрощения динамической модели технологической системы механической обработки // Машиностроение и автоматизация производства: Межвуз. сб. Вып. 14. СПб.: СЗПИ, 1998.-С.35−41.
- Васильков Д.В., Вейц В. Л., Максаров В. В. Моделирование процесса стружкообразования на основе кусочно-линейной аппроксимации // Академический вестник. Информатизация.Вып.1.-СПб.: СПбИмаш., 1998.-С.16−21.
- Васильков Д.В., Вейц В. Л., Шевченко B.C. Динамика технологической системы механической обработки. Санкт-Петербург: ТОО «Инвентекс», 1997−230с.
- Васильков Д.В., Козлова Е. Б. Обоснование выбора реологической моде-ли при решении нелинейных задач вязкоупругопластичности II Машинострое-ние и автоматизация производства: Межвуз. сб. Вып. 9. СПб.: СЗПИ, 1998. — С.36 — 44.
- Васильков Д.В., Козлова Е. Б., Максаров В. В. Анализ реологических уравнений для моделирования процесса резания// Машиностроение и автоматизация производства: Межвуз.сб.Вып.13.-СПб.: СЗПИ, 1999.-С.47−51.
- Васильков Д.В., Роменская Т. В. Вычислительные аспекты решения нелинейных задач динамики при исследовании контактных взаимодействий в технологической системе / Машиностроение иавтоматизация производства: межвуз.сб. Вып.8. СПб.: СЗПИ.1997. -С.74−82.
- Вейц В.Л. Динамика машинных агрегатов. Л.: Машиностроение, 1969.-370С.
- Вейц В.Л., Кочура А. Е. Динамика машинных агрегатов с двигателями внутреннего сгорания. Л.: Машиностроение, 1976. -384 с.
- Вейц В.Л.ч Кочура А. Е., Мартыненко A.M. Динамические расчеты приводов машин.-Л.: Машиностроение, 1971. 352 с.
- Вейц В.Л., Дондошанский В. К., Чиряев В. И. Вынужденные колебания в металлорежущих станках.-М.-Л.:Машгиз, 1959.-288с.
- Вейц В.Л., Журавлева Е. Ю. Метод кусочно-линейной аппроксимации в решении нелинейных задач динамики машин: часть I / Машиностроение и автоматизация производства.-СПб:Изд-во СЗПИ, 1999.-С.140−145.
- Вейц В. Л., Журавлева Е. Ю. Метод кусочно-линейной аппроксимации в решении нелинейных задач динамики машин: часть1./ Машиностроение и автоматизация производства.-СПб:Изд-во СЗПИ, 1999-С.126−134.
- Вейц В.Л., Козлова Е. Б., Максаров В. В. Математическое моделирование процесса стружкообразования при лезвийной обработке// Машиностроение и автоматизация производства: Межвуз.сб.Вып.14.-СПб.: СЗПИ, 1999. -С.70−77.
- Вейц В.Л., Максаров В. В. Локализация и неустойчивость пластической деформации в процессе стружкообразования при резании металлов // Машиностроение и автоматизация производства. Межвуз.сб.Вып.13.-СПб.: СЗПИ, 1999. С.39−43.
- Вейц В.Л., Максаров В. В. Модель формирования локализированных полос сдвига в зоне пластической деформациисрезаемого слоя // Машиностроение и автоматизация производства: Межвуз.сб. Вып. 14.-СПб.:СЗПИ, 1999. С.32−34.
- Вейц В.Л., Максаров В. В. Повышение устойчивости технологической системы при управлении реологическими параметрами процесса стружкообра-зования // Машиностроение и автоматизация производства: Межвуз. сб. Вып. 16. СПб: СЗПИ, 1999. -С.19 -29.
- Вейц В.Л.Ч Максаров В. В. Физические основы моделирования стружкообразования в процессе резания // Машиностроение и автоматизация производства: Межвуз. сб. Вып. 13. СПб.: СЗПИ, 1999.-С.44−46.
- Вейц В.Л., Максаров В. В. Динамическое моделирование стружкообразования в процессе резания // Машиностроение и автоматизация производства: Межвуз.сб. Вып.14.-СПШ:СЗПИ, 1999.-С. 15−20.
- Вейц В.Л., Мартыненко A.M. Автоколебания в механических кусочно линейных системах/ Нелинейные колебания и переходные процессы в машинах.-М.:Наука, 1972.- С.283−294.
- Глаговский Б.А., Московенко И. Б. Низкочастотные акустические методы контроля в машиностроении. Л.: Машиностроение, 1977. -208 с.
- Городецкий Ю.И. К теории возбуждения вибраций при резании металлов/ Сб.: Динамика машин. М.:Машгиз, 1963 — С.158−169
- Городецкий Ю.И. О колебаниях при резании металлов // Динамика систем. Межвуз. сб. Вып. 3. Горький: ГГУ, 1995. — С.58 — 89.
- Грановский Г. И., Грановский В. Г. Резание металлов. М.: Высшая шко-ла, 1985. — 304 с.
- Давиденков H.H. Динамические испытания металлов. М.: ОНТИ, 1936.-395 с.
- Детали и механизмы металлорежущих станков. Т1./ Под ред.Д. Н. Решетова. М. Машиностроение, 1972. — 664с.
- Джонсон У., Меллор П. Б. Теория пластичности для инженеров/ Пер. с англ.-М.:Машиностроение, 1979−567с.
- Диткин В.А., Придников А. П. Операционное исчисление. М.: Высшая школа, 1975. — 407 с.
- Дроздов H.A. К вопросу о вибрациях станка при токарной обработке//Станки и инструмент/, № 12, 1937
- Жарков И.Г. Вибрации при обработке лезвийным инструментом.-Л Машиностроение, 1986−184с.
- Заковоротный В.Л. Исследование динамической характеристики резания при автоколебаниях инструмента/Изв.техн.науки .- Ростов-на-Дону: Ростов. ин-тс.-х.машиностроения, 1976-С.37−44.
- Заковоротный В.Л. К теории управления динамикой резания// Системы управления металлорежущими станками и технологическими процессами. Ростов-на-Дону:СКНЦ, 1976.-С.22−27.
- Заре В.В. Вопросы самовозбуждения вибраций металлорежущих станков / Дис. докт. техн. наук. Л.: ЛПИ, 1972. — 238 с.
- Зорев H.H. Вопросы механики процесса резания металлов.- М.: Машгиз, 1956. 367 с.
- Ильницкий И.И. Колебания в металлорежущих станках и пути их устранения.-Москва-Свердловск:Машгиз, 1958−142с.
- Качанов Л.М. Основы теории пластичности.-М.:Наука, 1969.-480 с.
- Каширин А.И. Исследование вибраций при резании металлов. М.-Л.:АНСССР, 1944.-282с.
- Кедров С.С. Колебания металлорежущих станков. М.: Машинострое-ние, 1978. — 200 с.
- Колмогоров А.Н., Фомин C.B. Элементы теории функций и функционального анализа. М.:Наука, 1968. — 496с.
- Конструкционное демпфирование в неподвижных соединениях/ Калинин Н. Г., Лебедев Ю. А., Лебедева В. И. и др.- Под ред. Пановко Я.Г.- Рига: Изд-во АН ЛатССР, 1960.-169с.
- Конструкционное демпфирование в узлах вибрационных машин/Хвингия М.В., Цулая Г. Г., Гогилашвили В. Н. и др. -Тбилиси: Изд-во ГрПИ, 1973.-138с.
- Кочнева Л.Ф. Внутреннее трение в твердых телах при колебаниях. -М.: Наука, ^ 979.-96 с.
- Кривоухов В.А., Воронов А. Л. Высокочастотные вибрации резца при точении.М: Оборонгиз, 1956.
- Кудинов В.А. Автоколебания на низких и высоких частотах (устойчивость движения) при резании // СТИН, 1997, № 2. С. 16 — 22.
- Кудинов В.А. Динамика станков, М.:Машиностроение, 1967. — 359с.
- Кудинов В.А. Колебания в станках/Вибрации в технике: Справочник в 6-ти т. Ред. Челомей В.Н.- М.:Машиностроение, 1980.- т.З. Колебание машин, конструкций и их элементов/ Под ред. Диментберга Ф. М. и Колесникова К. С., 1980−544с.
- Кудинов В.А. Схема стружкообразования (динамическая модель про-цесса резания) // Станки и инструмент, 1992, № 10. С. 14 — 17, № 11. — С.26−29.
- Кудинов В.А. Теория вибраций при резании / Передовая технология машиностроения, АН СССР, 1955.
- Кучеряев Б.В. Механика сплошных сред. М.: МИСИС, 1999. -320с.
- Кучма Л.К. Устранение вибраций при обработке металлов резанием/ В кн.: Исследование колебаний металлорежущих станков при резании металлов.-М.: Машгиз, 1958.-С.158−219.
- Кучма Л.К. Учет сил сопротивления в автоколебательной системе станок-деталь-инструмент/Исследование колебаний металлорежущих станков при резании металлов.-М.: Машгиз, 1958.-С.220−227.
- Левин А.И. Математическое моделирование в исследовании и проектировании станков.-М.:Машиностроение, 1978.-184с.
- Левина З.М., Решетов Д. Н. Контактная жесткость машин.-М.Машиностроение, 1971.-264с.
- Максаров В.В. Реологическое представление при моделировании стружкообразования в процессе резания //Машиностроение и автоматизация производства: Межвуз.сб.Вып.14.-СПб.:СЗПИ, 1999.-С.21−24. ч
- Максаров В.В. Теория и методы моделирования и управления процессом стружкообразования при лезвийной механической обработке/СПб: ГТУ, 1999.-331 с.
- Манжос Г. А. Исследование вибраций в условиях скоростного точения и изыскание методов борьбы с ними/ Тоность механической обработки и пути ее повышения. Машгиз, 1951
- Мартынюк А.А. Техническая устойчивость в динамике. Киев: Техни-ка, 1973. — 188 с.
- Меркин Д.Р. Введение в теорию устойчивости движения. -М.:Наука, 1987. -304 с.
- Мокрицкий Б.Я. Особенности стружкообразования и устойчивости си-стемы при нестационарном резании // Изв. вузов. Машиностроение, 1996, № 4−6.-С.82 88.
- Мурашкин Л.С., Мурашкин С. Л. Прикладная нелинейная механика станков.-Л .".Машиностроение, 1977.-192с.
- Мурашкин С.Л. Вынужденные колебания самовозбуждающихся систем при вибрационной обработке материалов / Автоматизация и технология машиностроения. Труды ЛПИ, 1969. № 309. — С.234 -239.
- Мурашкин С.Л. Колебания и устойчивость движения систем станков с нелинейными характеристиками процесса резания / Дис. докт. техн. наук. Л.: ЛПИ, 1980. — 548 с.
- Неймарк Ю.И. Метод точечных отображений в теории нелинейных колебаний. М.: Наука, 1972. -472 с.
- Неймарк Ю.И. Об условиях самовозбуждения и стабилизации линеаризованных систем/Уч.зап. ГГУ, 1955-Вып.28
- Никитин Б.В. Расчет динамических характеристик металлорежущих станков.-М. :Машгиз, 1962.-112с.
- Остафьев В. А. Расчет динамической прочности режущего инструмента.^- М.: Машиностроение, 1979. 168с.
- Павлов А.Г. Эффективность снижения колебаний в станках // Вестник машиностроения, 1981 .-№ 7. С. 16 -18.
- Панин В.Е. Структурные уровни пластической деформации и разрушения. Новосибирск: Наука, 1990. — 251 с.
- Пановко Я.Г. Внутреннее трение при колебаниях упругих систем. -М.: ФИЗМАТГИЗ, 1960.-193с.
- Писаренко Г. С., Лебедев A.A. Сопротивление материалов деформированию и разрушению при сложном напряженном состоянии. Киев: Наукова Думка, 1969. — 208 с.
- Подураев В.Н. Автоматически регулируемые и комбинированные процессы резания.-М.Машиностроение, 1977.-304С.
- Пуш В. Э. Малые перемещения в станках.-М.:Машгиз, 1961- 351с.
- Расчет машиностроительных конструкций на прочность и жесткость/ Н. Н. Шапошников, Н. Д. Тарабасов, В. Б. Петров, В. И. Мяченков. М.:Машиностроение, 1981. — 333с.
- Решетов Д.Н., Левина З. М. Демпфирование колебаний в деталях станков/ Исследование колебаний металлорежущих станков при резании металлов. М.:Машгиз, 1958.-С.45−86.
- Роменская Т.В. Повышение производительности и точности при обработке резанием крупногабаритных маложестких заготовок / Дисс. канд. техн. наук. СПб.: ИМАШ, 1998. — 151с.
- Рыбин В.В. Большие пластические деформации и разрушение металлов. М.: Металлургия, 1986. — 224 с.
- Рыжков Д.И. Вибрации при резании металлов и методы их устранения.М.: Машгиз, 1961.
- Савин Г. Н., Рушицкий Я. Я. Элементы механики наследственных сред. К.: Вища школа, 1976. — 252 с.
- Садовский В.М. Разрывные решения в задачах динамики упругопластических сред. М.: Наука, 1997. — 208 с.
- Санкин Ю.Н. Динамика несущих систем металлорежущих станков. М. Машиностроение, 1986.-96 с.
- Скучик Е. Простые и сложные колебательные системы. М.: Мир, 1971.-557с.
- Соколовский А.П. Вибрации при работе на металлорежущих станках/ Исследование колебаний при резании металлов,-М.:Машгиз, 1958.-е. 15−18.
- Соколовский А.П., Жесткость в технологии машиностроения, Машгиз, 1946
- Соколовский В.В. Теория пластичности. М.: Высшая школа, 1969.-608с.
- Справочник по теории автоматического управления/ Под.ред.ААКрасовского. М: Наука, 1987. — 712с.
- Старков В.К. Обработка резанием. Управление стабильностью в автоматизированном производстве.- М. Машиностроение, 1972.-544С.
- Степанов М.Н. Статистические методы обработки результатов механических испытаний: Справочник.-М.: Машиностроение, 1985,-232с.
- Стренг Г. Линейная алгебра и ее применения/ Пер. с англ. -М.:Мир, 1980.-454с.
- Талантов Н.В. Физические основы процесса резания, изнашивания и разрушения инструмента. М.: Машиностроение, 1992. — 240 с.
- Тихонов А.Н., Васильева А.Б, Свешников А. Г. Дифференциальные уравнения.- М.: Наука, 1980. 232с.
- Тлустый И. Автоколебания в металлорежущих станках/ Пер. с чешек. М.:Машгиз, 1956.-395 с.
- Томович Р., Вукобратич М. Общая теория чувствительности / Пер, с сербск. И англ. М.:Изд-во «Советское радио», 1972. — 240с.
- Успенский В.А. Лекции о вычислимых функциях. М.:Физматгиз, 1960. -492с. v
- Хайрер Э., Нерсет С., Ваннер Г. Решение обыкновенных дифференциальных уравнений. Нежесткие задачи/ Пер. с англ. -М:Мир, 1995.-358с.
- Харкевич A.A. Избранные труды. В 3-х т. Т.2. Линейные и нелинейные системы. — М.:Наука, 1973.-С.253−379.
- Хильд А. Математическая статистика с техническими приложениями. М.:, 1936. — 385с.
- Циглер Г. Экспериментальные принципы термодинамики необратимых процессов и механика сплошной среды. М.: Мир, 1966. -135 с.
- Цыпкин Я.3. Основы теории автоматических систем. М.:Наука, 1977.-559с.
- Шенк X. Теория инженерного эксперимента / Пер. с англ. М.: Мир, 1972.-381 с.
- Шильман C.B. Метод производящих функций в теории динамических систем.-М.: Наука, 1978−336с.
- Штейнберг И.С., Устранение вибраций, возникающих при резании на токарном станке, Машгиз, 1947.
- Эльсгольц Л.Э. Дифференциальные уравнения и вариационное исчисление. М.:Наука, 1965.-424с.
- Эльясберг М.Е. Автоколебания металлорежущих станков: Теория и практика. С-Пб.:ОКБС, 1993.-182с.1.PLICIT REAL*8 (A-H.O-Z) !1. REAL*8 Y (3)
- EXTERNAL FCN, FCN1, FCN2,FCN3
- COMMON /MOD/XTEK, VS, BC, VK, P
- OPEN (8,FILE-DAN.DAT1,MODE-READ', STATUS-OLD')
- OPEN (9,FILE='REZ1.DAT', MODE=WRITE', STATUS='OLD')
- OPEN (7,FILE-REZ2.DAT1,MODE=WRITE', STATUS-OLD')1. XTEK=0.D01. X=0.D01. XEND=5.D01. H=0.01D01. HMAX=0.02D01. EPS=1.D-61. N=3
- Y (1)=1.D0 DO 1 1=2,3 Y (l)=0.D01.CONTINUE KUK=1
- READ (8,333)VS READ (8,333)BC 333 FORMAT (D15.8) P=0.D01. DO 111 111=1,500 N=31. A=P+0.13D01. X=0.0001DO
- CALL ODEX (N, FCN, X1Y, XEND, EPS, HMAX, H) IF (XTEKGE.XEND) GOTO 222 IF (XTEKGE.A) GOTO 21 X=0.0001D0
- CALL ODEX1 (N, FCN1, X, Y, XEND, EPS, HMAX, H) IF (XTEKGE.XEND) GOTO 222 IF (XTEKLE.A) GOTO 1121 P=XTEK A=P*2.D-222 X=0.0001D0
- CALL ODEX2 (N, FCN2, X)Y, XEND, EPS, HMAX, H) IF (XTEKGE.XEND) GOTO 222 IF (XTEKGE.A) GOTO 23 X=0.0001D0
- CALL ODEX3 (N, FCN3, X, Y, XEND, EPS, HMAX, H) IF (XTEKGE.XEND) GOTO 222 IF (XTEKLT.A) GOTO 2223 P=XTEK1. X=0.0001D0 !111 CONTINUE 222 STOP ENDC
- SUBROUTINE ODEX (N, FCN, X, Y, XEND, EPS1HMAX1H)1.PLICIT REAL*8 (A-H.O-Z)1.GICAL REJECT, LAST1. REAL*8 Y (N)1. EXTERNAL FCN
- COMMON/STAT/NFCN, NSTEP, NACCPT, NREJCT COMMON /EXJABL7 DZ (51), T (9,51), NJ (9)1HH (9), W (9)1ERR, FAC1 1 A (9), EPSD4, UROUND, FAC1, FAC2, SAFE2 DATA NJ/2,4,6,8,10,12,14,16,18/
- DATA A/3.D0,7.D0,13.DO, 21. DO, 31. DO, 43. DO, 57. DO, 73. DO, 91. DO/
- DATA NMAX/30 000/, KM/9/, UROUND/1.73D-18/
- DATA FAC1 /2. D-2/, FAC2/4. DO/, FAC3/. 9D0/, FAC4/. 8D0/
- DATA SAFE1/.65D0/.SAFE2/.94D0/1. EPSD4=EPS*SAFE 11. NSTEP=01. NREJCT=01. NACCPT=01. NFCN=0
- CALL MIDEX (KC, X, Y, HrHMAX, N, FCN) IF (ERR.GT.EPS) GO TO 100 60 X=X+H
- DO 70 1=1,N 70 Y (I)=T (1,I)
- H=HH (KC)*A (KOPT)/A (KC) END IF END IF K=KOPT GO TO 10 100 K=MIN0(K, KC)1. (KGT.2.AND.W (K-1).LT.W (K)*FAC3) K=K-1 NREJCT=NREJCT+1 H=HH (K) REJECT=.TRUE. GO TO 30 v 110 CONTINUE RETURN
- WRITE (6,*)' MORE THAN ', NMAX,' STEPS ' RETURN ENDC
- SUBROUTINE MIDEX (J, X, Y, H, HMAX, N, FCN) IMPLICIT REAL*8 (A-H.O-Z) EXTERNAL FCN
- CALL FCN (N, X+HJ*DFLOAT (MM), YH2, DY) DO 35 1=1, N YS=YH1(I) YH1(I)=YH2(I) 35 YH2(I)=YS+2.D0*HJ*DY (I) CALL FCN (N, X+H, YH2, DY) DO 40 1=1, N 40 T (J, I)=(YH1(I)+YH2(I)+HJ*DY (I))/2.D0 NFCN=NFCN+NJ (J) IF (J.EQ.1) RETURN DO 60 L=J, 2,-1
- FAC=(DFLOAT (NJ (J))/DFLOAT (NJ (L-1)))**2−1.DO DO 60 1=1,N
- T (L-1,I)=T (L, I)+(T (L, I)-T (L-1,I))/FAC
- CONTINUE ERR=O.DO DO 65 1=1,N
- SCAL=DMAX1 (DABS (Y (I)), DABS (T (1,1)), 1. D-6, UROUND/EPSD4) 65 ERR=ERR+((T (1, l)-T (2, l))/SCAL)**2 ERR=DSQRT (ERR/DFLOAT (N)) EXPO=1. DO/DFLOAT (2*J-1) FACMI N=FAC 1 **EXPO
- FAC=DMIN1(FAC2/FACMIN, DMAX1(FACMIN,(ERR/EPSD4)**EXPO/SAFE2))1. FAC=1.D0/FAC1. HH (J)=DMIN1(H*FAC, HMAX)1. W (J)=A (J)/HH (J)1. RETURN1. ENDC1. SUBROUTINE FCN (N, X, Y, F)1.PLICIT REAL*8 (A-H.O-Z)1. REAL*8 Y (N), F (N)
- COMMON /MOD/XTEK1VS, BC, VK, P1. A=16.D-61. B=2.4D-41. C=(6.D0A/S+1. D0)*1 .D-41. D=-0.053D0*BC1. F (1)=Y (2)
- F (2)=-Y (1)/A-B*Y (2)/A+Y (3)/A F (3)=D*Y (1)/C-Y (3)/C RETURN ENDC
- SUBROUTINE SOLOUT (NR, X, Y, N) IMPLICIT REAL*8 (A-H.O-Z) REAL*8 Y (N)
- COMMON /MOD/XTEK, VS, BC, VK, P KUK=1
- WRITE (6,*) (X+XTEK), Y (1), KUK WRITE (p,*) (X+XTEK), Y (1), KUK WRITE (7,*) Y (2), Y (3) VP=VS+Y (2) IF (VP.GE.VS) THEN XTEK=XTEK+X X=10.D0 END IF RETURN ENDC
- SUBROUTINE ODEX1 (N, FCN1, X, Y, XEND, EPS, HMAX, H)1.PLICIT REAL*8 (A-H.O-Z)1.GICAL REJECT, LAST1. REAL*8 Y (N)1. EXTERNAL FCN1
- COMMON /STAT/NFCN, NSTEP, NACCPT, NREJCT COMMON /EXTABL1/ DZ^IJ.T^.SIJ.NJ^.HH^J.W^.ERR.FAC, 1 A (9), EPSD4, UROUND, FAC1, FAC2. SAFE2 DATA NJ/2,4,6,8,10,12,14,16,18/
- DATA A/3.DO, 7. DO, 13. DO, 21. DO, 31. DO, 43. D0,57.DO, 73. DO, 91. DO/
- DATA NMAX/30 000/, KM/9/, UROUND/1.73D-18/
- DATA FAC ½. D-2/, FAC2/4. DO/, FAC3/.9D0/, FAC4/. 8D0/
- DATA SAF E 1/.65D0/, SAFE2/. 94D0/1. EPSD4=EPS*SAFE11. NSTEP=01. NREJCT=01. NACCPT=01. NFCN=0
- CALL MIDEX1(J, X, Y, H, HMAX, N, FCN1) 20 IF (J.GT. 1 .AND.ERR.LE.EPS) GO TO 60 GO TO 55 END IF 30 CONTINUE NSTEP=NSTEP+1 IF (NSTEP.GE.NMAX) GOTO 120 KC=K-1 DO 40 J=1,KC
- CALL MiDEX1(KC, X, Y, H, HMAX, N, FCN1) IF (ERR.GT.EPS) GO TO 100 60 X=X+H N
- DO 70 1=1, N 70 Y (I)=T (1,I)
- H=HH (KC)*A (KOPT)/A (KC) END IF END IF
- ERR=ERR+((T (1,I)-T (2,I))/SCAL)**2 !
- ERR=DSQRT (ERR/DFLOAT (N)) EXP0=1. DO/DFLOAT (2*J-1) F ACM I N=F AC 1 **EXPO
- FAC=DMIN1(FAC2/FACMIN, DMAX1(FACMIN,(ERR/EPSD4)**EXPO/SAFE2))1. FAC=1.D0/FAC1. HH (J)=DMIN1(H*FAC, HMAX)1. W (J)=A (J)/HH (J)1. RETURN1. ENDC1. SUBROUTINEfCN1(N, X, Y, F)1.PLICIT REAL*8 (A-H.O-Z)1. REAL*8 Y (N), F (N)
- COMMON /MOD/XTEK, VS, BC, VK, P1. A=15.5D-61. B=4.9D-4
- C=(6. D0/VS+1. DO)* 1. D-41. D=-0.051D0*BC1. F (1)=Y (2)
- F (2)=-Y (1)/A-B*Y (2)/A+Y (3)/A F (3)=D*Y (1)/C-Y (3)/C RETURN ENDC
- SUBROUTINE SOLOUT1(NR, X, Y, N) IMPLICIT REAL*8 (A-H.O-Z) REAL*8 Y (N)
- COMMON /MOD/XTEK, VS, BC, VK, P1. KUK=21. A=P+0.13DO
- WRITE (6,*) (X+XTEK), Y (1), KUK WRITE (9,*) (X+XTEK), Y (1), KUK WRITE (7,*) Y (2), Y (3) VP=VS+Y (2) IF (VP.LT.VS) THEN XTEK=XTEK+X X=10.D0 END IF RETURN END
- SUBROUTINE ODEX2 (N, FCN21X1Y, XEND, EPS, HMAX1H) IMPLICIT REAL*8 (A-H.O-Z) LOGICAL REJECT, LAST1. REAL*8 Y (N) !1. EXTERNAL FCN2
- COMMON /STAT/NFCN, NSTEP, NACCPT, NREJCT COMMON /EXTABL2/
- A (9), EPSD4, UROUND, FAC1, FAC2, SAFE2 DATA NJ/2,4,6,8,10,12,14,16,18/
- DATA A/3.D0,7.00,13.00,21.00,31.00,43.00,57.00,73.00,91. DO/
- DATA N MAX/30 000/, KM/9/, UROUND/1.73D-18/
- DATA F AC½. D-2/, F AC2/4. DO/, F AC3/. 9D0/, FAC4/. 8D0/
- DATA SAFE1/.65D0/, SAFE2/.94D0/1. EPSD4=EPS*SAFE11. NSTEP=0 v1. NREJCT=01. NACCPT=01. NFCN=0
- H=HH (K) REJECT=.TRUE. GO TO 30 110 CONTINUE RETURN
- WRITE (6,*)' MORE THAN ', NMAX,' STEPS ' RETURN END
- SUBROUTINE MIDEX2(J1X, Y, H, HMAX, N1FCN2) IMPLICIT RE^L*8 (A-H, 0-Z) EXTERNAL FCN2
- CALL FCN2(N, X+HJ*DFLOAT (MM), YH2, DY) DO 35 1=1, N YS=YH1(I) YH1(I)=YH2(I) 35 YH2(I)=YS+2. D0*HJ*DY (I) CALL FCN2(N, X+H, YH2, DY) DO 40 1=1, N 40 T (J, l)=(YH1 (I)+YH2(l)+HJ*DY (l))/2. DO NFCN=NFCN+NJ (J) IF (J.EQ.1) RETURN DO 60 L=J, 2,-1
- FAC=(DFLOAT (NJ (J))/DFLOAT (NJ (L-1)))**2−1.DO DO 60 l=1,N
- T (L-1,l)=T (L, l)+(T (L, l)-T (L-1,l)yFAC 60 CONTINUE ERR=0.D0 DO 65 1=1,N
- SCAL=DM AX1 (DABS (Y (I)), DABS (T (1,1)), 1, D-6,UROUND/EPSD4) 65 ERR=ERR+((T (1, l)-T (2, l))/SCAL)"2 ERR=DSQRT (ERR/DFLOAT (N)) EXPO=1. DO/DFLOAT (2*J-1) FACMIN=FAC1**EXPO
- FAC=DMIN1(FAC2/FACMIN, DMAX1(FACMINl (ERR/EPSD4)**EXPO/SAFE2))1. FAC=1.D0/FAC1. HH (J)=DMIN1(H*FAC, HMAX)1. W (J)=A (J)/HH (J)1. RETURN1. END
- SUBROUTINE FCN2(N, X, Y, F) IMPLICIT REAL*8 (A-H.O-Z) REAL*8 Y (N), F (N)
- COMMON /MOD/XTEK, VS, BC, VK, P1. A=16.D-6 s1. B=2.4D-41. C=(16.D0/VS+0.75D0)*1.D-41. D=-0.089D0*BC1. F (1)=Y (2)
- F (2)=-Y (1)/A-B*Y (2)/A+Y (3)/A F (3)=D*Y (1)/C-Y (3)/C RETURN END
- SUBROUTINE SOLOUT2(NR, X, Y, N) IMPLICIT REAL*8 (A-H.O-Z) REAL*8 Y (N)
- COMMON /MOD/XTEK, VS, BC, VK, P KUK=3
- WRITE (6,*) (X+XTEK), Y (1), KUK WRITE (9,*) (X+XTEK), Y (1), KUK WRITE (7,*) Y (2), Y (3) VP=VS+Y (2) IF (VP.GE.VS) THEN XTEK=XTEK+X X=10.D0 END IF RETURN END
- SUBROUTINE ODEX3 (N, FCN3, X, Y, XEND, EPS, HMAX, H)1.PLICIT REAL*8 (A-H.O-Z)1.GICAL REJECT, LAST1. REAL*8 Y (N)1. EXTERNAL FCN3
- COMMON/STAT/NFCN, NSTEP, NACCPT, NREJCT COMMON /EXTABL3/ DZ (51), T (9,51), NJ (9), HH (9), W (9), ERR, FAC,
- A (9), EPSD4, UROUND, FAC1, FAC2,SAFE2 DATA NJ/2,4,6,8,10,12,14,16,18/
- DATA A/3. DO, 7. DO, 13. DO, 21. DO, 31. DO, 43. DO, 57. DO, 73. DO, 91. DO/
- DATA NM AX/30 000/, KM/9/, UROUND/1.73D-18/
- DATA FAC½. D-2/, FAC2/4. DO/, FAC3/. 9D0/, FAC4/. 8D0/
- DATA SAFE 1/.65D0/, SAFE2/.94D0/1. EPSD4=EPS*SAFE11. NSTEP=01. NREJCT=0 '1. NACCPT=01. NFCN=0
- COMMON /MOD/XTEK, VS, BC, VK, P1. A=15.5D-61. B=4.9D-4
- C=(16. D0A/S+0.75D0)*1. D-4 D=-0.086D0*BC F (1)=Y (2) v
- F (2)=-Y (1)/A-B*Y (2)/A+Y (3)/A F (3)=D*Y (1)/C-Y (3)/C RETURN END
- SUBROUTINE SOLOUT3(NR, X, Y, N) IMPLICIT REAL*8 (A-H.O-Z) REAL*8 Y (N)
- COMMON /MOD/XTEK, VS, BC, VK, P KUK=4
- WRITE (6,*) (X+XTEK), Y (1), KU K WRITE (9,*) (X+XTEK), Y (1), KUK WRITE (7,*) Y (?), Y (3) VP=VS+Y (2) IF (VP.LT.VS) THEN XTEK=XTEK+X X=10.D0 END IF RETURN END
- УТВЕРЖДАЮ Директор НИИ «Энергостальп1. ВА. Повышевг
- От СПбИМаш: ' От НИИ Знергосталь:
- Д. т. н., профессор К.т.н. > <С5>. научн. сотр.1. К.т.н. доцент1. В.Л.Вейц1. В.В.Цуканов1. Д.В.Васильков
- УТВЕРЖДАЮ» ^^Рчальныи директор АО «НИТИ к. э. н1. Б.И.Катенев1. ТЕХНИЧЕСКИЙ АКТ
- Главный инженер АО «НИТИ Энергомаш"1. Н.Н.Ревин1. Инженерная академияотделение Машиностроение и инженерная механикаN
- Проректору по научной работе СПбИМаш академическому советнику РИА Ю. А. Державцуг 1
- Об использовании результатов научных исследований
- Перечисленные разработки включены в перечень передовых НИОКР для внедрения на предприятиях Санкт-Петербурга и Северо-Западного региона РФ под эгидой отделения «Машиностроение и инженерная механика».
- Главный ученый секретарь, академик
- Академик-секретарь отделения «Машиностроение и инженерная механика» Санкт-Петербургской Инженерной академии1. Ю.И. Мазуренко14.641. В. Смирнов1. Ректору
- Санкт-Петербургского института машиностроения М.А.Мартынову
- Генеральный директор фирмы «БИ ПИТРОН», к.т.н.