125 results about "Strain energy" patented technology

According to the bearing rigidity optimization and vibration reduction design method for the rotor and the bearing system, under the condition that the structural form of the rotor and the bearing system is not subversively adjusted, the elastic bearing and the extrusion oil film damper are arranged and designed, and the bearing rigidity of the rotor and the bearing system is optimized. The fulcrum rigidity is reduced to change the strain energy distribution of the rotor and adjust the critical rotating speed, and the vibration level of the rotor, especially vibration caused by sudden unbalance and transcritical of the multi-fulcrum rotor, is suppressed.

The invention discloses a heavy-duty gas turbine blade creep fatigue life design method considering an anisotropic effect, and relates to the technical field of gas turbine key component creep fatigue life evaluation.The method comprises the steps that orientation factor parameters on characteristic orientation are determined based on tensile data on [001] crystal orientation; constructing a modified stress relaxation formula; constructing an anisotropic failure strain energy density equation based on the pure creep performance data and the orientation factor parameters, and determining an upper platform of failure inelastic strain energy density on [001] crystal orientation; further determining an upper platform of the failure inelastic strain energy density on the feature orientation; constructing a creep damage prediction model based on a modified stress relaxation formula and an anisotropic failure strain energy density equation; and adopting the creep damage prediction model to obtain the creep fatigue life of the anisotropic turbine blade material sample. The method meets the actual requirements of service life evaluation of the high-temperature part of the heavy-duty gas turbine.

The invention discloses a method for predicting mechanical properties of ultra-high performance concrete, particularly relates to the technical field of dynamic property analysis of building materials, and aims to solve the problem that sudden brittle failure of materials cannot be accurately predicted under impact and explosion loads in the conventional method. The method comprises the following steps: calculating a macroscopic stress-strain response based on a continuum mechanical equation by acquiring a material ratio parameter and a dynamic load condition; constructing topological features through the strain energy entropy change rate to locate a damage area; determining an interface debonding energy threshold value and a crack growth rate of the damaged area in combination with interface fracture mechanics; utilizing the law of conservation of energy to establish dynamic correlation between interface debonding energy and strain energy density; when the crack growth rate exceeds an interface debonding energy threshold value, correcting macroscopic response data of a sudden stress drop interval based on the dynamic association relationship; the cross-scale coupling of a microscopic damage mechanism and a macroscopic mechanical behavior is realized, and the prediction precision under a dynamic load is remarkably improved.

The invention relates to the technical field of aviation part manufacturing and metal punch forming, in particular to an aviation volute spiral spring punch forming method which comprises the following steps: a data acquisition step: acquiring real-time modulus data and a target geometric trajectory equation of a to-be-processed plate; an energy mapping calculation step: performing energy mapping on the trajectory equation, and determining a target strain energy density distribution diagram representing the energy dissipation gradient; a potential barrier parameter determination step: determining a deformation energy potential barrier parameter of the current forming position in combination with the real-time modulus; a rheology compensation resolving step: resolving a rheology compensation control instruction set containing ultrasonic excitation and piezoelectric driving parameters based on potential barrier parameters; a dynamic forming execution step: driving the dynamic rheological compensation die system to execute stamping operation; according to the method, dimension raising from geometric shape copying to energy flow control is achieved, the fatigue life of the spring is remarkably prolonged, and geometric consistency of the spring is remarkably improved.

The invention belongs to the technical field of aircraft thin-wall structure design, and particularly relates to a method for correcting a quadrilateral shell unit in aircraft thin-wall structure finite element analysis, which comprises the following steps of: 1, establishing a unit coordinate system for the quadrilateral shell unit; step 2, under the unit coordinate system, establishing an equilibrium equation of a plane stress unit related to in-plane deformation, and calculating a rigidity matrix of the plane stress unit; step 3, under the unit coordinate system, establishing a balance equation of a bending plate unit related to out-of-plane deformation, and calculating a rigidity matrix of the bending plate unit; 4, applying an in-plane rotation penalty term on the rotational stiffness of the quadrilateral shell unit, and calculating the influence of the in-plane rotation penalty term on the strain energy function of the quadrilateral shell unit under the unit coordinate system to obtain an in-plane rotation stiffness matrix; 5, calculating a stiffness matrix of the quadrilateral shell unit under the unit coordinate system; and step 6, coordinate transformation is carried out to obtain a stiffness matrix of the quadrilateral shell unit under the overall coordinate system.

