Openfast: Динамическое изменение массовых свойств лопатки в рамках заданного моделирования в OpoenFAST и ElastoDyn

Мне кажется странным, что вы используете p_ED в FAST_Prog.f90. Я ожидал, что это будет Turbine(1)%ED%p . Если вы сами добавили тип p_ED , ни одна из переменных не будет инициализирована, поэтому она либо не будет проходить через циклы, либо приведет к ошибкам сегментации.

Дорогая Бонни,

Большое спасибо за быстрый ответ.
Я сам добавляю тип p_ED в заголовок FAST_prog.f90, потому что я думал, что аргумент в подпрограмме Fortran работает так же, как в функции Matlab, но, очевидно, это не так.
После того, как я использовал Turbine(1)%ED%p в FAST_Prog.f90, он заработал, прошел через петли и изменил массу лопаточного элемента.
Теперь я проверяю, изменяются ли формы режима лезвия, а также жесткость закрылка и кромки с изменением массы элемента лезвия.
Не могли бы вы оставить этот вопрос открытым? Потому что в будущем у меня наверняка появятся дополнительные вопросы по этой теме.

С наилучшими пожеланиями
Лоуренс

Уважаемый @bjonkman ,

Со времени последней публикации я внес некоторые изменения в исходный код ElstoDyn.f90, так что не только масса элемента лезвия p%BElmntMass(J,K) изменяется в ходе моделирования, но также и все параметры в SUBROUTINE Coeff(p,InputFileData, ErrStat, ErrMsg) . Для этого я создал SUBROUTINE UpdateCoeff(p, ErrStat, ErrMsg) , который в конечном итоге является копией подпрограммы Coeff () только с одним отличием - подпрограмма UpdateCoeff () не имеет типа данных ED_InputFile , потому что будет вызываться подпрограмма UpdateCoeff внутри подпрограммы ED_CalcOutput в середине моделирования, тогда как InputFileData уже уничтожается после инициализации. вызов подпрограммы UpdateCoeff работает следующим образом:

SUBROUTINE ED_CalcOutput( t, u, p, x, xd, z, OtherState, y, m, ErrStat, ErrMsg )
!..................................................................................................................................
!...............................................................................................................................
   ! set the values in the AllOuts array:
   !...............................................................................................................................
      ! Define the output channel specifying the current simulation time:

   m%AllOuts(  Time) = REAL( t, ReKi )

      ! Change Blade Element Masses:

    IF ( t == 0.50 ) THEN

        p%MassB(K,J) = 1.5*p%MassB(K,J)

        ! start with recalculation of the Parameters in subroutine Ceoff()

        CALL UpdateCoeff(p, ErrStat, ErrMsg)

    END IF
! This part is done to generate the blade mass element SpnMb as an output to see if the masses changed or not, and it did work  
   DO K = 1,p%NumBl
      DO I = 1, p%NBlGages

                BlElmMass = p%BElmntMass(I,K)

        m%AllOuts(  SpnMb(I,K)) = REAL( BlElmMass, R8Ki)

      END DO !I
   END DO !K

Но теперь, из-за отсутствия типа данных ED_Inputfile в подпрограмме UpdateCoeff, я сгенерировал новые параметры p% чтобы заменить InputFileData% в подпрограмме Coeff и использовать их после этого в подпрограмме UpdateCoeff следующим образом:

SUBROUTINE Coeff(p,InputFileData, ErrStat, ErrMsg)
!..................................................................................................................................
   ELSE                    ! 3-blader
      p%Hubg1Iner = InputFileData%HubIner
      p%HubInerX = InputFileData%HubIner                ! Line 4799 where the 1st error occured 
      p%Hubg2Iner = 0.0
   ENDIF

