Все выпуски

Гибкие тканые композиты относят к классу высокотехнологичных инновационных материалов. Благодаря совмещению различных компонентов наполнителя и элементов армирования такие материалы применяют в строительстве, оборонной промышленности, судо- и авиастроении и др. В отечественной литературе уделено недостаточное внимание к тканым композитам, изменяющим свою геометрическую структуру армирующего слоя в процессе деформирования. В настоящей работе приводится анализ предложенного ранее комплексного подхода к моделированию поведения гибких тканых композитов при статическом одноосном растяжении для дальнейшего обобщения подхода на двухосное растяжение. Работа нацелена на качественное и количественное описание механических деформационных процессов, протекающих в структуре исследуемых материалов при растяжении, к которым относится распрямления нитей армирующего слоя и увеличение величины взаимного надавливания накрест лежащих нитей армирования. В начале процесса деформирования распрямление нитей и увеличение взаимного надавливания нитей наиболее интенсивны. С увеличением уровня нагрузки изменение указанных параметров замедляется. Например, изгиб нитей армирования переходит в центральное растяжение, а величина нагрузки от взаимного надавливания более не увеличивается (стремится к константе). Для моделирования описанных процессов вводятся основные геометрические и механические параметры материала, влияющие на процесс формоизменения, приводятся необходимая терминология и описание характеристик. В связи с высокой геометрической нелинейностью все процессы описаны в приращениях, так как на начальных значениях нагрузки происходит значительное формоизменение армирующего слоя. Для количественного и качественного описания механических деформационных процессов, протекающих в армирующем слое, выведены аналитические зависимости, позволяющие определить приращение угла распрямления нитей армирования и нагрузки, вызванной взаимным надавливанием накрест лежащих нитей на каждом шаге приращения нагрузки. Для апробации выведенных зависимостей приведен пример их применения для гибких тканых композиционных материалов марок VP4126, VP6131 и VP6545. Результаты моделирования подтвердили предположения о процессах выпрямления нитей и замедления увеличения взаимного надавливания нитей. Приведенные в данной работе результаты и зависимости имеют непосредственное отношение к дальнейшему обобщению предложенных ранее аналитических моделей для двухосного растяжения, так как растяжение в двух направлениях существенно уменьшит выпрямление нитей и увеличит величину взаимного надавливания при аналогичных нагрузках.

Flexible woven composites are classified as high-tech innovative materials. Due to the combination of various components of the filler and reinforcement elements, such materials are used in construction, in the defense industry, in shipbuilding and aircraft construction, etc. In the domestic literature, insufficient attention is paid to woven composites that change their geometric structure of the reinforcing layer during deformation. This paper presents an analysis of the previously proposed complex approach to modeling the behavior of flexible woven composites under static uniaxial tension for further generalization of the approach to biaxial tension. The work is aimed at qualitative and quantitative description of mechanical deformation processes occurring in the structure of the studied materials under tension, which include straightening the strands of the reinforcing layer and increasing the value of mutual pressure of the cross-lying reinforcement strands. At the beginning of the deformation process, the straightening of the threads and the increase in mutual pressure of the threads are most intense. With the increase in the level of load, the change of these parameters slows down. For example, the bending of the reinforcement strands goes into the Central tension, and the value of the load from the mutual pressure is no longer increased (tends to constant). To simulate the described processes, the basic geometrical and mechanical parameters of the material affecting the process of forming are introduced, the necessary terminology and description of the characteristics are given. Due to the high geometric nonlinearity of the all processes described in the increments, as in the initial load values there is a significant deformation of the reinforcing layer. For the quantitative and qualitative description of mechanical deformation processes occurring in the reinforcing layer, analytical dependences are derived to determine the increment of the angle of straightening of reinforcement filaments and the load caused by the mutual pressure of the cross-lying filaments at each step of the load increment. For testing of obtained dependencies shows an example of their application for flexible woven composites brands VP4126, VP6131 and VP6545. The simulation results confirmed the assumptions about the processes of straightening the threads and slowing the increase in mutual pressure of the threads. The results and dependences presented in this paper are directly related to the further generalization of the previously proposed analytical models for biaxial tension, since stretching in two directions will significantly reduce the straightening of the threads and increase the amount of mutual pressure under similar loads.

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

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

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