The invention relates to the technical field of oil and gas reservoir development and geological modeling, and discloses a multi-scale factor constraint-based tectonic network fracture modeling method, which comprises the following steps of: synthesizing a Riemannian metric tensor field on the basis of earthquake, logging and geomechanics data, and defining non-Euclidean distance cost of a fracture expanded in an anisotropic medium; initial seed points are screened according to the elastic strain energy density, and initial growth potential energy in a limited range is distributed; solving the eikonal equation by using an anisotropic fast marching algorithm to carry out wavefront competitive growth, and dividing a grid region into generalized Voronoi units; identifying a wavefront contact interface, and extracting gradient features to judge a fusion or truncation type so as to establish fracture topological connection; and finally, tracking a geodesic line path along an anti-gradient direction to generate a three-dimensional discrete fracture network. According to the method, macro and micro constraints are unified through Riemannian geometry, clear physical significance is given to the fracture by utilizing an energy mechanism, and automatic and accurate construction of the complex fracture network topology structure is realized.

The invention discloses a roadway support digital design method and system based on surrounding rock fracture evolution simulation, particularly relates to the technical field of roadway engineering digital twin modeling and optical fiber stress sensing, and is used for solving the problems that in an existing support design method, the evolution process of a digital model and a physical entity is asynchronous, and the simulation precision is insufficient. Real-time surrounding rock strain data are obtained through a distributed optical fiber sensor, actual ground stress field distribution is obtained based on strain data inversion, the dominant direction of a potential fracture surface is identified, an anisotropic strength parameter field of a mechanical model is updated, and the updated model is operated to simulate the surrounding rock fracture evolution process. And identifying a composite high-risk area by combining strain energy density distribution and crustal stress field characteristics, and finally generating accurate digital support design parameters to realize dynamic optimization and accurate design of a roadway support scheme.

This invention belongs to the field of aerospace composite material manufacturing technology, specifically relating to a method and equipment for forming thermoplastic composite strain energy rods. First, prepreg is laid up to form a preform, which is then pressed by a flatbed pre-pressing device to form a preform sheet intermediate. Subsequently, the preform intermediate is fed into a hot-pressing module with a central heating and two-end cooling gradient temperature field. It undergoes melting and molding in a high-temperature zone and controlled cooling in a low-temperature zone to achieve final shaping. Finally, the formed rod is cooled and shaped, then pulled out at a uniform speed by a traction device and coiled and collected by a winding device, achieving continuous production. The hot-pressing module of the supporting hot-pressing equipment is equipped with a zoned heating and integrated cooling system, which can precisely construct a gradient temperature field and control the cooling rate. This invention, by using gradient temperature control to precisely regulate the material's crystallization behavior and release internal stress, effectively ensures the straightness, dimensional stability, and consistent deployment performance of the strain energy rod, perfectly meeting the manufacturing requirements of strain energy rods for large deployable structures in spacecraft.

This application provides a method for coating positioning detection and control of sponge substrate, relating to the field of flexible composite material manufacturing, including the following steps: feeding and conveying the sponge substrate, and controlling the tension and lateral position of the sponge substrate; unwinding the coating material, applying adhesive, and preheating; acquiring the full-field deformation field of the sponge substrate surface through a binocular high-speed vision unit, calculating the strain energy density distribution based on nonlinear strain analysis, and generating a deformation risk digital map with energy release rate characteristics; based on the energy release rate characteristics in the deformation risk digital map, calculating control parameters through a dynamic thermodynamic model, and performing local thermal adaptation control of the coating material to match the energy field; aligning and pressing the coating material, which has undergone local thermal adaptation control to match the energy field, with the sponge substrate; and curing and quality verification of the pressed composite material.

The invention relates to an engine component fatigue life evaluation method based on a structural simulation part, and belongs to the technical field of simulation data processing, and the method comprises the following steps: building a tensile strain energy life prediction model based on elastic strain energy and plastic strain energy, and predicting the fatigue life of a structural part based on the model; optimizing the design parameters of the simulation part of the engine component, manufacturing the simulation part based on the design parameters, performing a low-cycle fatigue test on the simulation part, and analyzing the fatigue life of the simulation part; comparing the fatigue life of the simulation part with the predicted fatigue life of the structural part, and if the requirements are not met, re-optimizing the design parameters of the simulation part; and if the requirement is met, outputting the fatigue life of the structural member. According to the method, the fatigue life of the simulation part is compared with the predicted fatigue life of the structural part, the model is further improved, a closed-loop research method is formed, the accuracy and efficiency of fatigue life prediction are improved, the limitation of a traditional hysteresis energy method is overcome, and the research cost and period are reduced.

