Поиск сил при помощи перемещений

[rasta89 113]

Всем привет.

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

как бы все ок, максимальную реакцию я увидел 194Н, чтобы проверить пошел от обратного, задал силу 200Н на внутреннюю часть лепестка слева, мол выдернет, а вот нифига, не хватает силы чтобы вырвать из паза, почему такое происходит, кто подскажет уважаемые?

[karachun 2 038]

@rasta89 А если приложить и перемещение и силу к одному узлу?

 

P.S. Сделайте лепесток гексой.

[rasta89 113]

3 минуты назад, karachun сказал:

@rasta89 А если приложить и перемещение и силу к одному узлу?

хм, но зачем? я же ищу при какой силе сорвется, что-то я не въехал(

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

@rasta89

P.S. Сделайте лепесток гексой.

ок, сделаю

[Nikon6260 20]

Руками не пробовали посчитать для начала? 

На 200 Н совсем не вылезла или почти не вылезла? Может побольше силу дать? Невязки все же есть в МКЭ. 

[rasta89 113]

3 минуты назад, Nikon6260 сказал:

Руками не пробовали посчитать для начала? 

На 200 Н совсем не вылезла или почти не вылезла? Может побольше силу дать? Невязки все же есть в МКЭ. 

1. руками считал конечно и есть данные испытаний, в принципе расчет в МКЭ близок к реальному

2. взял один из 22 лепестков для ускорения расчета, считал а потом складывал нагрузки, от испытательных отличается на ~100Н

3. лепесток почти вылез миллиметра не хватило, вот и думаю о погрешности 10-30Н не хватает от рассчитанного при помощи перемещений

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

[Nikon6260 20]

Отличие от опыта может быть из-за коэффициента трения в контакте. Как правило это самый неизвестный параметр. Вот и подберёте за одно и коэффициент трения. 

[karachun 2 038]

25 минут назад, rasta89 сказал:

а там меняются свойства стали

нет, модуль упругости не меняется

34 минуты назад, rasta89 сказал:

что-то я не въехал(

вьезжать не надо, делайте

[Fedor 1 597]

Цитата

А если приложить и перемещение и силу к одному узлу?

Уравнения  (е) с заданным перемещениями  (м)  вычеркиваются ввиду излишности. В остальных подставляются значения и  уходят в правую часть. То есть силу можете любую приложить и ничего не изменится. Логично приложить нулевую :) 

[Борман 2 375]

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

[karachun 2 038]

Вот по этому и надо было тянуть за один узел но мне было лень объяснять.

Будет интересно если это не поможет.

[Борман 2 375]

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

но мне было лень объяснять

Тогда нет смысла что-либо писать на форум. А умничать и перед зеркалом можно.

[rasta89 113]

5 часов назад, karachun сказал:

Вот по этому и надо было тянуть за один узел но мне было лень объяснять.

Будет интересно если это не поможет.

Ну я и тянул за один узел, за бортик слева, трение задано 0.1, получается лепесток в пазе соприкасается создаётся усиле чтобы его оттуда дёрнуло...а в обратном расчете в этот же узел, в этом же месте прикладывал силу

[karachun 2 038]

Тогда перейдем от слов к делу.

Прикладывайте модель.

Изменено 14 октября 2021 пользователем karachun

[rasta89 113]

17 часов назад, karachun сказал:

Тогда перейдем от слов к делу.

Прикладывайте модель.

Приложил проект.

Получается так. Что путем поиска силы срыва при помощи перемещений показало что нужно для срыва 100Н. При обратной проверке прикалдыванием силы и последовательных расчетах сорвало только на 150Н, с разбиением на 10 шагов. Получилось что погрешность 50%((

export1.wbpz

[rasta89 113]

Перемещение протестил с разделением на 10 шагов, получилось 127Н на срыв, итого разница:

При перемещении 127Н

При нагружении 150Н

[AlexKaz 1 821]

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

При нагружении 150Н

Не, между 128 и 144.

 

 Уровень по силам выше т.к. и невязка в решении растёт с ростом НДС и к 128 достигает примерно 4 Н. Чем больше узлов, тем хуже сходимость и выше невязка.

 

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

    FORCE CONVERGENCE VALUE  =   540.7      CRITERION=  0.5614   
     DISP CONVERGENCE VALUE   =  0.1939E-05  CRITERION=  0.1537E-04 <<< CONVERGED
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1939E-05
     DISP CONVERGENCE VALUE   =  0.1939E-05  CRITERION=  0.1569E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC =  0.1939E-05
     FORCE CONVERGENCE VALUE  =  0.9115E-01  CRITERION=  0.5725     <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   1
 *** LOAD STEP     1   SUBSTEP     2  COMPLETED.    CUM ITER =      8
 *** TIME =  0.200000         TIME INC =  0.100000   

     FORCE CONVERGENCE VALUE  =   535.6      CRITERION=  0.8413   
     DISP CONVERGENCE VALUE   =  0.8428E-05  CRITERION=  0.1537E-04 <<< CONVERGED
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.8428E-05
     DISP CONVERGENCE VALUE   =  0.8414E-05  CRITERION=  0.1569E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =  0.9983     SCALED MAX DOF INC =  0.8414E-05
     FORCE CONVERGENCE VALUE  =   1.218      CRITERION=  0.8518   
    EQUIL ITER   2 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2045E-05
     DISP CONVERGENCE VALUE   =  0.2045E-05  CRITERION=  0.1601E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC = -0.2045E-05
     FORCE CONVERGENCE VALUE  =  0.1654      CRITERION=  0.8715     <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   2
 *** LOAD STEP     1   SUBSTEP     3  COMPLETED.    CUM ITER =     10
 *** TIME =  0.300000         TIME INC =  0.100000   

     FORCE CONVERGENCE VALUE  =   515.1      CRITERION=   1.113   
     DISP CONVERGENCE VALUE   =  0.3069E-05  CRITERION=  0.1537E-04 <<< CONVERGED
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.3069E-05
     DISP CONVERGENCE VALUE   =  0.3069E-05  CRITERION=  0.1569E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC = -0.3069E-05
     FORCE CONVERGENCE VALUE  =  0.2479      CRITERION=   1.139     <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   1
 *** LOAD STEP     1   SUBSTEP     4  COMPLETED.    CUM ITER =     11
 *** TIME =  0.400000         TIME INC =  0.100000   

     FORCE CONVERGENCE VALUE  =   523.8      CRITERION=   1.395   
     DISP CONVERGENCE VALUE   =  0.3313E-04  CRITERION=  0.1683E-04
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.3313E-04
     DISP CONVERGENCE VALUE   =  0.2996E-04  CRITERION=  0.1717E-04
     LINE SEARCH PARAMETER =  0.9042     SCALED MAX DOF INC = -0.2996E-04
     FORCE CONVERGENCE VALUE  =   52.89      CRITERION=   1.452   
    EQUIL ITER   2 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.2482E-04
     DISP CONVERGENCE VALUE   =  0.2482E-04  CRITERION=  0.1752E-04
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC =  0.2482E-04
     FORCE CONVERGENCE VALUE  =   3.316      CRITERION=   1.457   
    EQUIL ITER   3 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1740E-04
     DISP CONVERGENCE VALUE   =  0.1326E-04  CRITERION=  0.1788E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =  0.7617     SCALED MAX DOF INC =  0.1326E-04
     FORCE CONVERGENCE VALUE  =   2.262      CRITERION=   1.474   
    EQUIL ITER   4 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.7478E-07
     DISP CONVERGENCE VALUE   =  0.7478E-07  CRITERION=  0.1824E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC = -0.7478E-07
     FORCE CONVERGENCE VALUE  =  0.1555      CRITERION=   1.504     <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   4
 *** LOAD STEP     1   SUBSTEP     5  COMPLETED.    CUM ITER =     15
 *** TIME =  0.500000         TIME INC =  0.100000   

     FORCE CONVERGENCE VALUE  =   494.4      CRITERION=   1.658   
     DISP CONVERGENCE VALUE   =  0.2113E-04  CRITERION=  0.1683E-04
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.2113E-04
     DISP CONVERGENCE VALUE   =  0.1974E-04  CRITERION=  0.1717E-04
     LINE SEARCH PARAMETER =  0.9343     SCALED MAX DOF INC =  0.1974E-04
     FORCE CONVERGENCE VALUE  =   33.69      CRITERION=   1.675   
    EQUIL ITER   2 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2309E-04
     DISP CONVERGENCE VALUE   =  0.2309E-04  CRITERION=  0.1752E-04
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC = -0.2309E-04
     FORCE CONVERGENCE VALUE  =   3.184      CRITERION=   1.732   
    EQUIL ITER   3 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.1356E-04
     DISP CONVERGENCE VALUE   =  0.9824E-05  CRITERION=  0.1788E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =  0.7247     SCALED MAX DOF INC = -0.9824E-05
     FORCE CONVERGENCE VALUE  =   2.076      CRITERION=   1.777   
    EQUIL ITER   4 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.9206E-07
     DISP CONVERGENCE VALUE   =  0.9206E-07  CRITERION=  0.1824E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC =  0.9206E-07
     FORCE CONVERGENCE VALUE  =  0.1425      CRITERION=   1.813     <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   4
 *** LOAD STEP     1   SUBSTEP     6  COMPLETED.    CUM ITER =     19
 *** TIME =  0.600000         TIME INC =  0.100000   

     FORCE CONVERGENCE VALUE  =   542.0      CRITERION=   1.957   
     DISP CONVERGENCE VALUE   =  0.1316E-03  CRITERION=  0.2213E-04
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.1316E-03
     DISP CONVERGENCE VALUE   =  0.7702E-04  CRITERION=  0.2258E-04
     LINE SEARCH PARAMETER =  0.5853     SCALED MAX DOF INC = -0.7702E-04
     FORCE CONVERGENCE VALUE  =   242.1      CRITERION=   2.069   
    EQUIL ITER   2 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.3019E-04
     DISP CONVERGENCE VALUE   =  0.3019E-04  CRITERION=  0.2304E-04
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC =  0.3019E-04
     FORCE CONVERGENCE VALUE  =   6.139      CRITERION=   2.078   
    EQUIL ITER   3 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.8557E-04
     DISP CONVERGENCE VALUE   =  0.8423E-04  CRITERION=  0.2352E-04
     LINE SEARCH PARAMETER =  0.9843     SCALED MAX DOF INC =  0.8423E-04
     FORCE CONVERGENCE VALUE  =   23.51      CRITERION=   2.039   
    EQUIL ITER   4 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.2336E-06
     DISP CONVERGENCE VALUE   =  0.2080E-06  CRITERION=  0.2399E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =  0.8904     SCALED MAX DOF INC =  0.2080E-06
     FORCE CONVERGENCE VALUE  =   2.637      CRITERION=   2.081   
    EQUIL ITER   5 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2534E-04
     DISP CONVERGENCE VALUE   =  0.2534E-04  CRITERION=  0.2448E-04
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC = -0.2534E-04
     FORCE CONVERGENCE VALUE  =   3.534      CRITERION=   2.149   
    EQUIL ITER   6 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.6476E-05
     DISP CONVERGENCE VALUE   =  0.3299E-05  CRITERION=  0.2498E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =  0.5094     SCALED MAX DOF INC = -0.3299E-05
     FORCE CONVERGENCE VALUE  =   2.087      CRITERION=   2.196     <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   6
 *** LOAD STEP     1   SUBSTEP     7  COMPLETED.    CUM ITER =     25
 *** TIME =  0.700000         TIME INC =  0.100000   

     FORCE CONVERGENCE VALUE  =   512.7      CRITERION=   2.219   
     DISP CONVERGENCE VALUE   =  0.1489E-03  CRITERION=  0.2213E-04
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.1489E-03
     DISP CONVERGENCE VALUE   =  0.9240E-04  CRITERION=  0.2258E-04
     LINE SEARCH PARAMETER =  0.6205     SCALED MAX DOF INC =  0.9240E-04
     FORCE CONVERGENCE VALUE  =   219.3      CRITERION=   2.177   
    EQUIL ITER   2 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2054E-04
     DISP CONVERGENCE VALUE   =  0.2054E-04  CRITERION=  0.2304E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC = -0.2054E-04
     FORCE CONVERGENCE VALUE  =   3.308      CRITERION=   2.246   
    EQUIL ITER   3 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.7983E-04
     DISP CONVERGENCE VALUE   =  0.7983E-04  CRITERION=  0.2352E-04
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC = -0.7983E-04
     FORCE CONVERGENCE VALUE  =   21.38      CRITERION=   2.369   
    EQUIL ITER   4 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.6662E-05
     DISP CONVERGENCE VALUE   =  0.2147E-05  CRITERION=  0.2399E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =  0.3223     SCALED MAX DOF INC = -0.2147E-05
     FORCE CONVERGENCE VALUE  =   14.64      CRITERION=   2.420   
    EQUIL ITER   5 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=  0.3045E-05
     DISP CONVERGENCE VALUE   =  0.3045E-05  CRITERION=  0.2448E-04 <<< CONVERGED
     LINE SEARCH PARAMETER =   1.000     SCALED MAX DOF INC =  0.3045E-05
     FORCE CONVERGENCE VALUE  =  0.5574      CRITERION=   2.466     <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   5
 *** LOAD STEP     1   SUBSTEP     8  COMPLETED.    CUM ITER =     30
 *** TIME =  0.800000         TIME INC =  0.100000   

     FORCE CONVERGENCE VALUE  =   554.7      CRITERION=   2.513   
     DISP CONVERGENCE VALUE   =  0.2173E-02  CRITERION=  0.1262E-03
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC= -0.2173E-02
     DISP CONVERGENCE VALUE   =  0.1796E-02  CRITERION=  0.1288E-03
     LINE SEARCH PARAMETER =  0.8266     SCALED MAX DOF INC = -0.1796E-02
     FORCE CONVERGENCE VALUE  =   4696.      CRITERION=   3.808   

 *** ERROR ***                           CP =     385.198   TIME= 10:07:27
 The value of UX at node 6671 is 6580940.52.  It is greater than the    
 current limit of 1000000 (which can be reset on the NCNV command).     
 This generally indicates rigid body motion as a result of an           
 unconstrained model.  Verify that your model is properly constrained.  

 

Изменено 15 октября 2021 пользователем AlexKaz

[forums.sis 46]

Если не с ложно, можно просто геометрию. А то у меня до сих пор 2019 стоит, чет не обновлю все никак.

[AlexKaz 1 821]

@forums.sis

Экспорт из студенческого WB, возможно откроется в Design Modeler.

geometry.rar

Geom.rar

Изменено 15 октября 2021 пользователем AlexKaz

[AlexKaz 1 821]

29 минут назад, AlexKaz сказал:

Не, между 128 и 144.

Сходимость при задании перемещений отличная, счётная погрешность низкая, и найденная реакция вполне себе попадает в диапазон выше

[rasta89 113]

12 часов назад, AlexKaz сказал:

Сходимость при задании перемещений отличная, счётная погрешность низкая, и найденная реакция вполне себе попадает в диапазон выше

Спасибо что посмотрели.

Это был один сегмент, в цанге их 22, к примеру целая цанга срывается на 3000N, проверочный расчет еще не делал, считает 5 часов полную модель. Скорее всего при проверке тоже будет примерно 3100-3200N.

Какая нагрузка все же ближе к истине, так что была при перемещении найдена, или на проверочном этапе?