설명

저널 무료 슬롯의 편심 e가 주어지면 압력을 구하십시오.
무료 슬롯에서 샤프트를 분리하는 윤활유의 분포.

이 모델은 COPS 벤치마킹 제품군에서 나온 것입니다.
참조http://www-unix.mcs.anl.gov/~more/cops/.내부 그리드 포인트 수는 다음 명령을 사용하여 지정할 수 있습니다
라인 매개변수 --nx 및 --ny. COPS 성능 테스트가 완료되었습니다.
nx-1 = 50, ny-1 = 25, 50, 75, 100에 대해 보고됨

대형 모델 유형 :NLP

카테고리 : 무료 슬롯 모델 라이브러리

메인 파일 : jbearing.gms

$title 저널 무료 슬롯 COPS 2.0 #16 (JBEARING,SEQ=244)

$onText
저널 무료 슬롯의 편심률 e가 주어지면 압력을 구하십시오.
무료 슬롯에서 샤프트를 분리하는 윤활유의 분포.

이 모델은 COPS 벤치마킹 제품군에서 나온 것입니다.
http://www-unix.mcs.anl.gov/~more/cops/를 참조하세요.

내부 격자점 수는 다음 명령을 사용하여 지정할 수 있습니다.
라인 매개변수 --nx 및 --ny. COPS 성능 테스트가 완료되었습니다.
nx-1 = 50, ny-1 = 25, 50, 75, 100에 대해 보고됨

Dolan, E D 등, JJ, 벤치마킹 최적화
COPS가 포함된 소프트웨어. 기술. 대표, 수학과 컴퓨터
과학부, 2000.

Capriz, G 및 Cimatti, G, 자유 경계 문제
유체역학적 윤활 이론: 조사. 파사노에서는 A,
및 Primicerio, M, Eds, 수학 연구 노트. 피트맨,
1983, pp. 613-635.

Averick, B M, Carter, R G, More, J J 및 Xue, G L, The
MINPACK-2 테스트 문제집. 기술. 대표, 수학과
아르곤 국립 연구소(Argonne National Laboratory) 컴퓨터 과학부,
1992.

키워드: 비선형 프로그래밍, 엔지니어링, 저널 무료 슬롯
$offText

$설정되지 않은 경우 nx $set nx 51
$설정되지 않은 경우 $set ny 51

세트
   nx / 0*%nx% /
   뉴욕 / 0*%ny% /;

별칭 (nx,i), (ny,j);

스칼라
   b '그리드는 (0,2*pi)x(0,2*b)' / 10 /
   e '이심률' / 0.1 /
   hx '그리드 간격'
   hy '그리드 간격'
   면적 '삼각형의 면적';

hx = 2*pi/%nx%;
hy = 2*b/%ny%;
면적 = 0.5*hx*hy;

매개변수 wq(nx);
wq(nx) = power(1 + e*cos((ord(nx) - 1)*hx),3);

양수 변수 v(nx,ny);

변수 obj;

방정식 defobj '목적 함수';

디포브..
   obj =e= (hx*hy/12)*sum(nx(i+1),ny(j+1)), (wq[i] + 2*wq[i+1])
                         *(sqr((v[i+1,j] - v[i,j])/hx) + sqr((v[i,j+1] - v[i,j])/hy))
        + (hx*hy/12)*sum(nx(i+1),ny(j+1)), (2*wq[i+1] + 2*wq[i])
                         *(sqr((v[i,j+1] - v[i+1,j+1])/hx) + sqr((v[i+1,j] - v[i+1,j+1])/hy))
        - hx*hy*sum(nx,ny), e*sin((ord(nx) - 1)*hx)*v[nx,ny];

* 출발점
v.l[nx,ny] = max(sin((ord(nx) - 1)*hx),0);
v.fx[nx, '0'] = 0;
v.fx[nx,'%ny%'] = 0;
v.fx[ '0',ny] = 0;
v.fx['%nx%',ny] = 0;

모델 무료 슬롯 / 모두 /;

$if set workSpace 무료 슬롯.workSpace = %workSpace%

Bearing.workFactor = 2;

nlp를 사용하여 obj를 최소화하는 무료 슬롯을 해결합니다.