Dynamic models for impact-initiated stress waves through snow columns

The objective of this research is to model snow’s response to dynamic, impact loading. Two constitutive relationships are considered: elastic and Maxwell-viscoelastic. These material models are applied to laboratory experiments consisting of 1000 individual impacts across 22 snow column configurations. The columns are 60 cm tall with a 30 cm by 30 cm cross-section. The snow ranges in density from 135-428 kg m-3 and is loaded with both short-duration (~1 ms) and long-duration (~10 ms) impacts. The Maxwell-viscoelastic model more accurately describes snow’s response because it contains a mechanism for energy dissipation, which the elastic model does not. Furthermore, the ascertained model parameters show a clear dependence on impact duration; shorter duration impacts resulted in higher wave speeds and greater damping coefficients. The stress wave’s magnitude is amplified when it hits a stiffer material because of the positive interference between incident and reflected waves. This phenome..., Collected in the Subzero Research Laboratory at Montana State University and processed with various Python and Matlab scripts., , # Dynamic models for impact-initiated stress waves through snow columns

This dataset contains supporting information for the Journal of Glaciology paper titled \"Dynamic models for impact-initiated stress waves through snow columns\" authored by S.V. Verplanck and E.E. Adams.

## Description of the data and file structure

The scripts to run the models are in two folders: \"FD Model\" and \"FE Model\". The finite difference (FD) model is scripted and executed Matlab R2023A. The finite element (FE) model is scripted in Python and executed in Abaqus 2022.

The \"FD Model\" folder contains a pdf with detailed documentation of the finite difference method implementation. This file is called \"FD_Documention.pdf\". There is also a top level matlab script titled \"runModel.m\" which calls \"wave_1D.m\".

The \"FE Model\" folder contains 12 python scripts. The filenaming convention is \"loop_\" + \"longDur\" or \"shortDur\" + \"E\" or \"VE\" + \"low\" or \"mid\" or \"high\". The permutations of this naming convention is s...

本研究的目标是构建雪体对动态冲击荷载的响应模型。本次研究考虑两种本构关系（constitutive relationships）：弹性模型与麦克斯韦粘弹性模型（Maxwell-viscoelastic）。将上述两种材料模型应用于实验室实验，实验包含针对22种雪柱构型的1000次独立冲击。雪柱高度为60 cm，截面尺寸为30 cm × 30 cm；雪样密度范围为135~428 kg·m⁻³，分别接受短时长（约1 ms）与长时长（约10 ms）的冲击加载。
麦克斯韦粘弹性模型能够更准确地描述雪体的响应特性，因为其内置了能量耗散机制，而弹性模型不具备该机制。此外，所确定的模型参数清晰展现出对冲击时长的依赖性：短时长冲击会产生更高的波速与更大的阻尼系数。当应力波遭遇刚度更高的材料时，由于入射波与反射波之间的正干涉效应，其幅值会被放大，该现象……本数据集采集自蒙大拿州立大学（Montana State University）低温研究实验室（Subzero Research Laboratory），并通过多种Python与Matlab脚本进行处理。

# 雪柱冲击诱发应力波的动态模型
本数据集为发表于《冰川学杂志》（Journal of Glaciology）、题为《雪柱冲击诱发应力波的动态模型》的论文提供支撑信息，该论文作者为S.V. Verplanck与E.E. Adams。

## 数据与文件结构说明
用于运行模型的脚本分为两个文件夹："FD Model"（有限差分模型）与"FE Model"（有限元模型）。有限差分（finite difference, FD）模型通过Matlab R2023A编写并执行；有限元（finite element, FE）模型使用Python编写，依托Abaqus 2022完成运行。
"FD Model"文件夹包含一份详细阐述有限差分法实现细节的PDF文档，文件名为"FD_Documention.pdf"；同时设有顶层Matlab脚本"runModel.m"，该脚本会调用"wave_1D.m"。
"FE Model"文件夹包含12个Python脚本，其命名规则为"loop_"+"longDur"或"shortDur"+"E"或"VE"+"low"或"mid"或"high"，该命名规则的排列组合为……