-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmOP2block_OQG.m
181 lines (164 loc) · 6.31 KB
/
mOP2block_OQG.m
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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
classdef mOP2block_OQG < mOP2block
%OP2BLOCK_OQG Reads table of single or multipoint constraint forces
% OP2BLOCK_OQG reads the table of single or multipoint constraint forces
% from a NASTRAN binary output 2 file.
%
% BLOCKINFO is the structure describing the block.
% DATA is a structure, containing the following fields:
%
% DATA(I), with I = 1:numel(DATA)
% |- ACODE, Device code + 10*approach code
% |- TCODE, Table code
% |- SUBCASE, Subcase or random identification number
% |- MODE*, Mode number
% |- EIGN*, Eigenvalue
% |- MODECYCL*, Mode or cycle
% |- TIMESTEP*, Time step
% |- LSDVMN*, Load set number
% |- LFTSFQ*, Load step
% |- RCODE, Random code
% |- FCODE, Format code
% |- NUMWDE, Number of words per entry in DATA record
% |- THERMAL, =1 for heat transfer and 0 otherwise
% |- TITLE, Title
% |- SUBTITL, Subtitle
% |- LABEL, Label
% |- num_entries, Number of entries in the table
% |- ElementID, Element ID for each entry
% |- Value, Output values for each entry (refer to Programmer's Guide)
%
% Note: fields marked with (*) are present only in the appropriate cases.
%
% Block described in MD Nastran 2006 DMAP Programmer's Guide, page 605.
% See also MOP2READ, MOP2INFO, MOP2BLOCK, MOP2BLOCK_*.
% $Revision: 1.1.0 $ $Date: 2022/09/08 $
%% Disclaimer
% Copyright (c) 2022 Giulio Molinari
%
% This file is part of mOP2.
% mOP2 is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% any later version. Also see the file "LICENSE".
properties
end % properties
properties(SetAccess=private)
end % properties(SetAccess=protected)
properties(Access=private)
data_copy
data_set = false
end % properties(Access=private)
properties(SetAccess=private, Dependent = true)
data
end % properties(SetAccess=private, Dependent = true)
methods
function block = mOP2block_OQG(blockInfo, fid, tableFmt)
block = block@mOP2block(blockInfo, fid, tableFmt);
block.interpreter = 'OQG';
end % mOP2block_OQG
function data = get.data(block)
if block.data_set == true
data = block.data_copy;
return
end
fid = block.fid;
blockInfo = block.blockInfo;
if isempty(strmatch('OQG', blockInfo.name)) && isempty(strmatch('OPG', blockInfo.name)) && isempty(strmatch('OQM', blockInfo.name))
error('mOP2block_OQG: Wrong block type')
end
for recordID = 1:(length(blockInfo.record)-1)/2
%% RECORD 1: IDENT
fseek(fid, blockInfo.record(recordID*2 - 1).TablesStart + 4, -1);
data(recordID).ACODE = fread(fid, 1, 'int32');
data(recordID).TCODE = fread(fid, 1, 'int32');
fseek(fid, 4, 0);
data(recordID).SUBCASE = fread(fid, 1, 'int32');
switch mOP2block.functionCodes(data(recordID).TCODE, 1)
case 1
switch mOP2block.functionCodes(data(recordID).ACODE, 4)
case {1 3 4 7 11} % Statics, PreBuckling, old NLStatic
data(recordID).LSDVMN = fread(fid, 1, 'int32');
fseek(fid, 2*4, 0);
case {2 9} % Real and Complex Eigenvalues
data(recordID).MODE = fread(fid, 1, 'int32');
data(recordID).EIGN = fread(fid, 1, 'float32');
data(recordID).MODECYCL = fread(fid, 1, 'float32');
case {8} % PostBuckling
data(recordID).MODE = fread(fid, 1, 'int32');
data(recordID).EIGN = fread(fid, 1, 'float32');
fseek(fid, 4, 0);
case {10} % Nonlinear statics (e.g. SOL400)
data(recordID).LFTSFQ = fread(fid, 1, 'int32');
fseek(fid, 2*4, 0);
case {6} % Transient
data(recordID).TIMESTEP = fread(fid, 1, 'float32');
fseek(fid, 2*4, 0);
otherwise
warning('mOP2block_OQG: ACODE %d not implemented', data(recordID).ACODE)
fseek(fid, 3*4, 0);
end
case 2
data(recordID).LSDVMN = fread(fid, 1, 'int32');
fseek(fid, 2*4, 0);
end
data(recordID).RCODE = fread(fid, 1, 'int32');
data(recordID).FCODE = fread(fid, 1, 'int32');
data(recordID).NUMWDE = fread(fid, 1, 'int32');
fseek(fid, 12*4, 0);
data(recordID).THERMAL = fread(fid, 1, 'int32');
fseek(fid, 27*4, 0);
data(recordID).TITLE = sprintf('%c', fread(fid, 32*4, 'char'));
data(recordID).SUBTITL = sprintf('%c', fread(fid, 32*4, 'char'));
data(recordID).LABEL = sprintf('%c', fread(fid, 32*4, 'char'));
%% RECORD 2: DATA
fseek(fid, blockInfo.record(recordID*2).TablesStart(1), -1);
% data(recordID).num_nodes = fread(fid, 1, 'int32') / (data(recordID).NUMWDE*4);
fseek(fid, 1*4, 0);
TablesLength = blockInfo.record(recordID*2).TablesLength;
data(recordID).num_nodes = sum(TablesLength) / (data(recordID).NUMWDE);
data(recordID).num_entries = data(recordID).num_nodes;
% Read all data
FloatData = block.readMatrix(blockInfo.record(recordID*2), data(recordID), 'float32');
IntData = block.readMatrix(blockInfo.record(recordID*2), data(recordID), 'int32');
CharData = nan(size(IntData)); % not used
switch mOP2block.functionCodes(data(recordID).TCODE, 1)
case 1
EKEY = IntData(:,1);
data(recordID).NodeID = floor(EKEY*.1);
case 2
switch mOP2block.functionCodes(data(recordID).ACODE, 4)
case {1 2 3 4 7 8 9 11 12}
EKEY = IntData(:,1);
data(recordID).NodeID = floor(EKEY*.1);
case {5 6 10}
data(recordID).FQ_TS = IntData(:,1);
otherwise
error('mOP2block_OQG: wrong ACODE')
end
end
% Define dataType
dataType = '';
switch mOP2block.functionCodes(data(recordID).TCODE, 7)
case {0 2}
% data(recordID).Value = FloatData(:,3:8);
dataType = ['nn' repmat('f', 1, 8-3+1)];
case 1
% data(recordID).Value = FloatData(:,3:14);
dataType = ['nn' repmat('f', 1, 14-3+1)];
otherwise
error('mOP2block_OQG: wrong TCODE')
end
% Create output data
if numel(dataType) > 0
data(recordID).Value = block.buildTable(dataType, IntData, FloatData, CharData);
end
end
% If the record is empty, initialize "data" as []
if length(blockInfo.record)-1 == 0
data = [];
end
block.data_set = true;
block.data_copy = data;
end % get.data
end % methods
end