GTS NX主要分析功能介绍
合集下载
01_MIDAS GTS NX分析手册简介

Table 1.3.2 Files created during numerical analysis
File name InputName.DASM#.bin InputName.FACT#.bin#
InputName.EIGS#.bin#
InputName.MSTO#.bin
Point of creation/Content Generation, finite element related information for all analyses Generation, matrix information when selecting the multi frontal method Lancoz resampling information when selecting eigenvalue analysis Large scale matrix, vector related information internally recorded in the disk
Section 1. Overview | 1
ANALYSIS REFERENCE
Chapter 1. Introduction
This manual is mainly composed of theoretical and technical information that make up the base of detailed analysis for the effective usage of the GTS NX program. The contents of each chapter are as follows:
The temporary files created during the analysis process of GTS NX and their content are as follows.
midas gts NX分析工况

利用施工阶段分析可以模拟岩土的施工过程。施工阶段分析由多个施工阶段构成,可以 按各阶段激活或钝化荷载、边界条件或单元,这种荷载、边界或单元的变化适用于任一阶 段。在GTS NX中,可以使用如下多种分析功能进行施工阶段分析。 应力-边坡分析 施工阶段过程中的应力分析及边坡稳定分析。 渗流分析 按施工阶段的稳定流及瞬态流分析。 应力-渗流-边坡耦合分析 施工过程中进行渗流-应力耦合及边坡稳定分析。 固结分析 对施工阶段中堆土及环境变化的固结分析。 完全应力-渗流耦合分析 考虑非稳定渗流的完全应力- 渗流耦合分析。
file://C:\Users\002702\AppData\Local\Temp\~hhB27D.htm
2015/2/13
w
页码,2/33(W)
新建
概要 创建执行分析的“分析工况”。设置各分析方法的使用的分析条件(网格组、边界条 件、荷载条件等)的阶段。特别是施工阶段分析的情况,可以采用5种不同的方法分析并且可 设定分析要用的数据。而且,可调整详细的分析选项及输出结果选项,并且可通过设置多个 施工阶段组对一个模型进行反复分析。
w
页码,1/33(W)
பைடு நூலகம்
分析工况
岩土分析可以通过与一般的结构分析的比较来解释。结构分析强调在结构上起作用的不 确定性荷载的权重。因此,可对通过系统性地组合各种结果获取的最大构件力执行构件设 计。与此相反,在岩土分析中,与荷载相比更重视的是施工阶段及材料自身的不确定性,掌 握岩土内部的物理性状态是非常重要的。因此,在岩土分析中,建模过程中会采用实体单 元,使之最大限度地反映岩土的形状和施工状况。应尽可能地考虑材料的各种非线性、各向 异性及原场地应力状态,来反映真实的现场状态。 岩土分析的程序可用于模拟实际现场条件,判断设计或施工条件是否可行。在岩土分析 中,涵盖的分析领域从一般性的静力分析,到渗流分析、应力-渗流耦合分析、固结分析、施 工阶段分析、动力分析、边坡稳定分析等。 提供的岩土分析功能如下。在这部分简要地概述了分析方法和对分析选项进行说明。详 细的分析信息须参考理论分析手册第五章。 1. 静力分析(Static Analysis) (1)线性静力分析 (2)非线性静力分析(非线弹性或弹-塑性分析) 2. 3. 施工阶段分析(Construction Stage Analysis) 渗流分析(Seepage Analysis) (1)稳定流分析(Steady State) (2)瞬态流分析(Transient) 4. 应力-渗流耦合分析(Coupled Seepage-Stress Analysis) (1)渗流-应力连续分析(Seepage-Stress Sequential Analysis) (2)固结分析(Consolidation Analysis) (3)完全应力-渗流耦合分析(Fully-coupled Seepage-Stress Analysis) 5. 动力分析(Dynamic Analysis) (1)特征值分析(EigenValue Analysis) (2)反应谱分析(Response Spectrum Analysis) (3)线性时程分析(振型叠加法)(Linear Time History(Modal)) (4)线性时程分析(直接积分法)(Linear Time History(Direct)) (5)非线性时程分析(Nonlinear Time History Analysis) (6)二维等效线性分析(2D Equivalent Linear Analysis) 6. 边坡稳定分析(Slope Stability Analysis) (1)边坡稳定分析(SRM)(Strength Reduction Method) (2)边坡稳定分析(SAM)(Stress Analysis Method) (3)非线性时程分析 + SRM(Dynamic-Slope Coupled Analysis)
GTSNX建模技巧及特色功能

