- 3장 - 동적 시스템의 시간역 모델링.

Slides:



Advertisements
Similar presentations
산업시스템분석 임성수 차수길 장연식 주혜림 7조7조.
Advertisements

이삼영 윤철호 김종국 박헌수 목련마을 화성빌리지 입구 예수 평강교회 앞 비타민 할인마트 옆 버스정류장 매화마을 관세청 입구 매화마을 건너편 탑마을 쪽문.
가. 3 일 취업완성 ! 성취반 _ 성공 취업 반 1 일차 시간 프로그램 내용비고 17:00 ~ 17:10 출석 및 일정 안내 17:10 ~ 18:10 취업 전략 및 채용 프로세스 이해 18:10 ~ 19:30 입사지원서 작성법 및 주요 기업 에세이 작성법 19:30 ~
구 분현존 무창계사 사육장 (1,000 평기준 ) 신개념 가금류사육장 (1,000 평기준 특허보유유럽에서 약 50 여년전 개발 2008 년 특허개발 ( 송백영농조합 ) 계사구조 별도 독립된 단층계사 500 평ⅹ 2 동 건축 많은 사육장면적 확보시 계사를 추가로 신축 500.
영남본부 현안사항 노동조합. -1- 현황 및 문제점 ( 단위 : 개 ) 가. 문제점 1) 현장구역과 회수업무 병행 - 공중전화유비보수, 모뎀회수 업무 2) 토, 일, 공휴일, 야간 근무 - KT 에서는 야간 및 휴일근무 부정적.
폭력. 폭력이란 무엇인가 우상의 눈물 물리적인 폭력 ( 최기표 ) VS 지능적인 폭력 ( 임형우, 담임선생님 )
일본주식시장의 신 고레가와긴조 투자전략 6 조 안승권. 신문수 발표자 : 신 문 수. 출 생 : 1897 효고현에서 출생 학 력 : 초등학교졸업, 사업가 1992 년 95 세 사망 유일한 자서전 1981 년 스미토모 금속광산 주식매매 200 억엔 벌다⇒ 일본 소득세 납세.
도민 기대수명 연장을 위한 도민 기대수명 연장을 위한. 함 양 군 보 건 소 추 진 배 경추 진 배 경 1 목적 및 목표 2 추진현황과 실적 3 향후 추진계획 4 목 차.
1 박 2 일 !!! 인천마장초등학교 유수아. 1 박 2 일 멤버 인기순 위 1 위 이승기 2 위 엄태웅 3 위 은지원 4 위 김종민, 이수근 ※인터넷에서 본것이기 때문에 사람에따라 서 다를 수 있다. ※
2009개정 중등 국어과 교육과정 울산광역시교육청 교육과정 컨설팅단 : 정일진.
효과적인 금연법 산재의료관리원 동해병원 건강관리센타.
우리나라 전통의 무술, 태권도 5학년 8반 김유승.
석관중앙교회 5남전도회 석 관 중 앙 교 회 회원 소식 통권 05-04호 발행일 : 2005년 04월 회 장 : 장진호 집사
남성에게 있어 조루증은 정말 부끄러우면서도 자신감을 잃게 만들게 합니다
한반도 통일 방식과 전망 10조 발표자료 발표: 변희섭 이종배 조영진.
달라지는 노동법 개정 내용 노무법인 正道 잠시나마… 주요 노동관계법 개정내용 3. 마무리 Contents
지역사회복지론 1조. 요양보호시설에 대해서 황성국 임재형 이동영
자살 사례 분석 경영학과 백승용 경영학부 하수정 경영학부 이은옥
박근혜 2년? 행복하셨습니까? 박근혜 2년은 이다 살의를 느끼게 한 세월 멈춰버린 대한민국 답답한 이년.
Ⅴ. 건설 기술과 환경.
* 그룹 상시 연락망 : 각사 조직도 기준 연락망으로 대체함
돼지가격 대표 기준 ‘탕박’변경 관련 설명자료
학습 주제 p 용해도 차이로 물질 분리하기.
막대의 영률 측정 1.실 험 목 적 Ewing 장치를 이용하여 휨에 의한 금속 막대의 영률을 계산한다.
I 문학의 개념과 역할 1. 문학의 개념 (1) 언어 예술로서의 문학 (2) 소통 활동으로서의 문학
직업재활실습 김은비.
4. 목적론적 윤리와 의무론적 윤리 01. 경험주의와 이성주의 01. 경험주의와 이성주의 02. 결과론적 윤리와 공리주의
수 계 소 화 설 비 론 - 소방시스템 기초이론 -.
열역학 Fundamentals of Thermodynamics(7/e) RICHARD E
순환&면역 6조 박아름 이명동 최제춘.
유체역학 (Fluid Mechanics)
운동에너지와 일에너지 에너지 물체의 상태와 관련된 스칼라 양 = 운동에너지 + 위치(퍼텐셜)에너지 + 열에너지 +…….
5장 기계공학 문제의 해결도구 및 해결절차.
에너지 운동량 방법: 일과 에너지법칙 1. 상자들이 초기속도 vo로 컨베이어 벨트로 운반되어 A에서 미끄러져서 B에서 떨어진다. μk= 0.40이고, 상자가 2.4m/s로 B점에서 떨어질 때 컨베이어 벨트의 속도를 구하라.
회전(Rotation) 각 변위 (Angular displacement) 각 속도 (Angular Velocity)
2장 수학적 기초 2.1 서 론 2.2 복소변수의 개념 2.3 미분방정식 2.4 라플라스변환
유독물 및 취급제한∙금지물질 관리자 교육 취급시설별 관리기준 2014 한강유역환경청 화학물질관리과.
개항기 조선과 동아시아 박 범 한국역사입문Ⅱ.
-9장- 디지털 제어시스템.
Chapter 8 운동량과 충격량, 충돌.
GS건설 하자보수 안내 GS건설에서 공용부 하자보수공사를 아래와 같이 실시합니다. 관리소의 공고
화재예방 안전교육 2013년 이현119안전센터장 김 용 래 작성
총인(Total Phosphorus ) 2조 김민구.
대한민국 No1. 고속열차 KTX 광고매체 제안서.
대구의 부도심 대구의 주요축 동대구 부도심 4조 강민석 / 박성균 / 최은지/ 황재현/김예지.
재생기술 regeneration technique
우리나라의 수자원 물 보기를 금같이 우리나라의 수자원 현황 우리나라의 수자원 이용 현황.
과학적 관리론 과 목 명 : 유아교육기관운영관리 교 수 명 : 홍성훈 교수님 학 과 : 유아교육과 학 번 :
문서의 제목 나눔고딕 45pt 작성자 | 소속팀 / 상위부서 | 이 문서는 나눔글꼴로 작성되었습니다. 설치하기.
지역의 자연 환경과 인문환경 조사 사회 1학년 1학기 Ⅰ.지역과 사회 탐구>1.지역사회의 지리적 환경(3/6
사도행전 13장 22절 말씀 –아멘 다 윗 을 왕 으 로 세 우 시 고 증 언 하 여 이 르 시 되 내 가 이 새 의 아 들
영률의 측정 Young’s modulus 이상호 수정일자 안선우 수정일자
CONTENTS Ⅰ. 대회목적 Ⅱ. 대회개요 Ⅲ. 대회요강 Ⅳ. 대회규정 Ⅴ. 운영계획 Ⅵ. 홍보계획 Ⅶ. 예산계획.
Ch8.기본적인 RL, RC 회로 자연응답, 강제응답, 시정수, 계단입력과 스위치 회로
평생 저축해도 강남 아파트 못산다 학 과 : 회계학과 1학년 B반 과 목 : 회계학원론 담당교수: 박성환 교수님
제목 입력란 그림 340mm x 170mm 글자 크기 30pt, 자간 52pt
악취저감 / 친환경 축산(농장)을 위한 시험보고서
문서의 제목 문서의 개요 작성자 이름 소속팀 소속팀 작성년월일
1월 교회학교 진급예배 및 성탄절 음악예배 찬 양 기 도 교 회 소 식 특 순 성 경 봉 독 말 씀 찬 양 축 도 인 도 자
인천경제자유구역 영종지구 영종하늘도시 개발사업
경찰행정과 세미나 결과를 공개해야한다. VS 비공개로 해야한다. 경찰의 근무성적평정 제도.
뜨거운 햇살을 받으며 양 손에 도시락 두 개를 들고, 콧 노래를 부르며, 시골 길을 걷고 있는 한 아이가 있었어요
Ⅲ. 선로전환기 청소 근거규정 및 점검요령.
유체역학 마이크로마노미터의 이론과 공식을 설명하라. 환경공학과 김기복.
우리나라의 최신무기 천 천 초 등 학 교 6학년 1반 16번 현 승 호.
제21강 교류회로와 전자기파 21.1 교류 회로에서의 저항 저항기: 전류와 전압이 같은 시간에 최대값 → 동일한 위상
근골격계 질환 예방교육.
위치 에너지 일과 중력 위치 에너지 물체끼리 주고받는 힘 때문에 생기는 에너지 물체의 배치(상대적 위치)에 따라 정해짐
2009개정 중등 국어과 교육과정.
쥐의 시간, 코끼리의 시간 (체중과 물질대사) 청원고등학교 배상기.
Presentation transcript:

- 3장 - 동적 시스템의 시간역 모델링

Contents 3.1 서론 3.2 동적 시스템 : 3.3 상태공간 모델식 3.4 비선형 시스템의 선형화 : 3.1 서론 3.2 동적 시스템 : 기계시스템, 유체시스템, 열시스템, 전기시스템 3.3 상태공간 모델식 3.4 비선형 시스템의 선형화 : Taylor 급수를 이용한 선형화, 피드백을 이용한 선형화 3.5 MATLAB을 이용한 동적 시스템의 시간역 모델링

3.1 서론  기술함수기법 이용 - 동적 시스템의 수학적 모델식 - 동적 시스템의 수학적 모델식 주파수역 모델식 : 선형 시불변 시스템에 적용, 전달함수로 표현 시간역 모델식 : 일반 시스템에 적용, 상태공간 모델식(상태방정식 + 출력방정식)으로 표현 - 대표적인 공학 시스템 병진 및 회전 기계시스템, 압축성 및 비압축성 유체시스템, 열시스템, 전기시스템 동적 시스템 간의 유사성 이용  다에너지역 시스템 조직적인 모델링 가능 동적 시스템의 기본 요소 : 에너지 발생, 저장, 발산, 변환 요소 - 비선형 시스템의 선형화 기법 소프트(soft) 비선형 시스템 : 비선형 요소의 입출력 관계 미분 가능한 시스템 Taylor 급수 이용(작동범위 작을 때) 피드백 선형화 이용(작동범위 클 때) 하드(hard) 비선형 시스템 : 비선형 요소의 입출력 관계 미분 불가능한 점 존재  기술함수기법 이용

3.2 동적 시스템 시스템의 각 요소가 가지고 있는 에너지 특성에 따라 에너지 발생, 저장, 발산, 변환 요소로 구분되며 3.2 동적 시스템 시스템의 각 요소가 가지고 있는 에너지 특성에 따라 에너지 발생, 저장, 발산, 변환 요소로 구분되며 이 요소들은 서로 접합되어 상호작용을 한다. 그림 3.1 동적 시스템의 에너지 흐름선도

파워(에너지) 보존 법칙 (동적 시스템의 기본 방정식) (3.1) 여기서 : 모든 에너지 소스(source)로부터 시스템에 가해진 파워(power) : 에너지 저장요소에 저장된 에너지의 시간변화율 : 에너지 발산요소에서 소멸된 파워 또는 (3.2) 여기서 P : 작용력(effort)변수 e와 흐름(flow)변수 f 의 곱으로 정의 (3.3)

표 3.1 대표적인 공학 시스템에 대한 작용력변수와 흐름변수

- 동적 시스템의 물리 변수 파워변수 : 작용력변수 e 와 흐름변수 f 에너지변수 : 일반화된 운동량 p 와 일반화된 변위 q - 동적 시스템의 기본 방정식 또는 여기서 C : 커패시턴스(위치에너지 저장요소) 여기서 I : 이너턴스(운동에너지 저장요소) 여기서 R : 저항(에너지 발산요소) 또는 (위치에너지) 또는 (운동에너지)

시스템의 접합 구조 - 0-접합 : 작용력변수가 같은 경우 (3.8) - 1-접합 : 흐름변수가 같은 경우 (3.9) - 0-접합 : 작용력변수가 같은 경우 (3.8) - 1-접합 : 흐름변수가 같은 경우 (3.9) 표 3.2 대표적인 공학 시스템에서 접합 구조의 물리적 의미

일반화된 커패시턴스 - 작용력변수 e가 일반화된 변위 q의 함수로 표현되는 위치에너지 저장요소 (3.10) - 위치에너지 저장요소의 예 스프링 : 토션바 : 어큐뮬레이터 : 커패시터 : - 일반적인 C 형태의 에너지 저장요소의 파워변수 관계식 또는 (3.12)

일반화된 이너턴스 - 흐름변수 f 가 일반화된 운동량 p의 함수로 표현되는 운동에너지 저장요소 (3.13) - 운동에너지 저장요소의 예 질량 : 회전관성 : 유체관성 : 인덕턴스 : - 일반적인 I 형태의 에너지 저장요소의 파워변수 관계식 또는 (3.15)

일반화된 저항 - 작용력변수 e가 흐름변수 f의 함수로 표현되는 에너지 발산요소 (3.16) - 선형 시스템에서의 에너지 발산요소 (3.17) - 에너지 발산요소의 예 감쇠 : 마찰 : 저항 :

에너지 변환요소 - 에너지 손실이 없다고 가정하면 에너지 변환요소의 입력파워 Pi 와 출력파워 Po 는 같다. (3.18) - 에너지 변환방법에 따른 에너지 변환 요소의 분류 트랜스포머(transformer) : (3.19) 여기서 Kt : 트랜스포머계수 자이레이터(gyrator) : (3.20) 여기서 Kg : 자이레이터계수

표 3.3 대표적인 공학 시스템에 대한 에너지 발생, 저장 및 발산 요소

표 3.4 대표적인 공학 시스템에 대한 에너지 변환요소

상사(similarity) 시스템 개념 기계, 유체, 열, 전기 시스템 등 플랜트의 종류가 다르더라도 각 시스템이 상사 시스템이면 수학적 모델식이 전적으로 같다. 사용하는 물리변수 및 파라미터 등의 용어만 다르다. 예) 기계시스템 : 힘, 속도, 변위, 운동량, 스프링상수, 질량, 감쇠계수 전기시스템 : 전압, 전류, 전하량, 자속쇄교, 커패시턴스, 인덕턴스, 저항 상사 시스템 개념을 적용하면 다에너지역 동적 시스템의 상태공간 모델식 또는 전달함수를 구하는 과정을 체계적으로 수행할 수 있다. 시스템 각 요소간의 접합 조건 및 파워 보존 법칙 적용 (2) 시스템의 각 요소방정식 정리