В работе предлагается метод определения состава входных переменных нейросетевых моделей для стационарных и нестационарных временных рядов, базирующийся на построении и анализе автокорреляционных функций. На основе предложенного метода разработаны алгоритм и программа в среде программирования Python, определяющая состав входных переменных прогнозной нейросетевой модели — персептрона, а также строящая саму модель. Осуществлена экспериментальная апробация предложенного метода на примере построения прогнозной нейросетевой модели временного ряда, отражающего потребление электроэнергии в разных регионах США, открыто опубликованной компанией PJM Interconnection LLC (PJM) — региональной сетевой организацией в Соединенных Штатах. Данный временной ряд является нестационарным и характеризуется наличием как тренда, так и сезонности. Прогнозирование очередных значений временного ряда на ос- нове предыдущих значений и построенной нейросетевой модели показало высокую точность аппроксимации, что доказывает эффективность предлагаемого метода.

This article studies a method of constructing a predictive neural network model of a time series based on determining the composition of input variables, constructing a training sample and training itself using the back propagation method. Traditional methods of constructing predictive models of the time series are: the autoregressive model, the moving average model or the autoregressive model — the moving average allows us to approximate the time series by a linear dependence of the current value of the output variable on a number of its previous values. Such a limitation as linearity of dependence leads to significant errors in forecasting.

Mining Technologies using neural network modeling make it possible to approximate the time series by a nonlinear dependence. Moreover, the process of constructing of a neural network model (determining the composition of input variables, the number of layers and the number of neurons in the layers, choosing the activation functions of neurons, determining the optimal values of the neuron link weights) allows us to obtain a predictive model in the form of an analytical nonlinear dependence.

The determination of the composition of input variables of neural network models is one of the key points in the construction of neural network models in various application areas that affect its adequacy. The composition of the input variables is traditionally selected from some physical considerations or by the selection method. In this work it is proposed to use the behavior of the autocorrelation and private autocorrelation functions for the task of determining the composition of the input variables of the predictive neural network model of the time series.

In this work is proposed a method for determining the composition of input variables of neural network models for stationary and non-stationary time series, based on the construction and analysis of autocorrelation functions. Based on the proposed method in the Python programming environment are developed an algorithm and a program, determining the composition of the input variables of the predictive neural network model — the perceptron, as well as building the model itself. The proposed method was experimentally tested using the example of constructing a predictive neural network model of a time series that reflects energy consumption in different regions of the United States, openly published by PJM Interconnection LLC (PJM) — a regional network organization in the United States. This time series is non-stationary and is characterized by the presence of both a trend and seasonality. Prediction of the next values of the time series based on previous values and the constructed neural network model showed high approximation accuracy, which proves the effectiveness of the proposed method.

Деформирование сплошных сред из материалов с памятью формы под влиянием возрастающей нагрузки и при постоянной температуре протекает обычным для металлов идеальным упругопластическим образом. При этом величина максимальных упругих деформаций много меньше предельных пластических. Восстановление формы происходит при повышенной температуре и невысоком уровне напряжений. Феноменологически «обратное» деформирование аналогично с точностью до знака изменению формыпри активном загружении силами. Так как в неупругом процессе решающую роль играет пластическая деформация, то анализ механического поведения целесообразно провести в рамках идеальной жесткопластической модели с двумя поверхностями нагружения. В этой модели поверхностям нагружения отвечают два физических состояния материала: пластическое течение при высоких напряжениях и плавление при сравнительно невысокой температуре. Во втором параграфе формулируется задача деформирования жесткопластических сред при постоянной температуре в двух формах: в виде принципа виртуальных скоростей с условием текучести Мизеса и как требование минимальности диссипативного функционала. Доказываются равносильность принятых формулировок и существование обобщенных решений в обоих принципах. В третьем параграфе изучается жесткопластическая модель сплошной среды при изменяющейся температуре с двумя поверхностями нагружения. Для принятой модели формулируются два оптимальных принципа, связывающих внешние нагрузки и скорости перемещений точек среды как при активном нагружении, так и в процессе восстановления формыпр и нагревании. Доказано существование обобщенных скоростей для широкого класса трехмерных областей. Связь вариационных принципов и изменяющейся температуры обеспечивается включением в расчетную схему первого и второго начал термодинамики. Существенно, что в процессе доказательств используется только феноменологическое описание явления. Аустенитно-мартенситные превращения сплавов, которые часто являются основными при объяснении механического поведения материалов с памятью формы, не используются. В четвертом параграфе дано определение материалов с памятью формы как сплошных сред с двумя поверхностями нагружения, доказано существование решений в принятых ограничениях. Показана адекватность модели и опытов по деформированию материалов с памятью формы. В заключении формулируются математические задачи, которые представляются интересными в будущих исследованиях.

