Hello
I want to use the ABAQUS2MATLAB tool for analysis in Abaqus / Explicit mode. The .inp file uses Abaqus / CAE to generate, I also run the analysis in Abaqus / Explicit without error. But I use Matlab 2015a and Abqus2matlab to analyze the following error:
>> Script_tn4_A2M
Simulation Started
Simulation Finished
Obtaining desired output data by Abaqus2Matlab
Abaqus Error: No input file was specified.The Abaqus analysis cannot be executed.
Abaqus ascfil execution exited with error(s).
NO NEED TO UPGRADE. THE OUTPUT DATABASE FILE IS ALREADY UP-TO-DATE
Abaqus License Manager checked out the following license(s):
"cae" release 6.14 from Flexnet server localhost
<1022 out of 1024 licenses remain available>.
NO NEED TO UPGRADE. THE OUTPUT DATABASE FILE IS ALREADY UP-TO-DATE
Abaqus License Manager checked out the following license(s):
"cae" release 6.14 from Flexnet server localhost
<1022 out of 1024 licenses remain available>.
Warning: The specified element set was not found in the output database tn4_A2M.odb. All element sets will be processed
NO NEED TO UPGRADE. THE OUTPUT DATABASE FILE IS ALREADY UP-TO-DATE
Abaqus License Manager checked out the following license(s):
"cae" release 6.14 from Flexnet server localhost
<1022 out of 1024 licenses remain available>.
Warning: The specified node set was not found in the output database tn4_A2M.odb. All node sets will be processed
NO NEED TO UPGRADE. THE OUTPUT DATABASE FILE IS ALREADY UP-TO-DATE
Abaqus License Manager checked out the following license(s):
"cae" release 6.14 from Flexnet server localhost
<1022 out of 1024 licenses remain available>.
Error using readFil (line 75)
The file tn4_A2M.fil does not exist
Error in Script_tn4_A2M (line 87)
S=readFil([ Inp_file '.fil' ],11);%Obtain the Stress(S);
>>
And here is the code from the Abaqus2Matlab I used to run
% If using this code for research or industrial purposes please cite:
% G. Papazafeiropoulos, M. Muniz-Calvente, E. Martinez-Paneda
% Abaqus2Matlab: a suitable tool for finite element post-processing.
% Advances in Engineering Software. Vol. 105, pp 9-16 (2017)
% DOI:10.1016/j.advengsoft.2017.01.006
% % G. Papazafeiropoulos, M. Muniz Calvente, E. Martinez-Paneda
% % % Abaqus2Matlab@gmail.com
% www.abaqus2matlab.com
%close all
clear
% Change the current directory
S = mfilename('fullpath');
f = filesep;
ind=strfind(S,f);
S1=S(1:ind(end)-1);
cd(S1) ;
%% 1st STEP - Run one FEM model
Inp_file= 'tn4_A2M';
disp('Simulation Started')
% Run the input file with Abaqus
% Initialize sw (boolean switch) as true
system(['abaqus job=' Inp_file]);
sw=true;
tic;
while sw
% Pause Matlab execution in order for the lck file to be created
pause(0.5);
% While the lck file exists, pause Matlab execution. If it is
% deleted, exit the while loop and proceed.
while exist([Inp_file '.lck'],'file')==2
pause(0.1)
% the lck file has been created and Matlab halts in this loop.
% Set sw to false to break the outer while loop and continue
% the code execution.
sw=false;
end
% In case that the lck file cannot be detected, then terminate
% infinite execution of the outer while loop after a certain
% execution time limit (5 sec)
if sw && (toc>5)
sw=false;
end
end
% NOTE: Alternatively, you can replace lines 27 to 49 by system(['abaqus job=' Inp_file ' interactive'])
disp('Simulation Finished')
%% 2st STEP - Postprocess Abaqus results file with Abaqu2Matlab
% Obtain the desired output data
disp('Obtaining desired output data by Abaqus2Matlab')
% NOTE: Some output variables are Matlab Cells, If you are not confortable working with Cells, you can use cell2mat()
% Convert fil file to fin (ASCII format)
!abaqus ascfil job=tn4_A2M
% Give Abaqus enough time for the file conversion
pause(2)
% Explore the field output contained in AbaqusInputFile.odb
odb_name='tn4_A2M.odb';
out=exploreFieldOdb(odb_name);
% You have selected "output,field" but any specific variable related to elements has been selected.
% It is an example of how to read from Odb the element outputs (Stresses).
odb_name='tn4_A2M.odb';
stepName='Step-1';
elsetName= 'False' ;
instanceName= 'False' ;
indOut='False'; %(no dat file)
indOut='True'; %(create dat file)
outputVar='S';% Stresses
[odbOut_S,odbDat_S,rpyOut_S]=readElementFieldOdb(odb_name,stepName,instanceName,elsetName,outputVar,indOut);
% You have selected "output,field" but any specific variable related to nodes has been selected.
% It is an example of how to read from Odb the node outputs (Coordinates).
odb_name='tn4_A2M.odb';
stepName='Step-1';
nsetName= 'False' ;
instanceName= 'False' ;
indOut='False'; %(no dat file)
indOut='True'; %(create dat file)
outputVar='COORD';%Coordinates
[odbOut_COORD,odbDat_COORD,rpyOut_COORD]=readNodeFieldOdb(odb_name,stepName,instanceName,nsetName,outputVar,indOut);
% Explore the history output contained in AbaqusInputFile.odb
odb_name='tn4_A2M.odb';
out=exploreHistoryOdb(odb_name);
S=readFil([ Inp_file '.fil' ],11);%Obtain the Stress(S);
SINV=readFil([ Inp_file '.fil' ],12);%Obtain the Stress Invariant(SINV);
This is the .inp file I try to analyze
*Heading
*Preprint, echo=NO, model=NO, history=NO, contact=NO
*Part, name=PART-1-1
*Node
1, 80., 0.
2, 80., 5.4545455
3, 80., 10.909091
4, 80., 16.363636
5, 80., 21.818182
6, 80., 27.272728
7, 80., 32.727272
8, 80., 38.1818199
9, 80., 43.636364
10, 80., 49.0909081
11, 80., 54.5454559
12, 80., 60.
13, 55., 0.
14, 55., 5.4545455
15, 55., 10.909091
16, 55., 16.363636
17, 55., 21.818182
18, 55., 27.272728
19, 55., 32.727272
20, 55., 38.1818199
21, 55., 43.636364
22, 55., 49.0909081
23, 55., 54.5454559
24, 55., 60.
25, 50., 0.
26, 50., 5.4545455
27, 50., 10.909091
28, 50., 16.363636
29, 50., 21.818182
30, 50., 27.272728
31, 50., 32.727272
32, 50., 38.1818199
33, 50., 43.636364
34, 50., 49.0909081
35, 50., 54.5454559
36, 50., 60.
37, 45., 0.
38, 45., 5.4545455
39, 45., 10.909091
40, 45., 16.363636
41, 45., 21.818182
42, 45., 27.272728
43, 45., 32.727272
44, 45., 38.1818199
45, 45., 43.636364
46, 45., 49.0909081
47, 45., 54.5454559
48, 45., 60.
49, 40., 0.
50, 40., 5.4545455
51, 40., 10.909091
52, 40., 16.363636
53, 40., 21.818182
54, 40., 27.272728
55, 40., 32.727272
56, 40., 38.1818199
57, 40., 43.636364
58, 40., 49.0909081
59, 40., 54.5454559
60, 40., 60.
61, 35., 0.
62, 35., 5.4545455
63, 35., 10.909091
64, 35., 16.363636
65, 35., 21.818182
66, 35., 27.272728
67, 35., 32.727272
68, 35., 38.1818199
69, 35., 43.636364
70, 35., 49.0909081
71, 35., 54.5454559
72, 35., 60.
73, 30., 0.
74, 30., 5.4545455
75, 30., 10.909091
76, 30., 16.363636
77, 30., 21.818182
78, 30., 27.272728
79, 30., 32.727272
80, 30., 38.1818199
81, 30., 43.636364
82, 30., 49.0909081
83, 30., 54.5454559
84, 30., 60.
85, 25., 0.
86, 25., 5.4545455
87, 25., 10.909091
88, 25., 16.363636
89, 25., 21.818182
90, 25., 27.272728
91, 25., 32.727272
92, 25., 38.1818199
93, 25., 43.636364
94, 25., 49.0909081
95, 25., 54.5454559
96, 25., 60.
97, 19.7261105, 0.
98, 19.751009, 5.4545455
99, 19.7759094, 10.909091
100, 19.800808, 16.363636
101, 19.8257065, 21.818182
102, 19.850605, 27.272728
103, 19.8755054, 32.727272
104, 19.900404, 38.1818199
105, 19.9253025, 43.636364
106, 19.950201, 49.0909081
107, 19.9751015, 54.5454559
108, 20., 60.
109, 14.7945824, 0.
110, 14.8132572, 5.4545455
111, 14.8319311, 10.909091
112, 14.850606, 16.363636
113, 14.8692799, 21.818182
114, 14.8879547, 27.272728
115, 14.9066286, 32.727272
116, 14.9253025, 38.1818199
117, 14.9439774, 43.636364
118, 14.9626513, 49.0909081
119, 14.9813261, 54.5454559
120, 15., 60.
121, 9.86305523, 0.
122, 9.87550449, 5.4545455
123, 9.88795471, 10.909091
124, 9.90040398, 16.363636
125, 9.91285324, 21.818182
126, 9.92530251, 27.272728
127, 9.93775272, 32.727272
128, 9.95020199, 38.1818199
129, 9.96265125, 43.636364
130, 9.97510052, 49.0909081
131, 9.98755074, 54.5454559
132, 10., 60.
133, 4.93152761, 0.
134, 4.93775225, 5.4545455
135, 4.94397736, 10.909091
136, 4.95020199, 16.363636
137, 4.95642662, 21.818182
138, 4.96265125, 27.272728
139, 4.96887636, 32.727272
140, 4.97510099, 38.1818199
141, 4.98132563, 43.636364
142, 4.98755026, 49.0909081
143, 4.99377537, 54.5454559
144, 5., 60.
145, 0., 0.
146, 0., 5.4545455
147, 0., 10.909091
148, 0., 16.363636
149, 0., 21.818182
150, 0., 27.272728
151, 0., 32.727272
152, 0., 38.1818199
153, 0., 43.636364
154, 0., 49.0909081
155, 0., 54.5454559
156, 0., 60.
*Element, type=CPS4R
1, 1, 2, 14, 13
2, 2, 3, 15, 14
3, 3, 4, 16, 15
4, 4, 5, 17, 16
5, 5, 6, 18, 17
6, 6, 7, 19, 18
7, 7, 8, 20, 19
8, 8, 9, 21, 20
9, 9, 10, 22, 21
10, 10, 11, 23, 22
11, 11, 12, 24, 23
12, 13, 14, 26, 25
13, 14, 15, 27, 26
14, 15, 16, 28, 27
15, 16, 17, 29, 28
16, 17, 18, 30, 29
17, 18, 19, 31, 30
18, 19, 20, 32, 31
19, 20, 21, 33, 32
20, 21, 22, 34, 33
21, 22, 23, 35, 34
22, 23, 24, 36, 35
23, 25, 26, 38, 37
24, 26, 27, 39, 38
25, 27, 28, 40, 39
26, 28, 29, 41, 40
27, 29, 30, 42, 41
28, 30, 31, 43, 42
29, 31, 32, 44, 43
30, 32, 33, 45, 44
31, 33, 34, 46, 45
32, 34, 35, 47, 46
33, 35, 36, 48, 47
34, 37, 38, 50, 49
35, 38, 39, 51, 50
36, 39, 40, 52, 51
37, 40, 41, 53, 52
38, 41, 42, 54, 53
39, 42, 43, 55, 54
40, 43, 44, 56, 55
41, 44, 45, 57, 56
42, 45, 46, 58, 57
43, 46, 47, 59, 58
44, 47, 48, 60, 59
45, 49, 50, 62, 61
46, 50, 51, 63, 62
47, 51, 52, 64, 63
48, 52, 53, 65, 64
49, 53, 54, 66, 65
50, 54, 55, 67, 66
51, 55, 56, 68, 67
52, 56, 57, 69, 68
53, 57, 58, 70, 69
54, 58, 59, 71, 70
55, 59, 60, 72, 71
56, 61, 62, 74, 73
57, 62, 63, 75, 74
58, 63, 64, 76, 75
59, 64, 65, 77, 76
60, 65, 66, 78, 77
61, 66, 67, 79, 78
62, 67, 68, 80, 79
63, 68, 69, 81, 80
64, 69, 70, 82, 81
65, 70, 71, 83, 82
66, 71, 72, 84, 83
67, 73, 74, 86, 85
68, 74, 75, 87, 86
69, 75, 76, 88, 87
70, 76, 77, 89, 88
71, 77, 78, 90, 89
72, 78, 79, 91, 90
73, 79, 80, 92, 91
74, 80, 81, 93, 92
75, 81, 82, 94, 93
76, 82, 83, 95, 94
77, 83, 84, 96, 95
78, 85, 86, 98, 97
79, 86, 87, 99, 98
80, 87, 88, 100, 99
81, 88, 89, 101, 100
82, 89, 90, 102, 101
83, 90, 91, 103, 102
84, 91, 92, 104, 103
85, 92, 93, 105, 104
86, 93, 94, 106, 105
87, 94, 95, 107, 106
88, 95, 96, 108, 107
89, 97, 98, 110, 109
90, 98, 99, 111, 110
91, 99, 100, 112, 111
92, 100, 101, 113, 112
93, 101, 102, 114, 113
94, 102, 103, 115, 114
95, 103, 104, 116, 115
96, 104, 105, 117, 116
97, 105, 106, 118, 117
98, 106, 107, 119, 118
99, 107, 108, 120, 119
100, 109, 110, 122, 121
101, 110, 111, 123, 122
102, 111, 112, 124, 123
103, 112, 113, 125, 124
104, 113, 114, 126, 125
105, 114, 115, 127, 126
106, 115, 116, 128, 127
107, 116, 117, 129, 128
108, 117, 118, 130, 129
109, 118, 119, 131, 130
110, 119, 120, 132, 131
111, 121, 122, 134, 133
112, 122, 123, 135, 134
113, 123, 124, 136, 135
114, 124, 125, 137, 136
115, 125, 126, 138, 137
116, 126, 127, 139, 138
117, 127, 128, 140, 139
118, 128, 129, 141, 140
119, 129, 130, 142, 141
120, 130, 131, 143, 142
121, 131, 132, 144, 143
122, 133, 134, 146, 145
123, 134, 135, 147, 146
124, 135, 136, 148, 147
125, 136, 137, 149, 148
126, 137, 138, 150, 149
127, 138, 139, 151, 150
128, 139, 140, 152, 151
129, 140, 141, 153, 152
130, 141, 142, 154, 153
131, 142, 143, 155, 154
132, 143, 144, 156, 155
*Nset, nset=SET-1, generate
1, 156, 1
*Elset, elset=SET-1, generate
1, 132, 1
*Solid Section, elset=SET-1, material=MATERIAL-1
,
*End Part
*Assembly, name=Assembly
*Instance, name=PART-1-1, part=PART-1-1
*End Instance
*Nset, nset=SET-1, instance=PART-1-1, generate
1, 85, 12
*Elset, elset=SET-1, instance=PART-1-1, generate
1, 67, 11
*Nset, nset=SET-2, instance=PART-1-1
1,
*Nset, nset=SET-5, instance=PART-1-1, generate
12, 156, 12
*Elset, elset=SET-5, instance=PART-1-1, generate
11, 132, 11
*Nset, nset=SET-6, instance=PART-1-1
12, 156
*Elset, elset=_SURF-1_S2, internal, instance=PART-1-1, generate
11, 132, 11
*Surface, type=ELEMENT, name=SURF-1
_SURF-1_S2, S2
*End Assembly
*Material, name=MATERIAL-1
*Density
7850.,
*Elastic
200000., 0.3
*Boundary
SET-1, 2, 2
*Boundary
SET-2, 1, 1
*Boundary
SET-2, 2, 2
*Step, name=Step-1, nlgeom=YES
*Dynamic, Explicit
, 1.
*Bulk Viscosity
0.06, 1.2
*Dsload
SURF-1, P, -10000.
*Restart, write, number interval=1, time marks=NO
*FILE OUTPUT
*Output,field,variable=PRESELECT
*Output,history,variable=PRESELECT
*EL FILE
S,SINV
*End Step
Please help me to change where needed in the inp file so that it can run without error. thank you very much