The application discloses a rotor coupling dynamics calculation method of complex three-dimensional shape layout. A coordinate conversion relationship caused by a front protrusion / after-sweeping / downward-reversing configuration is established; nonlinear elastic deformation of elastic blades waving, oscillating and twisting and nonlinear coupling relationship among them are described; nonlinear motion caused by the complex three-dimensional shape layout is described based on the coordinate conversion relationship and the nonlinear coupling relationship; and the relationship among strain energy, kinetic energy and virtual work of external force of the rotor is established based on the Hamilton principle according to the description of the nonlinear motion, and a finite element method of a 15-degree-of-freedom beam element is used to establish a rotor dynamics equation. The application can effectively reflect the structural coupling characteristics and coupling degree of the waving-torsion direction caused by the front / rear sweep of the blade tip, effectively reflect the structural coupling characteristics and coupling degree of the oscillation-torsion direction caused by the upward / downward reversal of the blade tip, and effectively reflect the rotor dynamics characteristic behaviors caused by the three-dimensional shape of the blade.

The embodiment of the application discloses a bridge anti-ship collision device impact force-impact depth model construction method, comprising the following steps: determining the total potential energy of the system in the process of the barge impacting the anti-collision device, wherein the total potential energy comprises the plastic strain energy of the outer panel of the anti-collision device, the compression energy consumption of the internal structure of the anti-collision device, the energy absorption of the bow of the barge and the work of the external impact force; each item in the total potential energy is expressed as a function of the impact force and / or the impact depth of the anti-collision device; each function is substituted into the minimum potential energy principle equation satisfied by the total potential energy, and an impact force-impact depth model for representing the relationship between the impact force and the impact depth of the anti-collision device is obtained. The embodiment can establish a clear mechanism model for the dynamic response of the anti-collision device under the ship collision.

This invention discloses an online monitoring method for vibration damage in building structures, specifically relating to the field of engineering structural health monitoring technology. It addresses the problem of insufficient damage location accuracy in existing vibration monitoring methods due to sparse measurement points. The method involves collecting structural vibration response data, first analyzing the overall input impedance spectrum to determine stiffness degradation, then calculating the strain energy density distribution and dividing the area to analyze synergy indicators for preliminary damage location. Next, the propagation path is constructed using the structural symmetry characteristics to verify damage reliability. Finally, precise damage location is achieved through the analysis of abrupt changes in local curvature mode distribution.

The application provides a simulation analysis method for compression deformation characteristics of a rubber pad, and relates to the field of rubber pad product testing. The application constructs a finite element grid model based on geometric parameters and initial material parameters of the rubber pad, identifies geometric singularity regions with displacement gradient mutation; applies a test loading load on the model, extracts displacement path curvature and performs gradient correction on material parameters, and constructs a local material constitutive relation; performs compression deformation simulation based on the local constitutive relation, analyzes strain energy accumulation sequence and determines an energy transmission control region; determines a dominant constraint mode in combination with the spatial corresponding relationship between the energy transmission control region and the geometric singularity region, and performs incremental correction on the simulation result, so as to realize accurate prediction of non-uniform compression deformation characteristics of the rubber pad.

The application discloses a kind of non-saturated solidified soil fatigue cracking damage prediction methods based on energy-mechanics method.The application introduces the effective stress formula proposed by Lytton and gives the solution method of relevant parameters in combination with energy-mechanics method, obtains the damage density evolution trend of non-saturated solidified soil material, finally based on the theory of fracture mechanics and continuous damage mechanics, on the basis of Paris equation, introduce damage density function, and based on dissipation strain energy DSE, the expression form of J integral is obtained, so as to establish the damage evolution prediction model capable of characterizing the multi-crack cracking process.The application solves the problem that the previous solidified soil fatigue cracking model cannot effectively reflect the influence of unsaturated characteristics on soil material cracking, and the model is based on the basis of mechanics theory, can better predict the multi-crack fatigue cracking damage evolution process from the mechanism level of non-saturated solidified soil material fatigue cracking, and provides basis and guidance for the design of solidified soil material and solidified soil road structure.