Under increasing loading and at a constant temperature shape memory solids become deformed in an ideal elastic plastic way as other metals, and the maximum elastic strains are much less than the ultimate plastic ones. The shape is restored at the elevated temperature and low stress level. Phenomenologically, the «reverse» deformation is equivalent to the change in shape under active loading up to sign. Plastic deformation plays a leading role in a non-elastic process; thus, the mechanical behavior should be analyzed within the ideal rigid-plastic model with two loading surfaces. In this model two physical states of the material correspond to the loading surfaces: plastic flow under high stresses and melting at a relatively low temperature. The second section poses a problem of deformation of rigid-plastic bodies at the constant temperature in two forms: as a principle of virtual velocities with the von Mises yield condition and as a requirement of the minimum dissipative functionаl. The equivalence of the accepted definitions and the existence of the generalized solutions is proved for both principles. The third section studies the rigid-plastic model of the solid at the variable temperature with two loading surfaces. For the assumed model two optimal principles are defined that link the external loads and the displacement velocities of the solid points both under active loading and in the process of shape restoration under heating. The existence of generalized velocities is proved for the wide variety of 3D domains. The connection between the variational principles and the variable temperature is ensured by inclusion of the first and second principles of thermodynamics in the calculation model. It is essential that only the phenomenological description of the phenomenon is used in the proving process. The austenite-tomartensite transformations of alloys, which are often the key elements in explanations of the mechanical behavior of shape memory materials, are not used here. The fourth section includes the definition of the shape memory materials as solids with two loading surfaces and proves the existence of solutions within the accepted restrictions. The adequacy of the model and the experiments on deformation of shape memory materials is demonstrated. In the conclusion mathematical problems that could be interesting for future research are defined.

Одной из особенностей беспилотных автомобильных транспортных средств является их способность к организованному движению в форме кластеров: последовательности движущихся с единой скоростью транспортных средств. Влияние образования и движения этих кластеров на динамику транспортных потоков представляет большой интерес. В настоящей работе предложена качественная имитационная модель кластерного движения беспилотных транспортных средств в гетерогенной транспортной системе, состоящей из двух типов агентов (транспортных средств): управляемых человеком и беспилотных. В основу описания временной эволюции системы положены правила 184 и 240 для элементарных клеточных автоматов. Управляемые человеком транспортные средства перемещаются по правилу 184 с добавлением случайного торможения, вероятность которого зависит от расстояния до находящегося впереди транспортного средства. Для беспилотных транспортных средств используется комбинация правил, в том числе в зависимости от типа ближайших соседей, в некоторых случаях независимо от расстояния до них, что привносит в модель нелокальное взаимодействие. При этом учтено, что группа последовательно движущихся беспилотных транспортных средств может сформировать организованный кластер. Исследовано влияние соотношения типов транспортных средств в системе на характеристики транспортного потока при свободномд вижении на круговой однополосной и двухполосной дорогах, а также при наличии светофора. Результаты моделирования показали, что эффект образования кластеров имеет существенное влияние при свободномдвиж ении, а наличие светофора снижает положительный эффект приблизительно вдвое. Также исследовано движение кластеров из беспилотных автомобилей на двухполосных дорогах с возможностью перестроения. Показано, что учет при перестроении беспилотными транспортными средствами типов соседних транспортных средств (беспилотное или управляемое человеком) положительно влияет на характеристики транспортного потока.

