Сражения с ветряными мельницами

[ДОБРЯК 603]

Для 10-ти элементов

  1        11   0.0000E+00   1.3458E-03   1.0312E-03   0.0000E+00  -1.4308E-03   2.0348E-03
Результаты разные только потому, что разная нагрузка. Опорные реакции я уже приводил. Попробуйте найти в этой помойке

[piden 2 874]

30 minutes ago, Борман said:

реши в сисе эту задачу.

А то я что-то не осилил правильный ответ. 

Та хз, вроде так же:

Spoiler

 

Здесь балка отображается, как ступенчатая, но все там tapered)

 

 

 

Spoiler

fini
/clear

/prep7  
k,  
k,,1000 

l,1,2
ET,1,BEAM188
 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2e5
MPDATA,PRXY,1,,0.3    

SECTYPE,   1, BEAM, RECT, big, 0
SECOFFSET, CENT 
SECDATA,100,100,

SECTYPE,   2, BEAM, RECT, small, 0  
SECOFFSET, CENT 
SECDATA,100,50,

SECTYPE,3,TAPER, ,taper 

SECDATA,   1,0,,,   
SECDATA,   2,1000,,,

secnum,3

/ESHAPE,1.0 
ESIZE,0,10, 		! ESIZE, SIZE, NDIV

LMESH,1

DK,1,ALL

gplot


!*** line pressure application

*get,my_elem_min,elem,,num,min
*get,my_elem_max,elem,,num,max

delta=my_elem_max-(my_elem_min-1)


*DO,i,my_elem_min,my_elem_max,1

esel,s,,,i
gplot
! SFBEAM,all,1,PRES, 13-(i-1)/delta*3, 13-(i)/delta*3, 		! formula loading, reversed
SFBEAM,all,1,PRES, 10+(i-1)/delta*3, 10+(i)/delta*3, 		! formula loading, normal
! /wait,0.2
*enddo

allsel
gplot

sfeli


/solu
solve

/post1
set,last
plnsol,u,z

 

 

[karachun 2 038]

1 час назад, piden сказал:

Но (как выяснилось!) в справке мелкими буквами написано, что второе сечение будет расположено не на другом конце балки, а на другом конце элемента (!!!). И если в балке больше одного элемента, то у нее будет пилообразная форма.

Да, В Femap/Nastran так же, если замешить потом кривую то получится пила.

@LazyBitch Но в Фемап, в Custom Tools, есть инструмент по правильному разрезанию таких балок. Видимо кто-то из пользователей часто такие балки считал и ему надоело по десять раз создавать сечения.

Вот, снял небольшое видео. https://fex.net/ru/s/pyrtaac

Инструмент находится по пути Custom Tools->Element Update -> Split Tapered Beams.

Изменено 29 мая 2022 пользователем karachun

[Борман 2 375]

@piden

спасибо. С точностью до капса исходник такой же.

 

@piden

а вспомнил.. у меня что-то не рисует поля на такой /ESHAPE,ON - балке 188.

 

Ты то как нарисовал ?

[piden 2 874]

7 minutes ago, Борман said:

Ты то как нарисовал ?

/GRAPHICS,POWER

13 minutes ago, Борман said:

/ESHAPE,ON

У меня срабатывает /eshape, 1

 

3 hours ago, ДОБРЯК said:

Чтобы получить реакцию 11 500 нужно приложить силу 11 500. В этом случае момент будет = 11 500, а у вас 6 250.

Распределенная сила - она такая)

 

Постойте! Так может, в терминологии все причины? Мы же тут имеем дело не абы с кем! 

 

@ДОБРЯК, а в твоем мега-комплексе можно к балке векторное давление прикладывать? Если можно, то вот у тебя был шанс, и ты им не воспользовался)

[Борман 2 375]

@piden

Не хочет чото.

Давно заметил,  что на 188 элементе какие то проблемы с этим. Какая то хитрость видимо есть.

[статист 609]

1 час назад, ДОБРЯК сказал:

Для 10-ти элементов

 

 

                                       Нагрузка вдоль Y                                                 Нагрузка вдоль Z                           

                            (вдоль меняющегося размера)                             (поперек меняющегося размера)

Аналитика                       0.0013401                                                                      0.00102756

C_A                                  0.0013412                                                                       0.00102836

Отклонение                    0.082%                                                                            0.078%

ISPA                                 0.0013458                                                                       0.0010312

Отклонение ISPA            0.425%                                                                           0.354

Lira                                  0.0013393                                                                       0.0010269

Отклонение Lira             0.06%                                                                              0.064%

Ansys                               0.0013592

Отклонение Ansys          1.425%

Abaqus                             0.0013566

Отклонение Abaqus        1.23%

 

Пока получается вот такая картина. Лира лучше всех себя на этом тесте ведет.

 

[LazyBitch 275]

2 часа назад, karachun сказал:

@LazyBitch Но в Фемап, в Custom Tools, есть инструмент по правильному разрезанию таких балок. Видимо кто-то из пользователей часто такие балки считал и ему надоело по десять раз создавать сечения.

Вот, снял небольшое видео. https://fex.net/ru/s/pyrtaac

К сожалению, видео не смогла посмотреть, сайт не грузит даже с vpn. @karachun, скажите, пожалуйста, инструментом Split Tapered Beams можно за один проход порезать балку на N-ое количество одинаковых интервалов?
Кстати, если работать этим инструментом с несколькими балками, для которых изначально была назначена одинаковая жесткость, то инструмент багует. Чтобы этого не происходило надо каждой балке продублировать жёсткость под уникальным id.

Скрытый текст

Мухин, держи результаты Настрана. Поскольку смогла резать балки только пополам, вместо 10-ти эл. у меня 8-мь, а вместо 30-ти соответственно - 32.

Скрытый текст

 

Изменено 29 мая 2022 пользователем LazyBitch

[статист 609]

 

В Autodesk Nastran получается вот такая вот картинка

 

 

 

0,00134972 м

 

 

 

0,001039 м

 

Но там для меняющихся сечений используются только beam, а они основаны на балке Тимошенко, а не Эйлера. Господа, может эти проценты отклонений с этим связаны?

 

 

                                      Нагрузка вдоль Y                                                 Нагрузка вдоль Z                           

                            (вдоль меняющегося размера)                             (поперек меняющегося размера)

Аналитика                       0.0013401                                                                      0.00102756

C_A                                  0.0013412                                                                       0.00102836

Отклонение                    0.082%                                                                            0.078%

ISPA                                 0.0013458                                                                       0.0010312

Отклонение ISPA            0.425%                                                                           0.354%

Lira                                  0.0013393                                                                       0.0010269

Отклонение Lira             0.06%                                                                              0.064%

Ansys                               0.0013592

Отклонение Ansys          1.425%

Abaqus                             0.0013566

Отклонение Abaqus        1.23%

Autodesk Nastran            0.00134972                                                                    0.001039

Отклонение AN               0,717%                                                                           1,11%

Изменено 29 мая 2022 пользователем статист

[ДОБРЯК 603]

1 час назад, статист сказал:

Лира лучше всех себя на этом тесте ведет.

Вы устроили очередное шоу. И считаете, что в Ансис работают глупые разработчики и не только в Ансис. 

Продолжайте и дальше свое шоу. Это ваш звездный час.

[статист 609]

34 минуты назад, ДОБРЯК сказал:

Продолжайте и дальше свое шоу. Это ваш звездный час.

 

[karachun 2 038]

2 часа назад, LazyBitch сказал:

Split Tapered Beams можно за один проход порезать балку на N-ое количество одинаковых интервалов?

Получается что за один проход нельзя. Не очень удобно но хоть что-то. А что касается видео то я там никаких секретов не открыл, просто показал как этот инструмент работает.

Изменено 30 мая 2022 пользователем karachun

[ДОБРЯК 603]

@статист для продолжения шоу сравните результаты решения для одного КЭ) в Ансис и строительном софте. И потом сравните с аналитическим решением.

И напряжения сравните и силовые факторы...

Полное сравнение проведите...

[Борман 2 375]

Вообще тут хитрость есть.

У балки статиста осевая линия имеет наклон 0,025/1. А значит если моделировать осевую линию, то вертикальная нагрузка действует не поперек оси балки со всеми вытекающими.

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

 

Понятно, что если все сводить к простой балке с переменным сечением, то все эти вопросы отпадают :)

 

Мутная тема какая-то.

[Fedor 1 597]

Цитата

За идею ничего не прошу

А она ничего и не стоит. Копаться в чужих кодах то еще занятие... Вон в линуксе копаются десятилетия, а все равно не особо чего хорошего получается. Лучше уж государству дать денег автору, с условием чтобы сделал коды открытыми и ПО бесплатным. Может кто и захочет покопаться   :) 

Изменено 30 мая 2022 пользователем Fedor

[piden 2 874]

7 hours ago, статист said:

там для меняющихся сечений используются только beam, а они основаны на балке Тимошенко, а не Эйлера. Господа, может эти проценты отклонений с этим связаны?

Похоже. В абакусе есть элемент B23 (beam 2d, cubic), он уже Euler-Bernulli, с ним прогиб на 10 элементах -1.3452. Но он доступен только для 2d.

 

19 minutes ago, Борман said:

со всеми вытекающими.

Вроде как да, но в 3д для солидов оба варианта сходятся к одному ответу. Да и у тебя для трех вариантов (+ ступенчатое сечение) тоже перемещения совпали.

Spoiler

 

[Борман 2 375]

@piden

я когда считал солидами, то получил ответ как у статиста. Отсюда и непонятки.

[ДОБРЯК 603]

6 минут назад, piden сказал:

В абакусе есть элемент B23 (beam 2d, cubic), он уже Euler-Bernulli

Теория стержня Бернулли для постоянного сечения. Для переменного сечения нет такой теории.

Или опять очередная новизна на форуме.

Так в этой новизне строительный софт впереди планеты всей. Хочешь 10 узлов в элементе задавай, можно и двадцать, можно и сто. 

А конечный элемент все равно 2-х узловой.)

8 минут назад, Борман сказал:

Отсюда и непонятки.

А вы задайте моменты инерции для одного стержня, а не геометрию для сечения. Результат не изменится. 

И не забывайте про h в кубе.

И в Лире задайте моменты инерции для одного КЭ и сравните результаты. 

Тогда перестанут может) кричать, что строительный софт впереди планеты всей.

[Борман 2 375]

@ДОБРЯК реши уже что-нибудь. Что тебе из задачи непонятно то ???

[piden 2 874]

3 minutes ago, Борман said:

реши уже что-нибудь.

Да он вроде сделал:

10 hours ago, ДОБРЯК said:

Для 10-ти элементов

  1        11   0.0000E+00   1.3458E-03   1.0312E-03   0.0000E+00  -1.4308E-03   2.0348E-03

 

Только, как обычно в его практических задачах, все как-то тихонько, без картинок, "джентльмены верят на слово" итп.