CN105912753B - Seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method - Google Patents

Seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method Download PDF

Info

Publication number
CN105912753B
CN105912753B CN201610206829.6A CN201610206829A CN105912753B CN 105912753 B CN105912753 B CN 105912753B CN 201610206829 A CN201610206829 A CN 201610206829A CN 105912753 B CN105912753 B CN 105912753B
Authority
CN
China
Prior art keywords
abaqus
file
seabed
dimensional
analyzed
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610206829.6A
Other languages
Chinese (zh)
Other versions
CN105912753A (en
Inventor
来向华
张恒
胡涛骏
黄潘阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Second Institute of Oceanography SOA
Original Assignee
Second Institute of Oceanography SOA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Second Institute of Oceanography SOA filed Critical Second Institute of Oceanography SOA
Priority to CN201610206829.6A priority Critical patent/CN105912753B/en
Publication of CN105912753A publication Critical patent/CN105912753A/en
Application granted granted Critical
Publication of CN105912753B publication Critical patent/CN105912753B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design 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 seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method.This method includes that S1. selectes region to be analyzed;S2. it treats analyzed area and carries out pre-treatment;S3. by derived excel file, ABAQUS is generated by ABAQUS numerical value interface routine and imports file .inp;S4. self-weight stress field is constructed;S5. stability analysis is carried out;S6. logarithm calculated result is analyzed.The present invention by ABAQUS numerical value interface routine, enable not only model import in ABAQUS, and can correctly carry out operation and obtain corresponding result of study.

Description

Seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method
Technical field
The present invention relates to a kind of Method for Slope Stability Analysis, especially a kind of seabed side slope three based on Strength Reduction Method Stability maintenance method for qualitative analysis.
Background technique
With China coastal seas sea area engineering development, seabed potential geological disaster especially seabed slope instability has become prestige An important factor for coercing offshore engineering.Submarine landslide, which refers to, not to be consolidated to obtain soft-sediment on submarine slope or has weak structural face Rock, there is a phenomenon where slide for the weak structural face in slope under external force.Seabed Method for Slope Stability Analysis Research be still within developing stage at present.With the generation and development of finite element, using point of stringent ess-strain method Possibility is stablized and be deformed into analysis structure, is gradually applied in geotechnical engineering slope stability analysis.
Most of ABAQUS finite element stability of slope Journal of Sex Research is all based on traditional two dimension and three-dimensional finite element at present It is that model obtains as a result, be no longer satisfied the needs of research and engineering precision, need one kind that can embody submarine geomorphy Finite element modeling calculate work to complete the numerical value of seabed stability of slope early period.
GOCAD (Geological Object Computer Aided Design) Geologic modeling software is French Nancy The Visualization Modeling software for being mainly used in geology field of university's exploitation.GOCAD software have powerful three-dimensional modeling, It visualizes, the function of GEOLOGICAL INTERPRETATION and analysis.Both surface modeling can be carried out, solid modelling can be carried out;Both it can design Geometric objects, can also be with expressive space property distribution.
ABAQUS is the finite element software of a set of powerful engineering simulation, and the range solved the problems, such as is from relatively easy Linear analysis to many complexity nonlinear problem.ABAQUS includes abundant, analog random geometry a list First library.ABAQUS is in addition to can solve a large amount of structures (stress/displacement)) topic, many problems in Other Engineering field can also be simulated, Such as heat transfer, mass diffusion, crack preventing, acoustic analysis, rock-soil mechanics analysis (fluid permeability/stress coupling analysis) And piezoelectric dielectric analysis.
Summary of the invention
Technical problem to be solved by the invention is to provide one kind can satisfy engineering precision need, be able to achieve GOCAD and The seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method of interface between ABAQUS.
It is as follows that the present invention solves technical solution used by above-mentioned technical problem:
Seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method, includes the following steps:
S1. region to be analyzed is selected;
S2. it treats analyzed area and carries out pre-treatment:
S21. the data for being analysed to region, which import, forms three-dimensional landform figure in GOCAD;
S22. the three-dimensional landform figure of formation is imported in Surfer, and forms seabed exterior view after interpolation processing;
S23. seabed exterior view is imported in GOCAD, by sequentially forming grid by point-line-face-body, composition is not advised Three-dimensional bottom relief map under the conditions of upper surface then;
S24. the three-dimensional bottom relief map in GOCAD is exported as into excel file;
S3. by derived excel file, ABAQUS is generated by ABAQUS numerical value interface routine and imports file .inp, number Value interface routine comprises the following processes:
S31. excel file is imported into numerical value interface routine;
S32. to the three-dimensional coordinate parameter of each node of generation unit grid in excel file, the i.e. three-dimensional of p0 to p7 Coordinate parameters are modified in the way of following table:
S33. ABAQUS is exported as from numerical value interface routine imports file .inp;
S4. self-weight stress field is constructed:
S41. it writes and imports file .inp, the direction x that boundary condition is limitation front and back two sides is displaced, the y on left and right two sides Direction is unique, and bottom surface limits the displacement on its two sides xy;
S42. after the completion of writing, it will import in file .inp importing ABAQUS and generate model;
S43. the corresponding module in ABAQUS is executed, by six components of stress of the unit center point in self-weight stress field F11, F22, F33, F12, F13, F23 are output to external file .txt and export;
S44. the file in external file .txt other than data line is deleted, and increases model before the data of each unit Name, format are " -1. serial number of model name ", and each be about to six components of stress belonging to the same unit and separated with comma;
S5. stability analysis is carried out:
S51. the gravitational field file after conversion is led back among former ABAQUS model, surface (table is generated at the top of model Face);
S52. defined analysis step-length, and reduction operation is carried out to cohesive strength and internal friction angle, and be input to input ABAQUS Among program, formula is as follows:
Wherein, C,Respectively cohesive strength and internal friction angle, CmCohesive strength and internal friction angle after being respectively reduced, Fm For strength reduction factor, numberical range 0.5-10;
S53. in ABAQUS, analysis seaway load is at t=0, t=0.25T, t=0.5T and t=0.75T moment, respectively Apply the wave pressure of vertical direction different moments on sea bed surface, according to small amplitude wave theory, the dynamic pressure of water body can be by formula (3) it indicates:
As z=-h, seabed surface pressure are as follows:
Wherein, z is corrugated relative altitude, and H is wave height, and σ is fluctuation circular frequency, and ρ is fluid density, and k is wave number, and h is quiet The water surface to seabed vertical range;
S54. model is integrally calculated using ABAQUS, output variable is safety coefficient FsWith sea bottom surface selected point Displacement;
S6. logarithm calculated result is analyzed:
S61. reduction integral unit displacement isogram analysis: according in t=0, t=0.25T, t=0.5T and t= When tetra- moment of 0.75T calculate termination, the reduction integral unit displacement isogram being respectively formed in ABAQUS, if occurring from more A sliding is developed towards single circle sliding surface, then slides the trend separated towards multiple sliding surfaces by single, then carries out S62;
S62. seabed slope stability is judged by safety coefficient:
S621. by safety coefficient FSDisplacement with sea bottom surface selected point is provided as output variable using ABAQUS Combine function obtains the safety coefficient F of different momentsSWith the variation relation figure of displacement;
S622. pass through safety coefficient FSThe corresponding relationship at value of safety factor value and moment is obtained with the variation relation figure of displacement, If the wave trend of dipping and heaving is presented with variation constantly for value of safety factor value, the minimum value of value of safety factor value is taken, and utilizes " building Slope project technical specification " (GB 50330-2013) judge whether side slope is stable.
Preferably, the data in region to be analyzed are with cad Software Create and with the preservation of .dxf format wait divide in S21 Analyse the bottom contour bathymetric data in region.
Preferably, analysis step-length is at least a cycle T, and cycle T is wave cycle in S53, it is bent to complete once sinusoidal The time of line.
Preferably, in S53, the calculation method of wave number k are as follows:
S531. in MATLAB environment, if wave number is unknown number x, calculation formula are as follows:
F=inline (' x- (/ (g*tanh (x*h))) ')
S532. the range for taking f is -100 to 100;
S533. zero point of the f between -20 20 is searched;
S534. after successive ignition, the value of wave number k is obtained.
Matlab function in the present invention: fzero, for searching the zero point of unitary continuous function;Ezplot, picture letter Number, it is not necessarily to data preparation, directly draws functional digraph;Inline is used to define the Inline Function of a class.
DXF in the present invention is autodesk, inc.'s exploitation for carrying out CAD data between AutoCAD and other softwares The CAD data file format of exchange.DXF is a kind of vector data form of opening.
The present invention compared with the existing technology has the following advantages that and effect:
1, since the present invention is by ABAQUS numerical value interface routine, enable not only model import in ABAQUS, but also And it can correctly carry out operation and obtain corresponding result of study.
2, it is connect due to the present invention by the application between Visual C# language platform development GOCAD and ABAQUS software Mouth program, realizing only need to be using CAD bathymetric chart, complicated seabed can be generated in finite element analysis software ABAQUS Shape and landform model.
3, it since the present invention is by ABAQUS output file conversion programs of Visual C# language platform development, solves Gravity field data format is converted under the conditions of big data quantity, and the field data file after this document is converted complies fully with leads back to Next step static(al) is carried out in ABAQUS software and dynamic load applies the requirement calculated.
4, since the present invention forms unique seabed side slope three-dimensional method for analyzing stability and technology path, utilization are limited Meta analysis means have combined influence of the Wave power load as external loads to seabed side slope stress state, improve The precision of calculating simplifies the process of seabed stability analysis.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with It obtains other drawings based on these drawings.
Fig. 1 is technology path schematic diagram of the invention.
Fig. 2 is model conversion control technology route schematic diagram of the invention.
Fig. 3 is value of safety factor value of the invention with changing schematic diagram constantly.
Fig. 4 is each node location diagram of unit grid of the invention.
Specific embodiment
The present invention will be further described in detail below with reference to the embodiments, following embodiment be explanation of the invention and The invention is not limited to following embodiments.
Embodiment 1:
Since ABAQUS procedural model construction ability is not strong, it is even more unable to do what one wishes especially to handle complicated seabed surface. GOCAD program can satisfy this requirement.GOCAD can generate landform face quickly through CAD, if will by Surfer program After landform carries out kriging analysis processing with finite difference calculus, three-dimensional landform will be more smooth.GOCAD grid generated with It mesh shape required by ABAQUS consistent (hexahedral mesh) and can be exported with Excel file, this is just that the realization of interface is created It may.
Based on above-mentioned discussion, stability analysis technology path are as follows: 1. by original CAD contour map import in GOCAD into Row pre-treatment forms basic contour point;2. will generation contour import Surfer in using kriging analysis method into Row data interpolating, then lead back to formation seabed surface in GOCAD program;3. treated, contour carries out three-dimensional geological modeling, Export modeling grid data;4. grid data, which is imported application programming interfaces, carries out data processing, and generates new ABAQUS meter Input file .inp is calculated, imports in ABAQUS and carries out operation.Its technology path schematic diagram is as shown in Figure 1.
In addition, existing research is achieved the interface between GOCAD and ABAQUS, but due to its grid interpretation sequence It goes wrong, leads to that model can be built out, but be unable to operation, program reports an error.The present invention is based on Visual C# language platform, Realize the grid data file that need to only import the derived Excel format after two program processing of GOCAD and Surfer, it is automatic to solve Data are translated, the generation of the new grid of ABAQUS and subsidiary boundary condition can be completed.Not only model can import the present invention In ABAQUS, and it simultaneously can correctly carry out operation and obtain corresponding result of study.The model conversion of the process controls skill Art route map is as shown in Figure 2.
The seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method of the present embodiment, includes the following steps:
S1. region to be analyzed is selected, in the present embodiment, the long 1610m in region, wide 890m, wavelength 36.5m to be analyzed, week Phase T=5s, wave height H=4m;
S2. it treats analyzed area and carries out pre-treatment:
S21. the data for being analysed to region, which import, forms three-dimensional landform figure in GOCAD;
S22. the three-dimensional landform figure of formation is imported in Surfer, and forms seabed exterior view after interpolation processing;
S23. seabed exterior view is imported in GOCAD, by being sequentially formed grid by point-line-face-body and being formed irregularly Upper surface under the conditions of three-dimensional bottom relief map;
S24. the three-dimensional bottom relief map in GOCAD is exported as into excel file;
S3. by derived excel file, ABAQUS is generated by ABAQUS numerical value interface routine and imports file .inp, number Value interface routine comprises the following processes:
S31. excel file is imported into numerical value interface routine;
S32. to the three-dimensional coordinate parameter of each node of generation unit grid in excel file, i.e. p0 to p7 (such as Fig. 4 It is shown), it modifies in the way of following table:
S33. ABAQUS is exported as from numerical value interface routine imports file .inp;
S4. self-weight stress field is constructed:
S41. it writes and imports file .inp, the direction x that boundary condition is limitation front and back two sides is displaced, the y on left and right two sides Direction is unique;Bottom surface limits the displacement on its two sides xy;
S42. after the completion of writing, it will import in file .inp importing ABAQUS and generate model;
S43. the corresponding module in ABAQUS is executed, by six components of stress of the unit center point in self-weight stress field F11, F22, F33, F12, F13, F23 are output to external file .txt and export;
S44. the file in external file .txt other than data line is deleted, and increases model before the data of each unit Name, format are " -1. serial number of model name ", and each be about to six components of stress belonging to the same unit and separated with comma;
S5. stability analysis is carried out:
S51. the gravitational field file after conversion is led back among former ABAQUS model, surface is generated at the top of model;
S52. defined analysis step-length, and reduction operation is carried out to cohesive strength and internal friction angle, and be input to input ABAQUS Among program, formula is as follows:
Wherein, C,Respectively cohesive strength and internal friction angle, Cm,Cohesive strength and internal friction angle after being respectively reduced, Fm is strength reduction factor, numberical range 0.5-10;
S53. in ABAQUS, analysis seaway load is at t=0, t=0.25T, t=0.5T and t=0.75T moment, respectively Apply the wave pressure of vertical direction different moments on sea bed surface, according to small amplitude wave theory, the dynamic pressure of water body can be by formula (3) it indicates:
As z=-h, seabed surface pressure are as follows:
Wherein, z is corrugated relative altitude, and H is wave height, and σ is fluctuation circular frequency, and ρ is fluid density, and k is wave number, and h is quiet The water surface to seabed vertical range;
S54. model is integrally calculated using ABAQUS, output variable is safety coefficient FSWith sea bottom surface selected point Displacement;
S6. logarithm calculated result is analyzed:
S61. reduction integral unit displacement isogram analysis: according in t=0, t=0.25T, t=0.5T and t= When tetra- moment of 0.75T calculate termination, the reduction integral unit displacement isogram being respectively formed in ABAQUS, if occurring from more A sliding is developed towards single circle sliding surface, then slides the trend separated towards multiple sliding surfaces by single, then carries out S62;
S62. seabed slope stability is judged by safety coefficient:
S621. by safety coefficient FSDisplacement with sea bottom surface selected point is provided as output variable using ABAQUS Combine function obtains the safety coefficient F of different momentsSWith the variation relation figure of displacement;
S622. pass through safety coefficient FSThe corresponding relationship at value of safety factor value and moment is obtained with the variation relation figure of displacement, If the wave trend of dipping and heaving is presented with variation constantly for value of safety factor value, the minimum value of value of safety factor value is taken, and utilizes " building Slope project technical specification " (GB 50330-2013) judge whether side slope is stable.
In above-mentioned S21, the data in region to be analyzed are with cad Software Create and with the region to be analyzed of .dxf format preservation Bottom contour bathymetric data.
In above-mentioned S53, analysis step-length is at least a cycle T, and cycle T is wave cycle, completes once sinusoidal curve Time.
In above-mentioned S53, the calculation method of wave number k are as follows:
S531. in MATLAB environment, if wave number is unknown number x, calculation formula are as follows:
F=inline (' x- (/ (g*tanh (x*h))) ')
S532. the range for taking f is -100 to 100;
S533. zero point of the f between -20 to 20 is searched;
S534. after successive ignition, the value of wave number k is obtained.
The calculation method of wave number k specifically:
In above-mentioned S43, in order to extract the field data of self-weight stress field, need to initially enter Visualization Post-processing module opens corresponding calculated result file.Again by six components of stress of the unit center point in self-weight stress field It is output in external file .txt and exports, just can be carried out in the model being then imported with field data under acting by external hydrostatic pressure into one The reduction of step is analyzed.But derived field data cannot be directly read by ABAQUS/CAE, it is necessary to by passing through after processing ABAQUS/command module is read.
In above-mentioned S44, it will increase before the data of file and each unit other than deletion data line in external file .txt Add model name, format is " -1. serial number of model name ", and each stress number for being about to different stress directions belonging to the same unit According to that is, six components of stress F11, F12, F13, F12, F13, F23 are separated with comma, such as former export data are as follows:
It is translated the file by field file converter are as follows:
Model name -1.1, F11, F12, F13, F12, F13, F23
Model name -1.2, F11, F12, F13, F12, F13, F23
Model name -1.3, F11, F12, F13, F12, F13, F23
In above-mentioned S622, " Technique Code for Building Slope Engineering " (GB 50330-2013) regulation: 5.3.1 is except check operating condition Outside, slope stability state is divided into stable, basicly stable, understable and unstable four kinds of states, can be according to slope stability system Number is determined by table 5.3.1.
Table 5.3.1 slope stability state demarcation
Note: Fst--- safety factor of slope.
5.3.2 safety factor of slope Fst should be determined by table 5.3.2, when stability factor of slope is less than stability of slope Side slope is answered to be handled when safety coefficient.
Table 5.3.3 safety factor of slope Fm
Note: when 1. earthquake operating condition, safety coefficient is only applicable to the side slope without important buildings or structures in slumping area;
2. pair geological conditions is very responsible or the CR Critical slope project of consequence of failure, buckling safety factor should be mentioned suitably It is high.
Above-mentioned FmFor strength reduction factor, numberical range is chosen between 0.5-10, in the present embodiment using
In addition, it should be noted that, the specific embodiments described in this specification, the shape of parts and components are named Title etc. can be different.The equivalent or simple change that all structure, feature and principles described according to the invention patent design are done, is wrapped It includes in the scope of protection of the patent of the present invention.Those skilled in the art can be to described specific implementation Example is done various modifications or additions or is substituted in a similar manner, and without departing from structure of the invention or surmounts this Range as defined in the claims, is within the scope of protection of the invention.