The gradual incorporation of automated vehicles into the global transport networks leads to the need to develop tools to assess the impact of this process on various aspects of traffic. This implies a more organized movement of automated vehicles which can form uniformly moving platoons. The influence of the formation and movement of these platoons on the dynamics of traffic flow is of great interest. The currently most developed traffic flow models are based on the cellular automaton approach. They are mainly developed in the direction of increasing accuracy. This inevitably leads to the complication of models, which in their modern form have significantly moved away from the original philosophy of cellular automata, which implies simplicity and schematicity of models at the level of evolution rules, leading, however, to a complex organized behavior of the system. In the present paper, a simulation model of connected automated vehicles platoon dynamics in a heterogeneous transport system is proposed, consisting of two types of agents (vehicles): human-driven and automated. The description of the temporal evolution of the system is based on modified rules 184 and 240 for elementary cellular automata. Human-driven vehicles move according to rule 184 with the addition of accidental braking, the probability of which depends on the distance to the vehicle in front. For automated vehicles, a combination of rules is used depending on the type of nearest neighbors, regardless of the distance to them, which brings non-local interaction to the model. At the same time, it is considered that a group of sequentially moving connected automated vehicles can form an organized platoon. The influence of the ratio of types of vehicles in the system on the characteristics of the traffic flow during free movement on a circular one-lane and two-lane roads, as well as in the presence of a traffic light, is studied. The simulation results show that the effect of platoon formation is significant for a freeway traffic flow; the presence of a traffic light reduces the positive effect by about half. The movement of platoons of connected automated vehicles on two-lane roads with the possibility of lane changing was also studied. It is shown that considering the types of neighboring vehicles (automated or human-driven) when changing lanes for automated vehicles has a positive effect on the characteristics of the traffic flow.

Волны возбуждения являются прообразом самоорганизующихся динамических структур в неравновесных системах. Они характеризуются своей собственной внутренней динамикой, приводящей к формированию бегущих волн различных типов и форм. Яркие примеры — это вращающиеся спирали и скрученные свитки. Интересная и сложная задача — найти способы управления их поведением, применяя внешние сигналы, влияющие на распространяющиеся волны. В качестве такого воздействия мы используем внешние электрические поля, наложенные на возбудимую реакцию Белоусова–Жаботинского (БЖ). Существенные эффекты влияния полей на волны включают изменение скорости волны, обращение направления распространения, взаимное уничтожение вращающихся в противоположных направлениях спиральных волн и переориентацию нитей скрученных свитков. Эти эффекты могут быть объяснены в численных экспериментах, при этом существенную роль играет отрицательно заряженный ингибиторбромид. Эффекты электрического поля также были исследованы в биологических возбудимых средах, таких как социальные амебы Dictyostelium discoideum. Совсем недавно мы начали исследовать влияние электрического поля на реакцию БЖ, протекающую в водно-масляной микроэмульсии. Удалось наблюдать дрейф сложных структур, а также изменение вязкости и электрической проводимости. Мы обсуждаем предположение, что эта система может выступать в качестве модели для дальнодействующего взаимодействия между нейронами.

Excitation waves are a prototype of self-organized dynamic patterns in non-equilibrium systems. They develop their own intrinsic dynamics resulting in travelling waves of various forms and shapes. Prominent examples are rotating spirals and scroll waves. It is an interesting and challenging task to find ways to control their behavior by applying external signals, upon which these propagating waves react. We apply external electric fields to such waves in the excitable Belousov–Zhabotinsky (BZ) reaction. Remarkable effects include the change of wave speed, reversal of propagation direction, annihilation of counter-rotating spiral waves and reorientation of scroll wave filaments. These effects can be explained in numerical simulations, where the negatively charged inhibitor bromide plays an essential role. Electric field effects have also been investigated in biological excitable media such as the social amoebae Dictyostelium discoideum. Quite recently we have started to investigate electric field effect in the BZ reaction dissolved in an Aerosol OT water-in-oil microemulsion. A drift of complex patterns can be observed, and also the viscosity and electric conductivity change. We discuss the assumption that this system can act as a model for long range communication between neurons.

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

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

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

