Какие ед.изм. КИН при PRCINT,PLCINT ? Скриншот внутри

[4estor 1]

Здравствуйте. В ANSYS APDL моделирую трещину, вывожу К1 через команды PRCINT и PLCINT, значения отличаются в 10^6 раз!

В каких единицах измерения какая команда выводит? Понимаю, что это зависит от единиц измерения давления и длины в моей задаче, но какова общая формула? Например типа давление*корень(длина).

Вот на моем скриншоте PLCINT вывел распределение КИН вдоль трещины (для 6го контура) порядка 1000 в среднем, а PRCINT порядка 0,001.

[Борман 2 375]

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

вдоль трещины

Какие у вас приняты единицы длины в модели ? Метры ? 

 

[4estor 1]

В 24.02.2018 в 08:37, Борман сказал:

В 23.02.2018 в 15:11, 4estor сказал:

вдоль трещины

Какие у вас приняты единицы длины в модели ? Метры ? 

миллиметры.

[soklakov 1 475]

@Борман , таки чего? В метры переходить?

[Борман 2 375]

@soklakov

ХЗ. Вопрос, я так понял, в том, что диаграмма и листинг показывают разные числа.

Я не решал таких задач.

 

 

[soklakov 1 475]

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

Вопрос, я так понял, в том, что диаграмма и листинг показывают разные числа.

да, пожалуй) перейдя в метры, можно узнать, сохранится ли сия странность. если да, то интересно. если нет - еще интереснее))

[4estor 1]

Я перевел модель в метры и в Па (геометрию, свойства, нагрузки скейлил) и вот что получилось:

 

 

Различие в 10^11 раз...

[Nikon6260 20]

@4estor 

Есть возможность посмотреть на расчетную базу?

И в какой версии ANSYS считаете?

[4estor 1]

@Nikon6260 что такое расчетная база? я не против показать, только что именно?)

14.0 

[piden 2 874]

49 minutes ago, 4estor said:

что такое расчетная база? я не против показать, только что именно?)

Выполните команды

/prep7
cdwrite,db,model,cdb

и потом выложите файл model.cdb, поместив в архив.

[4estor 1]

@piden ок, залил cdb https://yadi.sk/d/W2crofd93SqSz3

Хотелось бы узнать, в каких вообще единицах измерения листинг и диаграммы показывают, независимо от того, в чем у меня длины/силы. Ну то есть общий вид типа давление*корень из длины.

[Fedor 1 597]

https://en.wikipedia.org/wiki/Fracture_mechanics   Может метод Паскаля - заменяй определяемое определением - поможет ... :) 

[Fedor 1 597]

Там же корень из длины в    Irwin's modifications, отсюда и 10**6 видимо ... :)  

[Борман 2 375]

В 23.02.2018 в 16:11, 4estor сказал:

ANSYS APDL моделирую трещину

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

[Fedor 1 597]

Цитата

 

https://en.wikipedia.org/wiki/Stress_intensity_factor   тут еще посмотреть

 

17.7. POST1 - Crack Analysis  тут в хелпах много написано. Скорее всего единицы как обычно должны быть согласованы :)

Изменено 28 февраля 2018 пользователем Fedor

[Nikon6260 20]

1. Ваша модель выполнена в ANSYS 14.0. и похоже экспортирована из Workbench. Не используйте версию 14.0, она кривая. Стабильные версии 14.5, 17.2 и выше. Посчитал в 18.0 все сходиться (скрин прилагаю).

2. Для корректного расчета КИНов во фронте трещины, вокруг фронта строится тор с регулярной сеткой. Количество слоев 5-6 (это будут контура интегрирования). Представьте где находится 6-й контур для Вашей текущей сетки.

3. Не понял чего Вам 9 одинаковых расчетных шагов?

4. Также прикрепляю небольшой макрос для считывания и экспорта в текстовый файл КИНов.

SIF_OUT.mac

[Борман 2 375]

Только что, Борман сказал:

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

А бл.. и правда трещина.

[Борман 2 375]

Оценка на пальцах показала число порядка 43МПа*м^0.5 в самой глубокой точке. Была бы мельче сетка - было бы точнее.

[4estor 1]

@Fedor Прочитал, но я и так знаю определение КИН. Он пропорционален Давлению*корень(длина) или в системе СИ: Па*корень(м). Главное, по каким формулам ее вычисляет Ансис? Нашел только для J-интеграла. И неужели, для диаграммы и листинга эти формулы разные...

Когда перешел от МПа&мм в Па&м, то диаграммы увеличились в 3*10^4 раз - так и должно быть, т.к. 1 МПа*корень(мм) = 3,16*10^4 Па*корень(м). А значения в листингах наоборот уменьшились в 30 раз. 

 

@Борман ну это по порядку величины совпало с диаграммами.

 

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

9 подшагов не одинаковые, там давление линейно возрастает от 0 до 9 МПа (kbc,0).

Насчет тора с рег.сеткой (сингулярными элементами?) - знаю, но как быть, если 1 поверхность плоская (плоскость симметрии в моем случае), а вторая (в моем случае это впадина) не плоская, ведь тор подразумевает плоские торцы? Или будет немного торчать из впадины? Да, представляю, где мой 6й контур, возможно просто сгустить сетку помогло бы?

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

Спасибо за макрос. Удобней оцифровывать, чем стандартный листинг.

К листингу теперь доверия меньше... а ведь раньше всегда из него брал результаты K1 и J, значит возможно многие прошлые мои работы были ошибочными, стремно осознавать это.

[soklakov 1 475]

17 минут назад, 4estor сказал:

значит возможно многие прошлые мои работы были ошибочными, стремно осознавать это.

лайк  в каком-то конкретном месте поломку ждать?