-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathtornadobatch_ge.asv
147 lines (91 loc) · 3.18 KB
/
tornadobatch_ge.asv
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
function [CLa]=tornadobatch(ac_name,state_name);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% tornadobatch, subsidary function to TORNADO
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% ac_name = name of geometry file (string)
% state_name=name of state file (string)
% JID=job identifier, output filename (string)
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
results=[];
i=0;
j=0;
settings=config('startup')
cd(settings.acdir)
load('Z1');
cd(settings.hdir)
cd(settings.sdir)
load('E4');
cd(settings.hdir)
quest=1; %Simple state solution, could be something else
%but then you'll have to change below too,
%especially the "Load data" section.
JID='batchjob';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% CHANGE GEOMETRY
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Here you may enter whatever changes to the geometry you like
%
% Examples:
% Change number of panels for easy grid convergence.
% Change a rudder setting, or tailplane twist for trim computations
% Change the geometry in any other way.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
geo=geo;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% CHANGE STATE
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Here you may enter whatever changes to the state you like
%
% Examples:
% Change angle of attack or sideslip.
% Change altitude, airspeed or ainything else
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
state.PGcorr=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
latticetype=1;%Standard VLM
%latticetype=0;%Tornado freestream following wake VLM
alpha=2*pi/180;
[Alattice,ref]=fLattice_setup2(geo,state,latticetype);
Alattice=wingrotation2(1,geo,Alattice,[0 1 0],[0 0 0],alpha)
[s,void]=size(Alattice.N);
cd(settings.acdir)
load('ground1');
cd(settings.hdir)
geo.startz=-5;
[Glattice,void]=fLattice_setup2(geo,state,latticetype);
cd(settings.acdir)
load('Z1');
cd(settings.hdir)
lattice.XYZ=[Alattice.XYZ;Glattice.XYZ];
lattice.N=[Alattice.N;Glattice.N];
lattice.COLLOC=[Alattice.COLLOC;Glattice.COLLOC];
lattice.VORTEX=[Alattice.VORTEX;Glattice.VORTEX];
%lattice.XYZ=[Alattice.XYZ Glattice.XYZ];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% CHANGE Reference values.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Here you may enter whatever changes to the state Reference
% you values like.
%
% Examples:
% Set the reference spane to a constant.
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ref=ref;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[results]=solver9(results,state,geo,lattice,ref);
Aresults.F=results.F(1:s,:,:);
Aresults.M=results.M(1:s,:,:);
Aresults.gamma=results.gamma(1:s,:)
Aresults.FORCE=sum(Aresults.F,1); %Total force
Aresults.MOMENTS=sum(Aresults.M,1); %Summing up moments
[results]=coeff_create3(Aresults,Alattice,state,ref,geo);
%cd(settings.odir)
% load batchjob-Cx;
%cd(settings.hdir)