tailieunhanh - Introduction to Optimum Design phần 10

Tham khảo tài liệu 'introduction to optimum design phần 10', công nghệ thông tin, kỹ thuật lập trình phục vụ nhu cầu học tập, nghiên cứu và làm việc hiệu quả | IF AU - AL .LE. EPSLON GO TO 80 C C IMPLEMENT STEPS 4 5 OR 6 OF THE ALGORITHM C IF FA - FB 40 60 50 C C FA IS LESS THAN FB STEP 4 40 AU AB FU FB AB AA FB FA AA AL AU - AL - GR CALL FUNCT AA FA NCOUNT GO TO 30 C C FA IS GREATER THAN FB STEP 5 C 50 AL AA FL FA AA AB FA FB AB AL AU - AL GR CALL FUNCT AB FB NCOUNT GO TO 30 C C FA IS EQUAL TO FB STEP 6 C 60 AL AA FL FA AU AB FU FB AA AL - GR AU - AL CALL FUNCT AA FA NCOUNT AB AL AU - AL GR CALL FUNCT AB FB NCOUNT GO TO 30 C C MINIMUM IS FOUND C 80 ALFA AU AL CALL FUNCT ALFA F NCOUNT RETURN END FIGURE D-2 Continued 662 Appendix D Sample Computer Programs C THE MAIN PROGRAM FOR STEEPEST DESCENT METHOD C ------------------------------------------------------------ C DELTA INITIAL STEP LENGTH FOR LINE SEARCH C EPSLON LINE SEARCH ACCURACY C EPSL STOPPING CRITERION FOR STEEPEST DESCENT METHOD C NCOUNT NO. OF FUNCTION EVALUATIONS C NDV NO. OF DESIGN VARIABLES C NOC NO. OF CYCLES OF THE METHOD C X DESIGN VARIABLE VECTOR C D DIRECTION VECTOR C G GRADIENT VECTOR C WK WORK ARRAY USED FOR TEMPORARY STORAGE C ------------------------------------------------------------ IMPLICIT DOUBLE PRECISION A-H O-Z DIMENSION X 4 D 4 G 4 WK 4 C C DEFINE INITIAL DATA C DELTA EPSLON EPSL NCOUNT 0 NDV 3 NOC 100 C C STARTING VALUES OF THE DESIGN VARIABLES C X 1 x X 3 CALL GRAD X G NDV WRITE 10 10 FORMAT NO. COST FUNCT STEP SIZE NORM OF GRAD DO 20 K 1 NOC CALL SCALE G D NDV CALL GOLDM X D WK ALFA DELTA EPSLON F NCOUNT NDV CALL SCALE D D ALFA NDV CALL PRINT K X ALFA G F NDV CALL ADD X D X NDV CALL GRAD X G NDV IF TNORM G NDV .LE. EPSL GO TO 30 20 CONTINUE FIGURE D-3 Computer program for steepest descent method. Appendix D Sample Computer Programs 663 WRITE WRITE LIMIT ON NO. OF CYCLES HAS EXCEEDED WRITe the current design variables ARE write X CALL EXIT 30 WRITE WRITe THE OPTIMAL DESIGN VARIABLES ARE WRITe 40 X 40 FORMAT CALL FUNCT X F NCOUNT NDV WRITE 50 the optimum .