3.2.1 기계시스템 기계적(mechanical) : 고체 물체에 작용하는 힘이나 토크 그리고 물체의 운동을 해석함을 의미함 표 3.5 기계시스템의 기본적인 물리변수

표 3.6 대표적인 위치에너지 저장요소들에 대한 스프링상수

표 3.7 집중질량 등가 기계시스템 표 3.8 대표적인 물체의 회전관성

표 3.9 대표적인 기계적 에너지 발산요소들에 대한 감쇠마찰계수

질량-스프링-댐퍼 시스템의 동적 방정식 유도 예제 3.2 질량-스프링-댐퍼 시스템의 동적 방정식 유도 그림 3.12 질량-스프링-댐퍼 시스템 - 접합 구조를 고려한 파워 보존 법칙 - 에너지 저장요소들의 동적 방정식 또는 - 에너지 발산요소의 요소방정식 - 3차 시스템(상태변수 : 질량의 속도 , 그리고 스프링 힘 )

3.2.2 유체시스템 - 유체시스템의 분류 유공기압시스템 : 정적 압력만 고려 일반적인 유체시스템 : 정적 및 동적 압력 모두 고려 유압시스템 견고하고 고속이며 큰 동력을 발생시킬 수 있음 공작기계, 굴삭기, 동력전달장치 그리고 항공기 제어용 서보기구 등에 이용 공기압시스템 전기 모터보다 큰 파워를 얻을 수 있음 화재로부터 위험이 없음 응답속도가 유압시스템보다 느림 작동유체의 압축성 때문에 모델링이 복잡함 그림 3.13 일반적인 유압시스템

