| Правила | Регистрация | Пользователи | Поиск | Сообщения за день | Все разделы прочитаны |  Справка по форуму | Файлообменник |

Вернуться   Форум DWG.RU > Программное обеспечение > Расчетные программы > ANSYS > Спектр ответа - RESP и FILLDATA

Спектр ответа - RESP и FILLDATA

Ответ
Поиск в этой теме
Непрочитано 24.04.2019, 11:30 #1
Спектр ответа - RESP и FILLDATA
fasol
 
Регистрация: 30.09.2010
Сообщений: 370

Нужно построить спектр ответа в заданной точке конструкции при действии акселерограммы (полноценный нестационарный расчет).
Спектры ответа строю функцией RESP, с ней все понятно:
RESP, IR, номер_таблицы_частот, LDTAB, 3, 0.05, 0.02, 0, 15

На вход она, помимо всего прочего, принимает таблицу частот, заполненной FILLDATA или DATA. С этим сложнее:
FILLDATA, номер_таблицы_частот, LSTRT, LSTOP, LINC, VALUE, DVAL
Пробовал оставлять последние 5 параметров пустые, пробовал заполнять - разницы особой нет.
Построенный спектр очень сильно отличается от результатов работы скадовского Редактора Акселерограмм, если скормить тому поэтажную акселерограмму.
Что я делаю не так?
Просмотров: 3394
 
Непрочитано 24.04.2019, 12:27
#2
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


трудно без контекста советовать..

ktm=1000! точек на графике (можно запросить *get)
fmax=20! от 0.1 до 20 гц частоты учитываем
dfr=fmax/ktm! шаг по частотам
plvar,2! график ускорений
filldata,10,,,,0.1,dfr
plvar,10!проверьте таблицу

dt=0.02! какой шаг у вас ?
xvar,10! по абсциссе дб частоты

demp=0.3/6.28! при декременте 0.30

resp,3,10,2,3,demp,dt,,,1! спектр ускорений
plvar,3

Последний раз редактировалось СергейД, 24.04.2019 в 12:35. Причина: проверка
СергейД вне форума  
 
Непрочитано 24.04.2019, 17:31
#3
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


Эта задача решается элементарно в екселе. Строчек на бейсике будет не больше, чем в ансисе.
__________________
Сон разума рождает чудовищ. (испанская пословица)
Бахил вне форума  
 
Непрочитано 24.04.2019, 18:17
#4
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


тут несколько строк просто для читабельности
в принципе всего- три строки.

filldata,10,,,,0.1,1000
xvar,10! по абсциссе дб частоты
resp,3,10,2,3,0.05,0.02,,,1! спектр ускорений

все троллишь, Бахил?
СергейД вне форума  
 
Непрочитано 25.04.2019, 08:08
#5
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


Цитата:
Сообщение от СергейД Посмотреть сообщение
все троллишь, Бахил?
Да. На бейсике побольше будет.

зато более понятно
__________________
Сон разума рождает чудовищ. (испанская пословица)

Последний раз редактировалось Бахил, 26.04.2019 в 15:41.
Бахил вне форума  
 
Непрочитано 25.04.2019, 08:28
#6
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


кому как...
ты уверен, что это то что нужно???

у человека макрос под ансис со всей процедурой расчета. модель-расчет-обработка
речь идет о том, чтобы в нее добавить несколько строк. точнее, поправить.

зачем тут что-то передавать в эксель. там программировать, потом обратно в ансис загружать???
я понимаю еще питон применить.
СергейД вне форума  
 
Непрочитано 25.04.2019, 08:31
#7
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


Offtop: Да шучу я. Не бери в голову.
__________________
Сон разума рождает чудовищ. (испанская пословица)
Бахил вне форума  
 
Непрочитано 25.04.2019, 11:28
#8
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


кстати, планируем организовать обучение ansys на курсах.
но я только на базе apdl смогу научить..
WB потом самому можно изучить.
СергейД вне форума  
 
Автор темы   Непрочитано 26.04.2019, 11:49
#9
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Цитата:
Сообщение от СергейД Посмотреть сообщение
filldata,10,,,,0.1,1000
xvar,10! по абсциссе дб частоты
resp,3,10,2,3,0.05,0.02,,,1! спектр ускорений
Спасибо, стало намного лучше.
Как оказалось, на спектр сильно влияет практически все, а особенно максимальная частота и количество точек в графике, см. вложение.
Причем при увеличении количество точек (больше 1000), становится еще хуже - пик переезжает в другую половину графика.
В итоге наиболее точно ансис со скадом сходится при максимальной частоте 25 Гц и количестве точек как в акселерограмме (756).
Там же просто разложение Фурье, откуда такая чувствительность ко входным данным?
Миниатюры
Нажмите на изображение для увеличения
Название: ansys_scad.jpg
Просмотров: 64
Размер:	170.3 Кб
ID:	213595  
fasol вне форума  
 
Непрочитано 26.04.2019, 12:00
#10
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


Разный масштаб абсциссы.
__________________
Сон разума рождает чудовищ. (испанская пословица)
Бахил вне форума  
 
Автор темы   Непрочитано 26.04.2019, 12:00
#11
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Цитата:
Сообщение от Бахил Посмотреть сообщение
Разный масштаб абсциссы.
Лишь бы ляпнуть, да?
fasol вне форума  
 
Непрочитано 26.04.2019, 12:06
#12
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


По твоим графикам только одно можно сказать: СКАД считает правильно.
А с ансисом ты не разобрался.

----- добавлено через ~2 мин. -----
Да, и судя по графикам, по абсциссе отложена не частота, а период.
Можешь проверить по элементарной процедуре. Ну если разберёшься.
__________________
Сон разума рождает чудовищ. (испанская пословица)
Бахил вне форума  
 
Непрочитано 26.04.2019, 13:11
#13
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


число точек = число шагов решения=756
по горизонтали должны быть частоты.

вышли в личку акселерограмму.

Последний раз редактировалось СергейД, 26.04.2019 в 13:23.
СергейД вне форума  
 
Автор темы   Непрочитано 26.04.2019, 13:29
#14
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Выслал.
fasol вне форума  
 
Непрочитано 26.04.2019, 14:32
#15
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


Offtop: ...
__________________
Сон разума рождает чудовищ. (испанская пословица)

Последний раз редактировалось Бахил, 26.04.2019 в 15:42.
Бахил вне форума  
 
Непрочитано 29.04.2019, 13:48
#16
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


все очень хорошо совпало со скадовскими графиками.
СергейД вне форума  
 
Автор темы   Непрочитано 12.05.2019, 14:12
#17
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Методом тыка было установлено, что максимальная частота никак не влияет на график спектра.
Заметно влияет шаг RESP, IR, номер_таблицы_частот, LDTAB, 3, 0.05, 0.02.
Причем при совпадении шага с дискретностью акселерограммы (0.02) сходимость со СКАДом наименьшая.

fasol вне форума  
 
Непрочитано 12.05.2019, 16:45
#18
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


в личном письме напсиал, что шаг 0.02 грубый для прямой обработки акселерограммы в 10 гц. она негладкая. с изломами.
при уменьшении шага до 0.005с и менее точность разложения в ансис возрастает до необходимой. проверено на идеальной синусоиде.
СергейД вне форума  
 
Автор темы   Непрочитано 13.05.2019, 08:42
#19
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Цитата:
Сообщение от СергейД Посмотреть сообщение
при уменьшении шага до 0.005с и менее точность разложения в ансис возрастает до необходимой
Я как раз это на графике и привел.
0.005 почти не отличается от 0.01, фаза практически совпадает со скадовской, а вот амплитуда на 20% меньше.
Дальнейшее уменьшение шага ухудшает картину.
fasol вне форума  
 
Непрочитано 13.05.2019, 09:03
#20
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


Интересно взглянуть на акселерограмму.
Шаг вообще-то не должен сильно влиять.
__________________
Сон разума рождает чудовищ. (испанская пословица)
Бахил вне форума  
 
Непрочитано 13.05.2019, 10:20
#21
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


