[Getfem-commits] r4956 - in /trunk/getfem: doc/sphinx/source/tutorial/ interface/src/scilab/demos/ interface/tests/matlab/ interface/tests/py...

[Yves . Renard]

Author: renard
Date: Sat Apr 18 15:06:48 2015
New Revision: 4956

URL: http://svn.gna.org/viewcvs/getfem?rev=4956&view=rev
Log:
minor fix

Modified:
    trunk/getfem/doc/sphinx/source/tutorial/thermo_coupling.rst
    
trunk/getfem/interface/src/scilab/demos/demo_thermo_elasticity_electrical_coupling.sce
    
trunk/getfem/interface/tests/matlab/demo_thermo_elasticity_electrical_coupling.m
    
trunk/getfem/interface/tests/python/demo_thermo_elasticity_electrical_coupling.py
    trunk/getfem/tests/thermo_elasticity_electrical_coupling.cc

Modified: trunk/getfem/doc/sphinx/source/tutorial/thermo_coupling.rst
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/doc/sphinx/source/tutorial/thermo_coupling.rst?rev=4956&r1=4955&r2=4956&view=diff
==============================================================================
--- trunk/getfem/doc/sphinx/source/tutorial/thermo_coupling.rst (original)
+++ trunk/getfem/doc/sphinx/source/tutorial/thermo_coupling.rst Sat Apr 18 
15:06:48 2015
@@ -75,7 +75,7 @@
 
 .. math::
  
-  \bar{\sigma}(u) = \lambda^* \mbox{div}(u) I + 2\mu \bar{\varepsilon}(u) + 
\beta\theta I,
+  \bar{\sigma}(u) = \lambda^* \mbox{div}(u) I + 2\mu \bar{\varepsilon}(u) + 
\beta(T_0-\theta) I,
 
 :math:`\bar{\varepsilon}(u) = (\nabla u + (\nabla u)^T)/2` being the 
linearized strain tensor, :math:`I` the identity second order tensor and 
:math:`\lambda^*, \mu` being the |Lame| coefficients defined  by
 
@@ -87,7 +87,7 @@
 
 from :math:`E` the Young modulus and :math:`\nu` the Poisson ratio of the 
material.
 
-The term :math:`\beta\theta I` correspond to the thermal expansion term, where 
:math:`\beta = \alpha_{th} E/(1-2\nu)`, :math:`\alpha_{th}` being the thermal 
expansion coefficient.
+The term :math:`\beta(T_0-\theta) I` corresponds to the thermal expansion 
term, where :math:`\beta = \alpha_{th} E/(1-2\nu)`, :math:`\alpha_{th}` being 
the thermal expansion coefficient.
 
 
 The weak formulation
@@ -537,7 +537,7 @@
 
                 md.add_initialized_scalar_data("beta", alpha_th*E/(1-2*nu));
                 getfem::add_linear_generic_assembly_brick
-                  (md, mim, "beta*theta*Div_Test_u");
+                  (md, mim, "beta*(T0-theta)*Div_Test_u");
 ---------- 
----------------------------------------------------------------------------------------------------------------
 **Python** .. code-block:: python                                     
 
@@ -550,7 +550,7 @@
                 md.add_source_term_brick(mim, 'u', 'Fdata', RIGHT_BOUND)
 
                 md.add_initialized_data('beta', [alpha_th*E/(1-2*nu)])
-                md.add_linear_generic_assembly_brick(mim, 
'beta*theta*Div_Test_u')
+                md.add_linear_generic_assembly_brick(mim, 
'beta*(T0-theta)*Div_Test_u')
 ---------- 
----------------------------------------------------------------------------------------------------------------
 **Scilab** .. code-block:: matlab
 
@@ -563,7 +563,7 @@
                 gf_model_set(md, 'add source term brick', mim, 'u', 'Fdata', 
RIGHT_BOUND);
 
                 gf_model_set(md, 'add initialized data', 'beta', 
[alpha_th*E/(1-2*nu)]);
-                gf_model_set(md, 'add linear generic assembly brick', mim, 
'beta*theta*Div_Test_u');
+                gf_model_set(md, 'add linear generic assembly brick', mim, 
'beta*(T0-theta)*Div_Test_u');
 ---------- 
----------------------------------------------------------------------------------------------------------------
 **Matlab** .. code-block:: matlab
 