표 3.10 유공기압시스템의 기본적인 물리변수

(1) 위치에너지 저장요소 일정 면적 유체 탱크의 유체 커패시터 - 일정 면적 유체 탱크의 기본 방정식 (3.40) (3.41a) (3.41b) 여기서 P : 유체 탱크 바닥에서의 압력 V : 유체 탱크의 체적 (3.42) 여기서 g : 중력가속도상수 A : 유체 탱크의 단면적 그림 3.14(a) 일정 면적 유체 탱크의 유체 커패시터 - 유체 탱크의 커패시턴스 C (3.43)

강체 파이프의 유체 커패시터 dP : 유체 체적 V가 dV 만큼 증가할 때의 - 강체 파이프 내의 유체의 기본 방정식 (3.44) 여기서 B: 체적탄성계수(bulk modulus) dP : 유체 체적 V가 dV 만큼 증가할 때의 압력변화 또는 (3.45) 여기서 : 원래의 유체 체적 ( ) 그림 3.14(b) 강체 파이프의 유체 커패시터 - 강체 파이프의 유체 커패시턴스 C (3.47)

유연 파이프의 유체 커패시터  - 압력변화 P에 의해 파이프에서 발생되는 응력 (3.48) - 원주 방향의 스트레인 (3.49) 여기서 E : 탄성계수 - 반지름 증가에 의한 체적변화량 V (3.50) - 미소 스트레인 그림 3.14(c) 유연 파이프의 유체 커패시터 (3.51) - 유연 파이프의 커패시턴스  (3.52) (3.59) 여기서 공칭체적

