Экономико-математическое моделирование автоматизированного проектирования трудоемкости производственнотехнологических процессов организации сложных изделий
Автор: Чоракаев Олег Эдуардович, Попов Илья Валерьевич, Верушкин Олег Александрович
Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc
Рубрика: Механика и машиностроение
Статья в выпуске: 4-3 т.14, 2012 года.
Бесплатный доступ
В статье авторы разбирают численный пример экономико-математического моделирования автоматизированного проектирования системы и производственно – технологических процессов организации производства сложных изделий, в том числе элементов авиационной техники и представляет графоаналитическую интерпретацию состояния системы в зависимости от вложения средств по функции V = V(C), где V – вложения по стоимости С; приводят доказательство по предложенному методу и точности этого метода технико-экономического и экономико–математического моделирования производственно – технологических и организационных процессов организации производства САПР технико–экономического планирования (ТЭП) (САПРТЭП).
Экономико-математическое моделирование, автоматизированное проектирование, производственно-технологические процессы
Короткий адрес: https://sciup.org/148201306
IDR: 148201306
Текст научной статьи Экономико-математическое моделирование автоматизированного проектирования трудоемкости производственнотехнологических процессов организации сложных изделий
Адаптацию метода технико-экономического и экономико-математического моделирования (описания) САПРТЭП по методологии ФСИ проведем на численном примере (элементах авиационной техники) вложения средств (затрат) с получением дохода в развитие и совершенствование названной системы.
Итак: 1. Пусть даны две пары чисел C 0 = 874, V 0 = 175 (начальная точка) и C = 103423 д.е., V = 4370 д.е. (конечная точка), определяющие соответственно связь между минимальными и максимальными величинами затрат и доходов. Считая, что в САПРТЭП поддерживается потенциальный уровень средней цены сборочных работ (и вычислительных услуг) при наличии критерия оптимальности
C
V
1 0 = j x ln( к 1 x ) dx + j y ln( к 2 y ) dy , то есть име-
C
V 0
ется система, когда устанавливается зависимость V = V ( C ) , описывающая непрерывное изменение системы при переходе от начальной к конечной точке, обеспечивающих наиболее быстрое увеличение значений численной характеристики. Далее показано, что искомая зависимость будет иметь степенной вид.
Поскольку заданы две точки, степень У и
константа А определяются однозначно: у = ln( V 0 / V )/ln ( C 0/ C ) = 0,674, A = V 0 / Е 0 у = 1,821. Тогда из соотношения вытекает, что к 2 = 0,0084 определяется как корень уравнения к 2 ln у ( к 2 V 0 ) = у У ]v 0 и далее к 1 = у/ C 0 ln ( к 2 V 0 ) = 0,002 Поэтому искомая зависимость – степенная и имеет вид V 1 = 1,821 • C 0,674 , а соответствующие средние цены равны C 1 = V к 1 = 495,7, C 2 = 1 к 2 = 119,3.
Далее, 2. Пусть в добавление к нижней и верхней точкам к численному примеру по п. 1 имеется набор точек (табл. 1), где величины эффекта С и соответствующего ему вознаграждения V; Vп запланированные величины, а V1 и V2 – величины, рассчитанные по аппроксимирующим V зависимостям V 1 = 1,821 • C 0,674 и V 2 = 2,83 • C 0,64. В этом случае не существует степенной зависимости, на графике которой лежат данные точки (рис. 1). Далее найдем степенную зависимость, которая была бы по возможности наиболее близка к набору точек, приведенных в табл. 1. Выбирая в качестве меры близости сумму квадратичных отклонений по методу наименьших квадратов, получим V 2 = 2,83 • C 0,64. Значит, C 0 = 874, V 0 = 216, и далее к 2 = 0,0067 определяется как корень уравнения к 2 ln у ( к 2 V 0 ) = уУ Iv 0 ( при у = 0,64 ) Наконец, к j = у C^ 0 ln ( к 2 V 0 ) = 0,002 , а соответствующие средние цены равны C 1 = ^ к 1 = 495,7 , а C 2 = 1/ к 2 = 150,3.
В частности, оптимальное отношение цены механосборочной и слесарной работы к цене
Таблица 1. Набор точек величины эффекта C и вознаграждения V
С=100000 соответственно. На рис. 1 задается множество кривых, каждая из которых определяет множество состояний САПРТЭП организации производства машиностроительного предприятия, имеющих одно и то же значение технико-экономической информационной характеристики системы. При увеличении средних значений затрат примерно в 1,5 раза и дохода в 3 раза [нижняя и верхняя линии уровня – “потенциальные возможности” САПРТЭП], ее “сложность” увеличивается примерно в 5 раз (рис. 1, б).
Таким образом, в рамках рассматриваемой модели САПРТЭП по схеме “затраты-доход” полученные результаты позволяют определять оптимальные пропорции между ее составляющими (в комплексе с промышленными АСУ и САПРТЭП) и оценивать “потенциальные возможности” состояния САПРТЭП в смысле характеристики, связанной в силу способа ее построения с множеством реализаций. Следует отметить, что хотя величины затрат, дохода (см. табл. 1), работ (сборочных), вычислительных услуг и измеряются в условных единицах (д.е.), эти данные основаны на реальном исследовании сборочного производства авиастроительного предприятия. Величины С и VП имеют существенное значение, поскольку при их разработке явно или неявно использовались многочисленные содержательные свойства исследуемой организации производства машиностроительного предприятия как большой экономической системы при выполнении госбюджетных НИР и ОКР. При этом в примере по п.1 соответствует случай когда, по некоторым неформальным соображениям, максимальным и минимальным уровням затрат и доходов придается на авиастроительном предприятии решающее значение. В примере по п.2 отражен случай, когда все дан-

