% AbaqusPlot
%-----------------------------------------------------------------------
% An example of plotting data from Abaqus. The data are taken from the
% Abaqus-tutorial "AbaqusSimpleImperfektionLatTorsBuckling_2010.pdf"
% which is found on the Abaqus-page.
%-----------------------------------------------------------------------
% Johan Clausen
% Department of Civil Engineering
% Ålborg University
% October 2010
%-----------------------------------------------------------------------
clear all

end_moment = 19000 ; % The moment that is applied in Abaqus
lt_moment = 15663 ; % The lateral torsional buckling moment from "AbaqusSimpleLatTorsBuckling_2010.pdf"
%----- Data From Abaqus ---------
% Linear case, e = 0.1
Lin_e01 = [
0.05	0.00309835
0.1		0.0061967
0.175	0.0108442
0.275	0.0170409
0.375	0.0232376
0.475	0.0294343
0.575	0.035631
0.675	0.0418277
0.775	0.0480244
0.875	0.0542211
0.975	0.0604178
1	0.061967 ] ;

% Non-linear case, e = 0.1
Nlin_e01 = [
0	0
0.05	0.00343593
0.1	0.00760954
0.175	0.0155215
0.275	0.0301202
0.375	0.0515239
0.475	0.0836734
0.5	0.0940187
0.525	0.105488
0.5625	0.124911
0.576563	0.132898
0.597656	0.145573
0.629297	0.166023
0.641162	0.174068
0.65896	0.186449
0.685657	0.205528
0.725702	0.234561
0.765747	0.262963
0.805792	0.290194
0.845837	0.315555
0.885883	0.338699
0.925928	0.359344
0.965973	0.377168
1	0.389824 ] ;

% Non-linear case, e = 0.05
Nlin_e005 = [
0	0
0.05	0.00172565
0.1	0.00382762
0.175	0.00783285
0.275	0.0153019
0.375	0.0265026
0.475	0.0442188
0.575	0.0752461
0.6	0.0864575
0.625	0.0997907
0.65	0.115402
0.675	0.133691
0.7125	0.166425
0.75	0.202923
0.7875	0.241045
0.825	0.278613
0.8625	0.312781
0.9	0.34304
0.9375	0.369649
0.975	0.392594
0.98125	0.396044
0.9875	0.399383
0.996875	0.404187
1	0.405732 ] ;
% Non-linear case, e = 0.03
Nlin_e003 = [
0	0
0.05	0.00103601
0.1	0.00229882
0.175	0.0047066
0.2875	0.00990418
0.45625	0.0242974
0.519531	0.0340343
0.614453	0.0589126
0.638184	0.0685462
0.661914	0.0806096
0.685645	0.0955384
0.709375	0.113989
0.733105	0.136972
0.756836	0.164015
0.762769	0.171136
0.768701	0.178536
0.7776	0.189857
0.790949	0.207124
0.810971	0.233295
0.815977	0.239253
0.820982	0.245518
0.828491	0.254782
0.839754	0.268305
0.856648	0.287661
0.873542	0.305862
0.890436	0.322875
0.90733	0.338767
0.924224	0.353585
0.941118	0.367371
0.958012	0.38016
0.974906	0.391977
0.9918	0.402829
1	0.40775 ] ;

% % Non-linear case, e = 0.02
% Nlin_e002 = [
% 0	0
% 0.05	0.000690737
% 0.1	0.00153286
% 0.175	0.0031388
% 0.2875	0.00660781
% 0.45625	0.016263
% 0.519531	0.0228083
% 0.614453	0.0396436
% 0.650049	0.050619
% 0.703442	0.077048
% 0.756836	0.126879
% 0.76351	0.135201
% 0.770184	0.144189
% 0.780196	0.158416
% 0.795213	0.180975
% 0.810229	0.204303
% 0.813984	0.209844
% 0.817738	0.215547
% 0.823369	0.224038
% 0.831816	0.236524
% 0.844487	0.254569
% 0.857157	0.271694
% 0.869828	0.28785
% 0.882499	0.303063
% 0.895169	0.317376
% 0.90784	0.330842
% 0.92051	0.343515
% 0.933181	0.355438
% 0.945851	0.366656
% 0.958522	0.377206
% 0.971192	0.387108
% 0.990198	0.400794
% 1	0.407318 ] ;

% Non-linear case, e = 0.01
Nlin_e001 = [
0	0
0.05	0.000345382
0.1	0.000766508
0.175	0.00156958
0.2875	0.00330522
0.45625	0.00815555
0.709375	0.0424053
0.725195	0.0499267
0.741016	0.0598032
0.756836	0.0730728
0.772656	0.0904643
0.788477	0.114152
0.794409	0.125642
0.800342	0.136827
0.806274	0.148497
0.812207	0.160469
0.821106	0.178683
0.824443	0.185313
0.829449	0.195308
0.836957	0.209878
0.844466	0.223895
0.851974	0.237287
0.859483	0.250041
0.866991	0.262175
0.874499	0.273709
0.882008	0.284681
0.889516	0.295132
0.897025	0.305089
0.908288	0.319186
0.91955	0.332353
0.930813	0.344692
0.942076	0.356266
0.953338	0.367141
0.964601	0.377366
0.975864	0.386975
0.992758	0.400301
1	0.405632 ] ;
%---------------------------------

%----- Plotting -------------------------
figure(1) ; clf ; hold on ;
plot([0 0.39],[lt_moment lt_moment],'k--') ; % Buckling load of the perfectly straight beam
plot(Lin_e01(:,2),end_moment*Lin_e01(:,1),'k--','linewidth',1.5) ;
plot(Nlin_e01(:,2),end_moment*Nlin_e01(:,1),'b-','linewidth',1.5) ;
plot(Nlin_e005(:,2),end_moment*Nlin_e005(:,1),'c-','linewidth',1.5) ;
plot(Nlin_e003(:,2),end_moment*Nlin_e003(:,1),'r-','linewidth',1.5) ;
plot(Nlin_e001(:,2),end_moment*Nlin_e001(:,1),'m-','linewidth',1.5) ;
legend('Bifurcation buckling','Linear, {\it e} = 0.1','Non-linear, {\it e} = 0.1','Non-linear, {\it e} = 0.05','Non-linear, {\it e} = 0.03','Non-linear, {\it e} = 0.01','Location','Best')
xlabel('{\it y}-displacement','fontname','times','fontsize',13)
ylabel('Moment','fontname','times','fontsize',13)
%----------------------------------------