The application discloses a kind of early warning methods of easy argillization rock instability mutation, including making rock sample, respectively on rock sample under the condition of uniaxial compression acoustic emission test, and obtain the displacement, load and acoustic emission parameter numerical value of each rock sample;Calculate and analyze the mutation characteristics of the evolution of elastic strain energy curve of easy argillization rock;Construct the early warning model of elastic strain energy strain sequence of easy argillization rock instability mutation, when the simultaneous difference set equation Δ>0 indicates that easy argillization rock is in stable state, Δ<0 when easy argillization rock is in unstable state, Δ=0 when easy argillization rock is in critical state.The method of the application can comprehensively, quickly and accurately predict the instability mutation failure of easy argillization rock, and provide scientific reference basis for early warning prevention and treatment of rock instability mutation of mine exploitation.

The invention provides a path evolution method for energy absorption performance prediction and structural response simulation of an aramid fiber fabric, and relates to the technical field of performance simulation, and the method comprises the steps: obtaining three-dimensional weaving structure data and material characteristic parameters of the aramid fiber fabric, discretizing the data into a node unit set, and establishing a connection topology; and identifying the direct bearing unit group according to the high-speed impact load, and generating an evolution path sequence of energy propagation according to the topological coupling strength. And calculating the stress state of each unit along the path, marking a failure unit according to a composite failure criterion and blocking energy transfer, and meanwhile, accumulating the elastic strain energy release amount and the interface friction dissipation energy, and taking the sum of the elastic strain energy release amount and the interface friction dissipation energy as an energy absorption performance prediction value. According to the method, efficient and high-precision simulation of energy absorption and damage evolution of the aramid fabric under impact is realized.

The invention belongs to the field of engine design, and particularly relates to an aero-engine whole machine dynamic structure design parameter optimization method, which utilizes an APDL language of ANSYS software to carry out parametric modeling on an aero-engine whole machine (including a rotor, a support and a stator casing), and applies additional centrifugal force and additional gyroscopic moment generated by maneuvering flight. The stress level of the bearing case, the strain energy distribution of a rotor system and the rotor and stator gap change are obtained through bearing case stress analysis, critical speed analysis and unbalanced response analysis. Based on the analysis results, multi-objective optimization is carried out on complete machine dynamic design parameters, so that the stress level of a force bearing casing is lowest, rotor strain energy distribution meets rotor dynamic design requirements, rotor and stator gap changes are minimum, the rub-impact risk is reduced, optimization of the structural design parameters is completed through software, the design efficiency is improved, and the design cost is reduced. The design period is shortened.

The invention relates to the technical field of mechanical property testing, in particular to a laminated plate mechanical property evaluation method and system based on strain data. The method comprises the following steps: performing a biaxial tensile experiment on a laminated plate sample, synchronously acquiring load time-history data and a full-field speckle image, constructing a full-field strain matrix through a mapping algorithm, dispersing the full-field strain matrix into nodes, dividing a core pure shear region and a stress diffusion region, and then obtaining first principal strain through principal component analysis; a load transmission efficiency attenuation spectrum and full-field strain energy entropy density are calculated respectively, a dynamic correction model is constructed in combination with weighted coupling of preset shape characteristic parameters, and a real constitutive curve is obtained by using the model to calculate correction coefficient inversion nominal stress; the thermodynamic entropy characteristic is introduced, so that the problem that the stress field distortion at the nonlinear damage stage of the material is difficult to sense by a traditional method is solved, and the biaxial mechanical property test precision of the laminated plate is effectively improved.

The invention provides an engine composite material blade vibration characteristic prediction method. The method comprises the following steps: establishing a blade geometric model; establishing a coordinate system according to the blade geometric model, and obtaining coordinates of each point of the blade geometric model; calculating the curvature radius and deflection of the blade; constructing a displacement field based on a first-order shear deformation theory according to the coordinate of any point on the blade; establishing a strain-displacement relationship of the blade model; calculating the strain energy, the kinetic energy and the centrifugal force potential energy of the blade according to the curvature radius, the deflection and the strain-displacement relation; defining a vibration mode function by adopting a Chebyshev-ritz method, setting boundary conditions of the blade, and obtaining an undetermined coefficient; and calculating a characteristic value and a characteristic vector through a driving coefficient and a generalized characteristic value equation, and outputting a vibration characteristic. According to the method, the three-dimensional geometric model of the blade is established, and the curvature and deflection of the blade in the rotating process are calculated by adopting the Chebyshev-Ritz method through the first-order shear deformation theory, so that the vibration characteristics of the blade in the rotating process are predicted.