б)
а)

Рис. 1. Графоаналитическая интерпретация состояния системы от вложений по функциональной зависимости V = V ( С ) : а – линия уровня функции V=V(C); б – график функции V=V(C)
ные таблицы 1 имеют значение, причем значимость каждой точки учитывается в той или иной мере в зависимости от поставленной цели за счет выбора соответствующей аппроксимации.
На основании вышеизложенного, приведем соответствующее аналитическое доказательство по п.п. 1 и 2 примера.
1. По п. 1 примера доказательство вытекает из свойств интеграла, ограничений
X —■ +^ , то у ( x ) -Л X при x — +^ . Докажем.
Доказательство. Пусть заключение предложения неверно. Тогда найдется некоторое положительное 5 и последовательности X — +“ У п — +^ У — у ( x„ \ x „ — x ( у „ )] такие, что
xn yn
1 > 5, y^ -1 > 5. xn
1 ( 1 ^
C ( S ) < C , J ln f ( S ) dS < ln j f ( S ) dS I , где 0 1 0 )
f ( S ) почти всюду неотрицательная суммируемая на отрезке [0,1] функция. Покажем справедливость этого неравенства. Если f ( S ) равна нулю на некотором множестве положительной меры, то левая часть неравенства равна - ^ , и тогда неравенство верно. В противном случае для почти всех S е [ 0,1 ] в силу вогнутости логарифма получим, (касательная лежит выше графика вогнутой функции рис. 1):
ln f ( S ) < ln J f ( S ) dS