The article contains the description of developed mathematical models of fractures which can be used for numerical solution of exploration seismology problems with use of grid-characteristic method on unstructured triangular and tetrahedral meshes. The base of developed models is the concept of infinitely thin fracture. This fracture is represented by contact boundary. Such approach significantly reduces the consumption of computer resources by the absence of the mesh definition inside of fracture necessity. By the other side it lets state the fracture discretely in integration domain, therefore one can observe qualitative new effects which are not available to observe by use of effective models of fractures, actively used in computational seismic.

The main target in the development of models have been getting the most accurate result. Developed models thet can receive the response close to the actual response of the existing fracture in geological environment. We considered fluid-filled fractures, glued and partially glued fractures, and also fractures with dynamical friction force. Fracture behavior determinated by the nature of condition on the border.

Empty fracture was represented as free boundary condition. This condition give us opportunity for total reflection of wave fronts from fracture. Fluid-filling provided the condition for sliding on the border. Under this condition, there was a passage of longitudinal and total reflection of converted waves. For the real fractures, which has unequal distance between the borders has been proposed the model of partially glued fracture. At different points of the fracture's boundary were sat different conditions. Almost the same effect is achieved by using a fracture model of dynamic friction condition. But its disadvantage is the inabillity to specify the proportion of fracture's glued area due to the friction factor can take values from zero to infinity. The model of partially glued fracture is devoid of this disadvantage.

В данной работе исследованы статистические взаимосвязи между размерными и продукционными характеристиками фито- и зоопланктона, обитающего в водах Вислинского и Куршского заливов Балтийского моря. Исследования фито- и зоопланктона в пределах российских частей акваторий Вислинского и Куршского заливов проводили ежемесячно (с апреля по ноябрь) в рамках программы многолетнего мониторинга состояния экосистем заливов. Размерная структура планктонных организмов — основа понимания развития продукционных процессов, механизмов формирования видового разнообразия планктона и функционирования экосистем заливов. По результатам работы установлено, что максимальная скорость фотосинтеза и величина интегральной первичной продукции меняются по степенному закону с изменением среднего ценотического объема клеток фитопланктона. Полученный результат показывает, что чем меньше размер клеток водорослей в фитопланктонных сообществах, тем активнее в них протекают процессы метаболизма и тем эффективнее усваивается солнечная энергия. Показано, что формирование видового разнообразия планктона в экосистемах заливов самым тесным образом связано и с размерной структурой планктонных сообществ, и с особенностями развития продукционных процессов. Предложена структура пространственно однородной математической модели планктонной трофической цепи для экосистем заливов, учитывающая размерные спектры и характеристики фито- и зоопланктона. Параметры модели — размерно-зависимые показатели, аллометрически связанные со средними объемами клеток и организмов в разных диапазонах их размеров. В модели предложен алгоритм изменения во времени коэффициентов предпочтения в питании зоопланктонных организмов. Разработанная размерно-зависимая математическая модель водных экосистем позволяет учесть воздействие турбулентного обмена на размерную структуру и временную динамику планктонной пищевой цепи Вислинского и Куршского заливов. Модель может быть использована для исследования различных режимов динамического поведения планктонной системы в зависимости от изменений значений ее параметров и внешних воздействий, а также для количественной оценки перераспределения потоков вещества в экосистемах заливов.

