jaabell

Animate transient results with gmsh2opensees

Another tutorial on leveraging the power of OpenSeesPy with gmsh and my library gmsh2opensees. In this one we also look into performing transient analysis in OpenSeesPy with a simple model.

Using OpenSeesPy with gmsh in python

In this tutorial, I teach you how to create wonderful continuum FE models in OpenSeesPy using gmsh as pre and post-processor. We'll be using my new module gmsh2opensees.

Advanced seismic analysis with H5DRM OpenSees load pattern

Check out the new H5DRM load pattern for DRM-based seismic analysis on OpenSees. You can find the documentation and all the example files here.

DRM stands for the Domain Reduction Method. It is a game-changer for seismic analysis. The DRM allows the analyst to subject their soil-structure systems to earthquake input featuring rich 3D motions. That is, you're no-longer restricted to assuming that the in-coming seismic wave-field is a planar wave, which is what 99% of analysis nowadays still assume. This is very important for several types of structures, for example as I showed was true for nuclear power plants subjected to near-field earthquake motions.

Building OpenSees and OpenSeesPy on Ubuntu 22.04 using CMAKE

In this short video I teach you how to compile the latest OpenSees and OpenSeesPy on Ubuntu Linux 22.04.

Compilation is setup with no "bells and whistles", meaning its not optimized in any way and does not include any extra features.

The history file mentioned in the video can be found here and this is the CMakeLists.tex file.

EOSD2022 Recap

The Eurasian OpenSees days 2022 conference was a lot of fun! Especially because I had the opportunity to go with two students: Omar Oyarce and Alberto Hurtado (shown in the picture above with professor Pedro Arduino from University of Washington, Seattle) and show them what this little academic niche is all about.

EOSD2022 Summer School SSI Modeling

Short applicative lecture given at Eurasian OpenSees Days 2022 in Turin, Italy on July 6th.

Supplementary material:

• Slides
• mshopstools.py - A toy library to interface gmsh with openseespy.
• workshop-ssi1.geo - SSI gmsh .geo (geometry) file for the application.
• workshop-ssi1.msh - SSI gmsh .msh (mesh) file for the application, generated from the above file without any modifications.
• workshop-ssi1.py - Python script that reads the .msh file and creates and runs the OpenSees model using openseespy.

Buidling a computer cluster for OpenSees

As part of my recently-awarded research grant where I'll study the effects of near-field earthquakes in Santiago Chile (more on that in the future, maybe) I'm building a small cluster of computers (a beowulf cluster).

New 2021 additions to OpenSees

2021 has almost gone, I haven't written on my blog almost all year, but there has been progress! I've made a bunch of small contributions this year to OpenSees (see all contributions here) but two big ones: a new material model for high-cyclic analysis of sands (SANISAND-MS with my collaborators at TU Delft, NGI and Bristol, see here) and added isogeometric analysis (IGA) of shells capabilities. This last contribution is very recent and is still (as of writing this) in pull request mode. I'll let some time pass, for possible review, before I definitely merge into main OpenSees.

Temas de Magister/Memoria 2022

Al final de la licenciatura (terminando "Métodos Computacionales en Obras Civiles") el alumno UANDES de la especialidad de Obras Civiles debe considerar si se va por la mención de "Estructuras" o por la de "Diseño de Proyecto de Infraestructuras". Además, si extenderá su carrera con un magister de continuidad o cerrará con un proyecto de título (memorias en la jerga antigua).

Esta es mi oferta de temas de magister o proyecto de título:

Magister

Se busca alumno de Magister para el siguiente tema:

"Modelación de alta-fidelidad para análisis de estructuras sometidas a sismos de campo cercano "

Temas de Magister/Memoria 2021

Al final de la licenciatura el alumno UANDES debe comenzar a considerar si extenderá su carrera con un magister de continuidad o cerrará con un proyecto de título. En cualquier caso se deben considerar posibles temas a estudiar en esta etapa.

Esta es mi oferta de temas de magister o proyecto de título. Para más detalles y/o ayuda a discernir esta etapa, no dudes en contactarme directamente para conversar y aclarar dudas.

Magister

Se busca alumno de Magister para el siguiente tema:

"Modelación de alta-fidelidad para análisis de ciclo de vida de puentes chilenos"

Compiling OpenSees 3.2.0 on Ubuntu 20.10: a tutorial and commentary

In this short video I teach you how to compile OpenSees 3.2.0 on Ubuntu Linux 20.10.

Compilation is setup with no "bells and whistles", meaning its not optimized in any way and does not include any extra features.

The history file mentioned in the video can be found here and this is the Makefile.def.

Check out Silvia's Brainery and Michael Scott's blog for more awesome OpenSees stuff!

STKO OpenSees review

A created a short #review of STKO: a pre/post processor for #OpenSees which I think is really cool and brings fresh new stuff into OpenSees to benefit all. https://t.co/EnMAVAnRya

— José Antonio Abell (@RealJoseAbell) November 30, 2020

Nicolás Galano defends his MS Thesis

Congratulations Nico!

Remember this tweet?

I'm proud to announce the MS defense of my (co)-MS student Nicolás Galano, who worked on extending #PFEM on #OpenSees for scour modeling.

We stand on the shoulder of giants. Original #PFEM implementation by @mikusscott and @MinjieZhuOS. Thanks guys! https://t.co/K6rfp8FjPj pic.twitter.com/AiwGGnSFQH

— José Antonio Abell (@RealJoseAbell) September 9, 2020

Great results! pic.twitter.com/L6NikXiqh3

— José Antonio Abell (@RealJoseAbell) September 9, 2020

Paper04 - Modelling the cyclic ratcheting of sands through memory-enhanced bounding surface plasticity.

This is the first paper fruit of the on-going collaboration with Federico Pisanò at TU Delft. Here we tackle the problem of formulating a constitutive model, derived from the SAniSand framework by Yannis Dafalias, which can replicate the ratcheting phenomenon on sands.

What is ratcheting?

Large earthquake in Mexico triggers tsunami warnings

The magnitude I gave came from the Tsunami threat record. pic.twitter.com/2h5dppVLLQ

— José A. Abell 🇻🇦🇨🇱 (@jose_abell) June 23, 2020

#EngineeringWizardry | 00 | new series

Stuck teaching during COVID-19 lockdowns in Chile has reminded me one crucial thing about statics (the class I'm teaching this first semester in the southern hemisphere): it can be very frustrating for students. This frustriation comes from unmet expectations. They're just finishing their early maths and physics formation (Calculus, Linear Algebra, Differential Equations, physics etc. ), an ideal world where problems have nice answers such as \(\sqrt{2} \pi\) or no answer at all. Then, along comes statics, and the world seems upside down.

Using my gmshtranslator python tool to interface gmsh with opensees.

I wrote the gmshtranslator tool a while back during my …

Paper03 - Earthquake Soil Structure Interaction of Nuclear Power Plants, differences in response to 3-D, 3 × 1-D, and 1-D excitations

Summary: Finally! First paper since PhD was accepted for publication …

Video of 2-D SSI model of an RC Shear Wall Building

The RC shear-wall building is modeled in OpenSees using non-linear fiber based beam column elements. Soil is modeled as a continuum using quad elements and linear stress-strain relationship. Soil shear wave-speed is varied in depth such as to obtain a \(V_{s30}\) consistent with a class B site according to chilean seismic code. Lysmer-Kulhemeyer dashpots are used along the soil boundary to model seismic radiation and earthquake wave-field input.

The performance of the building will be assessed for varying site fundamental periods. This is an aspect of SSI that is not covered by the chilean seismic code, and has been shown to be a problem in past earthquakes.

This is part of an on-going study with Prof. Carolina Magna from Adolfo Ibañez Unversity and her MS student Miguel Ángel Rodriguez from UDP.

New cable element in OpenSees

With my student, Pablo Ibañez, we implemented a new catenary cable element in OpenSees. This element is based on the work by Salehi et al[1].

The stiffness of this element is obtained using a flexibility formulation. Basically the shape of the cable is determined by the integral:

Where \(\vect{x_1}\) is the position of the first node of the cable, \(\vect{w}\) is the weight vector in each direction, \(EA\) is the stiffness, \(\alpha \Delta T\) is the change in strain due to temperature and \(\vect{f}\) is the force vector at the start node. This equation is iterated (with the forces as variable) upon until the the shape of the cable matches the nodal postiions imposed by the finite element program (trial displacements). Then it is used to derive a stiffness matrix.

The element, as is, passes all our static verification tests. With the additional assumption of a lumped-mass matrix, we're currently working on a dynamic verification suite as well as some validation experiments.

The animation above was created using OpenSees to simulate the cable and Blender to render it.

Reference

[1] Salehi Ahmad Abad, M., Shooshtari, A., Esmaeili, V., & Naghavi Riabi, A. (2013). Nonlinear analysis of cable structures under general loadings. Finite Elements in Analysis and Design, 73, 11–19. https://doi.org/10.1016/j.finel.2013.05.002

OpenSees coming to python!

One of the gripes a lot of people have with …

Finite Elements for Shallow Water Equations

Talk Given at IngeoKring 2016 - Physics-Based Earthquake-Soil-Structure Interaction for Near-Field Induced Seismicity

This is a talk given at the IngeoKring 2016 Autumn …

NTS-02. On Rayleigh damping coefficients for FE analysis

Testing IPython Notebook integration into Pelican

SNE # 01. Example of ESSI simulation and visualization with visitESSI

These are the results of a simple elastic-domain simulation. Mesh …

NTS (Note-To-Self): Creating dynamically linked libraries

A nice extensive tutorial can be found here.

In a …

SNE # 0. Stochastic inverse pendulum

This is the first installment of "Small Numerical Experiments" (SNE …