[Kupec]

Сделал модель основание+фундамент.
Использовал для грунта элементы Solid 45, для плиты shell 63.
Грунт задал по модели Друкера-Прагера.
В результате получилась ерунда, после 2,5 часового расчета (машина не сильная).

Может быть дело в типе элемента?
Граничные условия: симметрия по торцам solid 45 и заделка внизу по плоскости.

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

Заранее благодарю.

[Фролов Дмитрий]

Сообщение от MrWhite Потому как в той статье, что приводил Дмитрий Фролов (настройки DP и EDP) многое неправильно.

это которая Имаоки статья?

[MrWhite]

Сообщение от Фролов Дмитрий это которая Имаоки статья?

Ага. Ни в коем случае Вам ни в упрек, статья хорошая только грешит принципиальными неточностями, из-за которых я очень много времени убил

[Фролов Дмитрий]

Сообщение от MrWhite статья хорошая только грешит принципиальными неточностями

может тогда написать ему, а то люди очень часто на эту статью выходят

[MrWhite]

Сообщение от Фролов Дмитрий EDP: C=1000 Fi=35 Del=35 Alpha1=(6*SIN(Fi))/(3-SIN(Fi)) Alpha2=(6*SIN(Del))/(3-SIN(Del)) Sigma=(6*C*COS(Fi))/(SQRT(3)*(3-SIN(Fi))) TB,EDP,matid,,,LYFUN TBDATA,1,Alpha1,Sigma TB,EDP,matid,,,LFPOT TBDATA,1,Alpha2

Ну вот вы по статье делали определение материала. Красным я отметил что не верно.
Наклон поверхности потенциала текучести определяется НЕ только углом дилатации. По сути при угле дилатации = 0, должно быть Alpha1 = Alpha2.
Sigma определяется исходя из того что поверхность DP вписана в поверхность Кулона-Мора, что в тексте хелпа и написано. Но на илюстрации хелпа и по сравнительным расчетам в материала DP Ansys поверхность описана вокруг Кулона-Мора.
И что-то еще по приведению шатровой модели EDP, но сейчас не на работе - не помню.
Как видите даже по русски путано выглядит

[Фролов Дмитрий]

да уж...а сильно ли отличаются правильные формулы? ну или подождем вашей статьи

[eraserheadd]

Сообщение от MrWhite Сообщение от Фролов Дмитрий EDP: C=1000 Fi=35 Del=35 Alpha1=(6*SIN(Fi))/(3-SIN(Fi)) Alpha2=(6*SIN(Del))/(3-SIN(Del)) Sigma=(6*C*COS(Fi))/(SQRT(3)*(3-SIN(Fi))) TB,EDP,matid,,,LYFUN TBDATA,1,Alpha1,Sigma TB,EDP,matid,,,LFPOT TBDATA,1,Alpha2

Сообщение от MrWhite Ну вот вы по статье делали определение материала. Красным я отметил что не верно. Наклон поверхности потенциала текучести определяется НЕ только углом дилатации. По сути при угле дилатации = 0, должно быть Alpha1 = Alpha2. Sigma определяется исходя из того что поверхность DP вписана в поверхность Кулона-Мора, что в тексте хелпа и написано. Но на илюстрации хелпа и по сравнительным расчетам в материала DP Ansys поверхность описана вокруг Кулона-Мора. И что-то еще по приведению шатровой модели EDP, но сейчас не на работе - не помню. Как видите даже по русски путано выглядит

В DP поверхность задается в пространстве √(J2 )-I1
В справке ANSYS написаны формулы для нахождения коэффициентов уравнения поверхности текучести:
β=(2*SIN(Fi))/(SQRT(3)*(3-SIN(Fi)))

σy=(6*C*COS(Fi))/(SQRT(3)*(3-SIN(Fi)))

Введенные в модель с и ϕ используется в этих формулах. С помощью них строится поверхность текучести DP, которая описывает поверхность Мора-Кулона с такими же значениями с и ϕ.
Для того, чтобы построить вписанную поверхность DP, необходимо использовать формулы :
β=(2*SIN(Fi))/(SQRT(3)*(3+SIN(Fi)))

σy=(6*C*COS(Fi))/(SQRT(3)*(3+SIN(Fi)))
Ну и, как правильно заметил MrWhite, с и ϕ тогда нужно задавать иначе.

В справке ANSYS написано, что закон пластического течения ассоциированный в том случае, если ϕ= ϕf (угол дилатансии) Тогда и будет Alpha1 = Alpha2. А пластический потенциал имеет такое же уравнение, как и поверхность текучести, только вместо ϕ - ϕf. При ϕf=0 течение неассоциированное будет, а не наоборот.

В EDP Linear, Power low и Hyperbolic поверхность задается в пространстве √(3J2 )-I1/3
Поэтому и формулы для коэффициентов уравнения отличаются:
α=(3*2*SIN(Fi))/(3-SIN(Fi))

σy=(6*C*COS(Fi))/(3-SIN(Fi))

По-моему, в EDP закон пластического течения можно задавать отличным от поверхности текучести. Например, поверхность текучести Linear, а пластический потенциал от Hyperbolic. Но в целом там тоже самое - уравнение поверхности текучести Linear соответствует пластическому потенциалу Linear и т.п.

В EDP Cap поверхность задается в пространстве √(J2 )-I1. Т.е. параметры α и σс в EDP Cap по сути те же, что и β с σy в DP.
НО у Имаоки верно сказано, что связь c параметрами других моделей Друкера-Прагера сохраняется до тех пор, пока идет линейное представление поверхности текучести. Если не учитывать шатры, то соответствие будет при параметре шатровой модели A=0.