Commit b65aed8f authored by Turnhout, M.C. van's avatar Turnhout, M.C. van
Browse files

update to current version and rerun sigma(18) simulations

parent a3d6c6cb
......@@ -129,6 +129,5 @@ octave-workspace
*.sig
documentation/*.f
documentation/paramspace/*.f
!documentation/param.inc
!documentation/pics/*.pdf
!documentation/paramspace/*.7z
......@@ -40,10 +40,11 @@ noi = [lebn; tebn; rebn; lebn(1)]; % nodes of interest
for m = 1:numel(mus)
param.expname = sprintf('lfmu_%02i_param', m);
% fid = fopen(sprintf('lfmu_%02i_results.txt', m), 'w');
% param.fiblab = [mus(m) pi/9 .75];
% TF_FOA2paraminc(param);
%
% fid = fopen(sprintf('lfmu_%02i_results.txt', m), 'w');
% param.fiblab = [mus(m) pi/9 .75];
% TF_FOA2paraminc(param);
% TF_writeFEM(param, 1, 1);
% for s = 2:numel(sigma)
% param.sigmax(1, 1) = sigma(s);
% TF_writeFEMinc(param, 1, 1);
......@@ -68,14 +69,14 @@ for m = 1:numel(mus)
% fclose(fid);
%
% if m ~= 9
% % rerun sigma(18), last correct for entire series
% param.sigmax(1, 1) = sigma(18);
% TF_writeFEMinc(param, 1, 1);
% % rerun sigma(18), last correct for entire series
% param.sigmax(1, 1) = sigma(18);
% TF_writeFEMinc(param, 1, 1);
%
% [status, res] = TF_runFEM(param, 1, 1);
% % print job status
% idx = regexp(res, '\n');
% fprintf('%s', res(idx(end-1)+1:end));
% [status, res] = TF_runFEM(param, 1, 1);
% % print job status
% idx = regexp(res, '\n');
% fprintf('%s', res(idx(end-1)+1:end));
% end
figure(1)
......@@ -186,7 +187,7 @@ for m = 1:numel(mus)
set(gca, 'view', [90 -1])
figure(12)
kl = sqrt(sum(kappa(:, 1:3).^2, 2));
kr = sqrt(sum(kappa(:, 4:6).^2, 2));
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 11:29:58
% @ 17-Nov-2018 11:50:33
param.expname = 'lfmu_01_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 11:38:46
% @ 17-Nov-2018 11:57:09
param.expname = 'lfmu_02_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 11:47:37
% @ 17-Nov-2018 12:05:33
param.expname = 'lfmu_03_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 11:56:53
% @ 17-Nov-2018 12:15:08
param.expname = 'lfmu_04_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 12:06:37
% @ 17-Nov-2018 12:25:56
param.expname = 'lfmu_05_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 12:14:59
% @ 17-Nov-2018 12:35:08
param.expname = 'lfmu_06_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 12:23:18
% @ 17-Nov-2018 12:41:48
param.expname = 'lfmu_07_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 12:32:38
% @ 17-Nov-2018 12:50:05
param.expname = 'lfmu_08_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 12:39:55
% @ 17-Nov-2018 12:56:11
param.expname = 'lfmu_09_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
% written by TF_writeparam
% http://mategit.wfw.wtb.tue.nl/STEM/TFlab
% @ 13-May-2018 12:44:54
% @ 17-Nov-2018 13:07:16
param.expname = 'lfmu_10_param';
param.nfilm = 1;
param.nframe = 3;
param.refframe = 3;
param.fibdef = {'fiblab'};
param.Efilm = 1.520000e+00;
param.vfilm = 4.900000e-01;
param.tfilm = 6.780000e-03;
param.Ecells = 7.000000e-04;
param.vcells = 3.000000e-01;
param.tcells = 2.300000e-03;
param.FEMbuild = 'matlab';
param.nelx = 10;
param.nely = 20;
param.nelz = 1;
param.mininc = 1.000000e-05;
param.maxinc = 1000;
param.nCPU = 4;
param.submit = 0;
param.curlcheck = 1;
param.inisigs = [7.500000e-04, 1.500000e-03];
param.cutback = 3.000000e-01;
......@@ -19,7 +28,7 @@ param.sigesttol = 1.000000e-06;
param.simfailtol = 5;
param.fancyest = 5.000000e-01;
param.instpath = '/home/mark/git/TFlab/';
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc'};
param.intvars = {'refframe', 'nfilm', 'nframe', 'stack', 'curlcheck', 'simfailtol', 'maxinc', 'nelx', 'nely', 'nelz', 'nCPU', 'submit'};
param.TFlabUrl = 'http://mategit.wfw.wtb.tue.nl/STEM/TFlab';
param.TFlabBranch = 'testing';
param.TFlabVersion = '0.1+';
......
Ep = 1.520E+00 ! [MPa] film modulus
vp = 0.49 ! [-] film Poisson ration
Ec = 7.000E-04 ! [MPa] cell layer modulus
vc = 0.30 ! [-] cell layer Poisson ration
sigmax = 1.000000E-02 ! [MPa] max cell stress
phi_sf(1) = 3.333E-02 ! fraction of fibres in direction pi*1/30
phi_sf(2) = 3.333E-02 ! fraction of fibres in direction pi*2/30
phi_sf(3) = 3.333E-02 ! fraction of fibres in direction pi*3/30
phi_sf(4) = 3.333E-02 ! fraction of fibres in direction pi*4/30
phi_sf(5) = 3.333E-02 ! fraction of fibres in direction pi*5/30
phi_sf(6) = 3.333E-02 ! fraction of fibres in direction pi*6/30
phi_sf(7) = 3.333E-02 ! fraction of fibres in direction pi*7/30
phi_sf(8) = 3.333E-02 ! fraction of fibres in direction pi*8/30
phi_sf(9) = 3.333E-02 ! fraction of fibres in direction pi*9/30
phi_sf(10) = 3.333E-02 ! fraction of fibres in direction pi*10/30
phi_sf(11) = 3.333E-02 ! fraction of fibres in direction pi*11/30
phi_sf(12) = 3.333E-02 ! fraction of fibres in direction pi*12/30
phi_sf(13) = 3.333E-02 ! fraction of fibres in direction pi*13/30
phi_sf(14) = 3.333E-02 ! fraction of fibres in direction pi*14/30
phi_sf(15) = 3.333E-02 ! fraction of fibres in direction pi*15/30
phi_sf(16) = 3.333E-02 ! fraction of fibres in direction pi*16/30
phi_sf(17) = 3.333E-02 ! fraction of fibres in direction pi*17/30
phi_sf(18) = 3.333E-02 ! fraction of fibres in direction pi*18/30
phi_sf(19) = 3.333E-02 ! fraction of fibres in direction pi*19/30
phi_sf(20) = 3.333E-02 ! fraction of fibres in direction pi*20/30
phi_sf(21) = 3.333E-02 ! fraction of fibres in direction pi*21/30
phi_sf(22) = 3.333E-02 ! fraction of fibres in direction pi*22/30
phi_sf(23) = 3.333E-02 ! fraction of fibres in direction pi*23/30
phi_sf(24) = 3.333E-02 ! fraction of fibres in direction pi*24/30
phi_sf(25) = 3.333E-02 ! fraction of fibres in direction pi*25/30
phi_sf(26) = 3.333E-02 ! fraction of fibres in direction pi*26/30
phi_sf(27) = 3.333E-02 ! fraction of fibres in direction pi*27/30
phi_sf(28) = 3.333E-02 ! fraction of fibres in direction pi*28/30
phi_sf(29) = 3.333E-02 ! fraction of fibres in direction pi*29/30
phi_sf(30) = 3.333E-02 ! fraction of fibres in direction pi*30/30
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment