CN112711884B - Simulation prediction method for quasi-static damage of brittle material by combination of implicit algorithm - Google Patents
Simulation prediction method for quasi-static damage of brittle material by combination of implicit algorithm Download PDFInfo
- Publication number
- CN112711884B CN112711884B CN202011605539.1A CN202011605539A CN112711884B CN 112711884 B CN112711884 B CN 112711884B CN 202011605539 A CN202011605539 A CN 202011605539A CN 112711884 B CN112711884 B CN 112711884B
- Authority
- CN
- China
- Prior art keywords
- brittle material
- damage
- loading
- quasi
- implicit
- 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.)
- Active
Links
- 239000000463 material Substances 0.000 title claims abstract description 97
- 230000006378 damage Effects 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000004088 simulation Methods 0.000 title claims abstract description 20
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 59
- 238000004364 calculation method Methods 0.000 claims abstract description 53
- 239000011521 glass Substances 0.000 claims description 18
- 238000010586 diagram Methods 0.000 description 10
- 230000000007 visual effect Effects 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000013001 point bending Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 239000005340 laminated glass Substances 0.000 description 1
- 238000012916 structural analysis 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The invention provides a simulation prediction method for quasi-static damage of brittle materials by combining implicit algorithm. The method comprises the following steps: establishing a finite element analysis model, and setting an implicit calculation analysis step; writing a subprogram, and importing the subprogram into an implicit calculation simulation process to obtain the time from the loading start to the damage of the brittle material; stopping a subroutine, adjusting the phenomenon time to be the time from the loading start to the time when the brittle material is about to be damaged, re-performing an implicit calculation simulation process, outputting a load curve of a loading point, determining the damage load of the brittle material under the quasi-static loading condition, storing the result of the implicit calculation of the loading, transmitting the result to an explicit analysis, and performing analysis of the quasi-static loading damage of the brittle material; and predicting the quasi-static damage characteristic of the brittle material according to the analysis result. The method can solve the problem that the implicit analysis method is difficult to converge in the brittle material damage stage, and simultaneously greatly improves the calculation efficiency.
Description
Technical Field
The invention relates to the field of analysis of quasi-static loading damage of brittle materials, in particular to a simulation prediction method of quasi-static damage of brittle materials by combining implicit algorithm.
Background
The quasi-static loading state of the brittle material refers to a state that a certain point or a certain area of the brittle material is subjected to load changing along with time and the loading speed of the load is relatively gentle. The quasi-static properties of brittle materials are primarily determined by the peak load of failure, and this variable is affected by the brittle material structure. At present, in the aspect of quasi-static loading analysis of brittle materials and related structures, no prediction analysis of the quasi-static failure characteristics of the brittle materials exists. The method for analyzing the quasi-static loading damage of the brittle material by using the finite element method is of great significance for further and deeply researching and predicting the quasi-static damage characteristics of the brittle material with high efficiency.
With the gradual maturity of finite element analysis technology, quasi-static loading and damage simulation technology for brittle materials are becoming more mature, and the whole quasi-static loading process of the brittle materials can be analyzed by using an implicit algorithm suitable for acceleration with little change; for the simulation of the fracture phenomenon of the brittle material, a method of inserting a cohesive force model into the brittle material to simulate a crack, namely, when the fracture strength of the cohesive force model, that is, the fracture strength of the brittle material is achieved, cohesive force units are failed and deleted to simulate the crack has been used, and good results are obtained.
However, there are two problems in purely implicit analysis of the quasi-static fracture of a brittle material, one is that after the brittle material reaches the fracture load, the fracture phase of the brittle material, especially glass, is often not more than 1ms (Giulio Castor. Structural analysis of failure behavior of laminated glass [ J ]. Composite Part B.2017), and the computation of the fracture phase is extremely difficult, if not lengthy; meanwhile, in ABAQUS, although the calculation of the whole breaking phase can be converged by setting the contact stability control, since the mechanism of the contact stability control of ABAQUS is to converge the calculation of the breaking phase by adjusting the contact damping of the contact pair, the breaking phase is a nonlinear phase with extremely short time, and the contact stability control automatically enlarges the contact damping at this phase, which leads to serious distortion of the analysis result of the breaking phase.
Meanwhile, if the quasi-static loading damage process of the brittle material is analyzed by using a purely explicit method, the time of the loading stage is often several seconds or even tens of seconds, so that the calculation time is very much, and the calculation efficiency is very low.
The ABAQUS finite element analysis software can realize data transmission between an implicit calculation step and an explicit calculation step, and the explicit calculation is more suitable for analysis with extremely short phenomenon time, so that for analysis of the quasi-static loading damage problem of the brittle material, the loading stage of the brittle material can be analyzed by using an implicit algorithm, and the result of the loading process is transmitted to the explicit calculation to analyze the damage stage of the brittle material; meanwhile, the ABAQUS is provided with various subprogram interfaces, the UVAARM subprogram can extract the stress of the units, the UVAARM subprogram can be used for extracting the stress of the cohesive units of the brittle material, and when the brittle material is about to be damaged by loading, the promoter subprogram stops calculation and performs data transmission work.
Disclosure of Invention
The invention provides a simulation method for analyzing quasi-static loading conditions of brittle materials by using a combination of implicit algorithm, wherein the implicit algorithm is used for analysis in the loading stage analysis of the brittle materials, a subprogram is inserted at the same time, when the brittle materials are about to be damaged, the subprogram is started to stop calculation, the phenomenon time of an implicit analysis step is adjusted to the phenomenon time when the previous calculation is stopped, the subprogram is stopped, and calculation data are transmitted to the explicit calculation step for analysis of the glass damage stage after the loading stage of the brittle materials is completely analyzed.
The object of the invention is achieved by at least one of the following technical solutions.
The simulation prediction method for the quasi-static damage of the brittle material by combining the implicit algorithm comprises the following steps:
s1, establishing a finite element analysis model of brittle material quasi-static loading damage in ABAQUS, setting an implicit calculation analysis step, and associating the ABAQUS, visual Studio 2012 and Intel Fortran XE 2013;
s2, writing an ABAQUS subroutine for judging whether the brittle material is about to be damaged, and introducing the ABAQUS subroutine into an implicit calculation simulation process of quasi-static loading of the brittle material, and starting a command for stopping calculation by the subroutine when the brittle material is about to be damaged, so as to obtain the time from the loading to the about to be damaged of the brittle material;
s3, stopping the ABAQUS subroutine in the step S2, adjusting the phenomenon time in the implicit calculation analysis step set in the step S1 to the time from the loading start to the about-to-break of the brittle material, which is obtained in the step S2, in the ABAQUS, re-performing the implicit calculation simulation process of the quasi-static loading of the brittle material, analyzing the loading condition of the brittle material in the time period, outputting a load curve of a loading point, determining the breaking load of the brittle material under the quasi-static loading condition, and simultaneously storing the result of the loading implicit calculation;
s4, transmitting result data of the implicit calculation to explicit analysis, and analyzing the quasi-static loading damage of the brittle material;
s5, predicting the quasi-static damage characteristic of the brittle material according to the analysis result in the step S4.
In step S1, an ABAQUS software is used for establishing a three-dimensional model of the brittle material quasi-static damage, dividing finite element grids, endowing the three-dimensional model with material properties, related contact and boundary conditions, and completing establishment of a finite element analysis model of the brittle material quasi-static loading damage.
Further, in step S1, the phenomenon time in the implicit calculation analysis step, i.e., the physical time of the loading process, is sufficiently long.
Further, in step S2, the method for determining the time from the start of loading to the point of breaking the brittle material is specifically as follows:
in the implicit calculation process of loading, obtaining the section equivalent tensile stress of a cohesive unit of the brittle material every time a time step is calculated, wherein the mechanism of simulating cracks by the cohesive unit is based on fracture mechanics theory; in fracture mechanics, the damage of the material is type I damage, type II damage and type III damage; the damage of the brittle material is typical tensile fracture, the I-type damage is open type damage, and the damage is the tensile damage caused by the section equivalent tensile stress, so that the section equivalent tensile stress of the glass cohesive force unit obtained by completing the calculation of one time step can be compared with the fracture stress standard set in the ABAQUS subroutine, and when the section equivalent tensile stress of the obtained brittle material cohesive force unit is greater than or equal to the fracture stress standard set by the subroutine, the calculation is stopped, and the time from the loading start to the brittle material to be damaged can be determined.
Further, the final result of the implicit loading calculation stored in step S3 is the final state of the brittle material calculated in the implicit loading stage, and is also the initial state of the analysis of the brittle material quasi-static loading destruction performed in step S4.
Further, in step S4, analysis of brittle material quasi-static loading damage is performed using an explicit algorithm carried by ABAQUS.
Further, in step S5, after the simulation analysis result in step S4 is obtained, the simulation result and the test result are compared and analyzed, so that the quasi-static fracture characteristics of the brittle material can be accurately predicted, and further, the quasi-static fracture characteristics of the brittle material can be further and deeply studied and predicted with high efficiency.
Compared with the prior art, the invention has the advantages that:
the invention uses the explicit algorithm to analyze the brittle material damage stage, not only can solve the convergence problem of the pure implicit analysis brittle material quasi-static loading damage stage, but also greatly improves the calculation efficiency and the calculation precision of the brittle material damage analysis.
Drawings
FIG. 1 is a schematic diagram of a three-point bending model of pure glass according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of the setting of equivalent tensile stress and shear stress in a cross section of a cohesive force model of glass according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of an associated Intel Fortran XE 2013 and Visual Studio 2012 in accordance with an embodiment of the present invention;
FIG. 4 is a schematic diagram of associating ABAQUS 2016 with Intel Fortran XE 2013 in an embodiment of the present invention;
FIG. 5 is a schematic diagram of a variable that needs to be called by a setup subroutine in an embodiment of the present invention;
FIG. 6 is a schematic diagram of a subroutine inserted in creating Job in an embodiment of the present invention;
FIG. 7 is a schematic diagram showing a stop of calculation when glass is about to break in an embodiment of the present invention;
FIG. 8 is a schematic diagram illustrating the setting of delete subroutine call variables in an embodiment of the present invention;
FIG. 9 is a stress cloud diagram of the loading phase ending in an embodiment of the present invention;
FIG. 10 is a graph of a stress cloud at the beginning of a failure phase in an embodiment of the invention;
FIG. 11 is a schematic diagram of a simulated crack in an embodiment of the invention;
FIG. 12 is a graph showing the load-displacement curve according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, a detailed description of the specific implementation of the present invention will be given below with reference to the accompanying drawings and examples.
Examples:
in this embodiment, taking three-point bending of pure glass with the size of 120mm×20mm×3mm under ABAQUS 2016 environment as an example, the simulation prediction method of quasi-static damage of brittle material combined by implicit algorithm comprises the following steps:
step 1: modeling a pure glass three-point bending model.
In this embodiment, as shown in fig. 1, the pure glass three-point bending finite element model is symmetrical about the xy plane, so only 1/2 model is taken to improve the calculation efficiency; the equivalent stress setting of the cohesive force model is shown in fig. 2.
Step 2: ABAQUS 2016, intel Fortran XE 2013, visual Studio 2012 are correlated.
The ABAQUS 2016 needs to be compiled before invoking a subroutine, meanwhile, all the subroutine interfaces of the ABAQUS use Fortran 77 language, and the Visual Studio 2012 is used for writing the subroutine, so that the ABAQUS 2016, the Intel Fortran XE 2013 and the Visual Studio 2012 have to be associated, and the ABAQUS 2012 can be normally compiled when the subroutine is inserted.
The association steps are as follows:
in this embodiment, the system is a 64-bit system, and all programs are opened, namely, intel Parallel Studio XE 2013, command template, parallel Studio XE with Intel Compiler XE v 14.0.0 Update 1, intel 64Visual Studio 2012mode, and after opening, command%com%/k (position of vcvarsal. Bat file) "x86_amd64 is input and run to associate Intel Fortran XE 2013 with Visual Studio 2012, as shown in FIG. 3;
as shown in fig. 4, the ABAQUS 2016 mount directory is edited in the abq2016.Bat and launcher. Bat files, and two lines of commands are added on the original basis:
call "(location of vcvarsall. Bat file)" x86_amd64;
call "(location of ifortvars. Bat file)" intel64 vs2012.
The role of these two lines of commands is to associate ABAQUS with Visual Studio 2012 and Intel Fortran XE 2013, respectively.
Step 3: and writing a UVARISC subroutine.
In this embodiment, since glass is a brittle material and breakage of the glass is in most cases tensile breakage, the first strength theory can be used as a criterion for judging whether or not calculation is stopped, and the method of programming the subroutine is as follows: and writing a frame of the subroutine according to a format specified by a UVARISC subroutine of ABAQUS, writing a command for acquiring the equivalent positive stress of the cross section of the cohesive unit in the x direction in the frame, and judging a statement for stopping calculation when the equivalent positive stress of the cohesive unit is greater than or equal to a failure stress standard.
In this embodiment, since the syntax of the subroutine of ABAQUS is Fortran 77 syntax, the written subroutine is saved as a file for, and ABAQUS can call the subroutine.
Step 4: a phenomenon time long enough to load implicit analysis is added while inserting subroutines.
The range of the subroutine acquisition model unit must be determined before inserting the UVARM subroutine, in this embodiment, the range of the subroutine acquisition model unit is the cohesive force unit of the whole glass, so the setting of "User Output Variable" is added to the cohesive force material property defining portion to determine that the subroutine acquisition model unit is the cohesive force unit, and the variables to be called are set, in this embodiment, since the subroutine only needs to call one variable, the number of variables is set to 1, as shown in fig. 5.
Importing the subroutine file generated in step 3 in the Job module before submitting the calculation example, as shown in fig. 6; the command to stop the calculation when the glass is about to break is initiated as shown in fig. 7.
Step 5: and (3) adjusting the phenomenon time for loading the implicit analysis to the phenomenon time when the calculation is stopped in the step (4), deleting the setting of the subroutine call variable, starting the restarting setting, and calculating again.
In ABAQUS, implicit and explicit computations cannot be seamlessly interfaced, so step 5 is an indispensable step in order to save the computation results of the load phase as an initial state of the destroy phase. The settings of the subroutine call variables need to be deleted before this, as shown in fig. 8.
Step 6: and (5) transmitting the loading implicit analysis calculation data in the step 5 to explicit analysis, and performing glass damage explicit analysis.
The equivalent stress calculation result after the implicit analysis of the loading process, the initial equivalent stress distribution after the explicit analysis of the destruction process, and the simulated crack at the destruction stage in this embodiment are shown in fig. 9, 10, and 11. Comparing fig. 9 and fig. 10, it can be known that the equivalent stress distribution after the implicit analysis of the loading process is finished and the equivalent stress distribution at the initial stage of the explicit analysis of the destruction process are the same, and the error of the numerical value of the stress interval is far less than 1%, which indicates that the calculated data of the implicit analysis of the loading process is successfully transmitted; the simulated failure mode shown in fig. 11 is a crack at the centerline of the glass, and corresponds to a failure mode in which the brittle material is broken under three-point bending loading, and theoretically only a crack is broken at the place where the equivalent stress is maximum.
Load-displacement curves for load-destruction processes for purely implicit and implicit-explicit analysis are shown in fig. 12. In fig. 12, the load-displacement curve of the load-destruction process and the load-displacement curve of the test using implicit-explicit analysis have high consistency, and in particular, the consistency of the destruction process curve is higher than that of the load-displacement curve of the pure implicit analysis, indicating that the implicit-explicit analysis method is higher in calculation accuracy than that of the pure implicit analysis method.
The following table shows the time comparisons of pure implicit method analysis and implicit-explicit combination method analysis of glass quasi-static loading damage:
TABLE 1
It can be seen that the analysis efficiency using implicit-explicit analysis of glass quasi-static loading damage is greatly improved over pure implicit analysis.
Claims (1)
1. The simulation prediction method for the quasi-static damage of the brittle material by combining the implicit algorithm is characterized by comprising the following steps of:
s1, establishing a finite element analysis model of brittle material quasi-static loading damage in ABAQUS, setting a phenomenon time, namely an implicit calculation analysis step with a physical time long enough in a loading process, and associating ABAQUS, visualStudio 2012 with Intel Fortran XE 2013;
s2, writing an ABAQUS subroutine for judging whether the brittle material is about to be damaged, and introducing the ABAQUS subroutine into an implicit calculation simulation process of quasi-static loading of the brittle material, and starting a command for stopping calculation by the ABAQUS subroutine when the brittle material is about to be damaged, so as to obtain the time from the loading to the about to be damaged of the brittle material;
s3, stopping the ABAQUS subroutine in the step S2, adjusting the phenomenon time in the implicit calculation analysis step set in the step S1 to the time from the loading start to the about-to-break of the brittle material obtained in the step S2 in the ABAQUS subroutine, re-performing the implicit calculation simulation process of the quasi-static loading of the brittle material, analyzing the loading condition of the brittle material in the time, outputting a load curve of a loading point, determining the breaking load of the brittle material under the quasi-static loading condition, and simultaneously storing the result of the loading implicit calculation;
s4, transmitting result data of the implicit calculation to explicit analysis, and analyzing the quasi-static loading damage of the brittle material;
s5, predicting the quasi-static damage characteristic of the brittle material according to the analysis result in the step S4;
in step S1, the establishing a finite element analysis model of the brittle material quasi-static loading failure specifically includes: establishing a three-dimensional model of the quasi-static damage of the brittle material by using ABAQUS, dividing a finite element grid, endowing the model with material properties, related contact and boundary conditions, and completing the establishment of a finite element analysis model of the quasi-static loading damage of the brittle material;
in step S2, the method for determining the time from the start of loading to the point of breaking the brittle material is specifically as follows:
in the implicit calculation process of loading, obtaining the section equivalent tensile stress of a cohesive unit of the brittle material every time a time step is calculated, wherein the mechanism of simulating cracks by the cohesive unit is based on fracture mechanics theory; in fracture mechanics, the damage of the material is type I damage, type II damage and type III damage; the damage of the brittle material is typical tensile fracture, the I-type damage is open type damage, and the damage is the tensile damage caused by the section equivalent tensile stress, so that the section equivalent tensile stress of the glass cohesive unit obtained by completing the calculation of one time step can be compared with the fracture stress standard set in the ABAQUS subroutine, and when the section equivalent tensile stress of the obtained brittle material cohesive unit is greater than or equal to the fracture stress standard set by the ABAQUS subroutine, the calculation is stopped, and the time from the loading start to the brittle material to be damaged can be determined;
in the step S3, the final result of the stored implicit loading calculation is the final state of the brittle material calculated in the implicit loading stage, and is also the initial state of the analysis of the quasi-static loading damage of the brittle material in the step S4;
in the step S4, an explicit algorithm carried by ABAQUS is used for analyzing the quasi-static loading damage of the brittle material;
in step S5, after the simulation analysis result in step S4 is obtained, the simulation result and the test result are compared and analyzed, so that the quasi-static damage characteristic of the brittle material can be accurately predicted, and further, the quasi-static damage characteristic of the brittle material can be further and deeply studied and predicted with high efficiency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011605539.1A CN112711884B (en) | 2020-12-29 | 2020-12-29 | Simulation prediction method for quasi-static damage of brittle material by combination of implicit algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011605539.1A CN112711884B (en) | 2020-12-29 | 2020-12-29 | Simulation prediction method for quasi-static damage of brittle material by combination of implicit algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112711884A CN112711884A (en) | 2021-04-27 |
| CN112711884B true CN112711884B (en) | 2024-01-23 |
Family
ID=75547096
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011605539.1A Active CN112711884B (en) | 2020-12-29 | 2020-12-29 | Simulation prediction method for quasi-static damage of brittle material by combination of implicit algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112711884B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104820731A (en) * | 2015-04-01 | 2015-08-05 | 华东建筑设计研究院有限公司 | Method for analyzing anti-progressive collapse performance of large-span spatial structure |
| JP2016031314A (en) * | 2014-07-29 | 2016-03-07 | 新日鐵住金株式会社 | Prediction method and device of ductile brittle fracture characteristic of member made of thin steel sheet, program therefor and recording medium |
| CN109558692A (en) * | 2018-12-25 | 2019-04-02 | 中国石油大学(华东) | Predict the finite element method of metal component residual stress and martensitic traoformation under particle high speed impact |
| WO2019238451A1 (en) * | 2018-06-13 | 2019-12-19 | Danmarks Tekniske Universitet | A method and a system for modelling and simulating a fractured geological structure |
| CN112084664A (en) * | 2020-09-14 | 2020-12-15 | 北京卫星环境工程研究所 | Sawtooth-shaped cohesion model, construction method and composite material damage simulation method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9274036B2 (en) * | 2013-12-13 | 2016-03-01 | King Fahd University Of Petroleum And Minerals | Method and apparatus for characterizing composite materials using an artificial neural network |
-
2020
- 2020-12-29 CN CN202011605539.1A patent/CN112711884B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016031314A (en) * | 2014-07-29 | 2016-03-07 | 新日鐵住金株式会社 | Prediction method and device of ductile brittle fracture characteristic of member made of thin steel sheet, program therefor and recording medium |
| CN104820731A (en) * | 2015-04-01 | 2015-08-05 | 华东建筑设计研究院有限公司 | Method for analyzing anti-progressive collapse performance of large-span spatial structure |
| WO2019238451A1 (en) * | 2018-06-13 | 2019-12-19 | Danmarks Tekniske Universitet | A method and a system for modelling and simulating a fractured geological structure |
| CN109558692A (en) * | 2018-12-25 | 2019-04-02 | 中国石油大学(华东) | Predict the finite element method of metal component residual stress and martensitic traoformation under particle high speed impact |
| CN112084664A (en) * | 2020-09-14 | 2020-12-15 | 北京卫星环境工程研究所 | Sawtooth-shaped cohesion model, construction method and composite material damage simulation method |
Non-Patent Citations (4)
| Title |
|---|
| Experimental determination of failure strength in automotive windscreens using acoustic emission and fractography;C. Brokmann et al;《Springer》;20181030;第229-236页 * |
| Structural analysis of failure behavior of laminated glass;Giulio Castori et al;《ElSEVIER:Composites》;20170323;第89-95页 * |
| 理查德·布洛克利 等.航空航天科技出版工程 3 结构技术.北京理工大学出版社,2016,第344-345页. * |
| 结构应力对飞机壁板冲击损伤影响的仿真分析;侯满义等;《弹箭与制导学报》;20101015;第30卷(第05期);第224-226页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112711884A (en) | 2021-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109684693B (en) | Method for predicting post-buckling of reinforced wallboard based on finite element analysis | |
| Dahmani et al. | Crack identification in reinforced concrete beams using ANSYS software | |
| CN111753373B (en) | Collision simulation method of automobile foamed plastic material | |
| CN113111438B (en) | On-orbit operation and maintenance method and system for digital twin satellite | |
| JP2004069638A (en) | Method for predicting crack development of elasto-plastic body and deformation predicting method | |
| CN112881174A (en) | Rock double-torsion test piece crack subcritical propagation fracture energy simulation method | |
| CN102799728B (en) | The method for making of leaf spring Dynamics Simulation Model | |
| CN108919670B (en) | Semi-physical simulation test system and test method for waste incineration control system | |
| CN112711884B (en) | Simulation prediction method for quasi-static damage of brittle material by combination of implicit algorithm | |
| CN112597684A (en) | Method for extracting static load of power assembly support 28 under working condition | |
| CN103605068A (en) | System and method for multifunctional and automated testing of keyboard | |
| CN112464401B (en) | Accurate modeling method for metal material welding spot | |
| CN110502783A (en) | Satellite vapour image analogue system and method | |
| CN113532884A (en) | Intelligent driving and ADAS system test platform and test method | |
| CN110362885B (en) | Simulation analysis method for crushing top of truck cab | |
| CN102183896B (en) | A kind of system utilizing Matlab to test nuclear power station relay protective scheme | |
| CN112648017B (en) | Numerical method for simulating serious bias large deformation of tunnel | |
| CN113591219A (en) | Multi-working-condition CAE transient fatigue analysis loading method for whole vehicle | |
| CN110472307B (en) | Optical cement stripping simulation method for display system | |
| CN112100873A (en) | Method for determining bearing capacity of hydraulic building | |
| CN117407031B (en) | Patch package upgrade file generation method, device, equipment and storage medium | |
| Liu et al. | A robust interface finite element formulation for modeling brittle material failure problems | |
| CN114912316A (en) | Method for optimizing thickness of automobile part under fixed working condition | |
| CN115758823A (en) | Target equipment simulation method, target equipment simulation device, target equipment generation system and storage medium | |
| CN117454684A (en) | Modeling method and system for numerical simulation of tunnel concrete cracking characteristics |
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 |