압축가스 어큐뮬레이터 - 등 엔트로피 과정으로 가정 일정 (3.60) 여기서 첨자 o : 초기값, Υ : 일정한 체적에서의 비열 에 대한 일정한 압력에서의 비열 의 비, 즉 (3.61) 여기서  ,  로 치환하면 (3.62) 그림 3.15 압축가스 어큐뮬레이터 - 압축가스의 커패시턴스 (3.63)

(2) 운동에너지 저장요소 - 유체관성 : 유체시스템에서 운동에너지를 저장하는 요소 - 유체관성 : 유체시스템에서 운동에너지를 저장하는 요소 - 유체관성 I는 작용력변수 P와 흐름변수 Q 사이의 관계식으로부터 유도 (3.68) (3.69) (3.70) - 유체관성 I 계산 (3.71) 여기서 : 유체의 가속도 그림 3.17 일정 면적 파이프 (3.72) (3.73)

(3) 에너지 발산요소 - 유체저항식 층류(laminar) 유체 경우( ) (3.75) 여기서 (3.77) 난류(turbulent) 유체 경우( ) (3.79) 그림 3.18 유체저항: (a) 다공성 마개, (b) 모세관 튜브, (c) 긴 파이프, (d) 오리피스, (e) 밸브 여기서 (3.80)