Claims (4)

1. a kind of seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method, it is characterized in that: including the following steps:
S1. region to be analyzed is selected;
S2. it treats analyzed area and carries out pre-treatment:
S21. the data for being analysed to region, which import, forms three-dimensional landform figure in GOCAD;
S22. the three-dimensional landform figure of formation is imported in Surfer, and forms seabed exterior view after interpolation processing;
S23. seabed exterior view is imported in GOCAD, by sequentially forming grid by point-line-face-body, is formed irregular Three-dimensional bottom relief map under the conditions of upper surface;
S24. the three-dimensional bottom relief map in GOCAD is exported as into excel file;
S3. by derived excel file, ABAQUS is generated by ABAQUS numerical value interface routine and imports file .inp, numerical value connects Mouth program comprises the following processes:
S31. excel file is imported into numerical value interface routine;
S32. to the three-dimensional coordinate parameter of each node of generation unit grid in excel file, the i.e. three-dimensional coordinate of p0 to p7 Parameter is modified in the way of following table:
S33. ABAQUS is exported as from numerical value interface routine imports file .inp;
S4. self-weight stress field is constructed:
S41. it writes and imports file .inp, the direction x that boundary condition is limitation front and back two sides is displaced, the direction y on left and right two sides Uniquely, bottom surface limits the displacement on its two sides xy;
S42. after the completion of writing, it will import in file .inp importing ABAQUS and generate model;
S43. execute ABAQUS in corresponding module, by six components of stress F11 of the unit center point in self-weight stress field, F22, F33, F12, F13, F23 are output to external file .txt and export;
S44. the file in the external file .txt other than data line is deleted, and increases model before the data of each unit Name, format are " -1. serial number of model name ", and each be about to six components of stress belonging to the same unit and separated with comma;
S5. stability analysis is carried out:
S51. the gravitational field file after conversion is led back among former ABAQUS model, surface is generated at the top of model;
S52. defined analysis step-length, and reduction operation is carried out to cohesive strength and internal friction angle, and be input to input ABAQUS program Among, formula is as follows:
Wherein, C,Respectively cohesive strength and internal friction angle, Cm,Cohesive strength and internal friction angle after being respectively reduced, FmFor Strength reduction factor, numberical range 0.5-10;
S53. in ABAQUS, seaway load is analyzed at t=0, t=0.25T, t=0.5T and t=0.75T moment, respectively in sea Bed surface applies the wave pressure of vertical direction different moments, and according to small amplitude wave theory, the dynamic pressure of water body can be by formula (3) table Show:
As z=-h, seabed surface pressure are as follows:
Wherein, z is corrugated relative altitude, and H is wave height, and σ is fluctuation circular frequency, and ρ is fluid density, and k is wave number, and h is standing level To the vertical range in seabed, T is time cycle unit, and each period is 5 seconds;
S54. model is integrally calculated using ABAQUS, output variable is the position of safety coefficient Fs and sea bottom surface selected point It moves;
S6. logarithm calculated result is analyzed:
S61. reduction integral unit displacement isogram analysis: according in t=0, t=0.25T, t=0.5T and t=0.75T tetra- When a moment calculates termination, the reduction integral unit displacement isogram being respectively formed in ABAQUS, if occurring from multiple slidings It is developed towards single circle sliding surface, then slides the trend separated towards multiple sliding surfaces by single, then carry out S62, wherein T is the time Unit;
S62. seabed slope stability is judged by safety coefficient:
S621. using the displacement of safety coefficient FS and sea bottom surface selected point as output variable, the Combine provided using ABAQUS Function obtains the safety coefficient FS of different moments with the variation relation figure of displacement;
S622. the corresponding relationship at value of safety factor value and moment is obtained with the variation relation figure of displacement by safety coefficient FS, if peace Complete set numerical value carves the wave trend that dipping and heaving is presented in variation at any time, takes the minimum value of value of safety factor value, and utilize " building slope Engineering legislation " (GB50330-2013) judge whether side slope is stable.
2. the seabed side slope three-dimensional method for analyzing stability according to claim 1 based on Strength Reduction Method, it is characterized in that: In the S21, the data in region to be analyzed are with cad Software Create and with the seabed in the region to be analyzed of .dxf format preservation Isobath bathymetric data.
3. the seabed side slope three-dimensional method for analyzing stability according to claim 1 based on Strength Reduction Method, it is characterized in that: In the S53, analysis step-length is at least a cycle T, the cycle T be wave cycle complete once sinusoidal curve when Between.
4. the seabed side slope three-dimensional method for analyzing stability according to claim 1 based on Strength Reduction Method, it is characterized in that: In the S53, the calculation method of wave number k are as follows:
S531. in MATLAB environment, if wave number is unknown number x, calculation formula are as follows:
F=inline (' x- (/ (g*tanh (x*h))) ')
S532. the range for taking f is -100 to 100;
S533. zero point of the f between -20 to 20 is searched;
S534. after iteration several times, the value of wave number k is obtained.
CN201610206829.6A 2016-04-01 2016-04-01 Seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method Expired - Fee Related CN105912753B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610206829.6A CN105912753B (en) 2016-04-01 2016-04-01 Seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610206829.6A CN105912753B (en) 2016-04-01 2016-04-01 Seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method

Publications (2)

Publication Number Publication Date
CN105912753A CN105912753A (en) 2016-08-31
CN105912753B true CN105912753B (en) 2019-05-14

Family

ID=56745338

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610206829.6A Expired - Fee Related CN105912753B (en) 2016-04-01 2016-04-01 Seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method

Country Status (1)

Country Link
CN (1) CN105912753B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106855637A (en) * 2016-12-10 2017-06-16 国家海洋局第二海洋研究所 Underwater Slope method for analyzing stability
CN106872081A (en) * 2017-03-28 2017-06-20 北京科技大学 Stress monitoring method for early warning based on the natural frequency of vibration
CN107330182B (en) * 2017-06-28 2020-08-14 西北农林科技大学 Method for calculating slope stability safety coefficient by reducing strength based on humidification condition
CN107391886A (en) * 2017-08-29 2017-11-24 广州海洋地质调查局 Submarine geological disaster evaluation method and system
CN108121853B (en) * 2017-11-29 2021-03-26 吴先富 System and method for calculating excavation and filling engineering quantity based on AutoCAD
CN110261573B (en) * 2019-05-16 2021-09-03 同济大学 Dynamic evaluation method for stability of high-position rocky landslide
CN111581836B (en) * 2020-05-14 2023-12-15 国网电力空间技术有限公司 Power transmission line landslide body stability calculation method
CN112257140B (en) * 2020-09-16 2022-10-25 南京工业大学 Safety coefficient calculation method for seabed slope stability
CN114036785B (en) * 2021-09-24 2024-04-05 自然资源部第二海洋研究所 Deep-sea hydrothermal area deep-structure stress field simulation method based on ocean shell speed structure

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003213591A (en) * 2002-01-18 2003-07-30 Dexter Corp Improvement of dry flexural strength of cold injection packaging material
CN105224704A (en) * 2014-06-25 2016-01-06 中国科学院海洋研究所 A kind of evaluation method of seabed erosion

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003213591A (en) * 2002-01-18 2003-07-30 Dexter Corp Improvement of dry flexural strength of cold injection packaging material
CN105224704A (en) * 2014-06-25 2016-01-06 中国科学院海洋研究所 A kind of evaluation method of seabed erosion

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Iterative velocity model building for 3-D depth migration by integrating GOCAD and AVS;Robert G.Clapp等;《Stanford Exploration Project》;20010515;1-9
基于ABAQUS的内河航道岸坡稳定性分析;程曦 等;《水道港口》;20120430;第33卷(第2期);第152-158页

