KR101901352B1 - Predicting System For Seismic Observation - Google Patents
Predicting System For Seismic Observation Download PDFInfo
- Publication number
- KR101901352B1 KR101901352B1 KR1020170083683A KR20170083683A KR101901352B1 KR 101901352 B1 KR101901352 B1 KR 101901352B1 KR 1020170083683 A KR1020170083683 A KR 1020170083683A KR 20170083683 A KR20170083683 A KR 20170083683A KR 101901352 B1 KR101901352 B1 KR 101901352B1
- Authority
- KR
- South Korea
- Prior art keywords
- value
- test structure
- measurement value
- acceleration sensor
- model
- Prior art date
Links
- 238000012360 testing method Methods 0.000 claims abstract description 63
- 238000005259 measurement Methods 0.000 claims abstract description 27
- 230000001133 acceleration Effects 0.000 claims abstract description 18
- 238000004088 simulation Methods 0.000 claims description 6
- 238000000034 method Methods 0.000 description 9
- 230000003111 delayed effect Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G01V1/01—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/008—Earthquake measurement or prediction
-
- E—FIXED CONSTRUCTIONS
- E04—BUILDING
- E04B—GENERAL BUILDING CONSTRUCTIONS; WALLS, e.g. PARTITIONS; ROOFS; FLOORS; CEILINGS; INSULATION OR OTHER PROTECTION OF BUILDINGS
- E04B1/00—Constructions in general; Structures which are not restricted either to walls, e.g. partitions, or floors or ceilings or roofs
- E04B1/62—Insulation or other protection; Elements or use of specified material therefor
- E04B1/92—Protection against other undesired influences or dangers
- E04B1/98—Protection against other undesired influences or dangers against vibrations or shocks; against mechanical destruction, e.g. by air-raids
-
- E—FIXED CONSTRUCTIONS
- E04—BUILDING
- E04H—BUILDINGS OR LIKE STRUCTURES FOR PARTICULAR PURPOSES; SWIMMING OR SPLASH BATHS OR POOLS; MASTS; FENCING; TENTS OR CANOPIES, IN GENERAL
- E04H9/00—Buildings, groups of buildings or shelters adapted to withstand or provide protection against abnormal external influences, e.g. war-like action, earthquake or extreme climate
- E04H9/02—Buildings, groups of buildings or shelters adapted to withstand or provide protection against abnormal external influences, e.g. war-like action, earthquake or extreme climate withstanding earthquake or sinking of ground
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
- F16F15/00—Suppression of vibrations in systems; Means or arrangements for avoiding or reducing out-of-balance forces, e.g. due to motion
- F16F15/02—Suppression of vibrations of non-rotating, e.g. reciprocating systems; Suppression of vibrations of rotating systems by use of members not moving with the rotating systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D1/00—Measuring arrangements giving results other than momentary value of variable, of general application
- G01D1/02—Measuring arrangements giving results other than momentary value of variable, of general application giving mean values, e.g. root means square values
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D1/00—Measuring arrangements giving results other than momentary value of variable, of general application
- G01D1/16—Measuring arrangements giving results other than momentary value of variable, of general application giving a value which is a function of two or more values, e.g. product or ratio
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P15/00—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration
- G01P15/006—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration by making use of fluid seismic masses
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G06N3/0481—
Abstract
Description
본 개시는 지진 관측에 대한 예측 시스템 및 방법에 관한 것으로 보다 상세하게는 테스트 구조물을 이용한 지진 관측에 대한 예측 시스템 및 방법에 관한 것이다.This disclosure relates to a system and method for predicting seismic observations, and more particularly to a system and method for predicting seismic observations using test structures.
고층빌딩 또는 브릿지 등의 경우 지진이나 바람을 견딜 수 있는 내진 설계가 필수적이다. 고층빌딩을 예로 들면, 동조 질량 댐퍼(TMD, tuned mass damper)와 같은 건물의 진동을 억제하기 위한 기기가 배치되는 것이 일반적이다.In the case of high-rise buildings or bridges, seismic design that can withstand earthquakes and wind is essential. Taking high-rise buildings as an example, devices for suppressing vibration of buildings such as tuned mass damper (TMD) are generally disposed.
그리고, 건축 공사가 시행되기 이전에 건물의 내진을 시험하기 위한 테스트 구조물이 이용될 수 있다. 테스트 구조물에 인위적인 충격, 진동이 가해지는 경우에 모니터링된 테스트 구조물의 거동(behavior) 데이터를 수집하여 건출물의 설계 및 시공시에 유용하게 이용될 수 있다.Test structures can be used to test the building's earthquake prior to the construction work. It can be useful in designing and constructing dry matter by collecting the behavior data of the monitored test structure when artificial impact and vibration are applied to the test structure.
이런 경우, 테스트 구조물의 거동 및 나아가 건축물의 거동을 용이하게 예측하는 기술의 대두가 요청된다. In this case, there is a demand for a technique for easily predicting the behavior of the test structure and further the behavior of the building.
한편, 상기와 같은 정보는 본 발명의 이해를 돕기 위한 백그라운드(background) 정보로서만 제시될 뿐이다. 상기 내용 중 어느 것이라도 본 발명에 관한 종래 기술로서 적용 가능할지 여부에 관해, 어떤 결정도 이루어지지 않았고, 또한 어떤 주장도 이루어지지 않는다On the other hand, the above information is only presented as background information to help understand the present invention. No decision has been made as to whether any of the above contents can be applied as the prior art related to the present invention and no claim is made
본 발명은 상술한 문제점을 해결하기 위해 안출된 것으로, 본 발명의 일 실시 예는 테스트 구조물을 이용한 지진 관측에 대한 예측 시스템 및 방법을 제안한다.The present invention has been made to solve the above-mentioned problems, and one embodiment of the present invention proposes a prediction system and method for a seismic observation using a test structure.
본 발명에서 이루고자 하는 기술적 과제들은 이상에서 언급한 기술적 과제들로 제한되지 않으며, 언급하지 않은 또 다른 기술적 과제들은 아래의 기재로부터 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 명확하게 이해될 수 있을 것이다.It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory and are not restrictive of the invention, unless further departing from the spirit and scope of the invention as defined by the appended claims. It will be possible.
상기한 과제를 실현하기 위한 본 발명의 일 실시 예와 관련된 지진 관측에 대한 예측 시스템은 복수의 플로어(floor)를 포함하고, 상기 복수의 플로어 중 최상층에 가속도센서가 구비된 테스트 구조물; 상기 테스트 구조물에 발생되는 진동을 감쇄시키는 TMD(tuned mass damper); 상기 테스트 구조물을 지지하면서 상기 테스트 구조물에 지진 신호를 트리거하는 진동 테이블; 및 상기 진동 테이블을 제어하는 제어부를 포함하며, 상기 제어부는, 상기 가속도센서를 통해 측정되는 실제 측정값 및 적어도 하나의 예측 모델을 이용하여 상기 실제 측정값을 예측하는 예측값을 산출하고, 소정의 조건을 만족하는 예측 모델을 선택할 수 있다.In order to achieve the above object, there is provided a system for predicting seismic observations according to an embodiment of the present invention, comprising: a test structure including a plurality of floors and having an acceleration sensor on the top of the plurality of floors; A tuned mass damper (TMD) for attenuating vibration generated in the test structure; A vibration table for supporting the test structure and triggering an earthquake signal on the test structure; And a control unit for controlling the vibration table. The control unit calculates a predicted value for predicting the actual measurement value using the actual measurement value measured through the acceleration sensor and at least one prediction model, Can be selected.
본 발명의 다양한 실시예에 따르면 아래와 같은 효과가 도출될 수 있다.According to various embodiments of the present invention, the following effects can be obtained.
첫째, 테스트 구조물을 이용한 지진 관측에 대한 예측 시스템이 제공됨으로써, 보다 정확하고 신속하게 테스트 구조물 및 실제 건축물의 거동이 예측될 수 있다.First, by providing a prediction system for seismic observations using test structures, the behavior of test structures and real buildings can be predicted more accurately and quickly.
둘째, 테스트 구조물의 거동을 측정에만 의지하지 않고 수학적으로 정확하게 예측하는 예측 시스템이 제공됨으로써, 장치 효율성 및 비용 절감성이 향상될 수 있다.Second, device efficiency and cost savings can be improved by providing a prediction system that accurately predicts the behavior of a test structure mathematically without resorting to measurement alone.
셋째, 이미 건축된 구조물의 내진 성능을 평가하는데 상기 시스템이 이용될 수 있는 바, 재난 예측 및 방지에 도움이 될 수 있다.Third, the system can be used to evaluate the seismic performance of already built structures, which can help predict and prevent disasters.
넷째, 고층빌딩 또는 긴 다리(long span bridge)의 안전 성능을 판단할 때, 상기 시스템이 이용되는 경우 다양한 상황을 시험해볼 수 있다.Fourth, when judging the safety performance of a high-rise building or a long span bridge, various situations can be tested when the system is used.
본 발명에서 얻은 수 있는 효과는 이상에서 언급한 효과들로 제한되지 않으며, 언급하지 않은 또 다른 효과들은 아래의 기재로부터 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 명확하게 이해될 수 있을 것이다.The effects obtained by the present invention are not limited to the above-mentioned effects, and other effects not mentioned can be clearly understood by those skilled in the art from the following description will be.
도 1은 실시 예에 따른 테스트 구조물을 이용한 지진 관측에 대한 예측 시스템을 나타낸다.
도 2는 실시 예에 따른 테스트 구조물에 배치된 TMD를 설명하기 위한 도면이다.
도 3은 실시 예에 따른 지진 관측에 대한 예측 시스템의 구성을 나타내는 블록도이다.
도 4(a) 내지 도 4(c)는 실시 예에 따른 NARX 모델을 나타낸다.
도 5는 실시 예에 따른 다양한 예측 모델을 나타내는 테이블이다.
도 6은 실시 예에 따른 진동 측정값을 특정 예측 모델에 기초하여 예측하는 그래프이다.
도 7은 실시 예에 따른 예측 시스템의 구동 방법을 나타낸다.Figure 1 shows a prediction system for seismic observations using a test structure according to an embodiment.
2 is a view for explaining a TMD disposed in a test structure according to an embodiment.
3 is a block diagram showing a configuration of a prediction system for seismic observation according to an embodiment.
4 (a) to 4 (c) show a NARX model according to an embodiment.
5 is a table showing various prediction models according to the embodiment.
6 is a graph for predicting vibration measurement values according to an embodiment based on a specific prediction model.
7 shows a method of driving the prediction system according to the embodiment.
첨부되는 도면들을 참조하는 하기의 상세한 설명은 청구항들 및 청구항들의 균등들로 정의되는 본 개시의 다양한 실시 예들을 포괄적으로 이해하는데 있어 도움을 줄 것이다. 하기의 상세한 설명은 그 이해를 위해 다양한 특정 구체 사항들을 포함하지만, 이는 단순히 예로서만 간주될 것이다. 따라서, 해당 기술 분야의 당업자는 여기에서 설명되는 다양한 실시 예들의 다양한 변경들 및 수정들이 본 개시의 범위 및 사상으로부터 벗어남이 없이 이루어질 수 있다는 것을 인식할 것이다. 또한, 공지의 기능들 및 구성들에 대한 설명은 명료성 및 간결성을 위해 생략될 수 있다.The following detailed description, which refers to the accompanying drawings, will serve to provide a comprehensive understanding of the various embodiments of the present disclosure, which are defined by the claims and the equivalents of the claims. The following detailed description includes various specific details for the sake of understanding, but will be considered as exemplary only. Accordingly, those skilled in the art will recognize that various changes and modifications of the various embodiments described herein may be made without departing from the scope and spirit of this disclosure. Furthermore, the descriptions of well-known functions and constructions may be omitted for clarity and conciseness.
도 1은 본 발명의 일 실시 예에 따른 테스트 구조물(10)을 이용한 지진 관측에 대한 예측 시스템(100)을 나타낸다.1 shows a
상기 시스템(100)은 테스트 구조물(10), 진동 테이블(110), 진동 컨트롤 박스(120), 제어 컴퓨터(PC)를 포함할 수 있다. The
테스트 구조물(10)은 세 개의 층(11-1~11-3)을 포함하고, 각 층마다(11-1~11-3) 가속도 센서(accelerometer, 130-1~130-3)가 배치될 수 있다. 다만, 가속도 센서 중 가장 상층(11-1)의 가속도 센서(130-1)의 진동이 테스트 구조물(10)의 배치상 가장 크게 관측되므로 가장 중요한 측정 센서로 사용될 수 있으나, 실시 예는 이에 국한되지 않는다.The
여기서, 테스트 구조물(10)의 제원을 예시적으로 살펴보면, 각 플로어(11)는 평평하고, 테스트 구조물(10)의 높이는 930 mm이고, 각 플로어(11-1~11-3)은 300 mm X 300 mm 의 사이즈에 3.4 Kg를 가지며, 각 층의 기둥(15-1~15-4)는 0.2 Kg 일 수 있다. 또한, 테스트 구조물(10)은 진동 테이블(110)과 4개의 접점(13-1~13-4)을 포함한다. 상기 테스트 구조물(10)은 실제 건축하고자 하는 건축물을 시뮬레이션한 구조물일 수 있으나, 실시 예는 이에 국한되지 않는다.The
테스트 구조물(10)은 진동 테이블(110) 상에 배치될 수 있다. 상기 진동 테이블(110)은 힘(F1, F2)에 의해 좌우로 움직이게 된다. 상기 힘은 제어 컴퓨터(PC)의 제어에 따라 컨트롤 박스(120)에서 발생될 수 있다. 그러면, 테스트 구조물(10) 역시 상기 힘에 의해 움직이게 된다. 이와 같이, 상기 시스템(100)은 인위적인 지진발생 상황을 연출할 수 있다.The
테스트 구조물(10)의 상부에는 TMD(20, tuned mass damper)가 배치될 수 있는데, 도 2를 참고하여 보다 상세하게 설명하기로 한다.A tuned mass damper (TMD) 20 may be disposed on top of the
여기서, 테스트 구조물(10) 상에는 수동 진동 제어기(passive vibration controller)인 TMD(20, tuned mass damper, 동조질량댐퍼)가 배치될 수 있다. TMD(20)는 구조물의 진동을 감쇄하기 위해서 사용될 수 있다. TMD(20)는 다소 무거운 질량의 무게추(M)를 포함하고, 상기 무게추(M)는 벽(21)에 고정된 용수철 등(22, 24)과 연결될 수 있다. 상기 무게추(M)는 1.80, 2.39 Kg가 이용될 수 있다. 상기 무게추(M)의 무게는 테스트 구조물(10)의 구조에 따라 다르게 변경될 수 있다. 도 2에 도시된 TMD(20)는 예시적으로 도시되었으므로 다양한 형태, 모양으로 구성될 수 있다.Here, a passive vibration controller (TMD) 20, a tuned mass damper, may be disposed on the
이하에서는 도 3의 블록도를 참고하여 지진 관측에 대한 예측 시스템(100)의 구성을 설명하기로 한다.Hereinafter, the configuration of the
상기 시스템(100)은 진동 테이블(110), 진동 컨트롤러(120), 가속도 센서(130), 저장부(140) 및 제어부(150)를 포함할 수 있다. 도 3에 도시된 구성요소들은 시스템(100)을 구현하는데 있어서 필수적인 것은 아니어서, 본 명세서 상에서 설명되는 시스템(100)은 위에서 열거된 구성요소들 보다 많거나, 또는 적은 구성요소들을 가질 수 있다The
진동 테이블(110)은 진동 컨트롤러(120)에 의해 제어될 수 있다. 진동 테이블(110)은 진동 컨트롤러(120)로부터 제공되는 필요전류에 따라 좌우로 진동될 수 있다. 이에 따라 진동 테이블(110) 위의 테스트 구조물(10)도 함께 진동하게 된다.The vibration table 110 may be controlled by the
진동 컨트롤러(120)는 진동 테이블(110)로부터 위치에 대해 인코딩된 피드백 신호를 수신하여 이를 제어부(150)로 전달한다. 큰 의미로 제어부(150)는 진동 컨트롤러(120)를 포함하는 제어 모듈인 바, 이하에서는 제어부(150)가 진동 컨트롤러(120)를 포함하는 것으로 설명하기로 한다.The
가속도 센서(130)는 이동하는 물체의 가속도나 출력의 세기를 측정할 수 있다. 가속도 센서(130)는 테스트 구조물(10) 각 층마다 또는 최상층에 배치되어 테스트 구조물(10)의 지진(진동) 신호를 센싱할 수 있다.The
저장부(140)는 제어부(150)의 제어에 따라 다양한 정보를 저장할 수 있다. 저장부(140)는 시스템(100)에서 구동되는 다수의 응용 프로그램, 시스템(100)의 동작을 위한 데이터들, 명령어들을 저장할 수 있다. 이러한 응용 프로그램 중 적어도 일부는, 무선 통신을 통해 외부 서버로부터 다운로드 될 수 있다. The
제어부(150)는 상기 시스템(100)을 전반적으로 제어하는 구성요소이다. 제어부(150)는 테스트하고자하는 진동신호를 생성하여 진동 컨트롤러(120)를 통해 진동 테이블(110)에 전달할 수 있다. 제어부(150)는 가속도 센서(130)로부터 수신되는 측정값을 무선 또는 유선으로 수집할 수 있다.The
이하에서는 본 발명의 실시 예에 따른 지진 관측에 대한 예측을 수행하는 내용을 자세히 설명하기로 한다.Hereinafter, a description will be made in detail of the prediction of the seismic observation according to the embodiment of the present invention.
상기 시스템(100)은 테스트 구조물(10)의 거동을 모델링할 수 있다. 여기서, 거동이란 테스트 구조물(10)이 지진(진동) 신호에 따라 움직이는 것을 말한다. 상기 시스템(100)은 테스트 구조물(10)에 발생된 진동 신호를 가속도 센서를 통해 측정하고, 아울러 측정된 진동 신호를 이용하여 수학적으로 상기 측정된 진동 신호가 얼마나 정확한지 예측하는 예측 모델을 제공할 수 있다.The
즉, 시스템(100)은 테스트 구조물(10)의 거동을 모델링하고, 실제의 테스트 구조물(10)의 거동과 수학적으로 예측한 테스트 구조물(10)의 예측 거동을 도출할 수 있다. 나아가, 상기 시스템(100)은 실제 구조물을 테스트 구조물(10)을 통해 거동을 예측할 수 있다. That is, the
이런 경우, 매번 측정을 통해 테스트 구조물(10)의 거동을 평가하지 않고, 수학적으로 테스트 구조물(10)의 거동을 예측하여 평가하면서 실제 건축물의 거동을 에측할 수 있다.In this case, the behavior of the actual structure can be observed while estimating the behavior of the
2016년 10월 18일에 공개된 “Identification of the Response of a controlled building structure subjected to seismic load by using nonlinear system models”를 참고하면, 2개의 시스템 예측/식별(system predictino/identification) 모델이 이용될 수 있다. 상기 모델은 nonlinear autoregresive with exogenous inputs(NARX)와 adaptive neuro-fuzzy inference system(ANFIS)을 포함한다. 상기 논문을 참고하면 NARX 모델이 테스트 구조물(10)의 거동에 대해 더 우수한 식별성을 갖는 바 본 명세서는 NARX 모델을 이용하여 시스템 거동을 예측하고자 하나, ANFIS 모델이 이용될 수도 있다.Two system predictor / identification models can be used, referring to "Identification of the Response to a Controlled Building Structure Subjected to Seismic Load by Using Nonlinear System Models", published Oct. 18, 2016. have. The model includes nonlinear autoregresive with exogenous inputs (NARX) and an adaptive neuro-fuzzy inference system (ANFIS). With reference to the above paper, since the NARX model has better discrimination against the behavior of the
NARX 모델의 패러미터는 입력값 및 출력값이 포함될 수 있다. 도 4(a) 내지 도 4(c)를 참고하여 설명하기로 한다. S(t)는 시간에 따른 지진 신호(제어부(150)에 의해 트리거됨)를 나타내고, a(t)는 가속도 센서에서 측정하는 측정값이고, a(t-1)은 1회 시간 지연된 측정값, a(t-2)는 2회 시간 지연된 측정값을 뜻한다. 여기서 1회 지연된 측정값이란 가속도 센서에서 측정되는 측정값 중에서 첫번째 지진신호가 지연된 측정값이고, 2회 지연된 측정값이란 첫번째 및 두번째 지진신호가 지연된 측정값을 말한다.The parameters of the NARX model can include input values and output values. Will be described with reference to Figs. 4 (a) to 4 (c). A (t) is a measurement value measured by an acceleration sensor, and a (t-1) is a measurement value of a time delayed one time , and a (t-2) denotes the measurement value delayed twice. Here, the measurement value delayed one time is a measurement value in which the first earthquake signal is delayed from the measurement values measured by the acceleration sensor, and the measurement value delayed twice is the measurement value in which the first and second earthquake signals are delayed.
도 4(a)에 따르면, 상기 모델은 SISO(single input single output)로 입력 패러미터가 하나이고, 출력 패러미터가 하나이다. 도 4(b) 및 도 4(c)는 MISO(multiple input single output)로, 입력 패러미터가 복수이고, 출력 패러미터가 하나이다.According to Fig. 4 (a), the model is a single input single output (SISO) with one input parameter and one output parameter. 4 (b) and 4 (c) show MISO (multiple input single output), which has a plurality of input parameters and one output parameter.
NARX의 예측 모델은 아래의 식 1으로 구조물의 예측 거동으로 이용할 수 있다.The predictive model of NARX can be used as the predictive behavior of the structure as shown in
<식 1><
여기서, f는 비선형함수(nonlinear function), u(t) 및 y(t)는 입출력 시퀀스, nu 및 ny는 입출력 측정의 최대 개수, nk는 입력 딜레이를 나타낸다. y(t)는 신경망 알고리즘을 적용한 비선형 함수에 해당된다. NARX 예측 모델은 입출력 패러미터들을 신경망 알고리즘을 통해 다양하게 추출하여 NARX 예측 모델에 적용할 수 있다.Where f is the nonlinear function, u (t) and y (t) are the input and output sequences, nu and ny are the maximum number of input and output measurements, and nk is the input delay. y (t) corresponds to a nonlinear function applying the neural network algorithm. The NARX prediction model can be applied to the NARX prediction model by variously extracting the input and output parameters through the neural network algorithm.
상기 논문에 개시된 바와 같이, sigmoid neurons과 wavelet neurons 가 적용될 수 있으나, 평가 결과가 우수한 sigmoid neurons 을 이용하여 nonlinear model 디자인을 평가하기 위해 sigmoid 함수를 이용될 수 있다.As described in the paper, sigmoid neurons and wavelet neurons can be applied, but the sigmoid function can be used to evaluate nonlinear model designs using sigmoid neurons with good evaluation results.
도 5는 NARX 모델의 평가에 사용되는 통계적 패러미터들을 나타낸다.Figure 5 shows the statistical parameters used in the evaluation of the NARX model.
상기 논문에 나타난 것처럼 모델 넘버는 1 내지 9일 수 있고, 모델은 모두 NARX가 선택될 수 있다. NARX(1)은 SISO, NARX(2)는 MISO이면서 입력을 S(t) 및 a(t-1)인 경우를, NARX(3)은 MISO 이면서 입력을 S(t), a(t-1), a(t-2)인 경우를 상정한 것이다.As shown in the article, the model number can be from 1 to 9, and all of the models can be selected as NARX. NARX (3) is MISO and input is S (t), a (t-1), NARX (1) is SISO and NARX (2) is MISO and input is S ) and a (t-2), respectively.
상기 NARX 알고리즘은 모델 패러미터들을 식별하고, 조정하는데 사용될 수 있다. 모델 패러미터를 식별하는 단계는 i) 패러미터 예측(estimation) 및 ii) 구조 탐지(structure detection)을 포함하고, 상기 단계는 모델 오더 선택 및 선택된 모델 오더를 적용하는 것으로 나뉠 수 있다. 모델 오더 선택은 구조 탐지의 일부로 판단될 수 있다.모델 오더를 결정하는 것은 trail and order method 에 의해 결정될 수 있다.The NARX algorithm can be used to identify and adjust model parameters. Identifying the model parameters includes i) parameter estimation and ii) structure detection, which can be divided into model order selection and applying a selected model order. Model order selection can be judged as part of structural detection. Determining the model order can be determined by the trail and order method.
모델 오더는 아래 식 2로 정해질 수 있다.The model order can be defined by
<식 2><
상술환 nu 및 ny는 입출력 측정의 최대 개수, nk는 입력 딜레이를 나타낸다. 특별히 모델별로 측정값과의 상관도를 판단하는 방법은 R2, MAE(mean absolute error), RMSE(root mean squre error)를 계산하는 것이다. 아래 식 3과 같이 R2(correlation coefficient)는 관측값 및 예측값 사이의 선행 의존 정보를 제공한다.The terms nu and ny represent the maximum number of input and output measurements, and nk represents the input delay. In particular, the method of determining the degree of correlation with measured values by model is to calculate R 2 , mean absolute error (MAE), and root mean square error (RMSE). As shown in
<식 3><
여기서, y, ys는 관측된 값이고, n은 관측 횟수이고, 타임 t =1 에서 t=n까지 예측된 성능이 나타난다. 더블도트를 포함하는 y, ys는 관측값의 평균이다.Where y and y s are observed values, n is the number of observations, and the performance predicted from time t = 1 to t = n appears. Y and y s , including double dots, are the average of the observations.
MAE(mean absolute error, 아래 식 4)는 에러의 평균 크기(magnitude)를 나타내는 지표로, 실제 측정값 및 산출된 예측값 간의 차이의 평균 크기를 나타낸다.The MAE (mean absolute error,
<식 4><
RMSE(root mean squre error, 아래 식 5)는 큰 에러에 큰 가중치를 두는 평균 크기(magnitude)를 나타내는 지표로, 실제 측정값 및 산출된 예측값 간의 차이에 따라 가중치를 다르게 부여하는 평균크기를 나타낸다.RMSE (root mean square error, Equation 5) represents an average magnitude that gives a large weight to a large error. The RMSE represents an average size giving different weights according to the difference between the actual measured value and the calculated predicted value.
<식 5>≪ EMI ID =
이때, R2의 값이 높을수록, MAE 및 RMSE 의 값이 낮을수록 측정값과 유사한 예측값이 산출된 것을 나타낸다.At this time, the higher the value of R2, the lower the values of MAE and RMSE, the more predicted values similar to the measured values are calculated.
도 5와 같이, model 8인 NARX(3)이 예측값이 측정값에 근사값을 형성하고 있다. 이와 같이 테스트 구조물(10)에 지진 신호를 부여하면서 지진 관측의 정확도가 예측될 수 있다. 이를 보여주는 그래프가 도 6이다.As shown in Fig. 5, the predicted value of NARX (3) of
도 6을 참고하면, X축이 측정되는 진동신호의 측정값을 나타내고, Y축이 예측모델(NARX 모델)을 이용한 예측값을 나타낸다. 도 6과 같이 NARX 예측값이 테스트 구조물(10)의 거동을 정확하게 예측하는데 이용될 수 있다.Referring to FIG. 6, the X-axis represents the measured value of the vibration signal to be measured, and the Y-axis represents the predicted value using the predictive model (NARX model). The NARX predicted value can be used to accurately predict the behavior of the
상기의 알고리즘을 이용하면 건물의 측정값과 예측값을 가장 효과적으로 도출하여 테스트 구조물에 계속적으로 진동값을 올려 테스트 구조물이 무너질 때까지 시험을 하지 않더라도 테스트 구조물의 다양한 패러미터들을 예측하여 측정할 수 있다.By using the above algorithm, it is possible to predict the various parameters of the test structure even without testing until the test structure is collapsed by deriving the measured value and the predicted value of the building most effectively and continuously increasing the vibration value to the test structure.
도 7은 실시 예에 따른 예측 시스템의 구동 방법을 나타낸다.7 shows a method of driving the prediction system according to the embodiment.
일단, 시뮬레이션 단계(S710내지S750)와 테스트 단계(S760 내지 S780)로 나뉠 수 있다. 시뮬레이션 단계는 보다 낮은 질량의 M을 구비한 TMD가 사용되는 단계이고, 테스트 단계는 보다 높은 질량의 M을 구비한 TMD가 사용되는 단계를 말한다. 이렇게 단계적으로 무게를 높여가며 테스트 구조물(10) 및 나아가 실제 건축물의 거동이 평가될 수 있다.The simulation steps S710 to S750 and the test steps S760 to S780 can be performed. The simulation step is a step in which a TMD with a lower mass M is used and a test step is a step in which a TMD with a higher mass M is used. The behavior of the
시뮬레이션 단계를 먼저 설명하기로 한다. The simulation step will be described first.
TMD 질량이 1.8 Kg에 대한 지진 신호인 입력과 가속도 센서 측정값인 출력을 선택한다(S710).An input, which is an earthquake signal for 1.8 Kg mass and an acceleration sensor measurement value, is selected (S710).
상기 1.8Kg 은 테스트 구조물(10)의 제원에 따라 달리 설정될 수 있다.The 1.8 Kg may be set differently depending on the specifications of the
그 다음으로, 상기 시스템(100)은 신호의 노이즈를 제어한다(S720). Next, the
이에 따라, 관측된 측정값의 노이즈가 제거될 수 있다.Thereby, the noise of the observed measurement value can be removed.
그 후에, 모델의 패러미터 및 오더를 입력하고(S730), 모델의 패러미터의 조정을 수행한다(S740).Thereafter, parameters and orders of the model are input (S730), and the parameters of the model are adjusted (S740).
본 단계는 복수의 예측 모델에 대해 복수의 패러미터에 따라 연산을 수행하고, 가장 최적의 패러미터가 무엇인지, 가장 최적의 예측 모델이 무엇인지 결정할 수 있다.In this step, calculations are performed on a plurality of prediction models according to a plurality of parameters, and it is possible to determine what is the most optimal parameter and which is the most optimal prediction model.
만약 개선이 있는 경우(S750), 시뮬레이션 단계를 마치고 테스트 단계로 진행한다(S760S 내지 S780).If there is an improvement (S750), the simulation step is completed and the process proceeds to the test step (S760S to S780).
테스트 단계에서, TMD 질량을 2.3 Kg에 대한 구조물 출력을 시험한다(S760). In the test phase, the structure output for 2.3 Kg of TMD mass is tested (S760).
개선이 있는 경우(S770), 최적 모델로 선택한다(S780).If there is an improvement (S770), the optimum model is selected (S780).
개선을 판단하는 것은 상술한 바와 같이, R2, MAE, RMSE가 이용될 수 있으나, 실시 예는 이에 국한되지 않는다.R 2 , MAE, and RMSE may be used to determine the improvement, but the embodiments are not limited thereto.
전술한 본 발명은, 프로그램이 기록된 매체에 컴퓨터가 읽을 수 있는 코드로서 구현하는 것이 가능하다. 컴퓨터가 읽을 수 있는 매체는, 컴퓨터 시스템에 의하여 읽혀질 수 있는 데이터가 저장되는 모든 종류의 기록장치를 포함한다. 컴퓨터가 읽을 수 있는 매체의 예로는, HDD(Hard Disk Drive), SSD(Solid State Disk), SDD(Silicon Disk Drive), ROM, RAM, CD-ROM, 자기 테이프, 플로피 디스크, 광 데이터 저장 장치 등이 있으며, 또한 캐리어 웨이브(예를 들어, 인터넷을 통한 전송)의 형태로 구현되는 것도 포함한다. 또한, 상기 컴퓨터는 시스템(100)의 제어부(150)를 포함할 수도 있다. 따라서, 상기의 상세한 설명은 모든 면에서 제한적으로 해석되어서는 아니되고 예시적인 것으로 고려되어야 한다. 본 발명의 범위는 첨부된 청구항의 합리적 해석에 의해 결정되어야 하고, 본 발명의 등가적 범위 내에서의 모든 변경은 본 발명의 범위에 포함된다.The present invention described above can be embodied as computer-readable codes on a medium on which a program is recorded. The computer readable medium includes all kinds of recording devices in which data that can be read by a computer system is stored. Examples of the computer readable medium include a hard disk drive (HDD), a solid state disk (SSD), a silicon disk drive (SDD), a ROM, a RAM, a CD-ROM, a magnetic tape, a floppy disk, , And may also be implemented in the form of a carrier wave (e.g., transmission over the Internet). The computer may also include a
Claims (6)
복수의 플로어(floor)를 포함하고, 상기 복수의 플로어 중 최상층에 가속도센서가 구비된 테스트 구조물;
상기 테스트 구조물에 발생되는 진동을 감쇄시키는 TMD(tuned mass damper);
상기 테스트 구조물을 지지하면서 상기 테스트 구조물에 지진 신호를 트리거하는 진동 테이블; 및
상기 진동 테이블을 제어하는 제어부를 포함하며,
상기 제어부는,
상기 가속도센서를 통해 측정되는 실제 측정값 및 적어도 하나의 예측 모델을 이용하여 상기 실제 측정값을 예측하는 예측값을 산출하고, 소정의 조건을 만족하는 예측 모델을 선택하되,
상기 소정의 조건은,
상기 실제 측정값 및 산출된 예측값 간 차이의 평균 크기를 나타내는 평균 절대치 에러(mean absolute error, MAE)의 값, 또는 상기 실제 측정값 및 산출된 예측값 간 차이에 따라 가중치를 다르게 부여하는 평균 제곱근 에러(root mean square error, RMSE) 값이 최소인 경우에 만족되며,
상기 소정의 조건을 만족하는 예측 모델은,
상기 지진 신호를 나타내는 S(t), 상기 실제 측정값의 제1 연기(delay) 측정값 및 상기 실제 측정값의 제2 연기 측정값 중 적어도 하나를 입력값으로, 상기 가속도센서의 예측되는 거동(behavior)값을 출력값으로 적용하고,
상기 소정의 조건을 만족하는 예측 모델은,
신경망 알고리즘을 적용하는 NARX(nonlinear autoregresive with exogenous inputs) 함수를 사용하며,
상기 TMD에 포함된 무게추의 무게는 시뮬레이션 단계 및 테스트 단계에서 서로 다르게 적용되며, 상기 시뮬레이션 단계에서 최적의 예측 모델이 결정되면, 상기 테스트 단계에서 상기 무게추의 무게가 증가되는, 지진 관측에 대한 예측 시스템.A prediction system for seismic observations,
A test structure including a plurality of floors and having an acceleration sensor on the top of the plurality of floors;
A tuned mass damper (TMD) for attenuating vibration generated in the test structure;
A vibration table for supporting the test structure and triggering an earthquake signal on the test structure; And
And a control unit for controlling the vibration table,
Wherein,
Calculating a predictive value for predicting the actual measurement value using the actual measurement value measured through the acceleration sensor and at least one predictive model, selecting a predictive model satisfying a predetermined condition,
The predetermined condition is that,
A mean square error (MAE) value indicating a mean magnitude of the difference between the actual measured value and the calculated predicted value or a mean square error (MAE) value giving a different weight according to the difference between the actual measured value and the calculated predicted value root mean square error (RMSE) is the minimum,
The predictive model that satisfies the predetermined condition,
At least one of a first delay measurement value of the actual measurement value and a second delay measurement value of the actual measurement value as an input value, S (t) representing the seismic signal, a predicted behavior of the acceleration sensor behavior) value as an output value,
The predictive model that satisfies the predetermined condition,
We use nonlinear autoregressive with exogenous inputs (NARX)
The weight of the weight included in the TMD is applied differently in the simulation step and the test step, and when the optimal prediction model is determined in the simulation step, the weight of the weight is increased in the test step. Prediction system.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020170083683A KR101901352B1 (en) | 2017-06-30 | 2017-06-30 | Predicting System For Seismic Observation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020170083683A KR101901352B1 (en) | 2017-06-30 | 2017-06-30 | Predicting System For Seismic Observation |
Publications (1)
Publication Number | Publication Date |
---|---|
KR101901352B1 true KR101901352B1 (en) | 2018-09-21 |
Family
ID=63720816
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020170083683A KR101901352B1 (en) | 2017-06-30 | 2017-06-30 | Predicting System For Seismic Observation |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101901352B1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257934A (en) * | 2020-01-17 | 2020-06-09 | 哈尔滨工业大学 | Seismic oscillation peak acceleration prediction method based on second-order neuron deep neural network |
KR20220011334A (en) * | 2020-07-21 | 2022-01-28 | 한양대학교 산학협력단 | System and method for the damage location detection of composite structures based on a convolutional neural network |
CN115099523A (en) * | 2022-07-20 | 2022-09-23 | 哈尔滨工业大学 | Earthquake motion parameter prediction method based on artificial intelligence and earthquake motion empirical model mixing |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150241300A1 (en) | 2014-02-25 | 2015-08-27 | Industry Foundation Of Konyang University | Structural State Evaluation System for Smart Structure Using Multi-Sensing |
KR101598417B1 (en) | 2014-12-04 | 2016-03-02 | 순천향대학교 산학협력단 | System identification method for Nonlinear behavior of smart construction structure |
-
2017
- 2017-06-30 KR KR1020170083683A patent/KR101901352B1/en active IP Right Grant
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150241300A1 (en) | 2014-02-25 | 2015-08-27 | Industry Foundation Of Konyang University | Structural State Evaluation System for Smart Structure Using Multi-Sensing |
KR101598417B1 (en) | 2014-12-04 | 2016-03-02 | 순천향대학교 산학협력단 | System identification method for Nonlinear behavior of smart construction structure |
Non-Patent Citations (1)
Title |
---|
김규진, 오병관, 박효선. 건축 구조물의 지진 응답 예측 모델 제안. 2016년 대한건축학회 추계학술발표대회논문집. 제36권 제2호(통권 제66집), 2016년 10월 4-6일. 744-745쪽 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257934A (en) * | 2020-01-17 | 2020-06-09 | 哈尔滨工业大学 | Seismic oscillation peak acceleration prediction method based on second-order neuron deep neural network |
CN111257934B (en) * | 2020-01-17 | 2022-03-11 | 哈尔滨工业大学 | Seismic oscillation peak acceleration prediction method based on second-order neuron deep neural network |
KR20220011334A (en) * | 2020-07-21 | 2022-01-28 | 한양대학교 산학협력단 | System and method for the damage location detection of composite structures based on a convolutional neural network |
KR102391064B1 (en) * | 2020-07-21 | 2022-04-27 | 한양대학교 산학협력단 | System and method for the damage location detection of composite structures based on a convolutional neural network |
CN115099523A (en) * | 2022-07-20 | 2022-09-23 | 哈尔滨工业大学 | Earthquake motion parameter prediction method based on artificial intelligence and earthquake motion empirical model mixing |
CN115099523B (en) * | 2022-07-20 | 2023-10-24 | 哈尔滨工业大学 | Earthquake motion parameter prediction method based on artificial intelligence and earthquake motion experience model mixing |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101901352B1 (en) | Predicting System For Seismic Observation | |
WO2014109324A1 (en) | Building safety verification system and building safety verification method | |
CA3018790C (en) | Seismic response assessment of man-made structures | |
Lin et al. | Semiactive control of building structures with semiactive tuned mass damper | |
Karami et al. | Developing a smart structure using integrated subspace‐based damage detection and semi‐active control | |
WO2017056426A1 (en) | Soil quality determination device, soil quality determination method, and recording medium having program stored thereon | |
JP6614165B2 (en) | Threshold decision method, threshold decision device, and elevator control system | |
JP6001740B1 (en) | High precision evaluation of structure transfer functions, earthquake response prediction, deterioration diagnosis system and method | |
EP2733103B1 (en) | Elevator operation control method and operation control device | |
KR101891175B1 (en) | Safety evaluation method for structure using modified dynamaic characteristics | |
Cara | Computing the modal mass from the state space model in combined experimental–operational modal analysis | |
KR102209749B1 (en) | Method and system for generating earthquake acceleration time history | |
JP5281475B2 (en) | Building health diagnostic method, diagnostic device and diagnostic program based on microtremor measurement | |
JP2015161657A (en) | Damage evaluation system and damage evaluation method for building | |
JP5495684B2 (en) | Elevator earthquake damage prediction device | |
Chu et al. | Application of the on‐line recursive least‐squares method to perform structural damage assessment | |
Xie et al. | Innovative substructure approach to estimating structural parameters of shear structures | |
Hashemi et al. | Vibration base identification of impact force using genetic algorithm | |
JP5799183B2 (en) | Building safety verification system, building safety verification method and program | |
JP2020004052A (en) | Earthquake analysis device, earthquake analysis system, and earthquake analysis method | |
JP2016017848A (en) | Structure verification system, structure verification device, and structure verification program | |
JP6499832B2 (en) | Structure safety verification system, structure safety verification method and program | |
JP5760268B2 (en) | Earthquake motion duration prediction system | |
JP2014102231A (en) | Safety evaluation device in seismic isolation architectural structure | |
Beskhyroun et al. | Ambient and forced vibration testing of a 13-story reinforced concrete building |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant |