CN110069843A - Based on the considerations of the Structural Topology Optimization Design method of ANSYS row wave effect - Google Patents
Based on the considerations of the Structural Topology Optimization Design method of ANSYS row wave effect Download PDFInfo
- Publication number
- CN110069843A CN110069843A CN201910293303.XA CN201910293303A CN110069843A CN 110069843 A CN110069843 A CN 110069843A CN 201910293303 A CN201910293303 A CN 201910293303A CN 110069843 A CN110069843 A CN 110069843A
- Authority
- CN
- China
- Prior art keywords
- rod piece
- wave effect
- row wave
- finite element
- stress
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of Structural Topology Optimization Design method based on the considerations of ANSYS row wave effect, is related to arbitrary excitation flowering structure topological optimization technology field.Steps are as follows by the present invention: step 1: considering the space earthquake motion time-history curves of row wave effect with obtaining non-smooth uniform field, be translated into arbitrary excitation acceleration value；Step 2: establishing finite element structure model；Input stimulus acceleration on the mold obtains rod piece internal force and rod piece stress；Step 3: establishing topological optimization model: rod piece stress being met into rod piece allowable stress be set to constraint condition finite element structure model is optimized；Step 4: repeating step 2- step 3, traverse the data in sectional area set, obtain new bar cross section productive set and close；Step 5: time-history analysis being carried out to data obtained in step 4, repeats step 3- step 4, until bar cross section product reaches and meets row wave effect and structure gross mass reaches minimum；This method saves steel, and the purpose for protecting structure not to be destroyed.
Description
Technical field
The present invention relates to arbitrary excitation flowering structure topological optimization technology fields, more particularly to one kind is based on the considerations of ANSYS
The Structural Topology Optimization Design method of row wave effect.
Background technique
Earthquake is a kind of great Sudden Natural Disasters of destructiveness, and when earthquake occurs, ground motion process is complicated, earthquake
The wave of generation is influenced to pass to bottom surface by seismic source mechanism and ground geological conditions, and each particle in ground is by different shake
Emotionally condition, or due to having certain distance between ground particle and focus, cause to generate on the seismic wave time in the presence of certain poor
It is different, there is different reaction results.The variation of the variability major embodiment of complicated earthquake motion over time and space, including traveling wave
Effect, local site effect and arbitrary coherency etc., it is the most significant with row wave effect." seismic design provision in building code "
(GB50011-2010) it is also indicated that in, calculates analysis when to carry out Random seismic field more than a certain range of span degree space structure
When, row wave effect should be specialized in.Currently, many scholars further investigate row wave effect, traveling wave is all proved
Effect is with the increase of view velocity of wave on the very big influence of the tools such as structural internal force, displacement, and proposition will consider row wave effect, but not
It is specific to propose to resist or reduce in the effective measures of row wave effect.Only consider to increase according to Traditional Thinking to be influenced by row wave effect
The problems such as big bar cross section product will cause Internal Force Redistribution, low so as to cause rod piece utilization efficiency, therefore, in structure design
Not only consider row wave effect, and specifically proposes economically and efficiently to solve due to a series of adverse effects of row wave effect bring
Measure becomes a urgent problem.
Summary of the invention
The technical problem to be solved by the present invention is in view of the above shortcomings of the prior art, provide a kind of examining based on ANSYS
Consider the Structural Topology Optimization Design method of row wave effect, this method uses time history analysis method by manually generated seismic wave time-history curves
It is converted to power excitation to be applied in structure, according to fully stressed criterion, structure reaches in each rod piece stress of row wave effect flowering structure
When fully stress, architecture quality reaches most light, arrived i.e. saving steel, and the purpose for protecting structure not to be destroyed.
In order to solve the above technical problems, the technical solution used in the present invention is:
The present invention provides a kind of Structural Topology Optimization Design method based on the considerations of ANSYS row wave effect, including following
Step:
Step 1: according to the classification of providing fortification against earthquakes of earthquake resistant code and local building, seismic fortification intensity, site condition, utilizing
Trigonometric series method considers coherent function and the artificial synthesized non-smooth uniform field of envelope function factor is added the sky of row wave effect
Between earthquake motion time history curve, and convert it into the arbitrary excitation acceleration value with row wave effect.
Step 2: establishing finite element structure model using finite element analysis software；It is required according to earthquake resistant code and engineering
Actual requirement is arranged bar cross section productive set and closesWherein i represents rod piece number, is first rod piece selection
Arbitrary section accumulates Ai, apply static load for finite element structure model and carry out static analysis, judge finite element structure model
Reasonability；Unreasonable reason is searched if unreasonable and re-establishes finite element structure model, if rationally calling this limited
Meta-model simultaneously reads in the step 1 arbitrary excitation acceleration generated with row wave effect in ANSYS, utilizes time history analysis method
Rod piece internal force and rod piece stress σ are obtained to the seismic response analysis that finite element structure model accounts for row wave effecti。
Step 3: topological optimization model is established, with bar cross section product for design variable, according to quality of materials density p and rod piece
Length l calculates the minimum gross mass minW of structural member, the formula of minimum gross mass minW are as follows:
The rod piece stress that consideration row wave effect obtains is met into rod piece allowable stress and is set to constraint condition to finite element structure
Model optimizes；Specific steps are as follows: bar is obtained by time history analysis method using in ANSYS post processor extraction step 2
Part stress value σi, rod piece stress value σ will be obtainediIt is compared with rod piece allowable stress value σ ', obtained ratioIf ratio?
Within the scope of 0.8-0.9, then the bar area be not changed, if ratioNot in this range, then by ratioMultiplied by bar
The sectional area A of part initial setting upi, new bar cross section product A ' is obtained, using the EMODIF order in finite element software by rod piece
Sectional area is by AiIt changes into sectional area AiBar cross section productive set close in be greater than A ' and A adjacent theretom。
Step 4: repeat step 2 to step 3, traverse bar cross section productive set close in all data, after traverse newly
Bar cross section productive set closes
Step 5: to a obtained in step 4mTime-history analysis is carried out, new rod piece internal force and rod piece stress is obtained；It repeats to walk
Rapid 3 and step 4, until reaching structural member sectional area meets row wave effect and structure gross mass reaches minimum；Structure is total
It is the optimum results for considering row wave effect that quality, which reaches the smallest rod piece,.
The beneficial effects of adopting the technical scheme are that one kind provided by the invention is based on the considerations of ANSYS
The Structural Topology Optimization Design method of row wave effect, this method can obtain structural section after the iterative cycles within 20 times
Product optimum results, this result not only consider that row wave effect bring influences, and provide and subtracted by way of changing bar cross section product
The measure that light row wave effect bring influences, obtained structure gross mass are still the smallest under the conditions of meeting structural stress.This
Manually generated seismic wave time-history curves are converted to power excitation using time history analysis method and are applied in structure by method, according to fully stress
Criterion, for structure when each rod piece stress of row wave effect flowering structure reaches fully stress, architecture quality reaches most light, arrived and saves
Steel are saved, and structure can be protected not to be destroyed under seismic response especially row wave effect effect；It is examined in structure optimization
Consider, increasing is influenced serious bar cross section product by row wave effect, effectively avoids row wave effect, solves the considerations of saying in specification row
Wave effect.
Detailed description of the invention
Fig. 1 is method flow diagram provided in an embodiment of the present invention.
Specific embodiment
With reference to the accompanying drawings and examples, specific embodiments of the present invention will be described in further detail.Implement below
Example is not intended to limit the scope of the invention for illustrating the present invention.
The method of the present embodiment is as described below.
The present invention provides a kind of Structural Topology Optimization Design method based on the considerations of ANSYS row wave effect, such as Fig. 1 institute
Show, comprising the following steps:
Step 1: according to the classification of providing fortification against earthquakes of earthquake resistant code and local building, seismic fortification intensity, site condition, utilizing
Trigonometric series method considers coherent function and the artificial synthesized non-smooth uniform field of envelope function factor is added the sky of row wave effect
Between earthquake motion time history curve, and convert it into the arbitrary excitation acceleration value with row wave effect.
The specific space earthquake motion time-history curves mode that generates is as follows, employs acceleration power spectrum according to ground each point earthquake
Spending Jacobian matrix S (v ω) indicates are as follows:
In formula, Svv(ω) and Svj(ω) is respectively Power spectral density and cross-spectral density function, wherein v ∈
N, j ∈ n；S (v ω) is subjected to Cholesky decomposition, obtains triangle complex matrix L (v ω) and Hermite matrix LH(v ω) multiplies
Product:
S (v ω)=L (v ω) LH(vω)
Then generate the steady space earthquake motion time history u in the time domain scale at n-th of pointK(t) such as following formula:
In formula, h=1,2 ..., n, ω are circular frequency, and t is time, ωl=l △ ω, △ ω=ωN/ N, △ ω are angular frequency
Rate, ωlIt is expressed as first of discrete frequency, ωNIndicate upper cut-off frequency；Be generally evenly distributed in [0,2 π) random phase
Parallactic angle, h and l represent the points of the earthquake motion on board different directions, Aγ(ωl) and θγ(ωl) be time-history curves amplitude and phase
Angle can be represented by the formula:
Wherein LγAnd LγThe value obtained after Cholesky decomposition is carried out for S (v ω)；
Step 2: establishing finite element structure model using finite element analysis software；It is required according to earthquake resistant code and engineering
Actual requirement is arranged bar cross section productive set and closesWherein i represents rod piece number, is first rod piece choosing
Take arbitrary section product Ai, apply static load for finite element structure model and carry out static analysis, judge finite element structure model
Reasonability；Unreasonable reason is searched if unreasonable and re-establishes finite element structure model, if rationally in finite element
The arbitrary excitation acceleration generated with row wave effect of input step 1 on structural model, using time history analysis method to limited
The seismic response analysis that meta structure model accounts for row wave effect obtains rod piece internal force and rod piece stress σi；
Step 3: topological optimization model is established, with bar cross section product for design variable, according to quality of materials density p and rod piece
Length l calculates the minimum gross mass W of structural member, the formula of minimum gross mass minW are as follows:
The rod piece stress that consideration row wave effect obtains is met into rod piece allowable stress and is set to constraint condition to finite element structure
Model optimizes；Specific steps are as follows: bar is obtained by time history analysis method using in ANSYS post processor extraction step 2
Part stress value σi, rod piece stress value σ will be obtainediIt is compared with rod piece allowable stress value σ ', obtained ratioIf ratio
Within the scope of 0.8-0.9, then the bar area be not changed, if ratioNot in this range, then by ratioMultiplied by bar
The sectional area A of part initial setting upi, new bar cross section product A ' is obtained, using the EMODIF order in finite element software by rod piece
Sectional area is by AiIt changes into sectional area AiBar cross section productive set close in be greater than A ' and A adjacent theretom.By bar cross section product
Sectional area in set is arranged according to ascending sequence, finds sectional area that is close with sectional area A ' and being greater than A '；
Step 4: repeat step 2 to step 3, traverse bar cross section productive set close in all data, after traverse newly
Bar cross section productive set closes
Step 5: to a obtained in step 4mTime-history analysis is carried out, new rod piece internal force and rod piece stress is obtained；It repeats to walk
Rapid 3 and step 4, until reaching structural member sectional area meets row wave effect and structure gross mass reaches minimum；Structure is total
It is the optimum results for considering row wave effect that quality, which reaches the smallest rod piece,.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify to technical solution documented by previous embodiment, or some or all of the technical features are equal
Replacement；And these are modified or replaceed, model defined by the claims in the present invention that it does not separate the essence of the corresponding technical solution
It encloses.
Claims (1)
1. a kind of Structural Topology Optimization Design method based on the considerations of ANSYS row wave effect, it is characterised in that: including following step
It is rapid:
Step 1: according to the classification of providing fortification against earthquakes of earthquake resistant code and local building, seismic fortification intensity, site condition, utilizing triangle
Series Method considers row wave effect spatially coherent function and the artificial synthesized non-smooth uniform field of envelope function factor is added
Time-history curves are shaken, and convert it into the arbitrary excitation acceleration value with row wave effect；
Step 2: establishing finite element structure model using finite element analysis software；According to earthquake resistant code requirement and the reality of engineering
It is required that setting bar cross section productive set closesWherein i represents rod piece number, first chooses and appoints for rod piece
Anticipate sectional area Ai, apply static load for finite element structure model and carry out static analysis, judge the conjunction of finite element structure model
Rationality；Unreasonable reason is searched if unreasonable and re-establishes finite element structure model, if rationally calling the finite element
Model simultaneously reads in the step 1 arbitrary excitation acceleration generated with row wave effect in ANSYS, utilizes time history analysis method pair
The seismic response analysis that finite element structure model accounts for row wave effect obtains rod piece internal force and rod piece stress σi；
Step 3: topological optimization model is established, with bar cross section product for design variable, according to quality of materials density p and rod length
L calculates the minimum gross mass W of structural member, the formula of minimum gross mass minW are as follows:
The rod piece stress that consideration row wave effect obtains is met into rod piece allowable stress and is set to constraint condition to finite element structure model
It optimizes；Specific steps are as follows: answered using rod piece is obtained by time history analysis method in ANSYS post processor extraction step 2
Force value σi, rod piece stress value σ will be obtainediIt is compared with rod piece allowable stress value σ ', obtained ratioIf ratioIn 0.8-
In 0.9 range, then the bar area be not changed, if ratioNot in this range, then by ratioIt is initial multiplied by rod piece
The sectional area A of settingi, the long-pending A ' of new bar cross section is obtained, it is using the EMODIF order in finite element software that bar cross section is long-pending
By AiIt changes into sectional area AiBar cross section productive set close in be greater than A ' and A adjacent theretom；
Step 4: repeating step 2 to step 3, traverse all data in the conjunction of bar cross section productive set, new rod piece after being traversed
Sectional area set
Step 5: to a obtained in step 4mTime-history analysis is carried out, new rod piece internal force and rod piece stress is obtained；Repeat step 3 and
Step 4, until reaching structural member sectional area meets row wave effect and structure gross mass reaches minimum；Structure gross mass reaches
It is the optimum results for considering row wave effect to the smallest rod piece.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910293303.XA CN110069843B (en) | 2019-04-12 | 2019-04-12 | ANSYS-based structural topology optimization design method considering traveling wave effect |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910293303.XA CN110069843B (en) | 2019-04-12 | 2019-04-12 | ANSYS-based structural topology optimization design method considering traveling wave effect |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110069843A true CN110069843A (en) | 2019-07-30 |
| CN110069843B CN110069843B (en) | 2023-04-18 |
Family
ID=67367578
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910293303.XA Active CN110069843B (en) | 2019-04-12 | 2019-04-12 | ANSYS-based structural topology optimization design method considering traveling wave effect |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110069843B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110565500A (en) * | 2019-08-30 | 2019-12-13 | 中铁大桥勘测设计院集团有限公司 | Intelligent design method for section of steel truss girder rod piece |
| CN112528375A (en) * | 2020-12-10 | 2021-03-19 | 扬州大学 | Method for analyzing seismic response of long-span bridge by considering traveling wave effect based on large-mass method |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103559383A (en) * | 2013-09-30 | 2014-02-05 | 上海交通大学 | Method for predicting and evaluating strong earthquake resistance of nuclear power station breakwater |
-
2019
- 2019-04-12 CN CN201910293303.XA patent/CN110069843B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103559383A (en) * | 2013-09-30 | 2014-02-05 | 上海交通大学 | Method for predicting and evaluating strong earthquake resistance of nuclear power station breakwater |
Non-Patent Citations (8)
| Title |
|---|
| S.Z.SHEN: "A Review of the Development of Spatial Structures in China" * |
| 全伟;李宏男;: "大跨斜拉桥多维多点地震激励减震控制方法分析" * |
| 吴长;王秀丽;周锟;: "开孔式双钢管约束屈曲支撑在双层网壳的减震控制分析" * |
| 张渊;赵海艳;宁志刚;刘巍;: "考虑地震行波效应的145m跨三心圆网壳结构时程分析" * |
| 李正坤;李晓文;孙晓岭;程元虎;: "超大跨长距离输煤栈桥考虑地震动行波效应影响分析研究" * |
| 李永华;朱芊;文明;: "竖向多点激励下输电塔线体系导(地)线及绝缘子动力响应分析" * |
| 江宜城;杨德喜;李黎;胡亮;: "LRB隔震桥梁空间变异性地震随机响应分析" * |
| 赵博;王元清;陈志华;石永久;: "行波效应下对称多跨大跨结构的随机地震响应研究" * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110565500A (en) * | 2019-08-30 | 2019-12-13 | 中铁大桥勘测设计院集团有限公司 | Intelligent design method for section of steel truss girder rod piece |
| CN110565500B (en) * | 2019-08-30 | 2021-05-25 | 中铁大桥勘测设计院集团有限公司 | Intelligent design method for section of steel truss girder rod piece |
| CN112528375A (en) * | 2020-12-10 | 2021-03-19 | 扬州大学 | Method for analyzing seismic response of long-span bridge by considering traveling wave effect based on large-mass method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110069843B (en) | 2023-04-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jones et al. | Convection-driven dynamos in a rotating plane layer | |
| CN110069843A (en) | Based on the considerations of the Structural Topology Optimization Design method of ANSYS row wave effect | |
| Kupka et al. | Total-variation-diminishing implicit–explicit Runge–Kutta methods for the simulation of double-diffusive convection in astrophysics | |
| CN104268322B (en) | A kind of BORDER PROCESSING technology of WENO difference methods | |
| CN109460622B (en) | Complete explicit power time course analysis method for large-scale building structure | |
| Grigoriadis et al. | LES of the flow past a rectangular cylinder using the immersed boundary concept | |
| Lee | Application of Chebyshev-tau method to the free vibration analysis of stepped beams | |
| Myers et al. | A 4th-Order Particle-in-Cell Method with Phase-Space Remapping for the Vlasov--Poisson Equation | |
| Miesch | Numerical modeling of the solar tachocline. II. Forced turbulence with imposed shear | |
| Hill et al. | Subharmonic resonance of short internal standing waves by progressive surface waves | |
| Hossain et al. | Three-dimensional compressible hydrodynamic convection in the Sun and stars | |
| Drui et al. | Experimenting with the p4est library for AMR simulations of two-phase flows | |
| Yan et al. | Large-scale three-dimensional phase-field simulations for phase coarsening at ultrahigh volume fraction on high-performance architectures | |
| Kurgansky | Inertial instability of the Kolmogorov flow in a rotating stratified fluid | |
| Sukhinov et al. | Parallel algorithms for modelling of suspended particles motion in channel for large peclet number | |
| Woodward et al. | The piecewise-parabolic Boltzmann advection scheme (PPB) applied to multifluid hydrodynamics | |
| Doltsinis et al. | Studies on parallel processing for coupled field problems | |
| Shashkov et al. | A composite scheme for gas dynamics in Lagrangian coordinates | |
| Audusse et al. | Parallelization of a relaxation scheme modelling the bedload transport of sediments in shallow water flow | |
| Brewin | Einstein-Bianchi system for smooth lattice general relativity. I. The Schwarzschild spacetime | |
| Tzanos | Higher-order differencing method with a multigrid approach for the solution of the incompressible flow equations at high Reynolds numbers | |
| Shitrit et al. | An algebraic multigrid solver for transonic flow problems | |
| Gholizadeh et al. | Seismic design of double layer grids by neural networks | |
| Vlădulescu | Design optimization of a civil structure using numerical simulations | |
| Corbà et al. | On the adiabatic subtraction of cosmological perturbations |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |