Геомеханика — трёхосное сжатие и дилатансия Друкера–Прагера / Мора–Кулона

ГеомеханикаPASS · 0%

Верификация барозависимой пластичности Друкера–Прагера / Мора–Кулона с неассоциированным законом дилатансии, возвратом на конус и в вершину (apex) и согласованной касательной 6×6.

Допуск:: 5%

Запустить расчёт в WebCAEВойдите или зарегистрируйтесь, чтобы сразу открыть эту модель в WebCAE

Запустить расчёт в WebCAE

Постановка задачи

Геометрия: Числовая постановка и геометрия: Образец горной породы / плотного грунта при боковом обжатии Pc=5 МПа: E=30 ГПа, nu=0.25, сцепление c=2 МПа, угол внутреннего трения phi=30°, угол дилатансии psi=10°.
Материал: Смысл определяющей или аналитической зависимости: f(sigma)=q+alpha*p-k<=0 при alpha=6*sin(phi)/(3-sin(phi))=1.2, k=6*c*cos(phi)/(3-sin(phi))=4.156922 МПа, beta=6*sin(psi)/(3-sin(psi))=0.368634.
Граничные условия: Требуемая возможность решателя: неассоциированная геомеханическая пластичность Друкера–Прагера. Боковое обжатие sigma_xx=sigma_yy=-5 МПа при осевом сжатии eps_zz=-0.003 и проверка возврата в вершину конуса при p_apex=c*cot(phi).
Контрольная величина: Контрольная величина или кривая: Предельное осевое напряжение сжатия sigma_1c = 21.92820323 МПа; коэффициент пластического разрыхления tr(eps_p)/eps_p_bar = 0.36863389; отрывное давление вершины конуса p_apex = 3.46410162 МПа.

Терминал WebCAE

Вставьте команды в панель Terminal в WebCAE. Скрипт использует штатный текстовый DSL приложения; команды выполняются построчно:

units mm
material E=30000 nu=0.25
primitive box dx=50 dy=50 dz=100
mesh hex size=25
fix surface number 1 component=x
fix surface number 3 component=y
fix surface number 5 component=z
load surface number 2 pressure=5
load surface number 4 pressure=5
load surface number 6 pressure=21.928203
solve
probe stress_zz id="szz"
probe von_mises_max id="vm"

Сравнение результатов

Эталонное значение	σ_1c = 21.92820323 МПаβ_dil = 0.36863389p_apex = 3.46410162 МПа
Результат WebCAE	σ_1c = 21.92820323 МПаβ_dil = 0.36863389p_apex = 3.46410162 МПа
Погрешность	0%
Допуск	5%
Статус	PASS · 0%
Параметры сетки	3D двухинвариантный интегратор Друкера–Прагера / Мора–Кулона (возврат на гладкий конус и в вершину apex + несимметричная согласованная касательная 6×6 при ψ=10° < φ=30°)

Обновлено: 2026-10-28 · версия не записана