! Calculate the tower shape functions (all derivatives):

      p%TwrFASF(1,J,2) = SHP( p%HNodesNorm(J), p%TwrFlexL, InputFileData%TwFAM1Sh(:), 2, ErrStat, ErrMsg )
      p%TwrFASF(2,J,2) = SHP( p%HNodesNorm(J), p%TwrFlexL, InputFileData%TwFAM2Sh(:), 2, ErrStat, ErrMsg )
      p%TwrFASF(1,J,1) = SHP( p%HNodesNorm(J), p%TwrFlexL, InputFileData%TwFAM1Sh(:), 1, ErrStat, ErrMsg )
      p%TwrFASF(2,J,1) = SHP( p%HNodesNorm(J), p%TwrFlexL, InputFileData%TwFAM2Sh(:), 1, ErrStat, ErrMsg )
      p%TwrFASF(1,J,0) = SHP( p%HNodesNorm(J), p%TwrFlexL, InputFileData%TwFAM1Sh(:), 0, ErrStat, ErrMsg )
      p%TwrFASF(2,J,0) = SHP( p%HNodesNorm(J), p%TwrFlexL, InputFileData%TwFAM2Sh(:), 0, ErrStat, ErrMsg )

      p%TwFAM1ShX(:) = InputFileData%TwFAM1Sh(:)     ! Line 5195 where the 2nd error occured 
      p%TwFAM2ShX(:) = InputFileData%TwFAM2Sh(:)

END SUBROUTINE Coeff

!**********************************************************************************************************************************

   SUBROUTINE UpdateCoeff(p, ErrStat, ErrMsg)
!..................................................................................................................................
    ELSE                    ! 3-blader
      p%Hubg1Iner = p%HubInerX
      p%Hubg2Iner = 0.0
   ENDIF

      ! Calculate the tower shape functions (all derivatives):

      p%TwrFASF(1,J,2) = SHP( p%HNodesNorm(J), p%TwrFlexL, p%TwFAM1ShX(:), 2, ErrStat, ErrMsg )

Я сделал это для всех InputFileData% в подпрограмме Coeff (), а также выделил новый определенный параметр в SetOtherParameters () следующим образом:

SUBROUTINE SetOtherParameters( p, InputFileData, ErrStat, ErrMsg )
!..................................................................................................................................
  CALL AllocAry(p%TwFAM1ShX,  PolyOrd,          'TwFAM1ShX' , ErrStat, ErrMsg ); IF ( ErrStat /= ErrID_None ) RETURN

Но когда я начал компилировать, я получил такую ​​ошибку:

Я думаю, что, возможно, я ошибаюсь в защите новых параметров.
Я надеюсь, что смогу объяснить проблему, и был бы рад, если бы вы помогли мне решить эту ошибку.

С наилучшими пожеланиями
Лоуренс

PS: Я взял только некоторые разделы из всего кода каждой подпрограммы, чтобы проиллюстрировать проблему.

Следует иметь в виду, что если вы изменяете какие-либо входные файлы реестра, вам может потребоваться выполнить сборку дважды (первый раз, когда Visual Studio регенерирует файлы * _Types.f90, но код, который зависит от файлов этих типов, не обязательно сборка с правильной версией файла * _Types.f90 до второй сборки). Это может быть причиной первой отображаемой ошибки (если что-то не так).

Одна проблема, которую я вижу, заключается в том, что ваш массив p%TwFAM1SHX не имеет того же размера, что и массив InputFileData%TwFAM1Sh . Вы должны выделить их одинакового размера, если вы собираетесь установить их равными друг другу. Вот распределение, которое я вижу для существующего массива:

CALL AllocAry  ( InputFileData%TwFAM1Sh,  PolyOrd-1, 'TwFAM1Sh'  , ErrStat, ErrMsg )

Ваш отводится на один элемент больше.

Если вы используете текущую ветвь dev , процесс сборки для файлов _Types.f90 изменился: они больше не являются частью стандартного процесса сборки. Вместо этого посмотрите это: https://raf-openfast.readthedocs.io/en/feature-registry_rework/source/dev/types_files.html~~

изменить: см. комментарий @rafmudaf ниже для ссылки.

На самом деле, это лучшая ссылка на документы: https://openfast.readthedocs.io/en/dev/source/dev/types_files.html.

Уважаемый OpenFAST Devloper,
@jjonkman
@bjonkman
Я мог бы изменить массу элементов лопасти в процессе моделирования, и я хотел бы обсудить с вами результат сложения массы.
Я увеличил массу лопаточных элементов, умножив массовую плотность каждого лопаточного элемента в 1,5 раза. Массовое добавление происходило как часть подпрограммы ED-CalcOutput путем вызова подпрограммы UpdateCoeff в заданное время моделирования (в моем примере при t = 0,5 [с]). Я прикрепил файл ElastoDyn.f90, чтобы увидеть все изменения в моем коде.

