CN109116753A - A kind of applying method of high-voltage power transmission tower space rigid finite element model wind load - Google Patents
A kind of applying method of high-voltage power transmission tower space rigid finite element model wind load Download PDFInfo
- Publication number
- CN109116753A CN109116753A CN201811047744.3A CN201811047744A CN109116753A CN 109116753 A CN109116753 A CN 109116753A CN 201811047744 A CN201811047744 A CN 201811047744A CN 109116753 A CN109116753 A CN 109116753A
- Authority
- CN
- China
- Prior art keywords
- wind
- unit
- windward side
- sample
- power transmission
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
A kind of applying method of high-voltage power transmission tower space rigid finite element model wind load, the application of wind load is divided into first accurate definition windward side component by the method, then sequencing extraction unit data and calculates the area of each windward side unit and the coordinate of unit center point；Wind speed simulation program is worked out using MATLAB, simulates Wind Velocity History sample, calculates wind pressure time-histories sample, and combines the area of windward side unit and the coordinate of unit center point, concentrated force time-histories sample is calculated, is packaged into independent interface document.The present invention works out wind speed simulation program using MATLAB, simulates Wind Velocity History sample, calculates wind pressure time-histories sample, and combine the area of windward side unit and the coordinate of unit center point, concentrated force time-histories sample is calculated, is packaged into independent interface document, it is ensured that result is accurate.
Description
Technical field
The present invention relates to a kind of applying methods of high-voltage power transmission tower space rigid finite element model wind load, and it is defeated to belong to high pressure
Ferroelectric tower technology field.
Background technique
The accurate application of space rigid Model load is the base that electric power pylon mechanical property carries out Deterministic Finite meta analysis
Plinth.And the advantage of ANSYS platform is finite element analysis, difficult to the graphic operation of large complicated shape, electric power pylon is in kind
The tower body rod piece of model composition is thousands of easily, and quantity is quite big, and sectional dimension is more, to the accurate application difficulty of its wind load compared with
Greatly.
Electric power pylon mainly has three class models, i.e. space truss, trusses mixing and space rigid model.A large amount of research and
Engineering practice shows in three class models that Space Truss Modeling and calculating are all most simple, but precision is minimum, and Space truss method is big
Engineering design requirements are fully met in most cases, but when load case complexity, calculated result and agreement with experimental data are obtained not
It reaches.Space rigid model accuracy is high, computational accuracy highest closest with actual loading situation, but its accurate modeling and outer load
Lotus, which accurately applies, to be also most difficult to.
Summary of the invention
The object of the present invention is in order to realize that high-voltage power transmission tower space rigid finite element model wind load is quick and precisely applied
Add, the present invention proposes a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load.
In order to solve the above technical problems, the technical solution adopted by the present invention is that: a kind of high-voltage power transmission tower space rigid has
The applying method for limiting meta-model wind load, is divided into first accurate definition windward side component for the application of wind load, then sequencing mentions
It takes cell data and calculates the area of each windward side unit and the coordinate of unit center point；Wind speed simulation is worked out using MATLAB
Program simulates Wind Velocity History sample, calculates wind pressure time-histories sample, and combine the area and unit center point of windward side unit
Coordinate, calculate concentrated force time-histories sample, be packaged into independent interface document.
A kind of applying method of high-voltage power transmission tower space rigid finite element model wind load, comprising the following steps:
(1) on the basis of ANSYS platform completion electric power pylon space rigid finite element model establishes work, basic step is such as
Under:
(1-1) defines windward side unit block；
(1-2) extraction unit data simultaneously calculate the area of each windward side unit and the coordinate of unit center point.
(2) in the work of MATLAB platform, basic step is as follows:
(2-1) works out wind speed simulation program, simulates each list according to the coordinate data of each unit central point in step (1-2)
First Wind Velocity History sample；
(2-2) calculates each unit wind pressure time-histories sample according to each unit Wind Velocity History sample in step (2-1)；
(2-3) is according to the area of each windward side unit in each unit wind pressure time-histories sample in step (2-2) and step (1-2)；Meter
Calculate the wind load of each windward side unit.
(3) in the work of ANSYS platform, basic step is as follows:
It is applied in step (2-1) using the wind load for meeting each wind face unit in step (2-3) as concentrated force load corresponding each
The central point of unit.
The high-voltage power transmission tower is voltage class 110kV and the above various types angle steel tower.
Step (2-1) the wind speed simulation program, using linear filtering method or harmonic and reactive detection method, wave spectrum is used
Davenport wind speed spectrum or Kaimal wind speed spectrum.
The step (2) is packaged into independent interface document in the course of work of MATLAB platform.
Described program extraction unit data, to two extreme coordinates values of each unit of windward side unit record of selection
And sectional dimension, to calculate the area of each unit center point coordinate, element length, unit.
It is described to simulate Wind Velocity History sample, it is programmed with MATLAB, using linear filtering method or harmonic and reactive detection method, wave
Spectrum uses Davenport wind speed spectrum or Kaimal wind speed spectrum analog wind speed, directly obtains a series of numbers that wind speed changes over time
According to as Wind Velocity History sample.
The wind pressure time-histories sample out can calculate Wind Velocity History sample according to the conversion formula of wind pressure and wind speed
Wind pressure time-histories sample.
The concentrated force time-histories sample sits the wind load load of each windward side unit at the center of windward side unit
It is concentrated force on punctuate, wind load changes over time, and wind pressure time-histories sample is obtained multiplied by corresponding unit area
The concentrated force time-histories sample of unit.
The wind load of the windward side unit, the wind load of each windward side unit be equal to wind pressure multiplied by it is each windward
Face cellar area.
The invention has the advantages that consider ANSYS platform difficult to the graphic operation of large complicated shape by the present invention,
And electric power pylon tower body rod piece quantity is quite big, sectional dimension is more, it is larger to the accurate application difficulty of its wind load and it is easy go out
Existing mistake, is divided into first accurate definition windward side component for the application of wind load, then sequencing extraction unit data and calculates
The area of each windward side unit and the coordinate of unit center point, reduce graphic operation difficulty, while guaranteeing the accuracy of data.
The present invention works out wind speed simulation program using MATLAB, simulates Wind Velocity History sample, calculates wind pressure time-histories sample
This, and the area of windward side unit and the coordinate of unit center point are combined, concentrated force time-histories sample is calculated, is packaged into independent
Interface document, it is ensured that result is accurate.
The present invention realizes wind load and applies sequencing, improves speed and accuracy that wind load applies, modifies.
Detailed description of the invention
Fig. 1 is step flow chart of the invention；Fig. 2 is ZMS2_39 type high-voltage power transmission tower rigid frame model；
Fig. 3 is ZMS2_39 type high-voltage power transmission tower windward side unit block；
Fig. 4 is Wind Velocity History sample；Fig. 5 is wind pressure time-histories sample；
Fig. 6 is that ZMS2_39 type high-voltage power transmission tower wind applies illustraton of model after load.
Specific embodiment
A specific embodiment of the invention is as shown in Figure 1.
A kind of applying method of high-voltage power transmission tower space rigid finite element model wind load of the present embodiment, steps are as follows:
Step 1, as shown in Fig. 2, definition is met on the basis of ANSYS platform establishes ZMS2_39 sections tower rigid frame finite element model
Wind face unit block, as shown in figure 3, extraction unit data and calculating the area of each windward side unit and the seat of unit center point
Mark.
A component is created in ANSYS platform, manually picks up windward side unit, this step operation must carefully, and be sure to
The correlation of lifting unit and geometry, Fig. 3 are the windward side unit different angle view finally defined.
How sequencing extraction unit data: to two extreme coordinates values of each unit of windward side unit record of selection
And sectional dimension, so as to calculate the area of each unit center point coordinate, element length, unit.Step 2, in MATLAB
Platform, the coordinate data of combining unit central point work out wind speed simulation program, simulate Wind Velocity History sample, be illustrated in figure 4
The Wind Velocity History sample (i.e. wind speed change with time relationship) obtained with Davenport wind speed spectrum harmonic and reactive detection method；According to wind
Fast time-histories sample calculates wind pressure time-histories sample (i.e. wind pressure change with time relationship), as shown in Figure 5.
Wind Velocity History sample is simulated, is programmed with MATLAB, is adopted using linear filtering method or harmonic and reactive detection method, wave spectrum
With Davenport wind speed spectrum or Kaimal wind speed spectrum analog wind speed, the volume of data that wind speed changes over time is directly obtained, i.e.,
For Wind Velocity History sample.
Wind Velocity History sample can be calculated outlet air according to the conversion formula of wind pressure and wind speed by the calculating of wind pressure time-histories sample
Press time-histories sample.
Step 3, in ANSYS platform, in conjunction with the windward side unit in step 1 area and step 2 in wind pressure time-histories
Sample calculates concentrated force time-histories sample, and is applied to the central point of corresponding unit, completes wind load as shown in FIG. 6 and applies；
When Fig. 6 is the 200th second after the wind load of ZMS2_39 type high-voltage power transmission tower space rigid model windward side unit application concentrated force
The schematic diagram at quarter.
The calculating of the wind load of windward side unit: the wind load of each windward side unit is equal to wind pressure and meets multiplied by each
Wind face cellar area.
Claims (10)
1. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load, which is characterized in that
The application of wind load is divided into first accurate definition windward side component by the method, then sequencing extraction unit data and is counted
Calculate the area of each windward side unit and the coordinate of unit center point；Wind speed simulation program is worked out using MATLAB, simulates outlet air
Fast time-histories sample calculates wind pressure time-histories sample, and combines the area of windward side unit and the coordinate of unit center point, calculates
Concentrated force time-histories sample, is packaged into independent interface document.
2. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load, which is characterized in that its feature exists
In the described method comprises the following steps:
(1) on the basis of ANSYS platform completion electric power pylon space rigid finite element model establishes work, basic step is such as
Under:
(1-1) defines windward side unit block；
(1-2) extraction unit data simultaneously calculate the area of each windward side unit and the coordinate of unit center point；
(2) in the work of MATLAB platform, basic step is as follows:
(2-1) works out wind speed simulation program, simulates each list according to the coordinate data of each unit central point in step (1-2)
First Wind Velocity History sample；
(2-2) calculates each unit wind pressure time-histories sample according to each unit Wind Velocity History sample in step (2-1)；
(2-3) is according to the area of each windward side unit in each unit wind pressure time-histories sample in step (2-2) and step (1-2)；Meter
Calculate the wind load of each windward side unit；
(3) in the work of ANSYS platform, basic step is as follows:
It is applied in step (2-1) using the wind load for meeting each wind face unit in step (2-3) as concentrated force load corresponding each
The central point of unit.
3. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load according to claim 1,
It is characterized in that, the high-voltage power transmission tower is voltage class 110kV and the above various types angle steel tower.
4. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load according to claim 1,
It is characterized in that, step (2-1) the wind speed simulation program, using linear filtering method or harmonic and reactive detection method, wave spectrum is used
Davenport wind speed spectrum or Kaimal wind speed spectrum.
5. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load according to claim 1,
It is characterized in that, the step (2) is packaged into independent interface document in the course of work of MATLAB platform.
6. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load according to claim 1,
It is characterized in that, described program extraction unit data, sit two endpoints of each unit of windward side unit record of selection
Scale value and sectional dimension, to calculate the area of each unit center point coordinate, element length, unit.
7. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load according to claim 1,
It is characterized in that, described simulate Wind Velocity History sample, programmed with MATLAB, using linear filtering method or harmonic and reactive detection
Method, wave spectrum use Davenport wind speed spectrum or Kaimal wind speed spectrum analog wind speed, directly obtain the system that wind speed changes over time
Column data, as Wind Velocity History sample.
8. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load according to claim 1,
It is characterized in that, the wind pressure time-histories sample out can calculate Wind Velocity History sample according to the conversion formula of wind pressure and wind speed
Wind pressure time-histories sample out.
9. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load according to claim 1,
It is characterized in that, the concentrated force time-histories sample, the wind load of each windward side unit is loaded in the unit of windward side
It is concentrated force in heart coordinate points, wind load changes over time, and is multiplied by corresponding unit area by wind pressure time-histories sample
Obtain the concentrated force time-histories sample of unit.
10. a kind of applying method of high-voltage power transmission tower space rigid finite element model wind load according to claim 2,
It is characterized in that, the wind load of the windward side unit, the wind load of each windward side unit is equal to wind pressure multiplied by each
Windward side cellar area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811047744.3A CN109116753A (en) | 2018-09-10 | 2018-09-10 | A kind of applying method of high-voltage power transmission tower space rigid finite element model wind load |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811047744.3A CN109116753A (en) | 2018-09-10 | 2018-09-10 | A kind of applying method of high-voltage power transmission tower space rigid finite element model wind load |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109116753A true CN109116753A (en) | 2019-01-01 |
Family
ID=64858282
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811047744.3A Pending CN109116753A (en) | 2018-09-10 | 2018-09-10 | A kind of applying method of high-voltage power transmission tower space rigid finite element model wind load |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109116753A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109885955A (en) * | 2019-03-01 | 2019-06-14 | 上海振华重工(集团)股份有限公司 | Wind load calculation method and computing system |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102968542A (en) * | 2012-12-12 | 2013-03-13 | 河北省电力勘测设计研究院 | Method for performing structure analysis on steel pylon by using ANSYS software |
| JP2014199635A (en) * | 2013-03-11 | 2014-10-23 | 株式会社リコー | Calculation device and simulation system |
| CN105354377A (en) * | 2015-10-29 | 2016-02-24 | 中国电力科学研究院 | Method for determining fluctuation wind induced vibration load of power transmission tower |
| CN105738143A (en) * | 2016-03-10 | 2016-07-06 | 国网福建省电力有限公司 | Cement electric pole equivalent wind load loading device and experiment method |
| CN105740548A (en) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | Method for calculating wind vibration of power transmission line under random wind load |
| CN105740549A (en) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | Wind vibration simulation method for 220kV power transmission line tower wire coupling system |
| CN108073758A (en) * | 2017-10-17 | 2018-05-25 | 广东电网有限责任公司电力科学研究院 | A kind of Transmission Tower-line wind-induced vibration response emulation mode and device |
-
2018
- 2018-09-10 CN CN201811047744.3A patent/CN109116753A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102968542A (en) * | 2012-12-12 | 2013-03-13 | 河北省电力勘测设计研究院 | Method for performing structure analysis on steel pylon by using ANSYS software |
| JP2014199635A (en) * | 2013-03-11 | 2014-10-23 | 株式会社リコー | Calculation device and simulation system |
| CN105354377A (en) * | 2015-10-29 | 2016-02-24 | 中国电力科学研究院 | Method for determining fluctuation wind induced vibration load of power transmission tower |
| CN105740548A (en) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | Method for calculating wind vibration of power transmission line under random wind load |
| CN105740549A (en) * | 2016-02-01 | 2016-07-06 | 西安交通大学 | Wind vibration simulation method for 220kV power transmission line tower wire coupling system |
| CN105738143A (en) * | 2016-03-10 | 2016-07-06 | 国网福建省电力有限公司 | Cement electric pole equivalent wind load loading device and experiment method |
| CN108073758A (en) * | 2017-10-17 | 2018-05-25 | 广东电网有限责任公司电力科学研究院 | A kind of Transmission Tower-line wind-induced vibration response emulation mode and device |
Non-Patent Citations (1)
| Title |
|---|
| 贺博，等: "强台风下高压输电线路塔_线耦联体系的力学行为仿真分析二：动力响应分析", 《高压电器》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109885955A (en) * | 2019-03-01 | 2019-06-14 | 上海振华重工(集团)股份有限公司 | Wind load calculation method and computing system |
| CN109885955B (en) * | 2019-03-01 | 2023-06-02 | 上海振华重工(集团)股份有限公司 | Wind load calculation method and calculation system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103020426B (en) | A kind of short-cut method of rectangular slab Plate with Inclined Center Crack propagation life of fatigue prediction | |
| CN102968811B (en) | Real-time motion simulation method for guide wire in minimally invasive vessel interventional operation | |
| CN106202732B (en) | A kind of axial modification method of involute spur gear pair and special parameters CAD system matched with its | |
| CN103353907A (en) | Calculating method of ultimate strength checking of connection of flange and bolt | |
| CN103280137B (en) | 220KV power transformer 3D modelling system | |
| CN106650046A (en) | Method for obtaining unsteady characteristic of air flow field in ship | |
| CN108106594A (en) | Building deformation method for real-time measurement | |
| CN109583377B (en) | Control method and device for pipeline model reconstruction and upper computer | |
| CN108073758B (en) | Simulation method and device for wind-induced vibration response of power transmission tower line | |
| CN103577651A (en) | Method for generating and processing isoline cloud charts with sections in any shapes intelligently in batch mode | |
| CN104268322A (en) | Boundary processing technology of WENO difference method | |
| CN105302994A (en) | Finite element simulation method of beam_plate shell structure | |
| CN106650086A (en) | Indirect air cooling system numerical simulation platform based on FLUENT software | |
| CN108756292A (en) | A kind of building synchronization of jacking up control method and system based on technology of Internet of things | |
| CN109116753A (en) | A kind of applying method of high-voltage power transmission tower space rigid finite element model wind load | |
| CN107422376B (en) | A kind of non-stationary non-gaussian earthquake motion time history analogy method | |
| CN103927409A (en) | Finite element analysis based concrete chimney stress analysis system | |
| CN106354929B (en) | Bearing structure load transfer path method for visualizing based on stiffness variation principle | |
| CN116226728A (en) | Floating structure frequency response function identification method based on single regular wave excitation | |
| CN104298806A (en) | Hydropower station dynamic property computer-assisted testing method | |
| CN109446658A (en) | A kind of determination method of titanium alloy taper impression constraint factor and representative strain | |
| CN113609734A (en) | Visualization method for thickness of simulated oil film of sliding bearing water turbine | |
| CN103399974B (en) | Quantize the method comparing random vibration emulated data and experimental data | |
| CN104699912A (en) | Strength calculation method of hub and spindle connecting bolt of wind generator set | |
| CN109559377A (en) | A method of it is cut using longitude and latitude polygon and generates three-dimensional map model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190101 |