molpit
Login:
Password:
remember

После анализа моего автореферата PIT и ZIM пришли к выводу, что я должен привести все рассчитанные диаграммы с указанием размерностей энергии, времени и начальных системы. Поэтому пришлось вспомнить и изучить новые приемы обезразмеривания и восстановления размерностей DIA00095, а также прикинуть реальные коэффициенты упругости DIA00096.

Итак встречайте:

Пример расчета

Для примера промоделируем образование связи между двумя объектами атомных масштабов. Возьмем два атома углерода, и из условия, что при сближении появляется объект организации, который притягивает атомы.

Динамика многоуровневой системы описывается следующим гамильтонианом DIA00069:

 \left\lbrace \begin{aligned}<br/>
	H = & \sum_{i,j} \mathcal{T}_j^i \left( \frac{(p^i_j)^2}{2m^i_j}<br/>
	+ \sum_{k \neq i} \mathcal{T}_j^k \frac 1 2 U_j^L (q^i_j, q^k_j) +\right.\\<br/>
	& + \sum_{k \not\in D^i_j} \mathcal{T}_{j+1}^k  \frac 1 2 U_j^{L+1} (q^i_j, q^k_{j+1}) 	<br/>
	+ \sum_{k \not\in P^i_j} \mathcal{T}_{j-1}^k	\frac 1 2 U_j^{L-1} (q^i_j, q^k_{j-1}) + \\<br/>
	& \left. + \sum_{k \in P^i_j} \mathcal{T}_{j-1}^k  \frac 1 2 U_j^P (q^i_j, q^k_{j-1})<br/>
	+ \sum_{k \in D^i_j} \mathcal{T}_{j+1}^k \frac 1 2 U_j^D (q^i_j, q^k_{j+1}) \right),<br/>
\end{aligned}\right.

где j — номер уровня организации; i — номер объекта организации на уровне; q_j^i — координаты i-ого объекта организации на уровне j; p_j^i — импульс i-го объекта на уровне j; m_j^i — масса i-го объекта на уровне j; U_j^L — потенциал взаимодействия между объектами уровня j; U_j^{L+1} — объектами уровней j и j+1; U_j^{L-1} — объектами уровней j и j-1; U_j^P — объектов уровня j c родительскими объектами уровня j-1; U_j^D — объектов уровня j c дочерними объектами уровня j+1; P^i_j — множество родительских объектов для i-го объекта уровня j; D^i_j — для i-го объекта уровня j.

Выберем потенциалы, триггер и сукцессор соответствующие постановке задачи:

 \left\lbrace<br/>
\begin{aligned}<br/>
U_j^L &= 1000 \left(\left( (r_i^j + r_i^k)/(q_i^j - q_i^k)\right)^{12}-2\left((r_i^j + r_i^k)/(q_i^j - q_i^k) \right)^{6}\right),\\<br/>
U_j^P &= U_j^D = 500 \, (q_i^j - q_i^k) / 2, \qquad U_j^{L-1}  = U_j^{L+1} = 0,\\<br/>
\mathcal{T}_j^i & = \Theta\left(\sum_{m,n \neq m, \lbrace n, m \rbrace \nsubseteq D^n_{j-1}, \lbrace n, m \rbrace \nsubseteq D^m_{j-1} } \sqcap\left( \frac{q^n_{j-1} - q^m_{j-1}}{r^n_{j-1} + r^m_{j-1}}\right) \right.\\<br/>
& - \left. \sum_{m \in P^n_{j},n \in P^n_{j} } \sqcap\left( \frac{q^n_{j-1} - q^m_{j-1}}{1.1(r^n_{j-1} + r^m_{j-1})}\right) \right),\\<br/>
\mathcal{S}_j & = (\mathcal{O}_j^n + \mathcal{O}_j^m) \left( \begin{array}{cc} 1/2 & 0 \\ 0 & 1/2  \end{array} \right), \quad \mathcal{O}_j^i = \left(\begin{array}{cc} m_j^i & r_j^i \end{array}\right).<br/>
\end{aligned}\right.

В результате расчета в приложении Levels 1.5 получим результат, изображенный на рисунке 1. Слева без диссипации энергии, справа с учетом диссипации 0.1 Па с.

double (70Кб)
Рис. 1. Расчет согласно уравнению выше, с массой объектов 12 а.е.м., радиусом — 1 ангстрем, начальным расстоянием между центрами 3 ангстрема, коэффициент жесткости для возникающего объекта 500 кг/c²

Ivan Denisov 08 Dec 2013 14:52
© International Open Laboratory for Advanced Science and Technology — MOLPIT, 2009–2021