%--------------------------------------------------------------------------
% PURPOSE
% TIMO_Multifiber - DYNAMIC - DYSBAC
%
%--------------------------------------------------------------------------
% REFERENCES
% Chaymaa LEJOUAD
% 05-02-2018
%--------------------------------------------------------------------------
% COMMENTS
%
%
%
%--------------------------------------------------------------------------
% MIT License
%
% Copyright (c) 2018 Benjamin Richard
%
% Permission is hereby granted, free of charge, to any person obtaining a
% copy of this software and associated documentation files (the "Software
% "), to deal in the Software without restriction, including without
% limitation the rights to use, copy, modify, merge, publish, distribute,
% sublicense, and/or sell copies of the Software, and to permit persons
% to whom the Software is furnished to do so, subject to the following
% conditions:
%
% The above copyright notice and this permission notice shall be included
% in all copies or substantial portions of the Software.
%
% THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
% OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
% MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
% IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY
% CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT,
% TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
% SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
%--------------------------------------------------------------------------
%% Clearing off
%% Global variables
global ME TP
%% Definition of the section
FILE = '103_1.mail';
ME = INPUT.ACQU(FILE,'MAIL');
%% Definition of the model - cross section scale
MO1 = MODEL('SBET' ,'MECHANICS','ELASTICITY','ISOTROPIC','DAMAGE','MAZARS','QUAS');
MO2 = MODEL('SAINF','MECHANICS','ELASTICITY','ISOTROPIC','PLASTICITY','VONMISES','POJS');
MO3 = MODEL('SASUP','MECHANICS','ELASTICITY','ISOTROPIC','PLASTICITY','VONMISES','POJS');
MOT = [MO1 [MO2 MO3]];
%% Definition of the materials - cross section scale
MA1 = CHAMELEM.MATE(MO1,'youn',36000e6,'nu',0.2,'rho',...
2500,'ay',0.83,'az',0.83,'at',0.8,'k0',1.0e-4,'bt',12000,...
'beta',1.06,'ac',1.2,'bc',1500);
MA2 = CHAMELEM.MATE(MO2,'youn',200000e6,'nu',0.3,'rho',...
7100,'sect',pi*0.012^2/4,'sigy',500e6,'K',0.01*210000e6,'ay',0.83);
MA3 = CHAMELEM.MATE(MO3,'youn',200000e6,'nu',0.3,'rho',...
7100,'sect',pi*0.012^2/4,'sigy',500e6,'K',0.01*210000e6,'ay',0.83);
MAT = [MA1 [MA2 MA3]];
%% Definition of the beam
FILE = '103_2.mail';
ME = INPUT.ACQU(FILE,'MAIL');
%% Definition of the model
MOB1 = MODEL('LPOUTRE','MECHANICS','ELASTICITY','SECTION','PLASTICITY',...
'SECTION','TIMO');
MOBT = MOB1;
%% Topology
TP = TOPOLOGY(MOBT);
%% Definition of the material
MAB1 = CHAMELEM.MATE(MOB1,'mods',MOT,'mats',MAT);
MABT = MAB1;
%% Stiffness matrix
RIG1 = MATRICE('STIFF',MOB1,MAB1);
%% Boundary conditions
CL1 = MATRICE('DIRI','P1G',1,2);
CL2 = MATRICE('DIRI','P1D',2);
CL3= MATRICE('DIRI','LPOUTRE',3,4,5);
CLT = [CL1 [CL2 CL3]];
%% Computation of the mass matrix
MAS = MATRICE('MASS',MOBT,MABT);
MAD1 = MATRICE('MADD','PVERIN2',[1 2 3],180./3);
MAD2 = MATRICE('MADD','PVERIN1',[1 2 3],180./3);
MAD = [MAD1 MAD2];
%% MODAL
PB2 = PROBLEM('model',MOBT,'mater',MABT,'diric',CLT,'solve_type','MODAL','lumped_mass',MAD);
SOL2 = SOLVE(PB2);
%%MASS_TOT
MASTO = [MAS [MAD1 MAD2]];
%% Definition of a prescribed displacement
FOU = CHPOINT('LABEL','LPOUTRE',2,-9.81/5.);
%% Equivalent force field
FO1 = MASTO * FOU;
%% Loading
EV1 = EVOL([0 0.1 0.2 0.3 20],[0 1 1 0 0],'Time','Displacement (m)');
CHT = TIMELOAD(FO1,EV1,'MECA');
%% Static analysis
PB1 = PROBLEM('model',MOBT,'mater',MABT,'diric',CLT,'loadt',CHT,'comp_time',0:0.01:2,...
'solve_type','DYNAMIC','lumped_mass',MAD);
SOL = SOLVE(PB1);
EV_OUT_D = EVOL.TIME(SOL,'displacement','P0',2);
plot(EV_OUT_D);
% Reference results (Cast3M)
ref1 = ...
[ 0.0000 0.0000
1.00000E-02 -4.65572E-06
2.00000E-02 -2.33076E-05
3.00000E-02 -5.55792E-05
4.00000E-02 -8.70562E-05
5.00000E-02 -1.04188E-04
6.00000E-02 -1.08119E-04
7.00000E-02 -1.13623E-04
8.00000E-02 -1.33805E-04
9.00000E-02 -1.66719E-04
0.10000 -1.97302E-04
0.11000 -2.08241E-04
0.12000 -1.92958E-04
0.13000 -1.67136E-04
0.14000 -1.57370E-04
0.15000 -1.73647E-04
0.16000 -1.99288E-04
0.17000 -2.07876E-04
0.18000 -1.90625E-04
0.19000 -1.65242E-04
0.20000 -1.57838E-04
0.21000 -1.71377E-04
0.22000 -1.77786E-04
0.23000 -1.51714E-04
0.24000 -1.01151E-04
0.25000 -5.93267E-05
0.26000 -5.04301E-05
0.27000 -6.48419E-05
0.28000 -6.89427E-05
0.29000 -3.97324E-05
0.30000 1.15371E-05
0.31000 4.62983E-05
0.32000 3.34478E-05
0.33000 -1.37875E-05
0.34000 -4.68495E-05
0.35000 -3.17419E-05
0.36000 1.60063E-05
0.37000 4.72901E-05
0.38000 2.99598E-05
0.39000 -1.81852E-05
0.40000 -4.76219E-05
0.41000 -2.81037E-05
0.42000 2.03178E-05
0.43000 4.78442E-05
0.44000 2.61790E-05
0.45000 -2.24011E-05
0.46000 -4.79539E-05
0.47000 -2.41931E-05
0.48000 2.44331E-05
0.49000 4.79478E-05
0.50000 2.21532E-05
0.51000 -2.64108E-05
0.52000 -4.78253E-05
0.53000 -2.00636E-05
0.54000 2.83279E-05
0.55000 4.75891E-05
0.56000 1.79264E-05
0.57000 -3.01769E-05
0.58000 -4.72426E-05
0.59000 -1.57440E-05
0.60000 3.19514E-05
0.61000 4.67878E-05
0.62000 1.35215E-05
0.63000 -3.36481E-05
0.64000 -4.62239E-05
0.65000 -1.12668E-05
0.66000 3.52662E-05
0.67000 4.55489E-05
0.68000 8.98815E-06
0.69000 -3.68042E-05
0.70000 -4.47632E-05
0.71000 -6.69136E-06
0.72000 3.82578E-05
0.73000 4.38699E-05
0.74000 4.37955E-06
0.75000 -3.96208E-05
0.76000 -4.28744E-05
0.77000 -2.05513E-06
0.78000 4.08873E-05
0.79000 4.17809E-05
0.80000 -2.77393E-07
0.81000 -4.20544E-05
0.82000 -4.05909E-05
0.83000 2.61051E-06
0.84000 4.31220E-05
0.85000 3.93041E-05
0.86000 -4.93556E-06
0.87000 -4.40904E-05
0.88000 -3.79214E-05
0.89000 7.24567E-06
0.90000 4.49577E-05
0.91000 3.64466E-05
0.92000 -9.53698E-06
0.93000 -4.57196E-05
0.94000 -3.48862E-05
0.95000 1.18072E-05
0.96000 4.63711E-05
0.97000 3.32464E-05
0.98000 -1.40529E-05
0.99000 -4.69099E-05
1.0000 -3.15308E-05
1.0100 1.62675E-05
1.0200 4.73368E-05
1.0300 2.97409E-05
1.0400 -1.84426E-05
1.0500 -4.76538E-05
1.0600 -2.78780E-05
1.0700 2.05711E-05
1.0800 4.78616E-05
1.0900 2.59463E-05
1.1000 -2.26485E-05
1.1100 -4.79580E-05
1.1200 -2.39527E-05
1.1300 2.46729E-05
1.1400 4.79395E-05
1.1500 2.19050E-05
1.1600 -2.66421E-05
1.1700 -4.78046E-05
1.1800 -1.98087E-05
1.1900 2.85511E-05
1.2000 4.75548E-05
1.2100 1.76667E-05
1.2200 -3.03927E-05
1.2300 -4.71937E-05
1.2400 -1.54810E-05
1.2500 3.21597E-05
1.2600 4.67245E-05
1.2700 1.32553E-05
1.2800 -3.38477E-05
1.2900 -4.61472E-05
1.3000 -1.09969E-05
1.3100 3.54554E-05
1.3200 4.54604E-05
1.3300 8.71411E-06
1.3400 -3.69818E-05
1.3500 -4.46632E-05
1.3600 -6.41406E-06
1.3700 3.84240E-05
1.3800 4.37578E-05
1.3900 4.10078E-06
1.4000 -3.97765E-05
1.4100 -4.27492E-05
1.4200 -1.77661E-06
1.4300 4.10328E-05
1.4400 4.16425E-05
1.4500 -5.54997E-07
1.4600 -4.21890E-05
1.4700 -4.04404E-05
1.4800 2.88760E-06
1.4900 4.32442E-05
1.5000 3.91431E-05
1.5100 -5.21264E-06
1.5200 -4.41989E-05
1.5300 -3.77509E-05
1.5400 7.52227E-06
1.5500 4.50526E-05
1.5600 3.62663E-05
1.5700 -9.81150E-06
1.5800 -4.58018E-05
1.5900 -3.46953E-05
1.6000 1.20778E-05
1.6100 4.64415E-05
1.6200 3.30446E-05
1.6300 -1.43184E-05
1.6400 -4.69683E-05
1.6500 -3.13190E-05
1.6600 1.65282E-05
1.6700 4.73819E-05
1.6800 2.95208E-05
1.6900 -1.86992E-05
1.7000 -4.76845E-05
1.7100 -2.76511E-05
1.7200 2.08233E-05
1.7300 4.78776E-05
1.7400 2.57127E-05
1.7500 -2.28952E-05
1.7600 -4.79602E-05
1.7700 -2.37119E-05
1.7800 2.49123E-05
1.7900 4.79292E-05
1.8000 2.16565E-05
1.8100 -2.68728E-05
1.8200 -4.77821E-05
1.8300 -1.95532E-05
1.8400 2.87732E-05
1.8500 4.75192E-05
1.8600 1.74059E-05
1.8700 -3.06068E-05
1.8800 -4.71440E-05
1.8900 -1.52166E-05
1.9000 3.23662E-05
1.9100 4.66602E-05
1.9200 1.29881E-05
1.9300 -3.40458E-05
1.9400 -4.60692E-05
1.9500 -1.07265E-05
1.9600 3.56435E-05
1.9700 4.53701E-05
1.9800 8.44005E-06
1.9900 -3.71584E-05
2.0000 -4.45614E-05];
hold on
plot(ref1(:,1),ref1(:,2),'go')
legend('CastLab-TIMO','Cast3M','Location','Best')
%close all
EV_OUT_V = EVOL.TIME(SOL,'velocity','P0',2);
plot(EV_OUT_V);
%reference Cast3M
ref2 = ...
[
0.0000 0.0000
1.00000E-02 -9.31144E-04
2.00000E-02 -2.79923E-03
3.00000E-02 -3.65510E-03
4.00000E-02 -2.64030E-03
5.00000E-02 -7.86043E-04
6.00000E-02 -2.48291E-07
7.00000E-02 -1.10054E-03
8.00000E-02 -2.93581E-03
9.00000E-02 -3.64691E-03
0.10000 -2.46973E-03
0.11000 2.81961E-04
0.12000 2.77452E-03
0.13000 2.38989E-03
0.14000 -4.36572E-04
0.15000 -2.81882E-03
0.16000 -2.30943E-03
0.17000 5.91766E-04
0.18000 2.85845E-03
0.19000 2.21816E-03
0.20000 -7.37408E-04
0.21000 -1.97024E-03
0.22000 6.88435E-04
0.23000 4.52593E-03
0.24000 5.58672E-03
0.25000 2.77804E-03
0.26000 -9.98717E-04
0.27000 -1.88363E-03
0.28000 1.06346E-03
0.29000 4.77861E-03
0.30000 5.47530E-03
0.31000 1.47694E-03
0.32000 -4.04705E-03
0.33000 -5.40001E-03
0.34000 -1.21238E-03
0.35000 4.23391E-03
0.36000 5.31573E-03
0.37000 9.41035E-04
0.38000 -4.40710E-03
0.39000 -5.22189E-03
0.40000 -6.65448E-04
0.41000 4.56908E-03
0.42000 5.11523E-03
0.43000 3.90040E-04
0.44000 -4.72309E-03
0.45000 -4.99292E-03
0.46000 -1.17644E-04
0.47000 4.86981E-03
0.48000 4.85543E-03
0.49000 -1.52492E-04
0.50000 -5.00644E-03
0.51000 -4.70634E-03
0.52000 4.23427E-04
0.53000 5.12893E-03
0.54000 4.54936E-03
0.55000 -6.97135E-04
0.56000 -5.23540E-03
0.57000 -4.38526E-03
0.58000 9.72115E-04
0.59000 5.32760E-03
0.60000 4.21147E-03
0.61000 -1.24418E-03
0.62000 -5.40908E-03
0.63000 -4.02485E-03
0.64000 1.50970E-03
0.65000 5.48170E-03
0.66000 3.82491E-03
0.67000 -1.76837E-03
0.68000 -5.54379E-03
0.69000 -3.61468E-03
0.70000 2.02288E-03
0.71000 5.59148E-03
0.72000 3.39836E-03
0.73000 -2.27595E-03
0.74000 -5.62211E-03
0.75000 -3.17796E-03
0.76000 2.52725E-03
0.77000 5.63660E-03
0.78000 2.95189E-03
0.79000 -2.77318E-03
0.80000 -5.63848E-03
0.81000 -2.71692E-03
0.82000 3.00961E-03
0.83000 5.63067E-03
0.84000 2.47162E-03
0.85000 -3.23518E-03
0.86000 -5.61276E-03
0.87000 -2.21820E-03
0.88000 3.45199E-03
0.89000 5.58143E-03
0.90000 1.96098E-03
0.91000 -3.66320E-03
0.92000 -5.53352E-03
0.93000 -1.70300E-03
0.94000 3.86967E-03
0.95000 5.46901E-03
0.96000 1.44376E-03
0.97000 -4.06871E-03
0.98000 -5.39114E-03
0.99000 -1.18026E-03
1.0000 4.25608E-03
1.0100 5.30359E-03
1.0200 9.10270E-04
1.0300 -4.42944E-03
1.0400 -5.20727E-03
1.0500 -6.34962E-04
1.0600 4.59012E-03
1.0700 5.09971E-03
1.0800 3.58406E-04
1.0900 -4.74147E-03
1.1000 -4.97748E-03
1.1100 -8.44243E-05
1.1200 4.88548E-03
1.1300 4.83964E-03
1.1400 -1.86312E-04
1.1500 -5.02060E-03
1.1600 -4.68882E-03
1.1700 4.56314E-04
1.1800 5.14286E-03
1.1900 4.52911E-03
1.2000 -7.28382E-04
1.2100 -5.24922E-03
1.2200 -4.36266E-03
1.2300 1.00246E-03
1.2400 5.34010E-03
1.2500 4.18803E-03
1.2600 -1.27508E-03
1.2700 -5.41875E-03
1.2800 -4.00186E-03
1.2900 1.54196E-03
1.3000 5.48811E-03
1.3100 3.80235E-03
1.3200 -1.80134E-03
1.3300 -5.54792E-03
1.3400 -3.59126E-03
1.3500 2.05499E-03
1.3600 5.59484E-03
1.3700 3.37278E-03
1.3800 -2.30604E-03
1.3900 -5.62536E-03
1.4000 -3.15009E-03
1.4100 2.55555E-03
1.4200 5.63896E-03
1.4300 2.92293E-03
1.4400 -2.80100E-03
1.4500 -5.63849E-03
1.4600 -2.68832E-03
1.4700 3.03805E-03
1.4800 5.62755E-03
1.4900 2.44377E-03
1.5000 -3.26398E-03
1.5100 -5.60717E-03
1.5200 -2.19009E-03
1.5300 3.47970E-03
1.5400 5.57493E-03
1.5500 1.93115E-03
1.5600 -3.68841E-03
1.5700 -5.52716E-03
1.5800 -1.67090E-03
1.5900 3.89220E-03
1.6000 5.46241E-03
1.6100 1.41033E-03
1.6200 -4.08972E-03
1.6300 -5.38288E-03
1.6400 -1.14708E-03
1.6500 4.27693E-03
1.6600 5.29251E-03
1.6700 8.78229E-04
1.6800 -4.45044E-03
1.6900 -5.19355E-03
1.7000 -6.03512E-04
1.7100 4.61019E-03
1.7200 5.08469E-03
1.7300 3.26160E-04
1.7400 -4.75913E-03
1.7500 -4.96246E-03
1.7600 -5.05363E-05
1.7700 4.90019E-03
1.7800 4.82465E-03
1.7900 -2.21266E-04
1.8000 -5.03327E-03
1.8100 -4.67259E-03
1.8200 4.90746E-04
1.8300 5.15502E-03
1.8400 4.51025E-03
1.8500 -7.61033E-04
1.8600 -5.26163E-03
1.8700 -4.34090E-03
1.8800 1.03344E-03
1.8900 5.35205E-03
1.9000 4.16449E-03
1.9100 -1.30568E-03
1.9200 -5.42873E-03
1.9300 -3.97806E-03
1.9400 1.57338E-03
1.9500 5.49516E-03
1.9600 3.77883E-03
1.9700 -1.83350E-03
1.9800 -5.55251E-03
1.9900 -3.56718E-03
2.0000 2.08657E-03];
hold on
plot(ref2(:,1),ref2(:,2),'go')
legend('CastLab-TIMO','Cast3M','Location','Best')
EV_OUT_A = EVOL.TIME(SOL,'acceleration','P0',2);
plot(EV_OUT_A);
%reference Cast3M
ref3 = ...
[ 0.0000 0.0000
1.00000E-02 -0.18623
2.00000E-02 -0.18739
3.00000E-02 1.62138E-02
4.00000E-02 0.18675
5.00000E-02 0.18410
6.00000E-02 -2.69455E-02
7.00000E-02 -0.19311
8.00000E-02 -0.17394
9.00000E-02 3.17214E-02
0.10000 0.20371
0.11000 0.34663
0.12000 0.15189
0.13000 -0.22881
0.14000 -0.33648
0.15000 -0.13997
0.16000 0.24185
0.17000 0.33839
0.18000 0.11494
0.19000 -0.24300
0.20000 -0.34811
0.21000 0.10155
0.22000 0.43019
0.23000 0.33731
0.24000 -0.12515
0.25000 -0.43659
0.26000 -0.31876
0.27000 0.14178
0.28000 0.44763
0.29000 0.29540
0.30000 -0.15606
0.31000 -0.64361
0.32000 -0.46118
0.33000 0.19059
0.34000 0.64693
0.35000 0.44233
0.36000 -0.22596
0.37000 -0.64897
0.38000 -0.42065
0.39000 0.25769
0.40000 0.65359
0.41000 0.39331
0.42000 -0.28408
0.43000 -0.66096
0.44000 -0.36167
0.45000 0.30770
0.46000 0.66735
0.47000 0.33014
0.48000 -0.33301
0.49000 -0.66857
0.50000 -0.30222
0.51000 0.36224
0.52000 0.66372
0.53000 0.27739
0.54000 -0.39330
0.55000 -0.65600
0.56000 -0.25165
0.57000 0.42168
0.58000 0.64979
0.59000 0.22130
0.60000 -0.44453
0.61000 -0.64660
0.62000 -0.18638
0.63000 0.46322
0.64000 0.64369
0.65000 0.15071
0.66000 -0.48207
0.67000 -0.63658
0.68000 -0.11850
0.69000 0.50433
0.70000 0.62319
0.71000 9.05330E-02
0.72000 -0.52916
0.73000 -0.60571
0.74000 -6.35254E-02
0.75000 0.55236
0.76000 0.58869
0.77000 3.31827E-02
0.78000 -0.57012
0.79000 -0.57489
0.80000 1.82795E-03
0.81000 0.58248
0.82000 0.56282
0.83000 -3.86115E-02
0.84000 -0.59320
0.85000 -0.54816
0.86000 7.26495E-02
0.87000 0.60626
0.88000 0.52778
0.89000 -0.10189
0.90000 -0.62220
0.91000 -0.50263
0.92000 0.12857
0.93000 0.63754
0.94000 0.47700
0.95000 -0.15713
0.96000 -0.64792
0.97000 -0.45458
0.98000 0.19009
0.99000 0.65209
1.0000 0.43518
1.0100 -0.22568
1.0200 -0.65298
1.0300 -0.41496
1.0400 0.25939
1.0500 0.65507
1.0600 0.38994
1.0700 -0.28803
1.0800 -0.66023
1.0900 -0.35974
1.1000 0.31254
1.1100 0.66607
1.1200 0.32791
1.1300 -0.33708
1.1400 -0.66811
1.1500 -0.29874
1.1600 0.36510
1.1700 0.66393
1.1800 0.27338
1.1900 -0.39613
1.2000 -0.65536
1.2100 -0.24880
1.2200 0.42611
1.2300 0.64691
1.2400 0.22062
1.2500 -0.45103
1.2600 -0.64159
1.2700 -0.18714
1.2800 0.47052
1.2900 0.63824
1.3000 0.15098
1.3100 -0.48814
1.3200 -0.63260
1.3300 -0.11671
1.3400 0.50804
1.3500 0.62121
1.3600 8.67640E-02
1.3700 -0.53117
1.3800 -0.60459
1.3900 -5.92766E-02
1.4000 0.55433
1.4100 0.58680
1.4200 2.98861E-02
1.4300 -0.57309
1.4400 -0.57169
1.4500 4.19482E-03
1.4600 0.58584
1.4700 0.55943
1.4800 -4.15328E-02
1.4900 -0.59522
1.5000 -0.54633
1.5100 7.76891E-02
1.5200 0.60573
1.5300 0.52823
1.5400 -0.10918
1.5500 -0.61957
1.5600 -0.50434
1.5700 0.13659
1.5800 0.63466
1.5900 0.47796
1.6000 -0.16392
1.6100 -0.64650
1.6200 -0.45351
1.6300 0.19488
1.6400 0.65228
1.6500 0.43252
1.6600 -0.22940
1.6700 -0.65346
1.6800 -0.41228
1.6900 0.26366
1.7000 0.65435
1.7100 0.38839
1.7200 -0.29349
1.7300 -0.65822
1.7400 -0.35884
1.7500 0.31817
1.7600 0.66421
1.7700 0.32593
1.7800 -0.34104
1.7900 -0.66814
1.8000 -0.29426
1.8100 0.36640
1.8200 0.66627
1.8300 0.26658
1.8400 -0.39554
1.8500 -0.65872
1.8600 -0.24140
1.8700 0.42555
1.8800 0.64932
1.8900 0.21440
1.9000 -0.45192
1.9100 -0.64212
1.9200 -0.18249
1.9300 0.47263
1.9400 0.63766
1.9500 0.14670
1.9600 -0.48997
1.9700 -0.63250
1.9800 -0.11130
1.9900 0.50837
2.0000 0.62238];
hold on
plot(ref3(:,1),ref3(:,2),'go')
legend('CastLab-TIMO','Cast3M','Location','Best')
%% Non regression test
if abs((EV_OUT_D.ordo(end) - ref1(end,2))/ref1(end,2)) > 1.0e-2||...
abs((EV_OUT_V.ordo(end) - ref2(end,2))/ref2(end,2)) > 4.0e-2||...
abs((EV_OUT_A.ordo(end) - ref3(end,2))/ref3(end,2)) > 1.0e-2||...
abs((SOL2.eigenfrequency(1) - 18.723)/18.723) > 1.0e-3 ||...
abs((SOL2.eigenfrequency(2) - 73.305)/73.305)> 1.0e-3 ||...
abs((SOL2.eigenfrequency(3) - 189.77)/189.77) > 1.0e-3
error('TEST IS NOT SUCCESSFUL')
else
disp('---------------------------------')
disp('0103_INCR_DYNAMIC_MULTIFIBER_TIMO')
disp('TEST IS SUCCESSFUL')
disp('---------------------------------')
end