[Cfytrr]

Siemens PLM Femap with NX Nastran "Real FEA made easy"




Тема посвящается вопросам, возникающим при работе в среде программного комплекса Femap

Примеры выполнения некоторых расчетов можно посмотреть:
на домашней страничке

на сайте www.iberisa.com http://iberisa.wordpress.com,
Форум (англ.)
Линейный анализ; Нелинейный анализ
Форум на домашней страничке

[100k]

Сообщение от Cfytrr Так в том то и дело, что методы отбора элементов (узлов) в группу не содержат никаких возможностей использовать свойства View.

Да, программно это сделать не сложно (используя API), как мне кажется на текущий момент.
Попытаюсь если будет время, набросать основу программы.

[Cfytrr]

100k,
Будет отлично, если получится...

Pilot25,
Не надо арматуру доводить до грани бетона, она вся должна быть внутри балки с отступом от граней. И как уже тут говорилось, солиды прикрепляются к балкам (пластинам) только шарнирно, и если арматурный стержень будет прямым, то он будет крутится вокруг своей оси превращая модель в механизм,следовательно, всегда делаем загибы.

p.s. ту выяснился еще один момент, в режиме Solid from Elements фемап строит очень грубую сетку:

Видимо, надо копаться в настройках, но как вариант можно самому добавить "центры кристализации" скопировав элементы арматуры:

и потом их удалив получить нормальную сетку:


Результат в файле, ( что то у Вас там с материалами не хорошо, для стали Е=2,12E+11, а для "бетона" Е=2,8E+11)

[100k]

както так должно выглядеть, в первом приближении, можно протестировать.
не нашел переменную, которая отвечает за выбор Limit Mode.
поэтому пока зашил 4.Outside. выбираются значения которые меньше нижней границы и больше верхней.
Создается группа под номером 1101, куда заносятся все нужные элементы.
Итересно, на сколько быстро работает, если элементов сотни тысяч.

Код:

[Выделить все]

Sub Main

	' переменные
    Dim App As femap.model
    Set App = feFemap()

    Dim v As femap.View
    Set v = App.feView
	Dim vId As Long

	'номер активного вида vid
	rc =App.feAppGetActiveView(vId)

	App.feAppMessage (FCL_BLACK, "Active view = " + Format (vId) )

	'ссылка на объект активного вида
	v.Get(vId)

	Dim os As femap.OutputSet
	Set os = App.feOutputSet
	Dim oput As femap.Output
	Dim osId As Long
	Dim vctId As Long
	Dim cMax, cMin As Double
	Dim vctMode As Long


	osId = v.OutputSet
	App.feAppMessage (FCL_BLACK, "OutPut id = " + Format (osId))
	vctId = v.ContourData
	App.feAppMessage (FCL_BLACK, "Vector countour = " + Format (vctId))
	cMax =v.CriteriaMax
	cMin =v.CriteriaMin
	App.feAppMessage (FCL_BLACK, "Crit min = " + Format(cMin)+ "; max = " + Format(cMax))
	vctMode = v.ContourVecMode
	App.feAppMessage (FCL_BLACK, "Crit mode = " + Format (vctMode))


	os.Get(osId)
	Set oput = os.Vector(vctId)
	Dim count As Long
	Dim el As femap.Elem
	Set el = App.feElem
	Dim nd As femap.Node
	Set nd = App.feNode

	Dim elemSET As femap.Set
	Set elemSET = App.feSet

	Dim gr As femap.Group
	Set gr= App.feGroup
	Dim vl As Double
	Dim entId As Long
	gr.ID = 1101
	gr.title = "Crit"
	'gr.Put(1101)

	'проверка вектора на принадлежность к узлам или элементам
	' помоему если >3000 то это элементы, если меньше узлы.
	'необходимо уточнить
	If(vctId<3000)Then
		count = App.Info_Count(FT_NODE)
	Else
		count = App.Info_Count(FT_ELEM)
	End If

	Debug.Clear
	Debug.Print os.title
	App.feAppMessage (FCL_BLACK, os.title)
	App.feAppMessage (FCL_BLACK, oput.title)


 	For i =1 To count
 		If(vctId<3000)Then
			nd.Next()
			entId = nd.ID
		Else
			el.Next()
			entId =el.ID
		End If
		vl =oput.Value(i)
		'App.feAppMessage (FCL_BLACK, Format(el.ID) +" = " + Format (vl))

		'If(oput.Value(i)> cMin And oput.Value(i)<cMax)Then 'BETWEEN
		If(oput.Value(i)<= cMin Or oput.Value(i)>=cMax)Then 'OUTSIDE
			elemSET.Add(entId)
			gr.Add(FT_ELEM, entId)
			'App.feAppMessage (FCL_BLACK, Format(entId) +" = " + Format(oput.Value(i)))
		End If

 	Next i
	App.feAppMessage (FCL_BLACK, "Selected elems = " + Format(elemSET.Count))
	'gr.Add(FT_ELEM,
	gr.Put(1101)
End Sub

Изменил немного. Должно работать с результатами узлов, проверьте работе ли, с балками, и пластинами.

с балочными элем работает, почему не хочет с солидами ума не приложу

[Pilot25]

Cfytrr
Благодарю! А "бетон" я задал просто наугад

[Cfytrr]

100k,
Что то c солидамии и эквивалентными напряжениями не работает, выбирает в группу все элементы, независимо от ндс и лимитов.

[100k]

Сообщение от Cfytrr 100k, Что то c солидамии и эквивалентными напряжениями не работает, выбирает в группу все элементы, независимо от ндс и лимитов.

файл задачи можно?
с включенными критериями.
я пробовал с балочными элементами.

[Cfytrr]

100k,
Файл

[selega]

Сообщение от Cfytrr 100k, Не надо арматуру доводить до грани бетона, она вся должна быть внутри балки с отступом от граней.

В принципе можно и доводить, только в этом случае придется вручную сместить ближайший узел сетки на поверхности в узел арматуры.
Но у меня возникла другая проблема.
Сделал все как в ролике для создания солидов на базе разбинения поверхностей бетонной части и КЭ арматуры. В итоге у меня и бетон и стержни получились бетонными.
Подскажите, пожалуйста, что я делаю не так? В ролике за операцией создания солидов сразу последовало наложение граничных условий и расчет.
Можно, конечно, потом обратно корректировать свойства КЭ арматуры, но это уже лишние операции, которые сильно затруднят подготовку относительно сложной задачи.

И еще, из личного опыта. Проверку модели на наличие совпадающих узлов лучше проводить дважды - после разбиения кривых, моделирующих арматуру, и после создания сетки на поверхности бетонной части. Без этого FEMAP все время плевался на этапе создания солидов.

[Pilot25]

Согласен по поводу процедуры слияния узлов - если после создания Plot сетки сшить узлы, то фемап успешнее завершит процедуру Solid from Elements. А без этого даже на простейшей сетке иногда завершается неудачей...
А вот по поводу того, что арматуру можно и на поверхность выводить, нужно только совпадение узлов на Plot сетке с арматурными - тут уж все равно не получается... Я делал командой Mesh-Mesh control-Mesh points on surface узлы на поверхностях, но все равно сетка заканчивалась неудачей...

[Cfytrr]

Сообщение от selega В итоге у меня и бетон и стержни получились бетонными.

Не знаю. Такой замены у меня не появляется, арматура остается "железной" и после создания бетона. Возможно, у Вас вкралась ошибка в момент назначения арматуре материала...

[selega]

Сообщение от Pilot25 А вот по поводу того, что арматуру можно и на поверхность выводить, нужно только совпадение узлов на Plot сетке с арматурными - тут уж все равно не получается... Я делал командой Mesh-Mesh control-Mesh points on surface узлы на поверхностях, но все равно сетка заканчивалась неудачей...

Нет, я это делал через Move To. Тогда получается нормально. А без этой процедуры FEMAP писал, что узел расположен на поверхности элемента (по-английски писал, разумеется).

Сообщение от Cfytrr Возможно, у Вас вкралась ошибка в момент назначения арматуре материала..

Так и было. Перебрал модель, и все получилось! Спасибо!

[Pilot25]

Сообщение от selega Нет, я это делал через Move To. Тогда получается нормально. А без этой процедуры FEMAP писал, что узел расположен на поверхности элемента (по-английски писал, разумеется).

Странно, ведь смысл получается тот же - узлы лежат на арматуре... Но ладно - буду дальше разбираться )

[selega]

Сообщение от Pilot25 Странно, ведь смысл получается тот же - узлы лежат на арматуре... Но ладно - буду дальше разбираться )

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

А вообще, даже с учетом нелинейности бетона и арматуры, подобная модель имеет ограниченное применение, потому что в железобетоне есть еще один такой крайне важный и так же крайне сложно учитываемый фактор, как сцепление арматуры с бетоном. Как ни бился над ним М.М.Холмянский и другие исследователи, дальше теоретических выкладок в отдельных монографиях и статьях дело не пошло. В нормах всех стран это явление учитывается эмпирическии коэффициентами. Поэтому с помощью подобной модели можно пытаться учесть, например, деформативность железобетонной бакли с учетом непругой работы бетона, да и то модель в обязательном порядке нужно протестировать на известных примерах. А определять напряжение в бетоне и арматуре в районе трещин в растянутой зоне изгибаемого элемента при абсолютно жесткой связи в узлах сетки стержневых КЭ и солидов - дело попросту бесперспективное. Погрешность может быть очень большой и при этом совершенно непредсказуемой.

[ьфлы]

Добрый день.
Не хочется показаться занудой, ну все-таки …., сделал на бумаги и в фемапе расчеты (простейшая задача)– корректно ли выполнены.
ПС смещение получились разные (сначала на бумаге – более 0,001625 см., потом посчитал в фемапе – 0,000713 см., пересчитал вручную еще раз на бумаге – 0,00875 см.)
я так понимаю 0,000713 и 0,00875 – правильный результат.

[100k]

Сообщение от Cfytrr 100k, Что то c солидамии и эквивалентными напряжениями не работает, выбирает в группу все элементы, независимо от ндс и лимитов.

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

внимание!!!!
'If(oput.Value(i)> cMin And oput.Value(i)<cMax)Then 'BETWEEN
If(oput.Value(i)<= cMin Or oput.Value(i)>=cMax)Then 'OUTSIDE

[Voyager2007]

Сообщение от ьфлы я так понимаю 0,000713 и 0,00875 – правильный результат.

Здравствуйте, я не смотрел модель в Femap. У меня замечание по поводу Вашего решения.
Формула, по которой Вы определяете растяжение/сжатие балки записана так, будто вся балка равномерно сжата,
а деформация то по высоте балки неодинаковая. Задачка оказалась с подковыркой.
Думаю этим объясняется разница в перемещениях с Nastranом. Удачи!

[ьфлы]

Voyager2007 спасибо)
я тоже об этом думал, поэтому и пересчитывал и ответ получился почти единтичный
хотя .... наверное надо было при расчете деформ. в сечении 1 и 2 брать силу не 200 и 450, она ведь не постоянно на данном участке (при 150 и 350 кгс - 0,0007125 см) -
ВСЕ ОК!. хуххх)
Voyager2007 а Вы можете посмотреть сам расчет в фемапе, все ли я корректно оформил)

[Cfytrr]

100k,
Что то, не получается.
Попытался сделать группу для пластин, вот так выглядит limit 0-100


а вот так получившаяся группа:


С балками, еще хуже собирает все в группу не смотря на лимиты
Может от того, что в одном файле и пластины и балки?

[Voyager2007]

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

[ьфлы]

Voyager2007, спасибо.