(4) 에너지 변환요소   트랜스포머 형태의 에너지 변환요소 - 용적형(positive-displacement type) 펌프 (3.90) - 출력 체적 와 축의 각변위 사이의 관계식 (3.94)  (3.95) - 파워 보존 법칙 (3.96)  (3.97) 트랜스포머 형태의 에너지 변환요소 그림 3.23 흐름소스로 모델링된 일정-변위 펌프

2차 시스템 (상태변수 : 펌프에서 공급되는 유량 Q1, 유체저항 통과 후의 압력 P3 ) 예제 3.3 유압시스템의 동적 방정식 유도 그림 3.25 유연하고 긴 파이프를 갖는 유압시스템 유체관성 I와 유체커패시터 C에 대한 동적 방정식 유체저항 R과 부하저항 RL에 대한 요소방정식 2차 시스템 (상태변수 : 펌프에서 공급되는 유량 Q1, 유체저항 통과 후의 압력 P3 )

3.2.3 열시스템 - 열시스템의 특징 열 에너지와 온도효과를 고려해야 함 운동에너지 저장요소가 필요 없음 - 열 시스템의 기본 방정식(Gibbs 방정식) (3.98) 여기서 u, s, v : 각각 단위 질량당 내부에너지, 엔트로피, 체적 표 3.11 열시스템의 물리적인 기본변수

열전달 시스템 - 엔트로피 개념은 추상적이므로 흐름변수로 엔트로피유량 대신 열(에너지)유량 Q를 사용 - 열역학적 온도는 절대온도뿐만 아니라 섭씨, 화씨 온도도 사용 - 표 3.12는 물리적 보다는 수학적인 유사성(analogy)을 갖는다. 표 3.12 열전달 시스템의 기본적인 물리변수

열전도에 의한 열저항 열대류에 의한 열저항  전도 열저항 Rk와 대류 열저항 Rc (3.101) 여기서 k : 열전도도, A : 열전달 면적, l : 두께 열대류에 의한 열저항 (3.102) 여기서 h : 대류 열전달계수, A : 열전달 면적, T1 , T2 : 벽과 유체의 온도  전도 열저항 Rk와 대류 열저항 Rc (3.104) (3.105) 그림 3.27 열전달 및 대류에 의한 열저항 요소

