Loading

   

Solid and Structural Mechanics Demo 1: A linear elasticity problem with planar deformations

This code solves linearized elastostatics problems with planar deformations.

  

An example

Let us run this program with the data presented in the following project file:

<?xml version="1.0" encoding="ISO-8859-1" ?>
<OFELI_File>
<info>
   <title></title>
   <date></date>
   <author></author>
</info>
<Project name="beam">
   <verbose>1</verbose>
   <output>1</output>
   <save>1</save>
   <plot>1</plot>
   <init>1</init>
   <mesh_file>beam.m</mesh_file>
   <plot_file>beam.d</plot_file>
   <prescription_file>beam.pr</prescription_file>
</Project>
</OFELI_File>

The only thing that may differ from other seen project files is that prescription data is given in a separate file: beam.pr. Let us give this one:

<?xml version="1.0" encoding="ISO-8859-1" ?>
<OFELI_File>
<info>
   <title></title>
   <date></date>
   <author></author>
</info>
<Prescription>
   <BoundaryCondition code="1">0.0</BoundaryCondition>
   <Traction code="1" dof="2">-1000.0</Traction>
</Prescription>
</OFELI_File>

Clearly, we have prescribed by this file a (Dirichlet) boundary condition that enforces a value of 0.0 for degrees of freedom with code 1. In addition, we have prescribed a traction of value -1000 for the second degrees of freedom (y-traction) having code 1.

 
Copyright © 1998-2016 Rachid Touzani