1,模型尽量大一些,接近实际为最佳;边界施加围岩预压力为其次,边界施加无反射边界为最后得方法;精度依次降低;
2,模型不进行切分,整体一个模型(研究振动,不要因为切分模型处理网格而导致传递力出现不连续)(爆破开挖失效,需要六面体,只是研究振动,如果模型很难处理为六面体就用10号纯四面体即可--WB划分10号算法得纯四面体)
3,如果使用四面体网格,材料不能使用jc本构,可以使用PK本构;
4,添加阻尼;阻尼分为刚度阻尼及振动阻尼;
1)当定义为all part时,自动使用质量加权阻尼,对于低频震荡十分有效;
*DAMPING_GLOBAL
0 0.005
2)指定了阻尼曲线LCID时,自动使用质量加权阻尼,对于低频震荡十分有效【可以只用全局阻尼】;
*DAMPING_PART_MASS
$# pid lcid sf flag
4 2 1000.0 0
*DAMPING_PART_MASS
$# pid lcid sf flag
1 2 1000.0 0
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp
2,0,1,1,0,0,0
0,0.000
500,1
501,0
1000,0
3)不使用阻尼曲线号且制定了阻尼系数COEF,刚度加权阻尼被使用于特定part,对于高频震荡有效;
*DAMPING_PART_STIFFNESS
$# pid coef
3 0.1
参考k文件如下:
$
*MAT_PLASTIC_KINEMATIC
1,2.65,0.374,0.25,0.0002452,0.045,0.00
0,0,0
$
*PART
$ partname
rocky.1
$ ID secid mid eosid hgid grav adopt tmid
1 1 1 0 0 0 0 0
*SECTION_SOLID
$ ID elform aet unused1
2 10 0
*MAT_PLASTIC_KINEMATIC
2,2.65,0.374,0.25,0.0002452,0.045,0.00
0,0,0
*PART
$ partname
xieshuidongchuzhi
$ ID secid mid eosid hgid grav adopt tmid
2 2 2 0 0 0 0 0
*SECTION_SOLID
$ ID elform aet unused1
3 10 0
*MAT_PLASTIC_KINEMATIC
3,2.65,0.374,0.25,0.0002452,0.045,0.00
0,0,0
*PART
$ partname
hengtongdaochuzhi
$ ID secid mid eosid hgid grav adopt tmid
3 3 3 0 0 0 0 0
*SECTION_SOLID
$ ID elform aet unused1
4 1 0
*MAT_PLASTIC_KINEMATIC
4,2.65,0.374,0.25,0.0002452,0.045,0.00
0,0,0
*PART
$ partname
tianchong
$ ID secid mid eosid hgid grav adopt tmid
4 4 4 0 0 0 0 0
*SECTION_SOLID
$ ID elform aet unused1
5 1 0
*MAT_PLASTIC_KINEMATIC
5,2.65,0.374,0.25,0.0002452,0.045,0.00
0,0,0
*PART
$ partname
erchen
$ ID secid mid eosid hgid grav adopt tmid
5 5 5 0 0 0 0 0
$
$
$
*DAMPING_GLOBAL
0 0.005
*DAMPING_PART_STIFFNESS
$# pid coef
1 0.1
*DAMPING_PART_STIFFNESS
$# pid coef
2 0.1
*DAMPING_PART_STIFFNESS
$# pid coef
3 0.1
*DAMPING_PART_STIFFNESS
$# pid coef
4 0.1
*DAMPING_PART_STIFFNESS
$# pid coef
5 0.1
$
$
*ALE_MULTI-MATERIAL_GROUP
6 1
7 1
8 1
9 1
10 1
11 1
$
*SET_PART_LIST
1
1
*SET_PART_LIST
2
6,7,8,9,10,11
*CONSTRAINED_LAGRANGE_IN_SOLID
1 2 0 0 2 5 2 0
0 0 0.15 0.1
2 0.3
$
*CONTROL_ALE
-1 1 2 -1.000000 0.000000 0.000000
0.000000 0.000000 0.000000
$
$
*SECTION_SOLID_ALE
$ ID elform aet unused1
6 11
*MAT_NULL
6 1.2929e-3 0 0
*EOS_LINEAR_POLYNOMIAL
6 -1e-6 0.00 0.00 0.00 0.4 0.4 0.00
2.068e-6 1.0
*MAT_HIGH_EXPLOSIVE_BURN
10 1.63 0.693 0.27 0.00E 0
*EOS_JWL
10 3.71 7.43E-02 4.15 0.950 0.30 7.0E-02 1.00
*initial_detonation
$10为part号
10,-1239,200,0,0
$
*PART
$ partname
air.1
$ ID secid mid eosid hgid grav adopt tmid
6 6 6 6 0 0 0 0
*PART
$ partname
air.2
$ ID secid mid eosid hgid grav adopt tmid
7 6 6 6 0 0 0 0
*PART
$ partname
air.3
$ ID secid mid eosid hgid grav adopt tmid
8 6 6 6 0 0 0 0
*PART
$ partname
air.4
$ ID secid mid eosid hgid grav adopt tmid
9 6 6 6 0 0 0 0
*PART
$ partname
tnt
$ ID secid mid eosid hgid grav adopt tmid
10 6 10 10 0 0 0 0
*PART
$ partname
air.5
$ ID secid mid eosid hgid grav adopt tmid
11 6 6 6 0 0 0 0
结果如下,高频震荡基本消失,阻尼起到应有的作用。