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

[ДОБРЯК 602]

4 минуты назад, Fedor сказал:

Плюс от ваших изысканий в том, что теперь не сомневаетесь в правильности написанного мной. :) 

@Fedor вы второй неугомонный Турта. Который то же что-то изобрел и себе это доказывает. 

Эта технология была опубликована в 1964 году. 

Никто догадками до вас не получал базисные функции.

А в целом клоунада про "новую" технологию построения базисных функций была смешная.

Многие даже верили что вы что-то новое изобрели... 

[Fedor 1 597]

Цитата

Эта технология была опубликована в 1964 году

Жаль что этого никто не заметил. Даже вы десятилетия помалкивали. Но ни одной системы базисных функций по ней и не получили.   Но молодец, доказали что все у меня правильно, а то что что-то похожее было раньше, так это мелкий вопрос. Остается  вам найти где опубликовыны и все другие мои результаты   :) 

[ДОБРЯК 602]

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

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

Я постоянно говорил, что вы переписали Зенкевича. А вы обижались. Никто не угадывал функции формы. А использовали эту технологию. 

Вы забываете то что вам говорят. Поэтому я много раз просил вас открыть тему и в ней доказывать что вы что-то открыли.

А вы ровным слоем размазываете своё "изобретение" практически по всем темам. Как Турта...

Вы и этот разговор забудете. 

@Fedor вы же не слышите что вам говорят

.

@Fedor соберите свое изобретение  в одну тему. 

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

 

[Fedor 1 597]

Есть единственный способ сказать что-то - сделать это. Зенкевич тоже переписывал Куранта и Хренникова если уж на то пошло . Не говоря об Эйлере с его вариационным исчислением :)

[Fedor 1 597]

Цитата

 Вы же 20 лет на этом форуме говорили, что я переписалнемецкие хелпы. 

А разве не так ?  :)

[Fedor 1 597]

Цитата

А использовали эту технологию.

Какую эту ?   Попробуйте что-нибудь получить по этим наметкам :)

[ДОБРЯК 602]

31 минуту назад, Fedor сказал:

Какую эту ?

У вас явные проблемы с памятью.

Я же только что дал вам ссылку.

Цитата

 for (i = 0 ; i < 8 ; ++i) {

       switch (i) {

       case  0 : ksi_r = -1 ; nu_r = -1 ; dzeta_r = -1 ; break ;
       case  1 : ksi_r =  1 ; nu_r = -1 ; dzeta_r = -1 ; break ;
       case  2 : ksi_r =  1 ; nu_r =  1 ; dzeta_r = -1 ; break ;
       case  3 : ksi_r = -1 ; nu_r =  1 ; dzeta_r = -1 ; break ;

       case  4 : ksi_r = -1 ; nu_r = -1 ; dzeta_r = 1 ; break ;
       case  5 : ksi_r =  1 ; nu_r = -1 ; dzeta_r = 1 ; break ;
       case  6 : ksi_r =  1 ; nu_r =  1 ; dzeta_r = 1 ; break ;
       case  7 : ksi_r = -1 ; nu_r =  1 ; dzeta_r = 1 ; break ;

       }

 

      uk_qk [i][0] = 1. ;
      uk_qk [i][1] = ksi_r ;
      uk_qk [i][2] = nu_r ;
      uk_qk [i][3] = dzeta_r ;
      

      uk_qk [i][4] = ksi_r * nu_r ;
      uk_qk [i][5] = ksi_r * dzeta_r ;
      uk_qk [i][6] = nu_r * dzeta_r ;

      uk_qk [i][7] = ksi_r * nu_r * dzeta_r ;

   }

   matr_obr (8, uk_qk, uk_m)  ;
 

Вы считаете, что вы это открыли первый в мире. По этой технологии получены функции формы в этих координатах. 

И вы постоянно будете доказывать себе что вы это открыли первый, а до вас угадывали.

Это очень смешная клоунада. И главное бесконечная. Как у Турта.

[Fedor 1 597]

Это вы пытаетесь изобразить .  Ну вызвали и обратили матрицу.  И где базисные функции ?  Хоть бы почитали более позднего Зенкевича  4.4 .  :) 

И исходные базисы из треугольника (пирамиды) Паскаля далеко не всегда срабатывают и позволяют обращать матрицу . Для этого пришлось придумать окаймляющий исходный базис пространства.  :)  

Изменено 19 сентября 2021 пользователем Fedor

[ДОБРЯК 602]

33 минуты назад, Fedor сказал:

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

Вы можете придумывать и изобретать все что угодно. Для треугольника (тетраэдра) функции формы в L-координатах. Ссылку на технологию я вам давал в декабре прошлого года.

281 оболочка.

Интеграл берут по трем точкам. Такая же химия как для 8-ми узлов.

Для этой оболочки условия Кирхгофа не выполняются. 

[Fedor 1 597]