там частота 10 гц. период 0.1 с. при 0.02 с пять точек на две полуволны. естественно, график с изломами. пики пропускаются.

при разложении в спектр разумеется будет ошибка.оттого и максимумы меньше.


при 0.005 20 точек на период- график более менее гладкий.
я просто для наглядности сгенерировал синусоиду 10 гц и сохранил с различными шагами по времени.
все очевидно.

в скад видимо работали с исходным вариантом с меньшим шагом или создавали сплайн-аппроксимацию между точками.
СергейД вне форума  
 
Автор темы   Непрочитано 13.05.2019, 10:43
#22
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Цитата:
Сообщение от СергейД Посмотреть сообщение
я просто для наглядности сгенерировал синусоиду 10 гц и сохранил с различными шагами по времени.
Покажите, если не трудно, как у вас сходятся спектры акселерограммы, что высылал вам.
Аппроксимации в скаде точно нет.
fasol вне форума  
 
Непрочитано 13.05.2019, 13:48
#23
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


про скад откуда вы знаете что у них внутри?
исходный график должен быть похож без изломов и не пропускать экстремумы (5 точек на две полуволны- это же недопустимо) , тогда и частота точнее получится в спектре.
файл в личку вышлю. сейчас занят.
СергейД вне форума  
 
Автор темы   Непрочитано 13.05.2019, 18:14
#24
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Цитата:
Сообщение от СергейД Посмотреть сообщение
про скад откуда вы знаете что у них внутри?
Потому что точек на графике там заметно меньше.

Синусоида у меня тоже более-менее сходится (шаг акселерограммы 0.02).
fasol вне форума  
 
Автор темы   Непрочитано 17.05.2019, 08:37
#25
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Есть у кого-нибудь возможность прогнать акселерограмму в Одисее или схожем ПО?
Появилось подозрение, что Скад временами привирает.
Акселерограмма - https://pastebin.com/DzeVrphU
Результирующий скадовский спектр 5% - https://pastebin.com/TKMX3nXP
fasol вне форума  
 
Непрочитано 17.05.2019, 09:10
#26
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


опять 0.02...
при таком грубом шаге будут все те же коллизии

затребуйте с меньшим шагом или сами ее интерполируйте сплайном
СергейД вне форума  
 
Автор темы   Непрочитано 17.05.2019, 09:51
#27
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Цитата:
Сообщение от СергейД Посмотреть сообщение
опять 0.02...
Это поэтажная акселерограмма с шагом как у исходной.
Все интерполяции уже при обработке, иначе смысл.
fasol вне форума  
 
Непрочитано 18.05.2019, 22:08
#28
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


11,6 Гц резонанс 302
Бахил вне форума  
 
Непрочитано 19.05.2019, 04:36
#29
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


в ансис обрабатывается график с линейным изменением между точками.
если шаг 0.02с (50 гц) а функция изменяется в основном 10-12 гц,
то как вы две полуволны правильно опишете пятью точками? будут изломы графика.

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

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


2 бахил. вы по что сейчас?. своей программой считали или какой?
создайте синусоиду10 гц, запишите ее с шагом 0.02с и определите ее спектр при нулевой и ненулевом демпф. что получится?

Последний раз редактировалось СергейД, 19.05.2019 в 04:41.
СергейД вне форума  
 
Автор темы   Непрочитано 19.05.2019, 11:52
#30
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Цитата:
Сообщение от СергейД Посмотреть сообщение
получите той или иной интерполяцией перед ансис дополнительные точки на графике и считайте себе спектр в ансис.
Попробую.

Цитата:
Сообщение от СергейД Посмотреть сообщение
в ноль остаточные перемещения приходят?
Нет, акселерограмма не сбалансированная, как и любая синтезированная. Но это никак не влияет на спектр, я проверял.
fasol вне форума  
 
Непрочитано 20.05.2019, 08:29
#31
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


Offtop: Кнопкотыки...
0,02 - шаг периода для которого считается максимальная реакция линейного осциллятора.
Для данной акселерограммы максимальная реакция осциллятора= 302 при его частоте 11,6 Гц. Для периода 0,1с примерно 297.
Странная акселерограмма. Бункер считаешь?
__________________
Сон разума рождает чудовищ. (испанская пословица)

