혜 르. 101; &#E양 씨 계수,-나.축 단면 의 직경 방향 관성 모멘트,-k 초기 화mh -균열 을 완전히 벌 리 는 등가 링크 에 대응 하 는 모멘트 강도 계수. 이 방정식 은 균열 과 관련 된 유연성 과 무 균열 보 의 유연성 두 가지 가 있다.그것

 주의해 야 할 것 은 들보 의 모양,경계 조건 과 균열 의 위치 가 어떻든 간 에 일반 양 유도 에 대한 공헌 은 계수 에 만 달 려 있다. k 초기 화mh ,i.e.첫 번 째 총화. 들보 가 회전 할 때 균열 단면 의 유연성 에 변화 가 생 긴 다.우 리 는 다음 과 같은 정 보 를 얻 었 다.혜 르. 101; 

 g o&#들보 의 유연성 균열 없 음(  q) = -kmh부하 에 주어진 모멘트 강도 전류 계수q. -적분 두 방정식(5)과(7)을 얻 을 수 있 습 니 다.혜 르.101; 방정식(6)을 고려 하여 우 리 는 모멘트 의 강도 가 위상 차 에 따라 변화 하 는 규칙 을 얻 었 다.

 아인타이늄. q 

 :등가 링크 의 힘 모멘트 강 도 는 축 직경,재료 특성 과 균열 깊이 에 만 달 려 있다.축 속 균열 위치의 변화,축의 특성(지지 장치 포함)은 등가 링크 의 힘 모멘트 강 도 를 바 꾸 지 않 습 니 다.&#균열 이 있 는 단면 은 균열 이 있 는 단면 의 모든 부분 에서 같은 조건 을 유지한다.획득이전에 얻 은 강도 행렬 로 돌아 가면 강도 계 수 는 다음 과 같이 쓸 수 있 습 니 다. 혜 르. 

101; k

초기 화 

 ee ,  k  초기 화 

 응. 완전 균열 대응 축 상의 모멘트 강도 초기 값&#. 과업 k초기 화 ee, k초기 화  -응.계산 은 두 가지 방법 으로 해결 할 수 있다.하 나 는 유한 원 프로그램의 커 브 강도 계산 이다.둘 째 는 단열 역학 이론 을 운용 하 는 것 이다.개 구 부 균열 의 기하학 적 모양,축 직경 과 재료 특성 을 알 고 있다 면 개 구 부 균열 의 국부 유도 계 수 치 를 계산 할 기회 가 있다[7,10].

 균열 로터 시스템 시 뮬 레이 션 알고리즘 다음 절 차 를 통 해 균열 의 모멘트 강도 계 수 를 얻어 야 한다.1.전문 적 인 프로그램 에서 완전한 로터 모델 을 만들어 로터 동력학(예 를 들 어:예컨대 동력학 R4 에서).2.균열 이 있 는 로터 부분 을 강조 한다.  3.균열 은 축 을 두 개의 서브 시스템 으로 나눈다.변수 매트릭스 에서 설명 하 는 링크강도 계수[K （q） 

 , 

j） ]

빌리다

사이즈 6x6 은 서브 시스템 사이 에 있 습 니 다.

4.모멘트 강도 의 초기 계수

k

초기 화

ee , k초기 화응.제 시 된 방법 을 통 해 열 린 균열 의 해석 을 얻 었 다.위.이 데이터 들 은 초기 계산 데이터 이다.적분 운동 방정식 과 동시에 아 날로 그 균열 링크 의 강도 행렬 계 수 를 계산 했다. 

로터 시스템 의 모든q. 비 선형 동력학 모델 의 비 선형 상태 매트릭스 방정식 을 묘사 하 다.로터 시스템 의 구 조 는 다음 과 같다. 혜 르.101; [M ]

-관성 계수 행렬;

 [ C ]

-댐퍼 와 자 이 로 계수 행렬;[K

]  

 행렬 강도 계수; { &#u},{ u},{u} -해당 하 는 진동 가속도,진동 속도 와 진동 변위 열;{-F(t)} &&아무 동적 부하 유형-내부 와 외부.등가 링크 의 강도 행렬 은 두 부분 으로 나 눌 수 있다.상수 와 변 수 는 다음 과 같다.[&K C]구조의 일반 강도 행렬 에 포함[K] 계통행렬（）[q j]n, K계산 에 사용 하 다비 선형 링크 의 반응:-혜 르.101;  

 u 접수 하 다. ,  

 u 천둥.1 묘단면 이 상응하는 축의 실제 회전 을 감다.시스템 의 최종 운동 방정식 은 다음 과 같다.:주어진 방정식 은 수치 방법(예 를 들 어 용 격 법)을 통 해 풀 수 있다.쿠 타 법,뉴 트 라 마크 법 등.두 가지 알고리즘 의 유연성 을 비교 함으로써 이 알고리즘 의 충분 성 을 증명 하 였 다버팀목유한 원 분석 에서 얻 은 균열동력학 R4 에서 가정 한 알고리즘 에 따라 단원 시스템 을 최적화 시킨다.본 논문 의 임 무 는 들보 가 단위 력 의 작용 하에 균열 과 받 는 힘 의 서로 다른 단계 에서 균열 단면 의 플 렉 시 를 계산 하 는 것 이다.그림 3 은 권장 알고리즘 검사 결 과 를 보 여 준다.세 가지 결 과 를 비교 했다.유한 원 법(FEM)을 사용 하여 유 도 를 계산 하 다.유한 원 에서 균열 보 모형 의 직경 방향 유 도 를 계산 하 였 다.

전체 균열 각도 위치 범위 의 시스템;&#동력학 R4 를 사용 하여 유연성 을 계산 하고 FEM 을 사용 하여 초기 데 이 터 를 얻 습 니 다.모멘트 초기 값완전히 벌 어 진 균열 의 강도k 초기 화ee,  k초기 화

응.방정식 풀이(7)를 통 해 얻 을 수 있다.

k-초기 화

mh-균열 보 의 직경 방향 유연성

gc-상응하는 방향 에서 유한 원 법 을 사용 하여 계산 하 다.(10)법칙 에 따라 중간 각 균열 위치의 모멘트 강 도 는 최소 에서 최대 로 변 한다.

동력학 R4 를 사용 하여 유 도 를 계산 하 다.완전히 벌 어 진 균열 의 구부러짐 유도 초기 데이터

-k-ee

 초기 화 ,  k-응.  

초기 화단열 역학의 알고리즘 을 사용 하여[7,10]을 해석 하여 얻 었 다.모멘트 유연성 의 가치중간 각 균열 의 위 치 는 법칙(10)에 따라 최소 에서 최대 로 변화 한다.그림 3 한 바퀴 돌 때마다 균열 단면 에서 들보 의 유연성 변화 유한 원 계산 결과모델 은 동력학 R4 의 모델 과 매우 가깝다.초기 조건 에서 의 계산 결과 와 유한 원 결과 의 차 이 는 1%보다 적다.이 동시에 초기 강도 의 분석 계산 은 유한 원 의 계산 보다 훨씬 빠 르 고 필요 한 작업 시간 이 적어 서 실현 하기 쉽다.도착 하 다사용  균열 로터 의 기하학 적 매개 변수 균열 회전자 가 있 는 기하학 적 도형 을 선택 하여 알고리즘 작업 의 가장 좋 은 장점 을 나타 낸다.표 1 과 같다.중심 판 을 가 진 로터,받침 대 는 축 끝 에 놓 여 있다.② ②②②

 ② ②-②②②②②②②②②  ② 

②

 ② ② ②-②-②-② 

②

 ② ② ②