Geotechnical boundary conditions, initial conditions, and loading for Abaqus
Expert guidance on setting up boundary conditions, initial conditions, and loading sequences for geotechnical Abaqus simulations. Covers geostatic stress initialization, groundwater conditions, staged construction loading, and infinite boundaries.
岩土 Abaqus 仿真中的边界条件、初始条件和加载序列设置指导。覆盖地应力初始化、地下水条件、分步施工加载和无限边界。
*STEP, NAME=Geostatic
*GEOSTATIC
*DLOAD
ALL_SOIL, GRAV, 9.81, 0., -1., 0.
*END STEP
*INITIAL CONDITIONS, TYPE=STRESS, GEOSTATIC
SOIL_SET, -100., 10., -30., 0., 0.8
| 土体类型 | K₀ 范围 | 推荐值 |
|---|---|---|
| 正常固结砂土 | 0.35-0.50 | 1 - sin φ |
| 正常固结黏土 | 0.50-0.70 | 1 - sin φ |
| 超固结黏土 | 0.50-2.50 | K₀_nc × OCR^(sin φ) |
| 软岩 | 0.30-0.50 | 0.40 |
| 硬岩 | 0.10-0.50 | 按实测 |
*STEP, NAME=VerifyGeostatic
*GEOSTATIC
*OUTPUT, FIELD
*ELEMENT OUTPUT
S, E, COORD
*NODE OUTPUT
U, RF
*END STEP
*INITIAL CONDITIONS, TYPE=PORE PRESSURE
SOIL_SET, 0., 10., 100., 0.
*BOUNDARY
BOTTOM_NODES, 8, 8, 100.
SURFACE_NODES, 8, 8, 0.
*BOUNDARY
WELL_NODES, 8, 8, -50.
*BOUNDARY
BOTTOM, 1, 3, 0. ** 底面固定 (x,y,z)
LEFT, 1, 1, 0. ** 左侧法向约束 (x)
RIGHT, 1, 1, 0. ** 右侧法向约束 (x)
FRONT, 2, 2, 0. ** 前面法向约束 (y) - 3D
BACK, 2, 2, 0. ** 后面法向约束 (y) - 3D
** 顶面自由（地表）
*BOUNDARY
SYMM_PLANE, YSYMM ** y=0 对称面：u2=ur1=ur3=0
*ELEMENT, TYPE=CIN3D8, ELSET=INFINITE
...
*DLOAD
ALL_ELEMENTS, GRAV, 9.81, 0., -1., 0.
*DSLOAD
SURFACE_TOP, P, 100. ** 均布压力 100 kPa
*CLOAD
NODE_SET, 2, -50. ** 竖向集中力
*STEP, NAME=Layer1
*STATIC
*MODEL CHANGE, ADD
FILL_LAYER1
*DLOAD
FILL_LAYER1, GRAV, 9.81, 0., -1., 0.
*END STEP
*STEP, NAME=Layer2
*STATIC
*MODEL CHANGE, ADD
FILL_LAYER2
*DLOAD
FILL_LAYER2, GRAV, 9.81, 0., -1., 0.
*END STEP
*STEP, NAME=Earthquake
*DYNAMIC, IMPLICIT
0.01, 10., 1e-6, 0.05
*BOUNDARY, TYPE=ACCELERATION
BASE_NODES, 1, 1, 1.0
*AMPLITUDE, NAME=EQ_MOTION, INPUT=earthquake.txt
*BOUNDARY, AMPLITUDE=EQ_MOTION
BASE_NODES, 1, 1, 1.0
*END STEP
Step 1: GEOSTATIC（地应力平衡）
Step 2: 开挖/施工（MODEL CHANGE）
Step 3: 固结（SOILS, CONSOLIDATION）
Step 4: 长期运营/地震（视需求）
*STEP, NAME=Consolidation
*SOILS, CONSOLIDATION, END=PERIOD
0.001, 100., 1e-8, 10.
*END STEP
*STEP, NLGEOM=YES
| 问题 | 原因 | 解决方案 |
|---|---|---|
| 地应力不平衡 | K₀ 与弹性参数不匹配 | K₀ = ν/(1-ν) 或使用 *GEOSTATIC |
| 孔压初始不对 | 水位线设置错误 | 检查 z 坐标方向和参考水位 |
| 无限元不起作用 | 极点方向反了 | 检查节点编号顺序 |
| 分步施工应力跳变 | 新加单元无初始应力 | 在 MODEL CHANGE 后加 *GEOSTATIC 子步 |
| 固结步发散 | 初始增量太大 | 减小初始增量至 0.001 或更小 |
| 问题类型 | 水平范围 | 竖向深度 |
|---|---|---|
| 浅基础 | ≥ 5B (B=基础宽) | ≥ 3B |
| 深基坑 | ≥ 3~5 倍开挖深度 | ≥ 2 倍开挖深度 |
| 隧道 | ≥ 3~5D (D=洞径) | 洞底下 ≥ 3D |
| 边坡 | 坡脚外 ≥ 2H (H=坡高) | ≥ 1.5H |
| 桩基 | ≥ 20D (D=桩径) | 桩底下 ≥ 5D |
abaqusgpt ask "深基坑开挖边界条件怎么设置" --domain geotechnical
abaqusgpt ask "如何做地应力平衡" --domain geotechnical
from abaqusgpt.skills.geotech_boundary import GeotechBoundarySkill
skill = GeotechBoundarySkill()
result = skill.execute({
"query_type": "boundary_setup",
"problem_type": "excavation",
"dimensions": "3D"
})