f ( S ) - J f ( S ) dS
Интегрируя это неравенство по отрезку [0,1], получим искомое неравенство, что и требовалось доказать.
-
2. По п. 2 примера 2.1. Действительно, рассмотрим, следующее семейство функций ( e > 0): C E ( S ) — ( 2 e ) - 1 ■ C , B E — ( 2 e ) ■ B , u f ( S ) = ( 2 e )- 1 ■ U , V ( S ) — ( 2 e )- 1 ■ V при S 0 - E < S < S 0 + E и C E ( S ) — 0, B E ( S ) — 0, U E ( S ) — 0, V ( S ) — 0 для остальных S е [0,1], где S 0 - некоторое число из интервала (0,1). Тогда непосредственно проверяется, что функции этого семейства удовлетворяют соответствующим ограничениям и Ф [ C e ( S )> B E ( S )> V E ( S )> U E ( S )] — “ при E — 0.
-
2.2. Сформируем и докажем два вспомогательных утверждения:
-
2.2Л . Пусть f ( t ) – монотонная положительная непрерывная функция времени, определен-
- z
ная при t > X0. Положим F(z) — J f (t) dt при x 0
-
z > X 0 и пусть
F(z) — +^ при z — +^ .
Пусть кроме того, некоторого R > 0 при всех z > X 0 выполняется неравенство F ( z ) < R ■ zf ( z ). Предположим далее, что каждому х однозначно сопоставлен некоторый y — y ( x ) [соответственно обратная зависимость: x — f ( у ) , причем y ( x ) —^ при X —+^ ]. Тогда, если F ( x ) ^ F [ у ( x )] при
По теореме о среднем для интеграла и в силу f (t) получим
монотонности
- |< | F (у„)-F (x„ )
" Х" min [ f ( x„ ) , f ( у „ )] . (2)
Если f ( x n )— min [ f ( x n ), f ( у п )] , то неравенство (2) можно переписать в виде
I Уп - x„ | < F ^ y V^ V JLXnJ . ■ F ( x „ )
Xn ” F ( Xn ) Xnf ( Xn ) . (3)
Если f ( y „ ) — min [ f ( X„ ), f ( y „ ) ] , то неравенство (2) принимает вид
I У п - Xn | < I F ( У п )- F ( Xn )| . F ( У п ) y „ F ( У п ) У п У ( У п ) . ( )
Но по условию п. 2.2.1 правые части неравенства (2)-(4) стремятся к нулю при п —> +^ , что противоречит определимости от нуля левых частей этих неравенств в силу соотношений (1).
-
2. 2.2. Для функции f ( t ) — , некоторо- V ln t
-
3. Пусть B — B ( t ), C — C ( t ) , U — U ( t ) и V — v ( t ) – решение системы и после интегрирования этих уравнений получим соотношения:
го X 0 > 1 и построенной по ним функции F ( z ) , выполняются сделанные в п. 2.2.1 допущения. Докажем это.
Доказательство. Свойства монотонности, положительности и непрерывности функции f ( X ) очевидны, а условие F ( z )——+^ при z —— +^ выполняется в силу расходимости соответствующего интеграла. Требуемая константа R > 0 имеет место в силу непрерывности и существования конечного предела соответствующей функции, который (интеграл) вычисляется непосредственно с помощью правила Лапиталя:
F ( z ) г f ( z )
lim —— lim / x / x =1 z—+“ zf (z) z—+“ f (z)+ zf'(z)
C2 — B2 lnB -1 в2 + R1,(5)
V2 — U2 lnU -1U2 + R2, где Ci(i = 1,2)- постоянные, определяемые по формулам
2 2 12 2 2 1 2
R1 C 0 B 0 ln B 0 2 B 0 , R 2 V 0 U 0 ln U 0 2 U 0 "
Используя (5)-(6) получим:
B xdx в0 x2 In x — 0,5x2 + R1
= t,
P 2 = 4 z ln ( kz ) - 4 z + 4 R , (11) где R – постоянная, в силу начальных условий и уравнения (9), задается формулой:
R = ( C 0 + V 3 ) 2 - U 2ln ( kU 0 2 ) + U 0 2" (12)
Следовательно, из (11) вытекает, что:
U 2
J
dz
z ln( kz)- z + R
= 2 1 .
U xdx
J I =^ = t"
и 0 V x 2 In x — 0,5 x 2 + R 2
Соотношения (5)-(7) полностью задают искомое решение. Установим требуемые свойства этого решения. Из (7) вытекает, что
B xdx U xdx
B 0 J x 2 ln x - 0,5 x 2 + R 1 U 0 x 2 ln x - 0,5 x 2 + R 2 (8)
Здесь нетрудно видеть, что подкоренные выражения слева и справа в (8) – возрастающие непрерывные функции переменной х при x > B0 и x > U0, равные x = B0 и x = U0 значениям C0 и V0 соответственно. Поэтому левая и правая части выражения (8) – непрерывные возрастающие функции переменных B и U. Множествами значений этих функций являются все неотрицательные числа в силу несуммируемости подынтегральных выражений в (8) (они эквивалентны функции 1/д/lnx при x ——+^ ). Следовательно соотношение (8) взаимооднозначно определяет переменные В и U, причем
В силу условия k > 1/ U 0 подкоренное выражение слева (13) возрастающая непрерывная функция и при z > U 0 , равная при z = U 0 значению ( C 0 + V 0 ) 2. Поэтому (13) однозначно определяет (см. аналогичные рассуждения в п. 2.3) искомое решение U = U ( t ) , причем U ( t ) — +^ при t —— +^. Далее получим решение C = C"tt ) , V = V ( t ) в виде
C ( t ) = C 0 + J ln [ U ( x ) ] dx , V ( t ) = V 0 + J ln [ kU ( x ) ] dx . (14) 00
Соотношения (13)-(14) полностью задают искомое решение. Из (14) вытекает, что C ( t ) — +^ и V ( t ) — +^ при t — +~ , а в силу эквивалентности ln U Л ln ( kU ) при U — +^ получим C Л V .
-
5. Из системы уравнений вытекает, что искомое решение C = C ( t ), V = V ( t ) определяется соотношениями:
C dx V dx
—7---x--- = t , ----x--- = t . (15)
C 0ln(k 1x )+ 1 Jln(k 2 x )+ 1
в — +^ и U — +^
Из (8) вытекает:
B
I
max( C 0 , B ) ) + 1