ElastoDyn.txt

После компиляции OpenFAST я выполнил тест 18 (5MW_Land_DLL_WTurb) и построил массы элементов лезвия для списка узлов лезвий, которые имеют тензодатчики BldGagNd на рисунке 1, эти сигналы (SpnMb) я сгенерировал сам, чтобы увидеть сложение физической массы. На рисунке 2 я сравнил скорость ротора каждой версии OpenFAST и функции OpenFAST + AddMass.

Рисунок_1.pdf

Рисунок_2.pdf

Фактически, я искупил постоянное снижение скорости ротора после добавления массы в 1,5 раза для каждого элемента лопасти, но цифры этого не показывают!
Мои ожидания верны? Правильно ли моя реализация массового сложения?

С наилучшими пожеланиями
Лоуренс

@jjonkman @bjonkman
Еще один вопрос, я хотел бы расширить входной файл ElastoDyn.dat, чтобы использовать увеличивающий коэффициент.
_AddMassF_ и Time _AddMassT_ для массового добавления в качестве InputData, поэтому мне не нужно компилировать OpenFAST каждый раз, когда мне нужно изменить время или коэффициент массового добавления:

Я уже добавил эти входные данные и соответствующие им параметры в ElastoDyn_Registry.txt и в ElastoDyn_IO.f90.
Я мог бы после этого скомпилировать OpenFAST, но когда я запускаю симуляцию, у меня возникает такая ошибка:

Вы уверены, что правильно читаете два новых ввода из файла? Вы добавили две новые строки CALL ReadVar() и поместили их после пробела, где написано PtfmYIner ?

Привет,

Можно ли использовать сигналы от ServoDyn в ElstoDyn? Я вижу, что в ElastoDyn есть сигналы, которые можно использовать в ServoDyn, например:
`` ''
! ................................................. .................................................. ............................
! Выходы, необходимые для ServoDyn
! ................................................. .................................................. ............................

y% Yaw = x% QT (DOF_Yaw)
y% YawRate = x% QDT (DOF_Yaw)
y% YawAngle = x% QT (DOF_Yaw) + x% QT (DOF_Y)! грубое приближение ошибки рыскания ... (без вычитания ее из направления ветра)
y% BlPitch = u% BlPitchCom! OtherState% BlPitch
y% LSS_Spd = x% QDT (DOF_GeAz)
y% HSS_Spd = ABS (p% GBRatio) * x% QDT (DOF_GeAz)
y% RotSpeed ​​= x% QDT (DOF_GeAz) + x% QDT (DOF_DrTr)

Now I’m trying to do exact the opposite. I want to generate a signal in ServoDyn to use it in ElasoDyn.
To do that I defined first an output signal 
`y%BlPitchED(K) = p%BlPitchF(K)`
Then I allocated this signal in SrvD_OutputType data structure (y)

РАСПРЕДЕЛЕНИЕ ВЫЗОВОВ (y% BlPitchED, p% NumBl, 'BlPitchED', ErrStat2, ErrMsg2)
CALL CheckError (ErrStat2, ErrMsg2)
ЕСЛИ (ErrStat> = AbortErrLev) ВОЗВРАТ
`` Afterthat I defended this output in ServoDyn_Registry.txt typedef ^ OutputType ReKi BlPitchED {:} - - "Заданные углы наклона лопастей (ElastoDyn)" радианы`
Но теперь я не знаю, что мне делать, используя этот сигнал в ElstoDyn?
Спасибо за помощь

Уважаемый @LaurenceWETI ,

В структуре модуляции FAST / OpenFAST, когда модули связаны вместе, входы и выходы должны быть согласованными, т. Е. Входы для данного модуля должны быть производными от выходов других модулей. Итак, если вы хотите добавить выход в ServoDyn, который становится входом в ElastoDyn, вы должны определить выход в ServoDyn, определить вход в ElastoDyn, а затем изменить связующий код OpenFAST, чтобы передать выход ServoDyn как вход в ElastoDyn. Например, посмотрите, как выходные данные ServoDyn BlPitchCom передаются в качестве входных данных в ElastoDyn.

С наилучшими пожеланиями,

@jjonkman спасибо за быстрый ответ