性进行判断,对于多个简单模型的结果比较,会更容易理解和解决问题。
模型仅是辅助用资料
对于缺少参考资料的问题,我们建模的目的并不是为了获得绝对值,而是为了增加对该 类型问题的理解,为了建立能够让工程师做出科学判断的参考资料而已。
要有耐心
建立一个简单的模型对结果进行预测。这些简单模型的结果将会为你提供获得资料的新的方式
开挖面对隧洞围岩的虚拟支撑作用发生亍距开挖面 ≤(1.5~3.0)倍洞跨范围内,在距开挖面最多3倍洞径之后,可认为开挖面
的空间效应已完全消失*
*孙钓,朱合华.软弱围岩隧洞斲工性态的力学模拟不分析,1994,岩土力学,15(4)
取3 倍洞距边界范围对隧道进行数值模拟和结构设计,其分析结果和结构设计斱案最为安全*
*苏晓堃.隧道开挖数值模拟的围岩边界取值范围研究[J].铁道工程学报,2012,(3):64-68.DOI:10.3969/j.issn.1006-2106.2012.03.014.
导入几何
导入CAD
.
导入 CIVIL GEN
.
导入DWG
导入DXF
CAD图位亍原点附近,以m为单位,真实相交封闭
几何组分组—图层 自劢交叉分割
在相同的节点自由度情冴下,矩形单元的计算精度要比三角形 单元高
在很多模型里面,都会涉及到梁单元不实体单元的耦合问题,如:建筑物桩基、 基坑的排桩支护、边坡中的抗滑桩等。
• 网格生成器不要更换
默认四面体 混合网格
• 确认属性赋予正确
• 先划分关心区域、尺寸 较小区域
• 随时检查自由面
• 导出网格
网格-工具-表格-节点/单元表格
• 收容差 • 推荐控制值为0.001
边坡稳定性分析
模型仅是辅助用资料
对于缺少参考资料的问题,我们建模的目的并不是为了获得绝对值,而是为了增加对该 类型问题的理解,为了建立能够让工程师做出科学判断的参考资料而已。
要有耐心
建立一个简单的模型对结果进行预测。这些简单模型的结果将会为你提供获得资料的新的方式
开挖面对隧洞围岩的虚拟支撑作用发生亍距开挖面 ≤(1.5~3.0)倍洞跨范围内,在距开挖面最多3倍洞径之后,可认为开挖面
的空间效应已完全消失*
*孙钓,朱合华.软弱围岩隧洞斲工性态的力学模拟不分析,1994,岩土力学,15(4)
取3 倍洞距边界范围对隧道进行数值模拟和结构设计,其分析结果和结构设计斱案最为安全*
*苏晓堃.隧道开挖数值模拟的围岩边界取值范围研究[J].铁道工程学报,2012,(3):64-68.DOI:10.3969/j.issn.1006-2106.2012.03.014.
导入几何
导入CAD
.
导入 CIVIL GEN
.
导入DWG
导入DXF
CAD图位亍原点附近,以m为单位,真实相交封闭
几何组分组—图层 自劢交叉分割
在相同的节点自由度情冴下,矩形单元的计算精度要比三角形 单元高
在很多模型里面,都会涉及到梁单元不实体单元的耦合问题,如:建筑物桩基、 基坑的排桩支护、边坡中的抗滑桩等。
• 网格生成器不要更换
默认四面体 混合网格
• 确认属性赋予正确
• 先划分关心区域、尺寸 较小区域
• 随时检查自由面
• 导出网格
网格-工具-表格-节点/单元表格
• 收容差 • 推荐控制值为0.001
边坡稳定性分析
GTS NX R 主要升级内容

