weac 2.0.4

Weak layer anticrack nucleation model

WEAC  ·  Weak Layer Anticrack Nucleation Model

Implementation of closed-form analytical models for the analysis of dry-snow slab avalanche release.

View demo · Report bug · Request feature · Read the docs

Contents

1. About the project
2. Installation
3. Usage
4. Roadmap
5. Release history
6. How to contribute
7. License
8. Contact

About the project

WEAC implements closed-form analytical models for the mechanical analysis of dry-snow slabs on compliant weak layers, the prediction of anticrack onset, and, in particular, allows for the analysis of stratified snow covers. The model covers propagation saw tests (a), and uncracked (b) or cracked (c) skier-loaded buried weak layers.

Please refer to the companion papers for model derivations, illustrations, dimensions, material properties, and kinematics:

• Rosendahl, P. L., & Weißgraeber, P. (2020). Modeling snow slab avalanches caused by weak-layer failure – Part 1: Slabs on compliant and collapsible weak layers. The Cryosphere, 14(1), 115–130. https://doi.org/10.5194/tc-14-115-2020
• Rosendahl, P. L., & Weißgraeber, P. (2020). Modeling snow slab avalanches caused by weak-layer failure – Part 2: Coupled mixed-mode criterion for skier-triggered anticracks. The Cryosphere, 14(1), 131–145. https://doi.org/10.5194/tc-14-131-2020

Written in Python and built with Atom, GitKraken, and Jupyter. Note that release v1.0 was written and built in MATLAB.

Installation

Install globally using the pip Package Installer for Python

pip install weac

or clone the repo

git clone https://github.com/2phi/weac

for local use.

Needs

• Numpy for matrix operations
• Scipy for solving optimization problems
• Pandas for data handling
• Matplotlib for plotting

Usage

The following describes the basic usage of WEAC. Please refer to the demo for more examples and read the documentation for details.

Load the module.

import weac

Choose a profile from the database (see demo) or create your own.

myprofile = [[180, 100],
             [190,  40],
             [230, 130],
             [380,  20],
             [270, 100]]

Create a model instance with optional custom layering.

skier = weac.Layered(system='skier', layers=myprofile)

Calculate lists of segment lengths, locations of foundations, and position and magnitude of skier loads from the inputs total length L (mm), crack length a (mm), and skier weight m (kg). We can choose to analyze the situtation before a crack appears even if a crack length > 0 is set by replacing the 'crack' key thorugh the 'uncracked' key.

segments = skier.calc_segments(L=10000, a=300, m=80)['crack']

Assemble the system of linear equations and solve the boundary-value problem for the free constants C providing the inclination phi in degrees.

C = skier.assemble_and_solve(phi=38, **segments)

Prepare the output by rasterizing the solution vector at all horizontal positions xq. The result is returned in the form of the ndarray zq. We also get xb that only contains x-coordinates that lie on a foundation.

xq, zq, xb = skier.rasterize_solution(C=C, phi=38, **segments)

Visualize the results.

# Visualize deformations as a contour plot
weac.plot.contours(skier, x=xq, z=zq, window=200, scale=100)

# Plot slab displacements
weac.plot.displacements(skier, x=xq, z=zq, **segments)

# Plot weak-layer stresses
weac.plot.stresses(skier, x=xb, z=zq, **segments)

Roadmap

See the open issues for a list of proposed features and known issues.

v2.2

• Finite fracture mechanics implementation for layered snow covers

v2.1

• Flip coordinate system and order of layers of the slab
• Allow for cutting PSTs from left or right side

v2.0

• Completely rewritten in Python
• Coupled bending-extension ODE solver implemented
• Stress analysis of arbitrarily layered snow slabs
• FEM validation of
  • displacements
  • weak-layer stresses
  • energy release rates in weak layers
• Documentation
• Demo and examples

Release history

v1.0

• Written in MATLAB
• Deformation analysis of homogeneous snow labs
• Weak-layer stress prediction
• Energy release rates of cracks in weak layers
• Finite fracture mechanics implementation
• Prediction of anticrack nucleation

How to contribute

1. Fork the project
2. Create your feature branch (git checkout -b feature/AmazingFeature)
3. Commit your changes (git commit -m 'Add some AmazingFeature')
4. Push to the branch (git push origin feature/AmazingFeature)
5. Open a pull request

License

Copyright 2phi GbR, 2021.

We currently do not offer an open source license. Please contact us for private licensing options.

Contact

E-mail: [e-mail] · Web: https://2phi.de · Project Link: https://github.com/2phi/weac · Project DOI: http://dx.doi.org/10.5281/zenodo.5773113