In the paper the statistical relationships between the size and production characteristics of phytoplankton and zooplankton of the Vistula and Curonian lagoons, the Baltic Sea, were investigated. Research phytoplankton and zooplankton within the Russian part of the area of the Vistula and the Curonian lagoon was carried out on the monthly basis (from April to November) within the framework of long-term monitoring program on evaluating of ecological status of the lagoons. The size structure of plankton is the basis for understanding of the development of production processes, mechanisms of formation of the plankton species diversity and functioning of the lagoon ecosystems. As results of the work it was found that the maximum rate of photosynthesis and the integral value of the primary production with a change in cell volume of phytoplankton are changed according to a power law. The result shows that the smaller the size of algal cells in phytoplankton communities the more actively occur metabolism and the more effective they assimilate the solar energy. It is shown that the formation of plankton species diversity in ecosystems of lagoons is closely linked with the size structure of plankton communities and with features of development of the production processes. It is proposed the structure of a spatially homogenous mathematical model of the plankton food chain for the lagoon ecosystems taking into account the size spectrum and the characteristics of phytoplankton and zooplankton. The model parameters are the sizedependent indicators allometrically linked with average volumes of cells and organisms in different ranges of their sizes. In the model the algorithm for changes over time the coefficients of food preferences in the diet of zooplankton was proposed. Developed the size-dependent mathematical model of aquatic ecosystems allows to consider the impact of turbulent exchange on the size structure and temporal dynamics of the plankton food chain of the Vistula and Curonian lagoons. The model can be used to study the different regimes of dynamic behavior of plankton systems depending on the changes in the values of its parameters and external influences, as well as to quantify the redistribution of matter flows in ecosystems of the lagoons.

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

We present results of a study aimed at identification of a controlled object’s channels based on postprocessing of measurements with development of a model of a multiple-input controlled object and subsequent active modelling experiment. The controlled object model is developed using approximation of its behavior by a neural network model using trends obtained during a passive experiment in the mode of normal operation. Recurrent neural network containing feedback elements allows to simulate behavior of dynamic objects; input and feedback time delays allow to simulate behavior of inertial objects with pure delay. The model was taught using examples of the object’s operation with a control system and is presented by a dynamic neural network and a model of a regulator with a known regulation function. The neural network model simulates the system’s behavior and is used to conduct active computing experiments. Neural network model allows to obtain the controlled object’s response to an exploratory stimulus, including a periodic one. The obtained complex frequency response is used to evaluate parameters of the object’s transfer system using the least squares method. We present an example of identification of a channel of the simulated control system. The simulated object has two input ports and one output port and varying transport delays in transfer channels. One of the input ports serves as a controlling stimulus, the second is a controlled perturbation. The controlled output value changes as a result of control stimulus produced by the regulator operating according to the proportional-integral regulation law based on deviation of the controlled value from the task. The obtained parameters of the object’s channels’ transfer functions are close to the parameters of the input simulated object. The obtained normalized error of the reaction for a single step-wise stimulus of the control system model developed based on identification of the simulated control system doesn’t exceed 0.08. The considered objects pertain to the class of technological processes with continuous production. Such objects are characteristic of chemical, metallurgic, mine-mill, pulp and paper, and other industries.

В различных областях науки, техники, защиты окружающей среды, в строительстве актуальными являются вопросы изучения процессов взаимодействия пористых материалов с веществами, находящимися в различных агрегатных состояниях. Особенно актуальными с точки зрения экологии и защиты окружающей среды являются исследования процессов взаимодействия пористых материалов с водой в жидкой и газообразной фазе. Поскольку в одном моле воды содержится 6.022140857 · 10²³ молекул H₂O, для описания свойств, например, водяного пара в поре в основном используются макроскопические подходы, рассматривающие водяной пар как сплошную среду в рамках аэродинамики. В данной работе построена и использовалась для моделирования макроскопическая двумерная диффузионная модель [Bitsadze, Kalinichenko, 1980] поведения водяного пара внутри изолированной поры. Наряду с макроскопической моделью в работе предложена микроскопическая модель поведения водяного пара внутри изолированной поры, построенная в рамках молекулярно-динамического подхода [Gould et al., 2005]. В данной модели на основе классической механики Ньютона описывается движение каждой молекулы воды, взаимодействующей как с другими молекулами воды, так и со стенками поры. Рассматривается эволюция системы «водяной пар – пора» с течением времени. В зависимости от внешних по отношению к поре условий система эволюционирует к различным состояниям равновесия, которые характеризуются различными значениями макроскопических характеристик, таких как температура, плотность, давление. Сравнение результатов молекулярно-динамического моделирования с результатами вычислений на основе макроскопической диффузионной модели и экспериментальными данными позволяет сделать вывод о необходимости сочетания макроскопического и микроскопического подхода для адекватного и более точного описания процессов взаимодействия водяного пара с пористыми материалами.