복사에 의한 열저항 - Stefan-Boltzmann 법칙 (3.107) 여기서 (Stefan-Boltzmann 상수) 유효 방사율(emissivity) 형상계수 열전달 면적 - 식 (3.107)을 작동점 근처에서 선형화 (3.109) 여기서 (3.110)

열 커패시턴스 - 열 커패시턴스의 기본 방정식 여기서 c : 비열(specific heat) 그림 3.28 열 커패시턴스 - 비열 c는 온도에 따라 변하지만, 온도의 변화량 T 가 작다고 가정 또는 (3.114) (3.115) (3.116) 여기서 (3.117)

- 열 커패시턴스 C와 열저항 R를 고려한 열시스템에 대한 동적 모델식 예제 3.4 열시스템에 대한 동적 모델식 유도 그림 3.29 기본적인 열시스템의 예 단열재와 재료 블록 사이의 열저항 그리고 단열재와 주변의 대기온도 T0 사이의 열저항을 합한 등가 열저항을 R로 나타냄 - 대기온도는 T0 로 일정하다고 가정 - 열 커패시턴스 C와 열저항 R를 고려한 열시스템에 대한 동적 모델식 또는 - 1차 시스템(상태변수 : 온도 T )

3.2.4 전기시스템 표 3.13 전기시스템의 기본적인 물리변수

(1) 위치에너지 저장요소 - 위치에너지 저장요소의 기본 방정식 (3.118) 또는 (3.119) - 커패시턴스 C의 단위 : F(farad) = 쿨롬/전압(C/V) 그림 3.30 커패시턴스 C에 대한 기호

(2) 운동에너지 저장요소 - 운동에너지 저장요소의 기본 방정식 (3.120) 또는 (3.121) - 유도 효과는 자기 인덕턴스와 상호 인덕턴스로 구분 인덕턴스: 단일 코일이 코일전류에 의하여 만들어진 자기장이 코일 자체와 관련될 때 발생 인덕턴스: 두 인덕터가 한 전기회로 안에 포함되어 있을 때, 서로 미치는 자기장의 영향 인덕턴스의 단위: H(henry) = 웨버/암페어(Wb/A) 자기 상호 그림 3.31 인덕턴스 L에 대한 회로 기호 그림 3.32 상호 연성을 갖는 인덕터

(3) 에너지 발산요소 (4) 에너지 발생 및 변환요소 - 에너지 발산요소의 기본 방정식 (3.126) - 저항 R의 단위: - 전류 소스 : 시간의 함수로 전류를 발생시키는 에너지 발생요소 - 전압 소스 : 시간의 함수로 전압을 발생시키는 에너지 발생요소 전류소스 전압소스 그림 3.36 전류 및 전압소스에 대한 회로 기호

- 전기적 에너지 변환요소 트랜스포머 형태의 에너지 변환식 (3.127) 여기서 n은 트랜스포머계수 자이레이터 형태의 에너지 변환식 (3.128) 여기서 r는 자이레이터계수 (a) 트랜스포머 (b) 자이레이터 그림 3.37 전기적 에너지 변환요소

- 2차 시스템(상태변수 : 커패시터에 걸리는 전압 Vc, 인덕터 L에 흐르는 전류 iL) 예제 3.5 전기시스템에 대한 동적 모델식 유도 그림 3.44 직렬 및 병렬로 연결된 전기시스템 에너지 저장요소인 커패시터 C 와 인덕터 L에 대한 동적 방정식 저항 R1과 R2의 요소방정식 - 2차 시스템(상태변수 : 커패시터에 걸리는 전압 Vc, 인덕터 L에 흐르는 전류 iL)