Могу, у нас свободная страна. А такие координаты мне не нравятся. Лучше в обычных как это у Бенерджи с Баттерфильдом . Результат тот же, но единообразнее  с кубиками. 

Интересно глянуть на with,  а не without . Считает то по умолчанию с ключом 0. А у таких функций углы присутствуют дважды.     :) 

Изменено 19 сентября 2021 пользователем Fedor

[ДОБРЯК 602]

42 минуты назад, Fedor сказал:

Могу, у нас свободная страна. А такие координаты мне не нравятся.

Вы и говорите все что угодно на форуме. 

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

[Fedor 1 597]

Вы Зенкевичу не верите который о серендиповых элементах так писал ?  И Митчеллу с Уэйтом, которые писали - " Построение базисных функций является одним из самых важных и часто одним из самых трудных моментов в методе конечных элементов"стр 179 ?  :) 

Знали бы технологию, то о трудностях не писали бы :) 

[Fedor 1 597]

Свободой не прикрываются, в ней живут. Иначе бесплодие.  

Цитата

Петр пришел к пониманию того, что корни европейских технических достижений кроются в высвобождении человеческого разума

из книги Петр Великий Роберта Масси. Жить в Питере и не быть свободным не серьезно, это вам не Москва с ее начальниками   :) 

Изменено 21 сентября 2021 пользователем Fedor

[ДОБРЯК 602]

19.09.2021 в 19:45, Fedor сказал:

Ну вызвали и обратили матрицу.  И где базисные функции ?

На более простом примере. Для 4-ч узловой мембраны Заполнили матрицу и обратили.

Цитата

   for (i = 0 ; i < 4 ; ++i) {
       switch (i) {
          case 0 : ksi_r = -1. ; nu_r = -1. ; break ;
          case 1 : ksi_r =  1. ; nu_r = -1  ; break ;
          case 2 : ksi_r =  1. ; nu_r =  1. ; break ;
          case 3 : ksi_r = -1. ; nu_r =  1. ; break ;
       }

       uk_qk [i][0] = 1. ;
       uk_qk [i][1] = ksi_r ;
       uk_qk [i][2] = nu_r ;
       uk_qk [i][3] = ksi_r * nu_r ;
   }

   matr_obr_new (4, uk_qk, uk_m) ;
 

Далее вектор {1, ksi, nu, ksi  * nu} умножаем на обратную матрицу.

ff [0] = 0.250000  + -0.250000 * ksi + -0.250000 * nu  +  0.250000 * ksi  * nu ;   
ff [1] = 0.250000  +  0.250000 * ksi + -0.250000 * nu  + -0.250000 * ksi  * nu ;   
ff [2] = 0.250000  +  0.250000 * ksi +  0.250000 * nu  +  0.250000 * ksi  * nu ;   
ff [3] = 0.250000  + -0.250000 * ksi +  0.250000 * nu  + -0.250000 * ksi  * nu ;

И получаем функции формы. Также и с производными.

Вы видимо забыли то что придумали. Первый в мире.  
 

[Fedor 1 597]

Дошло и до вас наконец то на пятые сутки как до утки. Долго же соображали  :) 

[ДОБРЯК 602]

9 минут назад, Fedor сказал:

Дошло и до вас наконец то на пятые сутки как до утки. Долго же соображали  :) 

Очень смешная клоунада. Эта технология была известна. Только вы один доказываете, что до вас функции формы угадывали.

Функции формы для 38-ми узлового объемного элемента.

Цитата

ff [0] =   0.171875 + 0.140625 * ksi + 0.140625 * nu  + 0.140625 * dze  + -0.296875 * ksi * nu  + -0.296875 * ksi * dze   + -0.296875 * nu * dze   + 0.296875 * ksi * nu * dze  + -0.171875 * ksi2 + -0.171875 * nu2  + -0.171875 * dze2  + 0.015625 * ksi2 * nu  + 0.015625 * ksi2 * dze   + 0.015625 * nu2 * ksi  + 0.015625 * nu2 * dze  + 0.015625 * dze2 * ksi  + 0.015625 * dze2 * nu  + 0.140625 * ksi2 * nu * dze  + 0.140625 * ksi * nu2 * dze  + 0.140625 * ksi * nu * dze2  + -0.156250 * ksi2 * nu2 * dze  + -0.156250 * ksi * nu2 * dze2  + -0.156250 * ksi2 * nu * dze2  + 0.156250 * ksi2 * nu2 + 0.156250 * nu2 * dze2  + 0.156250 * ksi2 * dze2  + -0.140625 * ksi3 + -0.140625 * nu3 + -0.140625 * dze3 + 0.140625 * ksi3 * nu  + 0.140625 * ksi3 * dze   + 0.140625 * nu3 * ksi  + 0.140625 * nu3 * dze  + 0.140625 * dze3 * ksi  + 0.140625 * dze3 * nu  + -0.140625 * ksi3 * nu * dze  + -0.140625 * ksi * nu3 * dze  + -0.140625 * ksi * nu * dze3  ;

 

ff [1] =   0.171875 + -0.140625 * ksi + 0.140625 * nu  + 0.140625 * dze  + 0.296875 * ksi * nu  + 0.296875 * ksi * dze   + -0.296875 * nu * dze   + -0.296875 * ksi * nu * dze  + -0.171875 * ksi2 + -0.171875 * nu2  + -0.171875 * dze2  + 0.015625 * ksi2 * nu  + 0.015625 * ksi2 * dze   + -0.015625 * nu2 * ksi  + 0.015625 * nu2 * dze  + -0.015625 * dze2 * ksi  + 0.015625 * dze2 * nu  + 0.140625 * ksi2 * nu * dze  + -0.140625 * ksi * nu2 * dze  + -0.140625 * ksi * nu * dze2  + -0.156250 * ksi2 * nu2 * dze  + 0.156250 * ksi * nu2 * dze2  + -0.156250 * ksi2 * nu * dze2  + 0.156250 * ksi2 * nu2 + 0.156250 * nu2 * dze2  + 0.156250 * ksi2 * dze2  + 0.140625 * ksi3 + -0.140625 * nu3 + -0.140625 * dze3 + -0.140625 * ksi3 * nu  + -0.140625 * ksi3 * dze   + -0.140625 * nu3 * ksi  + 0.140625 * nu3 * dze  + -0.140625 * dze3 * ksi  + 0.140625 * dze3 * nu  + 0.140625 * ksi3 * nu * dze  + 0.140625 * ksi * nu3 * dze  + 0.140625 * ksi * nu * dze3  ;

 

ff [2] =   0.171875 + -0.140625 * ksi + -0.140625 * nu  + 0.140625 * dze  + -0.296875 * ksi * nu  + 0.296875 * ksi * dze   + 0.296875 * nu * dze   + 0.296875 * ksi * nu * dze  + -0.171875 * ksi2 + -0.171875 * nu2  + -0.171875 * dze2  + -0.015625 * ksi2 * nu  + 0.015625 * ksi2 * dze   + -0.015625 * nu2 * ksi  + 0.015625 * nu2 * dze  + -0.015625 * dze2 * ksi  + -0.015625 * dze2 * nu  + -0.140625 * ksi2 * nu * dze  + -0.140625 * ksi * nu2 * dze  + 0.140625 * ksi * nu * dze2  + -0.156250 * ksi2 * nu2 * dze  + 0.156250 * ksi * nu2 * dze2  + 0.156250 * ksi2 * nu * dze2  + 0.156250 * ksi2 * nu2 + 0.156250 * nu2 * dze2  + 0.156250 * ksi2 * dze2  + 0.140625 * ksi3 + 0.140625 * nu3 + -0.140625 * dze3 + 0.140625 * ksi3 * nu  + -0.140625 * ksi3 * dze   + 0.140625 * nu3 * ksi  + -0.140625 * nu3 * dze  + -0.140625 * dze3 * ksi  + -0.140625 * dze3 * nu  + -0.140625 * ksi3 * nu * dze  + -0.140625 * ksi * nu3 * dze  + -0.140625 * ksi * nu * dze3  ;
и т. д.

И вы доказываете, что такие функции формы угадывали до вас.

[Fedor 1 597]

Теперь и вы научились. Остается научиться придавать им общепринятый вид в произведениях. Не догадываетесь почему обычно так представляют в книжках ?   :) 

[ДОБРЯК 602]

16 минут назад, Fedor сказал:

Остается научиться придавать им общепринятый вид в произведениях. Не догадываетесь почему обычно так представляют в книжках ? 

Так это вы великий изобретатель. 

Читайте Зенкевича.

Для программирования удобнее такая форма записи.

ff [0] = 0.250000  + -0.250000 * ksi + -0.250000 * nu  +  0.250000 * ksi  * nu ;   
ff [1] = 0.250000  +  0.250000 * ksi + -0.250000 * nu  + -0.250000 * ksi  * nu ;   
ff [2] = 0.250000  +  0.250000 * ksi +  0.250000 * nu  +  0.250000 * ksi  * nu ;   
ff [3] = 0.250000  + -0.250000 * ksi +  0.250000 * nu  + -0.250000 * ksi  * nu ;

[Fedor 1 597]

Ну вот, познакомились с формулой Лагранжа из алгебры :) 

[AlexKaz 1 821]

21.09.2021 в 12:11, ДОБРЯК сказал:

Для программирования удобнее такая форма записи.

Слишком много умножений. Можно сократить до пяти.

Но успехам Вашим можно только радоваться. Сегодня не каждый дипломированный магистрант имеет совесть заглянуть в Зенкевича.

Осталось добраться до главы с переходом из ЛСК в ГСК.