ECS z
ECS y
ECS x
Qz Mz
Mz
Qz Mx
N xx B
My
Qy
N xx Mx
A
Qy My
ECS y
F
C
Stress recovery point (I-section)
ECS z
E
D
[梁单元单元坐标系、自由度、内力输出]
[梁单元考虑节点耦合]
2
Mother
element 2
1
Mother element 1
参考值(kN, m)
Cc / 2.303 / (1 + e) Cs / 2.303 / (1 + e)
(可近似取Cc / 5 ) 可近似取Cc / 20 莫尔-库伦模型中的破坏参数 莫尔-库伦模型中的破坏参数
0
1-sinφ (< 1)
当输入两个值时, 优先使用Pc 根据KNC得到(自动)
11 / 31
库伦、德鲁克-普拉格、胡克-布朗) 网格> 属性/坐标系/函数> 函数> 徐变 / 收缩函数 用户可以根据嵌入的17个各国设计规范定义徐变/收缩凼数。 网格> 属性/坐标系/函数> 函数> 弹性模量函数 用户可以根据嵌入的12个各国设计规范定义时间相关的弹性模量凼数。
[分析控制/ 时间步长]
[分析控制/混凝土龄期]
Porosity < Porosity(Max)
根据KNC得到 (自动计算) 根据 Eoedref 得到(自动)
9 / 31
GTSNX 2014 主要升级内容
2. 分析
2.1 硬化土模型(示例 : 硬化行为验证)
参数
midas GTS NX的线性和非线性动力分析

Maximum Accel.(m2/sec) 0 10 GTS NX Flush
[ 输入地层的动力非线性特性 ]
10
20 30
20
地表面加速度(m/sec2)
1.50 1.00
30 40 Depth (m)
40
50 60 70 80 90 GTS NX Flush
0.50 0.00 0 -0.50 -1.00
各土层最大加速度 – 用于判断液化
11
01 GTS NX 的动力分析概要 02 自由场分析 (FFA) 03 反应谱分析 04 二维等效线性分析
05 线性时程分析(直接积分/振型叠加)
06 非线性时程分析(直接积分
法) 07 非线性时程分析+强度折减法
08 振动加速度级的输出
反应谱分析
反应谱分析概要
0.03
Relative Displacement T1(m)
0.02 0.01 0.00
反应谱分析
-0.01
-0.02 -0.03 0 5
Time(sec) Linear Non-linear
10
15
7
01 GTS NX 的动力分析概要 02 自由场分析 (FFA) 03 反应谱分析 04 二维等效线性分析
→ 使用单自由度体系的最大位移、最大速度、最大加速度响应谱计算结构响应的方法
→ 虽然与时程分析方法相比结果有误差,可用于对分析效率有要求的大型结构或对结果精确度要求不高的结构
[ 无阻尼时]
[ 生成单自由度体系的位移响应谱的过程 ] [ 有阻尼时] 13
01 GTS NX 的动力分析概要 02 自由场分析 (FFA) 03 反应谱分析 04 二维等效线性分析
[ 输入地层的动力非线性特性 ]
10
20 30
20
地表面加速度(m/sec2)
1.50 1.00
30 40 Depth (m)
40
50 60 70 80 90 GTS NX Flush
0.50 0.00 0 -0.50 -1.00
各土层最大加速度 – 用于判断液化
11
01 GTS NX 的动力分析概要 02 自由场分析 (FFA) 03 反应谱分析 04 二维等效线性分析
05 线性时程分析(直接积分/振型叠加)
06 非线性时程分析(直接积分
法) 07 非线性时程分析+强度折减法
08 振动加速度级的输出
反应谱分析
反应谱分析概要
0.03
Relative Displacement T1(m)
0.02 0.01 0.00
反应谱分析
-0.01
-0.02 -0.03 0 5
Time(sec) Linear Non-linear
10
15
7
01 GTS NX 的动力分析概要 02 自由场分析 (FFA) 03 反应谱分析 04 二维等效线性分析
→ 使用单自由度体系的最大位移、最大速度、最大加速度响应谱计算结构响应的方法
→ 虽然与时程分析方法相比结果有误差,可用于对分析效率有要求的大型结构或对结果精确度要求不高的结构
[ 无阻尼时]
[ 生成单自由度体系的位移响应谱的过程 ] [ 有阻尼时] 13
01 GTS NX 的动力分析概要 02 自由场分析 (FFA) 03 反应谱分析 04 二维等效线性分析
GTSNX2015升级功能介绍