Also Published As

Publication number Publication date
CN105912753A (en) 2016-08-31

Similar Documents

Publication Publication Date Title
CN105912753B (en) Seabed side slope three-dimensional method for analyzing stability based on Strength Reduction Method
Penalba et al. Using NEMOH for modelling wave energy converters: A comparative study with WAMIT
Armesto et al. Comparative analysis of the methods to compute the radiation term in Cummins’ equation
Capizzano Turbulent wall model for immersed boundary methods
Fuhrman et al. Tsunami generation, propagation, and run-up with a high-order Boussinesq model
Bletzinger et al. Efficient modeling in shape optimal design
CN108763625A (en) A kind of method and system of FLAC3D softwares pre-treatment
Maier On some issues in shakedown analysis
Pan et al. Numerical simulation for 2D shallow water equations by using Godunov-type scheme with unstructured mesh
Chan SUMMAC-A numerical model for water waves
Bai et al. Convergence of multilayer nonhydrostatic models in relation to Boussinesq-type equations
He et al. New speedup algorithms for nonlinear dynamic time history analysis of supertall building structures under strong earthquakes
Li et al. ISWFoam: A numerical model for internal solitary wave simulation in continuously stratified fluids
Hojjat et al. Fluid-structure interaction in the context of shape optimization and computational wind engineering
Harris et al. Simulation of wave forces on a gravity based foundation by a BEM based on fully nonlinear potential flow
Shin Simulation of two-dimensional internal waves generated by a translating and pitching foil
Lee et al. A time-dependent radiation condition for transient wave-structure interactions
Han et al. Transient analysis of three‐dimensional dynamic interaction between multilayered soil and rigid foundation
Cavallo et al. A Parallel Adaptation Package for Three-Dimensional Mixed-Element Unstructured Meshes
Chandar Development of a Parallel Overset Grid Framework for Moving Body Simulations in Open FOAM.
Spears et al. Nonlinear Time Domain Seismic Soil-Structure Interaction (SSI) Deep Soil Site Methodology Development
Kvamsdal et al. Fluid-structure interaction for structural design
CN111583400A (en) Three-dimensional geological model modeling method, system, device and storage medium
Rai A new integrated finite volume–finite volume hydrodynamic modelling framework for wave-structure interactions
Wei et al. Turbulence model verification and validation in an open source environment

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190514

Termination date: 20200401

CF01 Termination of patent right due to non-payment of annual fee