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>Finite Element Mesh of a beam</title>
<date>January 1, 2010</date>
<author>R. Touzani</author>
</info>
<Project name="proj">
<mesh_file>proj-5x10.m</mesh_file>
<plot_file>proj.pl</plot_file>
<time_step>0.01</time_step>
<max_time>1.0</max_time>
<verbose>0</verbose>
<output>0</output>
<save>1</save>
</Project>
<Prescription>
<BoundaryCondition code="1">-tanh(10)*(exp(t)-1)</BoundaryCondition>
<BoundaryCondition code="2"> tanh(10)*(exp(t)-1)</BoundaryCondition>
<Source>tanh(10*y)*(exp(t)+200*(exp(t)-1)/(cosh(10*y)*cosh(10*y)))</Source>
</Prescription>
</OFELI_File>
|
In summary, this file looks mainly like the one in the previous example. We show here
in addition how to prescribe a boundary condition by an algebraic expression. The
OFELI library is equipped with the expression parser fparser.
We use here the variables x, y, z,
and t for space and time variables.