In various areas of science, technology, environment protection, construction, it is very important to study processes of porous materials interaction with different substances in different aggregation states. From the point of view of ecology and environmental protection it is particularly actual to investigate processes of porous materials interaction with water in liquid and gaseous phases. Since one mole of water contains 6.022140857 · 10²³ molecules of H₂O, macroscopic approaches considering the water vapor as continuum media in the framework of classical aerodynamics are mainly used to describe properties, for example properties of water vapor in the pore. In this paper we construct and use for simulation the macroscopic two-dimensional diffusion model [Bitsadze, Kalinichenko, 1980] describing the behavior of water vapor inside the isolated pore. Together with the macroscopic model it is proposed microscopic model of the behavior of water vapor inside the isolated pores. This microscopic model is built within the molecular dynamics approach [Gould et al., 2005]. In the microscopic model a description of each water molecule motion is based on Newton classical mechanics considering interactions with other molecules and pore walls. Time evolution of “water vapor – pore” system is explored. Depending on the external to the pore conditions the system evolves to various states of equilibrium, characterized by different values of the macroscopic characteristics such as temperature, density, pressure. Comparisons of results of molecular dynamic simulations with the results of calculations based on the macroscopic diffusion model and experimental data allow to conclude that the combination of macroscopic and microscopic approach could produce more adequate and more accurate description of processes of water vapor interaction with porous materials.

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

The paper presents the physical and mathematical formulation of the problems of internal ballistics of an artillery shot for a charge consisting of a set of powder tubes and their stress-strain state. Combustion and movement of a bundle of powder tubes along the barrel channel is modeled by an equivalent tubular charge of all-round combustion. It is assumed that the equivalent tube moves along the axis of the bore. The speed of movement of an equivalent tubular charge and its current position are determined from Newton’s second law. When calculating the flow parameters, two-dimensional axisymmetric equations of gas dynamics were used, for the solution of which an axisymmetric orthogonalized difference grid is constructed, which adapts to the flow conditions. The control volume method is used to numerically solve the system of gas-dynamic equations. The gas parameters at the boundaries of the control volumes are determined using a self-similar solution to the Godunov’s problem of the decay of an arbitrary discontinuity. The stress-strain state is modeled for a separate burning powder tube located in the field of gas-dynamic parameters. The calculation of the gas-dynamic parameters of the shot is carried out without taking into account the deformed state of the powder elements. The behavior of powder elements during firing is considered under these conditions. The finite element method with the division of the calculation area into triangular elements is used to solve the problem of elasticity. In the process of powder tube burnout, the computational grid on each time layer of the dynamic problem is completely updated due to a change in the boundaries of the powder element due to combustion. The paper shows the time dependences of the parameters of the internal ballistics process and the stress-strain state of powder elements, as well as the distribution of the main parameters of the flow of combustion products at different points in time. It has been established that the tubular powder elements during the shot experience significant deformations, which must be taken into account when solving the basic problem of internal ballistics. The data obtained give an idea of the level of equivalent stresses acting at various points of the powder element. The results obtained indicate the relevance of the conjugate formulation of the problem of gas dynamics and the stress-strain state for charges consisting of tubular powders, since this allows a new approach to the design of tubular charges and opens up the possibility of determining the parameters on which the physics of the combustion process of gunpowder significantly depends, therefore, and the dynamics of the shot process.