The invention discloses a composite sealing structure stress inversion method based on CT in-situ observation and end-to-end deep learning, and relates to the technical field of stress inversion, and the method comprises the steps: obtaining measurement data, taking the measurement data as input, training a displacement network and a constitutive network under physical constraints, and obtaining a stress inversion result; a trained displacement network and a trained constitutive network are obtained, the displacement of each position point in the prediction loading step is determined by using the trained displacement network, the deformation coordinate of each position point in the prediction loading step is further determined, a displacement field is obtained, the strain invariant of each position point in the prediction loading step is further calculated, and the deformation coordinate of each position point in the prediction loading step is obtained. The trained constitutive network is used for determining the strain energy density of each position point in the prediction loading step, the stress of each position point in the prediction loading step is further determined, the stress field is obtained, and stress field inversion can be efficiently and accurately completed under the condition that the internal state of the composite sealing structure is invisible and a material constitutive model is unknown.

The invention relates to the technical field of numerical control, in particular to a real-time data driven cutting bed automatic control method and system and a storage medium. The method comprises the following steps: collecting dynamic stress data of a cutting bed cutter and microscopic deformation data of a fabric in real time through a multi-source sensing network, and constructing a physical coupling state space of the cutter and the fabric; according to the invention, through a collaborative optimization mechanism driven by multi-source sensing data, high-precision collaborative control of cutter dynamics and fabric deformation is realized in the field of cutting bed control, and through hardware synchronous triggering of the piezoelectric ceramic array sensor and the distributed optical fiber sensor, space-time alignment monitoring of cutter stress fluctuation and fabric interlayer shear deformation is realized. According to the design, the problem of time sequence mismatch of a traditional system in cross-scale data fusion is solved, the thermodynamic state of the cutter and fabric strain energy distribution form a computable incidence matrix, and a solvable space foundation of an optimization algorithm is laid.

This invention provides a method for determining the critical state of cracks in viscoelastic materials under dynamic interaction of cracking and healing. Specifically, it involves: establishing a geometric model of the energy change region under dynamic interaction of cracking and healing; constructing energy dissipation and recovery strain energy models for the apparent and actual configurations of the viscoelastic material sample during the loading stage; determining the size and net increase in damage density of the viscoelastic material sample under dynamic interaction of cracking and healing; constructing a model relating the recovery strain energy and the net increase in damage density during the healing stage; using the sum of the increase in volume free energy due to viscoelastic failure volume recovery and the decrease in surface energy due to crack healing as the total free energy increment; calculating the rate of change of the free energy increment; and determining the critical state of the crack based on this method. This invention provides a basis for determining the critical transition of crack instability propagation in materials with viscoelastic and healing properties.

A method and system for determining a critical inelastic strain rate during a creep-fatigue damage process are provided. The method includes: obtaining a plurality of high-temperature equipment materials; under a same temperature condition, conducting a creep test, a fatigue test, and a creep-fatigue interaction test on the high-temperature equipment materials, respectively; establishing a stress relaxation formula based on the creep test result; establishing a fitting equation for an inelastic strain energy density at a failure based on the fatigue test result; establishing a net tensile hysteresis energy-induced fatigue damage equation based on the creep-fatigue interaction test result; and according to a linear cumulative damage rule, determining a creep-fatigue life of the high-temperature equipment materials based on the stress relaxation formula, the fitting equation for an inelastic strain energy density at a failure, and the net tensile hysteresis energy-induced fatigue damage equation.

A shoe sole includes a shock absorber that buckles when the shock absorber receives compressive force applied in a normal direction to a ground contact surface. When a load is applied to the shoe sole, the shock absorber starts to buckle in a state in which a stress occurring in the shock absorber is within a range from 0.15 MPa to 0.80 MPa and a strain of the shock absorber is within a range from 10% to 60%. When a time point at which a strain energy density of the shock absorber reaches 0.157 J / cm3 is defined as a specific time point, a maximum value of the stress occurring in the shock absorber until the specific time point is 0.80 MPa or less, and a tangential elastic modulus of the shock absorber at the specific time point is 5.00 MPa or less.

