CN110096808A - A kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load - Google Patents
A kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load Download PDFInfo
- Publication number
- CN110096808A CN110096808A CN201910364967.0A CN201910364967A CN110096808A CN 110096808 A CN110096808 A CN 110096808A CN 201910364967 A CN201910364967 A CN 201910364967A CN 110096808 A CN110096808 A CN 110096808A
- Authority
- CN
- China
- Prior art keywords
- vault
- load
- value
- linear
- stability
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 52
- 238000004088 simulation Methods 0.000 title claims abstract description 28
- 238000012797 qualification Methods 0.000 claims abstract description 8
- 230000003321 amplification Effects 0.000 claims abstract description 5
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 5
- 238000006073 displacement reaction Methods 0.000 claims description 22
- 229910000831 Steel Inorganic materials 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 239000010959 steel Substances 0.000 claims description 4
- 238000000034 method Methods 0.000 abstract description 24
- 238000009826 distribution Methods 0.000 abstract description 10
- 238000013461 design Methods 0.000 description 14
- 230000000694 effects Effects 0.000 description 9
- 230000009286 beneficial effect Effects 0.000 description 5
- 240000002853 Nelumbo nucifera Species 0.000 description 4
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 4
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 4
- 150000001875 compounds Chemical class 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000002787 reinforcement Effects 0.000 description 3
- 238000005452 bending Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000005728 strengthening Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012913 prioritisation Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Classifications
-
- 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)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of spherical shell vault Finite Element Simulation Analysis methods with ribbing under multiple spot load, comprising the following steps: establishes analysis finite element geometrical model；Concentrfated load gradually is inputted on the vault of simulation, obtains its linear instability modes；It obtains vault nonlinear stability and analyzes initial load value；Concentrfated load is inputted on vault, obtains first linear stable property coefficient, and compare with assessed value；When linear stable property coefficient is less than assessed value, local enhancement component is added to it and is considered as qualification more than or equal to assessed value；Model when according to qualification by the concentrfated load equal proportion amplification setting multiplying power on original input vault and is analyzed, and the collapsing load value of vault is obtained；Calculate the coefficient of stability.The present invention simulates the practical mode and deformation distribution of spherical shell vault with ribbing under the conditions of different loads, analyzes the stability of spherical shell vault with ribbing under the conditions of different loads, and obtains the anti-unstability method of spherical shell vault with ribbing under the conditions of different loads.
Description
Technical field
The present invention relates to Structural Dynamics technical fields, and in particular to the spherical shell vault with ribbing under a kind of multiple spot load
Finite Element Simulation Analysis method.
Background technique
Spherical shell vault with ribbing is domestic distinctive tank structure, simple compared to truss-like or column brace type residual oil in tank,
Good economy performance.But its bearing capacity is limited, and current standard only accounts for spherical shell when calculating and bears uniform load and can not check
Additional multiple spot load bring influences thereon.In practical projects, vault will not only bear certain external pressure effect, certain feelings
Condition still suffers from biggish concentrated force effect, the even load of multiple spot distribution.Such as: since internal pressure storage tank occurs needing to maintain
When higher gas phase seals, the high-precision starting being arranged in tank deck is needed to adjust tank deck sealing system, which need to set in tank deck
It sets pipeline and valve group bring is up to the additional multiple spot load of more than ten tons, be typically designed item far more than national standard spherical shell vault with ribbing
Part, if careful consideration is not added, tank deck is being concentrated on and is being collapsed risk under external pressure collective effect there are overall collapse.
When carrying out spherical shell vault intensive analysis with ribbing, current design method has ignored tank where concentrfated load at present
The design strength distribution on top has ignored the great variety that roof stability generates after concentrfated load applies.For concentrfated load without
Method is checked, and stress, linear and nonlinear stability especially in the case where currently needing to concentrfated load lack safe guarantor
Barrier.
Summary of the invention
Technical problem to be solved by the invention is to provide the spherical shell vault finite element simulations with ribbing under a kind of multiple spot load
Analysis method simulates the practical mode and deformation distribution of spherical shell vault with ribbing under the conditions of different loads, analyzes different loads condition
Under spherical shell vault with ribbing stability, and obtain the anti-unstability method of spherical shell vault with ribbing under the conditions of different loads.
The technical scheme to solve the above technical problems is that
A kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load, comprising the following steps:
S1: analysis finite element geometrical model is established；
S2: gradually inputting concentrfated load on the vault to simulation on the basis of the analysis finite element geometrical model, obtains
The linear instability modes of the vault；
S3: the initial load value of the vault nonlinear stability analysis is obtained according to the linear instability modes；
S4: according to the initial load value, inputting the concentrfated load on the vault, obtains first linear stable
Property coefficient, and the linear stable property coefficient is compared with assessed value；
S5: when the linear stable property coefficient is less than institute's rating value, local enhancement component, repeating said steps are added
S4 is until the linear stable property coefficient is greater than or equal to institute's rating value and is considered as qualification；
S6: after the linear stable property coefficient is qualified, model when according to qualification will be on original input vault
The concentrfated load equal proportion amplification setting multiplying power carries out geometrical nonlinear analysis, obtains the collapsing load value of the vault；
S7: the coefficient of stability of the vault is calculated according to the initial load value and the collapsing load value.
The beneficial effects of the present invention are: the present invention initially sets up finite element geometrical model, simulation arch is obtained after preliminary analysis
The linear instability modes on top, and the instability status of vault under different concentrfated loads is analyzed, obtain initial load when vault collapses
Value；Then on the basis of vault linear instability modes, concentrfated load is inputted to vault step by step, obtains the Nonlinear Instability of vault
Mode, and analyze the collapsing load value for obtaining vault when Nonlinear Instability collapses；Finally according to the initial load of vault obtained
Charge values and collapsing load value calculate the coefficient of stability of simulation vault, so that analysis obtains the scheme of the anti-unstability of vault.The present invention
The practical mode and deformation distribution for simulating spherical shell vault with ribbing under the conditions of different loads, analyze spherical shell with ribbing under the conditions of different loads
The stability of vault, and obtain the anti-unstability method of spherical shell vault with ribbing under the conditions of different loads.
Based on the above technical solution, the present invention can also be improved as follows.
Further, the step S3 the following steps are included:
S31: according to the xyz three direction displacement for inputting the corresponding point of load of the concentrfated load, each point of load is calculated
Displacement vector values, and to multiple displacement vector values be ranked up obtain its maximum value, find maximum institute's displacement vector
It is worth corresponding strain maximum point；
S32: the maximum displacement vector values are converted after obtaining conversion factor according to 1/300 total span, root
Institute according to the conversion factor to cell nodes whole in the analysis finite element geometrical model under the linear instability modes
Displacement vector value is converted, and the shift value after obtaining all cell node conversions is as the initial load value.
Beneficial effect using above-mentioned further scheme is to calculate the position vector value of each point of load first, and to multiple
Displacement vector values are ranked up, and obtain maximum displacement vector value, and find the corresponding strain maximum point of maximum displacement vector value；So
Maximum displacement vector value is converted according to 1/300 total span afterwards, obtains conversion factor；Finally, according to above-mentioned conversion system
Several pairs of whole cell nodes are converted, and initial load value is obtained.
Further, the step S6 the following steps are included:
S61: first nonlinear stability property coefficient is obtained according to geometrical nonlinear analysis, and by the nonlinear stability
Coefficient and non-linear setting characteristic value compare, and judge whether the nonlinear stability property coefficient restrains, the nonlinear stability
Coefficient shows its convergence when being greater than or equal to the non-linear setting characteristic value；
S62: when the nonlinear stability property coefficient is less than the non-linear setting characteristic value, to the nonlinear stability
The corresponding point of load of property coefficient adds the local enhancement component and is reinforced, and repeating said steps S61 obtains described non-
Linear stable property coefficient, until the nonlinear stability property coefficient is greater than or equal to the non-linear setting characteristic value, obtaining should
The corresponding collapsing load value of the nonlinear stability property coefficient.
Beneficial effect using above-mentioned further scheme is to obtain first nonlinear stability according to geometrical nonlinear analysis
Property coefficient, and the nonlinear stability property coefficient and non-linear setting characteristic value are compared, judge the nonlinear stability system
Whether number restrains, the nonlinear stability property coefficient corresponding collapsing load value when obtaining convergence.
Further, the non-linear characteristic value that sets in the step S61 is 3.0-3.5.
Further, the coefficient of stability in the step S7 is the collapsing load value/initial load value.
Beneficial effect using above-mentioned further scheme is the coefficient of stability for calculating simulation vault, analyzes different loads item
The stability of vault under part, and obtain the anti-unstability method of spherical shell vault with ribbing under the conditions of different loads.
Further, institute's rating value in the step S4 is 4.2.
Further, the local enhancement component in the step S5 is steel rib.
Beneficial effect using above-mentioned further scheme is to increase the stability of vault.
Detailed description of the invention
Fig. 1 is one of the structural schematic diagram that vault is simulated in the present invention；
Fig. 2 is the second structural representation that vault is simulated in the present invention；
Fig. 3 is the simulation drawing that vault load is reinforced for the first time in the present invention；
Fig. 4 is second of simulation drawing reinforced of vault load in the present invention；
Fig. 5 is the multiple strengthened simulation drawing of vault load in the present invention；
Fig. 6 is the simulation drawing that deformation occurs after vault reinforcement in the present invention；
Fig. 7 is the simulation drawing that vault increases local reinforcement in the present invention；
Fig. 8 is the structural schematic diagram of Spherical Tank Shell Plates with ribbing in the prior art.
In attached drawing, parts list represented by the reference numerals are as follows:
1, vault, 2, local enhancement component, 3, top plate, 4, rib.
Specific embodiment
Principles and features of the present invention are described below in conjunction with drawings and the specific embodiments, example is served only for solving
The present invention is released, is not intended to limit the scope of the present invention.
As shown in Figures 1 to 7, the present invention provides the spherical shell vault with ribbing Finite Element Simulation Analysis side under a kind of multiple spot load
Method, comprising the following steps:
S1: analysis finite element geometrical model is established；
S2: gradually inputting concentrfated load on the basis of the analysis finite element geometrical model on the vault of simulation 1, obtains
The linear instability modes of the vault 1；
S3: obtaining the initial load value of 1 nonlinear stability of the vault analysis according to the linear instability modes, including
Following steps:
S31: according to the xyz three direction displacement for inputting the corresponding point of load of the concentrfated load, each point of load is calculated
Displacement vector values, and to multiple displacement vector values be ranked up obtain its maximum value, find maximum institute's displacement vector
It is worth corresponding strain maximum point；
S32: the maximum displacement vector values are converted to obtain conversion factor according to 1/300 total span (described
Absolute value/300 of conversion factor=maximum displacement vector values) after, according to the conversion factor to the analysis finite element
(cell node refers to is divided into limited unit, Mei Gedan for analysis finite element geometrical model to whole cell nodes in geometrical model
Member represents a cell node) displacement vector values under the linear instability modes are converted, it obtains described in whole
Shift value after cell node conversion is as the initial load value；
S4: according to the initial load value, inputting the concentrfated load on the vault 1, obtains first linearly surely
Qualitative coefficient, and the linear stable property coefficient is compared with assessed value, institute's rating value in the step S4 is
4.2；
S5: it when the linear stable property coefficient is less than institute's rating value, needs to continue to choose full model bending deformation most
General goal (i.e. the point of load of 1 unstability of vault) adds local enhancement component 2, and repeating said steps S4 is until the Linear Stability
Coefficient is greater than or equal to institute's rating value and is considered as qualification；
S6: after the linear stable property coefficient is qualified, model when according to qualification will be on original input vault 1
The concentrfated load equal proportion amplification setting multiplying power carries out geometrical nonlinear analysis, obtains the collapsing load value of the vault 1, wraps
Include following steps:
S61: first nonlinear stability property coefficient is obtained according to geometrical nonlinear analysis, and by the nonlinear stability
Coefficient and non-linear setting characteristic value compare, and the non-linear characteristic value that sets is 3.0-3.5, the non-linear setting characteristic value
3.0,3.3 and 3.5 can be taken, judges whether the nonlinear stability property coefficient restrains, the nonlinear stability property coefficient is greater than
Or show its convergence when equal to the non-linear setting characteristic value；
S62: when the nonlinear stability property coefficient is less than the non-linear setting characteristic value, to the nonlinear stability
The corresponding point of load of property coefficient adds the local enhancement component 2 and is reinforced, and repeating said steps S61 obtains described non-
Linear stable property coefficient, until the nonlinear stability property coefficient is greater than or equal to the non-linear setting characteristic value, obtaining should
The corresponding collapsing load value of the nonlinear stability property coefficient；
S7: the coefficient of stability of the vault 1, the stabilization are calculated according to the initial load value and the collapsing load value
Coefficient is the collapsing load value/initial load value.
Preferably, in the present invention, the local enhancement component 2 in the step S5 is steel rib, increases the strong of vault 1
Degree, to increase the stability of vault 1.
It should be noted that the setting multiplying power of the amplification of concentrfated load equal proportion described in step S6 is usually at 2.5~50 times
Between.It is selected when practical operation, such as is chosen according to the value of linear stable property coefficient × 1.25 times according to demand.
In above-mentioned steps S4 and step S5, the method that linear newton-La Fuxun iteration both can be used inputs load on vault 1
The method of the nonlinear Newtonian-La Fuxun iteration input load on vault 1 can also be used in lotus, carries out to the structure of vault 1 non-thread
Property analysis, nonlinear analysis engine use arc-length methods, with newton-La Fuxun iteration, when being multiplied to concentrfated load, geometry
The Nonlinear Cumulative of displacement is iterated, when finally finding out Nonlinear Instability and collapsing, locating load level in structure, and root
Ratio is loaded according to concentrfated load, obtains the coefficient of stability under final instability status.
In addition, number involved in Fig. 3 to Fig. 5 is the linear unstability characteristic value of vault 1 under the conditions of different loads.
The present invention initially sets up finite element geometrical model, and the linear instability modes of simulation vault 1 are obtained after preliminary analysis, and
The instability status of vault 1 under the conditions of analysis different loads, obtains initial load value when vault 1 collapses；Then linear in vault 1
On the basis of instability modes, concentrfated load is inputted to vault 1 step by step, obtains the Nonlinear Instability mode of vault, and analyze and obtain
The collapsing load value of vault 1 when Nonlinear Instability collapses；Finally according to the initial load value and collapsing load of vault 1 obtained
Value calculates the coefficient of stability of simulation vault 1, so that analysis obtains the scheme of the anti-unstability of vault 1.The present invention simulates different loads
Under the conditions of vault 1 practical mode and deformation distribution, analyze the stability of vault 1 under the conditions of different loads, and obtain different loads
The anti-unstability method of vault 1 under the conditions of lotus.
In the present invention, latticed vault under the conditions of different loads non-linear setting characteristic value be greater than or equal to 3.5 when, institute
The vault construction designed is more stable, can keep relative stability under the conditions of different loads.
Involved finite element method (FEA, Finite Element Analysis) is close using mathematics in the present invention
As method actual physical system (geometry and load working condition) is simulated.(i.e. using simple and interaction element
Unit), so that it may it is gone to approach the real system of unlimited unknown quantity with the unknown quantity of limited quantity.
In addition, finite element analysis is solved again after replacing challenge with better simply problem.It is regarded as domain is solved
It is made of many small interconnection subdomains for being known as finite element, suitable (better simply) approximate solution is assumed to each unit,
What then this domain of Derivation was total meets condition (equilibrium condition of such as structure), to obtain the solution of problem.Because actually asking
Topic is replaced by better simply problem, so this solution not instead of Exact Solutions, approximate solution.Since most of practical problems are difficult to
Exact Solutions are obtained, and not only computational accuracy is high for finite element, but also adapts to various complicated shapes, thus become effective work
Journey analysis means.
The present invention by finite element method, to the vault 1 of spherical shell with ribbing by larger asymmetric half across, concentrate and carry
Failure mode under lotus effect is analyzed, and according to the special failure mode, has been obtained total under larger concentrfated load effect
Body strain, local train, overall stress, overall stability and local stability obtain the overall collapse mould under geometrical nonlinear analysis
State, according to the mode, it is determined that local optimum scheme, and repeatedly can check and obtain rapidly that optimized reality is non-thread
Property invariant feature value, it is whether qualified under corresponding assessment criteria according to the value, to determine its safety.
Method advantage provided by the present invention is: first, it is combined using overall stability and local stability, finally
Geometrical nonlinear stability is analyzed, and controls unstability characteristic coefficient, analyzes the reality for locally applying concentrated force on spherical shell with ribbing
Border mode and deformation distribution, achieve the anti-unstability side to large-scale ribbed spherical shell in concentrfated load and uniform load while under acting on
Method.Second, this method gives local ribbing simultaneously enhances the Optimized Measures of anti-unstability ability, and passes through comparison Different Optimization
Influence of the measure to overall stability finally to obtain the optimization design for meeting bearing capacity.Third, this method consider linearly
Influence under non-uniformly distributed wind load effect simultaneously to 1 stability of vault simulates the vaults 1 under a variety of non-uniformly distributed wind pressure several
The coping behavior of what non-linear behavior, and by 1 wind resistance behavioral study of vault, it greatly strengthens common in Various Complex wind pressure
Structural stability under effect, makes structure have broad sense compliance characteristics to the oil gas field ground surface environment of all kinds of geologic features.
Method involved in the present invention yet forms a whole set of prioritization scheme to result.The program is according to practical geometry
The specific location that nonlinear stability occurs carries out specific aim reinforcement, gradually extracts nonlinear deformation maximum mesh, increases part
Radial direction carries out local strengthening to angular single rib or cross rib.The design scheme is realized sets about from most severe strain place, fills up
1 overall stiffness notch of vault, because caused by concentrfated load rigidity wooden pail effect be iterated analysis.
Such as the structural schematic diagram that Fig. 8 is Spherical Tank Shell Plates with ribbing, the traditional analysis of 1 stability of vault of spherical shell with ribbing is such as
Under:
The external load allowable of spherical shell with ribbing should be calculated as follows:
In formula: [P] --- the external load allowable (kPa) of spherical shell with ribbing；
E --- the elasticity modulus (MPa) of steel under design temperature；
Rs--- the radius of curvature (m) of spherical shell；
th--- the effective thickness (mm) of tank deck plate；
tm--- the converting thickness (mm) of spherical shell with ribbing.
The converting thickness of spherical shell with ribbing should be calculated according to the following formula.
In formula: t1m--- the converting thickness (mm) of latitudinal ribs and top plate compound section；
h1--- broadwise rib width (mm)；
b1--- latitudinal ribs effective thickness (mm)；
L1S--- latitudinal ribs through to spacing (mm)；
n1--- latitudinal ribs and top plate through to area conversion factor；
e1--- latitudinal ribs and top plate through to compound section centroid into top plate at a distance from face (mm)；
t2m--- through the converting thickness (mm) to rib and top plate compound section；
h2--- through to rib width (mm)；
b2--- through to rib effective thickness (mm)；
L2S--- through to rib broadwise spacing (mm)；
n2--- through to rib and top plate broadwise area conversion factor；
e2--- through to rib and top plate broadwise compound section centroid into top plate at a distance from face (mm).
In traditional design method to Intensity Design to stability by uniformly distributed design external loads lotus and external applied load allowable into
Row comparison obtains qualified criterion, and allowable load is according to limit collapsing load value.It is mainly calculated in the calculation formula of this method
The thickness of uniform load lower roof plate 3 and rib 4, the bending strength of shell thickness do not consider that the structure under concentrfated load is anti-
Curved performance.In fact, after carrying a certain number of concentrfated loads in structure, the anti-stability and integrated stress of structure entirety
Distribution situation is changed, and can be significantly changed in the instability modes of stability, original to occur geometrical non-linearity first
The system point of behavior can change.New maximum stress point also occurs due to the effect of concentrfated load in netted girder structure
Uneven distribution, assesses this geometrical non-linearity behavior and uneven distribution stress level has exceeded the design of existing technologies scheme
Scope, therefore can not be checked suitable for the stability and stress intensity of spherical shell vault 1 with ribbing under the conditions of concentrfated load.
Traditional design scheme is corresponding corresponding uniform load design, but concentrfated load can not be checked, especially
It is stress in the case where currently needing to concentrfated load, linear and nonlinear stability.Lack safety guarantee.Due to acting on arch
It pushes up concentrfated load on 1, all kinds of dead load parts and uniformly design external pressure three summation is larger.Therefore full payload is far more than self-supporting
The bearing capacity of vault 1.If easily causing 1 local buckling of vault or even overall collapse using conventional design, causing to damage.
In comparison, traditional design method can not analyze concentrfated load, can not assess geometrical non-linearity unstability,
It can not obtain local stability strengthened scheme.The method of the present invention realizes the geometrical nonlinear analysis under concentrfated load, to concentration
The analysis of load has comparable accuracy, and position of playing a game moves past general goal and carries out local strengthening, so that analysis result safely may be used
It leans on.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all in spirit of the invention and
Within principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (7)
1. a kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load, which comprises the following steps:
S1: analysis finite element geometrical model is established；
S2: gradually inputting concentrfated load on the vault to simulation on the basis of the analysis finite element geometrical model, described in acquisition
The linear instability modes of vault；
S3: the initial load value of the vault nonlinear stability analysis is obtained according to the linear instability modes；
S4: according to the initial load value, inputting the concentrfated load on the vault, obtains first Linear Stability system
Number, and the linear stable property coefficient is compared with assessed value；
S5: when the linear stable property coefficient is less than institute's rating value, local enhancement component is added, repeating said steps S4 is straight
It is greater than or equal to institute's rating value to the linear stable property coefficient and is considered as qualification；
S6: after the linear stable property coefficient is qualified, model when according to qualification will be described in original input on the vault
Concentrfated load equal proportion amplification setting multiplying power carries out geometrical nonlinear analysis, obtains the collapsing load value of the vault；
S7: the coefficient of stability of the vault is calculated according to the initial load value and the collapsing load value.
2. the spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load according to claim 1, feature exist
In, the step S3 the following steps are included:
S31: according to the xyz three direction displacement for inputting the corresponding point of load of the concentrfated load, the position of each point of load is calculated
Vector value is moved, and multiple displacement vector values are ranked up and obtain its maximum value, finds the maximum displacement vector values pair
The strain maximum point answered；
S32: the maximum displacement vector values are converted after obtaining conversion factor according to total span of l/300, according to institute
State institute rheme of the conversion factor to cell nodes whole in the analysis finite element geometrical model under the linear instability modes
It moves vector value to be converted, the shift value after obtaining all cell node conversions is as the initial load value.
3. the spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load according to claim 1, feature exist
In, the step S6 the following steps are included:
S61: first nonlinear stability property coefficient is obtained according to geometrical nonlinear analysis, and by the nonlinear stability property coefficient
It is compared with non-linear setting characteristic value, judges whether the nonlinear stability property coefficient restrains, the nonlinear stability property coefficient
Show its convergence when more than or equal to the non-linear setting characteristic value；
S62: when the nonlinear stability property coefficient is less than the non-linear setting characteristic value, to the nonlinear stability system
The corresponding point of load of number adds the local enhancement component and is reinforced, and repeating said steps S61 obtains described non-linear
It is described to obtain this until the nonlinear stability property coefficient is greater than or equal to the non-linear setting characteristic value for stability coefficient
The corresponding collapsing load value of nonlinear stability property coefficient.
4. the spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load according to claim 3, feature exist
In: the non-linear characteristic value that sets in the step S61 is 3.0-3.5.
5. the spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load according to claim 1, feature exist
In: the coefficient of stability in the step S7 is the collapsing load value/initial load value.
6. the spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load according to claim 1, feature exist
In: institute's rating value in the step S4 is 4.2.
7. the spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load according to claim 1, feature exist
In: the local enhancement component in the step S5 is steel rib.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910364967.0A CN110096808A (en) | 2019-04-30 | 2019-04-30 | A kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910364967.0A CN110096808A (en) | 2019-04-30 | 2019-04-30 | A kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110096808A true CN110096808A (en) | 2019-08-06 |
Family
ID=67446651
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910364967.0A Pending CN110096808A (en) | 2019-04-30 | 2019-04-30 | A kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110096808A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112784460A (en) * | 2021-01-28 | 2021-05-11 | 北京理工大学重庆创新中心 | Method for analyzing stability of mechanical metamaterial compression bar |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104765908A (en) * | 2015-03-12 | 2015-07-08 | 南车青岛四方机车车辆股份有限公司 | High-speed train body local structure instability finite element simulation device and method |
| CN105488280A (en) * | 2015-11-30 | 2016-04-13 | 上海宇航系统工程研究所 | Method for analyzing composite structure with complex configuration of lunar explorer |
| CN107169209A (en) * | 2017-05-19 | 2017-09-15 | 上海宇航系统工程研究所 | A kind of structure design and optimization method biased under concentrated force effect |
| CN108153981A (en) * | 2017-12-26 | 2018-06-12 | 中航沈飞民用飞机有限责任公司 | A kind of composite material fuselage Material Stiffened Panel Post-Buckling Analysis of Structures method based on finite element analysis |
| CN109684693A (en) * | 2018-12-12 | 2019-04-26 | 中国航空工业集团公司西安飞机设计研究所 | A method of based on the estimated Material Stiffened Panel post-buckling of finite element analysis |
-
2019
- 2019-04-30 CN CN201910364967.0A patent/CN110096808A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104765908A (en) * | 2015-03-12 | 2015-07-08 | 南车青岛四方机车车辆股份有限公司 | High-speed train body local structure instability finite element simulation device and method |
| CN105488280A (en) * | 2015-11-30 | 2016-04-13 | 上海宇航系统工程研究所 | Method for analyzing composite structure with complex configuration of lunar explorer |
| CN107169209A (en) * | 2017-05-19 | 2017-09-15 | 上海宇航系统工程研究所 | A kind of structure design and optimization method biased under concentrated force effect |
| CN108153981A (en) * | 2017-12-26 | 2018-06-12 | 中航沈飞民用飞机有限责任公司 | A kind of composite material fuselage Material Stiffened Panel Post-Buckling Analysis of Structures method based on finite element analysis |
| CN109684693A (en) * | 2018-12-12 | 2019-04-26 | 中国航空工业集团公司西安飞机设计研究所 | A method of based on the estimated Material Stiffened Panel post-buckling of finite element analysis |
Non-Patent Citations (3)
| Title |
|---|
| 尹晔昕: "大型拱顶储罐的有限元计算", 《油气储运》, 31 January 2003 (2003-01-31), pages 23 - 26 * |
| 朱明亮: "《弦支叉筒网壳结构体系研究》", 东南大学出版社 * |
| 西北地区五省一市力学学会: "《压力容器的安全与强度计算》", 天津科学技术出版社, pages: 177 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112784460A (en) * | 2021-01-28 | 2021-05-11 | 北京理工大学重庆创新中心 | Method for analyzing stability of mechanical metamaterial compression bar |
| CN112784460B (en) * | 2021-01-28 | 2023-05-16 | 北京理工大学重庆创新中心 | Stability analysis method for mechanical metamaterial compression bar |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Montuori et al. | Design criteria for diagrid tall buildings: Stiffness versus strength | |
| Enevoldsen et al. | Reliability-based optimization in structural engineering | |
| US20150019177A1 (en) | Method of Determining Prestressing Force of Cable Dome Based on Whole Process Analysis of Cable Dome Tensioning and Bearing | |
| Wadee et al. | Comparative studies of localized buckling in sandwich struts with different core bending models | |
| CN109299540A (en) | Plane frame structure Static Non-linear Analyisis based on rigid body criterion | |
| Mashhadiali et al. | Dynamic increase factor for investigation of progressive collapse potential in tall tube-type buildings | |
| CN106384193A (en) | ICS information safety assessment method based on analytic hierarchy method | |
| Snijder et al. | Lateral torsional buckling design imperfections for use in non‐linear FEA | |
| CN110096808A (en) | A kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load | |
| Shin et al. | Optimum retrofit strategy of FRP column jacketing system for non-ductile RC building frames using artificial neural network and genetic algorithm hybrid approach | |
| CN103924663B (en) | A kind of method for designing improving RC progressive collapse-resisting frame structure ability | |
| Meng et al. | A novel study of structural reliability analysis and optimization for super parametric convex model | |
| Ferreira et al. | Comparison of the General Method with the Overall Method for the out-of-plane stability of members with lateral restraints | |
| Pandey et al. | Modelling and design of cold-formed S960 steel brace-rotated tubular T-and X-joints | |
| CN112380606B (en) | Method for evaluating importance of components and vulnerability of system in civil structure system | |
| CN107066666A (en) | A kind of water supply network cascade robustness measuring method of dynamic emergency recovery mechanism | |
| CN106759373B (en) | A kind of calculation method of long-short pile combined support structure | |
| CN104008258A (en) | Steel structure fire disaster temperature field inverse presumption method based on displacement forms | |
| Al-Aukaily et al. | Sensitivity analysis for displacement-controlled finite-element analyses | |
| CN116070510A (en) | Slope foundation pile horizontal bearing reliability calculation method based on active learning Kriging model | |
| CN111931281B (en) | Method for searching critical failure path of gravity dam-foundation structure system | |
| Savsani | Implementation of modified artificial bee colony (ABC) optimization technique for minimum cost design of welded structures | |
| CN113343424A (en) | Method for calculating ultimate strength of thin-walled beam under action of combined load of bending, shearing and twisting | |
| Jafarian et al. | Tying resistance of reverse channel connection to concrete filled square and rectangular tubular sections | |
| CN110765538A (en) | Improved generalized α method for nonlinear dynamical analysis of complex structure |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20210223 Address after: 100120 No.6 liupukang street, Xicheng District, Beijing Applicant after: CHINA PETROLEUM ENGINEERING & CONSTRUCTION Corp. Applicant after: CHINA NATIONAL PETROLEUM Corp. Applicant after: BEIJING DWELL PETROLEUM & GAS TECHNOLOGY DEVELOPMENT Co.,Ltd. Applicant after: CNPC ENGINEERING DESIGN Co.,Ltd. Address before: 100085 room 301-303, 3 / F, block B, No.22, Shangdi Information Road, Haidian District, Beijing Applicant before: BEIJING DESIGN BRANCH OF CHINA PETROLEUM ENGINEERING CONSTRUCTION Co.,Ltd. |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190806 |