Есть ли способ не устранить результаты «NAN» для OpenSees NodereAction (i); где я = узловой? ⇐ Python

[Anonymous]

imports < /strong>
importseSpy.openseees as ops
import opsvis как opsv
import numpy как np
import matplotlib.pyplot как plt < /p>
Проблема < /strong>
i -nefiving [nan, nan, nan, nan, nan, nan, nan lonse, nan, nan, nan, lonseless, nan, nan, nan, nan, nan, nan, nan, nan, nan, lon Анализ трехмерной модели, когда я называю ops.nodereaction (i) , где i = Nodetag, даже когда ops.Analyze (1) успешно с 0, а модель была применена с нагрузкой. E0, b) и ops.uniaxialmaterial ('concrete07', 2, fcu, epsc, ec, ft, et, xp, xn, r) со всеми назначенными аргументами. 3, '-ndf', 6)
Срезы волокон использовались как для ops.dispbeamcolumn вертикальные элементы, так и для ops.nonlinearbeamcolumn горизонтальные элементы, использующие ops.section ('fiber', sectag, '-gri-emements, следовать ops.section (' Fiber ', decectag, eme-gri, следовал ops.section (' Fiber ', dececag, следовал ops.section. *ops.patch ('Quad', mattag, nfibcorezp, nfibcoreyp, verticescore) , где Sectag составляет 1 для вертикального элемента, и 2, 3, 4 для 3 горизонтальных элементов (они отражают о вертикальном элементе, когда они образуют форму T), другие аргументы - это то, что они называют их, и Mattag = 2 для 3 -х для 3 -х похват; За последующим *Ops.Layer ('Straight', 1, NBARZP, ABAR1, *start_botlayr, end_botlayr) для стальной арматуры. Все аргументы были объявлены соответствующим образом до вызова команды. Геометрические трансформации были определены с помощью *ops.geomtransf ('pdelta', transftag, transfvect) соответственно для всех разделов с трансфтаг 1, 2, 3 и 4 для вертикального элемента и трех горизонтальных элементов соответственно. Transfvect были объявлены как [0, -1, 0] для вертикального элемента, и [-1, 0, 0] для горизонтальных элементов перед вызовом команды *Ops.geomtransf ('pdelta', Transftag, Transfvect) .
Модель: модель -модели и 3 -й интеграции
Модель. Вертикальный консольный элемент с использованием элемента ops.dispbeamcolumn с одним нижним узлом (nodetag = 1), приложенным с ops.fix (1, 1,1,1,1,1,1,1) и еще один верхний узел (Nodetag = 2), соединенные с ops.nononlinearbeamcolumn emezontal samesementsemements (inodetag = 2). [3,4], [4,5], [5,6]) через *ops.rigidlink ('beam', [2,3]) и к ops.nonlinearbeamcolumn горизонтальные элементы луча (соединение Nodetags [7,8], [8,9], [9,10] через A *ops. [2,7]) с другой стороны, сделав T -образную консольную модель. Три типа элементов горизонтального луча имеют разные поперечные сечения.
ops.beamintegration ('lobatto', 1, sectag, n) использовали для вертикального ops.dispbeamcolumn с sectag = 1 и n = 4. For ops.nonlinearBeamColumn horizontal elements, numIntgrPts=4 was as well used in *ops.element('nonlinearBeamColumn', eleTag, ops.eleNodes, numIntgrPts, secTag, transfTag)
Load
The load was applied on Nodetag = 2 As Ops.load (2, 100,0, 0,0, 0,0, 0,0, 0,0, 0,0) после определения Ops.Timeseries ('Constant', 1) , pattern ('plain', 1, 1) . ops.system ('bandgen') , ops.numberer ('rcm') , ops.constraints ('transformation') , ops.integrator ('loadcontrol', 1) , Ops.Integrator ('lineControl', 1) , Ops. ops.analysis('Static').
I executed analyze(1) command and printed it and it was successful as it returned zero (0) output.
However, when I tried to execute the command ops.nodeReaction(1) after calling ops.reactions() in a preceding line, and printed the Результаты я получил [NAN, NAN, NAN, NAN, NAN, NAN]. Я даже построил модель в 3D с opsv.plot_extruded_shapes_3d (ele_shapes)
со всеми «ele_shapes dict», соответственно, следовал plt.axis ('ровно') и экстрадированная форма (3d), как я ожидал. Разделы выглядели именно так, как я ожидал при построении с помощью opsv.plot_fiber_section (fib_sec_list, fillflag = 1, matcolor = matcolor) и plt.axis ('ровно') со всеми параметрами, предварительно предварительно декелянными. Есть помощь? Вывод:
0
[a, b, c, d, e, f], где a, b, c, d, e, f являются некоторыми значениями поплавки < /p>
Фактический результат вывод:
0
[nan, nan, nan, nan, nan, nan] < /p>

Подробнее здесь: https://stackoverflow.com/questions/796 ... ioni-where