CN104849750B - Nuclear power floor based on the analysis of target waveform composes artificial ripple approximating method - Google Patents
Nuclear power floor based on the analysis of target waveform composes artificial ripple approximating method Download PDFInfo
- Publication number
- CN104849750B CN104849750B CN201510191268.2A CN201510191268A CN104849750B CN 104849750 B CN104849750 B CN 104849750B CN 201510191268 A CN201510191268 A CN 201510191268A CN 104849750 B CN104849750 B CN 104849750B
- Authority
- CN
- China
- Prior art keywords
- ripple
- target
- spectrum
- artificial
- floor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention belongs to the safety analysis of nuclear power generating equipment earthquake resistant engineering and assessment technique field, it is related to a kind of nuclear power floor based on the analysis of target waveform to compose Artificial Seismic Wave approximating method.It is characterized in that first analyzing selected floor spectrum target wave amplitude, phase angle, frequency and energy feature；Target spectrum is obtained by uncertain artificial treatment to target ripple response spectrum；Calculate and correspond to target spectrum Fourier amplitude spectrum；Linear combination target ripple and the Fourier amplitude spectrum from target spectrum try to achieve original manual ripple Fourier amplitude spectrum；Original manual ripple is defined with target wave phase spectrum and strength envelope；Multigroup artificial ripple is generated by control condition of the strength envelope of target ripple during original manual ripple iterative fitting, artificial ripple is composed by the preferred floor of index of itself linearly dependent coefficient and between target wave energy curve between time-history curves.The present invention can make floor compose artificial ripple has the wave character of certain destination time-histories reaction, can improve expression effect and engineering adaptability that floor composes artificial ripple.
Description
Technical field
Artificial Seismic Wave approximating method is composed the present invention relates to a kind of nuclear power floor based on the analysis of target waveform, belongs to core
The safety analysis of electric equipment earthquake resistant engineering and assessment technique field.
Background technology
Ensure that Nuclear Safety is to ensure that nuclear power plant construction smoothly implements the key with safe operation, ensure that nuclear power resists among these
The key for shaking safety is that Nuclear Power Station safe level engineering structure and equipment will meet structure complete under OBE and SSE grades of seismic condition
The basic demand of whole property and operational feasibility.Compared to the limit seismic bearing capacity of nuclear power engineering structure, safe level nuclear power generating equipment
One of control condition for often more paying close attention to can be normally run, and the aseismic analysis of nuclear power generating equipment are with nuclear power work where it
The floor earthquake response spectrum of journey structure is used as input foundation.
The principle of conservative design is followed, the aseismic analysis of nuclear power engineering structure often use various models, a plurality of input ground
Seismic wave carries out comparing calculation, and floor spectrum reaction result to obtaining takes largest enveloping.Meanwhile, consider to calculate ginseng by code requirement
The probabilistic influence of number, carries out expanding peak value frequency range and reduces the artificial treatment of peak height, finally to floor spectrum calculated value
Formation can be used for the floor response spectrum that nuclear power generating equipment antidetonation is checked.It can be seen that, eventually for the destination spectrum of nuclear power generating equipment antidetonation
It is corresponding to the statistics of nuclear power floor dynamic response under a certain security level.On this basis, intend by using artificial ripple
Conjunction technology forms floor and composes artificial ripple, for the time domain kinematic analysis of safe level nuclear power generating equipment.The building structure that compares input
Earth's surface design response spectrum in earthquake engineering, floor response spectrum contains the frequency characteristic of more structures, so as to often have
There is complex-shaped obvious characteristic, be mainly manifested in the aspects such as peak value is more, peak valley difference is big, give the fitting of artificial wave height precision, and
Waveform control brings not molehill.
From intuitively angle, engineer and nuclear facility personnel, in the design of equipment quake-resistant safety and evaluating
The time-histories feature that the artificial ripple of input for being used should as much as possible have place floor to compose true dynamic response, this is also the present invention
The starting point of consideration.Determine rational technical indicator and Technology Ways, the true dynamic response of structural floor is special where description
Levy, and these information are inherited as much as possible in the floor with statistical nature composes artificial ripple fitting, acquisition meets target ripple ripple
The time generated seismic wave of shape.
At present, the research in terms of artificial ripple fitting, is concentrated mainly on response spectrum iterative fitting method.Usual way is to set up
Mathematical Modeling regards seismic wave as the superposition of the trigonometrical number that frequency is different, phase angle is random, is pushed away using Maharaj K.Kaul
Derived response spectrum and power spectrum approximate relation obtain corresponding power spectrum, then go to synthesize artificial ripple by this power spectrum.Usual people
Work ripple fitting core procedure be：
(1) composed from aimed acceleration response spectrum corrected power, amplitude spectrum is obtained by power spectrum.
(2) original manual seismic wave is generated based on random phase angle.
(3) by the method such as peak regulation value and Overlapping intensities envelope, the shape of the time-history curves of generation is controlled.
(4) the iterative fitting process of goal response spectrum.
However, above-mentioned artificial ripple fit procedure is to obtain power spectrum i.e. averagely from the approximation relation of goal response spectrum
Fourier spectra, while using random phase angle, the wave character of target ripple is not embodied；Strength envelope in iterative process
General is also statistics, it is believed that factor is relevant with earthquake magnitude size, the distance of epicentral distance etc. for it, not for a certain target ripple.
The content of the invention
The present invention is a kind of holistic approach of the nuclear power floor spectrum Artificial Seismic Wave fitting based on the analysis of target waveform.For
Floor is set to compose intensity envelope, non stationary frequencies, the wave character of Energy distribution that artificial ripple has target ripple simultaneously, it is proposed that
Based on target waveform specificity analysis, the floor that can inherit target ripple basic waveform feature composes artificial ripple approximating method.
Innovation of the invention has：The Fourier amplitude spectrum of original manual ripple is based on the Fourier width that optimization takes target ripple
The linear combination of value spectrum and the Fourier amplitude spectrum from destination spectrum；Using target ripple in response spectrum iterative process
Phase angle and strength envelope, i.e., using the essential characteristic of target ripple；Based on the linear phase between energy curve between time-history curves
Closing property investigates the feature association that floor composes artificial ripple and target ripple.
The present invention comprises the following steps：
(1) it is target ripple from true floor response time-histories, analyzes its amplitude, phase angle, frequency, energy feature, it is determined that
The absolute value strength envelope of target ripple；By carrying out the Fourier amplitude spectrum F that Fourier transform obtains target ripple to target ripple1
(w) and phase spectrum θ1(w), energy distribution curve
E (t)=∫ β2(t)dt。
(2) response spectrum to target ripple obtains target spectrum by uncertain artificial treatment, calculates what target was composed by formula
Approximate power spectral, obtains Fourier amplitude spectrum F of the original manual ripple corresponding to target spectrum2(w)。
(3) by formula F (w)=α F1(w)+(1-α)F2W () tries to achieve the final Fourier amplitude that original manual ripple is used
Spectrum, 0≤α≤1, it is 0.01 that α can take step-length, generates Fourier's amplitude spectrum of multigroup alternative original manual ripple.
(4) phase angle and shape facility of original manual ripple are defined using the phase spectrum and strength envelope of target ripple, i.e.,
Using the essential characteristic of target ripple, several original manual ripples are generated.
(5) each bar original manual ripple obtained for step (4) is generated corresponding many based on response spectrum iterative fitting process
Bar time generated seismic wave, the strength envelope during iterative fitting with target ripple is as control condition.
(6) linearly dependent coefficient between the artificial ripple and target wave energy curve that are obtained with step (5) between time-history curves
It is Artificial Seismic Wave to be asked to recommend coefficient correlation artificial ripple high in the set of the artificial ripple obtained from step (5) for index.Its
Whether the energy processes that the coefficient correlation between middle energy curve reflects these artificial ripples are similar to target ripple；Time-history curves
Between coefficient correlation then more directly reflect these artificial ripples and relevance of the target ripple on waveform.
The present invention can make floor compose artificial ripple has the feature of certain destination time-histories reaction, and solve is carried out at present
Nuclear power generating equipment quake-resistant safety is designed composes artificial ripple with building structure floor response time history in frequency with the floor used when evaluating
The larger problem of difference on non-stationary, Energy distribution, wave character, can improve expression effect and engineering that floor composes artificial ripple
Applicability.
Brief description of the drawings
Fig. 1 is the nuclear power floor spectrum Artificial Seismic Wave approximating method flow chart based on the analysis of target waveform.
Fig. 2 a are that the floor obtained by certain transactional analysis in embodiment as target ripple responds time-histories.
Fig. 2 b be in embodiment as target ripple the floor response time-histories absolute value obtained by certain transactional analysis and
Its absolute value strength envelope.
Fig. 2 c are the bandwidth amplitude spectrums of Fourier 20% that embodiment floor composes target ripple.
Fig. 3 is the destination spectrum after the response spectrum and uncertain artificial treatment of embodiment target ripple.
Fig. 4 is the time generated seismic wave sample that generally artificial ripple fit procedure is obtained.
Fig. 5 a are that the floor synthesized by embodiment method composes time generated seismic wave.
Fig. 5 b be the floor synthesized by embodiment method composes artificial ripple, generally artificial ripple fit procedure is obtained artificial ripple with
Target wave energy distribution curve comparison diagram.
Fig. 5 c are the response spectrum of the artificial ripple synthesized by embodiment method, the artificial ripple of generally artificial ripple fit procedure acquisition
Response spectrum and target compose comparison diagram.
In figure：1 selected floor dynamic response time-histories (target ripple) curve；
2 target ripple acceleration absolute value time-histories；3 target ripple Acceleration time course absolute value strength envelopes；
The bandwidth amplitude spectrum of 4 target ripple Fourier 20%；The response spectrum of 5 target ripples；
6 targets obtained through uncertain artificial treatment are composed；
The time generated seismic wave sample that 7 generally artificial ripple fit procedures are obtained；
8 compose time generated seismic wave by the floor of embodiment method synthesis；The energy curve of 9 target ripples；
10 energy curves that artificial ripple is composed by the floor of embodiment method synthesis；
The energy curve of the artificial ripple that 11 generally artificial ripple fit procedures are obtained；
12 response spectrums that artificial ripple is composed by the floor of embodiment method synthesis；
The response spectrum of the artificial ripple that 13 generally artificial ripple fit procedures are obtained.
Specific embodiment
With reference to concrete scheme and accompanying drawing, specific embodiment of the invention is described in detail.
The present embodiment method is illustrated as a example by being fitted the destination of 5% damping ratio and composing artificial ripple.Including following step
Suddenly：
(1) target waveform specificity analysis.Time-histories is responded as target ripple, such as Fig. 2 a from true floor；Calculate its anti-
Ying Pu, such as Fig. 3；Determine envelope of the target ripple absolute value envelope as artificial ripple fit procedure, such as Fig. 2 b.
Calculate the Fourier phase spectrum θ of target ripple1(w) and amplitude spectrum F1(w), such as Fig. 2 c.
Fourier amplitude spectrum F1(w), phase spectrum θ1The definition of (w) such as formula 1：
In formula 1：A (w)-true floor responds the frequency domain presentation form of time-histories.
The real part of C (w)-A (w).
The imaginary part of S (w)-A (w).
The energy distribution curve of target ripple, such as Fig. 5 b are calculated by formula 2.
E (t)=∫ β2(t)dt (2)
In formula 2：The Acceleration time course of β (t)-true floor response.
(2) uncertain artificial treatment is carried out to the floor response spectrum that target ripple is calculated, obtains destination spectrum (mesh
Mark spectrum).Method is：2-3 obvious peak value is selected in the floor response spectrum that target ripple is calculated, to frequency at this peak value
Square upward continuation 15%, at the same peak value direction superiors drop 15%, it is determined that response spectrum be target spectrum, such as Fig. 3.
(3) composed by target with usual method and calculate time generated seismic wave.The corresponding power spectrum of target spectrum is calculated, using power
The approximation relation of spectrum and amplitude spectrum obtains Fourier amplitude spectrum F2(w), such as formula 3.
In formula 3：G(wk)-be power spectrum.
According to the phase spectrum θ for being uniformly distributed hypothesis each Fourier spectral line of generation at random phase angle2W (), is derived from base
In the artificial ripple sample that the usual method of destination spectrum is obtained, such as Fig. 4, its energy distribution curve such as Fig. 5 b.
(5) the final Fourier amplitude spectrum that original manual ripple is used is calculated.Using the Fourier amplitude spectrum F of target ripple1
The Fourier amplitude spectrum F of w artificial ripple that () and generally artificial ripple fit procedure are obtained2W the linear combination of () obtains fitting artificial
Fourier amplitude spectrum F (w) of the original manual ripple required for ripple, such as formula 4.
F (w)=α F1(w)+(1-α)F2(w) (4)
In formula 4, it is big that α reflects contribution of the Fourier amplitude spectrum of target ripple to the Fourier amplitude spectrum of original manual ripple
It is small.In the present embodiment to α in the range of 0~1 value, obtain multigroup original manual ripple Fourier amplitude spectrum.
(6) θ is composed with target wave phase1W () and target intensity of wave envelope obtain a plurality of original manual as control condition
Ripple.
(7) the artificial ripple of iterative fitting.Iterative fitting process based on response spectrum, is control with the strength envelope of target ripple
Condition, fitting is iterated to a plurality of original manual ripple obtained in the previous step, until precision meets requiring.
(8) artificial ripple is preferably recommended.The energy distribution curve of a plurality of artificial ripple obtained in the previous step is calculated by formula 2；Calculate
Linearly dependent coefficient between these artificial ripples and target wave energy distribution curve and between time-history curves, recommend wherein coefficient correlation compared with
Artificial ripple high.
Recommend the artificial ripple of α=0.25, such as Fig. 5 a, its energy distribution curve such as Fig. 5 b in the present embodiment.
Time generated seismic wave (Fig. 4) and floor spectrum target ripple (Fig. 2 a) that generally artificial ripple fit procedure is obtained are calculated by formula 5
The correlation coefficient ρ of time-histories1≈ 0.01, the floor for calculating the synthesis of embodiment method composes time generated seismic wave (Fig. 5 a) with floor spectrum target
The correlation coefficient ρ of ripple time-histories (Fig. 2 a)2≈0.4。
The response spectrum for calculating the work ripple to be asked for help that embodiment method is obtained is artificial with what generally artificial ripple fit procedure was obtained
The response spectrum of ripple is simultaneously contrasted, such as Fig. 5 c with target ripple response spectrum.
From the foregoing, it will be observed that the floor of embodiment method synthesis composes artificial ripple is meeting target spectrum security, conservative requirement
Meanwhile, compared with the artificial ripple that usual artificial ripple fit procedure is obtained, with the correlation between the energy distribution curve of target ripple and
Correlation between time-history curves is significantly increased, and is provided with the basic waveform feature that floor composes target ripple.
Claims (2)
1. a kind of nuclear power floor based on the analysis of target waveform composes Artificial Seismic Wave approximating method, it is characterised in that comprising following
Step：
(1) it is target ripple from true floor response time-histories, analyzes amplitude, phase angle, frequency and the energy feature of target ripple, really
Set the goal the absolute value strength envelope of ripple；By carrying out the Fourier amplitude spectrum that Fourier transform obtains target ripple to target ripple
F1(w) and phase spectrum θ1(w), energy distribution curve E (t)=∫ β2When () dt, wherein β (t) are the acceleration of true floor response t
Journey；
(2) response spectrum to target ripple obtains target spectrum by uncertain artificial treatment, and the approximate of target spectrum is calculated by formula
Power spectrum, obtains Fourier amplitude spectrum F of the original manual ripple corresponding to target spectrum2(w)；
(3) by formula F (w)=α F1(w)+(1-α)F2W () tries to achieve the final Fourier amplitude spectrum that original manual ripple is used, 0≤
α≤1, generates Fourier's amplitude spectrum of multigroup alternative original manual ripple, and wherein α is the Fourier amplitude spectrum of target ripple to initial people
The contribution of the Fourier amplitude spectrum of work ripple；
(4) phase angle and shape facility of original manual ripple are defined using the phase spectrum and strength envelope of target ripple, that is, is used
The essential characteristic of target ripple, generates several original manual ripples；
(5) each bar original manual ripple obtained for step (4) is based on response spectrum iterative fitting process and generates corresponding a plurality of people
Work ripple time-histories, the strength envelope during iterative fitting with target ripple is as control condition；
(6) linearly dependent coefficient between the artificial ripple and target wave energy curve that are obtained with step (5) between time-history curves is finger
It is Artificial Seismic Wave to be asked to recommend coefficient correlation artificial ripple high in the set of the artificial ripple that mark is obtained from step (5).
2. a kind of nuclear power floor based on the analysis of target waveform according to claim 1 composes Artificial Seismic Wave fitting side
Method, it is characterised in that comprise the steps of：It is 0.01 that the α takes step-length.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510191268.2A CN104849750B (en) | 2015-04-21 | 2015-04-21 | Nuclear power floor based on the analysis of target waveform composes artificial ripple approximating method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510191268.2A CN104849750B (en) | 2015-04-21 | 2015-04-21 | Nuclear power floor based on the analysis of target waveform composes artificial ripple approximating method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104849750A CN104849750A (en) | 2015-08-19 |
| CN104849750B true CN104849750B (en) | 2017-06-13 |
Family
ID=53849532
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510191268.2A Expired - Fee Related CN104849750B (en) | 2015-04-21 | 2015-04-21 | Nuclear power floor based on the analysis of target waveform composes artificial ripple approximating method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104849750B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107967240B (en) * | 2017-11-23 | 2020-12-11 | 大连理工大学 | Earthquake-resistant time-course acceleration optimization algorithm based on artificial earthquake motion generation process |
| CN109002678B (en) * | 2018-10-17 | 2020-08-18 | 中国水利水电科学研究院 | Artificial earthquake simulation method based on hydroelectric engineering hydraulic structure earthquake-resistant design |
| CN109508439B (en) * | 2018-11-15 | 2023-03-24 | 中国地震局工程力学研究所 | Time course analysis method for nuclear power plant floor response spectrum |
| CN112462424A (en) * | 2020-11-13 | 2021-03-09 | 大连理工大学 | Artificial wave generation method considering real wave phase characteristics |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101236256A (en) * | 2008-01-25 | 2008-08-06 | 北京工业大学 | Multi- damping ratio goal response spectrum compatible artificial earthquake wave synthesis method |
| CN101944140A (en) * | 2010-08-08 | 2011-01-12 | 顺特阿海珐电气有限公司 | Earthquake proof performance analysis method of dry type transformer for nuclear power |
| CN104008227A (en) * | 2014-04-18 | 2014-08-27 | 江苏大学 | Optimum design method for gap between movable part and static part of centrifugal pump |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100294384B1 (en) * | 1999-05-24 | 2001-07-12 | 최동환 | Synthesis methode of shock waveform for shock response spectrum test |
-
2015
- 2015-04-21 CN CN201510191268.2A patent/CN104849750B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101236256A (en) * | 2008-01-25 | 2008-08-06 | 北京工业大学 | Multi- damping ratio goal response spectrum compatible artificial earthquake wave synthesis method |
| CN101944140A (en) * | 2010-08-08 | 2011-01-12 | 顺特阿海珐电气有限公司 | Earthquake proof performance analysis method of dry type transformer for nuclear power |
| CN104008227A (en) * | 2014-04-18 | 2014-08-27 | 江苏大学 | Optimum design method for gap between movable part and static part of centrifugal pump |
Non-Patent Citations (4)
| Title |
|---|
| 人工地震波反应谱拟合技术的改进;翟希梅 等;《哈尔滨工业大学学报》;19951231;第27卷(第6期);第130-133页 * |
| 人工拟合地震动时程参数对场地反应的影响分析;王恒知 等;《西北地震学报》;20090630;第31卷(第2期);第126-130页 * |
| 基于相位差谱的多点激励人工波数值生成算法研究;李建波 等;《大连理工大学学报》;20090731;第49卷(第4期);第558-563页 * |
| 基于粒子群优化算法的人工地震波反应谱拟合技术;王雪妮 等;《水电能源科学》;20131231;第31卷(第12期);第98-101、142页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104849750A (en) | 2015-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104849750B (en) | Nuclear power floor based on the analysis of target waveform composes artificial ripple approximating method | |
| Mohamed | Aero-acoustics noise evaluation of H-rotor Darrieus wind turbines | |
| Yunus et al. | ARIMA-based frequency-decomposed modeling of wind speed time series | |
| CN104201700B (en) | Meter and the regional power grid thermoelectricity frequency modulation crew qiting method of the uncertain fluctuation of wind-powered electricity generation | |
| CN106709641A (en) | Monte-Carlo simulation based small interference probability risk analysis and simulation method | |
| Chen et al. | The investigation of tower height matching optimization for wind turbine positioning in the wind farm | |
| CN106873359A (en) | A kind of wind-powered electricity generation noise rating method based on cluster analysis and neutral net | |
| Gong et al. | Influence of non-Gaussian wind characteristics on wind turbine extreme response | |
| CN104283210B (en) | Based on the two-stage distributed emulation mode of wind-powered electricity generation collection region that Distributed Power Flow calculates | |
| CN101697014A (en) | Earthquake motion synthesis method and system | |
| CN107070568A (en) | A kind of frequency spectrum sensing method based on Hilbert-Huang transform | |
| CN108647415A (en) | The reliability estimation method of electric system for high proportion wind-electricity integration | |
| CN104899432A (en) | Kernel function combination-based PSO-LSSVM fluctuating wind speed prediction method | |
| CN104077478A (en) | Numerical simulation method for downburst non-stationary fluctuating wind speed | |
| CN109214119A (en) | Bridge Earthquake Resistance Design method based on response surface model | |
| CN107977492A (en) | Based on the non-linear windage yaw reliability degree calculation method of Monte Carlo insulator chain | |
| Chen et al. | Probabilistic analysis of wind-induced vibration mitigation of structures by fluid viscous dampers | |
| Chen et al. | The iteration method for tower height matching in wind farm design | |
| CN108169794B (en) | Adjusting method for matching near-fault artificial seismic waves with standard response spectrum | |
| Ke et al. | Non‐Gaussian characteristics and extreme distribution of fluctuating wind pressures on large cylindrical–conical steel cooling towers | |
| Gu | Wind-resistant studies on tall buildings and structures | |
| Viterna | The NASA-LeRC wind turbine sound prediction code | |
| CN110490215A (en) | A kind of modal identification method and system of wind power plant-power grid interaction | |
| Shuang et al. | Reliability analysis of wind turbines under non‐Gaussian wind load | |
| di Roseto et al. | Performance-based seismic design of a modular pipe-rack |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Li Jianbo Inventor after: Zhou Lei Inventor after: Mei Runyu Inventor after: Pan Dong Inventor after: Lin Gao Inventor after: Liu Jun Inventor before: Li Jianbo Inventor before: Zhou Lei Inventor before: Mei Runyu Inventor before: Lin Gao Inventor before: Liu Jun |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20170519 Address after: 116024 Liaoning, Dalian, Ganjingzi Ling Road, No. 2 Applicant after: Dalian University of Technology Applicant after: China International Water & Electric Corp. Address before: 116024 Liaoning, Dalian, Ganjingzi Ling Road, No. 2 Applicant before: Dalian University of Technology |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170613 Termination date: 20200421 |