主要内容

从USAF数字DATCOM文件导入

此示例演示如何使用航空航天工具箱将美国空军(USAF)数字DATCOM文件引入MATLAB®环境™ 软件。

美国空军数字数据通信文件示例

以下是美国空军数字数据通信的输入文件样本,用于机翼-机身-水平尾翼-垂直尾翼配置,运行超过5个阿尔法,2马赫数,2个高度,并计算静态和动态导数:

类型astdatcom.in
2.0,ALT(1)=5000.08000.0,BLREF=41.15美元。该参考参考参考=41.15美元。合成合成的X X=7.08,参考参考参考参考=41.15美元。合成的X=7.08,ZCTG=7.08,ZCG=0.0=0.0,ZCHCHD(1)为0.0.0,AlTCN=5,阿尔TCN=5,AlTCN=5,AlTCON=5,AlTCON=5,AlTCON=5,AlTCON=5,AlCHD(1,AlCHCHD(1),AlCHD(1),AlCHCHD(1=5,AlCHD(1=5),AlCHD(1=5),AlCHCHD(1=5),AlCHCHD(1=5),AlCHD(1=5),AlCHCHCHD(1=2.0=2.0,1=2.0,1=2.0,0,0 0,VERTUP=.TRUE。$$BODY NX=10.0,X(1)=-4.9,0.0,3.0,6.1,9.1,13.3,20.2,23.5,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,0.0,0.0.0,$WGPPLNF CHRDP=4.0,SSPNE=18.7,SSPNE=18.7,SSPNE=18.7,SSPNP=20.7,SSPN=20.6,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,0,DHDAD=0,DHDADID=3.0,DHDAD=3.0,DHDADO=3.0,D=3.0,D=3.0,类型=1.0,类型=1.0美元NACA-1.0美元NACA-1.0-NACA-6-6-W-6-6-6-6-6-6-6-6-6-6-6-64AATwistA=0.0,类型=1.0$NACA-H-4-0012$VTPLNF CHRDTP=2.7,SSPNE=5.0,SSPN=5.2,CHRDR=5.3,SAVSI=31.3,CHSTAT=0.25,TWISTA=0.0,类型=1.0$NACA-V-4-0012 CASEID SKYHOGG机身-机翼-水平尾翼-垂直尾翼-下一个案例

这是美国空军数字数据通信生成的相同的机翼-机身-水平尾翼-垂直尾翼配置的输出文件,运行超过5个阿尔法,2马赫数,2个高度:

类型astdatcom.out
本软件及随附文档发布的“AS IS”。美国政府不做任何明示或暗示的任何保证,以这软件和任何附带的文档,包括但不限于任何适销性或适合于特定用途。IN NO EVENT WILL THE U.S. GOVERNMENT BE LIABLE FOR ANY DAMAGES, INCLUDING LOST PROFITS, LOST SAVINGS OR OTHER INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING OUT OF THE USE, OR INABILITY TO USE, THIS SOFTWARE OR ANY ACCOMPANYING DOCUMENTATION, EVEN IF INFORMED IN ADVANCE OF THE POSSIBILITY OF SUCH DAMAGES. **************************************************** * USAF STABILITY AND CONTROL DIGITAL DATCOM * * PROGRAM REV. JAN 96 DIRECT INQUIRIES TO: * * WRIGHT LABORATORY (WL/FIGC) ATTN: W. BLAKE * * WRIGHT PATTERSON AFB, OHIO 45433 * * PHONE (513) 255-6764, FAX (513) 258-4054 * **************************************************** 1 CONERR - INPUT ERROR CHECKING 0 ERROR CODES - N* DENOTES THE NUMBER OF OCCURENCES OF EACH ERROR 0 A - UNKNOWN VARIABLE NAME 0 B - MISSING EQUAL SIGN FOLLOWING VARIABLE NAME 0 C - NON-ARRAY VARIABLE HAS AN ARRAY ELEMENT DESIGNATION - (N) 0 D - NON-ARRAY VARIABLE HAS MULTIPLE VALUES ASSIGNED 0 E - ASSIGNED VALUES EXCEED ARRAY DIMENSION 0 F - SYNTAX ERROR 0****************************** INPUT DATA CARDS ****************************** $FLTCON NMACH=2.0,MACH(1)=0.1,0.2$ $FLTCON NALT=2.0,ALT(1)=5000.0,8000.0$ $FLTCON NALPHA=5.,ALSCHD(1)=-2.0,0.0,2.0, ALSCHD(4)=4.0,8.0,LOOP=2.0$ $OPTINS SREF=225.8,CBARR=5.75,BLREF=41.15$ $SYNTHS XCG=7.08,ZCG=0.0,XW=6.1,ZW=-1.4,ALIW=1.1,XH=20.2, ZH=0.4,ALIH=0.0,XV=21.3,ZV=0.0,VERTUP=.TRUE.$ $BODY NX=10.0, X(1)=-4.9,0.0,3.0,6.1,9.1,13.3,20.2,23.5,25.9, R(1)=0.0,1.0,1.75,2.6,2.6,2.6,2.0,1.0,0.0$ $WGPLNF CHRDTP=4.0,SSPNE=18.7,SSPN=20.6,CHRDR=7.2,SAVSI=0.0,CHSTAT=0.25, TWISTA=-1.1,SSPNDD=0.0,DHDADI=3.0,DHDADO=3.0,TYPE=1.0$ NACA-W-6-64A412 $HTPLNF CHRDTP=2.3,SSPNE=5.7,SSPN=6.625,CHRDR=0.25,SAVSI=11.0, CHSTAT=1.0,TWISTA=0.0,TYPE=1.0$ NACA-H-4-0012 $VTPLNF CHRDTP=2.7,SSPNE=5.0,SSPN=5.2,CHRDR=5.3,SAVSI=31.3, CHSTAT=0.25,TWISTA=0.0,TYPE=1.0$ NACA-V-4-0012 CASEID SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG DAMP NEXT CASE 1 THE FOLLOWING IS A LIST OF ALL INPUT CARDS FOR THIS CASE. 0 $FLTCON NMACH=2.0,MACH(1)=0.1,0.2$ $FLTCON NALT=2.0,ALT(1)=5000.0,8000.0$ $FLTCON NALPHA=5.,ALSCHD(1)=-2.0,0.0,2.0, ALSCHD(4)=4.0,8.0,LOOP=2.0$ $OPTINS SREF=225.8,CBARR=5.75,BLREF=41.15$ $SYNTHS XCG=7.08,ZCG=0.0,XW=6.1,ZW=-1.4,ALIW=1.1,XH=20.2, ZH=0.4,ALIH=0.0,XV=21.3,ZV=0.0,VERTUP=.TRUE.$ $BODY NX=10.0, X(1)=-4.9,0.0,3.0,6.1,9.1,13.3,20.2,23.5,25.9, R(1)=0.0,1.0,1.75,2.6,2.6,2.6,2.0,1.0,0.0$ $WGPLNF CHRDTP=4.0,SSPNE=18.7,SSPN=20.6,CHRDR=7.2,SAVSI=0.0,CHSTAT=0.25, TWISTA=-1.1,SSPNDD=0.0,DHDADI=3.0,DHDADO=3.0,TYPE=1.0$ NACA-W-6-64A412 $HTPLNF CHRDTP=2.3,SSPNE=5.7,SSPN=6.625,CHRDR=0.25,SAVSI=11.0, CHSTAT=1.0,TWISTA=0.0,TYPE=1.0$ NACA-H-4-0012 $VTPLNF CHRDTP=2.7,SSPNE=5.0,SSPN=5.2,CHRDR=5.3,SAVSI=31.3, CHSTAT=0.25,TWISTA=0.0,TYPE=1.0$ NACA-V-4-0012 CASEID SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG DAMP NEXT CASE 0 INPUT DIMENSIONS ARE IN FT, SCALE FACTOR IS 1.0000 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM WING SECTION DEFINITION 0 IDEAL ANGLE OF ATTACK = 0.00000 DEG. ZERO LIFT ANGLE OF ATTACK = -3.09292 DEG. IDEAL LIFT COEFFICIENT = 0.40000 ZERO LIFT PITCHING MOMENT COEFFICIENT = -0.08719 MACH ZERO LIFT-CURVE-SLOPE = 0.09654 /DEG. LEADING EDGE RADIUS = 0.00993 FRACTION CHORD MAXIMUM AIRFOIL THICKNESS = 0.12000 FRACTION CHORD DELTA-Y = 2.46808 PERCENT CHORD 0 MACH= 0.1000 LIFT-CURVE-SLOPE = 0.09693 /DEG. XAC = 0.26404 0 MACH= 0.2000 LIFT-CURVE-SLOPE = 0.09811 /DEG. XAC = 0.26457 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM HORIZONTAL TAIL SECTION DEFINITION 0 IDEAL ANGLE OF ATTACK = 0.00000 DEG. ZERO LIFT ANGLE OF ATTACK = 0.00000 DEG. IDEAL LIFT COEFFICIENT = 0.00000 ZERO LIFT PITCHING MOMENT COEFFICIENT = 0.00000 MACH ZERO LIFT-CURVE-SLOPE = 0.09596 /DEG. LEADING EDGE RADIUS = 0.01587 FRACTION CHORD MAXIMUM AIRFOIL THICKNESS = 0.12000 FRACTION CHORD DELTA-Y = 3.16898 PERCENT CHORD 0 MACH= 0.1000 LIFT-CURVE-SLOPE = 0.09636 /DEG. XAC = 0.25854 0 MACH= 0.2000 LIFT-CURVE-SLOPE = 0.09761 /DEG. XAC = 0.25881 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM VERTICAL TAIL SECTION DEFINITION 0 IDEAL ANGLE OF ATTACK = 0.00000 DEG. ZERO LIFT ANGLE OF ATTACK = 0.00000 DEG. IDEAL LIFT COEFFICIENT = 0.00000 ZERO LIFT PITCHING MOMENT COEFFICIENT = 0.00000 MACH ZERO LIFT-CURVE-SLOPE = 0.09596 /DEG. LEADING EDGE RADIUS = 0.01587 FRACTION CHORD MAXIMUM AIRFOIL THICKNESS = 0.12000 FRACTION CHORD DELTA-Y = 3.16898 PERCENT CHORD 0 MACH= 0.1000 LIFT-CURVE-SLOPE = 0.09636 /DEG. XAC = 0.25854 0 MACH= 0.2000 LIFT-CURVE-SLOPE = 0.09761 /DEG. XAC = 0.25881 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM CHARACTERISTICS AT ANGLE OF ATTACK AND IN SIDESLIP WING-BODY-VERTICAL TAIL-HORIZONTAL TAIL CONFIGURATION SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG ----------------------- FLIGHT CONDITIONS ------------------------ -------------- REFERENCE DIMENSIONS ------------ MACH ALTITUDE VELOCITY PRESSURE TEMPERATURE REYNOLDS REF. REFERENCE LENGTH MOMENT REF. CENTER NUMBER NUMBER AREA LONG. LAT. HORIZ VERT FT FT/SEC LB/FT**2 DEG R 1/FT FT**2 FT FT FT FT 0 0.100 5000.00 109.70 1.7609E+03 500.843 6.1507E+05 225.800 5.750 41.150 7.080 0.000 0 -------------------DERIVATIVE (PER DEGREE)------------------- 0 ALPHA CD CL CM CN CA XCP CLA CMA CYB CNB CLB 0 -2.0 0.032 0.113 -0.0340 0.112 0.035 -0.304 8.926E-02 -2.105E-02 -3.458E-03 9.142E-04 -6.161E-04 0.0 0.035 0.296 -0.0752 0.296 0.035 -0.254 9.350E-02 -2.034E-02 -6.205E-04 2.0 0.042 0.487 -0.1153 0.488 0.025 -0.236 9.732E-02 -1.971E-02 -6.268E-04 4.0 0.052 0.685 -0.1541 0.687 0.004 -0.224 1.005E-01 -1.927E-02 -6.349E-04 8.0 0.084 1.098 -0.2304 1.099 -0.069 -0.210 1.059E-01 -1.890E-02 -6.554E-04 0 ALPHA Q/QINF EPSLON D(EPSLON)/D(ALPHA) 0 -2.0 1.000 0.953 0.571 0.0 1.000 2.094 0.583 2.0 1.000 3.284 0.606 4.0 1.000 4.520 0.610 8.0 1.000 6.897 0.594 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM DYNAMIC DERIVATIVES WING-BODY-VERTICAL TAIL-HORIZONTAL TAIL CONFIGURATION SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG ----------------------- FLIGHT CONDITIONS ------------------------ -------------- REFERENCE DIMENSIONS ------------ MACH ALTITUDE VELOCITY PRESSURE TEMPERATURE REYNOLDS REF. REFERENCE LENGTH MOMENT REF. CENTER NUMBER NUMBER AREA LONG. LAT. HORIZ VERT FT FT/SEC LB/FT**2 DEG R 1/FT FT**2 FT FT FT FT 0 0.100 5000.00 109.70 1.7609E+03 500.843 6.1507E+05 225.800 5.750 41.150 7.080 0.000 DYNAMIC DERIVATIVES (PER DEGREE) 0 -------PITCHING------- -----ACCELERATION------ --------------ROLLING-------------- --------YAWING-------- 0 ALPHA CLQ CMQ CLAD CMAD CLP CYP CNP CNR CLR 0 -2.00 9.739E-02 -8.918E-02 1.874E-02 -4.247E-02 -7.824E-03 -1.516E-03 -1.498E-04 -1.059E-03 6.334E-04 0.00 1.913E-02 -4.336E-02 -8.226E-03 -1.649E-03 -4.034E-04 -1.068E-03 1.240E-03 2.00 1.991E-02 -4.512E-02 -8.599E-03 -1.792E-03 -6.631E-04 -1.073E-03 1.878E-03 4.00 2.003E-02 -4.540E-02 -8.890E-03 -1.942E-03 -9.290E-04 -1.073E-03 2.542E-03 8.00 1.952E-02 -4.424E-02 -9.387E-03 -2.262E-03 -1.479E-03 -1.060E-03 3.926E-03 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM CHARACTERISTICS AT ANGLE OF ATTACK AND IN SIDESLIP WING-BODY-VERTICAL TAIL-HORIZONTAL TAIL CONFIGURATION SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG ----------------------- FLIGHT CONDITIONS ------------------------ -------------- REFERENCE DIMENSIONS ------------ MACH ALTITUDE VELOCITY PRESSURE TEMPERATURE REYNOLDS REF. REFERENCE LENGTH MOMENT REF. CENTER NUMBER NUMBER AREA LONG. LAT. HORIZ VERT FT FT/SEC LB/FT**2 DEG R 1/FT FT**2 FT FT FT FT 0 0.200 5000.00 219.39 1.7609E+03 500.843 1.2301E+06 225.800 5.750 41.150 7.080 0.000 0 -------------------DERIVATIVE (PER DEGREE)------------------- 0 ALPHA CD CL CM CN CA XCP CLA CMA CYB CNB CLB 0 -2.0 0.028 0.114 -0.0335 0.113 0.032 -0.297 9.000E-02 -2.124E-02 -3.465E-03 8.781E-04 -6.226E-04 0.0 0.031 0.298 -0.0751 0.298 0.031 -0.252 9.421E-02 -2.051E-02 -6.270E-04 2.0 0.038 0.491 -0.1155 0.492 0.021 -0.235 9.800E-02 -1.987E-02 -6.332E-04 4.0 0.048 0.690 -0.1546 0.692 0.000 -0.223 1.011E-01 -1.943E-02 -6.413E-04 8.0 0.081 1.105 -0.2316 1.106 -0.074 -0.209 1.065E-01 -1.906E-02 -6.614E-04 0 ALPHA Q/QINF EPSLON D(EPSLON)/D(ALPHA) 0 -2.0 1.000 0.957 0.573 0.0 1.000 2.103 0.585 2.0 1.000 3.297 0.609 4.0 1.000 4.537 0.612 8.0 1.000 6.923 0.596 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM DYNAMIC DERIVATIVES WING-BODY-VERTICAL TAIL-HORIZONTAL TAIL CONFIGURATION SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG ----------------------- FLIGHT CONDITIONS ------------------------ -------------- REFERENCE DIMENSIONS ------------ MACH ALTITUDE VELOCITY PRESSURE TEMPERATURE REYNOLDS REF. REFERENCE LENGTH MOMENT REF. CENTER NUMBER NUMBER AREA LONG. LAT. HORIZ VERT FT FT/SEC LB/FT**2 DEG R 1/FT FT**2 FT FT FT FT 0 0.200 5000.00 219.39 1.7609E+03 500.843 1.2301E+06 225.800 5.750 41.150 7.080 0.000 DYNAMIC DERIVATIVES (PER DEGREE) 0 -------PITCHING------- -----ACCELERATION------ --------------ROLLING-------------- --------YAWING-------- 0 ALPHA CLQ CMQ CLAD CMAD CLP CYP CNP CNR CLR 0 -2.00 9.840E-02 -8.993E-02 1.900E-02 -4.307E-02 -7.877E-03 -1.525E-03 -1.499E-04 -1.057E-03 6.448E-04 0.00 1.940E-02 -4.398E-02 -8.276E-03 -1.659E-03 -4.038E-04 -1.066E-03 1.264E-03 2.00 2.018E-02 -4.574E-02 -8.646E-03 -1.802E-03 -6.637E-04 -1.070E-03 1.915E-03 4.00 2.030E-02 -4.602E-02 -8.934E-03 -1.953E-03 -9.297E-04 -1.070E-03 2.593E-03 8.00 1.978E-02 -4.483E-02 -9.423E-03 -2.273E-03 -1.479E-03 -1.057E-03 4.003E-03 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM CHARACTERISTICS AT ANGLE OF ATTACK AND IN SIDESLIP WING-BODY-VERTICAL TAIL-HORIZONTAL TAIL CONFIGURATION SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG ----------------------- FLIGHT CONDITIONS ------------------------ -------------- REFERENCE DIMENSIONS ------------ MACH ALTITUDE VELOCITY PRESSURE TEMPERATURE REYNOLDS REF. REFERENCE LENGTH MOMENT REF. CENTER NUMBER NUMBER AREA LONG. LAT. HORIZ VERT FT FT/SEC LB/FT**2 DEG R 1/FT FT**2 FT FT FT FT 0 0.100 8000.00 108.52 1.5721E+03 490.151 5.6457E+05 225.800 5.750 41.150 7.080 0.000 0 -------------------DERIVATIVE (PER DEGREE)------------------- 0 ALPHA CD CL CM CN CA XCP CLA CMA CYB CNB CLB 0 -2.0 0.032 0.113 -0.0340 0.112 0.036 -0.305 8.926E-02 -2.106E-02 -3.458E-03 9.190E-04 -6.161E-04 0.0 0.035 0.296 -0.0753 0.296 0.035 -0.254 9.350E-02 -2.034E-02 -6.205E-04 2.0 0.042 0.487 -0.1154 0.488 0.025 -0.236 9.732E-02 -1.971E-02 -6.268E-04 4.0 0.052 0.685 -0.1541 0.687 0.004 -0.224 1.005E-01 -1.927E-02 -6.349E-04 8.0 0.085 1.098 -0.2304 1.099 -0.069 -0.210 1.059E-01 -1.891E-02 -6.554E-04 0 ALPHA Q/QINF EPSLON D(EPSLON)/D(ALPHA) 0 -2.0 1.000 0.953 0.571 0.0 1.000 2.094 0.583 2.0 1.000 3.284 0.606 4.0 1.000 4.520 0.610 8.0 1.000 6.897 0.594 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM DYNAMIC DERIVATIVES WING-BODY-VERTICAL TAIL-HORIZONTAL TAIL CONFIGURATION SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG ----------------------- FLIGHT CONDITIONS ------------------------ -------------- REFERENCE DIMENSIONS ------------ MACH ALTITUDE VELOCITY PRESSURE TEMPERATURE REYNOLDS REF. REFERENCE LENGTH MOMENT REF. CENTER NUMBER NUMBER AREA LONG. LAT. HORIZ VERT FT FT/SEC LB/FT**2 DEG R 1/FT FT**2 FT FT FT FT 0 0.100 8000.00 108.52 1.5721E+03 490.151 5.6457E+05 225.800 5.750 41.150 7.080 0.000 DYNAMIC DERIVATIVES (PER DEGREE) 0 -------PITCHING------- -----ACCELERATION------ --------------ROLLING-------------- --------YAWING-------- 0 ALPHA CLQ CMQ CLAD CMAD CLP CYP CNP CNR CLR 0 -2.00 9.739E-02 -8.918E-02 1.874E-02 -4.247E-02 -7.824E-03 -1.516E-03 -1.498E-04 -1.060E-03 6.334E-04 0.00 1.913E-02 -4.336E-02 -8.226E-03 -1.649E-03 -4.034E-04 -1.069E-03 1.240E-03 2.00 1.991E-02 -4.512E-02 -8.599E-03 -1.792E-03 -6.631E-04 -1.073E-03 1.878E-03 4.00 2.003E-02 -4.540E-02 -8.890E-03 -1.942E-03 -9.290E-04 -1.074E-03 2.542E-03 8.00 1.952E-02 -4.424E-02 -9.387E-03 -2.262E-03 -1.479E-03 -1.061E-03 3.926E-03 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM CHARACTERISTICS AT ANGLE OF ATTACK AND IN SIDESLIP WING-BODY-VERTICAL TAIL-HORIZONTAL TAIL CONFIGURATION SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG ----------------------- FLIGHT CONDITIONS ------------------------ -------------- REFERENCE DIMENSIONS ------------ MACH ALTITUDE VELOCITY PRESSURE TEMPERATURE REYNOLDS REF. REFERENCE LENGTH MOMENT REF. CENTER NUMBER NUMBER AREA LONG. LAT. HORIZ VERT FT FT/SEC LB/FT**2 DEG R 1/FT FT**2 FT FT FT FT 0 0.200 8000.00 217.04 1.5721E+03 490.151 1.1291E+06 225.800 5.750 41.150 7.080 0.000 0 -------------------DERIVATIVE (PER DEGREE)------------------- 0 ALPHA CD CL CM CN CA XCP CLA CMA CYB CNB CLB 0 -2.0 0.028 0.114 -0.0335 0.113 0.032 -0.297 9.000E-02 -2.124E-02 -3.465E-03 8.829E-04 -6.226E-04 0.0 0.031 0.298 -0.0751 0.298 0.031 -0.252 9.421E-02 -2.051E-02 -6.270E-04 2.0 0.038 0.491 -0.1156 0.492 0.021 -0.235 9.800E-02 -1.987E-02 -6.332E-04 4.0 0.049 0.690 -0.1546 0.692 0.000 -0.223 1.011E-01 -1.943E-02 -6.413E-04 8.0 0.081 1.105 -0.2316 1.106 -0.073 -0.209 1.065E-01 -1.906E-02 -6.614E-04 0 ALPHA Q/QINF EPSLON D(EPSLON)/D(ALPHA) 0 -2.0 1.000 0.957 0.573 0.0 1.000 2.103 0.585 2.0 1.000 3.297 0.609 4.0 1.000 4.537 0.612 8.0 1.000 6.923 0.596 1 AUTOMATED STABILITY AND CONTROL METHODS PER APRIL 1976 VERSION OF DATCOM DYNAMIC DERIVATIVES WING-BODY-VERTICAL TAIL-HORIZONTAL TAIL CONFIGURATION SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL TAIL CONFIG ----------------------- FLIGHT CONDITIONS ------------------------ -------------- REFERENCE DIMENSIONS ------------ MACH ALTITUDE VELOCITY PRESSURE TEMPERATURE REYNOLDS REF. REFERENCE LENGTH MOMENT REF. CENTER NUMBER NUMBER AREA LONG. LAT. HORIZ VERT FT FT/SEC LB/FT**2 DEG R 1/FT FT**2 FT FT FT FT 0 0.200 8000.00 217.04 1.5721E+03 490.151 1.1291E+06 225.800 5.750 41.150 7.080 0.000 DYNAMIC DERIVATIVES (PER DEGREE) 0 -------PITCHING------- -----ACCELERATION------ --------------ROLLING-------------- --------YAWING-------- 0 ALPHA CLQ CMQ CLAD CMAD CLP CYP CNP CNR CLR 0 -2.00 9.840E-02 -8.993E-02 1.900E-02 -4.307E-02 -7.877E-03 -1.525E-03 -1.499E-04 -1.057E-03 6.448E-04 0.00 1.940E-02 -4.398E-02 -8.276E-03 -1.659E-03 -4.038E-04 -1.066E-03 1.264E-03 2.00 2.018E-02 -4.574E-02 -8.646E-03 -1.802E-03 -6.637E-04 -1.071E-03 1.915E-03 4.00 2.030E-02 -4.602E-02 -8.934E-03 -1.953E-03 -9.297E-04 -1.071E-03 2.593E-03 8.00 1.978E-02 -4.483E-02 -9.424E-03 -2.273E-03 -1.479E-03 -1.057E-03 4.003E-03 1 THE FOLLOWING IS A LIST OF ALL INPUT CARDS FOR THIS CASE. 0 1 END OF JOB.

从DATCOM文件导入数据

使用datcomimport函数将数字DATCOM数据带入MATLAB。

alldata=datcomimport(“astdatcom.out”,对,0);

检查导入的DATCOM数据

datcomimport函数创建包含数字DATCOM输出文件中数据的结构单元数组。

data=alldata{1}
=数据结构体字段:案例:“SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL尾配置“马赫:[0.1000 - 0.2000]alt:α(5000 8000):[2 0 2 4 8]nmach: 2 nalt: 2 nalpha: 5 rnnub:[]大卖场:0循环:2 sref: 225.8000 cbar: 5.7500 blref: 41.1500暗淡:“英尺”引出:“度”stmach: 0.6000 tsmach: 1.4000保存:0药栓:[]修剪:0潮湿:1建立:1部分:0 highsym: 0 highasy: 0 highcon: 0 tjet: 0 hypeff: 0磅:0压水式反应堆:0接地:0 wsspn: 18.7000 hsspn: 5.7000 ndelta: 0三角洲:[]deltal: [] deltar:[]已:0 grndht:[]配置:[1 x1 struct]版本:1976 cd: [5 x2x2双重的]cl: [5 x2x2双重的]厘米:[5 x2x2双]的cn: [5 x2x2双]的ca: [5 x2x2双]的xcp: [5 x2x2双重的]cma: [5 x2x2双)的地方:[5x2x2 double] clp: [5x2x2 double] clp: [5x2x2 double] clp: [5x2x2 double] clr: [5x2x2 double] clp: [5x2x2 double] clr: [5x2x2 double] clp: [5x2x2 double] clr: [5x2x2 double] clp: [5x2x2 double] clr: [5x2x2 double] clp: [5x2x2 double] clr: [5x2x2 double] clp: [5x2x2 double] clr: [5x2x2 double] clp: [5x2x2 double] clr: [5x2x2 double] clr: [5x2x2 double] clr: [5x2x2 double]

填写丢失的DATCOM数据

默认情况下,如果不存在DATCOM方法或该方法不适用,缺失数据点设置为99999,数据点设置为NaN。

可以看到,在数字DATCOM输出文件和检查导入的数据

$$C{Y\beta},$$
;$$C{n\beta},$$和#xA;$$C{Lq},$$和
;$$C_{mq}$$

仅在第一个alpha值中包含数据。以下是导入的数据值。

data.cyb
ans(:,:,1)=1.0e+04*-0.0000-0.0000 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999 ans(:,:,2)=1.0e+04*-0.0000-0.0000 9.99999 9.99999 9.99999 9.99999 9.99999
data.cnb
ans(:,:,1)=1.0e+04*0.0000 0.0000 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999 ans(:,:,2)=1.0e+04*0.0000 0.0000 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999
data.clq
ans(:,:,1)=1.0e+04*0.0000 0.0000 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999 ans(:,:,2)=1.0e+04*0.0000 0.0000 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999
data.cmq
ans(:,:,1)=1.0e+04*-0.0000-0.0000 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999 9.99999 ans(:,:,2)=1.0e+04*-0.0000-0.0000 9.99999 9.99999 9.99999 9.99999 9.99999

丢失的数据点将用第一个alpha的值来填充,因为这些数据点意味着用于所有alpha值。

aerotab = {“地方”“cnb”“clq”“cmq的,};对于k = 1:长度(aerotab)对于m=1:data.nmach对于h = 1:数据。nalt数据。(aerotab {k}) (:, m h) =数据。(aerotab {k}) (1 m h);结束结束结束

以下是更新后的导入数据值:

data.cyb
ans(:,:,1)=-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035-0.0035
data.cnb
ans(:,:,1)=1.0e-03*0.9142 0.8781 0.9142 0.8781 0.9142 0.8781 0.9142 0.8781 0.9142 0.8781 ans(:,:,2)=1.0e-03*0.9190 0.8829 0.9190 0.8829 0.9190 0.8829 0.9190 0.8829 0.9190 0.8829
data.clq
ans(:,:,1)=0.0974 0.0984 0.0974 0.0984 0.0984 0.0974 0.0984 0.0974 0.0984 ans(:,:,2)=0.0974 0.0984 0.0974 0.0984 0.0974 0.0984 0.0984 0.0984 0.0974 0.0984
data.cmq
ans (:,: 1) = -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 ans (:,: 2) = -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899

绘制空气动力系数

绘制升力曲线、阻力极轴和俯仰力矩。

h1 =图;figtitle = {的升程曲线''};对于k=1:2子地块(2,1,k)图(data.alpha,permute(data.cl(:,k,:),[1,3,2]))网格标签(['升力系数(马赫数)'num2str(数据马赫数(k))')'])标题(figtitle{k});结束xlabel(‘攻角(度)’

h2=图;figtitle={“拖极”''};对于k=1:2子地块(2,1,k)图(排列(data.cd(:,k,:),[1 3 2]),排列(data.cl(:,k,:),[1 3 2])网格标签(['升力系数(马赫数)'num2str(数据马赫数(k))')'])标题(figtitle{k})结束xlabel(“阻力系数”

h3 =图;figtitle = {“俯仰力矩”''};对于k = 1:2次要情节(2,1,k)情节(排列(data.cm (k,::),[1 3 2]),排列(data.cl (k,::),[1 3 2]))网格ylabel (['升力系数(马赫数)'num2str(数据马赫数(k))')'])标题(figtitle{k})结束xlabel(“俯仰力矩系数”

关闭(h1, h2, h3);%#好的<*NOPTS>