1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
% Copyright (C) 2017-2020 Yves Renard.
%
% This file is a part of GetFEM++
%
% GetFEM++ is free software; you can redistribute it and/or modify it
% under the terms of the GNU Lesser General Public License as published
% by the Free Software Foundation; either version 3 of the License, or
% (at your option) any later version along with the GCC Runtime Library
% Exception either version 3.1 or (at your option) any later version.
% This program is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
% or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
% License and GCC Runtime Library Exception for more details.
% You should have received a copy of the GNU Lesser General Public License
% along with this program; if not, write to the Free Software Foundation,
% Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA.
% -*- matlab -*- (enables emacs matlab mode)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% parameters for program thermo_elasticity_electrical_coupling %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
epsilon = 1.; % Thickness of the plate (cm)
E = 21E6; % Young Modulus (N/cm^2)
nu = 0.3; % Poisson ratio
F = 100E2; % Force density at the right boundary (N/cm^2)
kappa = 4.; % Thermal conductivity (W/(cm K))
D = 10.; % Heat transfer coefficient (W/(K cm^2))
air_temp = 20.; % Temperature of the air in oC.
alpha_th = 16.6E-6; % Thermal expansion coefficient (/K).
T0 = 20.; % Reference temperature in oC.
rho_0 = 1.754E-8; % Resistance temperature coefficient at T0
alpha = 0.0039; % Second resistance temperature coefficient.
h = 1.; % Approximate mesh size
elements_degree = 2; % Degree of the finite element methods
export_mesh = 1; % Draw the mesh after mesh generation or not
solve_in_two_steps = 1; % Solve the elasticity pb separately or not
|