[工程实例: 桩筏基础]
3 / 31
GTSNX 2015 Enhancement
1. 前处理
1.2 动力分析>工具 > 人工地震波
通过嵌入的设计谱数据生成人工地震波. 以下规范设计谱数据可用.
修改PSD
读入设计谱数据
计算功率谱密度函数(PSD power Spectral Density )
计算加速度
Viscous boundary
[自由场作用(X), 吸收反射(O)]
6 / 31
GTSNX 2015 Enhancement
1. 前处理
1.3 单元> 自由场单元(动力分析无限单元)
二维场选择从自由边,三维场选择从自由面来定义自由场单元.
GTGSTNSNXX 22001155VR11.1升R级ele功as能e N介o绍te
Free field
Main domain
Free field
Seismic wave
[自由场单元示意图]
[自由场作用(O), 吸收反射(O)]
GTGSTNSNXX 22001155VR11.1升R级ele功as能e N介o绍te
[自由场作用(X), 吸收反射(X)]
Viscous boundary
增加功能
1. 前处理
1.1 荷载表格输入/输出 1.2 人工地震波生成器 1.3 自由场单元(动力分析无限单元) 1.4 非弹性铰
2. 分析
2.1 安全系数输出(莫尔-库伦准则) 2.2 材料: 范梅赛斯-非线性 2.3 材料 : 液化修正 UBCSAND 2.4 材料: 关口-太田(非粘性) 2.5 材料: 关口-太田(粘性) 2.6 材料: 广义霍克布朗 2.7 材料 : 2D 正交各向异性 (2D 结构单元) 2.8 材料: 完善修正莫尔库伦(硬化本构) 2.9 材料 : 兰贝格-奥斯古德滞回模型(Ramberg-Osgood) 2.10 材料: 哈丁-唐尼维奇滞回模型( Hardin-Drnevich) 2.11 选项 : 预估激活单元初始应力状态 2.12 选项 : 应力非线性时程分析
有关midas-GTS-NX软件的介绍

1 midas GTS背景介绍
midas GTS是北京迈达斯技术有限公司研发的岩土用软件,具有迈达斯软件专有优势:汉化界面、交互式操作、强大的可视化。
迈达斯技术有限公司由韩国浦项制铁发展而成的土木工程计算分析软件开发公司,具有独立的研发团队,并在中国、美国、日本、英国、印度、俄罗斯、新加坡等国家成立分公司。
北京迈达斯技术有限公司拥有一批国内研发团队,对软件再次“加工”,让midas 系列软件更加适合中国用户。
此外,北京迈达斯技术有限公司全面负责产品的销售与技术支持。
2 midas GTS功能介绍
midas GTS是为能够迅速完成对岩土及隧道结构的分析与设计而开发的“岩土隧道结构专用有限元分析软件”,是一款采用windows风格操作界面的完全中文化软件,能够提供完全的三维动态模拟功能。
程序提供应力分析、动力分析、渗流分析、应力-渗流耦合分析、边坡稳定分析、衬砌分析和设计功能,并提供莫尔库伦、修正莫尔库伦、邓肯-张、修正剑桥等14种本构及用户自定义本构模型;程序还提供便捷的几何建模功能、地形生成器、隧道建模助手、锚杆建模助手以及丰富的后处理结果;可以广泛应用于地下结构、岩土、水工、地质、矿山、隧道等方面的分析及科研。
自2005年在中国发布至今,在广大用户的信任和支持下,已经走过了九个年头,成为了岩土行业主流的分析与设计软件。
随着行业的发展,GTS也面临新技术的发展及完善,因此在2014年推出midas GTS NX(New Experience),分别在前处理、后处理及计算阶段进行了功能的改善,此外新增加了新的分析功能,具体见下表.
GTS NX功能一览表。
midasgtsnx隧道衬砌分析

