64713

Physics-Mathematical modeling

Лекция

Экономическая теория и математическое моделирование

Using the general laws of the thermodynamics such integrals can be written for any function we are interested in. For example, the equilibrium states can be described by the minimum of the potential energy...

Английский

2014-07-10

401 KB

0 чел.

Physics-Mathematical modeling is one of the basic research methods of approach to various phenomenon and processes. A co-ordinate space is associated to an every natural object here and the processes are described by differential equations. Commonly these equations couple a number of physical functions Fк(xi) of the independent variable xi, defined in. As example to an element of a cooler (fig.1а), we can associate co-ordinate space (xi, х=1, 2, 3) (fig.1в). In this space we can define functions of the stress ij(xi), of the temperature T(xi) and of the velocity vi(xi). In the courses on “Mechanics of materials” “Mechanics of fluids” and “Heat and mass transfer” is presented, that differential equations, connecting these functions, can be derived from equilibrium of the infinitesimal area around point in .

   

а)      b)

Фиг.1

Фиг.2

If we examine mechanical equilibrium of the volume with infinitesimal edges size dx1, dx2, dx3 (fig.2), we receive:

(1.1)  ,

          ,

    ,

were Хi  are the components of the dead load,  is the density and ui are the displacements (the equations will be discussed in the next parts of these materials).

For this volume the next boundary conditions are in force:

33=0 at x3=R  - there isn’t any stress on the open  surface,

 33=p (p - pressure) for x3 = r  -  the stress on the inner side of the pipe is equal to the pressure.

(1.1а)  ij=oij at x1=0  -  the stress on the left boundary of the investigated area is equal to the stress on the right side of the next area or is defined by fixing.

 ij=lij at x1=l  - the stress on the right boundary of the investigated area is equal to the stress on the left side of the next area or is defined by fixing.

In the same way we can derive from the thermal quantity balance the next:

(1.2)   ,

were Т is the temperature and а is the conduction coefficient. The boundary conditions are:

  Т=Т0 at х1=0 for the temperature,

(1.2а)    at х1→0 for thermal flux.

For the fluid stream the equation of the continuity is in force:      

(1.3)   .

In this case the boundary conditions are:

(1.3а)  V1 = 0  за x3 = r; V1 = V1(o) за x1 = 0; V1= V1(l) за x1 = l, V2 = 0,  V3 = 0.

So, in existence of sufficiency of boundary conditions 1а-3а, the equations 1-3 can be solved and the functions Fк(xi), which describe the object state, can be determined. But the solution of these equations (by finite differences as example), runs into big difficulties when complex geometries or complex boundary conditions exist. For solving the problem nowadays it’s increasingly used the minimum-principles (variation principles) of Hamilton and Lagrange. According to these principles some definite integrals, representative for the process under investigation, obtain the minimum value if the unknown integral function describes the process in the right way.   

Using the general laws of the thermodynamics such integrals can be written for any function we are interested in. For example, the equilibrium states can be described by the minimum of the potential energy, the motion of mechanical systems - by the minimum of the kinetic potential etc. In this way all equations we are interested in can be obtained from one origin – variation principals of the continuum.

Thus the solution of the problem can be obtained by a procedure in which an approximate continuous function is assumed to represent the solution and an above integral formulation is used to create a system of algebraic equations. Finite element method assumes that the investigated area is divided into smaller pieces (elements) and the procedure above mentioned is implemented for every element separately as a first step of the solution (discretization), the results are collated in a global matrix at the second step (assembling) and at the end – a system of algebraic equations for the unknown quantities in the all area is obtained (expanding). A brief description of this three-step procedure is presented in the next part of these materials.


 

А также другие работы, которые могут Вас заинтересовать

41964. Написать программу на языке C++, моделирующую поведение курицы (Hen) путём создания соответствующего класса 14.17 KB
  Листинг программы: include iostrem include cstring include cmth include cstdlib using nmespce std; clss Chickhen { privte: chr nme; double w h f; Кормление урожай норма кормления sttic int e; норма яйценосности public: Chickhenvoid; Chickhenchr double; Chickhenconst Chickhen ; virtul Chickhen; double hrvest; double feeddouble; }; int Chickhen::e=10; Chickhen::Chickhen { w=0; h=0; f=0; nme=new chr[7]; strcpy nme nonme ; } Chickhen::Chickhen chrndouble F { nme=new chr[strlenn1]; strcpynmen; f=F; h=0; w=0;...
41968. Дослідження стійкості ланки другого порядку 114.05 KB
  Для лінійних систем автоматичного керування, які описуються характеристичним рівнянням виду a0pn+a1pn-1+…+an-1p+an=0 стійкість не залежить від величини і вигляду збурення і визначається коренями характеристичного рівняння, яке залежить від параметрів системи Для зручності зафіксуємо L C та змінюватимемо R withinttrns; urovnenie:=TTpp2xiTp1; h:=k p urovnenie; l:=invlplcehpt; sol:=solveurovneniep: sol[1];sol[2]; Аперіодичний процес Вибираємо L=50мГн.05;C:=2010^6;R:=250;T:=sqrtLC;xi:=RsqrtC L 2;k:=1;p1:=sol[1];p2:=sol[2];задання параметрів для даного виду процесу l:=invlplcehpt;розрахунок зворотнього перетворення Лапласа plotlt=0.05;C:=2010^6;R:=100;T:=sqrtLC;xi:=RsqrtC L 2;k:=1;p1:=sol[1];p2:=sol[2]; l:=invlplcehpt:...
41969. ДОСЛIДЖЕННЯ ВЕКТОРНИХ ПЛОТТЕРІВ 78.68 KB
  Все рассматриваемые здесь команды находятся в основной части языка HPGL 2. Первыми идут команды ини рйализации для установки размера изображения и другие параметры после них следуют команды для прорисовки линий фигур и трок символов а также одна или две команды для завершения процесса. Некоторые команды имеющие числовые аргументы требуют целых значений в то время как другие команды допускают наличие чисел с десятичной точкой. Некоторые команды передают результаты обратно хосткомпьютеру: например 01 сообщает идентификацию модели...
41970. Функції введення/виведення printf(), scanf(). Лінійні обчислювальні процеси 14.14 KB
  Обладнання: ПКПО Borlnd C Хід роботи Вивчити теоретичні відомості Ознайомитися з форматом функцій printf і scnf.h void min { long ; double b; unsigned c; flot d; cout Вводите n ; cin b c d; cout Long n ; printf 16.5d ; printf n double n ; printf 16.
41971. Обчислювальний процес, що розгалужується, з різними логічними умовами: оператор if... else, умовна операція (?:), оператор switch, оператор break, оператор goto 23.72 KB
  else умовна операція : оператор switch оператор brek оператор goto Ціль роботи: Вивчити реалізацію в мові ветвящихся обчислювальних процесів . Навчитися писати програми використовуючи оператори: розгалуження if.else переключення switch у парі з оператором brek оператор переходу goto тернарную умовну операцію .
41972. Розробка програм з циклічними обчислювальними процесами 44.87 KB
  Розробка програм з циклічними обчислювальними процесами Ціль роботи: Вивчити написання програм мовою С, використовуючи ітераційні (циклічні) методи, освоїти основні оператори, що підтримують роботу з циклами (for, while, do... while). Навчитися писати програми, використовуючи дані оператори.