This invention relates to the field of electronic digital data processing technology, specifically to computer-aided engineering calculation and numerical simulation data processing technology. This invention discloses a design method for determining the maximum hanging weight of an ultra-thin lifting beam. In incremental iterative solutions, the equivalent stress gradient and strain energy density are extracted in real time to calculate local instability characteristic values. When the rate of change of element strain energy density exceeds a threshold, a mesh reconstruction mechanism is triggered. The nodal displacement field and stress field are extracted as initial boundary conditions. The coarse mesh is split into a fine mesh, and the coordinates of newly added nodes are smoothly interpolated based on the minimum potential energy principle. The local and global mesh degrees of freedom are coupled and assembled to continue nonlinear buckling analysis iterations until the load-displacement curve shows an extreme point. The corresponding load is then extracted as the maximum hanging weight. This invention avoids stiffness matrix singularities and computational divergence caused by distortion of the fixed mesh in local instability regions, ensuring the computational convergence of determining the maximum hanging weight.

This invention discloses a true triaxial energy evolution analysis method for coal and rock under high static load and disturbance coupling, relating to the field of mining engineering technology. The method includes the following steps: applying high static background stress to a coal and rock sample on a six-axis servo loading platform, and applying periodic disturbance loads in predetermined loading directions, generating loading time sequence records and disturbance effect lists corresponding to the high static background stress and periodic disturbance loads. This invention performs fractional analysis of energy in each principal stress direction under true triaxial loading conditions, achieving independent characterization of strain energy, dissipated energy, and disturbance input energy, and can identify the implicit energy accumulation process under high static load. Through energy density evolution, damage mapping, and three-dimensional visualization and curve fitting, the energy release law and dynamic stability threshold are extracted, providing quantitative basis for identifying coal and rock dynamic instability and improving the interpretability and application value of experimental analysis.

This invention discloses an in-situ self-generated TiBw and TiCp synergistically reinforced bi-matrix heterostructure titanium matrix composite material prepared by powder metallurgy, and its preparation method. This invention aims to solve the problems of "strength-plasticity mismatch" and insufficient dynamic impact energy absorption efficiency in traditional titanium alloys and their composites. The composite material uses pure Ti as the soft phase, TC4 as the hard phase, and B4C particles as the reinforcing precursor. It can be formed by low-energy ball milling, vacuum drying, rapid hot pressing sintering, and heat treatment, without subsequent rolling or forging. B4C reacts in-situ with the matrix to generate TiBw (TiBw) and TiCp (TiCp), constructing a discontinuous in-situ self-generated bi-scale TiBw and TiCp synergistically reinforced heterostructure at the soft-hard phase interface, achieving a significant synergistic effect of strengthening and toughening. The prepared composite material has a quasi-static tensile strength of 1085 MPa and a fracture strain of 16.3%; at 3000 s⁻¹, the tensile strength is 1085 MPa. ‑1 The dynamic compressive strength under strain rate reaches 1562 MPa, with a true strain of 27.8% and an energy absorption efficiency of 408 MJ / m³. This composite material achieves synergistic optimization of high strength, high plasticity, and high energy absorption through mechanisms such as twinned coordinated plasticity, multi-stage whisker fracture energy dissipation, and adiabatic shear band suppression. The invention features a simple process and low cost, making it suitable for the fabrication of lightweight, high-impact-resistant structural components for aerospace, armor protection, and other applications.

The invention belongs to the technical field of potential earthquake risk prediction analysis, and discloses an earthquake risk analysis method considering historical earthquake energy release, which comprises the following steps: constructing a strain energy three-dimensional grid of a research area; calculating historical earthquake energy release and a release range; calculating accumulated strain energy of each grid; calculating residual strain energy and identifying a potential earthquake generation area; and earthquake risk analysis. According to the method, an original analysis method of a numerical simulation result is broken through, energy release of historical earthquakes and energy accumulation between earthquakes are considered from a new perspective and method, and a new thought for carrying out earthquake risk analysis based on the numerical simulation result is provided.