@@ -576,7 +576,7 @@
                 gf_model_set(md, 'add source term brick', mim, 'u', 'Fdata', 
RIGHT_BOUND);
 
                 gf_model_set(md, 'add initialized data', 'beta', 
[alpha_th*E/(1-2*nu)]);
-                gf_model_set(md, 'add linear generic assembly brick', mim, 
'beta*theta*Div_Test_u');
+                gf_model_set(md, 'add linear generic assembly brick', mim, 
'beta*(T0-theta)*Div_Test_u');
 ========== 
================================================================================================================
 
 Electric potential problem

Modified: 
trunk/getfem/interface/src/scilab/demos/demo_thermo_elasticity_electrical_coupling.sce
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/interface/src/scilab/demos/demo_thermo_elasticity_electrical_coupling.sce?rev=4956&r1=4955&r2=4956&view=diff
==============================================================================
--- 
trunk/getfem/interface/src/scilab/demos/demo_thermo_elasticity_electrical_coupling.sce
      (original)
+++ 
trunk/getfem/interface/src/scilab/demos/demo_thermo_elasticity_electrical_coupling.sce
      Sat Apr 18 15:06:48 2015
@@ -176,7 +176,7 @@
 
 // Thermal expansion term
 gf_model_set(md, 'add initialized data', 'beta', [alpha_th*E/(1-2*nu)]);
-gf_model_set(md, 'add linear generic assembly brick', mim, 
'beta*theta*Trace(Grad_Test_u)');
+gf_model_set(md, 'add linear generic assembly brick', mim, 
'beta*(T0-theta)*Trace(Grad_Test_u)');
 
 
 //

Modified: 
trunk/getfem/interface/tests/matlab/demo_thermo_elasticity_electrical_coupling.m
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/interface/tests/matlab/demo_thermo_elasticity_electrical_coupling.m?rev=4956&r1=4955&r2=4956&view=diff
==============================================================================
--- 
trunk/getfem/interface/tests/matlab/demo_thermo_elasticity_electrical_coupling.m
    (original)
+++ 
trunk/getfem/interface/tests/matlab/demo_thermo_elasticity_electrical_coupling.m
    Sat Apr 18 15:06:48 2015
@@ -180,7 +180,7 @@
 
 % Thermal expansion term
 gf_model_set(md, 'add initialized data', 'beta', [alpha_th*E/(1-2*nu)]);
-gf_model_set(md, 'add linear generic assembly brick', mim, 
'beta*theta*Div_Test_u');
+gf_model_set(md, 'add linear generic assembly brick', mim, 
'beta*(T0-theta)*Div_Test_u');
 
 
 %

Modified: 
trunk/getfem/interface/tests/python/demo_thermo_elasticity_electrical_coupling.py
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/interface/tests/python/demo_thermo_elasticity_electrical_coupling.py?rev=4956&r1=4955&r2=4956&view=diff
==============================================================================
--- 
trunk/getfem/interface/tests/python/demo_thermo_elasticity_electrical_coupling.py
   (original)
+++ 
trunk/getfem/interface/tests/python/demo_thermo_elasticity_electrical_coupling.py
   Sat Apr 18 15:06:48 2015
@@ -184,7 +184,7 @@
 
 # Thermal expansion term
 md.add_initialized_data('beta', [alpha_th*E/(1-2*nu)])
-md.add_linear_generic_assembly_brick(mim, 'beta*theta*Trace(Grad_Test_u)')
+md.add_linear_generic_assembly_brick(mim, 'beta*(T0-theta)*Trace(Grad_Test_u)')
 
 
 #

Modified: trunk/getfem/tests/thermo_elasticity_electrical_coupling.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/tests/thermo_elasticity_electrical_coupling.cc?rev=4956&r1=4955&r2=4956&view=diff
==============================================================================
--- trunk/getfem/tests/thermo_elasticity_electrical_coupling.cc (original)
+++ trunk/getfem/tests/thermo_elasticity_electrical_coupling.cc Sat Apr 18 
15:06:48 2015
@@ -237,7 +237,7 @@
   // Thermal expansion term
   md.add_initialized_scalar_data("beta", alpha_th*E/(1-2*nu));
   getfem::add_linear_generic_assembly_brick
-    (md, mim, "beta*theta*Trace(Grad_Test_u)");
+    (md, mim, "beta*(T0-theta)*Trace(Grad_Test_u)");
 
 
   //