3.3 상태공간 모델식 - 상태공간 : 어떤 지정된 시간에서 시스템의 상태를 완전하게 서술하는 데 필요한 임의의 - 상태공간 : 어떤 지정된 시간에서 시스템의 상태를 완전하게 서술하는 데 필요한 임의의 최소 변수의 집합 - 상태변수 : 시스템의 상태를 정의할 수 있는 임의의 최소 변수 - 상태벡터 : 상태변수를 벡터 형태로 표현 - 상태공간 모델식 : 상태방정식과 출력방정식(선형 시불변 다변수 시스템의 경우) (3.140) (3.141) 여기서 은 상태벡터, 은 제어입력벡터, 는 출력벡터 - 상태변수 선정 시스템 내에 있는 독립된 에너지 저장요소의 파워변수[위치에너지 저장요소의 작용력 변수 e 와 운동에너지 저장요소의 흐름변수 f ] 또는 Lagrange 변수[위치에너지 저장 요소의 일반화된 변위변수 와 운동에너지 저장요소의 일반화된 속도 ] 시스템 내에 관심 있는 변수, 즉 시스템의 출력으로 운동에너지 저장요소의 흐름변수 를 적분한 일반화된 변수 가 포함되어 있을 때, 이를 상태변수로 추가

- 질량-스프링-댐퍼 기계시스템의 상태공간모델식 유도 그림 3.46 질량-스프링-댐퍼 기계시스템 - Newton의 제2법칙 적용 (3.142) 상태변수 선정 : ① Lagrange변수 : 위치에너지 저장요소인 스프링의 변위 x(t)와 운동에너지 저장요소인 질량의 속도 v(t) ② 파워변수 : 위치에너지 저장요소인 스프링의 작용력 Fs(t)와

- 상태방정식 및 출력방정식 유도 ① Lagrange 상태변수를 사용한 경우 상태방정식 (3.145) 또는 (3.146) 여기서 출력방정식 (3.147) 여기서

② 파워 상태변수를 사용한 경우 상태방정식 (3.150) 또는 (3.151) 여기서 출력방정식 (3.152) 여기서

상태공간 모델식으로부터 전달함수 유도 - 상태공간 모델식에 대해 Laplace 변환 수행 (3.156) 여기서 (3.157) - 파워변수에 의한 상태공간 모델식으로부터 전달함수 유도 (3.159) - 기계시스템의 동적 방정식 (3.142)에 대해 Laplace 변환 수행 (3.160) 상태변수의 선정이 다르더라도 전달함수는 동일하게 유도되고, 시스템행렬 A의 고윳값은 동일함

발전기: 자이레이터 형태의 에너지 변환요소 Kg 종합문제 3.1 수력발전시스템의 상태공간 모델식 유도 - 입력: 비에 의해 공급되는 유량 Q - 출력: 저수량 VR, V - 수력발전시스템의 구성요소 저수지: 위치에너지 저장요소 C 밸브: 저항요소 Rv 터빈: 트랜스포머 형태의 에너지 변환요소 Kt 축: 운동에너지 저장요소 J, 저항요소 R 발전기: 자이레이터 형태의 에너지 변환요소 Kg 그림 3.52 수력발전시스템 에너지 변환 : 유체의 위치에너지 기계에너지(회전) 전기에너지 터빈 발전기 ① 저수지와 축의 동적 방정식 여기서 P : 저수지 수압, : 축의 각속도, Q : 수로에서의 유량, J : 축토크  

② 에너지 발산요소(밸브저항 Rv, 축의 점성저항 cT, 전기저항 R)의 요소방정식 여기서 Pt : 터빈 측 수압, t : 터빈 측 토크, g : 발전기 측 토크, V : 발전기 전압, i : 전기저항에 흐르는 전류 ③ 에너지 변환요소(터빈 Kt , 발전기 Kg )의 요소방정식 여기서 Kt : 터빈의 트랜스포머계수, Kg: 발전기의 자이레이터계수

④ 각 요소방정식들을 정리하여 상태방정식 유도 또는 여기서 ⑤ 출력방정식 유도 여기서

3.4 비선형 시스템의 선형화 선형화 기법 비선형 시스템에 대해서도 가능하면 선형화하여 조직적인 선형 제어기법 적용 - 선형 시스템 중첩 및 동차성의 원리를 만족함 시스템의 특성을 파악할 때 시스템 입력의 크기를 고려할 필요 없음 - 비선형 시스템 중첩 및 동차성의 원리를 만족하지 않음 시스템 특성이 시스템 입력의 크기에 따라 변함 해석 및 설계기법이 어렵고 사용범위가 제한됨 비선형 시스템에 대해서도 가능하면 선형화하여 조직적인 선형 제어기법 적용 선형화 기법 Taylor 급수를 이용한 선형화 기법 피드백을 이용한 선형화 기법 또는 입출력 선형화 기법