Последний раз редактировалось Бахил, 20.05.2019 в 10:03.
Бахил вне форума  
 
Непрочитано 20.05.2019, 08:53
| 1 #32
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


акс-мму нужно интерполировать сплайном, введя точек 5 промежуточных.
тогда ее можно спокойно применять.

вопрос и в использовании этой акс-ммы для расчетов отклика.
какие частоты у нагружаемого оборудования?
волновать должны пики спектра нагрузки в зоне низших собственных частот.
если эти частоты 30-50 гц, то разницы при пиках спектров 9 или 11 гц особой не заметишь.
СергейД вне форума  
 
Непрочитано 20.05.2019, 09:32
#33
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


Да.
__________________
Сон разума рождает чудовищ. (испанская пословица)
Бахил вне форума  
 
Автор темы   Непрочитано 23.05.2019, 12:32
#34
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Цитата:
Сообщение от СергейД Посмотреть сообщение
получите той или иной интерполяцией перед ансис дополнительные точки на графике и считайте себе спектр в ансис.
Сделал линейную интерполяцию результирующей сейсмограммы, увеличив количество точек в два раза.
Потом через псевдозадачу как у вас снял NSOL и опять прогнал RESP - результат еще хуже, фаза на месте, а вот амплитуда просела на 20%.
Слишком трудоемкая задача оказалась, возвращаюсь к инструментам скада.


Цитата:
Сообщение от СергейД Посмотреть сообщение
какие частоты у нагружаемого оборудования?
волновать должны пики спектра нагрузки в зоне низших собственных частот.
если эти частоты 30-50 гц, то разницы при пиках спектров 9 или 11 гц особой не заметишь.
Такой подход меня не устраивает, это уже шаг назад.
fasol вне форума  
 
Непрочитано 23.05.2019, 12:35
#35
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


речь идет о нелинейной интерполяции графика. гладким сплайном. от линейной действительно только хуже будет.
СергейД вне форума  
 
Автор темы   Непрочитано 23.05.2019, 13:19
#36
fasol


 
Регистрация: 30.09.2010
Сообщений: 370


Цитата:
Сообщение от СергейД Посмотреть сообщение
гладким сплайном.
А функционал APDL позволяет это сделать?
Этож нужно ловить несколько точек до и несколько после рассматриваемого интервала и аппроксимировать полиномом?
fasol вне форума  
 
Непрочитано 23.05.2019, 14:53
#37
СергейД

расчеты МКЭ и CFD. ктн
 
Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,089


сделайте в матлабе или еще где-то.

в ансис можно изловчиться-извратиться= вариантов много.
например, можно цифры в координаты KP превратить и провести сплайн (можно до 200 точек в одной порции). разбить линии ldiv и сохранить координаты всех точек.
СергейД вне форума  
 
Непрочитано 23.05.2019, 15:51
#38
Бахил

?
 
Регистрация: 17.06.2014
Царицын
Сообщений: 8,145


Интерполировать исходную акселерограмму бесполезно
Бахил вне форума  
Ответ
Вернуться   Форум DWG.RU > Программное обеспечение > Расчетные программы > ANSYS > Спектр ответа - RESP и FILLDATA

Размещение рекламы
Опции темы Поиск в этой теме
Поиск в этой теме:

Расширенный поиск


Похожие темы
Тема Автор Раздел Ответов Последнее сообщение
Как построить спектр ответа конструкции при заданной акселерограмме(Abaqus)? ИванМоскалев Расчетные программы 55 02.12.2014 12:46
ExAc 0.3: Экспорт поэтажных акселерограмм и спектров ответа. Александр Бауск Robot 2 16.07.2011 03:32
Поэтажные спектры ответа от реальных акселеллограм (Lira/ SCAD) DKI Расчетные программы 2 13.07.2010 09:40
Программа для спктров ответа diesel-generator Расчетные программы 3 22.08.2009 00:20
Прибор диагностики свай СПЕКТР 2.0 Danger_pgs Разное 9 30.01.2009 06:01