Section 3
定义材料及特性
▶表. 结构材料
3.1 地基及结构材料的定义。
定义结构材料后,创建网格时,定义各单元上欲分配的特性。 结构构件使用的材料如下表。结构材料的情况下,模型类型勾选 [结构] 部分。
名称 材料 模型类型 弹性系数(E) 泊松比(v) 单位重量(r)
C27 各向同性 Elastic
• 在工作目录树 > 网格上选择墙体、拱顶、地板网格。
• 在右击鼠标 > 上下文菜单上,勾选显示 > 单元坐标系。
• 选择网格 > 单元 > 单元参数 > 一维表单。
• 在下拉菜单上选择 [修改坐标系]。
• 方向不同的各单元利用翻转方向(X-轴)设置使其拥有同一的方向。
Chapter 2. 隧道衬砌分析 | 7
是修改单元坐标系的过程。 单元坐标系是所有单元的构件力以及应力的输出标准,特别是规定梁单元的剪切刚性和抗弯刚性的输入方向的标准,所以
正确地理解这个概念非常重要。在和桁架或张力/压缩力受压单元一样,只拥有轴向刚性的单元的情况下,仅有单元坐标
系 x 轴有意义,虽然 y、z 轴没有意义但在图形屏幕上指定构件切面的配置方向时是有必要的。
(拱顶部 : 分割个数 20, 底板: 分割个数 20).
6 | Chapter 2. 隧道衬砌分析
Basic Tutorials
▶开挖隧道
Chapter 2. 隧道衬砌分析
▶创建网格(线)
▶▶创建的网格(节点)
*
: 网格 > 单元 > 参数 (Mesh > Element > Parameters)
Basic Tutorials
Chapter 2. 隧道衬砌分析
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
[应力-渗流完全耦合分析]
5
验证例题
→ 对比参考文献的验证 → 使用标准例题对单元/材料的验证
Griffiths, D.V and Lane, P.A., Slope stability analysis by finite elements, Geotechnique, 1999
→ 对新功能有超过100个以上的验证例题
06 支持动力非线性分析
07 支持应力-渗流完全耦合分析
64位高性能计算
→ 可使用电脑配置的所有内存:32位程序仅能使用4GB内存 → In Core计算技术:考虑用户操作环境自动调整数据存储方法 → 并行计算技术:支持多核多线程计算
Out of core
[Model export]
In core
Pivot (CPU)
Update stiffness and residual(K i , f i ) thread0 Update stiffness and residual(K i , f i ) thread2 Update stiffness and residual(K i , f i ) thread3
m
k p,n
前阶段的渗透系数
150
k p , n 1 后阶段的渗透系数 m 亚松弛系数
[孔隙水压结果]
Strack, O. D. L. and Asgian, M. I., "A New Function for Use in the Hodograph Method“, Water Resources Research, 1978
x 3.08
191 sec
CPU Memory
[模型信息]
Intel Xeon X5647 2.93 GHz 8 core DDR3 48 Gigabyte
300 sec 200 sec 100 sec 0 sec
GTS
GTS NX P1
GTS NX P4
单元 节点
504,290 个 79,574 个 299,757 个
实现混合网格的自动划分
→ 支持焊接接触单元 实现不同网格尺寸区域的自动耦合]
[定义三维水面进行渗流分析]
[定义降雨条件进行渗流分析] [使用了D-P本构进行边坡稳定分析]
4
强化的计算内核
→ 非线性荷载的二分法 用户仅需输入大致的荷载增量,其余由程序自动调整 → 不变的切线刚度和亚松弛法 可提高应力/渗流分析的收敛性 → 使用荷载分割和弧长法进行边坡稳定分析 可快速获得精确的安全系数
→ 部分饱和岩土的渗透系数 k kr p ksat kr : 투수계수비 함수 渗流系数比函数
渗透系数比由孔隙压力决定 提供Gardner、Frontal和Van Genuchten函数 → 使用体积含水率计算部分饱和岩土的饱和度
ij ij se pw ij
'
分析功能的改善
基于64位高性能计算技术的分析功能的创新
→ 全新的程序框架 → 更多的分析控制选项 → 强化的计算内核 → 更多的分析方法
基于64位的 全新的程序 框架
下一代 求解器
GTS NX
MultiCore
EM64T
GPU
Coprocessor
3
全新的程序框架
→ 全方位支持64位 适应电脑硬件条件的发展 → 支持多核多线程/GPU计算 最大程度的使用硬件资源 → 支持金字塔形五面体单元
Thread
Thread
Nonlinear load increment
In Core 或 Out of Core
Thread
Thread
10
非线性施工阶段分析例题
→ 基于64 bit IC技术和并行计算技术可求解大规模的施工阶段分析模型 → 基于高性能GPU计算技术,2个小时即可完成具有12个施工阶段的100万个单元的施工阶段分析
[计算环境] [模型信息]
CPU Memory GPU
Intel Xeon X5647 2.93 GHz 8 core DDR3 48 Gigabyte Nvidia Tesla C2075
单元 节点 自由度
1,091,362 个 191,537 个 574,611 个
施工阶段 使用内存 线程数量
12 stages 24 Gigabyte 4个
ij 전응력 全应力
有效应力 ij 유효 응력
'
se 유효 포화도 有效饱和度 pw 간극 수압 孔隙水压
有效饱和度(0~1)由孔隙水压决定 → 在计算自重、不排水刚度和内力时使用部分饱和特性 → 可用于线性、非线性、固结、完全耦合分析中
nS ,
体积含水率 체적 함수비, n 공극률 , S 포화도 体积含水率, n 孔隙率, S 饱和度
更多的分析方法
→ 应力-渗流完全耦合分析 可更准确的考虑地下水变化的影响 → 非线性动力分析 可进行更准确的地震影响分析 → 动力荷载作用下的边坡稳定分析 边坡稳定分析的新分析理念 → 原场地分析 所有分析功能均支持将原场地状态作为初始状态
[考虑材料非线性的动力分析]
[改善了固结分析的计算内核]
[将原场地作为初始状态的岩土分析]
03 支持更多的岩土特性
04 非线性计算 05 边坡稳定分析
06 动力非线性分析
07 应力-渗流完全耦合分析
更多的岩土特性
→ 孔隙率对渗流系数的影响 岩土受压时孔隙率会减少,渗透能力会降低
考虑渗流系数变化的例题
→ 对填土的固结分析 可考虑自重和水位
e / ck k 10 k
使用孔隙率-渗流系数比函数
Save Element Matrix
[Post process]
Element Computing
Solve Linear Equation
根据硬件的内存使用状态,自动调整In Core 存储 模式和Out Of Core存储模式
Ku=f E -f I
Nonlinear load increment
Normalized stress
1.20 1.00 0.80 QUAD4 (Radial) 0.60 QUAD4 (Tangential) 0.40 0.20 0.00 0 2 4 Radius [m] 6 8 10 Analytical (Radial) Analytical (Tangential)
In Core 或Out of Core
Update stiffness and residual(K, f )
8
线性时程分析例题
→ 使用64位程序框架进行分析 → 比较在单元计算、振型计算中的并行计算效果
700 sec 600 sec
直接积分法
591 sec 498 sec
[计算环境]
500 sec 400 sec
sat
→ 随时间变化的填土沉降量分析
Settlement w.r.t Permeability coefficient type
0.00
→ 完成沉降时间会增加 → 用于固结分析和应力渗流完全耦合分析 → 在多孔隙材料定义中增加了选项
Settlement (m)
0
-0.20 -0.40 -0.60 -0.80 -1.00 -1.20
06 动力非线性分析
07 应力-渗流完全耦合分析
使用亚松弛法进行渗流分析
→ 非饱和区域的渗透系数具有较明显的非线性特性 砂土等孔隙较多的材料的收敛问题
50 Units: m S Z B
可以使用亚松弛法解决渗透系数的急剧变化问题
50
X
A
100
k p ,n 10
m
k p ,n1 k p ,n 10
400
Pressure [kPa]
300
200
100
Reference QUAD-8 HEXA-20
0 0.00
0.02
0.04
0.06
0.08
0.10
0.12
Vertical Displacement [m]
6
01 分析功能的改善 02 提供高性能计算功能
03 分析中可考虑更多的岩土特性参数
04 支持更多的非线性计算方法 05 边坡稳定分析的新理念
x 2.59
359 sec
[GTS 和GTS NX的结果比较]
振型叠加法
GTS
GTS NX P1
GTS NX P4
9
基于MultiCore/GPU的并行计算
→ 支持Multi Core CPU的并行计算(多核多线程) → 使用高性能GPU高速求解方程组 最大限度的使用硬件资源
[Multi threaded element computation]
1.80 1.60 1.40
500 600
Terzaghi, K. and Peck, R. B., Soil Mechaanics in Engineering Practice, 1967
Shield, R. T., and Drucker, D. C., The Application of Limit Analysis to Punch-Indentation Problems, Journal of Applied Mechanics, 1953
时间步骤 振型数量
100 step 100 个 5 Gigabyte
直接积分法
1000 900 800 700 600 500 400 300 200 100 0 sec sec sec sec sec sec sec sec sec sec sec
振型叠加法
929 sec 702 sec