3.4.1 Taylor 급수를 이용한 선형화 - 입출력 관계에서 미분 불가능한 점이 존재하지 않는 비선형 시스템에 적용 (3.161) 여기서 는 출력 가 입력 의 비선형 함수임을 나타냄 - 작동점 근처에서 Taylor 급수 이용 (3.162) - 편차 가 작다고 가정 여기서 그림 3.48 비선형 요소에 대한 입출력 관계

작동점 vo근처에서 선형화한 상태방정식을 이용한 전달함수 G(s)=v(s)/f(s) 구하기 예제 3.6 (3.167) - 평형상태에서의 입력 fo (3.168) - 상태변수 v(t)와 입력 f(t)에 작은 섭동(perturbation)을 가했다고 가정 (3.169) - 식 (3.169)로부터 (3.170) - 식 (3.169)와 식 (3.170)을 식 (3.167)에 대입 또는 (3.171)

- 식 (3.168)을 식 (3.171)에 대입 - 시스템이 작동점 근처에서 운동한다고 가정하면 항은 무시가능 (3.172) - 초기조건이 0인 가정 아래서 식 (3.172)에 Laplace 변환 수행 여기서 f(s) : 입력 f(t)에 대한 Laplace 변환함수 v(s) : 출력 v(t)에 대한 Laplace 변환함수 - 선형화된 시스템의 전달함수 G(s)

m : 구의 질량, (x, i) : 자기력에 의해 발생되는 부상력, x : 구의 부상높이, 예제 3.7 한 개의 전자석이 구 베어링을 부상시키는 구 부상기의 선형화된 모델 유도 그림 3.49 구 부상기의 모델 및 실험적으로 구한 자기부상력 곡선 - Newton의 제2법칙을 이용한 구에 대한 비선형 운동방정식 유도 여기서 m : 구의 질량, (x, i) : 자기력에 의해 발생되는 부상력, x : 구의 부상높이, i : 코일에 흐르는 전류, 그리고 g : 중력가속도상수 9.81 m/sec2

평형상태에서의 부상속도 x2o= 0, 자기부상력 0 = mg - 상태방정식 여기서 제어입력 u : 전자석에 흐르는 전류 i - 평형상태에서의 상태벡터 xo=[x1o x2o]T이고, 제어입력 uo일 때의 평형조건식 평형상태에서의 부상속도 x2o= 0, 자기부상력 0 = mg - 선형화를 위해 상태변수와 제어입력을 평형값(x1o , x2o , u o)과 섭동량( )으로 표시

- 선형화된 상태방정식 - 평형상태에서 x2o = 0 , 0 = mg - 선형화된 상태방정식 여기서 고정된 x1 = 2.5 mm에서 전류의 변화량에 따른 자기부상력의 변화량 그림 3.49에 표시된 기울기 K(약 13 N/m) - 선형화된 상태방정식 여기서

3.4.2 피드백을 이용한 선형화 예제 3.8 유체 탱크 시스템에서 지정된 수위 hd가 유지되도록 하는 적절한 제어입력 u(t) 선정 그림 3.50 유체 탱크 시스템 (가정) 유체 탱크 및 파이프 단면적은 각각 A와 a로 일정 - 유체 탱크 시스템의 동적 모델식 (3.179) 수위 가 작동점 ho 근처에서 작동 시(즉, 섭동량 δh가 충분히 작을 때) Taylor 급수를 이용하여 선형화 가능 (3.180)

h(t)의 크기에 상관없이 항상 선형성 유지  피드백 선형화 기법 적용 (3.181) 여기서 v(t) : 피드백 선형화된 시스템의 등가 제어입력 - 피드백 선형화된 모델식 (3.182) h(t)의 크기에 상관없이 항상 선형성 유지 그림 3.51 피드백 선형화를 이용한 제어시스템의 구조

3.5 MATLAB을 이용한 동적 시스템의 시간역 모델링 예제 3.9 상태공간 모델식을 전달함수로 표현하고, 시스템의 극점 구하기 MATLAB 프로그램 3.1 - 전달함수 - G(s)의 극점

다변수 시스템의 상태공간 모델식을 전달함수행렬 G(s)로 표현 예제 3.10 다변수 시스템의 상태공간 모델식을 전달함수행렬 G(s)로 표현 MATLAB 프로그램 3.2 - 전달함수행렬

전달함수 G(s)를 상태공간 모델식으로 표현 예제 3.11 전달함수 G(s)를 상태공간 모델식으로 표현 - 상태공간 모델식 MATLAB 프로그램 3.3 여기서