- Л
U ( B )
I
max ( C 0 , B 0 ) + 1
dx ln x
при B — +^ , где —< - знак предопределения.
Тогда, по п.п. 2.2.1 и 2.2.2, получим B Л U (при этом B — +^, U — +^) и далее отсюда и из (5)-(6) устанавливаем, кроме того, что и C Л V (при C — +^, V — +^ ), что и требовалось доказать.
4. Пусть C = C ( t ) , U = U ( t ) и V = V ( t ) -решение системы, тогда имеем соответственно:
(C + V )= ln( kU )2, ^ U2 J = 2(C + V).
Отсюда
( v
V U 2 V
= 2ln ( kU 2 ) .
Полагая z = U2 и делая стандартную замену z = р ( z ) из (10) после понижения порядка уравнения и последующего интегрирования получим
Рассуждая так же, как в п. 3, с применением п.п. 2.2.1 и 2.2.2 получим, что C —— ^, V —— ^ и C Л V .
Следовательно, после анализа численного примера с доказательствами по п.п. 1, 2, 3, 4 и 5, имеется основание полагать, что разработанный экономико-математический метод описания САПРТЭП в составе комплексной системы опытного производства исследуемого предприятия, логичен, жизнеспособен и функционально применим для описания технико-экономических процессов опытного производства, а также может быть кооптирован на любые комплексные автоматизированные управляющие системы отрасли машиностроения и электронного приборостроения страны.
На основании вышеизложенного, необходимо отметить, что:
1. Существующие классические методы технико-экономического и экономико-математического моделирования производственно-технологических и организационных процессов организации производства с использованием САПРТЭП обеспечивают ряд вопросов модели-
рования , но не всегда эффективны из-за относительной сложности учета ограничений , налагаемых на изменение переменных , и необходимости решения двухточечных краевых задач с использованием программных комплексов и ЭВМ, что значительно затрудняется высокой стоимостью функции расчетов.
-
2. Проведенный анализ и исследование технико-экономических и экономико-математических методов описания САПРТЭП, а также разработанная методика технико-экономического и экономико-математического моделирования производственно-технологических и организационных процессов с использованием САПРТЭП, основанная на принципе технологии «затраты-до-ход» и «вложения-доход» по критериям функциональности и стоимости на уровне функции при верификации разработок, дают основание полагать, что разработанные методы и методика описания САПРТЭП наиболее полно обеспечивают взаимоувязку компонентов системы по функции (как полезному свойству, состоянию или действию) как с позиции технико-экономического и экономико-математического моделирования, так и математического интерпретирования и технико-экономического описания процессов организации производства технических систем.
-
3. Использование системы вложений по методике Канбан – «затраты-доход» и «вложения-
- доход», позволяет уяснить принцип «затратности» любой исследовательской работы, которая может принести доход, доход весьма весомый, а может и остаться в виде затрат, без положительного экономического эффекта! Это позволяет исследователям при проведении НИР и ОКР настраиваться в направлении получения обязательного экономического эффекта разработок и экспериментов, что является положительным моментом «затратной» системы исследований.
Список литературы Экономико-математическое моделирование автоматизированного проектирования трудоемкости производственнотехнологических процессов организации сложных изделий
- Ляшко Ф.Е. Методика проведения расчетов экономической эффективности работы сборочного производства самолетов на основе математического моделирования процессов в САПР, АСУТП, АСТПП комплексной АСУП//Проблемы машиностроения и автоматизации. Международный журнал. №2/2008. С.34-39
- Маркорян И.Н. Разработка моделей и методов снижения трудоемкости сборочных процессов авиационного производства на основе функционально-стоимостного анализа. Дисс. … канд. экон. наук. 08.00.13. -Математические и инструментальные методы экономики. Ульяновск, 2006. С.164.
- Попов П.М. Оптимизация технических решений автоматизированного проектирования и управления разработками для комплексного повышения эффективности подготовки машиностроительных производств. Дисс. … докт. техн. наук. 05.13.12 -«Системы автоматизации проектирования (по техническим наукам -промышленность)». Ульяновск. УлГТУ, 2001. С. 288.